Search

(Reduced) Using the NewHorizon simulation, we have studied ten gas-versus-star counter-rotating galaxies in field environments with a stellar mass of M*~[1-5]x10^10 Msun. For all of them, the retrograde accretion of gas either from gas stripping from a nearby companion or from the circumgalactic medium is the starting point of the formation process. Then follows the co-existence of two distinct disks of gas rotating in opposite directions, the pre-existing one in the inner parts and the accreted gas in the outer parts of the galaxy. The latter progressively replaces the former leading to the final gas-star kinetic misalignment configuration. During the process, the star formation is first enhanced and then progressively decreases. We roughly estimate that a higher fraction of the pre-existing gas is converted into stars rather than being expelled. We also found that the black hole activity (BH) tends to be enhanced during the removal of the pre-existing gas. Furthermore, our analysis suggests that the formation of a counter-rotating gas component is always accompanied with the formation of counter-rotating stellar disks. These stellar disks can have diverse properties but host in general a younger and more metal rich population of stars with respect to the main disc, depending on the star formation history and BH activity. The central part of counter-rotating disks tend also to be characterized by a younger population, an enhanced star formation rate and a higher metallicity than their outer parts. The high metallicity comes the progressive metal enrichment of the accreted gas by mixing with the pre-existing gas and by supernovae activity as it sinks toward the center of the galaxy. In case of major mergers, a large amount of accreted stars from the companion would be distributed at large distances from the remnant center due to conservation of the initial orbital angular momentum., Comment: accepted for publication in A&A

Among the various implications of the spiral arms, it has been observed in the Milky Way disc that the distribution of radial actions from the Gaia Data Release 3 exhibits structures that may be related to the spiral arms. Our goal is to investigate the relationship between regions of low radial action identified in simulated discs and the location of the spiral arms, such as that suggested in Palicio et al. (2023) for the Galaxy. For a sample of 23 simulated spiral galaxies, we modelled the axisymmetric component of their gravitational potential to compute the radial action of their stellar particles using the Staeckel fudge. The spatial distribution of the radial action was then compared to the location of the spiral arms, identified as overdensities in the stellar surface density using a kernel density estimator. Our analysis reveals a strong correlation between the radial action distribution and the spiral arms in 18 of 23 simulated galaxies. However, notable discrepancies are observed in the remaining five, since they are profoundly out-of-equilibrium systems, such as galaxies influenced by external interactions or spiral arms still in the process of winding up. We have confirmed that, in general, there is a tendency of spatial correlation between spiral arms and stellar populations featuring low values of the radial action, as discussed in Palicio et al. (2023). However, discrepancies between features in the radial action distribution and the spiral structure can be interpreted as signatures of recent disturbances, a scenario applicable to the Milky Way. Furthermore, populations at least as old as 3 Gyr trace the spiral arms with no significant misalignment across age bins. A linear relation between the maximum value of the radial action of the spiral arms and the vertical scale-length is found, which is also satisfied by the Milky Way., Comment: Accepted for publication in A&A. The manuscript has been updated with the new layout

Recent observations reveal a rapid dust build-up in high-redshift galaxies (z > 4), challenging current models of galaxy formation. While our understanding of dust production and destruction in the interstellar medium (ISM) is advancing, probing baryonic processes in the early Universe remains a complex task. We characterize the evolution of 98 z~5 star-forming galaxies observed as part of the ALPINE survey by constraining the physical processes underpinning the gas and dust production, consumption, and destruction in their ISM. We make use of chemical evolution models to simultaneously reproduce the observed dust and gas content. For each galaxy, we estimate initial gas mass, inflows and outflows, and efficiencies of dust growth and destruction. We test the models with the canonical Chabrier and top-heavy initial mass functions (IMFs), with the latter enabling rapid dust production on shorter timescales. Our models successfully reproduce gas and dust content in older galaxies (> 600 Myr) regardless of the IMF, with Type II SNe as the primary dust source and no dust growth in ISM with moderate inflow of primordial gas. In case of intermediate-age galaxies (300 - 600 Myr), we reproduce the gas and dust content through Type II SNe and dust growth in ISM, though we observe an over-prediction of dust mass in older galaxies, potentially indicating an unaccounted dust destruction mechanism and/or an overestimation of the observed dust masses. The number of young galaxies (< 300 Myr) reproduced, increases for models assuming top-heavy IMF but with maximal prescriptions of dust production. Galactic outflows are necessary to reproduce observed gas and dust masses. The Chabrier IMF models reproduce 65% of galaxies, while top-heavy IMF models improve this to 93%, easing tensions with observations. Upcoming JWST data will refine these models by resolving degeneracies in intrinsic galaxy properties., Comment: Accepted for publication in A&A

We introduce a new concept -- termed "planarity" -- which aims to quantify planar structure in galaxy satellite systems without recourse to the number or thickness of planes. We use positions and velocities from the Gaia EDR3 to measure planarity in Milky Way (MW) satellites and the extent to which planes within the MW system are kinematically supported. We show that the position vectors of the MW satellites exhibit strong planarity but the velocity vectors do not, and that kinematic coherence cannot, therefore, be confirmed from current observational data. We then apply our methodology to NewHorizon, a high-resolution cosmological simulation, to compare satellite planarity in MW-like galaxies in a {\Lambda}CDM-based model to that in the MW satellite data. We demonstrate that kinematically supported planes are common in the simulation and that the observed planarity of MW satellites is not in tension with the standard {\Lambda}CDM paradigm.

