We study the distribution of the Milky Way satellites stellar and dark matter debris. For the first time we address the question of the tidal disruption of satellites in simulations by utilising simultaneously a) a realistic set of orbits extracted from cosmological simulations, b) a three component host galaxy with live halo, disc and bulge components, and c) satellites from hydrodynamical simulations. We analyse the statistical properties of the satellite debris of all massive galaxies reaching the inner Milky Way on a timescale of 2 Gyr. Up to 80$\%$ of the dark matter is stripped from the satellites, while this happens for up to 30$\%$ of their stars. The stellar debris ends mostly in the inner Milky Way halo, whereas the dark matter debris shows a flat mass distribution over the full main halo. The dark matter debris follows a density profile with inner power law index $\alpha_{\rm DM}=-0.66$ and outer index $\beta_{\rm DM}=2.94$, while for stars $\alpha_{*}=-0.44$ and $\beta_{*}=6.17$. In the inner 25 kpc, the distribution of the stellar debris is flatter than that of the dark matter debris and the orientations of their short axes differ significantly. Changing the orientation of the stellar disc by 90$^{\rm{o}}$ has only a minor impact on the distribution of the satellite debris. Our results indicate that the dark matter is more easily stripped than stars from the Milky Way satellites. The structure of the debris is dominated by the satellite orbital properties. The radial profiles, the flattening and the orientation of the stellar and dark matter debris are significantly different, which prevents the prediction of the dark matter distribution from the observed stellar component.

Aims. We aim to show how encounters with low-mass satellite galaxies may alter the bar formation in a Milky Way-like disc galaxy. Methods. We use high-resolution N-body simulations of a disc galaxy prone to mild bar instability. For realistic initial conditions of satellites, we take advantage of cosmological simulations of Milky Way-like dark matter haloes. Results. The satellites may have a significant impact on the time of bar formation. Some runs with satellites demonstrate a delay, while others show an advancement in bar formation compared to the isolated run, with such time differences reaching $\sim$ 1 Gyr. Meanwhile, the final bar configuration, including its very appearance and the bar characteristics such as the pattern speed and the exponential growth rate of its amplitude are independent of the number of encounters and their orbits. The contribution of satellites with masses below $10^9 M_{\odot}$ is insignificant, unless their pericentre distances are small. We suggest that the encounters act indirectly via inducing perturbations across the disc that evolve to delayed waves in the central part and interfere with an emerging seed bar. The predicted effect for the present-day host galaxy is expected to be even more significant at redshifts $z \gtrsim 0.5$.

We present a high resolution study of the impact of realistic satellite galaxies, extracted from cosmological simulations of Milky Way haloes including 6 Aquarius suites and Via Lactea \rom{2}, on the dynamics of the galactic disc. The initial conditions for the multi-component Milky Way galaxy were generated using the GalIC code, to ensure a system in dynamical equilibrium state prior to addition of satellites. Candidate subhaloes that came closer than 25\,kpc to the centre of the host DM haloes with initial mass enclosed within the tidal radius, $M_\textrm{tid}$ $\ge$ 10$^{8} M_{\odot}$\,=\,0.003 $M_\textrm{disc}$, were identified, inserted into our high resolution N-body simulations and evolved for 2 Gyr. We quantified the vertical heating due to such impacts by measuring the disc thickness and squared vertical velocity dispersion $\sigma_{z}^{2}$ across the disc. According to our analysis the strength of heating is strongly dependent on the high mass end of the subhalo distribution from cosmological simulations. The mean increase of the vertical dispersion is $\sim$ 20\,km$^{2}$\,s$^{-2}$\,Gyr$^{-1}$ for R $>$ 4\,kpc with a flat radial profile while, excluding Aq-F2 results, the mean heating is $<$ 12\,km$^{2}$\,s$^{-2}$\,Gyr$^{-1}$, corresponding to 28\% and 17\% of the observed vertical heating rate in the solar neighbourhood. Taking into account the statistical dispersion around the mean we miss the observed heating rate by more than 3$\sigma$. We observed a general flaring of the disc height in the case of all 7 simulations in the outer disc.