Hydrodynamical cosmological simulations are increasing their level of realism by considering more physical processes and having greater resolution or larger statistics. However, usually either the statistical power of such simulations or the resolution reached within galaxies are sacrificed. Here, we introduce the NewHorizon project in which we simulate at high resolution a zoom-in region of $\sim(16\,\rm Mpc)^3$ that is larger than a standard zoom-in region around a single halo and is embedded in a larger box. A resolution of up to 34 pc is reached within galaxies; this allows the simulation to capture the multi-phase nature of the interstellar medium and the clumpy nature of the star formation process in galaxies. In this introductory paper, we present several key fundamental properties of galaxies and their black holes, including the galaxy mass function, cosmic star formation rate, galactic metallicities, the Kennicutt-Schmidt relation, the stellar-to-halo mass relation, galaxy sizes, stellar kinematics and morphology, gas content within galaxies and its kinematics, and the black hole mass and spin properties over time. The various scaling relations are broadly reproduced by NewHorizon with some differences with the standard observables. Owing to its exquisite spatial resolution, NewHorizon captures the inefficient process of star formation in galaxies, which evolve over time from being more turbulent, gas rich, and star bursting at high redshift. These high-redshift galaxies are also more compact, and they are more elliptical and clumpier until the level of internal gas turbulence decays enough to allow for the formation of discs. The NewHorizon simulation gives access to a broad range of galaxy formation and evolution physics at low-to-intermediate stellar masses, which is a regime that will become accessible in the near future through surveys such as the LSST.

Developing an efficient code for large, multiscale astrophysical simulations is crucial in preparing the upcoming era of exascale computing. RAMSES is an astrophysical simulation code that employs parallel processing based on the Message Passing Interface (MPI). However, it has limitations in computational and memory efficiency when using a large number of CPU cores. The problem stems from inefficiencies in workload distribution and memory allocation that inevitably occur when a volume is simply decomposed into domains equal to the number of working processors. We present RAMSES-yOMP, which is a modified version of RAMSES designed to improve parallel scalability. Major updates include the incorporation of Open Multi-Processing into the MPI parallelization to take advantage of both the shared and distributed memory models. Utilizing this hybrid parallelism in high-resolution benchmark simulations with full prescriptions for baryonic physics, we achieved a performance increase of a factor of 2 in the total run-time, while using 75% less memory and 30% less storage compared to the original code, when using the same number of processors. These improvements allow us to perform larger or higher-resolution simulations than what was feasible previously.

Dust plays a pivotal role in shaping the observed morphology of galaxies. While traditional cosmological simulations often assume a fixed dust-to-gas (DTG) or dust-to-metal (DTM) mass ratio to model dust effects, recent advancements have enabled on-the-fly (OTF) dust modeling that captures the spatial and temporal evolution of dust. In this work, we investigate the impact of dust modeling on galaxy morphology using the NewCluster simulation, which implements a detailed OTF dust model. We generate mock images of NewCluster galaxies under both OTF and fixed DTM models using the radiative transfer code SKIRT, and compare their morphology to JWST observations. We measure morphology indices and use the $G-M_{20}$ test to classify galaxies. We find that the OTF galaxy models exhibit brighter centers and more pronounced bulges than those of the fixed DTM models, resulting in a lower late-type galaxy (LTG) fraction, particularly at high redshifts. This central brightening is linked to a phenomenon we refer to as the DTM cavity, a localized depression in the DTM ratio driven by intense bulge starbursts. Our results highlight the importance of modeling dust evolution in a physically motivated manner, as fixed DTM models fail to capture key morphological features.

We present the Extreme-Horizon (EH) cosmological simulation: EH models galaxy formation with stellar and AGN feedback and uses a very high resolution in the intergalactic and circumgalactic medium. The high resolution in low-density regions results in smaller-size massive galaxies at redshift $z=2$, in better agreement with observations compared to other simulations. This results from the improved modeling of cold gas flows accreting onto galaxies. Besides, the EH simulation forms a population of particularly compact galaxies with stellar masses of $10^{10-11}$\,M$_\sun$ that are reminiscent of observed ultracompact galaxies at $z\simeq2$. These objects form mainly through repeated major mergers of low-mass progenitors, independently of baryonic feedback mechanisms. This formation process can be missed in simulations using a too low resolution in low-density intergalactic regions.