We study the stellar mass function (SMF) and the co-evolution with dark matter halos via abundance matching in the largest redshift range to date $0.25$, we find increased abundances of massive (log$\, M_{\star}/M_{\odot}>10.5$) implying integrated star formation efficiencies (SFE) $\epsilon_{\star}\equiv M_{\star}\, f_{\rm b}^{-1} M_{\rm halo}^{-1} \gtrsim 0.5$. We find a flattening of the SMF at the high-mass end that is better described by a double power law at $z>5.5$. At $z \lesssim 5.5$ it transitions to a Schechter law which coincides with the emergence of the first massive quiescent galaxies in the Universe. We trace the cosmic stellar mass density (SMD) and infer the star formation rate density (SFRD), which at $z>7.5$ agrees remarkably with recent \JWST{} UV luminosity function-derived estimates. However, at $z \lesssim 3.5$, we find significant tension ($\sim 0.3$ dex) with the cosmic star formation (SF) history from instantaneous SF measures, the causes of which remain poorly understood. We infer the stellar-to-halo mass relation (SHMR) and the SFE from abundance matching out to $z=12$, finding a non-monotonic evolution. The SFE has the characteristic strong dependence with mass in the range of $0.02 - 0.2$, and mildly decreases at the low mass end out to $z\sim3.5$. At $z\sim3.5$ the SFE increases sharply from $\sim 0.1$ to approach high SFE of $0.8-1$ by $z\sim 10$ for log$(M_{\rm h}/M_{\odot})\approx11.5$, albeit with large uncertainties. Finally, we use the SHMR to track the SFE and stellar mass growth throughout the halo history and find that they do not grow at the same rate -- from the earliest times up until $z\sim3.5$ the halo growth rate outpaces galaxy assembly, but at $z>3.5$ halo growth stagnates and accumulated gas reservoirs keep the SF going and galaxies outpace halos.

The stellar mass-star formation rate (M$_\star$ - SFR) plane is essential for distinguishing galaxy populations, but how galaxies move within this plane over cosmic time remains unclear. This study aims to describe galaxy migrations in the M$_\star$ - SFR plane by reconstructing star formation histories (SFHs) for a sample of galaxies out to redshift $z < 4$. This provides insights into the physical processes driving star formation. We use data from the COSMOS field, selecting 299131 galaxies at $z < 4$ with COSMOS-Web NIRCam data (m$_\mathrm{F444W} < 27$) over 0.54 deg$^2$. Using the SED modeling code CIGALE with non-parametric SFHs, we derive physical properties and migration vectors for these galaxies. These vectors describe the direction and velocity of evolutionary paths across the M$_\star$ - SFR plane. To assess the accuracy of these vectors, we compare them to results from the Horizon-AGN simulation. Galaxies within the main sequence show low migration amplitudes and dispersed movement directions, indicating oscillation within the main sequence. Most progenitors were already on the main sequence a billion years earlier. Starburst galaxies assembled half their mass in the last 350 Myr and originated from the main sequence. Passive galaxies show uniformly declining SFHs and include massive galaxies already in the passive region at $z = 3.5-4$, with increasing density over time. Using reconstructed SFHs up to $z < 4$, we propose a coherent picture of how galaxies migrate over cosmic time in the M$_\star$ - SFR plane, highlighting the connection between major phases in the star-formation history.

We propose a novel method to reconstruct the full posterior distribution of the star formation histories (SFHs) of galaxies from broad-band photometry. Our method combines simulation-based inference (SBI) using a neural network trained with SFHs and photometry from the {\sc Horizon-AGN} hydrodynamical cosmological simulation. We apply it to reconstruct SFHs using COSMOS2020 photometry at redshift $0 9$ had a first event of mass assembly around $z\sim 3$, independent of mass. This work represents a pilot study for the future analysis of the \textit{Euclid} Deep fields that will reach similar depths in alike set of photometric bands, but with over an order-of-magnitude larger area, opening the possibility of deriving SFHs for millions of galaxies in a robust manner.

