Jinsu Rhee, Sukyoung K. Yi, Jongwan Ko, Emanuele Contini, J. K. Jang, Seyoung Jeon, San Han, Christophe Pichon, Yohan Dubois, Katarina Kraljic, Sébastien Peirani We study star formation (SF) quenching of satellite galaxies with $M_{*} > 10^7\,M_{\odot}$ within two low-mass groups ($M_{\rm vir}=10^{12.9}$ and $10^{12.7} \,M_{\odot}$) using the NewHorizon simulation. We confirm that satellite galaxies ($M_{*}\lesssim10^{10}\,M_{\odot}$) are more prone to quenching than their field counterparts. This quenched fraction decreases with increasing stellar mass, consistent with recent studies. Similar to the findings in cluster environments, we note a correlation between the orbital motions of galaxies within these groups and the phenomenon of SF quenching. Specifically, SF is suppressed at the group center, and for galaxies with $M_{*} > 10^{9.1}\,M_{\odot}$, there is often a notable rejuvenation phase following a temporary quenching period. The SF quenching at the group center is primarily driven by changes in star formation efficiency and the amount of gas available, both of which are influenced by hydrodynamic interactions between the interstellar medium and surrounding hot gas within the group. Conversely, satellite galaxies with $M_{*} < 10^{8.2}\,M_{\odot}$ experience significant gas removal within the group, leading to SF quenching. Our analysis highlights the complexity of SF quenching in satellite galaxies in group environments, which involves an intricate competition between the efficiency of star formation (which depends on the dynamical state of the gas) on the one hand, and the availability of cold dense gas on the other hand. This challenges the typical understanding of environmental effects based on gas stripping through ram pressure, suggesting a need for a new description of galaxy evolution under mild environmental effects.
Sang Hyeok Im, Ho Seong Hwang, Jaehong Park, Jaehyun Lee, Hyunmi Song, Stephen Appleby, Yohan Dubois, C. Gareth Few, Brad K. Gibson, Juhan Kim, Yonghwi Kim, Changbom Park, Christophe Pichon, Jihye Shin, Owain N. Snaith, Maria Celeste Artale, Eric Gawiser, Lucia Guaita, Woong-Seob Jeong, Kyoung-Soo Lee, et al (4) We test whether Lyman alpha emitters (LAEs) and Lyman-break galaxies (LBGs) can be good tracers of high-z large-scale structures, using the Horizon Run 5 cosmological hydrodynamical simulation. We identify LAEs using the Ly\alpha emission line luminosity and its equivalent width, and LBGs using the broad-band magnitudes at z~2.4, 3.1, and 4.5. We first compare the spatial distributions of LAEs, LBGs, all galaxies, and dark matter around the filamentary structures defined by dark matter. The comparison shows that both LAEs and LBGs are more concentrated toward the dark matter filaments than dark matter. We also find an empirical fitting formula for the vertical density profile of filaments as a binomial power-law relation of the distance to the filaments. We then compare the spatial distributions of the samples around the filaments defined by themselves. LAEs and LBGs are again more concentrated toward their filaments than dark matter. We also find the overall consistency between filamentary structures defined by LAEs, LBGs, and dark matter, with the median spatial offsets that are smaller than the mean separation of the sample. These results support the idea that the LAEs and LBGs could be good tracers of large-scale structures of dark matter at high redshifts.
Feedback from stars and active galactic nuclei (AGNs) primarily affects the formation and evolution of galaxies and the circumgalactic medium, leaving some kind of imprint on larger scales. Based on the \sc Simba hydrodynamical simulation suite and using the full set of Minkowski functionals (MFs), this study systematically analyses the time evolution of the global geometry and topology of the gas temperature, pressure, density (total, HI, and H$_2$), and the metallicity fields between redshifts $z=5$ and $z=0$. The MFs show that small-scale astrophysical processes are persistent and manifest on larger, up to tens of Mpc scales, highlighting the specific morphological signatures of the relevant feedback mechanisms on these scales in the last $\sim12$~Gyr. In qualitative terms, we were able establish a ranking that varies according to the field considered: stellar feedback mostly determines the morphology of the pressure and density fields and AGN jets are the primary origin of the morphology of the temperature and metallicity fields, while X-ray heating and AGN winds play the second most important role in shaping the geometry and topology of all the gaseous fields, except metallicity. Hence, the cosmic evolution of the geometry and topology of fields characterising the thermodynamical and chemical properties of the cosmic web offers complementary, larger scale constraints to galaxy formation models.
The formation pathways of compact stellar systems (CSSs) are still under debate. We utilize the \NH simulation to investigate the origins of such objects in the field environment. We identified 55 CSS candidates in the simulation whose properties are similar to those of the observed ultra-compact dwarfs and compact ellipticals. All but two most massive objects (compact elliptical candidates) are a result of a short starburst. Sixteen are formed by tidal stripping, while the other 39 are intrinsically compact from their birth. The stripped objects originate from dwarf-like galaxies with a dark halo, but most of their dark matter is stripped through their orbital motion around a more massive neighbor galaxy. The 39 intrinsically compact systems are further divided into ``associated'' or ``isolated'' groups, depending on whether they were born near a massive dark halo or not. The isolated intrinsic compact objects (7) are born in a dark halo and their stellar properties are older and metal-poor compared to the associated counterparts (32). The stripped compact objects occupy a distinct region in the age-metallicity plane from the intrinsic compact objects. The associated intrinsic compact objects in our sample have never had a dark halo; they are the surviving star clumps of a massive galaxy.
Aaron R. Rowntree, Ankit Singh, Fiorenzo Vincenzo, Brad K. Gibson, Céline Gouin, Daniela Galárraga-Espinosa, Jaehyun Lee, Juhan Kim, Clotilde Laigle, Changbom Park, Christophe Pichon, Gareth Few, Sungwook E. Hong, Yonghwi Kim Metallicity offers a unique window into the baryonic history of the cosmos, being instrumental in probing evolutionary processes in galaxies between different cosmic environments. We aim to quantify the contribution of these environments to the scatter in the mass-metallicity relation (MZR) of galaxies. By analysing the galaxy distribution within the cosmic skeleton of the Horizon Run 5 cosmological hydrodynamical simulation at redshift $z = 0.625$, computed using a careful calibration of the T-ReX filament finder, we identify galaxies within three main environments: nodes, filaments and voids. We also classify galaxies based on the dynamical state of the clusters and the length of the filaments in which they reside. We find that the cosmic environment significantly contributes to the scatter in the MZR; in particular, both the gas metallicity and its average relative standard deviation increase when considering denser large-scale environments. The difference in the average metallicity between galaxies within relaxed and unrelaxed clusters is $\approx 0.1 \text{ dex}$, with both populations displaying positive residuals, $\delta Z_{g}$, from the averaged MZR. Moreover, the difference in metallicity between node and void galaxies accounts for $\approx 0.14 \, \text{dex}$ in the scatter of the MZR at stellar mass $M_{\star} \approx 10^{9.35}\,\text{M}_{\odot}$. Finally, both the average [O/Fe] in the gas and the galaxy gas fraction decrease when moving to higher large-scale densities in the simulation, suggesting that the cores of cosmic environments host, on average, older and more massive galaxies, whose enrichment is affected by a larger number of Type Ia Supernova events.
Jaewon Yoo, Changbom Park, Cristiano G. Sabiu, Ankit Singh, Jongwan Ko, Jaehyun Lee, Christophe Pichon, M. James Jee, Brad K. Gibson, Owain Snaith, Juhan Kim, Jihye Shin, Yonghwi Kim, Hyowon Kim One intriguing approach for studying the dynamical evolution of galaxy clusters is to compare the spatial distributions among various components, such as dark matter, member galaxies, gas, and intracluster light (ICL). Utilizing the recently introduced Weighted Overlap Coefficient (WOC) \citep2022ApJS..261...28Y, we analyze the spatial distributions of components within 174 galaxy clusters ($M_{\rm tot}> 5 \times 10^{13} M_{\odot}$, $z=0.625$) at varying dynamical states in the cosmological hydrodynamical simulation Horizon Run 5. We observe that the distributions of gas and the combination of ICL with the brightest cluster galaxy (BCG) closely resembles the dark matter distribution, particularly in more relaxed clusters, characterized by the half-mass epoch. The similarity in spatial distribution between dark matter and BCG+ICL mimics the changes in the dynamical state of clusters during a major merger. Notably, at redshifts $>$ 1, BCG+ICL traced dark matter more accurately than the gas. Additionally, we examined the one-dimensional radial profiles of each component, which show that the BCG+ICL is a sensitive component revealing the dynamical state of clusters. We propose a new method that can approximately recover the dark matter profile by scaling the BCG+ICL radial profile. Furthermore, we find a recipe for tracing dark matter in unrelaxed clusters by including the most massive satellite galaxies together with BCG+ICL distribution. Combining the BCG+ICL and the gas distribution enhances the dark matter tracing ability. Our results imply that the BCG+ICL distribution is an effective tracer for the dark matter distribution, and the similarity of spatial distribution may be a useful probe of the dynamical state of a cluster.
Rayleigh-Levy flights have played a significant role in cosmology as simplified models for understanding how matter distributes itself under gravitational influence. These models also exhibit numerous remarkable properties that enable the prediction of a wide range of characteristics. Here, we derive the one and two point statistics of extreme points within Rayleigh-Levy flights spanning one to three dimensions, stemming directly from fundamental principles. In the context of the mean field limit, we provide straightforward closed-form expressions for Euler counts and their correlations, particularly in relation to their clustering behaviour over long distances. Additionally, quadratures allow for the computation of extreme value number densities. A comparison between theoretical predictions in 1D and Monte Carlo measurements shows remarkable agreement. Given the widespread use of Rayleigh-Levy processes, these comprehensive findings offer significant promise not only in astrophysics but also in broader applications beyond the field.
The long-term relaxation of rotating, spherically symmetric globular clusters is investigated through an extension of the orbit averaged Chandrasekhar non-resonant formalism. A comparison is made with the long-term evolution of the distribution function in action space, measured from averages of sets of $N$-body simulations up to core collapse. The impact of rotation on in-plane relaxation is found to be weak. In addition, we observe a clear match between theoretical predictions and $N$-body measurements. For the class of rotating models considered, we find no strong gravo-gyro catastrophe accelerating core collapse. Both kinetic theory and simulations predict a reshuffling of orbital inclinations from overpopulated regions to underpopulated ones. This trend accelerates as the amount of rotation is increased. Yet, for orbits closer to the rotational plane, the non-resonant prediction does not reproduce numerical measurements. We argue that this mismatch stems from these orbits' coherent interactions, which are not captured by the non-resonant formalism that only addresses local deflections.
Sebastien Peirani, Yasushi Suto, Ricarda S. Beckmann, Marta Volonteri, Yen-Ting Lin, Yohan Dubois, Sukyoung K. Yi, Christophe Pichon, Katarina Kraljic, Minjung Park, Julien Devriendt, San Han, Wei-Huai Chen (Reduced) Using the recent cosmological high-resolution zoom-in simulations, NewHorizon and Galactica, in which the evolution of black hole spin is followed on the fly, we have tracked the cosmic history of a hundred of black holes (BHs) with a mass greater than 2x10^4 Ms. For each of them, we have studied the variations of the three dimensional angle (Psi) subtended between the BH spins and the angular momentum vectors of their host galaxies. The analysis of the individual evolution of the most massive BHs suggests that they are generally passing by three different regimes. First, for a short period after their birth, low mass BHs (<3x10^4 Ms) are rapidly spun up by gas accretion and their spin tends to be aligned with their host galaxy spin. Then follows a second phase in which the accretion of gas onto low mass BHs (<10^5 Ms) is quite chaotic and inefficient, reflecting the complex and disturbed morphologies of forming proto-galaxies at high redshifts. The variations of Psi are rather erratic during this phase and are mainly driven by the rapid changes of the direction of the galaxy angular momentum. Then, in a third and long phase, BHs are generally well settled in the center of galaxies around which the gas accretion becomes much more coherent (>10^5 Ms). In this case, the BH spins tend to be well aligned with the angular momentum of their host galaxy and this configuration is generally stable even though BH merger episodes can temporally induce misalignment. We have also derived the distributions of cos(Psi) at different redshifts and found that BHs and galaxy spins are generally aligned. Finally, based on a Monte Carlo method, we also predict statistics for the 2-d projected spin-orbit angles lambda. In particular, the distribution of lambda traces well the alignment tendency in the 3-d analysis. Such predictions provide an interesting background for future observational analyses.
We present LinearResponse.jl, an efficient, versatile public library written in julia to compute the linear response of self-gravitating (3D spherically symmetric) stellar spheres and (2D axisymmetric razor-thin) discs. LinearResponse.jl can scan the whole complex frequency plane, probing unstable, neutral and (weakly) damped modes. Given a potential model and a distribution function, this numerical toolbox estimates the modal frequencies as well as the shapes of individual modes. The libraries are validated against a combination of previous results for the spherical isochrone model and Mestel discs, and new simulations for the spherical Plummer model. Beyond linear response theory, the realm of applications of LinearResponse.jl also extends to the kinetic theory of self-gravitating systems through a modular interface.
Exclusion zones in the cross-correlations between critical points (peak-void, peak-wall, filament-wall, filament-void) of the density field define quasi-standard rulers that can be used to constrain dark matter and dark energy cosmological parameters. The average size of the exclusion zone is found to scale linearly with the typical distance between extrema. The latter changes as a function of the matter content of the universe in a predictable manner, but its comoving size remains essentially constant in the linear regime of structure growth on large scales, unless the incorrect cosmology is assumed in the redshift-distance relation. This can be used to constrain the dark energy parameters when considering a survey that scans a range of redshifts. The precision of the parameter estimation is assessed using a set of cosmological simulations, and is found to be a 4$\sigma$ detection of a change in matter content of 5%, or about 3.8$\sigma$ detection of 50% shift in the dark energy parameter using a full sky survey up to redshift 0.5.
San Han, Sukyoung K. Yi, Sree Oh, Mina Pak, Scott M. Croom, Julien Devriendt, Yohan Dubois, Taysun Kimm, Katarina Kraljic, Christophe Pichon, Marta Volonteri Observed and simulated galaxies exhibit a significant variation in their velocity dispersion profiles. We examine the inner and outer slopes of stellar velocity dispersion profiles using integral field spectroscopy data from two surveys, SAMI (for $z < 0.115$) and CALIFA (for $z < 0.03$), comparing them with results from two cosmological hydrodynamic simulations: Horizon-AGN (for $z = 0.017$) and NewHorizon (for $z\lesssim1$). The simulated galaxies closely reproduce the variety of velocity dispersion slopes and stellar mass dependence of both inner and outer radii ($0.5\,r_{50}$ and $3\,r_{50}$) as observed, where $r_{50}$ stands for half-light radius. The inner slopes are mainly influenced by the relative radial distribution of the young and old stars formed in-situ: a younger center shows a flatter inner profile. The presence of accreted (ex-situ) stars has two effects on the velocity dispersion profiles. First, because they are more dispersed in spatial and velocity distributions compared to in-situ formed stars, it increases the outer slope of the velocity dispersion profile. It also causes the velocity anisotropy to be more radial. More massive galaxies have a higher fraction of stars formed ex-situ and hence show a higher slope in outer velocity dispersion profile and a higher degree of radial anisotropy. The diversity in the outer velocity dispersion profiles reflects the diverse assembly histories among galaxies.
R. A. Jackson, S. Kaviraj, S. K. Yi, S. Peirani, Y. Dubois, G. Martin, J. E. G. Devriendt, A. Slyz, C. Pichon, M. Volonteri, T. Kimm, K. Kraljic The `core-cusp' problem is considered a key challenge to the LCDM paradigm. Halos in dark matter only simulations exhibit `cuspy' profiles, where density continuously increases towards the centre. However, the dark matter profiles of many observed galaxies (particularly in the dwarf regime) deviate strongly from this prediction, with much flatter central regions (`cores'). We use NewHorizon (NH), a hydrodynamical cosmological simulation, to investigate core formation, using a statistically significant number of galaxies in a cosmological volume. Halos containing galaxies in the upper (M* > 10^10.2 MSun) and lower (M* < 10^8 MSun) ends of the stellar mass distribution contain cusps. However, halos containing galaxies with intermediate (10^8 MSun < M* < 10^10.2 MSun) stellar masses are generally cored, with typical halo masses between 10^10.2 MSun and 10^11.5 MSun. Cores form through supernova-driven gas removal from halo centres, which alters the central gravitational potential, inducing dark matter to migrate to larger radii. While all massive (M* > 10^9.5 MSun) galaxies undergo a cored-phase, in some cases cores can be removed and cusps reformed. This happens if a galaxy undergoes sustained star formation at high redshift, which results in stars (which, unlike the gas, cannot be removed by baryonic feedback) dominating the central gravitational potential. After cosmic star formation peaks, the number of cores, and the mass of the halos they are formed in, remain constant, indicating that cores are being routinely formed over cosmic time after a threshold halo mass is reached. The existence of cores is, therefore, not in tension with the standard paradigm.
We build a model to predict from first principles the properties of major mergers. We predict these from the coalescence of peaks and saddle points in the vicinity of a given larger peak, as one increases the smoothing scale in the initial linear density field as a proxy for cosmic time. To refine our results, we also ensure, using a suite of $\sim 400$ power-law Gaussian random fields smoothed at $\sim 30$ different scales, that the relevant peaks and saddles are topologically connected: they should belong to a persistent pair before coalescence. Our model allows us to (a) compute the probability distribution function of the satellite-merger separation in Lagrangian space: they peak at three times the smoothing scale; (b) predict the distribution of the number of mergers as a function of peak rarity: haloes typically undergo two major mergers ($>$1:10) per decade of mass growth; (c) recover that the typical spin brought by mergers: it is of the order of a few tens of percent.
Katarina Kraljic, Florent Renaud, Yohan Dubois, Christophe Pichon, Oscar Agertz, Eric Andersson, Julien Devriendt, Jonathan Freundlich, Sugata Kaviraj, Taysun Kimm, Garreth Martin, Sébastien Peirani, Álvaro Segovia Otero, Marta Volonteri, Sukyoung K. Yi The scaling relations between the gas content and star formation rate of galaxies provide useful insights into processes governing their formation and evolution. We investigate the emergence and the physical drivers of the global Kennicutt-Schmidt (KS) relation at $0.25 \leq z \leq 4$ in the cosmological hydrodynamic simulation NewHorizon capturing the evolution of a few hundred galaxies with a resolution of $\sim$ 40 pc. The details of this relation vary strongly with the stellar mass of galaxies and the redshift. A power-law relation $\Sigma_{\rm SFR} \propto \Sigma_{\rm gas}^{a}$ with $a \approx 1.4$, like that found empirically, emerges at $z \approx 2 - 3$ for the most massive half of the galaxy population. However, no such convergence is found in the lower-mass galaxies, for which the relation gets shallower with decreasing redshift. At the galactic scale, the star formation activity correlates with the level of turbulence of the interstellar medium, quantified by the Mach number, rather than with the gas fraction (neutral or molecular), confirming previous works. With decreasing redshift, the number of outliers with short depletion times diminishes, reducing the scatter of the KS relation, while the overall population of galaxies shifts toward low densities. Using pc-scale star formation models calibrated with local Universe physics, our results demonstrate that the cosmological evolution of the environmental and intrinsic conditions conspire to converge towards a significant and detectable imprint in galactic-scale observables, in their scaling relations, and in their reduced scatter.
Sukyoung K. Yi, J. K. Jang, Julien Devriendt, Yohan Dubois, San Han, Taysun Kimm, Katarina Kraljic, Minjung Park, Sebastien Peirani, Christophe Pichon, Jinsu Rhee Thick disks are a prevalent feature observed in numerous disk galaxies including our own Milky Way. Their significance has been reported to vary widely, ranging from a few to 100% of the disk mass, depending on the galaxy and the measurement method. We use the NewHorizon simulation which has high spatial and stellar mass resolutions to investigate the issue of thick disk mass fraction. We also use the NewHorizon2 simulation that was run on the same initial conditions but additionally traced nine chemical elements. Based on a sample of 27 massive disk galaxies with M* > 10^10 M_⊙ in NewHorizon, the contribution of the thick disk was found to be 34 \pm 15% in r-band luminosity or 48 \pm 13% in mass to the overall galactic disk, which seems in agreement with observational data. The vertical profiles of 0, 22, and 5 galaxies are best fitted by 1, 2, or 3 sech2 components, respectively. The NewHorizon2 data show that the selection of thick disk stars based on a single [\alpha/Fe] cut is severely contaminated by stars of different kinematic properties while missing a bulk of kinematically thick disk stars. Vertical luminosity profile fits recover the key properties of thick disks reasonably well. The majority of stars are born near the galactic mid-plane with high circularity and get heated with time via fluctuation in the force field. Depending on the star formation and merger histories, galaxies may naturally develop thick disks with significantly different properties.
We propose a new method for finding galaxy protoclusters that is motivated by structure formation theory and also directly applicable to observations. We adopt the conventional definition that a protocluster is a galaxy group whose virial mass $M_{\rm vir} < M_{\rm cl}$ at its epoch, where $M_{\rm cl}=10^{14}\,M_{\odot}$, but would exceed that limit when it evolves to $z=0$. We use the critical overdensity for complete collapse at $z = 0$ predicted by the spherical top-hat collapse model to find the radius and total mass of the regions that would collapse at $z=0$. If the mass of a region centered at a massive galaxy exceeds $M_{\rm cl}$, the galaxy is at the center of a protocluster. We define the outer boundary of protocluster as the zero-velocity surface at the turnaround radius so that the member galaxies are those sharing the same protocluster environment and showing some conformity in physical properties. We use the cosmological hydrodynamical simulation Horizon Run 5 (HR5) to calibrate this prescription and demonstrate its performance. We find that the protocluster identification method suggested in this study is quite successful. Its application to the high-redshift HR5 galaxies shows a tight correlation between the mass within the protocluster regions identified according to the spherical collapse model and the final mass to be found within the clusters at $z=0$, meaning that the regions can be regarded as the bona fide protoclusters with high reliability. We also confirm that the redshift-space distortion does not significantly affect the performance of the protocluster identification scheme.
Juhan Kim, Jaehyun Lee, Clotilde Laigle, Yohan Dubois, Yonghwi Kim, Changbom Park, Christophe Pichon, Brad Gibson, C. Gareth Few, Jihye Shin, Owain Snaith We investigate the impact of the surface brightness (SB) limit on the galaxy stellar mass functions (GSMFs) using mock surveys generated from the Horizon Run 5 (HR5) simulation. We compare the stellar-to-halo-mass relation, GSMF, and size-stellar mass relation of the HR5 galaxies with empirical data and other cosmological simulations. The mean SB of simulated galaxies are computed using their effective radii, luminosities, and colors. To examine the cosmic SB dimming effect, we compute $k$-corrections from the spectral energy distributions of individual simulated galaxy at each redshift, apply the $k$-corrections to the galaxies, and conduct mock surveys based on the various SB limits. We find that the GSMFs are significantly affected by the SB limits at a low-mass end. This approach can ease the discrepancy between the GSMFs obtained from simulations and observations at $0.625\le z\le 2$. We also find that a redshift survey with a SB selection limit of $\left<\mu_r\right>^e =$ 28 mag arcsec${}^{-2}$ will miss 20% of galaxies with $M_\star^g=10^{9}~{\rm M_\odot}$ at $z=0.625$. The missing fraction of low-surface-brightness galaxies increases to 50%, 70%, and 98% at $z=0.9$, 1.1, and 1.9, respectively, at the SB limit.
Here we use the Horizon-AGN simulation to test whether the spins of SMBHs in merger-free galaxies are higher. We select samples using an observationally motivated bulge-to-total mass ratio of < 0.1, along with two simulation motivated thresholds selecting galaxies which have not undergone a galaxy merger since z = 2, and those SMBHs with < 10% of their mass due to SMBH mergers. We find higher spins (> 5\sigma ) in all three samples compared to the rest of the population. In addition, we find that SMBHs with their growth dominated by BH mergers following galaxy mergers, are less likely to be aligned with their galaxy spin than those that have grown through accretion in the absence of galaxy mergers (3.4\sigma ). We discuss the implications this has for the impact of active galactic nuclei (AGN) feedback, finding that merger-free SMBHs spend on average 91% of their lifetimes since z = 2 in a radio mode of feedback (88% for merger-dominated galaxies). Given that previous observational and theoretical works have concluded that merger-free processes dominate SMBH-galaxy co-evolution, our results suggest that this co-evolution could be regulated by radio mode AGN feedback.
J. K. Jang, Sukyoug K. Yi, Yohan Dubois, Jinsu Rhee, Christophe Pichon, Taysun Kimm, Julien Devriendt, Marta Volonteri, Sugata Kaviraj, Sebastien Peirani, Sree Oh, Scott Croom Based on the recent advancements in the numerical simulations of galaxy formation, we anticipate the achievement of realistic models of galaxies in the near future. Morphology is the most basic and fundamental property of galaxies, yet observations and simulations still use different methods to determine galaxy morphology, making it difficult to compare them. We hereby perform a test on the recent NewHorizon simulation which has spatial and mass resolutions that are remarkably high for a large-volume simulation, to resolve the situation. We generate mock images for the simulated galaxies using SKIRT that calculates complex radiative transfer processes in each galaxy. We measure morphological indicators using photometric and spectroscopic methods following observer's techniques. We also measure the kinematic disk-to-total ratios using the Gaussian mixture model and assume that they represent the true structural composition of galaxies. We found that spectroscopic indicators such as $V/{\sigma}$ and ${\lambda}_{R}$ closely trace the kinematic disk-to-total ratios. In contrast, photometric disk-to-total ratios based on the radial profile fitting method often fail to recover the true kinematic structure of galaxies, especially for small galaxies. We provide translating equations between various morphological indicators.
Cluster galaxies exhibit substantially lower star formation rates than field galaxies today, but it is conceivable that clusters were sites of more active star formation in the early universe. Herein, we present an interpretation of the star formation history (SFH) of group/cluster galaxies based on the large-scale cosmological hydrodynamic simulation, Horizon-AGN. We find that massive galaxies in general have small values of e-folding timescales of star formation decay (i.e., ``mass quenching'') regardless of their environment, whilst low-mass galaxies exhibit prominent environmental dependence. In massive host halos (i.e., clusters), the e-folding timescales of low-mass galaxies are further decreased if they reside in such halos for a longer period of time. This ``environmental quenching'' trend is consistent with the theoretical expectation from ram pressure stripping. Furthermore, we define a ``transition epoch'' as where cluster galaxies become less star-forming than field galaxies. The transition epoch of group/cluster galaxies varies according to their stellar and host cluster halo masses. Low-mass galaxies in massive clusters show the earliest transition epoch of $\sim 7.6$ Gyr ago in lookback time. However, it decreases to $\sim 5.2$ Gyr for massive galaxies in low-mass clusters. Based on our findings, we can describe cluster galaxy's SFH with regard to the cluster halo-to-stellar mass ratio.
Residuals between measured galactic radii and those predicted by the Fundamental Plane (FP) are possible tracers of weak lensing magnification. However, observations have shown these to be systematically correlated with the large-scale structure. We use the Horizon-AGN hydrodynamical cosmological simulation to analyse these intrinsic size correlations (ISCs) for both elliptical (early-type) and spiral (late-type) galaxies at $z=0.06$. We fit separate FPs to each sample, finding similarly distributed radius residuals, $\lambda$, in each case. We find persistent $\lambda\lambda$ correlations over three-dimensional separations $0.5-17\,h^{-1}{\rm{Mpc}}$ in the case of spiral galaxies, at $>3\sigma$ significance. When relaxing a mass-selection, applied for better agreement with galaxy clustering constraints, the spiral $\lambda\lambda$ detection strengthens to $9\sigma$; we detect a $5\sigma$ density-$\lambda$ correlation; and we observe intrinsically-large spirals to cluster more strongly than small spirals over scales $\lesssim10\,h^{-1}{\rm{Mpc}}$, at $>5\sigma$ significance. Conversely, and in agreement with the literature, we observe lower-mass, intrinsically-small ellipticals to cluster more strongly than their large counterparts over scales $0.5-17\,h^{-1}{\rm{Mpc}}$, at $>5\sigma$ significance. We model $\lambda\lambda$ correlations using a phenomenological non-linear size model, and predict the level of contamination for cosmic convergence analyses. We find the systematic contribution to be of similar order to, or dominant over the cosmological signal. We make a mock measurement of an intrinsic, systematic contribution to the projected surface mass density $\Sigma(r)$ and find statistically significant, low-amplitude, positive (negative) contributions from lower-mass spirals (ellipticals), which may be of concern for large-scale ($\gtrsim\,7\,h^{-1}$ Mpc) measurements.
We present LyMAS2, an improved version of the "Lyman-\alpha Mass Association Scheme" aiming at predicting the large-scale 3d clustering statistics of the Lyman-\alpha forest (Ly-\alpha) from moderate resolution simulations of the dark matter (DM) distribution, with prior calibrations from high resolution hydrodynamical simulations of smaller volumes. In this study, calibrations are derived from the Horizon-AGN suite simulations, (100 Mpc/h)^3 comoving volume, using Wiener filtering, combining information from dark matter density and velocity fields (i.e. velocity dispersion, vorticity, line of sight 1d-divergence and 3d-divergence). All new predictions have been done at z=2.5 in redshift-space, while considering the spectral resolution of the SDSS-III BOSS Survey and different dark matter smoothing (0.3, 0.5 and 1.0 Mpc/h comoving). We have tried different combinations of dark matter fields and found that LyMAS2, applied to the Horizon-noAGN dark matter fields, significantly improves the predictions of the Ly-\alpha 3d clustering statistics, especially when the DM overdensity is associated with the velocity dispersion or the vorticity fields. Compared to the hydrodynamical simulation trends, the 2-point correlation functions of pseudo-spectra generated with LyMAS2 can be recovered with relative differences of ~5% even for high angles, the flux 1d power spectrum (along the light of sight) with ~2% and the flux 1d probability distribution function exactly. Finally, we have produced several large mock BOSS spectra (1.0 and 1.5 Gpc/h) expected to lead to much more reliable and accurate theoretical predictions.
We investigate the long-term relaxation of one-dimensional (${1D}$) self-gravitating systems, using both kinetic theory and $N$-body simulations. We consider thermal and Plummer equilibria, with and without collective effects. All combinations are found to be in clear agreement with respect to the Balescu-Lenard and Landau predictions for the diffusion coefficients. Interestingly, collective effects reduce the diffusion by a factor ${\sim 10}$. The predicted flux for Plummer equilibrium matches the measured one, which is a remarkable validation of kinetic theory. We also report on a situation of quasi kinetic blocking for the same equilibrium.
G. Martin, A. E. Bazkiaei, M. Spavone, E. Iodice, J. C. Mihos, M. Montes, J. A. Benavides, S. Brough, J. L. Carlin, C. A. Collins, P. A. Duc, F. A. Gómez, G. Galaz, H. M. Hernández-Toledo, R. A. Jackson, S. Kaviraj, J. H. Knapen, C. Martínez-Lombilla, S. McGee, D. O'Ryan, et al (33) Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilising automated techniques and human visual classification in conjunction with realistic mock images produced using the NEWHORIZON cosmological simulation, we investigate the nature, frequency and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-year depth of the Legacy Survey of Space and Time (30-31 mag / sq. arcsec), falling to 60 per cent assuming a shallower final depth of 29.5 mag / sq. arcsec. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M*~10^11.5 Msun). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterisation of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimised, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z<0.2).
Changbom Park, Jaehyun Lee, Juhan Kim, Donghui Jeong, Christophe Pichon, Brad K. Gibson, Owain N. Snaith, Jihye Shin, Yonghwi Kim, Yohan Dubois, C. Gareth Few We investigate the formation and morphological evolution of the first galaxies in the cosmic morning ($10 \gtrsim z \gtrsim 4$) using the Horizon Run 5 (HR5) simulation. For galaxies above the stellar mass $M_{\star, {\rm min}} = 2\times 10^9\,M_{\odot}$, we classify them into disk, spheroid, and irregular types according to their asymmetry and stellar mass morphology. We find that about 2/3 of the galaxies have a Sérsic index $< 1.5$, reflecting the dominance of disk-type morphology in the cosmic morning. The rest are evenly distributed as incidental and transient irregulars or spheroids. These fractions are roughly independent of redshift and stellar mass up to $\sim10^{10}\,M_{\odot}$. Almost all the first galaxies with $M_{\star}> M_{\star, {\rm min}}$ at $z>4$ form at initial peaks of the matter density field. Large-scale structures in the universe emerge and grow like cosmic rhizomes as the underlying matter density fluctuations grow and form associations of galaxies in rare overdense regions and the realm of the galactic world is stretched into relatively lower-density regions along evolving filaments. The cosmic web of galaxies forms at lower redshifts when most rhizomes globally percolate. The primordial angular momentum produced by the induced tidal torques on protogalactic regions is correlated with the internal kinematics of galaxies and tightly aligned with the angular momentum of the total galaxy mass. The large-scale tidal field imprinted in the initial conditions seems responsible for the dominance of disk morphology, and for the tendency of galaxies to re-acquire a disk post-distortion.
We model the response of spherical, non-rotating Milky Way (MW) dark matter and stellar halos to the Large Magellanic Cloud (LMC) using the matrix method of linear response theory. Our computations reproduce the main features of the dark halo response from simulations. We show that these features can be well separated by a harmonic decomposition: the large scale over/underdensity in the halo (associated with its reflex motion) corresponds to the $\ell=1$ terms, and the local overdensity to the $\ell\geq2$ multipoles. Moreover, the dark halo response is largely dominated by the first order 'forcing' term, with little influence from self-gravity. This makes it difficult to constrain the underlying velocity distribution of the dark halo using the observed response of the stellar halo, but it allows us to investigate the response of stellar halo models with various velocity anisotropies: a tangential (respectively radial) halo produces a shallower (respectively stronger) response. We also show that only the local wake is responsible for these variations, the reflex motion being solely dependent on the MW potential. Therefore, we identify the structure (orientation and winding) of the in-plane quadrupolar ($m=2$) response as a potentially good probe of the stellar halo anisotropy. Finally, our method allows us to tentatively relate the wake strength and shape to resonant effects: the strong radial response could be associated with the inner Lindblad resonance, and the weak tangential one with corotation.
Globular clusters are dense stellar systems whose core slowly contracts under the effect of self-gravity. The rate of this process was recently found to be directly linked to the initial amount of velocity anisotropy: tangentially anisotropic clusters contract faster than radially anisotropic ones. Furthermore, initially anisotropic clusters are found to generically tend towards more isotropic distributions during the onset of contraction. Chandrasekhar's "non-resonant" (NR) theory of diffusion describes this relaxation as being driven by a sequence of local two-body deflections along each star's orbit. We explicitly tailor this NR prediction to anisotropic clusters, and compare it with $N$-body realisations of Plummer spheres with varying degrees of anisotropy. The NR theory is shown to recover remarkably well the detailed shape of the orbital diffusion and the associated initial isotropisation, up to a global multiplicative prefactor which increases with anisotropy. Strikingly, a simple effective isotropic prescription provides almost as good a fit, as long as the cluster's anisotropy is not too strong. For these more extreme clusters, accounting for long-range resonant relaxation may be necessary to capture these clusters' long-term evolution.
Katarina Kraljic, Clotilde Laigle, Christophe Pichon, Sebastien Peirani, Sandrine Codis, Junsup Shim, Corentin Cadiou, Dmitri Pogosyan, Stéphane Arnouts, Matthiew Pieri, Vid Iršič, Sean S. Morrison, Jose Oñorbe, Ignasi Pérez-Ràfols, Gavin Dalton The upcoming WEAVE-QSO survey will target a high density of quasars over a large area, enabling the reconstruction of the 3D density field through Lyman-$\alpha$ tomography over unprecedented volumes smoothed on intermediate scales ($\approx$ 16 Mpc/$h$). We produce mocks of the Lyman-$\alpha$ forest using LyMAS, and reconstruct the 3D density field between sightlines through Wiener filtering in a configuration compatible with the future WEAVE-QSO observations. The fidelity of the reconstruction is assessed by measuring one- and two-point statistics from the distribution of critical points in the cosmic web. In addition, initial Lagrangian statistics are predicted from first principles, and measurements of the connectivity of the cosmic web are performed. The reconstruction captures well the expected features in the auto- and cross-correlations of the critical points. This remains true after a realistic noise is added to the synthetic spectra, even though sparsity of sightlines introduces systematics, especially in the cross-correlations of points with mixed signature. Specifically, for walls and filaments, the most striking clustering features could be measured with up to 4 sigma of significance with a WEAVE-QSO-like survey. Moreover, the connectivity of each peak identified in the reconstructed field is globally consistent with its counterpart in the original field, indicating that the reconstruction preserves the geometry of the density field not only statistically, but also locally. Hence the critical points relative positions within the tomographic reconstruction could be used as standard rulers for dark energy by WEAVE-QSO and similar surveys.
Supermassive black holes dominate the gravitational potential in galactic nuclei. In these dense environments, stars follow nearly Keplerian orbits and see their orbital planes relax through the potential fluctuations generated by the stellar cluster itself. For typical astrophysical galactic nuclei, the most likely outcome of this vector resonant relaxation (VRR) is that the orbital planes of the most massive stars spontaneously self-align within a narrow disc. We present a maximum entropy method to systematically determine this long-term distribution of orientations and use it for a wide range of stellar orbital parameters and initial conditions. The heaviest stellar objects are found to live within a thin equatorial disk. The thickness of this disk depends on the stars' initial mass function, and on the geometry of the initial cluster. This work highlights a possible (indirect) novel method to constrain the distribution of intermediate mass black holes in galactic nuclei.
Galaxies form from the accretion of cosmological infall of gas. In the high redshift Universe, most of this gas infall is expected to be dominated by cold filamentary flows which connect deep down inside halos, and, hence, to the vicinity of galaxies. Such cold flows are important since they dominate the mass and angular momentum acquisition that can make up rotationally-supported disks at high-redshifts. We study the angular momentum acquisition of gas into galaxies, and in particular, the torques acting on the accretion flows, using hydrodynamical cosmological simulations of high-resolution zoomed-in halos of a few $10^{11}\,\rm M_\odot$ at $z=2$. Torques can be separated into those of gravitational origin, and hydrodynamical ones driven by pressure gradients. We find that coherent gravitational torques dominate over pressure torques in the cold phase, and are hence responsible for the spin-down and realignment of this gas. Pressure torques display small-scale fluctuations of significant amplitude, but with very little coherence on the relevant galaxy or halo-scale that would otherwise allow them to effectively re-orientate the gas flows. Dark matter torques dominate gravitational torques outside the galaxy, while within the galaxy, the baryonic component dominates. The circum-galactic medium emerges as the transition region for angular momentum re-orientation of the cold component towards the central galaxy's mid-plane.
J. Reddish, K. Kraljic, M. S. Petersen, K. Tep, Y. Dubois, C. Pichon, S. Peirani, F. Bournaud, H. Choi, J. Devriendt, R. Jackson, G. Martin, M. J. Park, M. Volonteri, S. K. Yi We use the NewHorizon simulation to study the redshift evolution of bar properties and fractions within galaxies in the stellar masses range $M_{\star} = 10^{7.25} - 10^{11.4} \ \rm{M}_{\odot}$ over the redshift range $z = 0.25 - 1.3$. We select disc galaxies using stellar kinematics as a proxy for galaxy morphology. We employ two different automated bar detection methods, coupled with visual inspection, resulting in observable bar fractions of $f_{\rm bar} = 0.070_{{-0.012}}^{{+0.018}}$ at $z\sim$ 1.3, decreasing to $f_{\rm bar} = 0.011_{{-0.003}}^{{+0.014}}$ at $z\sim$ 0.25. Only one galaxy is visually confirmed as strongly barred in our sample. This bar is hosted by the most massive disk and only survives from $z=1.3$ down to $z=0.7$. Such a low bar fraction, in particular amongst Milky Way-like progenitors, highlights a missing bars problem, shared by literally all cosmological simulations with spatial resolution $<$100 pc to date. The analysis of linear growth rates, rotation curves and derived summary statistics of the stellar, gas and dark matter components suggest that galaxies with stellar masses below $10^{9.5}-10^{10} \ \rm{M}_{\odot}$ in NewHorizon appear to be too dominated by dark matter to form a bar, while more massive galaxies typically have formed large bulges that prevent bar persistence at low redshift. This investigation confirms that the evolution of the bar fraction puts stringent constraints on the assembly history of baryons and dark matter onto galaxies.
Supermassive black holes in the centre of galaxies dominate the gravitational potential of their surrounding stellar clusters. In these dense environments, stars follow nearly Keplerian orbits, which get slowly distorted as a result of the potential fluctuations generated by the stellar cluster itself as a whole. In particular, stars undergo a rapid relaxation of their eccentricities through both resonant and non-resonant processes. An efficient implementation of the resonant diffusion coefficients allows for detailed and systematic explorations of the parameter space describing the properties of the stellar cluster. In conjunction with recent observations of the S-cluster orbiting SgrA*, this framework can be used to jointly constrain the distribution of the unresolved, old, background stellar cluster and the characteristics of a putative dark cluster. Specifically, we show how this can be used to estimate the typical mass and cuspide exponent of intermediate-mass black holes consistent with the relaxed state of the distribution of eccentricities in the observed S-cluster. This should prove useful in constraining super massive black hole formation scenarios.
Globular clusters contain a finite number of stars. As a result, they inevitably undergo secular evolution (`relaxation') causing their mean distribution function (DF) to evolve on long timescales. On one hand, this long-term evolution may be interpreted as driven by the accumulation of local deflections along each star's mean field trajectory -- so-called `non-resonant relaxation'. On the other hand, it can be thought of as driven by non-local, collectively dressed and resonant couplings between stellar orbits, a process termed `resonant relaxation'. In this paper we consider a model globular cluster represented by a spherical, isotropic isochrone DF, and compare in detail the predictions of both resonant and non-resonant relaxation theories against tailored direct $N$-body simulations. In the space of orbital actions (namely the radial action and total angular momentum), we find that both resonant and non-resonant theories predict the correct morphology for the secular evolution of the cluster's DF, although non-resonant theory over-estimates the amplitude of the relaxation rate by a factor ${\sim 2}$. We conclude that the secular relaxation of hot isotropic spherical clusters is not dominated by collectively amplified large-scale potential fluctuations, despite the existence of a strong ${\ell = 1}$ damped mode. Instead, collective amplification affects relaxation only marginally even on the largest scales. The predicted contributions to relaxation from smaller scale fluctuations are essentially the same from resonant and non-resonant theories.
Associating the formation sites of haloes with the maxima of the smoothed linear density field, we present non-perturbative predictions for the Lagrangian and evolved halo correlation functions that are valid at all separations. In Lagrangian space, we find significant deviations from the perturbative bias calculation at small scales, in particular, a pronounced exclusion region where $\xi=-1$ for maxima of unequal height. Our predictions are in good agreement with the Lagrangian clustering of dark matter proto-haloes reconstructed from N-body simulations. Our predictions for the mean infall and velocity dispersion of haloes, which differ from the local bias expansion, show a similar level of agreement with simulations. Finally, we displace the initial density peaks according to the Zeldovich approximation in order to predict the late-time clustering of dark matter haloes. While we are able to reproduce the early evolution of this conserved set of tracers, our approximation fails at the collapse epoch (z=0) on non-linear scales r<10Mpc/h, emphasizing the need for a non-perturbative treatment of the halo displacement field.
Feedback processes play a fundamental role in the regulation of the star formation (SF) activity in galaxies and, in particular, in the quenching of early-type galaxies (ETGs) as has been inferred by observational and numerical studies of Lambda CDM models. At z = 0, ETGs exhibit well-known fundamental scaling relations, but the connection between them and the physical processes shaping ETG evolution remains unknown.This work aims at studying the impact of the energetic feedback due to active galactic nuclei (AGN) on the formation and evolution of ETGs.We focus on assessing the impact of AGN feedback on the evolution of the mass-plane and the fundamental plane (FP, defined by using mass surface density) as well as on morphology, kinematics, and stellar age across the FP.The Horizon-AGN and Horizon-noAGN cosmological hydrodynamical simulations were performed with identical initial conditions and including the same physical processes except for the activation of the AGN feedback in the former. We select a sample of central ETGs from both simulations using the same criteria and exhaustively study their SF activity, kinematics, and scaling relations for z <= 3. We find that Horizon-AGN ETGs identified at z = 0 follow the observed fundamental scaling relations (mass-plane, FP, mass-size relation) and qualitatively reproduce kinematic features albeit conserving a rotational inner component with a mass fraction regulated by the AGN feedback. AGN feedback seems to be required to reproduce the bimodality in the spin parameter distribution reported by observational works and the mass-size relation (with more massive galaxies having older stellar populations (SPs), larger sizes, and being slower rotators). We study the evolution of the fundamental relations with redshift, finding .Abridged
There have been many studies aiming to reveal the origins of the star-gas misalignment found in galaxies, but there still is a lack of understanding of the contribution from each formation channel candidate. We aim to answer the question by investigating the misaligned galaxies in Horizon-AGN, a cosmological large-volume simulation of galaxy formation. There are 27,903 galaxies of stellar mass $M_* > 10^{10} M_\odot$ in our sample, of which 5,984 are in a group of the halo mass of $M_{200} > 10^{12} M_\odot$. We have identified four main formation channels of misalignment and quantified their level of contribution: mergers (35%), interaction with nearby galaxies (23%), interaction with dense environments or their central galaxies (21%), and secular evolution including smooth accretion from neighboring filaments (21%). We found in the simulation that the gas, rather than stars, is typically more vulnerable to dynamical disturbances; hence, misalignment formation is mainly due to the change in the rotational axis of the gas rather than stars, regardless of the origin. We have also inspected the lifetime (duration) of the misalignment. The decay timescale of the misalignment shows a strong anti-correlation with the kinematic morphology ($V/{\sigma}$) and the cold gas fraction of the galaxy. The misalignment has a longer lifetime in denser regions, which is linked with the environmental impact on the host galaxy. There is a substantial difference in the length of the misalignment lifetime depending on the origin, and it can be explained by the magnitude of the initial position angle offset and the physical properties of the galaxies.
Focusing on both small separations and Baryonic Acoustic Oscillation scales, the cosmic evolution of the clustering properties of peak, void, wall, and filament-type critical points is measured using two-point correlation functions in $\Lambda$CDM dark matter simulations as a function of their relative rarity. A qualitative comparison to the corresponding theory for Gaussian Random fields allows us to understand the following observed features: i) the appearance of an exclusion zone at small separation, whose size depends both on rarity and on the signature (\ie the number of negative eigenvalues) of the critical points involved; ii) the amplification of the Baryonic Acoustic Oscillation bump with rarity and its reversal for cross-correlations involving negatively biased critical points; iii) the orientation-dependent small-separation divergence of the cross-correlations of peaks and filaments (voids and walls) which reflects the relative loci of such points in the filament's (wall's) eigenframe. The most significant features of the correlations are tabulated. The (cross-) correlations involving the most non-linear critical points (peaks, voids) display significant variation with redshift, while those involving less non-linear critical points seem mostly insensitive to redshift evolution, which should prove advantageous to model. The relative distances to the maxima of the peak-to-wall and peak-to-void over that of the peak-to-filament cross-correlation are in ratios of $\sim\sqrt{2}$ and $\sim\sqrt{3}$, respectively which could be interpreted as an indication of the cosmic crystal being on average close to a cubic lattice. The insensitivity to redshift evolution suggests that the absolute and relative clustering of critical points could become a topologically robust alternative to standard clustering techniques when analyzing upcoming large scale surveys such as Euclid or LSST.
R. A. Jackson, S. Kaviraj, G. Martin, J. E. G. Devriendt, A. Slyz, J. Silk, Y. Dubois, S. K. Yi, C. Pichon, M. Volonteri, H. Choi, T. Kimm, K. Kraljic, S. Peirani In the standard Lambda-CDM paradigm, dwarf galaxies are expected to be dark-matter-rich, as baryonic feedback is thought to quickly drive gas out of their shallow potential wells and quench star formation at early epochs. Recent observations of local dwarfs with extremely low dark matter content appear to contradict this picture, potentially bringing the validity of the standard model into question. We use NewHorizon, a high-resolution cosmological simulation, to demonstrate that sustained stripping of dark matter, in tidal interactions between a massive galaxy and a dwarf satellite, naturally produces dwarfs that are dark-matter-deficient, even though their initial dark matter fractions are normal. The process of dark matter stripping is responsible for the large scatter in the halo-to-stellar mass relation in the dwarf regime. The degree of stripping is driven by the closeness of the orbit of the dwarf around its massive companion and, in extreme cases, produces dwarfs with halo-to-stellar mass ratios as low as unity, consistent with the findings of recent observational studies. ~30 per cent of dwarfs show some deviation from normal dark matter fractions due to dark matter stripping, with 10 per cent showing high levels of dark matter deficiency (Mhalo/M*<10). Given their close orbits, a significant fraction of dark-matter-deficient dwarfs merge with their massive companions (e.g. ~70 per cent merge over timescales of ~3.5 Gyrs), with the dark-matter-deficient population being constantly replenished by new interactions between dwarfs and massive companions. The creation of these galaxies is, therefore, a natural by-product of galaxy evolution and their existence is not in tension with the standard paradigm.
Minjung J. Park, Sukyoung K. Yi, Sebastien Peirani, Christophe Pichon, Yohan Dubois, Hoseung Choi, Julien Devriendt, Sugata Kaviraj, Taysun Kimm, Katarina Kraljic, Marta Volonteri Ever since a thick disk was proposed to explain the vertical distribution of the Milky Way disk stars, its origin has been a recurrent question. We aim to answer this question by inspecting 19 disk galaxies with stellar mass greater than $10^{10}\,\rm M_\odot$ in recent cosmological high-resolution zoom-in simulations: Galactica and NewHorizon. The thin and thick disks are reasonably reproduced by the simulations with scale heights and luminosity ratios as observed. We then spatially classify the thin and thick disks and find that the thick disk stars are older, metal-poorer, kinematically-hotter, and higher in accreted star fraction, while both disks are dominated by the stars formed in situ. Half of the in-situ stars in the thick disks are formed before the galaxies develop their disks, and the rest are formed in spatially and kinematically thinner disks and then thickened with time by heating. However, the 19 galaxies have various properties and evolutionary routes, highlighting the need for statistically-large samples to draw general conclusions. We conclude from our simulations that the thin and thick disk components are not entirely distinct in terms of formation processes, but rather markers of the evolution of galactic disks. Moreover, as the combined result of the thickening of the existing disk stars and the continued formation of young thin-disk stars, the vertical distribution of stars does not change much after the disks settle, pointing to the modulation of both orbital diffusion and star formation by the same confounding factor: the proximity of galaxies to marginal stability.
Yohan Dubois, Ricarda Beckmann, Frédéric Bournaud, Hoseung Choi, Julien Devriendt, Ryan Jackson, Sugata Kaviraj, Taysun Kimm, Katarina Kraljic, Clotilde Laigle, Garreth Martin, Min-Jung Park, Sébastien Peirani, Christophe Pichon, Marta Volonteri, Sukyoung K. Yi 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.
We examine how the mass assembly of central galaxies depends on their location in the cosmic web. The HORIZON-AGN simulation is analysed at z~2 using the DISPERSE code to extract multi-scale cosmic filaments. We find that the dependency of galaxy properties on large-scale environment is mostly inherited from the (large-scale) environmental dependency of their host halo mass. When adopting a residual analysis that removes the host halo mass effect, we detect a direct and non-negligible influence of cosmic filaments. Proximity to filaments enhances the build-up of stellar mass, a result in agreement with previous studies. However, our multi-scale analysis also reveals that, at the edge of filaments, star formation is suppressed. In addition, we find clues for compaction of the stellar distribution at close proximity to filaments. We suggest that gas transfer from the outside to the inside of the haloes (where galaxies reside) becomes less efficient closer to filaments, due to high angular momentum supply at the vorticity-rich edge of filaments. This quenching mechanism may partly explain the larger fraction of passive galaxies in filaments, as inferred from observations at lower redshifts.
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.
G. Martin, R. A. Jackson S. Kaviraj, H. Choi, J. E. G. Devriendt, Y. Dubois, T. Kimm, K. Kraljic, S. Peirani, C. Pichon, M. Volonteri, S. K. Yi Dwarf galaxies (M*<10^9 Msun) are key drivers of mass assembly in high mass galaxies, but relatively little is understood about the assembly of dwarf galaxies themselves. Using the \textscNewHorizon cosmological simulation (40 pc spatial resolution), we investigate how mergers and fly-bys drive the mass assembly and structural evolution of around 1000 field and group dwarfs up to z=0.5. We find that, while dwarf galaxies often exhibit disturbed morphologies (5 and 20 per cent are disturbed at z=1 and z=3 respectively), only a small proportion of the morphological disturbances seen in dwarf galaxies are driven by mergers at any redshift (for 10^9 Msun, mergers drive only 20 per cent morphological disturbances). They are instead primarily the result of interactions that do not end in a merger (e.g. fly-bys). Given the large fraction of apparently morphologically disturbed dwarf galaxies which are not, in fact, merging, this finding is particularly important to future studies identifying dwarf mergers and post-mergers morphologically at intermediate and high redshifts. Dwarfs typically undergo one major and one minor merger between z=5 and z=0.5, accounting for 10 per cent of their total stellar mass. Mergers can also drive moderate star formation enhancements at lower redshifts (3 or 4 times at z=1), but this accounts for only a few per cent of stellar mass in the dwarf regime given their infrequency. Non-merger interactions drive significantly smaller star formation enhancements (around two times), but their preponderance relative to mergers means they account for around 10 per cent of stellar mass formed in the dwarf regime.
R. A. Jackson, G. Martin, S. Kaviraj, M. Ramsøy, J. E. G. Devriendt, T. Sedgwick, C. Laigle, H. Choi, R. S. Beckmann, M. Volonteri, Y. Dubois, C. Pichon, S. K. Yi, A. Slyz, K. Kraljic, T. Kimm, S. Peirani, I. Baldry Low-surface-brightness galaxies (LSBGs) -- defined as systems that are fainter than the surface-brightness limits of past wide-area surveys -- form the overwhelming majority of galaxies in the dwarf regime (M* < 10^9 MSun). Using NewHorizon, a high-resolution cosmological simulation, we study the origin of LSBGs and explain why LSBGs at similar stellar mass show the large observed spread in surface brightness. New Horizon galaxies populate a well-defined locus in the surface brightness -- stellar mass plane, with a spread of ~3 mag arcsec^-2, in agreement with deep SDSS Stripe data. Galaxies with fainter surface brightnesses today are born in regions of higher dark-matter density. This results in faster gas accretion and more intense star formation at early epochs. The stronger resultant supernova feedback flattens gas profiles at a faster rate which, in turn, creates shallower stellar profiles (i.e. more diffuse systems) more rapidly. As star formation declines towards late epochs (z<1), the larger tidal perturbations and ram pressure experienced by these systems (due to their denser local environments) accelerate the divergence in surface brightness, by increasing their effective radii and reducing star formation respectively. A small minority of dwarfs depart from the main locus towards high surface brightnesses, making them detectable in past wide surveys. These systems have anomalously high star-formation rates, triggered by recent, fly-by or merger-driven starbursts. We note that objects considered extreme/anomalous at the depth of current datasets, e.g. `ultra-diffuse galaxies', actually dominate the predicted dwarf population and will be routinely visible in future surveys like LSST.
Solène Chabanier, Frédéric Bournaud, Yohan Dubois, Sandrine Codis, Damien Chapon, David Elbaz, Christophe Pichon, Olivier Bressand, Julien Devriendt, Raphael Gavazzi, Katarina Kraljic, Taysun Kimm, Clotilde Laigle, Jean-Baptiste Lekien, Garreth Martin, Nathalie Palanque-Delabrouille, Sébastien Peirani, Pierre-Franck Piserchia, Adrianne Slyz, Maxime Trebitsch, et al (1) 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}$\u2009M$_\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.
Melanie Habouzit, Yuan Li, Rachel S. Somerville, Shy Genel, Annalisa Pillepich, Marta Volonteri, Romeel Davé, Yetli Rosas-Guevara, Stuart McAlpine, Sébastien Peirani, Lars Hernquist, Daniel Anglés-Alcázar, Amy Reines, Richard Bower, Yohan Dubois, Dylan Nelson, Christophe Pichon, Mark Vogelsberger The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the M_BH-M_star relation and the BH mass function. All simulations predict tight M_BH-M_star relations, and struggle to produce the lowest (M_BH< 10^7.5 Msun) in galaxies of M_star~10^10.5-10^11.5 Msun. While the time evolution of the mean M_BH-M_star relation is mild (<1 dex in BH mass for 0<z<5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for M_star<10^10.5 Msun. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with the observations of the local universe, we find an excess of BHs with M_BH>10^9 Msun in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 f_Edd >-2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.
Jaehyun Lee, Jihye Shin, Owain N. Snaith, Yonghwi Kim, C. Gareth Few, Julien Devriendt, Yohan Dubois, Leah M. Cox, Sungwook E. Hong, Oh-Kyoung Kwon, Chan Park, Christophe Pichon, Juhan Kim, Brad K. Gibson, Changbom Park Horizon Run 5 (HR5) is a cosmological hydrodynamical simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1kpc. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of $1049\times119\times127\,{\rm cMpc}^3$.The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, UV background, star formation, supernova feedback, chemical evolution tracking the enrichment of oxygen and iron, the growth of supermassive black holes and feedback from active galactic nuclei (AGN) in the form of a dual jet-heating mode. For this simulation we implemented a hybrid MPI-OMP version of RAMSES, specifically targeted for modern many-core many thread parallel architectures. In addition to the traditional simulation snapshots, light-cone data was generated on the fly. For the post-processing, we extended the Friends-of-Friend (FoF) algorithm and developed a new galaxy finder PGalF to analyse the outputs of HR5. The simulation successfully reproduces observations, such as the cosmic star formation history and connectivity of galaxy distribution, We identify cosmological structures at a wide range of scales, from filaments with a length of several cMpc, to voids with a radius of ~100 cMpc. The simulation also indicates that hydrodynamical effects on small scales impact galaxy clustering up to very large scales near and beyond the baryonic acoustic oscillation (BAO) scale. Hence, caution should be taken when using that scale as a cosmic standard ruler: one needs to carefully understand the corresponding biases. The simulation is expected to be an invaluable asset for the interpretation of upcoming deep surveys of the Universe.
We study shapes and alignments of 45 dark matter (DM) haloes and their brightest cluster galaxies (BCGs) using a sample of 39 massive clusters from Hubble Frontier Field (HFF), Cluster Lensing And Supernova survey with Hubble (CLASH), and Reionization Lensing Cluster Survey (RELICS). We measure shapes of the DM haloes by strong gravitational lensing, whereas BCG shapes are derived from their light profiles in Hubble Space Telescope images. Our measurements from a large sample of massive clusters presented here provide new constraints on dark matter and cluster astrophysics. We find that DM haloes are on average highly elongated with the mean ellipticity of $0.482\pm 0.028$, and position angles of major axes of DM haloes and their BCGs tend to be aligned well with the mean value of alignment angles of $22.2\pm 3.9$ deg. We find that DM haloes in our sample are on average more elongated than their BCGs with the mean difference of their ellipticities of $0.11\pm 0.03$. In contrast, the Horizon-AGN cosmological hydrodynamical simulation predicts on average similar ellipticities between DM haloes and their central galaxies. While such a difference between the observations and the simulation may well be explained by the difference of their halo mass scales, other possibilities include the bias inherent to strong lensing measurements, limited knowledge of baryon physics, or a limitation of cold dark matter.
Marta Volonteri, Hugo Pfister, Ricarda S. Beckman, Yohan Dubois, Monica Colpi, Christopher J. Conselice, Massimo Dotti, Garreth Martin, Ryan Jackson, Katarina Kraljic, Christophe Pichon, Maxime Trebitsch, Sukyoung. K. Yi, Julien Devriendt, Sebastien Peirani 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.
Galaxy merger histories correlate strongly with stellar mass, largely regardless of morphology. Thus, at fixed stellar mass, spheroids and discs share similar assembly histories, both in terms of the frequency of mergers and the distribution of their mass ratios. Since mergers are the principal drivers of disc-to-spheroid morphological transformation, and the most massive galaxies typically have the richest merger histories, it is surprising that discs exist at all at the highest stellar masses (e.g. beyond the knee of the mass function). Using Horizon-AGN, a cosmological hydro-dynamical simulation, we show that extremely massive (M*> 10^11.4 MSun) discs are created via two channels. In the primary channel (accounting for ~70% of these systems and ~8% of massive galaxies) the most recent, significant merger (stellar mass ratio > 1:10) between a massive spheroid and a gas-rich satellite `spins up' the spheroid by creating a new rotational stellar component, leaving a massive disc as the remnant. In the secondary channel (accounting for ~30% of these systems and ~3% of massive galaxies), a system maintains a disc throughout its lifetime, due to an anomalously quiet merger history. Not unexpectedly, the fraction of massive discs is larger at higher redshift, due to the Universe being more gas-rich. The morphological mix of galaxies at the highest stellar masses is, therefore, a strong function of the gas fraction of the Universe. Finally, these massive discs have similar black-hole masses and accretion rates to massive spheroids, providing a natural explanation for why a minority of powerful AGN are surprisingly found in disc galaxies.
Marko Shuntov, J. Pasquet, S. Arnouts, O. Ilbert, M. Treyer, E. Bertin, S. de la Torre, Y. Dubois, D. Fouchez, K. Kraljic, C. Laigle, C. Pichon, D. Vibert Improving distance measurements in large imaging surveys is a major challenge to better reveal the distribution of galaxies on a large scale and to link galaxy properties with their environments. Photometric redshifts can be efficiently combined with the cosmic web (CW) extracted from overlapping spectroscopic surveys to improve their accuracy. We apply a similar method using a new generation of photometric redshifts based on a convolution neural network (CNN). The CNN is trained on the SDSS images with the main galaxy sample (SDSS-MGS, $r \leq 17.8$) and the GAMA spectroscopic redshifts up tor $\sim 19.8$. The mapping of the CW is obtained with 680,000 spectroscopic redshifts from the MGS and BOSS surveys. The redshift probability distribution functions (PDF), which are well calibrated (unbiased and narrow, $\leq 120$ Mpc), intercept a few CW structure along the line of sight. Combining these PDFs with the density field distribution provides new photometric redshifts, $z_{web}$, whose accuracy is improved by a factor of two (i.e.,${\sigma} \sim 0.004(1+z)$) for galaxies with $r \leq 17.8$. For half of them, the distance accuracy is better than 10 cMpc. The narrower the original PDF, the larger the boost in accuracy. No gain is observed for original PDFs wider than 0.03. The final $z_{web}$ PDFs also appear well calibrated. The method performs slightly better for passive galaxies than star-forming ones, and for galaxies in massive groups since these populations better trace the underlying large-scale structure. Reducing the spectroscopic sampling by a factor of 8 still improves the photometric redshift accuracy by 25%. Extending the method to galaxies fainter than the MGS limit still improves the redshift estimates for 70% of the galaxies, with a gain in accuracy of 20% at low $z$ where the resolution of the CW is the highest.
The merging rate of cosmic structures is computed, relying on the Ansatz that they can be predicted in the initial linear density field from the coalescence of critical points with increasing smoothing scale, used here as a proxy for cosmic time. Beyond the mergers of peaks with saddle points (a proxy for halo mergers), we consider the coalescence and nucleation of all sets of critical points, including wall-saddle to filament-saddle and wall-saddle to minima (a proxy for filament and void mergers respectively), as they impact the geometry of galactic infall, and in particular filament disconnection. Analytical predictions of the one-point statistics are validated against multiscale measurements in 2D and 3D realisations of Gaussian random fields (the corresponding code being available upon request) and compared qualitatively to cosmological $N$-body simulations at early times ($z\geq 10$) and large scales ($\geq 5\, \mathrm{Mpc}/h$). The rate of filament coalescence is compared to the merger rate of haloes and the two-point clustering of these events is computed, along with their cross-correlations with critical points. These correlations are qualitatively consistent with the preservation of the connectivity of dark matter haloes, and the impact of the large scale structures on assembly bias. The destruction rate of haloes and voids as a function of mass and redshift is quantified down to $z=0$ for a $\Lambda$CDM cosmology. The one-point statistics in higher dimensions are also presented, together with consistency relations between critical point and critical event counts.
Maxime Trebitsch, Yohan Dubois, Marta Volonteri, Hugo Pfister, Corentin Cadiou, Harley Katz, Joakim Rosdahl, Taysun Kimm, Christophe Pichon, Ricarda S. Beckmann, Julien Devriendt, Adrianne Slyz We present the Obelisk project, a cosmological radiation-hydrodynamics simulation following the assembly and reionization of a protocluster progenitor during the first two billions of years from the big bang, down to $z = 3.5$. The simulation resolves haloes down to the atomic cooling limit, and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The Obelisk project is designed specifically to study the coevolution of high redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to high redshift protoclusters. We find that hydrogen reionization happens inside-out and is completed by $z \sim 6$ in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at $z \sim 4$.
Donghyeon J. Khim, Sukyoung K. Yi, Yohan Dubois, Julia J. Bryant, Christophe Pichon, Scott M. Croom, Joss Bland-Hawthorn, Sarah Brough, Hoseung Choi, Julien Devriendt, Brent Groves, Matt S. Owers, Samuel N. Richards, Jesse van de Sande, Sarah M. Sweet 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.
It is known observationally that the major axes of galaxy clusters and their brightest cluster galaxies are roughly aligned with each other. To understand the origin of the alignment, we identify 40 cluster-sized dark matter (DM) haloes with masses higher than $5\times10^{13}~M_{\odot}$ and their central galaxies (CGs) at $z\approx 0$ in the Horizon-AGN cosmological hydrodynamical simulation. We trace the progenitors at 50 different epochs between $0<z<5$. We then fit their shapes and orientations with a triaxial ellipsoid model. While the orientations of both DM haloes and CGs change significantly due to repeated mergers and mass accretions, their relative orientations are well aligned at each epoch even at high redshifts, $z>1$. The alignment becomes tighter with cosmic time; the major axes of the CGs and their host DM haloes at present are aligned on average within $\sim 30^{\circ}$ in the three dimensional space and $\sim 20^{\circ}$ in the projected plane. The orientations of the major axes of DM haloes on average follow one of the eigen-vectors of the surrounding tidal field that corresponds to the \it slowest collapsing (or even stretching) mode, and the alignment with the tidal field also becomes tighter. This implies that the orientations of CGs and DM haloes at the present epoch are largely imprinted in the primordial density field of the Universe, whereas strong dynamical interactions such as mergers are important to explain their mutual alignment at each epoch.
Elliptical galaxies today appear aligned with the large-scale structure of the Universe, but it is still an open question when they acquire this alignment. Observational data is currently insufficient to provide constraints on the time evolution of intrinsic alignments, and hence existing models range from assuming that galaxies gain some primordial alignment at formation, to suggesting that they react instantaneously to tidal interactions with the large-scale structure. Using the cosmological hydrodynamical simulation Horizon-AGN, we measure the relative alignments between the major axes of galaxies and eigenvectors of the tidal field as a function of redshift. We focus on constraining the time evolution of the alignment of the main progenitors of massive $z=0$ elliptical galaxies, the main weak lensing contaminant at low redshift. We show that this population, which at $z=0$ has a stellar mass above $10^{10.4}$ M$_\odot$, transitions from having no alignment with the tidal field at $z=3$, to a significant alignment by $z=1$. From $z=0.5$ they preserve their alignment at an approximately constant level until $z=0$. We find a mass-dependence of the alignment signal of elliptical progenitors, whereby ellipticals that are less massive today ($10^{10.4}<M/{\rm M}_\odot<10^{10.7}$) do not become aligned till later redshifts ($z<2$), compared to more massive counterparts. We also present an extended study of progenitor alignments in the parameter space of stellar mass and galaxy dynamics, the impact of shape definition and tidal field smoothing.
J. Japelj, C. Laigle, M. Puech, C. Pichon, H. Rahmani, Y. Dubois, J. E. G. Devriendt, P. Petitjean, F. Hammer, E. Gendron, L. Kaper, S. Morris, N. Pirzkal, R. Sánchez-Janssen, A. Slyz, S. D. Vergani, Y. Yang 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]
Katarina Kraljic, Christophe Pichon, Sandrine Codis, Clotilde Laigle, Romeel Davé, Yohan Dubois, Ho Seong Hwang, Dmitri Pogosyan, Stéphane Arnouts, Julien Devriendt, Marcello Musso, Sébastien Peirani, Adrianne Slyz, Marie Treyer We investigate the impact of the number of filaments connected to the nodes of the cosmic web on the physical properties of their galaxies using the Sloan Digital Sky Survey. We compare these measurements to the cosmological hydrodynamical simulations Horizon-(no)AGN and Simba. We find that more massive galaxies are more connected, in qualitative agreement with theoretical predictions and measurements in dark matter only simulation. The star formation activity and morphology of observed galaxies both display some dependence on the connectivity of the cosmic web at fixed stellar mass: less star forming and less rotation supported galaxies also tend to have higher connectivity. These results qualitatively hold both for observed and virtual galaxies, and can be understood given that the cosmic web is the main source of fuel for galaxy growth. The simulations show the same trends at fixed halo mass, suggesting that the geometry of filamentary infall impacts galaxy properties beyond the depth of the local potential well. Based on simulations, it is also found that AGN feedback is key in reversing the relationship between stellar mass and connectivity at fixed halo mass. Technically, connectivity is a practical observational proxy for past and present accretion (minor mergers or diffuse infall).
C. Welker, J. Bland-Hawthorn, J. Van de Sande, C. Lagos, P. Elahi, D. Obreschkow, J. Bryant, C. Pichon, L. Cortese, S. N. Richards, S. M. Croom, M. Goodwin, J. S. Lawrence, S. Sweet, A. Lopez-Sanchez, A. Medling, M. S. Owers, Y. Dubois, J. Devriendt We present the first detection of mass dependent galactic spin alignments with local cosmic filaments with over 2 sigma confidence using IFS kinematics. The 3D network of cosmic filaments is reconstructed on Mpc scales across GAMA fields using the cosmic web extractor DisPerSe. We assign field galaxies from the SAMI survey to their nearest filament segment in 3D and estimate the degree of alignment between SAMI galaxies kinematic spin axis and their nearest filament in projection. Low-mass galaxies align their spin with their nearest filament while higher mass counterparts are more likely to display an orthogonal orientation. The stellar transition mass from the first trend to the second is bracketed between log stellar masses 10.4 and 10.9, with hints of an increase with filament scale. Consistent signals are found in the HorizonAGN cosmological hydrodynamic simulation. This supports a scenario of early angular momentum build-up in vorticity rich quadrants around filaments at low stellar mass followed by progressive flip of spins orthogonal to the cosmic filaments through mergers at high stellar mass. Conversely, we show that dark-matter only simulations post-processed with a semi-analytic model treatment of galaxy formation struggles to reproduce this alignment signal. This suggests that gas physics is key in enhancing the galaxy-filament alignment.
Understanding how rotationally-supported discs transform into dispersion-dominated spheroids is central to our comprehension of galaxy evolution. Morphological transformation is largely merger-driven. While major mergers can efficiently create spheroids, recent work has highlighted the significant role of other processes, like minor mergers, in driving morphological change. Given their rich merger histories, spheroids typically exhibit large fractions of `ex-situ' stellar mass, i.e. mass that is accreted, via mergers, from external objects. This is particularly true for the most massive galaxies, whose stellar masses typically cannot be attained without a large number of mergers. Here, we explore an unusual population of extremely massive (M* > 10^11 MSun) spheroids, in the Horizon-AGN simulation, which exhibit anomalously low ex-situ mass fractions, indicating that they form without recourse to significant merging. These systems form in a single minor-merger event (with typical merger mass ratios of 0.11 - 0.33), with a specific orbital configuration, where the satellite orbit is virtually co-planar with the disc of the massive galaxy. The merger triggers a catastrophic change in morphology, over only a few hundred Myrs, coupled with strong in-situ star formation. While this channel produces a minority (~5 per cent) of such galaxies, our study demonstrates that the formation of at least some of the most massive spheroids need not involve major mergers -- or any significant merging at all -- contrary to what is classically believed.
Hydrodynamical cosmological simulations have recently made great advances in reproducing galaxy mass assembly over cosmic time - as often quantified from the comparison of their predicted stellar mass functions to observed stellar mass functions from data. In this paper we compare the clustering of galaxies from the hydrodynamical cosmological simulated lightcone Horizon-AGN, to clustering measurements from the VIDEO survey observations. Using mocks built from a VIDEO-like photometry, we first explore the bias introduced into clustering measurements by using stellar masses and redshifts derived from SED-fitting, rather than the intrinsic values. The propagation of redshift and mass statistical and systematic uncertainties in the clustering measurements causes us to underestimate the clustering amplitude. We find then that clustering and halo occupation distribution (HOD) modelling results are qualitatively similar in Horizon-AGN and VIDEO. However at low stellar masses Horizon-AGN underestimates the observed clustering by up to a factor of ~3, reflecting the known excess stellar mass to halo mass ratio for Horizon-AGN low mass haloes, already discussed in previous works. This reinforces the need for stronger regulation of star formation in low mass haloes in the simulation. Finally, the comparison of the stellar mass to halo mass ratio in the simulated catalogue, inferred from angular clustering, to that directly measured from the simulation, validates HOD modelling of clustering as a probe of the galaxy-halo connection.
The intrinsic correlations of galaxy shapes and orientations across the large-scale structure of the Universe are a known contaminant to weak gravitational lensing. They are known to be dependent on galaxy properties, such as their mass and morphologies. The complex interplay between alignments and the physical processes that drive galaxy evolution remains vastly unexplored. We assess the sensitivity of intrinsic alignments (shapes and angular momenta) to Active Galactic Nuclei -AGN- feedback by comparing galaxy alignment in twin runs of the cosmological hydrodynamical Horizon simulation, which do and do not include AGN feedback respectively. We measure intrinsic alignments in three dimensions and in projection at z=0 and z=1. We find that the projected alignment signal of all galaxies with resolved shapes with respect to the density field in the simulation is robust to AGN feedback, thus giving similar predictions for contamination to weak lensing. The relative alignment of galaxy shapes around galaxy positions is however significantly impacted, especially when considering high-mass ellipsoids. Using a sample of galaxy "twins" across simulations, we determine that AGN changes both the galaxy selection and their actual alignments. Finally, we measure the alignments of angular momenta of galaxies with their nearest filament. Overall, these are more significant in the presence of AGN as a result of the higher abundance of massive pressure-supported galaxies.
We study the Intra-Halo Stellar Component (IHSC) of Milky Way-mass systems up to galaxy clusters in the Horizon-AGN cosmological hydrodynamical simulation. We identify the IHSC using an improved phase-space galaxy finder algorithm which provides an adaptive, physically motivated and shape-independent definition of this stellar component, that can be applied to halos of arbitrary masses. We explore the IHSC mass fraction-total halo's stellar mass, $f_{M*,IHSC}-M*$, relation and the physical drivers of its scatter. We find that on average the $f_{M*,IHSC}$ increases with $M_{*,tot}$, with the scatter decreasing strongly with mass from 2 dex at $M_{*,tot}\sim10^{11}M_\odot$ to 0.3 dex at group masses. At high masses, $M_{*,tot}>10^{11.5}M_\odot$, $f_{M*,IHSC}$ increases with the number of substructures, and with the mass ratio between the central galaxy and largest satellite, at fixed $M_{*,tot}$. From mid-size groups and systems below $M_{*,tot}<10^{12}M_\odot$, we find that the central galaxy's stellar rotation-to-dispersion velocity ratio, V/\sigma, displays the strongest (anti)-correlation with $f_{M*,IHSC}$ at fixed $M_{*,tot}$ of all the galaxy and halo properties explored, transitioning from $f_{M*,IHSC}$<0.1% for high V/\sigma, to $f_{M*,IHSC}\sim5$% for low V/\sigma galaxies. By studying the $f_{M*,IHSC}$ temporal evolution, we find that, in the former, mergers not always take place, but if they did, they happened early (z>1), while the high $f_{M*,IHSC}$ population displays a much more active merger history. In the case of massive groups and galaxy clusters, $M_{*,tot}>10^{12}M_\odot$, a fraction $f_{M*,IHSC}\sim$10-20% is reached at $z\sim1$ and then they evolve across lines of constant $f_{M*,IHSC}$ modulo some small perturbations. Because of the limited simulation's volume, the latter is only tentative and requires a larger sample of simulated galaxy clusters to confirm.
We present analysis of two- and three-point correlation functions of Ly$\alpha$ forest (at $2\le z\le 2.5$) using X-Shooter spectra of three background quasar triplets probing transverse separations of 0.5-1.6 pMpc. We present statistics based on transmitted flux and clouds identified using Voigt profile fitting. We show that the observed two-, three-point correlation functions and reduced three-point correlation (i.e Q) are well reproduced by our simulations. We assign probabilities for realising all the observed correlation properties simultaneously using our simulations. Our simulations suggest an increase in correlation amplitudes and Q with increasing $N_{\rm HI}$. We roughly see this trend in the observations too. We identify a concurrent gap of 17$\mathring{A}$ (i.e 14.2 $h^{-1}$cMpc, one of the longest reported) wide in one of the triplets. Such gap is realised only in 14.2% of our simulated sightlines and most of the time belongs to a void in the matter distribution. In the second triplet, we detect DLAs along all three sightlines (with spatial separations 0.64 to 1.6 pMpc) within a narrow redshift interval (i.e $\Delta z = 0.088$). Detection of a foreground quasar ($\sim$ 1 pMpc from the triplet sightlines) and excess partial Lyman Limit systems around these DLAs suggest that we may be probing a large over-dense region. We also report positive CIV- CIV correlations up to $\sim 500$ $km s^{-1}$ only in the longitudinal direction. Additionally, we conclude a positive CIV-Ly$\alpha$ correlations for higher $N_{\rm HI}$ thresholds up to a scale of $\sim 1000$ $km s^{-1}$ both in transverse and longitudinal directions.
We study the spin alignment of galaxies and halos with respect to filaments and walls of the cosmic web, identified with DisPerSE, using the SIMBA simulation from z=0-2. Massive halos' spins are oriented perpendicularly to their closest filament's axis and walls, while low mass halos tend to have their spins parallel to filaments and in the plane of walls. A similar mass-dependent spin flip is found for galaxies, albeit with a weaker signal particularly at low mass and low-z, suggesting that galaxies' spins retain memory of their larger-scale environment. Low-z star-forming and rotation-dominated galaxies tend to have spins parallel to nearby filaments, while quiescent and dispersion-dominated galaxies show preferentially perpendicular orientation; the star formation trend can be fully explained by the stellar mass correlation, but the morphology trend cannot. There is a strong dependence on HI mass, such that high-HI galaxies tend to have parallel spins while low-HI galaxies are perpendicular, which persists even when matching samples in stellar mass, suggesting that HI content traces anisotropic infall more faithfully than the stellar component. Finally, at fixed stellar mass, the strength of spin alignments correlates with the filament's density, with parallel alignment for galaxies in high density environments. These findings are consistent with conditional tidal torque theory, and highlight a significant correlation between galactic spin and the larger scale tides that are important e.g. for interpreting weak lensing studies. SIMBA allows us to rule out numerical grid locking as the cause of previously-seen low mass alignment.
Iary Davidzon, Clotilde Laigle, Peter L. Capak, Olivier Ilbert, Daniel C. Masters, Shoubaneh Hemmati, Nikolaos Apostolakos, Jean Coupon, Sylvain de la Torre, Julien Devriendt, Yohan Dubois, Daichi Kashino, Stephane Paltani, Christophe Pichon Using the Horizon-AGN hydrodynamical simulation and self-organising maps (SOMs), we show how to compress the complex data structure of a cosmological simulation into a 2-d grid which is much easier to analyse. We first verify the tight correlation between the observed 0.3$\!-\!5\mu$m broad-band colours of Horizon-AGN galaxies and their high-resolution spectra. The correlation is found to extend to physical properties such as redshift, stellar mass, and star formation rate (SFR). This direct mapping from colour to physical parameter space is shown to work also after including photometric uncertainties that mimic the COSMOS survey. We then label the SOM grid with a simulated calibration sample and estimate redshift and SFR for COSMOS-like galaxies up to $z\sim3$. In comparison to state-of-the-art techniques based on synthetic templates, our method is comparable in performance but less biased at estimating redshifts, and significantly better at predicting SFRs. In particular our "data-driven" approach, in contrast to model libraries, intrinsically allows for the complexity of galaxy formation and can handle sample biases. We advocate that obtaining the calibration for this method should be one of the goals of next-generation galaxy surveys.
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.
E. Darragh-Ford, C. Laigle, G. Gozaliasl, C. Pichon, J. Devriendt, A. Slyz, S. Arnouts, Y. Dubois, A. Finoguenov, R. Griffiths, K. Kraljic, H. Pan, S. Peirani, F. Sarron Cosmic filaments are the channel through which galaxy groups assemble their mass. Cosmic connectivity, namely the number of filaments connected to a given group, is therefore expected to be an important ingredient in shaping group properties. The local connectivity is measured in COSMOS around X-Ray detected groups between redshift 0.5 and 1.2. To this end, large-scale filaments are extracted using the accurate photometric redshifts of the COSMOS2015 catalogue in two-dimensional slices of thickness 120 comoving Mpc centred on the group's redshift. The link between connectivity, group mass and the properties of the brightest group galaxy (BGG) is investigated. The same measurement is carried out on mocks extracted from the lightcone of the hydrodynamical simulation Horizon-AGN in order to control systematics. More massive groups are on average more connected. At fixed group mass in low-mass groups, BGG mass is slightly enhanced at high connectivity, while in high mass groups BGG mass is lower at higher connectivity. Groups with a star-forming BGG have on average a lower connectivity at given mass. From the analysis of the Horizon-AGN simulation, we postulate that different connectivities trace different paths of group mass assembly: at high group mass, groups with higher connectivity are more likely to have grown through a recent major merger, which might be in turn the reason for the quenching of the BGG. Future large-field photometric surveys, such as Euclid and LSST, will be able to confirm and extend these results by probing a wider mass range and a larger variety of environment.
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.