Understanding the processes that trigger morphological transformation is central to understanding how and why the Universe transitions from being disc-dominated at early epochs to having the morphological mix that is observed today. We use Horizon-AGN, a cosmological hydrodynamical simulation, to perform a comprehensive study of the processes that drive morphological change in massive (M > 10^10 MSun) galaxies over cosmic time. We show that (1) essentially all the morphological evolution in galaxies that are spheroids at z=0 is driven by mergers with mass ratios greater than 1:10, (2) major mergers alone cannot produce today's spheroid population -- minor mergers are responsible for a third of all morphological transformation over cosmic time and are its dominant driver after z~1, (3) prograde mergers trigger milder morphological transformation than retrograde mergers -- while both types of events produce similar morphological changes at z>2, the average change due to retrograde mergers is around twice that due to their prograde counterparts at z~0, (4) remnant morphology depends strongly on the gas fraction of a merger, with gas-rich mergers routinely re-growing discs, and (5) at a given stellar mass, discs do not exhibit drastically different merger histories from spheroids -- disc survival in mergers is driven by acquisition of cold gas (via cosmological accretion and gas-rich interactions) and a preponderance of prograde mergers in their merger histories.

The origin of the disk and spheroid of galaxies has been a key open question in understanding their morphology. Using the high-resolution cosmological simulation, New Horizon, we explore kinematically decomposed disk and spheroidal components of 144 field galaxies with masses greater than $\rm 10^9\,M_{\odot}$at$z=0.7$. The origins of stellar particles are classified according to their birthplace (in situ or ex situ) and their orbits at birth. Before disk settling, stars form mainly through chaotic mergers between proto-galaxies and become part of the spheroidal component. When disk settling starts, we find that more massive galaxies begin to form disk stars from earlier epochs; massive galaxies commence to develop their disks at $z\sim1-2$, while low-mass galaxies do after $z\sim1$. The formation of disks is affected by accretion as well, as mergers can trigger gas turbulence or induce misaligned gas infall that prevents galaxies from forming co-rotating disk stars. The importance of accreted stars is greater in more massive galaxies, especially in developing massive spheroids. A significant fraction of the spheroids comes from the disk stars that are perturbed, which becomes more important at lower redshifts. Some ($\sim12.5\%$) of our massive galaxies develop counter-rotating disks from the gas infall misaligned with the existing disk plane, which can last for more than a Gyr until they become the dominant component, and flip the angular momentum of the galaxy in the opposite direction. The final disk-to-total ratio of a galaxy needs to be understood in relation to its stellar mass and accretion history. We quantify the significance of the stars with different origins and provide them as guiding values.

Gravitational lensing allows to quantify the angular distribution of the convergence field around clusters of galaxies to constrain their connectivity to the cosmic web. We describe in this paper the corresponding theory in Lagrangian space where analytical results can be obtained by identifying clusters to peaks in the initial field. We derive the three-point Gaussian statistics of a two-dimensional field and its first and second derivatives. The formalism allows us to study the statistics of the field in a shell around a central peak, in particular its multipolar decomposition. The peak condition is shown to significantly remove power from the dipolar contribution and to modify the monopole and quadrupole. As expected, higher order multipoles are not significantly modified by the constraint. Analytical predictions are successfully checked against measurements in Gaussian random fields. The effect of substructures and radial weighting is shown to be small and does not change the qualitative picture. The non-linear evolution is shown to induce a non-linear bias of all multipoles proportional to the cluster mass.We predict the Gaussian and weakly non-Gaussian statistics of multipolar moments of a two-dimensional field around a peak as a proxy for the azimuthal distribution of the convergence field around a cluster of galaxies. A quantitative estimate of this multipolar decomposition of the convergence field around clusters in numerical simulations of structure formation and in observations will be presented in two forthcoming papers.