The James Webb Space telescope has discovered an abundant population of broad line emitters, typical signposts for Active Galactic Nuclei (AGN). Many of these sources have red colors and a compact appearance that has led to naming them `Little Red Dots'. In this paper we develop a detailed framework to estimate the photometry of AGN embedded in galaxies extracted from the \Obelisk{} cosmological simulation to understand the properties of color-selected Little Red Dots (cLRDs) in the context of the full AGN and massive black hole population. We find that using realistic spectral energy distributions (SEDs) and attenuation for AGN we can explain the shape of the cLRD SED as long as galaxies host a sufficiently luminous AGN that is not too much or too little attenuated. When attenuation is too low or too high, AGN do not enter the cLRD selection, because the AGN dominates over the host galaxy too much in blue filters, or it does not contribute to photometry anywhere, respectively. cLRDs are also characterized by high Eddington ratios, possibility super-Eddington, and/or high ratios between black hole and stellar mass., Comment: Accepted for publication in A&A

Galaxy mergers represent a fundamental physical process under hierarchical structure formation, but their role in triggering AGNs is still unclear. We aim to investigate the merger-AGN connection using state-of-the-art observations and novel methods in detecting mergers and AGNs. We selected stellar mass-limited samples at redshift z<1 from KiDS, focusing on the KiDS-N-W2 field with a wide range of multi-wavelength data. Three AGN types, selected in the MIR, X-ray, and via SED modelling, were analysed. To identify mergers, we used convolutional neural networks trained on two cosmological simulations. We created mass and redshift-matched control samples of non-mergers and non-AGNs. We observe a clear AGN excess (a factor of 2-3) in mergers with respect to non-mergers for the MIR AGNs, and a mild excess for the X-ray and SED AGNs, indicating that mergers could trigger all 3 types but are more connected with the MIR AGNs. About half of the MIR AGNs are in mergers but it is unclear whether mergers are the main trigger. For the X-ray and SED AGNs, mergers are unlikely to be the dominant trigger. We also explore the relation using the continuous AGN fraction $f_{AGN}$ parameter. Mergers exhibit a clear excess of high $f_{AGN}$ values relative to non-mergers, for all AGNs. We unveil the first merger fraction $f_{merg}-f_{AGN}$ relation with two distinct regimes. When the AGN is not dominant, the relation is only mildly increasing or even flat, with the MIR AGNs showing the highest $f_{merg}$. In the regime of very dominant AGNs ($f_{AGN}\geq0.8$), $f_{merg}$ shows a steeply rising trend with increasing $f_{AGN}$ for all AGN types. These trends are also seen when plotted against AGN bolometric luminosity. We conclude that mergers are most connected with dust-obscured AGNs (linked to a fast-growing phase of the SMBH) and are the main or even the sole fuelling mechanism of the most powerful AGNs., Comment: 24 pages, 17 figures, 6 Appendix pages. Public available multi-wavelength catalogue of galaxies. Accepted for publication in A&A

The Kennicutt-Schmidt (KS) relation between the gas and the star formation rate (SFR) surface density ($\Sigma_{\rm gas}$-$\Sigma_{\rm SFR}$) is essential to understand star formation processes in galaxies. So far, it has been measured up to z~2.5 in main-sequence galaxies. In this letter, we aim to put constraints at z~4.5 using a sample of four massive main-sequence galaxies observed by ALMA at high resolution. We obtained ~0.3"-resolution [CII] and continuum maps of our objects, which we then converted into gas and obscured SFR surface density maps. In addition, we produced unobscured SFR surface density maps by convolving Hubble ancillary data in the rest-frame UV. We then derived the average $\Sigma_{\rm SFR}$ in various $\Sigma_{\rm gas}$ bins, and estimated the uncertainties using a Monte Carlo sampling. Our galaxy sample follows the KS relation measured in main-sequence galaxies at lower redshift and is slightly lower than predictions from simulations. Our data points probe the high end both in terms of $\Sigma_{\rm gas}$ and $\Sigma_{\rm gas}$, and gas depletion timescales (285-843 Myr) remain similar to z~2 objects. However, three of our objects are clearly morphologically disturbed, and we could have expected shorter gas depletion timescales (~100 Myr) similar to merger-driven starbursts at lower redshifts. This suggests that the mechanisms triggering starbursts at high redshift may be different than in the low- and intermediate-z Universe., Comment: 9 pages, 7 figures, 2 tables, accepted by A&A (letter)

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., Comment: 13 pages, 11 figures, Accepted to MNRAS

Star formation rate functions (SFRFs) give an instantaneous view of the distribution of star formation rates (SFRs) in galaxies at different epochs. They are a complementary and more stringent test for models than the galaxy stellar mass function, which gives an integrated view of the past star formation activity. However, the exploration of SFRFs has been limited thus far due to difficulties in assessing the SFR from observed quantities and probing the SFRF over a wide range of SFRs. We overcome these limitations thanks to an original method that predicts the infrared luminosity from the rest-frame UV/optical color of a galaxy and then its SFR over a wide range of stellar masses and redshifts. We applied this technique to the deep imaging survey HSC-CLAUDS combined with near-infrared and UV photometry. We provide the first SFR functions with reliable measurements in the high- and low-SFR regimes up to $z=2$ and compare our results with previous observations and four state-of-the-art hydrodynamical simulations., Comment: 29 pages, 19 figures

Massive black-hole (BH) mergers are predicted to be powerful sources of low-frequency gravitational waves (GWs). Coupling the detection of GWs with an electromagnetic (EM) detection can provide key information about merging BHs and their environments. We study the high-resolution cosmological radiation-hydrodynamics simulation OBELISK, run to redshift $z=3.5$, to assess the GW and EM detectability of high-$z$ BH mergers, modelling spectral energy distribution and obscuration. For EM detectability, we further consider sub-grid dynamical delays in postprocessing. We find that most of the merger events can be detected by LISA, except for high-mass mergers with very unequal mass ratios. Intrinsic binary parameters are accurately measured, but the sky localisation is poor generally. Only $\sim 40\%$ of these high-$z$ sources have a sky localisation better than $10\,\mathrm{deg}^2$. Merging BHs are hard to detect in the restframe UV since they are fainter than the host galaxies, which at high $z$ are star-forming. A significant fraction, $15-35\%$, of BH mergers instead outshine the galaxy in X-rays, and about $5-15\%$ are sufficiently bright to be detected with sensitive X-ray instruments. If mergers induce an Eddington-limited brightening, up to $30\%$ of sources can become observable. The transient flux change originating from such a brightening is often large, allowing $4-20\%$ of mergers to be detected as EM counterparts. A fraction, $1-30\%$, of mergers are also detectable at radio frequencies. Observable merging BHs tend to have higher accretion rates and masses and are overmassive at a fixed galaxy mass with respect to the full population. Most EM-observable mergers can also be GW-detected with LISA, but their sky localisation is generally poorer. This has to be considered when using EM counterparts to obtain information about the properties of merging BHs and their environment., Comment: 18 pages, 13 figures. Matches accepted version (A&A, 22/06/23)

Massive black hole (BH) mergers will be key targets of future gravitational wave and electromagnetic observational facilities. In order to constrain BH evolution with the information extracted from BH mergers, one must take into account the complex relationship between the population of merging BHs and the global BH population. We analysed the high-resolution cosmological radiation-hydrodynamics simulation OBELISK, run to redshift $z=3.5$, to study the properties of the merging BH population, and its differences with the underlying global BH population in terms of BH and galaxy properties. In post-processing, we calculated dynamical delays between the merger in the simulation at the resolution limit and the actual coalescence well below the resolution scale. We find that merging BHs are hosted in relatively massive galaxies with stellar mass $M_\ast\gtrsim10^9\,M_\odot$. Given that galaxy mass is correlated with other BH and galaxy properties, BH mergers tend to also have a higher total BH mass and higher BH accretion rates than the global population of main BHs. These differences generally disappear if the merger population is compared with a BH population sampled with the same galaxy mass distribution as merger hosts. Galaxy mergers can temporarily boost the BH accretion rate and the host's star formation rate, which can remain active at the BH merger if sub-resolution delays are not taken into account. When dynamical delays are taken into account, the burst has generally faded by the time the BHs merge. BH spins are followed self-consistently in the simulation under the effect of accretion and BH mergers. We find that merging BHs have higher spins than the global population, but similar or somewhat lower spins compared to a mass-matched sample. For our sample, mergers tend to decrease the spin of the final BH remnant., Comment: 13 pages, 9 figures. Matches accepted version (A&A, 28/03/23)

Recent observational and theoretical studies have suggested that supermassive black holes (SMBHs) grow mostly through non-merger (`secular') processes. Since galaxy mergers lead to dynamical bulge growth, the only way to observationally isolate non-merger growth is to study galaxies with low bulge-to-total mass ratio (e.g. B/T < 10%). However, bulge growth can also occur due to secular processes, such as disk instabilities, making disk-dominated selections a somewhat incomplete way to select merger-free systems. Here we use the Horizon-AGN simulation to select simulated galaxies which have not undergone a merger since z = 2, regardless of bulge mass, and investigate their location on typical black hole-galaxy scaling relations in comparison to galaxies with merger dominated histories. While the existence of these correlations has long been interpreted as co-evolution of galaxies and their SMBHs driven by galaxy mergers, we show here that they persist even in the absence of mergers. We find that the correlations between SMBH mass and both total mass and stellar velocity dispersion are independent of B/T ratio for both merger-free and merger-dominated galaxies. In addition, the bulge mass and SMBH mass correlation is still apparent for merger-free galaxies, the intercept for which is dependent on B/T. Galaxy mergers reduce the scatter around the scaling relations, with merger-free systems showing broader scatter. We show that for merger-free galaxies, the co-evolution is dominated by radio-mode feedback, and suggest that the long periods of time between galaxy mergers make an important contribution to the co-evolution between galaxies and SMBHs in all galaxies., Comment: RJS and RSB are joint first authors. 12 pages, 7 figures. Accepted 2023 February 13. Received 2023 February 2; in original form 2022 November 24

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., Comment: RSB and RJS are joint first authors. 11 pages, 7 figures. Accepted 2023 June 9. Received 2023 April 25; in original form 2022 November 24

