The advent of a new generation of large-scale galaxy surveys is pushing cosmological numerical simulations in an uncharted territory. The simultaneous requirements of high resolution and very large volume pose serious technical challenges, due to their computational and data storage demand. In this paper, we present a novel approach dubbed Dynamic Zoom Simulations -- or DZS -- developed to tackle these issues. Our method is tailored to the production of lightcone outputs from N-body numerical simulations, which allow for a more efficient storage and post-processing compared to standard comoving snapshots, and more directly mimic the format of survey data. In DZS, the resolution of the simulation is dynamically decreased outside the lightcone surface, reducing the computational work load, while simultaneously preserving the accuracy inside the lightcone and the large-scale gravitational field. We show that our approach can achieve virtually identical results to traditional simulations at half of the computational cost for our largest box. We also forecast this speedup to increase up to a factor of 5 for larger and/or higher-resolution simulations. We assess the accuracy of the numerical integration by comparing pairs of identical simulations run with and without DZS. Deviations in the lightcone halo mass function, in the sky-projected lightcone, and in the 3D matter lightcone always remain below 0.1%. In summary, our results indicate that the DZS technique may provide a highly-valuable tool to address the technical challenges that will characterise the next generation of large-scale cosmological simulations.

We investigate Lyman-$\alpha$ transmission spikes at $z > 5$ in synthetic quasar spectra, and discuss their connection to the properties of the intergalactic medium and their ability to constrain reionization models. We use state-of-the-art radiation-hydrodynamic simulations from the Cosmic Reionization On Computers series to predict the number of transmission spikes as a function of redshift, both in the ideal case of infinite spectral resolution and in a realistic observational setting. Transmission spikes are produced in highly-ionized underdense regions located in the vicinity of UV sources. We find that most of the predicted spikes are unresolved by current observations, and show that our mock spectra are consistent with observations of the quasar ULAS J1120+0641 in about 15% of the realizations. The spike height correlates with both the gas density and the ionized fraction, but the former link is erased when synthetic spectra are smoothed to realistically achievable spectral resolutions. There exists a linear relationship between spike width and the extent of the associated underdense region, with a slope that is redshift-dependent. In agreement with observations, the spike transmitted flux is suppressed at small distance from bright galaxies as these reside in overdense regions. We argue that this anti-correlation can be used to constrain large-scale density modes.

The detection of an unexpectedly large number of faint candidate quasars (QSO) at $z \gtrsim 4$ has motivated Madau and Haardt to investigate a reionization scenario in which active galactic nuclei dominate the photon budget at all times. Their analytical study reveals that this picture is in agreement with the evolution of the HI volume fraction in the intergalactic medium (IGM) and the optical depth to Thomson scattering of the cosmic microwave background. We employ a suite of hydrodynamical simulations post-processed with a radiative transfer code to further investigate the properties of such a reionization history. In particular, we generate synthetic absorption-line spectra for HI and HeII that we compare with observational data. Although we confirm the analytical results mentioned above, we also find that the evolution of the IGM temperature and of the HeII optical depth are at odds with current observational constraints. Nevertheless, the QSO-dominated scenario presents the attractive feature that it naturally generates an inhomogeneous IGM and thus produces an extended tail at high values in the distribution function of the HI effective optical depth, in agreement with recent observations. We find tentative evidence that considering some QSO emission at redshift $z \gtrsim 5$ could reconcile such observations with numerical simulations that have so far failed to reproduce this feature. Finally, we show that collecting more HeII-absorption spectra at $z \gtrsim 3$ and studying their distinctive characteristics will be key to precisely constraining the QSO contribution to reionization.

The radial acceleration measured in bright galaxies tightly correlates with that generated by the observed distribution of baryons, a phenomenon known as the radial acceleration relation (RAR). Dwarf spheroidal satellite galaxies have been recently found to depart from the extrapolation of the RAR measured for more massive objects but with a substantially larger scatter. If confirmed by new data, this result provides a powerful test of the theory of gravity at low accelerations that requires robust theoretical predictions. By using high-resolution hydrodynamical simulations, we show that, within the standard model of cosmology ($\Lambda$CDM), satellite galaxies are expected to follow the same RAR as brighter systems but with a much larger scatter which does not correlate with the physical properties of the galaxies. In the simulations, the RAR evolves mildly with redshift. Moreover, the acceleration due to the gravitational field of the host has no effect on the RAR. This is in contrast with the External Field Effect in Modified Newtonian Dynamics (MOND) which causes galaxies in strong external fields to deviate from the RAR. This difference between $\Lambda$CDM and MOND offers a possible way to discriminate between them.