Intrinsic alignments of galaxies are a significant astrophysical systematic affecting cosmological constraints from weak gravitational lensing. Obtaining numerical predictions from hydrodynamical simulations of expected survey volumes is expensive, and a cheaper alternative relies on populating large dark matter-only simulations with accurate models of alignments calibrated on smaller hydrodynamical runs. This requires connecting the shapes and orientations of galaxies to those of dark matter haloes and to the large-scale structure. In this paper, we characterise galaxy-halo alignments in the Horizon-AGN cosmological hydrodynamical simulation. We compare the shapes and orientations of galaxies in the redshift range $0

Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: Horizon-AGN with a large volume and low resolution which tracks the high-mass (> 1e7 Msun) MBH population, and NewHorizon with a smaller volume but higher resolution that traces the low-mass (< 1e7 Msun) MBH population. While Horizon-AGN can be used to estimate the rate of inspirals for Pulsar Timing Arrays, NewHorizon can investigate MBH mergers in a statistical sample of dwarf galaxies for LISA, which is sensitive to low-mass MBHs. We use the same method to analyse the two simulations, post-processing MBH dynamics to account for time delays mostly determined by dynamical friction and stellar hardening. In both simulations, MBHs typically merge long after the galaxies do, so that the galaxy morphology at the time of the MBH merger is no longer determined by the galaxy merger from which the MBH merger originated. These time delays cause a loss of high-z MBH coalescences, shifting the peak of the MBH merger rate to z~1-2. This study shows how tracking MBH mergers in low-mass galaxies is crucial to probing the MBH merger rate for LISA and investigate the properties of the host galaxies.

Recent integral field spectroscopy observations have found that about 11% of galaxies show star-gas misalignment. The misalignment possibly results from external effects such as gas accretion, interaction with other objects, and other environmental effects, hence providing clues to these effects. We explore the properties of misaligned galaxies using Horizon-AGN, a large-volume cosmological simulation, and compare the result with the result of the Sydney-AAO Multi-object integral field spectrograph (SAMI) Galaxy Survey. Horizon-AGN can match the overall misalignment fraction and reproduces the distribution of misalignment angles found by observations surprisingly closely. The misalignment fraction is found to be highly correlated with galaxy morphology both in observations and in the simulation: early-type galaxies are substantially more frequently misaligned than late-type galaxies. The gas fraction is another important factor associated with misalignment in the sense that misalignment increases with decreasing gas fraction. However, there is a significant discrepancy between the SAMI and Horizon-AGN data in the misalignment fraction for the galaxies in dense (cluster) environments. We discuss possible origins of misalignment and disagreement.

Mapping of the large-scale structure through cosmic time has numerous applications in the studies of cosmology and galaxy evolution. At $z > 2$, the structure can be traced by the neutral intergalactic medium (IGM) by way of observing the Ly$\alpha$, forest towards densely-sampled lines-of-sight of bright background sources, such as quasars and star forming galaxies. We investigate the scientific potential of MOSAIC, a planned multi-object spectrograph on the European Extremely Large Telescope (ELT), for the 3D mapping of the IGM at $z \gtrsim 3$. We simulate a survey of $3 \lesssim z \lesssim 4$ galaxies down to a limiting magnitude of $m_{r}\sim 25.5$ mag in an area of 1 degree$^2$ in the sky. Galaxies and their spectra (including the line-of-sight Ly$\alpha$ absorption) are taken from the lightcone extracted from the Horizon-AGN cosmological hydrodynamical simulation. The quality of the reconstruction of the original density field is studied for different spectral resolutions and signal-to-noise ratios of the spectra. We demonstrate that the minimum $S/N$ (per resolution element) of the faintest galaxies that such survey has to reach is $S/N = 4$. We show that a survey with such sensitivity enables a robust extraction of cosmic filaments and the detection of the theoretically-predicted galaxy stellar mass and star-formation rate gradients towards filaments. By simulating the realistic performance of MOSAIC we obtain $S/N(T_{\rm obs}, R, m_{r})$ scaling relations. We estimate that $\lesssim 35~(65)$ nights of observation time are required to carry out the survey with the instrument's high multiplex mode and with the spectral resolution of $R=1000~(2000)$. A survey with a MOSAIC-concept instrument on the ELT is found to enable the mapping of the IGM at $z > 3$ on Mpc scales, and as such will be complementary to and competitive with other planned IGM tomography surveys. [abridged]