While it is well established that supermassive black holes (SMBHs) co-evolve with their host galaxy, it is currently less clear how lower mass black holes, so-called intermediate mass black holes (IMBHs), evolve within their dwarf galaxy hosts. In this paper, we present results on the evolution of a large sample of IMBHs from the NewHorizon simulation. We show that occupation fractions of IMBHs in dwarf galaxies are at least 50 percent for galaxies with stellar masses down to 1E6 Msun, but BH growth is very limited in dwarf galaxies. In NewHorizon, IMBH growth is somewhat more efficient at high redshift z = 3 but in general IMBH do not grow significantly until their host galaxy leaves the dwarf regime. As a result, NewHorizon under-predicts observed AGN luminosity function and AGN fractions. We show that the difficulties of IMBH to remain attached to the centres of their host galaxies plays an important role in limiting their mass growth, and that this dynamic evolution away from galactic centres becomes stronger at lower redshift., Comment: 15 pages, submitted to MNRAS

Dynamically, cosmic rays with energies above about one GeV/nucleon may be important agents of galaxy evolution. Their pressures compare with the thermal and magnetic ones impacting galactic gas accretion, fountains and galactic outflows, and alter the mass cycling between the gas phases, its efficiency depends on the properties of CR transport in the different media. We aim to study the dynamical role of CRs in shaping the interstellar medium of a galaxy when changing their propagation mode. We perform MHD simulations with the AMR code RAMSES of the evolution of the same isolated galaxy (dwarf galaxy of $10^{11}$ M$_{\odot}$ down to 9-pc resolution) and compare the impact of the simplest cosmic-ray transport assumption of uniform diffusion. We have also updated the observational relation seen between the $\gamma$-ray luminosities and SFR of galaxies using the latest detection of Fermi LAT sources. We find that the radial and vertical distributions, and mass fractions of the gas in the different phases are marginally altered when changing CR transport. We observe positive feedback of CR on the amplification of the magnetic field in the inner half of the galaxy, except for fast isotropic diffusion. The increase in CR pressure for slow or anisotropic diffusion can suppress star formation by up to 50\%, but the dual effect of cosmic-ray pressure and magnetic amplification can reduce star formation by a factor 2.5. The $\gamma$-ray luminosities and SFR of the simulated galaxies are fully consistent with the trend seen in the observations in the case of anisotropic $10^{27.5-29}$ cm$^2$ s$^{-1}$ diffusion and for isotropic diffusion slower or equal to $3 \times 10^{28}$cm$^2$ s$^{-1}$. These results, therefore, do not confirm claims of very fast $10^{29-31}$ cm$^2$ s$^{-1}$ diffusion to match the Fermi LAT observations.

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., Comment: Submitted to MNRAS (one iteration with the referee) - 29 pages