We use zoom hydrodynamical simulations to investigate the properties of satellites within galaxy-sized dark-matter haloes with different assembly histories. We consider two classes of haloes at redshift $z=0$: `stalled' haloes that assembled at $z>1$ and `accreting' ones that are still forming nowadays. Previously, we showed that the stalled haloes are embedded within thick filaments of the cosmic web while the accreting ones lie where multiple thin filaments converge. We find that satellites in the two classes have both similar and different properties. Their mass spectra, radial count profiles, baryonic and stellar content, and the amount of material they shed are indistinguishable. However, the mass fraction locked in satellites is substantially larger for the accreting haloes as they experience more mergers at late times. The largest difference is found in the satellite kinematics. Substructures fall towards the accreting haloes along quasi-radial trajectories whereas an important tangential velocity component is developed, before accretion, while orbiting the filament that surrounds the stalled haloes. Thus, the velocity anisotropy parameter of the satellites ($\beta$) is positive for the accreting haloes and negative for the stalled ones. This signature enables us to tentatively categorize the Milky Way halo as stalled based on a recent measurement of $\beta$. Half of our haloes contain clusters of satellites with aligned orbital angular momenta corresponding to flattened structures in space. These features are not driven by baryonic physics and are only found in haloes hosting grand-design spiral galaxies, independently of their assembly history.

The clustering of dark matter haloes with fixed mass depends on their formation history, an effect known as assembly bias. We use zoom N -body simulations to investigate the origin of this phenomenon. For each halo at redshift z=0, we determine the time in which the physical volume containing its final mass becomes stable. We consider five examples for which this happens at z~1.5 and two that do not stabilize by z=0. The zoom simulations show that early-collapsing haloes do not grow in mass at z=0 while late-forming ones show a net inflow. The reason is that 'accreting' haloes are located at the nodes of a network of thin filaments feeding them. Conversely, each 'stalled' halo lies within a prominent filament that is thicker than the halo size. Infalling material from the surroundings becomes part of the filament while matter within it recedes from the halo. We conclude that assembly bias originates from quenching halo growth due to tidal forces following the formation of non-linear structures in the cosmic web, as previously conjectured in the literature. Also the internal dynamics of the haloes change: the velocity anisotropy profile is biased towards radial (tangential) orbits in accreting (stalled) haloes. Our findings reveal the cause of the yet unexplained dependence of halo clustering on the anisotropy. Finally, we extend the excursion-set theory to account for these effects. A simple criterion based on the ellipticity of the linear tidal field combined with the spherical collapse model provides excellent predictions for both classes of haloes.

The growth-rate and the internal dynamics of galaxy-sized dark-matter haloes depend on their location within the cosmic web. Haloes that sit at the nodes grow in mass till the present time and are dominated by radial orbits. Conversely, haloes embedded in prominent filaments do not change much in size and are dominated by tangential orbits. Using zoom hydrodynamical simulations including star formation and feedback, we study how gas accretes onto these different classes of objects that, for simplicity, we dub 'accreting' and 'stalled' haloes. We find that all haloes get a fresh supply of newly accreted gas in their inner regions, although this slowly decreases with time, in particular for the stalled haloes. The inflow of new gas is always higher than (but comparable with) that of recycled material. Overall, the cold-gas fraction increases (decreases) with time for the accreting (stalled) haloes. In all cases, a stellar disc and a bulge form at the centre of the simulated haloes. The total stellar mass is in excellent agreement with expectations based on the abundance-matching technique. Many properties of the central galaxies do not seem to correlate with the large-scale environment in which the haloes reside. However, there are two notable exceptions that characterise stalled haloes with respect to their accreting counterparts: i) the galaxy disc contains much older stellar populations; ii) its vertical scale-height is larger by a factor of two or more. This thickening is likely due to the heating of the long-lived discs by mergers and close flybys.

We perform for the first time high-resolution zoom-in re-simulations of individual halos in the context of the Multi-coupled Dark Energy (McDE) scenario, which is characterised by the existence of two distinct dark matter particle species with opposite couplings to a Dark Energy scalar field. We compare the structural properties of the simulated halos to the standard Lambda-CDM results. The zoomed-in initial conditions are set up using a specifically designed code called ZInCo that we publicly release along with the present paper. Our numerical results allow to investigate in detail and with unprecedented resolution the halo segregation process that characterises McDE cosmologies from its very early stages. In particular, we find that in contrast to what could be inferred from previous numerical analysis at lower resolution, the segregation process is already in place at redshifts as high as z ~ 7. Most remarkably, we find that the subsequent evolution of the segregation leads to the formation of cored total matter density profiles with a core size that progressively increases in time. The shape of the cored profiles can be accurately predicted as the superposition of two NFW profiles with an increasing offset, thereby confirming the interpretation of the simulations results in terms of the segregation of the two dark matter components of the halo as a consequence of their different coupling to the Dark Energy field.