Our statistical understanding of galaxy evolution is fundamentally driven by objects that lie above the surface-brightness limits of current wide-area surveys (mu ~ 23 mag arcsec^-2). While both theory and small, deep surveys have hinted at a rich population of low-surface-brightness galaxies (LSBGs) fainter than these limits, their formation remains poorly understood. We use Horizon-AGN, a cosmological hydrodynamical simulation to study how LSBGs, and in particular the population of ultra-diffuse galaxies (UDGs; mu > 24.5 mag arcsec^-2), form and evolve over time. For M* > 10^8 MSun, LSBGs contribute 47, 7 and 6 per cent of the local number, mass and luminosity densities respectively (~85/11/10 per cent for M* > 10^7 MSun). Today's LSBGs have similar dark-matter fractions and angular momenta to high-surface-brightness galaxies (HSBGs; mu < 23 mag arcsec^-2), but larger effective radii (x2.5 for UDGs) and lower fractions of dense, star-forming gas (more than x6 less in UDGs than HSBGs). LSBGs originate from the same progenitors as HSBGs at z > 2. However, LSBG progenitors form stars more rapidly at early epochs. The higher resultant rate of supernova-energy injection flattens their gas-density profiles, which, in turn, creates shallower stellar profiles that are more susceptible to tidal processes. After z ~ 1, tidal perturbations broaden LSBG stellar distributions and heat their cold gas, creating the diffuse, largely gas-poor LSBGs seen today. In clusters, ram-pressure stripping provides an additional mechanism that assists in gas removal in LSBG progenitors. Our results offer insights into the formation of a galaxy population that is central to a complete understanding of galaxy evolution, and which will be a key topic of research using new and forthcoming deep-wide surveys.

Context. Accurate model predictions including the physics of baryons are required to make the most of the upcoming large cosmological surveys devoted to gravitational lensing. The advent of hydrodynamical cosmological simulations enables such predictions on sufficiently sizeable volumes. Aims. Lensing quantities (deflection, shear, convergence) and their statistics (convergence power spectrum, shear correlation functions, galaxy-galaxy lensing) are computed in the past lightcone built in the Horizon-AGN hydrodynamical cosmological simulation, which implements our best knowledge on baryonic physics at the galaxy scale in order to mimic galaxy populations over cosmic time. Methods. Lensing quantities are generated over a one square degree field of view by performing multiple-lens plane ray-tracing through the lightcone, taking full advantage of the 1 kpc resolution and splitting the line of sight over 500 planes all the way to redshift z~7. Two methods are explored (standard projection of particles with adaptive smoothing, and integration of the acceleration field) to assert a good implementation. The focus is on small scales where baryons matter most. Results. Standard cosmic shear statistics are impacted at the 10% level by the baryonic component for angular scales below a few arcmin. The galaxy-galaxy lensing signal, or galaxy-shear correlation function, is consistent with measurements for the redshift z~0.5 massive galaxy population. At higher redshift z>1, the impact of magnification bias on this correlation is relevant for separations greater than 1 Mpc. Conclusions. This work is pivotal for all current and upcoming weak lensing surveys and represents a first step towards building a full end-to-end generation of lensed mock images from large cosmological hydrodynamical simulations.

Using the two large cosmological hydrodynamical simulations, Horizon-AGN (H-AGN) and Horizon-noAGN (H-noAGN, no AGN feedback), we investigate how a typical sub-grid model for AGN feedback affects the evolution of the total density profiles (dark matter + stars) at the effective radius of massive early-type galaxies (M*>10^11 Msun). We have studied the dependencies of the mass-weighted density slope gamma'_tot with the effective radius, the galaxy mass and the host halo mass at z~0.3 and found that the inclusion of AGN feedbackalways leads to a much better agreement with observational values and trends. Our analysis suggests also that the inclusion of AGN feedback favours a strong correlation between gamma'_tot and the density slope of the dark matter component while, in the absence of AGN activity, gamma'_tot is rather strongly correlated with the density slope of the stellar component. Finally, we find that gamma'_tot derived from our samples of galaxies increases from z=2 to z=0,in good agreement with the expected observational trend. The derived slopes are slightly lower than in the data when AGN is included because the simulated galaxies tend to be too extended, especially the least massive ones. However, the simulated compact galaxies without AGN feedback have gamma'_tot values that are significantly too high compared to observations.