We use the COSMOS2020 catalogue to measure the stellar-to-halo mass relation (SHMR) divided by central and satellite galaxies from $z=0.2$ to $z = 5.5$. Starting from accurate photometric redshifts we measure the near-infrared selected two-point angular correlation and stellar mass functions in ten redshift bins and fit them with an HOD-based model. At each redshift, we measure the ratio of stellar mass to halo mass, $M_*/M_h$, which shows the characteristic strong dependence of halo mass with a peak at $M_h^{\rm peak} \sim 10^{12}\, M_{\odot}$. Our results are in accordance with the scenario in which the peak of star-formation efficiency moves towards more massive halos at higher redshifts. We also measure the fraction of satellites as a function of stellar mass and redshift. For all stellar mass thresholds the satellite fraction decreases at higher redshifts. At a given redshift there is a higher fraction of low-mass satellites. The satellite contribution to the total stellar mass budget in halos becomes more important than centrals at halo masses of about $M_h > 10^{13} \, M_{\odot}$ and always stays below by peak, indicating that quenching mechanisms are present in massive halos that keep the star-formation efficiency low. Finally, we compare our results with three hydrodynamical simulations Horizon-AGN, Illustris-TNG-100 and EAGLE. We find that the most significant discrepancy is at the high mass end, where the simulations generally show that satellites have a higher contribution to the total stellar mass budget than the observations. This, together with the finding that the fraction of satellites is higher in the simulations, indicates that the feedback mechanisms acting in group-and cluster-scale halos appear to be less efficient in quenching the mass assembly of satellites, and/or that quenching occurs much later in the simulations.

In our hierarchical structure-formation paradigm, the observed morphological evolution of massive galaxies -- from rotationally-supported discs to dispersion-dominated spheroids -- is largely explained via galaxy merging. However, since mergers are likely to destroy discs, and the most massive galaxies have the richest merger histories, it is surprising that any discs exist at all at the highest stellar masses. Recent theoretical work by our group has used a cosmological, hydrodynamical simulation to suggest that extremely massive (M* > 10^11.4 MSun) discs form primarily via minor mergers between spheroids and gas-rich satellites, which create new rotational stellar components and leave discs as remnants. Here, we use UV-optical and HI data of massive galaxies, from the SDSS, GALEX, DECaLS and ALFALFA surveys, to test these theoretical predictions. Observed massive discs account for ~13% of massive galaxies, in good agreement with theory (~11%). ~64% of the observed massive discs exhibit tidal features, which are likely to indicate recent minor mergers, in the deep DECaLS images (compared to ~60% in their simulated counterparts). The incidence of these features is at least four times higher than in low-mass discs, suggesting that, as predicted, minor mergers play a significant (and outsized) role in the formation of these systems. The empirical star-formation rates agree well with theoretical predictions and, for a small galaxy sample with HI detections, the HI masses and fractions are consistent with the range predicted by the simulation. The good agreement between theory and observations indicates that extremely massive discs are indeed remnants of recent minor mergers between spheroids and gas-rich satellites., Comment: Accepted in MNRAS, 9 pages, 7 figures

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., Comment: 25 pages, 13 figures, 5 tables, accepted by Monthly Notices of Royal Astronomical Society (MNRAS) on 16/02/22

The excitation of the filamentary gas structures surrounding giant elliptical galaxies at the center of cool-core clusters, a.k.a BCGs (brightest cluster galaxies), is key to our understanding of active galactic nucleus feedback, and of the impact of environmental and local effects on star formation. We investigate the contribution of the thermal radiation from the cooling flow surrounding BCGs to the excitation of the filaments. We explore the effects of small levels of extra-heating (turbulence), and of metallicity, on the optical and infrared lines. Using the Cloudy code, we model the photoionization and photodissociation of a slab of gas of optical depth AV{\leq}30mag at constant pressure, in order to calculate self-consistently all of the gas phases, from ionized gas to molecular gas. The ionizing source is the EUV and soft X-ray radiation emitted by the cooling gas. We test these models comparing their predictions to the rich multi-wavelength observations, from optical to submillimeter. These models reproduce most of the multi-wavelength spectra observed in the nebulae surrounding the BCGs, not only the LINER-like optical diagnostics: [O iii]{\lambda} 5007 {\AA}/H\b{eta}, [N ii]{\lambda} 6583 {\AA}/H{\alpha} and ([S ii]{\lambda} 6716 {\AA}+[S ii]{\lambda} 6731 {\AA})/H{\alpha} but also the infrared emission lines from the atomic gas. The modeled ro-vib H2 lines also match observations, which indicates that near and mid-IR H2 lines are mostly excited by collisions between H2 molecules and secondary electrons produced naturally inside the cloud by the interaction between the X-rays and the cold gas in the filament. However, there is still some tension between ionized and molecular line tracers (i.e. CO), which requires to optimize the cloud structure and the density of the molecular zone., Comment: 16 pages, 9 figures, 3 tables. Accepted for publication in A&A

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, Comment: 18 pages, 17 figures, Accepted in A&A

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., Comment: Accepted for publication in MNRAS, 13 pages, 11 figures

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 \textsc{NewHorizon} 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., Comment: 22 pages, 19 figures, version accepted for publication in MNRAS. Some minor changes and additional appendices

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., Comment: 16 pages, 15 figures, Accepted for publication in MNRAS

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.

