  We study the vertical perturbations in the galactic disc of the Milky
Way-size high-resolution hydrodynamical cosmological simulation named GARROTXA.
We detect phase spirals in the vertical projection $Z- V_{Z}$ of disc's stellar
particles for the first time in this type of simulations. Qualitatively similar
structures were detected in the recent Gaia data, and their origin is still
under study. In our model the spiral-like structures in the phase space are
present in a wide range of times and locations across the disc. By accounting
for an evolving mix of stellar populations, we observe that, as seen in the
data, the phase spirals are better observed in the range of
younger-intermediate star particles. We measure the intensity of the spiral
with a Fourier decomposition and find that these structures appear stronger
near satellite pericenters. Current dynamical models of the phase spiral
considering a single perturber required a mass at least of the order of
10$^{10}$ M$_\odot$, but all three of our satellites have masses of the order
of $\sim$10$^8$ M$_\odot$. We suggest that there are other mechanisms at play
which appear naturally in our model such as the physics of gas, collective
effect of multiple perturbers, and a dynamically cold population that is
continuously renovated by the star formation Complementing collisionless
isolated N-body models with the use of fully-cosmological simulations with
enough resolution can provide new insights into the nature/origin of the phase
spiral.