Accurate cosmology from upcoming weak lensing surveys relies on knowledge of the total matter power spectrum at percent level at scales k &lt; 10 $h$/Mpc, for which modelling the impact of baryonic physics is crucial. We compare measurements of the total matter power spectrum from the Horizon cosmological hydrodynamical simulations: a dark matter-only run, one with full baryonic physics, and another lacking Active Galactic Nuclei (AGN) feedback. Baryons cause a suppression of power at $k\simeq 10$ $h/$Mpc of &lt;15\% at $z=0$, and an enhancement of a factor of a few at smaller scales due to the more efficient cooling and star formation. The results are sensitive to the presence of the highest mass haloes in the simulation and the distribution of dark matter is also impacted up to a few percent. The redshift evolution of the effect is non-monotonic throughout $z=0-5$ due to an interplay between AGN feedback and gas pressure, and the growth of structure. We investigate the effectiveness of an analytic `baryonic correction model' in describing our results. We require a different redshift evolution and propose an alternative fitting function with $4$ free parameters that reproduces our results within $5\%$. Compared to other simulations, we find the impact of baryonic processes on the total matter power spectrum to be smaller at $z=0$. Correspondingly, our results suggest that AGN feedback is not strong enough in the simulation. Total matter power spectra from the Horizon simulations are made publicly available at this https URL

As endpoints of the hierarchical mass-assembly process, the stellar populations of local early-type galaxies encode the assembly history of galaxies over cosmic time. We use Horizon-AGN, a cosmological hydrodynamical simulation, to study the merger histories of local early-type galaxies and track how the morphological mix of their progenitors evolves over time. We provide a framework for alleviating `progenitor bias' -- the bias that occurs if one uses only early-type galaxies to study the progenitor population. Early-types attain their final morphology at relatively early epochs -- by $z\sim1$, around 60 per cent of today's early-types have had their last significant merger. At all redshifts, the majority of mergers have one late-type progenitor, with late-late mergers dominating at z&gt;1.5 and early-early mergers becoming significant only at z&lt;0.5. Progenitor bias is severe at all but the lowest redshifts -- e.g. at $z\sim0.6$, less than 50 per cent of the stellar mass in today's early-types is actually in progenitors with early-type morphology, while, at $z\sim2$, studying only early-types misses almost all (80 per cent) of the stellar mass that eventually ends up in local early-type systems. At high redshift, almost all massive late-type galaxies, regardless of their local environment or star-formation rate, are progenitors of local early-type galaxies, as are lower-mass (M$_\star$ &lt; 10$^{10.5}$ M$_{\odot}$) late-types as long as they reside in high density environments. In this new era of large observational surveys (e.g. LSST, JWST), this study provides a framework for studying how today's early-type galaxies have been built up over cosmic time.

Using a suite of three large cosmological hydrodynamical simulations, Horizon-AGN, Horizon-noAGN (no AGN feedback) and Horizon-DM (no baryons), we investigate how a typical sub-grid model for AGN feedback affects the evolution of the inner density profiles of massive dark matter haloes and galaxies. Based on direct object-to-object comparisons, we find that the integrated inner mass and density slope differences between objects formed in these three simulations (hereafter, H_AGN, H_noAGN and H_DM) significantly evolve with time. More specifically, at high redshift (z~5), the mean central density profiles of H_AGN and H_noAGN dark matter haloes tend to be much steeper than their H_DM counterparts owing to the rapidly growing baryonic component and ensuing adiabatic contraction. By z~1.5, these mean halo density profiles in H_AGN have flattened, pummelled by powerful AGN activity ("quasar mode"): the integrated inner mass difference gaps with H_noAGN haloes have widened, and those with H_DM haloes have narrowed. Fast forward 9.5 billion years, down to z=0, and the trend reverses: H_AGN halo mean density profiles drift back to a more cusped shape as AGN feedback efficiency dwindles ("radio mode"), and the gaps in integrated central mass difference with H_noAGN and H_DM close and broaden respectively. On the galaxy side, the story differs noticeably. Averaged stellar profile central densities and inner slopes are monotonically reduced by AGN activity as a function of cosmic time, resulting in better agreement with local observations. As both dark matter and stellar inner density profiles respond quite sensitively to the presence of a central AGN, there is hope that future observational determinations of these quantities can be used constrain AGN feedback models.