In an effort to better understand the formation of galaxy groups, we examine the kinematics of a large sample of spectroscopically confirmed X-ray galaxy groups in the Cosmic Evolution Survey (COSMOS) with a high sampling of galaxy group members up to $z=1$. We compare our results with predictions from the cosmological hydrodynamical simulation of {\sc Horizon-AGN}. Using a phase-space analysis of dynamics of groups with halo masses of $M_{\mathrm{200c}}\sim 10^{12.6}-10^{14.50}M_\odot$, we show that the brightest group galaxies (BGG) in low mass galaxy groups ($M_{\mathrm{200c}}<2 \times 10^{13} M_\odot$) have larger proper motions relative to the group velocity dispersion than high mass groups. The dispersion in the ratio of the BGG proper velocity to the velocity dispersion of the group, $\sigma_{\mathrm{BGG}}/\sigma_{group}$, is on average $1.48 \pm 0.13$ for low mass groups and $1.01 \pm 0.09$ for high mass groups. A comparative analysis of the {\sc Horizon-AGN} simulation reveals a similar increase in the spread of peculiar velocities of BGGs with decreasing group mass, though consistency in the amplitude, shape, and mode of the BGG peculiar velocity distribution is only achieved for high mass groups. The groups hosting a BGG with a large peculiar velocity are more likely to be offset from the $L_x-\sigma_{v}$ relation; this is probably because the peculiar motion of the BGG is influenced by the accretion of new members., Comment: 9 pages, 5 figures

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], Comment: 28 pages, accepted to A&A; language edited

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., Comment: 23 pages, 18 figures, accepted in MNRAS

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., Comment: Accepted for publication in MNRAS, 12 pages, 6 figures

Extended filamentary H$\alpha$ emission nebulae are a striking feature of nearby galaxy clusters but the formation mechanism of the filaments, and the processes which shape their morphology remain unclear. We conduct an investigation into the formation, evolution and destruction of dense gas in the center of a simulated, Perseus-like, cluster under the influence of a spin-driven jet. We particularly study the role played by condensation of dense gas from the diffuse intracluster medium, and the impact of direct uplifting of existing dense gas by the jets, in determining the spatial distribution and kinematics of the dense gas. We present a hydrodynamical simulation of an idealised Perseus-like cluster using the adaptive mesh refinement code {\sc ramses}. Our simulation includes a supermassive black hole (SMBH) that self-consistently tracks its spin evolution via its local accretion, and in turn drives a large-scale jet whose direction is based on the black hole's spin evolution. We show that the formation and destruction of dense gas is closely linked to the SMBH's feedback cycle, and that its morphology is highly variable throughout the simulation. While extended filamentary structures readily condense from the hot intra-cluster medium, they are easily shattered into an overly clumpy distribution of gas during their interaction with the jet driven outflows. Condensation occurs predominantly onto infalling gas located 5 - 15 kpc from the center during quiescent phases of the central AGN, when the local ratio of the cooling time to free fall time falls below 20, i.e. when $t_{\rm cool}/t_{\rm ff} < 20$. We find evidence for both condensation and uplifting of dense gas, but caution that purely hydrodynamical simulations struggle to effectively regulate the cluster cooling cycle and produce overly clumpy distributions of dense gas morphologies, compared to observation., Comment: 18 pages, 17 Figures. Accepted for publication in A&A

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., Comment: 14 pages, 10 figures, 1 table; submitted to MNRAS; Comments are welcome!

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., Comment: Submitted to A&A, 14 pages, comments welcome

Using the ligthcone from the cosmological hydrodynamical simulation Horizon-AGN, we produced a photometric catalogue over $0