We investigate the rate of orbital orientation dilution of young stellar clusters in the vicinity of supermassive black holes. Within the framework of vector resonant relaxation, we predict the time evolution of the two-point correlation function of the stellar orbital plane orientations as a function of their initial angular separation and diversity in orbital parameters (semi-major axis, eccentricity). As expected, the larger the spread in initial orientations and orbital parameters, the more efficient the dilution of a given set of co-eval stars, with a characteristic timescale set up by the coherence time of the background potential fluctuations. A Markovian prescription which matches numerical simulations allows us to efficiently probe the underlying kinematic properties of the unresolved nucleus when requesting consistency with a given dilution efficiency, imposed by the observed stellar disc within the one arcsecond of Sgr A*. As a proof of concept, we compute maps of constant dilution times as a function of the semi major axis cusp index and fraction of intermediate mass black holes in the old background stellar cluster. This computation suggests that vector resonant relaxation should prove useful in this context since it impacts orientations on timescales comparable to the stars' age.

We present the first quantitative detection of large-scale filamentary structure at $z \simeq 0.7$ in the large cosmological volume probed by the VIMOS Public Extragalactic Redshift Survey (VIPERS). We use simulations to show the capability of VIPERS to recover robust topological features in the galaxy distribution, in particular the filamentary network. We then investigate how galaxies with different stellar masses and stellar activities are distributed around the filaments and find a significant segregation, with the most massive or quiescent galaxies being closer to the filament axis than less massive or active galaxies. The signal persists even after down-weighting the contribution of peak regions. Our results suggest that massive and quiescent galaxies assemble their stellar mass through successive mergers during their migration along filaments towards the nodes of the cosmic web. On the other hand, low-mass star-forming galaxies prefer the outer edge of filaments, a vorticity rich region dominated by smooth accretion, as predicted by the recent spin alignment theory. This emphasizes the role of large scale cosmic flows in shaping galaxy properties.

Context. Upcoming weak lensing surveys such as Euclid will provide an unprecedented opportunity to quantify the geometry and topology of the cosmic web, in particular in the vicinity of lensing clusters. Aims. Understanding the connectivity of the cosmic web with unbiased mass tracers, such as weak lensing, is of prime importance to probe the underlying cosmology, seek dynamical signatures of dark matter, and quantify environmental effects on galaxy formation. Methods. Mock catalogues of galaxy clusters are extracted from the N-body PLUS simulation. For each cluster, the aperture multipolar moments of the convergence are calculated in two annuli (inside and outside the virial radius). By stacking their modulus, a statistical estimator is built to characterise the angular mass distribution around clusters. The moments are compared to predictions from perturbation theory and spherical collapse. Results. The main weakly chromatic excess of multipolar power on large scales is understood as arising from the contraction of the primordial cosmic web driven by the growing potential well of the cluster. Besides this boost, the quadrupole prevails in the cluster (ellipsoidal) core, while at the outskirts, harmonic distortions are spread on small angular modes, and trace the non-linear sharpening of the filamentary structures. Predictions for the signal amplitude as a function of the cluster-centric distance, mass, and redshift are presented. The prospects of measuring this signal are estimated for current and future lensing data sets. Conclusions. The Euclid mission should provide all the necessary information for studying the cosmic evolution of the connectivity of the cosmic web around lensing clusters using multipolar moments and probing unique signatures of, for example, baryons and warm dark matter.