Multi-phase filamentary structures around Brightest Cluster Galaxies are likely a key step of AGN-feedback. We observed molecular gas in 3 cool cluster cores: Centaurus, Abell S1101, and RXJ1539.5 and gathered ALMA and MUSE data for 12 other clusters. Those observations show clumpy, massive and long, 3--25 kpc, molecular filaments, preferentially located around the radio bubbles inflated by the AGN (Active Galactic Nucleus). Two objects show nuclear molecular disks. The optical nebula is certainly tracing the warm envelopes of cold molecular filaments. Surprisingly, the radial profile of the H$\alpha$/CO flux ratio is roughly constant for most of the objects, suggesting that (i) between 1.2 to 7 times more cold gas could be present and (ii) local processes must be responsible for the excitation. Projected velocities are between 100--400 km s$^{-1}$, with disturbed kinematics and sometimes coherent gradients. This is likely due to the mixing in projection of several thin unresolved filaments. The velocity fields may be stirred by turbulence induced by bubbles, jets or merger-induced sloshing. Velocity and dispersions are low, below the escape velocity. Cold clouds should eventually fall back and fuel the AGN. We compare the filament's radial extent, r$_{fil}$, with the region where the X-ray gas can become thermally unstable. The filaments are always inside the low-entropy and short cooling time region, where t$_{cool}$/t$_{ff}$<20 (9 of 13 sources). The range t$_{cool}$/t$_{ff}$, 8-23 at r$_{fil}$, is likely due to (i) a more complex gravitational potential affecting the free-fall time (e.g., sloshing, mergers); (ii) the presence of inhomogeneities or uplifted gas in the ICM, affecting the cooling time. For some of the sources, r$_{fil}$ lies where the ratio of the cooling time to the eddy-turnover time, t$_{cool}$/t$_{eddy}$, is approximately unity., Comment: 33 Pages, 17 Figures, submitted to A&A

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

The strikingly anisotropic large-scale distribution of matter made of an extended network of voids delimited by sheets, themselves segmented by filaments, within which matter flows towards compact nodes where they intersect, imprints its geometry on the dynamics of cosmic flows, ultimately shaping the distribution of galaxies and the redshift evolution of their properties. The (filament-type) saddle points of this cosmic web provide a local frame in which to quantify the induced physical and morphological evolution of galaxies on large scales. The properties of virtual galaxies within the Horizon-AGN simulation are stacked in such a frame. The iso-contours of the galactic number density, mass, specific star formation rate (sSFR), kinematics and age are clearly aligned with the filament axis with steep gradients perpendicular to the filaments. A comparison to a simulation without feedback from active galactic nuclei (AGN) illustrates its impact on quenching star formation of centrals away from the saddles. The redshift evolution of the properties of galaxies and their age distribution are consistent with the geometry of the bulk flow within that frame. They compare well with expectations from constrained Gaussian random fields and the scaling with the mass of non-linearity, modulo the redshift dependent impact of feedback processes. Physical properties such as sSFR and kinematics seem not to depend only on mean halo mass and density: the residuals trace the geometry of the saddle, which could point to other environment-sensitive physical processes, such as spin advection, and AGN feedback at high mass., Comment: 28 pages, 27 figures, submitted to MNRAS

This work investigates the alignment of galactic spins with the cosmic web across cosmic time using the cosmological hydrodynamical simulation Horizon-AGN. The cosmic web structure is extracted via the persistent skeleton as implemented in the DISPERSE algorithm. It is found that the spin of low-mass galaxies is more likely to be aligned with the filaments of the cosmic web and to lie within the plane of the walls while more massive galaxies tend to have a spin perpendicular to the axis of the filaments and to the walls. The mass transition is detected with a significance of 9 sigmas. This galactic alignment is consistent with the alignment of the spin of dark haloes found in pure dark matter simulations and with predictions from (anisotropic) tidal torque theory. However, unlike haloes, the alignment of low-mass galaxies is weak and disappears at low redshifts while the orthogonal spin orientation of massive galaxies is strong and increases with time, probably as a result of mergers. At fixed mass, alignments are correlated with galaxy morphology: the high-redshift alignment is dominated by spiral galaxies while elliptical centrals are mainly responsible for the perpendicular signal. These predictions for spin alignments with respect to cosmic filaments and unprecendently walls are successfully compared with existing observations. The alignment of the shape of galaxies with the different components of the cosmic web is also investigated. A coherent and stronger signal is found in terms of shape at high mass. The two regimes probed in this work induce competing galactic alignment signals for weak lensing, with opposite redshift and luminosity evolution. Understanding the details of these intrinsic alignments will be key to exploit future major cosmic shear surveys like Euclid or LSST., Comment: accepted for publication in MNRAS

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

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

Understanding the processes that drive the formation of black holes (BHs) is a key topic in observational cosmology. While the observed $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ correlation in bulge-dominated galaxies is thought to be produced by major mergers, the existence of a $M_{\mathrm{BH}}$--$M_{\star}$ relation, across all galaxy morphological types, suggests that BHs may be largely built by secular processes. Recent evidence that bulge-less galaxies, which are unlikely to have had significant mergers, are offset from the $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ relation, but lie on the $M_{\mathrm{BH}}$--$M_{\star}$ relation, has strengthened this hypothesis. Nevertheless, the small size and heterogeneity of current datasets, coupled with the difficulty in measuring precise BH masses, makes it challenging to address this issue using empirical studies alone. Here, we use Horizon-AGN, a cosmological hydrodynamical simulation to probe the role of mergers in BH growth over cosmic time. We show that (1) as suggested by observations, simulated bulge-less galaxies lie offset from the main $M_{\mathrm{BH}}$--$M_{\mathrm{Bulge}}$ relation, but on the $M_{\mathrm{BH}}$--$M_{\star}$ relation, (2) the positions of galaxies on the $M_{\mathrm{BH}}$--$M_{\star}$ relation are not affected by their merger histories and (3) only $\sim$35 per cent of the BH mass in today's massive galaxies is directly attributable to merging -- the majority ($\sim$65 per cent) of BH growth, therefore, takes place gradually, via secular processes, over cosmic time., Comment: 12 pages, 9 figures, accepted for publication in MNRAS

We quantify the quenching impact of the group environment using the spectroscopic survey Galaxy and Mass Assembly (GAMA) to z=0.2. The fraction of red (quiescent) galaxies, whether in groups or isolated, increases with both stellar mass and large-scale (5 Mpc) density. At fixed stellar mass, the red fraction is on average higher for satellites of red centrals than of blue (star-forming) centrals, a galactic conformity effect that increases with density. Most of the signal originates from groups that have the highest stellar mass, reside in the densest environments, and have massive, red only centrals. Assuming a color-dependent halo-to-stellar-mass ratio, whereby red central galaxies inhabit significantly more massive halos than blue ones of the same stellar mass, two regimes emerge more distinctly: at log(Mhalo/Msol) < 13, central quenching is still ongoing, conformity is no longer existent, and satellites and group centrals exhibit the same quenching excess over field galaxies at all mass and density, in agreement with the concept of "group quenching"; at log(Mhalo/Msol) > 13, a cutoff that sets apart massive (log(M*/Msol) > 11), fully quenched group centrals, conformity is meaningless, and satellites undergo significantly more quenching than their counterparts in smaller halos. The latter effect strongly increases with density, giving rise to the density-dependent conformity signal when both regimes are mixed. The star-formation of blue satellites in massive halos is also suppressed compared to blue field galaxies, while blue group centrals and the majority of blue satellites, which reside in low mass halos, show no deviation from the color-stellar mass relation of blue field galaxies., Comment: 23 pages, 22 figures, published version with changes to text and figures (results unchanged)

Feedback by Active Galactic Nuclei is often divided into quasar and radio mode, powered by radiation or radio jets, respectively. Both are fundamental in galaxy evolution, especially in late-type galaxies, as shown by cosmological simulations and observations of jet-ISM interactions in these systems. We compare AGN feedback by radiation and by collimated jets through a suite of simulations, in which a central AGN interacts with a clumpy, fractal galactic disc. We test AGN of $10^{43}$ and $10^{46}$ erg/s, considering jets perpendicular or parallel to the disc. Mechanical jets drive the more powerful outflows, exhibiting stronger mass and momentum coupling with the dense gas, while radiation heats and rarifies the gas more. Radiation and perpendicular jets evolve to be quite similar in outflow properties and effect on the cold ISM, while inclined jets interact more efficiently with all the disc gas, removing the densest $20\%$ in $20$ Myr, and thereby reducing the amount of cold gas available for star formation. All simulations show small-scale inflows of $0.01-0.1$ M$_\odot$/yr, which can easily reach down to the Bondi radius of the central supermassive black hole (especially for radiation and perpendicular jets), implying that AGN modulate their own duty cycle in a feedback/feeding cycle., Comment: 21 pages, 15 figures, 2 tables

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

The role of the cosmic web in shaping galaxy properties is investigated in the GAMA spectroscopic survey in the redshift range $0.03 \leq z \leq 0.25$. The stellar mass, $u - r$ dust corrected colour and specific star formation rate (sSFR) of galaxies are analysed as a function of their distances to the 3D cosmic web features, such as nodes, filaments and walls, as reconstructed by DisPerSE. Significant mass and type/colour gradients are found for the whole population, with more massive and/or passive galaxies being located closer to the filament and wall than their less massive and/or star-forming counterparts. Mass segregation persists among the star-forming population alone. The red fraction of galaxies increases when closing in on nodes, and on filaments regardless of the distance to nodes. Similarly, the star-forming population reddens (or lowers its sSFR) at fixed mass when closing in on filament, implying that some quenching takes place. Comparable trends are also found in the state-of-the-art hydrodynamical simulation Horizon-AGN. These results suggest that on top of stellar mass and large-scale density, the traceless component of the tides from the anisotropic large-scale environment also shapes galactic properties. An extension of excursion theory accounting for filamentary tides provides a qualitative explanation in terms of anisotropic assembly bias: at a given mass, the accretion rate varies with the orientation and distance to filaments. It also explains the absence of type/colour gradients in the data on smaller, non-linear scales., Comment: 26 pages, 18 figures, 7 tables, accepted for publication in MNRAS

The mass, accretion rate and formation time of dark matter haloes near proto-filaments (identified as saddle points of the potential) are analytically predicted using a conditional version of the excursion set approach in its so-called "upcrossing" approximation. The model predicts that at fixed mass, mass accretion rate and formation time vary with orientation and distance from the saddle, demonstrating that assembly bias is indeed influenced by the tides imposed by the cosmic web. Starved, early forming haloes of smaller mass lie preferentially along the main axis of filaments, while more massive and younger haloes are found closer to the nodes. Distinct gradients for distinct tracers such as typical mass and accretion rate occur because the saddle condition is anisotropic, and because the statistics of these observables depend on both the conditional means and their covariances. The theory is extended to other critical points of the potential field. The response of the mass function to variations of the matter density field (the so-called large scale bias) is computed, and its trend with accretion rate is shown to invert along the filament. The signature of this model should correspond at low redshift to an excess of reddened galactic hosts at fixed mass along preferred directions, as recently reported in spectroscopic and photometric surveys and in hydrodynamical simulations. The anisotropy of the cosmic web emerges therefore as a significant ingredient to describe jointly the dynamics and physics of galaxies, e.g. in the context of intrinsic alignments or morphological diversity., Comment: 30 pages, 23 figures, Submitted to MNRAS