Showing 1–50 of 115 results for author: Slyz, A

Inferring the Ionizing Photon Contributions of High-Redshift Galaxies to Reionization with JWST NIRCam Photometry

Authors: Nicholas Choustikov, Richard Stiskalek, Aayush Saxena, Harley Katz, Julien Devrient, Adrianne Slyz

Abstract: JWST observations have the potential to provide unprecedented constraints on the history of reionization and the sources responsible for the ionizing photons due to the detection of large populations of faint galaxies at $z\gg6$. Modelling reionization requires knowing both the number of ionizing photons that are produced by galaxies and the fraction of those photons that escape into the intergala… ▽ More JWST observations have the potential to provide unprecedented constraints on the history of reionization and the sources responsible for the ionizing photons due to the detection of large populations of faint galaxies at $z\gg6$. Modelling reionization requires knowing both the number of ionizing photons that are produced by galaxies and the fraction of those photons that escape into the intergalactic medium. Observational estimates of these values generally rely on spectroscopy for which large samples with well-defined selection functions remain limited. To overcome this challenge, we present an implicit likelihood inference (ILI) pipeline trained on mock photometry to predict the escaped ionizing luminosity of individual galaxies ($\dot{n}_{\rm ion}$) based on photometric magnitudes and redshifts. Compared to traditional SED-fitting methods, the new ILI pipeline is consistently more accurate and significantly faster. We deploy the method on a sample of 4,559 high-redshift galaxies from the JADES Deep survey, finding a gentle redshift evolution of $\log_{10}(\dot{n}_{\rm ion}) = (0.08\pm0.01)z + (51.60\pm0.06)$, with late-time values for $\dot{N}_{\rm ion}$ consistent with theoretical models and observations. We measure the evolution of the volume-averaged ionized fraction and optical depth to find that observed populations of star-forming galaxies are capable of driving reionization to completion at $z\sim 5.3$ without the need for AGN or other exotic sources. The $20\%$ of UV-brightest galaxies ($M_{\rm UV}<-18.5$) in our sample can reionize only $\sim30\%$ of the survey volume, demonstrating that faint LyC emitters are crucial for reionization. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 17 pages, 11 figures, 2 tables, submitted to MNRAS

Dwarf galaxies as a probe of a primordially magnetized Universe

Authors: Mahsa Sanati, Sergio Martin-Alvarez, Jennifer Schober, Yves Revaz, Adrianne Slyz, Julien Devriendt

Abstract: The true nature of primordial magnetic fields (PMFs) and their role in the formation of galaxies still remains elusive. To shed light on these unknowns, we investigate their impact by varying two sets of properties: (i) accounting for the effect of PMFs on the initial matter power spectrum, and (ii) accounting for their magneto-hydrodynamical effects on the formation of galaxies. By comparing both… ▽ More The true nature of primordial magnetic fields (PMFs) and their role in the formation of galaxies still remains elusive. To shed light on these unknowns, we investigate their impact by varying two sets of properties: (i) accounting for the effect of PMFs on the initial matter power spectrum, and (ii) accounting for their magneto-hydrodynamical effects on the formation of galaxies. By comparing both we can determine the dominant agent in shaping galaxy evolution. We use the magneto-hydrodynamics code RAMSES, to generate multiple zoom-in simulations for eight different host halos of dwarf galaxies across a wide luminosity range of $10^3-10^6\,L_{\odot}$. We explore a variety of primordial magnetic field (comoving) strengths ranging from $0.05$ to $0.50\,\mathrm{nG}$. We find magnetic fields in the interstellar medium not only modify star formation in dwarf spheroidal galaxies but also completely prevent the formation of stars in less compact ultra-faints with halo mass and stellar mass below $\sim 2.5\cdot10^9$ and $3\cdot10^6\,M_{\odot}$, respectively. At high redshifts, the impact of PMFs on host halos of dwarf galaxies through the modification of the matter power spectrum is more dominant than the influence of magneto-hydrodynamics in shaping their gaseous structure. Through the amplification of small perturbations ranging in mass from $10^7$ to $10^9\,M_{\odot}$ in the $Λ$CDM$+$PMFs matter power spectrum, primordial fields expedite the formation of the first dark matter halos, leading to an earlier onset and a higher star formation rate at redshifts $z>12$. △ Less

Submitted 8 March, 2024; originally announced March 2024.

The Great Escape: Understanding the Connection Between Ly$α$ Emission and LyC Escape in Simulated JWST Analogues

Authors: Nicholas Choustikov, Harley Katz, Aayush Saxena, Thibault Garel, Julien Devriendt, Adrianne Slyz, Taysun Kimm, Jeremy Blaizot, Joki Rosdahl

Abstract: Constraining the escape fraction of Lyman Continuum (LyC) photons from high-redshift galaxies is crucial to understanding reionization. Recent observations have demonstrated that various characteristics of the Ly$α$ emission line correlate with the inferred LyC escape fraction ($f_{\rm esc}^{\rm LyC}$) of low-redshift galaxies. Using a data-set of 9,600 mock Ly$α$ spectra of star-forming galaxies… ▽ More Constraining the escape fraction of Lyman Continuum (LyC) photons from high-redshift galaxies is crucial to understanding reionization. Recent observations have demonstrated that various characteristics of the Ly$α$ emission line correlate with the inferred LyC escape fraction ($f_{\rm esc}^{\rm LyC}$) of low-redshift galaxies. Using a data-set of 9,600 mock Ly$α$ spectra of star-forming galaxies at $4.64 \leq z \leq 6$ from the SPHINX$^{20}$ cosmological radiation hydrodynamical simulation, we study the escape of Ly$α$ and LyC photons. We find that our mock Ly$α$ observations are representative of high-redshift observations and that typical observational methods tend to over-predict the Ly$α$ escape fraction ($f_{\rm esc}^{\rm Lyα}$) by as much as two dex. We investigate the correlations between $f_{\rm esc}^{\rm LyC}$ and $f_{\rm esc}^{\rm Lyα}$, Ly$α$ equivalent width ($W_λ({\rm Lyα})$), peak separation ($v_{\rm sep}$), central escape fraction ($f_{\rm cen}$), and red peak asymmetry ($A_f^{\rm red}$). We find that $f_{\rm esc}^{\rm Lyα}$ and $f_{\rm cen}$ are good diagnostics for LyC leakage, selecting for galaxies with lower neutral gas densities and less UV attenuation that have recently experienced supernova feedback. In contrast, $W_λ({\rm Lyα})$ and $v_{\rm sep}$ are found to be necessary but insufficient diagnostics, while $A_f^{\rm red}$ carries little information. Finally, we use stacks of Ly$α$, H$α$, and F150W mock surface brightness profiles to find that galaxies with high $f_{\rm esc}^{\rm LyC}$ have less extended Ly$α$ and F150W but larger H$α$ haloes than their non-leaking counterparts. This confirms that Ly$α$ spectral profiles and surface brightness morphology can be used to better understand the escape of LyC photons from galaxies during the Epoch of Reionization. △ Less

Submitted 24 June, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Comments: 22 pages, 18+3 figures, accepted by MNRAS. Accepted version includes a new appendix

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). VII. A Tomographic View of Far-infrared and Radio Polarimetric Observations through MHD Simulations of Galaxies

Authors: Sergio Martin-Alvarez, Enrique Lopez-Rodriguez, Tara Dacunha, Susan E. Clark, Alejandro S. Borlaff, Rainer Beck, Francisco Rodríguez Montero, S. Lyla Jung, Julien Devriendt, Adrianne Slyz, Julia Roman-Duval, Evangelia Ntormousi, Mehrnoosh Tahani, Kandaswamy Subramanian, Daniel A. Dale, Pamela M. Marcum, Konstantinos Tassis, Ignacio del Moral-Castro, Le Ngoc Tram, Matt J. Jarvis

Abstract: The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy t… ▽ More The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of $h_\text{I radio} \sim h_\text{PI radio} = 0.23 \pm 0.03$ kpc) relative to the total FIR and polarized emission that are concentrated in the disk midplane ($h_\text{I FIR} \sim h_\text{PI FIR} = 0.11 \pm 0.01$ kpc). Radio emission traces magnetic fields at scales of $\gtrsim 300$ pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution ($<300$ pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism. △ Less

Submitted 24 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Accepted for publication in the ApJ. Replaced to match final version. 35 pages, 17 figures, 4 tables

Evaluating the reconstruction of individual haloes in constrained cosmological simulations

Authors: Richard Stiskalek, Harry Desmond, Julien Devriendt, Adrianne Slyz

Abstract: Constrained cosmological simulations play an important role in modelling the local Universe, enabling investigation of the dark matter content of local structures and their formation histories. We introduce a method for determining the extent to which individual haloes are reliably reconstructed between constrained simulations, and apply it to the Constrained Simulations in BORG (CSiBORG) suite of… ▽ More Constrained cosmological simulations play an important role in modelling the local Universe, enabling investigation of the dark matter content of local structures and their formation histories. We introduce a method for determining the extent to which individual haloes are reliably reconstructed between constrained simulations, and apply it to the Constrained Simulations in BORG (CSiBORG) suite of $101$ high-resolution realisations across the posterior probability distribution of initial conditions from the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. The method is based on the overlap of the initial Lagrangian patch of a halo in one simulation with those in another, and therefore measures the degree to which the haloes' particles are initially coincident. By this metric we find consistent reconstructions of $M\gtrsim10^{14}~M_\odot / h$ haloes across the CSiBORG simulations, indicating that the constraints from the BORG algorithm are sufficient to pin down the masses, positions and peculiar velocities of clusters to high precision. The effect of the constraints tapers off towards lower mass however, and the halo spins and concentrations are largely unconstrained at all masses. We document the advantages of evaluating halo consistency in the initial conditions, describe how the method may be used to quantify our knowledge of the halo field given galaxy survey data analysed through the lens of probabilistic inference machines such as BORG, and describe applications to matched but unconstrained simulations. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Comments: 11 pages, 11 figures. To be submitted to MNRAS. Comments are welcome!

The formation of cores in galaxies across cosmic time -- the existence of cores is not in tension with the LCDM paradigm

Authors: 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

Abstract: 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),… ▽ More 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. △ Less

Submitted 8 January, 2024; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 13 pages, 11 figures, Accepted to MNRAS

The impact of cosmic rays on the interstellar medium and galactic outflows of Milky Way analogues

Authors: Francisco Rodríguez Montero, Sergio Martin-Alvarez, Adrianne Slyz, Julien Devriendt, Yohan Dubois, Debora Sijacki

Abstract: During the last decade, cosmological simulations have managed to reproduce realistic and morphologically diverse galaxies, spanning the Hubble sequence. Central to this success was a phenomenological calibration of the few included feedback processes, whilst glossing over higher complexity baryonic physics. This approach diminishes the predictive power of such simulations, preventing to further ou… ▽ More During the last decade, cosmological simulations have managed to reproduce realistic and morphologically diverse galaxies, spanning the Hubble sequence. Central to this success was a phenomenological calibration of the few included feedback processes, whilst glossing over higher complexity baryonic physics. This approach diminishes the predictive power of such simulations, preventing to further our understanding of galaxy formation. To tackle this fundamental issue, we investigate the impact of cosmic rays (CRs) and magnetic fields on the interstellar medium (ISM) and the launching of outflows in a cosmological zoom-in simulation of a Milky Way-like galaxy. We find that including CRs decreases the stellar mass of the galaxy by a factor of 10 at high redshift and $\sim 4$ at cosmic noon, leading to a stellar mass to halo mass ratio in good agreement with abundance matching models. Such decrease is caused by two effects: i) a reduction of cold, high-density, star-forming gas, and ii) a larger fraction of SN events exploding at lower densities, where they have a higher impact. SN-injected CRs produce enhanced, multi-phase galactic outflows, which are accelerated by CR pressure gradients in the circumgalactic medium of the galaxy. While the mass budget of these outflows is dominated by the warm ionised gas, warm neutral and cold gas phases contribute significantly at high redshifts. Importantly, our work shows that future JWST observations of galaxies and their multi-phase outflows across cosmic time have the ability to constrain the role of CRs in regulating star formation. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 20 pages, 11 figures, submitted MNRAS

The Physics of Indirect Estimators of Lyman Continuum Escape and their Application to High-Redshift JWST Galaxies

Authors: Nicholas Choustikov, Harley Katz, Aayush Saxena, Alex J. Cameron, Julien Devriendt, Adrianne Slyz, Joki Rosdahl, Jeremy Blaizot, Leo Michel-Dansac

Abstract: Reliable indirect diagnostics of LyC photon escape from galaxies are required to understand which sources were the dominant contributors to reionization. While multiple LyC escape fraction ($f_{\rm esc}$) indicators have been proposed to trace favourable conditions for LyC leakage from the interstellar medium of low-redshift ''analog'' galaxies, it remains unclear whether these are applicable at h… ▽ More Reliable indirect diagnostics of LyC photon escape from galaxies are required to understand which sources were the dominant contributors to reionization. While multiple LyC escape fraction ($f_{\rm esc}$) indicators have been proposed to trace favourable conditions for LyC leakage from the interstellar medium of low-redshift ''analog'' galaxies, it remains unclear whether these are applicable at high redshifts where LyC emission cannot be directly observed. Using a library of 14,120 mock spectra of star-forming galaxies with redshifts $4.64 \leq z \leq 10$ from the SPHINX$^{20}$ cosmological radiation hydrodynamics simulation, we develop a framework for the physics that leads to high $f_{\rm esc}$. We investigate LyC leakage from our galaxies based on the criteria that successful LyC escape diagnostics must i) track a high specific star formation rate, ii) be sensitive to stellar population age in the range $3.5-10~$Myr representing the times when supernova first explode to when LyC production significantly drops, and iii) include a proxy for neutral gas content and gas density in the interstellar medium. ${\rm O}_{32}$, $Σ_{\rm SFR}$, M$_{\rm UV}$, and H$β$ equivalent width select for one or fewer of our criteria, rendering them either necessary but insufficient or generally poor diagnostics. In contrast, UV slope ($β$), and ${\rm E(B-V)}$ match two or more of our criteria, rendering them good $f_{\rm esc}$ diagnostics (albeit with significant scatter). Using our library, we build a quantitative model for predicting $f_{\rm esc}$ based on direct observables. When applied to bright $z > 6$ Ly$α$ emitters observed with JWST, we find that the majority of them have $f_{\rm esc} \lesssim 10\%$. △ Less

Submitted 13 March, 2024; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: 16 pages, 15 figures, 1 table, amended to reflect version accepted to MNRAS

Boosting galactic outflows with enhanced resolution

Authors: Martin P. Rey, Harley B. Katz, Alex J. Cameron, Julien Devriendt, Adrianne Slyz

Abstract: We study how better resolving the cooling length of galactic outflows affect their energetics. We perform radiative-hydrodynamical galaxy formation simulations of an isolated dwarf galaxy ($M_{\star}=10^{8}\, \mathrm{M}_\odot$) with the Ramses-RTZ code, accounting for non-equilibrium cooling and chemistry coupled to radiative transfer. Our simulations reach a spatial resolution of… ▽ More We study how better resolving the cooling length of galactic outflows affect their energetics. We perform radiative-hydrodynamical galaxy formation simulations of an isolated dwarf galaxy ($M_{\star}=10^{8}\, \mathrm{M}_\odot$) with the Ramses-RTZ code, accounting for non-equilibrium cooling and chemistry coupled to radiative transfer. Our simulations reach a spatial resolution of $18 \, \mathrm{pc}$ in the interstellar medium (ISM) using a traditional quasi-Lagrangian scheme. We further implement a new adaptive mesh refinement (AMR) strategy to resolve the local gas cooling length, allowing us to gradually increase the resolution in the stellar-feedback-powered outflows, from $\geq 200 \, \mathrm{pc}$ to $18 \, \mathrm{pc}$. The propagation of outflows into the inner circumgalactic medium (CGM) is significantly modified by this additional resolution, but the ISM, star formation and feedback remain by and large the same. With increasing resolution in the diffuse gas, the hot outflowing phase ($T > 8 \times 10^{4} \, \mathrm{K}$) systematically reaches overall higher temperatures and stays hotter for longer as it propagates outwards. This leads to two-fold increases in the time-averaged mass and metal outflow loading factors away from the galaxy ($r=5\, \mathrm{kpc}$), a five-fold increase in the average energy loading factor, and a $\approx$50 per cent increase in the number of sightlines with $N_{\text{OVI}} \geq 10^{13}\, \mathrm{cm}^{-2}$. Such a significant boost to the energetics of outflows without new feedback mechanisms or channels strongly motivates future studies quantifying the efficiency with which better-resolved multiphase outflows regulate galactic star formation in a cosmological context. △ Less

Submitted 11 February, 2024; v1 submitted 16 February, 2023; originally announced February 2023.

Comments: Matching published version in MNRAS. Updated Fig.9 following bug correction in outflow loading factor normalisation, main conclusions unchanged. Additional robustness tests, explanations and visualisations throughout the paper following peer-review

PRISM: A Non-Equilibrium, Multiphase Interstellar Medium Model for Radiation Hydrodynamics Simulations of Galaxies

Authors: Harley Katz, Shenghua Liu, Taysun Kimm, Martin P. Rey, Eric P. Andersson, Alex J. Cameron, Francisco Rodriguez-Montero, Oscar Agertz, Julien Devriendt, Adrianne Slyz

Abstract: We introduce the PRISM interstellar medium (ISM) model for thermochemistry and its implementation in the RAMSES-RTZ code. The model includes a non-equilibrium primordial, metal, and molecular chemistry network for 115 species coupled to on-the-fly multifrequency radiation transport. PRISM accurately accounts for the dominant ISM cooling and heating processes in the low-density regime (i.e.… ▽ More We introduce the PRISM interstellar medium (ISM) model for thermochemistry and its implementation in the RAMSES-RTZ code. The model includes a non-equilibrium primordial, metal, and molecular chemistry network for 115 species coupled to on-the-fly multifrequency radiation transport. PRISM accurately accounts for the dominant ISM cooling and heating processes in the low-density regime (i.e. $ρ<10^5\ {\rm cm^{-3}}$), including photoheating, photoelectric heating, H$_2$ heating/cooling, cosmic-ray heating, H/He cooling, metal-line cooling, CO cooling, and dust cooling (recombination and gas-grain collisions). We validate the model by comparing 1D equilibrium simulations across six dex in metallicity to existing 1D ISM models in the literature. We apply PRISM to high-resolution (4.5 pc) isolated dwarf galaxy simulations that include state-of-the-art models for star formation and stellar feedback to take an inventory of which cooling and heating processes dominate each different gas phase of a galaxy and to understand the importance of non-equilibrium effects. We show that most of the ISM gas is either close to thermal equilibrium or exhibits a slight cooling instability, while from a chemical perspective, the non-equilibrium electron fraction is often more than three times higher or lower than the equilibrium value, which impacts cooling, heating, and observable emission lines. Electron enhancements are attributed to recombination lags while deficits are shown to be due to rapid cosmic-ray heating. The PRISM model and its coupling to RAMSES-RTZ is applicable to a wide variety of astrophysical scenarios, from cosmological simulations to isolated giant molecular clouds, and is particularly useful for understanding how changes to ISM physics impact observable quantities such as metallic emission lines. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 21 pages, 16 figures, submitted to MNRAS

Two Modes of LyC Escape From Bursty Star Formation: Implications for [C II] Deficits and the Sources of Reionization

Authors: Harley Katz, Aayush Saxena, Joki Rosdahl, Taysun Kimm, Jeremy Blaizot, Thibault Garel, Leo Michel-Dansac, Martin Haehnelt, Richard S. Ellis, Laura Penterrici, Julien Devriendt, Adrianne Slyz

Abstract: We use the SPHINX$^{20}$ cosmological radiation hydrodynamics simulation to study how Lyman Continuum (LyC) photons escape from galaxies and the observational signatures of this escape. We define two classes of LyC leaker: Bursty Leakers and Remnant Leakers, based on their star formation rates (SFRs) that are averaged over 10 Myr (SFR$_{10}$) or 100 Myr (SFR$_{100}$). Both have $f_{\rm esc}>20\%$… ▽ More We use the SPHINX$^{20}$ cosmological radiation hydrodynamics simulation to study how Lyman Continuum (LyC) photons escape from galaxies and the observational signatures of this escape. We define two classes of LyC leaker: Bursty Leakers and Remnant Leakers, based on their star formation rates (SFRs) that are averaged over 10 Myr (SFR$_{10}$) or 100 Myr (SFR$_{100}$). Both have $f_{\rm esc}>20\%$ and experienced an extreme burst of star formation, but Bursty Leakers have ${\rm SFR_{10}>SFR_{100}}$, while Remnant Leakers have ${\rm SFR_{10}<SFR_{100}}$. The maximum SFRs in these bursts were typically $\sim100$ times greater than the SFR of the galaxy prior to the burst, a rare $2σ$ outlier among the general high-redshift galaxy population. Bursty Leakers are qualitatively similar to ionization-bounded nebulae with holes, exhibiting high ionization parameters and typical HII region gas densities. Remnant Leakers show properties of density-bounded nebulae, having normal ionization parameters but much lower HII region densities. Both types of leaker exhibit [CII]$_{\rm 158μm}$ deficits on the [CII]-SFR$_{100}$ relation, while only Bursty Leakers show deficits when SFR$_{10}$ is used. We predict that [CII] luminosity and SFR indicators such as H$α$ and M$_{\rm 1500Å}$ can be combined to identify both types of LyC leaker and the mode by which photons are escaping. These predictions can be tested with [CII] observations of known $z=3-4$ LyC leakers. Finally, we show that leakers with $f_{\rm esc}>20\%$ dominate the ionizing photon budget at $z\gtrsim7.5$ but the contribution from galaxies with $f_{\rm esc}<5\%$ becomes significant at the tail-end of reionization. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Comments: 16 pages, 14 figures, accepted by MNRAS

The information on halo properties contained in spectroscopic observations of late-type galaxies

Authors: Tariq Yasin, Harry Desmond, Julien Devriendt, Adrianne Slyz

Abstract: Rotation curves are the key observational manifestation of the dark matter distribution around late-type galaxies. In a halo model context, the precision of constraints on halo parameters is a complex function of the properties of the measurements as well as properties of the galaxy itself. Forthcoming surveys will resolve rotation curves to varying degrees of precision, or measure their integrate… ▽ More Rotation curves are the key observational manifestation of the dark matter distribution around late-type galaxies. In a halo model context, the precision of constraints on halo parameters is a complex function of the properties of the measurements as well as properties of the galaxy itself. Forthcoming surveys will resolve rotation curves to varying degrees of precision, or measure their integrated effect in the HI linewidth. To ascertain the relative significance of the relevant quantities for constraining halo properties, we study the information on halo mass and concentration as quantified by the Kullback-Leibler divergence of the kinematics-informed posterior from the uninformative prior. We calculate this divergence as a function of the different types of spectroscopic observation, properties of the measurement, galaxy properties and auxiliary observational data on the baryonic components. Using the SPARC sample, we find that fits to the full rotation curve exhibit a large variation in information gain between galaxies, ranging from ~1 to ~11 bits. The variation is predominantly caused by the vast differences in the number of data points and the size of velocity uncertainties between the SPARC galaxies. We also study the relative importance of the minimum HI surface density probed and the size of velocity uncertainties on the constraining power on the inner halo slope, finding the latter to be significantly more important. We spell out the implications of these results for optimising the strategy of galaxy surveys aiming to constrain galaxies' dark matter distributions, highlighting spectroscopic precision as the most important factor. △ Less

Submitted 8 January, 2024; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: 13 pages, 7 figures. Revision to match MNRAS published version

Journal ref: MNRAS 2023 525 (4): 5066-5079

Impact of Radiation Feedback on the Formation of Globular Cluster Candidates during Cloud-Cloud Collisions

Authors: Daniel Han, Taysun Kimm, Harley Katz, Julien Devriendt, Adrianne Slyz

Abstract: To understand the impact of radiation feedback during the formation of a globular cluster (GC), we simulate a head-on collision of two turbulent giant molecular clouds (GMCs). A series of idealized radiation-hydrodynamic simulations is performed, with and without stellar radiation or Type II supernovae. We find that a gravitationally bound, compact star cluster of mass… ▽ More To understand the impact of radiation feedback during the formation of a globular cluster (GC), we simulate a head-on collision of two turbulent giant molecular clouds (GMCs). A series of idealized radiation-hydrodynamic simulations is performed, with and without stellar radiation or Type II supernovae. We find that a gravitationally bound, compact star cluster of mass $M_{\rm GC} \sim 10^5\,M_\odot$ forms within $\approx 3\,{\rm Myr}$ when two GMCs with mass $M_{\rm GMC}=3.6\times10^5\,M_\odot$ collide. The GC candidate does not form during a single collapsing event but emerges due to the mergers of local dense gas clumps and gas accretion. The momentum transfer due to the absorption of the ionizing radiation is the dominant feedback process that suppresses the gas collapse and photoionization becomes efficient once a sufficient number of stars form. The cluster mass is larger by a factor of $\sim 2$ when the radiation feedback is neglected, and the difference is slightly more pronounced (16\%) when extreme $\rm Lyα$ feedback is considered in the fiducial run. In the simulations with radiation feedback, supernovae explode after the star-forming clouds are dispersed, and their metal ejecta are not instantaneously recycled to form stars. △ Less

Submitted 12 July, 2022; originally announced July 2022.

Comments: 14 pages, 7 figures, Accepted for publication in the Astrophysical Journal

The Challenges of Identifying Population III Stars in the Early Universe

Authors: Harley Katz, Taysun Kimm, Richard S. Ellis, Julien Devriendt, Adrianne Slyz

Abstract: The recent launch of JWST has enabled the exciting prospect of detecting the first generation of metal-free, Population III (Pop. III) stars. Determining the emission line signatures that robustly signify the presence of Pop. III stars against other possible contaminants represents a key challenge for interpreting JWST data. To this end, we run high-resolution (sub-pc) cosmological radiation hydro… ▽ More The recent launch of JWST has enabled the exciting prospect of detecting the first generation of metal-free, Population III (Pop. III) stars. Determining the emission line signatures that robustly signify the presence of Pop. III stars against other possible contaminants represents a key challenge for interpreting JWST data. To this end, we run high-resolution (sub-pc) cosmological radiation hydrodynamics simulations of the region around a dwarf galaxy at $z\geq10$ to predict the emission line signatures of the Pop. III/Pop. II transition. We show that the absence of metal emission lines is a poor diagnostic of Pop. III stars because metal-enriched galaxies in our simulation can maintain low [OIII] 5007${\rm Å}$ emission that may be undetectable due to sensitivity limits. Combining spectral hardness probes (e.g. HeII 1640${\rm Å}$/H$α$) with metallicity diagnostics is more likely to probe the presence of metal-free stars, although contamination from Wolf-Rayet stars, X-ray binaries, or black holes may be important. The hard emission from Pop. III galaxies fades fast due to the short stellar lifetimes of massive Pop. III stars, which could further inhibit detection. Similarly, Pop. III stars may be detectable after they evolve off the main-sequence due to the cooling radiation from nebular gas or a supernova remnant; however, these signatures are also short-lived (i.e. few Myr), and contaminants such as flickering black holes might confuse this diagnostic. While JWST will provide a unique opportunity to spectroscopically probe the nature of the earliest galaxies, both the short timescales associated with pristine systems and ambiguities in interpreting key diagnostic emission lines may hinder progress. Special care will be needed before claiming the discovery of systems with pure Pop. III stars. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: 9 pages, 8 figures, submitted to MNRAS

Inferring dark matter halo properties for HI-selected galaxies

Authors: Tariq Yasin, Harry Desmond, Julien Devriendt, Adrianne Slyz

Abstract: We set constraints on the dark matter halo mass and concentration of ~22,000 individual galaxies visible both in HI (from the ALFALFA survey) and optical light (from the SDSS). This is achieved by combining two Bayesian models, one for the HI line width as a function of the stellar and neutral hydrogen mass distributions in a galaxy using kinematic modelling, and the other for the galaxy's total b… ▽ More We set constraints on the dark matter halo mass and concentration of ~22,000 individual galaxies visible both in HI (from the ALFALFA survey) and optical light (from the SDSS). This is achieved by combining two Bayesian models, one for the HI line width as a function of the stellar and neutral hydrogen mass distributions in a galaxy using kinematic modelling, and the other for the galaxy's total baryonic mass using the technique of inverse subhalo abundance matching. We hence quantify the constraining power on halo properties of spectroscopic and photometric observations, and assess their consistency. We find good agreement between the two sets of posteriors, although there is a sizeable population of low-line width galaxies that favour significantly smaller dynamical masses than expected from abundance matching (especially for cuspy halo profiles). Abundance matching provides significantly more stringent bounds on halo properties than the HI line width, even with a mass--concentration prior included, although combining the two provides a mean gain of 40% for the sample when fitting an NFW profile. We also use our kinematic posteriors to construct a baryonic mass--halo mass relation, which we find to be near power-law, and with a somewhat shallower slope than expected from abundance matching. Our method demonstrates the potential of combining photometric and spectroscopic observations to precisely map out the dark matter distribution at the galaxy scale using upcoming HI surveys such as the SKA. △ Less

Submitted 30 June, 2023; v1 submitted 30 June, 2022; originally announced June 2022.

Comments: 23 pages, 15 figures, accepted for publication in MNRAS

Journal ref: MNRAS 2023 526 (4): 5861-5882

MgII in the JWST Era: a Probe of Lyman Continuum Escape?

Authors: Harley Katz, Thibault Garel, Joakim Rosdahl, Valentin Mauerhofer, Taysun Kimm, Jérémy Blaizot, Léo Michel-Dansac, Julien Devriendt, Adrianne Slyz, Martin Haehnelt

Abstract: Limited constraints on the evolution of the Lyman Continuum (LyC) escape fraction represent one of the primary uncertainties in the theoretical determination of the reionization history. Due to the intervening intergalactic medium (IGM), the possibility of observing LyC photons directly in the epoch of reionization is highly unlikely. For this reason, multiple indirect probes of LyC escape have be… ▽ More Limited constraints on the evolution of the Lyman Continuum (LyC) escape fraction represent one of the primary uncertainties in the theoretical determination of the reionization history. Due to the intervening intergalactic medium (IGM), the possibility of observing LyC photons directly in the epoch of reionization is highly unlikely. For this reason, multiple indirect probes of LyC escape have been identified, some of which are used to identify low-redshift LyC leakers (e.g. O32), while others are primarily useful at $z>6$ (e.g. [OIII]/[CII] far infrared emission). The flux ratio of the resonant MgII doublet emission at 2796$\dot{\rm A}$ and 2803$\dot{\rm A}$ as well as the MgII optical depth have recently been proposed as ideal diagnostics of LyC leakage that can be employed at $z>6$ with JWST. Using state-of-the-art cosmological radiation hydrodynamics simulations post-processed with CLOUDY and resonant-line radiative transfer, we test whether MgII is indeed a useful probe of LyC leakage. Our simulations indicate that the majority of bright, star-forming galaxies with high LyC escape fractions are expected to be MgII emitters rather than absorbers at $z=6$. However, we find that the MgII doublet flux ratio is a more sensitive indicator of dust rather than neutral hydrogen, limiting its use as a LyC leakage indicator to only galaxies in the optically thin regime. Given its resonant nature, we show that MgII will be an exciting probe of the complex kinematics in high-redshift galaxies in upcoming JWST observations. △ Less

Submitted 23 May, 2022; originally announced May 2022.

Comments: 25 Pages, 20 Figures, MNRAS Accepted

On the Viability of Determining Galaxy Properties from Observations I: Star Formation Rates and Kinematics

Authors: Kearn Grisdale, Laurence Hogan, Dimitra Rigopoulou, Niranjan Thatte, Miguel Pereira-Santaella, Julien Devriendt, Adrianne Slyz, Ismael García-Bernete, Yohan Dubois, Sukyoung K. Yi, Katarina Kraljic

Abstract: We explore how observations relate to the physical properties of the emitting galaxies by post-processing a pair of merging $z\sim2$ galaxies from the cosmological, hydrodynamical simulation NewHorizon using LCARS (Light from Cloudy Added to RAMSES) to encode the physical properties of the simulated galaxy into H$α$ emission line. By carrying out mock observations and analysis on these data cubes… ▽ More We explore how observations relate to the physical properties of the emitting galaxies by post-processing a pair of merging $z\sim2$ galaxies from the cosmological, hydrodynamical simulation NewHorizon using LCARS (Light from Cloudy Added to RAMSES) to encode the physical properties of the simulated galaxy into H$α$ emission line. By carrying out mock observations and analysis on these data cubes we ascertain which physical properties of the galaxy will be recoverable with the HARMONI spectrograph on the European Extremely Large Telescope (ELT). We are able to estimate the galaxy's star formation rate and dynamical mass to a reasonable degree of accuracy, with values within a factor of $1.81$ and $1.38$ of the true value. The kinematic structure of the galaxy is also recovered in mock observations. Furthermore, we are able to recover radial profiles of the velocity dispersion and are therefore able to calculate how the dynamical ratio varies as a function of distance from the galaxy centre. Finally, we show that when calculated on galaxy scales the dynamical ratio does not always provide a reliable measure of a galaxy's stability against gravity or act as an indicator of a minor merger. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 20 pages, 17 figures, Accepted for Publication in MNRAS, Comments and questions welcome

The effect of local universe constraints on halo abundance and clustering

Authors: Maxwell L. Hutt, Harry Desmond, Julien Devriendt, Adrianne Slyz

Abstract: Cosmological $N$-body simulations of the dark matter component of the universe typically use initial conditions with a fixed power spectrum and random phases of the density field, leading to structure consistent with the local distribution of galaxies only in a statistical sense. It is, however, possible to infer the initial phases which lead to the configuration of galaxies and clusters that we s… ▽ More Cosmological $N$-body simulations of the dark matter component of the universe typically use initial conditions with a fixed power spectrum and random phases of the density field, leading to structure consistent with the local distribution of galaxies only in a statistical sense. It is, however, possible to infer the initial phases which lead to the configuration of galaxies and clusters that we see around us. We analyse the CSiBORG suite of 101 simulations, formed by constraining the density field within 155 Mpc$/h$ with dark matter particle mass $4.38\times10^9 M_\odot$, to quantify the degree to which constraints imposed on 2.65 Mpc$/h$ scales reduce variance in the halo mass function and halo-halo cross-correlation function on a range of scales. This is achieved by contrasting CSiBORG with a subset of the unconstrained Quijote simulations and expectations for the $Λ$CDM average. Using the FOF, PHEW and HOP halofinders, we show that the CSiBORG suite beats cosmic variance at large mass scales ($\gtrsim 10^{14}M_{\odot}/h$), which are most strongly constrained by the initial conditions, and exhibits a significant halo-halo cross-correlation out to $\sim30$ Mpc$/h$. Moreover, the effect of the constraints percolates down to lower mass objects and to scales below those on which they are imposed. Finally, we develop an algorithm to "twin" halos between realisations and show that approximately 50% of halos with mass greater than $10^{15}M_{\odot}/h$ can be identified in all realisations of the CSiBORG suite. We make the CSiBORG halo catalogues publicly available for future applications requiring knowledge of the local halo field. △ Less

Submitted 23 August, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 11 pages, 8 figures; minor revision, as accepted for publication by MNRAS

Simulating Jellyfish Galaxies: A Case Study for a Gas-Rich Dwarf Galaxy

Authors: Jaehyun Lee, Taysun Kimm, Jérémy Blaizot, Harley Katz, Wonki Lee, Yun-Kyeong Sheen, Julien Devriendt, Adrianne Slyz

Abstract: We investigate the formation of jellyfish galaxies using radiation-hydrodynamic simulations of gas-rich dwarf galaxies with a multi-phase interstellar medium (ISM). We find that the ram-pressure-stripped (RPS) ISM is the dominant source of molecular clumps in the near wake within 10 kpc from the galactic plane, while in-situ formation is the major channel for dense gas in the distant tail of the g… ▽ More We investigate the formation of jellyfish galaxies using radiation-hydrodynamic simulations of gas-rich dwarf galaxies with a multi-phase interstellar medium (ISM). We find that the ram-pressure-stripped (RPS) ISM is the dominant source of molecular clumps in the near wake within 10 kpc from the galactic plane, while in-situ formation is the major channel for dense gas in the distant tail of the gas-rich galaxy. Only 20% of the molecular clumps in the near wake originate from the intracluster medium (ICM); however, the fraction reaches 50% in the clumps located at $80\,{\rm kpc}$ from the galactic center since the cooling time of the RPS gas tends to be short due to the ISM--ICM mixing ($\lesssim$ 10 Myr). The tail region exhibits a star formation rate of $0.001-0.01\,{\rm M_{\odot}\,yr^{-1}}$, and most of the tail stars are born in the stripped wake within 10 kpc from the galactic plane. These stars induce bright H$α$ blobs in the tail, while H$α$ tails fainter than $6\times10^{38}\,{\rm erg\,s^{-1}\,kpc^{-2}}$ are mostly formed via collisional radiation and heating due to mixing. We also find that the stripped tails have intermediate X-ray to H$α$ surface brightness ratios (1.5$\lesssim F_{\rm X}/F_{\rm Hα}\lesssim$20), compared to the ISM ($\lesssim$1.5) or pure ICM ($\gg$20). Our results suggest that jellyfish features emerge when the ISM from gas-rich galaxies is stripped by strong ram pressure, mixes with the ICM, and enhances the cooling in the tail. △ Less

Submitted 15 February, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: 19 pages, 15 figures; Accepted for publication in ApJ

Towards convergence of turbulent dynamo amplification in cosmological simulations of galaxies

Authors: Sergio Martin-Alvarez, Julien Devriendt, Adrianne Slyz, Debora Sijacki, Mark L. A. Richardson, Harley Katz

Abstract: Our understanding of the process through which magnetic fields reached their observed strengths in present-day galaxies remains incomplete. One of the advocated solutions is a turbulent dynamo mechanism that rapidly amplifies weak magnetic field seeds to the order of ${\sim}μ$G. However, simulating the turbulent dynamo is a very challenging computational task due to the demanding span of spatial s… ▽ More Our understanding of the process through which magnetic fields reached their observed strengths in present-day galaxies remains incomplete. One of the advocated solutions is a turbulent dynamo mechanism that rapidly amplifies weak magnetic field seeds to the order of ${\sim}μ$G. However, simulating the turbulent dynamo is a very challenging computational task due to the demanding span of spatial scales and the complexity of the required numerical methods. In particular, turbulent velocity and magnetic fields are extremely sensitive to the spatial discretisation of simulated domains. To explore how refinement schemes affect galactic turbulence and amplification of magnetic fields in cosmological simulations, we compare two refinement strategies. A traditional quasi-Lagrangian adaptive mesh refinement approach focusing spatial resolution on dense regions, and a new refinement method that resolves the entire galaxy with a high resolution quasi-uniform grid. Our new refinement strategy yields much faster magnetic energy amplification than the quasi-Lagrangian method, which is also significantly greater than the adiabatic compressional estimate indicating that the extra amplification is produced through stretching of magnetic field lines. Furthermore, with our new refinement the magnetic energy growth factor scales with resolution following $\propto \Dres^{-1/2}$, in much better agreement with small-scale turbulent box simulations. Finally, we find evidence suggesting most magnetic amplification in our simulated galaxies occurs in the warm phase of their interstellar medium, which has a better developed turbulent field with our new refinement strategy. △ Less

Submitted 21 April, 2022; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: Replaced to match version accepted for publication in MNRAS. 20 pages, 17 figures and 2 appendices

Momentum deposition of supernovae with cosmic rays

Authors: Francisco Rodríguez Montero, Sergio Martin-Alvarez, Debora Sijacki, Adrianne Slyz, Julien Devriendt, Yohan Dubois

Abstract: The cataclysmic explosions of massive stars as supernovae are one of the key ingredients of galaxy formation. However, their evolution is not well understood in the presence of magnetic fields or cosmic rays (CRs). We study the expansion of individual supernova remnants (SNRs) using our suite of 3D hydrodynamical (HD), magnetohydrodynamical (MHD) and CRMHD simulations generated using RAMSES. We ex… ▽ More The cataclysmic explosions of massive stars as supernovae are one of the key ingredients of galaxy formation. However, their evolution is not well understood in the presence of magnetic fields or cosmic rays (CRs). We study the expansion of individual supernova remnants (SNRs) using our suite of 3D hydrodynamical (HD), magnetohydrodynamical (MHD) and CRMHD simulations generated using RAMSES. We explore multiple ambient densities, magnetic fields and fractions of supernova energy deposited as CRs ($χ_{\rm CR}$), accounting for cosmic ray anisotropic diffusion and streaming. All our runs have comparable evolutions until the end of the Sedov-Taylor phase. However, our CRMHD simulations experience an additional CR pressure-driven snowplough phase once the CR energy dominates inside the SNR. We present a model for the final momentum deposited by supernovae that captures this new phase: $p_{\rm SNR} = 2.87\times 10^{5} (χ_{\text{CR}} + 1)^{4.82}\left(\frac{n}{\text{cm}^{-3}}\right)^{-0.196} M_{\odot}$ km s$^{-1}$. Assuming a 10% fraction of SN energy in CRs leads to a 50% boost of the final momentum, with our model predicting even higher impacts at lower ambient densities. The anisotropic diffusion of CRs assuming an initially uniform magnetic field leads to extended gas and cosmic ray outflows escaping from the supernova poles. We also study a tangled initial configuration of the magnetic field, resulting instead in a quasi-isotropic diffusion of CRs and earlier momentum deposition. Finally, synthetic synchrotron observations of our simulations using the POLARIS code show that the local magnetic field configuration in the interstellar medium modifies the overall radio emission morphology and polarisation. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: 17 pages, 11 figures, submitted to MNRAS

Catalogues of voids as antihalos in the local Universe

Authors: Harry Desmond, Maxwell L. Hutt, Julien Devriendt, Adrianne Slyz

Abstract: A recently-proposed algorithm identifies voids in simulations as the regions associated with halos when the initial overdensity field is negated. We apply this method to the real Universe by running a suite of constrained simulations of the 2M++ volume with initial conditions inferred by the BORG algorithm, along with the corresponding inverted set. Our 101 inverted and uninverted simulations, spa… ▽ More A recently-proposed algorithm identifies voids in simulations as the regions associated with halos when the initial overdensity field is negated. We apply this method to the real Universe by running a suite of constrained simulations of the 2M++ volume with initial conditions inferred by the BORG algorithm, along with the corresponding inverted set. Our 101 inverted and uninverted simulations, spanning the BORG posterior, each identify ~150,000 "voids as antihalos" with mass exceeding $4.38\times10^{11} \: \mathrm{M_\odot}$ (100 particles) at $z=0$ in a full-sky sphere of radius 155 Mpc/h around the Milky Way. We calculate the size function, volume filling fraction, ellipticity, central density, specific angular momentum, clustering and stacked density profile of the voids, and cross-correlate them with those produced by VIDE on the same simulations. We make our antihalo and VIDE catalogues publicly available. △ Less

Submitted 24 January, 2022; v1 submitted 20 September, 2021; originally announced September 2021.

Comments: 5 pages, 5 figures; minor revision, matches MNRASL accepted version

Journal ref: MNRAS Letters 511 (1): L45-L49 (2022)

The Nature of High ${\rm [OIII]}_{\rm 88μm}$/${\rm [CII]}_{\rm 158μm}$ Galaxies in the Epoch of Reionization: Low Carbon Abundance and a Top-Heavy IMF?

Authors: Harley Katz, Joakim Rosdahl, Taysun Kimm, Thibault Garel, Jérémy Blaizot, Martin G. Haehnelt, Léo Michel-Dansac, Sergio Martin-Alvarez, Julien Devriendt, Adrianne Slyz, Romain Teyssier, Pierre Ocvirk, Nicolas Laporte, Richard Ellis

Abstract: ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μm}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenom… ▽ More ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μm}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenomenon, we employ SPHINX$^{20}$, a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization. We find that the observed $z>6$ [OIII]-SFR and [CII]-SFR relations can only be reproduced when the C/O abundance ratio is $\sim8\times$ lower than Solar and the total metal production is $\sim4\times$ higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [CII]-SFR and [CII] luminosity functions are predicted to converge to observed values at $z\sim4.5$. While we demonstrate that ionisation parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [OIII]/[CII], observed values at $z>6$ can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [CII] and [OIII] can be used to predict the values of ionisation parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed $z>6$ galaxies. Finally, we demonstrate that [OI]$_{\rm 63μm}$ can be used as a replacement for [CII] in high-redshift galaxies where [CII] is unobserved and argue that more observation time should be used to target [OI] at $z>6$. (Abridged) △ Less

Submitted 5 January, 2022; v1 submitted 2 August, 2021; originally announced August 2021.

Comments: 20 pages, 18 figures, Accepted by MNRAS

Introducing SPHINX-MHD: The Impact of Primordial Magnetic Fields on the First Galaxies, Reionization, and the Global 21cm Signal

Authors: Harley Katz, Sergio Martin-Alvarez, Joakim Rosdahl, Taysun Kimm, Jérémy Blaizot, Martin G. Haehnelt, Léo Michel-Dansac, Thibault Garel, Jose Oñorbe, Julien Devriendt, Adrianne Slyz, Omar Attia, Romain Teyssier

Abstract: We present the first results from SPHINX-MHD, a suite of cosmological radiation-magnetohydrodynamics simulations designed to study the impact of primordial magnetic fields (PMFs) on galaxy formation and the evolution of the intergalactic medium during the epoch of reionization. The simulations are among the first to employ multi-frequency, on-the-fly radiation transfer and constrained transport id… ▽ More We present the first results from SPHINX-MHD, a suite of cosmological radiation-magnetohydrodynamics simulations designed to study the impact of primordial magnetic fields (PMFs) on galaxy formation and the evolution of the intergalactic medium during the epoch of reionization. The simulations are among the first to employ multi-frequency, on-the-fly radiation transfer and constrained transport ideal MHD in a cosmological context to simultaneously model the inhomogeneous process of reionization as well as the growth of PMFs. We run a series of $(5\,\text{cMpc})^3$ cosmological volumes, varying both the strength of the seed magnetic field ($B_0$) and its spectral index ($n_B$). We find that PMFs that have $n_B > -0.562\log_{10}\left(\frac{B_0}{1{\rm n}G}\right) - 3.35$ produce electron optical depths ($τ_e$) that are inconsistent with CMB constraints due to the unrealistically early collapse of low-mass dwarf galaxies. For $n_B\geq-2.9$, our constraints are considerably tighter than the $\sim{\rm n}G$ constraints from Planck. PMFs that do not satisfy our constraints have little impact on the reionization history or the shape of the UV luminosity function. Likewise, detecting changes in the Lya forest due to PMFs will be challenging because photoionisation and photoheating efficiently smooth the density field. However, we find that the first absorption feature in the global 21cm signal is a sensitive indicator of the properties of the PMFs, even for those that satisfy our $τ_e$ constraint. Furthermore, strong PMFs can marginally increase the escape of LyC photons by up to 25\% and shrink the effective radii of galaxies by $\sim44\%$ which could increase the completeness fraction of galaxy surveys. Finally, our simulations show that surveys with a magnitude limit of ${\rm M_{UV,1500}=-13}$ can probe the sources that provide the majority of photons for reionization out to $z=12$. △ Less

Submitted 22 July, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: 29 pages, 27 figures, accepted MNRAS

Predicting the Observability of Population III Stars with ELT-HARMONI via the Helium $1640{\rmÅ}$ emission line

Authors: Kearn Grisdale, Niranjan Thatte, Julien Devriendt, Miguel Pereira-Santaella, Adrianne Slyz, Taysun Kimm, Yohan Dubois, Sukyoung K. Yi

Abstract: Population III (Pop. III) stars, as of yet, have not been detected, however as we move into the era of extremely large telescopes this is likely to change. One likely tracer for Pop. III stars is the HeII$\lambda1640$ emission line, which will be detectable by the HARMONI spectrograph on the European Extremely Large Telescope (ELT) over a broad range of redshifts ($2\leq z\leq14$). By post-process… ▽ More Population III (Pop. III) stars, as of yet, have not been detected, however as we move into the era of extremely large telescopes this is likely to change. One likely tracer for Pop. III stars is the HeII$\lambda1640$ emission line, which will be detectable by the HARMONI spectrograph on the European Extremely Large Telescope (ELT) over a broad range of redshifts ($2\leq z\leq14$). By post-processing galaxies from the cosmological, AMR-hydrodynamical simulation NewHorizon with theoretical spectral energy distributions (SED) for Pop. III stars and radiative transfer (i.e. the Yggdrasil Models and CLOUDY look-up tables respectively) we are able to compute the flux of HeII$\lambda1640$ for individual galaxies. From mock 10 hour observations of these galaxies we show that HARMONI will be able to detect Pop. III stars in galaxies up to $z\sim10$ provided Pop. III stars have a top heavy Initial Mass Function (IMF). Furthermore, we find that should Pop. III stars instead have an IMF similar to those of the Pop. I stars, the HeII$\lambda1640$ line would only be observable for galaxies with Pop. III stellar masses in excess of $10^{7}{\,M}_\odot\,{\rm yr}^{-1}$, average stellar age $<1{\rm Myr}$ at $z=4$. Finally, we are able to determine the minimal intrinsic flux required for HARMONI to detect Pop. III stars in a galaxy up to $z=10$. △ Less

Submitted 5 January, 2021; originally announced January 2021.

Comments: 23 Pages, 17 Figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

Rivers of Gas I.: Unveiling The Properties of High Redshift Filaments

Authors: Marius Ramsøy, Adrianne Slyz, Julien Devriendt, Clotilde Laigle, Yohan Dubois

Abstract: At high redshift, the cosmic web is widely expected to have a significant impact on the morphologies, dynamics and star formation rates of the galaxies embedded within it, underscoring the need for a comprehensive study of the properties of such a filamentary network. With this goal in mind, we perform an analysis of high-$z$ gas and dark matter (DM) filaments around a Milky Way-like progenitor si… ▽ More At high redshift, the cosmic web is widely expected to have a significant impact on the morphologies, dynamics and star formation rates of the galaxies embedded within it, underscoring the need for a comprehensive study of the properties of such a filamentary network. With this goal in mind, we perform an analysis of high-$z$ gas and dark matter (DM) filaments around a Milky Way-like progenitor simulated with the {\sc ramses} adaptive mesh refinement (AMR) code from cosmic scales ($\sim$1 Mpc) down to the virial radius of its DM halo host ($\sim$20 kpc at $z=4$). Radial density profiles of both gas and DM filaments are found to have the same functional form, namely a plummer-like profile modified to take into account the wall within which these filaments are embedded. Measurements of the typical filament core radius $r_0$ from the simulation are consistent with that of isothermal cylinders in hydrostatic equilibrium. Such an analytic model also predicts a redshift evolution for the core radius of filaments in fair agreement with the measured value for DM $(r_0 \propto (1+z)^{-3.18\pm 0.28})$. Gas filament cores grow as $(r_0 \propto (1+z)^{-2.72\pm 0.26})$. In both gas and DM, temperature and vorticity sharply drop at the edge of filaments, providing an excellent way to constrain the outer filament radius. When feedback is included the gas temperature and vorticity fields are strongly perturbed, hindering such a measurement in the vicinity of the galaxy. However, the core radius of the filaments as measured from the gas density field is largely unaffected by feedback, and the median central density is only reduced by about 20%. △ Less

Submitted 7 January, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: 19 pages, 11 figures

The role of AGN on the structure, kinematics and evolution of ETGs in the Horizon simulations

Authors: M. S. Rosito, S. E. Pedrosa, P. B. Tissera, N. E. Chisari, R. Dominguez-Tenreiro, Y. Dubois, S. Peirani, J. Devriendt, C. Pichon, A. Slyz

Abstract: 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 evolut… ▽ More 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 △ Less

Submitted 19 May, 2021; v1 submitted 18 December, 2020; originally announced December 2020.

Comments: 18 pages, 17 figures, Accepted in A&A

Journal ref: A&A 652, A44 (2021)

Unraveling the origin of magnetic fields in galaxies

Authors: Sergio Martin-Alvarez, Harley Katz, Debora Sijacki, Julien Devriendt, Adrianne Slyz

Abstract: Despite their ubiquity, there are many open questions regarding galactic and cosmic magnetic fields. Specifically, current observational constraints cannot rule out if magnetic fields observed in galaxies were generated in the Early Universe or are of astrophysical nature. Motivated by this we use our magnetic tracers algorithm to investigate whether the signatures of primordial magnetic fields pe… ▽ More Despite their ubiquity, there are many open questions regarding galactic and cosmic magnetic fields. Specifically, current observational constraints cannot rule out if magnetic fields observed in galaxies were generated in the Early Universe or are of astrophysical nature. Motivated by this we use our magnetic tracers algorithm to investigate whether the signatures of primordial magnetic fields persist in galaxies throughout cosmic time. We simulate a Milky Way-like galaxy in four scenarios: magnetised solely by primordial magnetic fields, magnetised exclusively by SN-injected magnetic fields, and two combined primordial + SN magnetisation cases. We find that once primordial magnetic fields with a comoving strength $B_0 >10^{-12}$ G are considered, they remain the primary source of galaxy magnetisation. Our magnetic tracers show that, even combined with galactic sources of magnetisation, when primordial magnetic fields are strong, they source the large-scale fields in the warm metal-poor phase of the simulated galaxy. In this case, the circumgalactic and intergalactic medium can be used to probe $B_0$ without risk of pollution by magnetic fields originated in the galaxy. Furthermore, whether magnetic fields are primordial or astrophysically-sourced can be inferred by studying local gas metallicity. As a result, we predict that future state-of-the-art observational facilities of magnetic fields in galaxies will have the potential to unravel astrophysical and primordial magnetic components of our Universe. △ Less

Submitted 28 July, 2021; v1 submitted 23 November, 2020; originally announced November 2020.

Comments: Accepted for publication in MNRAS. 19 pages, 13 figures and 1 appendix

Journal ref: MNRAS,504, 2517 (2021)

Dual Effects of Ram Pressure on Star Formation in Multi-phase Disk Galaxies with Strong Stellar Feedback

Authors: Jaehyun Lee, Taysun Kimm, Harley Katz, Joakim Rosdahl, Julien Devriendt, Adrianne Slyz

Abstract: We investigate the impact of ram pressure stripping due to the intracluster medium (ICM) on star-forming disk galaxies with a multi-phase interstellar medium (ISM) maintained by strong stellar feedback. We carry out radiation-hydrodynamics simulations of an isolated disk galaxy embedded in a 10^11 Msun dark matter halo with various ICM winds mimicking the cluster outskirts (moderate) and the centr… ▽ More We investigate the impact of ram pressure stripping due to the intracluster medium (ICM) on star-forming disk galaxies with a multi-phase interstellar medium (ISM) maintained by strong stellar feedback. We carry out radiation-hydrodynamics simulations of an isolated disk galaxy embedded in a 10^11 Msun dark matter halo with various ICM winds mimicking the cluster outskirts (moderate) and the central environment (strong). We find that both star formation quenching and triggering occur in ram pressure-stripped galaxies, depending on the strength of the winds. HI and H$_2$ in the outer galactic disk are significantly stripped in the presence of the moderate winds, whereas turbulent pressure provides support against ram pressure in the central region where star formation is active. Moderate ICM winds facilitate gas collapsing, increasing the total star formation rates by ~40% when the wind is oriented face-on or ~80% when it is edge-on. In contrast, strong winds rapidly blow away neutral and molecular hydrogen gas from the galaxy, suppressing the star formation by a factor of two within ~200 Myr. Dense gas clumps with N_H > 10 Msun pc^-2 are easily identified in extraplanar regions, but no significant young stellar populations are found in such clumps. In our attempts to enhance radiative cooling by adopting a colder ICM of T=10^6K, only a few additional stars are formed in the tail region, even if the amount of newly cooled gas increases by an order of magnitude. △ Less

Submitted 21 October, 2020; originally announced October 2020.

Comments: 26 pages, 17 figures, accepted for publication in ApJ

Dark-matter-deficient dwarf galaxies form via tidal stripping of dark matter in interactions with massive companions

Authors: 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

Abstract: 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 questio… ▽ More 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. △ Less

Submitted 7 January, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: Accepted for publication in MNRAS, 13 pages, 11 figures

Beyond halo mass: quenching galaxy mass assembly at the edge of filaments

Authors: Hyunmi Song, Clotilde Laigle, Ho Seong Hwang, Julien Devriendt, Yohan Dubois, Katarina Kraljic, Christophe Pichon, Adrianne Slyz, Rory Smith

Abstract: 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 residu… ▽ More 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. △ Less

Submitted 23 June, 2021; v1 submitted 31 August, 2020; originally announced September 2020.

Comments: 22 pages, 18 figures, accepted to MNRAS, comments are welcome

The origin of low-surface-brightness galaxies in the dwarf regime

Authors: 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

Abstract: 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 surfac… ▽ More 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. △ Less

Submitted 6 January, 2021; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: 16 pages, 15 figures, Accepted for publication in MNRAS

Formation of compact galaxies in the Extreme-Horizon simulation

Authors: 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, Christophe Yèche

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

Submitted 9 July, 2020; originally announced July 2020.

Comments: Submitted to A&A Link to a projection of the simulation at z = 2 in high definition: http://www.galactica-simulations.eu/EH_L50/index.html

Journal ref: A&A 643, L8 (2020)

Spatially offset black holes in the Horizon-AGN simulation and comparison to observations

Authors: Deaglan J. Bartlett, Harry Desmond, Julien Devriendt, Pedro G. Ferreira, Adrianne Slyz

Abstract: We study the displacements between the centres of galaxies and their supermassive black holes (BHs) in the cosmological hydrodynamical simulation Horizon-AGN, and in a variety of observations from the literature. The BHs in Horizon-AGN feel a sub-grid dynamical friction force, sourced by the surrounding gas, which prevents recoiling BHs being ejected from the galaxy. We find that i) the fraction o… ▽ More We study the displacements between the centres of galaxies and their supermassive black holes (BHs) in the cosmological hydrodynamical simulation Horizon-AGN, and in a variety of observations from the literature. The BHs in Horizon-AGN feel a sub-grid dynamical friction force, sourced by the surrounding gas, which prevents recoiling BHs being ejected from the galaxy. We find that i) the fraction of spatially offset BHs increases with cosmic time, ii) BHs live on prograde orbits in the plane of the galaxy with an orbital radius that decays with time but stalls near $z=0$, and iii) the magnitudes of offsets from the galaxy centres are substantially larger in the simulation than in observations. We attribute the stalling of the infall and excessive offset magnitudes to the fact that dynamical friction from stars and dark matter is not modelled in the simulation, and hence provide a way to improve the black hole dynamics of future simulations. △ Less

Submitted 3 December, 2020; v1 submitted 2 July, 2020; originally announced July 2020.

Comments: 19 pages, 15 figures, accepted for publication in MNRAS

How primordial magnetic fields shrink galaxies

Authors: Sergio Martin-Alvarez, Adrianne Slyz, Julien Devriendt, Carlos Gómez-Guijarro

Abstract: As one of the prime contributors to the interstellar medium energy budget, magnetic fields naturally play a part in shaping the evolution of galaxies. Galactic magnetic fields can originate from strong primordial magnetic fields provided these latter remain below current observational upper limits. To understand how such magnetic fields would affect the global morphological and dynamical propertie… ▽ More As one of the prime contributors to the interstellar medium energy budget, magnetic fields naturally play a part in shaping the evolution of galaxies. Galactic magnetic fields can originate from strong primordial magnetic fields provided these latter remain below current observational upper limits. To understand how such magnetic fields would affect the global morphological and dynamical properties of galaxies, we use a suite of high-resolution constrained transport magneto-hydrodynamic cosmological zoom simulations where we vary the initial magnetic field strength and configuration along with the prescription for stellar feedback. We find that strong primordial magnetic fields delay the onset of star formation and drain the rotational support of the galaxy, diminishing the radial size of the galactic disk and driving a higher amount of gas towards the centre. This is also reflected in mock UVJ observations by an increase in the light profile concentration of the galaxy. We explore the possible mechanisms behind such a reduction in angular momentum, focusing on magnetic braking. Finally, noticing that the effects of primordial magnetic fields are amplified in the presence of stellar feedback, we briefly discuss whether the changes we measure would also be expected for galactic magnetic fields of non-primordial origin. △ Less

Submitted 20 May, 2020; originally announced May 2020.

Comments: 23 pages, 15 figures, accepted for publication in MNRAS

New Methods for Identifying Lyman Continuum Leakers and Reionization-Epoch Analogues

Authors: Harley Katz, Dominika Ďurovčíková, Taysun Kimm, Joki Rosdahl, Jeremy Blaizot, Martin G. Haehnelt, Julien Devriendt, Adrianne Slyz, Richard Ellis, Nicolas Laporte

Abstract: Identifying low-redshift galaxies that emit Lyman Continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, it also remains uncertain whether the low-redshift LyC leakers are truly "analogues" of the sources that reionized the Universe. Here, we use high-… ▽ More Identifying low-redshift galaxies that emit Lyman Continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, it also remains uncertain whether the low-redshift LyC leakers are truly "analogues" of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission line properties to observed LyC leakers at $z\sim3$ and $z\sim0$. We find that the simulated galaxies with high LyC escape fractions ($f_{\rm esc}$) often exhibit high O32 and populate the same regions of the R23-O32 plane as $z\sim3$ LyC leakers. However, we show that viewing angle, metallicity, and ionisation parameter can all impact where a galaxy resides on the O32-$f_{\rm esc}$ plane. Based on emission line diagnostics and how they correlate with $f_{\rm esc}$, lower-metallicity LyC leakers at $z\sim3$ appear to be good analogues of reionization-era galaxies. In contrast, local [SII]-deficient galaxies do not overlap with the simulated high-redshift LyC leakers on the SII-BPT diagram; however, this diagnostic may still be useful for identifying leakers. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using IR and nebular emission lines. We show that our model using only [CII]$_{\rm 158μm}$ and [OIII]$_{\rm 88μm}$ can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS1149_JD1 at $z=9.1$ is the most likely galaxy to be actively contributing to the reionization of the Universe. △ Less

Submitted 5 August, 2020; v1 submitted 4 May, 2020; originally announced May 2020.

Comments: 18 pages, 11 figures, Accepted to MNRAS

Detecting the Cosmic Web: Lyα Emission from Simulated Filaments at z=3

Authors: Lydia M. Elias, Shy Genel, Amiel Sternberg, Julien Devriendt, Adrianne Slyz, Eli Visbal, Nicolas Bouché

Abstract: The standard cosmological model ($Λ$CDM) predicts the existence of the cosmic web: a distribution of matter into sheets and filaments connecting massive halos. However, observational evidence has been elusive due to the low surface brightness of the filaments. Recent deep MUSE/VLT data and upcoming observations offer a promising avenue for Ly$α$ detection, motivating the development of modern theo… ▽ More The standard cosmological model ($Λ$CDM) predicts the existence of the cosmic web: a distribution of matter into sheets and filaments connecting massive halos. However, observational evidence has been elusive due to the low surface brightness of the filaments. Recent deep MUSE/VLT data and upcoming observations offer a promising avenue for Ly$α$ detection, motivating the development of modern theoretical predictions. We use hydrodynamical cosmological simulations run with the AREPO code to investigate the potential detectability of large-scale filaments, excluding contributions from the halos embedded in them. We focus on filaments connecting massive ($M_{200c}\sim(1-3)\times10^{12} M_\odot$) halos at z=3, and compare different simulation resolutions, feedback levels, and mock-image pixel sizes. We find increasing simulation resolution does not substantially improve detectability notwithstanding the intrinsic enhancement of internal filament structure. By contrast, for a MUSE integration of 31 hours, including feedback increases the detectable area by a factor of $\simeq$5.5 on average compared with simulations without feedback, implying that even the non-bound components of the filaments have substantial sensitivity to feedback. Degrading the image resolution from the native MUSE scale of (0.2")$^2$ per pixel to (5.3")$^2$ apertures has the strongest effect, increasing the detectable area by a median factor of $\simeq$200 and is most effective when the size of the pixel roughly matches the width of the filament. Finally, we find the majority of Ly$α$ emission is due to electron impact collisional excitations, as opposed to radiative recombination. △ Less

Submitted 20 April, 2020; originally announced April 2020.

The Obelisk simulation: galaxies contribute more than AGN to HI reionization of protoclusters

Authors: 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

Abstract: 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 collisio… ▽ More 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$. △ Less

Submitted 26 July, 2022; v1 submitted 10 February, 2020; originally announced February 2020.

Comments: 24+2 pages, 18+1 figures, published in Astronomy & Astrophysics

Journal ref: A&A 653, A154 (2021)

Reionization history constraints from neural network based predictions of high-redshift quasar continua

Authors: D. Ďurovčíková, H. Katz, S. E. I. Bosman, F. B. Davies, J. Devriendt, A. Slyz

Abstract: Observations of the early Universe suggest that reionization was complete by $z\sim6$, however, the exact history of this process is still unknown. One method for measuring the evolution of the neutral fraction throughout this epoch is via observing the Ly$α$ damping wings of high-redshift quasars. In order to constrain the neutral fraction from quasar observations, one needs an accurate model of… ▽ More Observations of the early Universe suggest that reionization was complete by $z\sim6$, however, the exact history of this process is still unknown. One method for measuring the evolution of the neutral fraction throughout this epoch is via observing the Ly$α$ damping wings of high-redshift quasars. In order to constrain the neutral fraction from quasar observations, one needs an accurate model of the quasar spectrum around Ly$α$, after the spectrum has been processed by its host galaxy but before it is altered by absorption and damping in the intervening IGM. In this paper, we present a novel machine learning approach, using artificial neural networks, to reconstruct quasar continua around Ly$α$. Our QSANNdRA algorithm improves the error in this reconstruction compared to the state-of-the-art PCA-based model in the literature by 14.2% on average, and provides an improvement of 6.1% on average when compared to an extension thereof. In comparison with the extended PCA model, QSANNdRA further achieves an improvement of 22.1% and 16.8% when evaluated on low-redshift quasars most similar to the two high-redshift quasars under consideration, ULAS J1120+0641 at $z=7.0851$ and ULAS J1342+0928 at $z=7.5413$, respectively. Using our more accurate reconstructions of these two $z>7$ quasars, we estimate the neutral fraction of the IGM using a homogeneous reionization model and find $\bar{x}_\mathrm{HI} = 0.25^{+0.05}_{-0.05}$ at $z=7.0851$ and $\bar{x}_\mathrm{HI} = 0.60^{+0.11}_{-0.11}$ at $z=7.5413$. Our results are consistent with the literature and favour a rapid end to reionization. △ Less

Submitted 17 May, 2020; v1 submitted 2 December, 2019; originally announced December 2019.

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 493, Issue 3, 2020

When galaxies align: intrinsic alignments of the progenitors of elliptical galaxies in the Horizon-AGN simulation

Authors: James Bate, Nora Elisa Chisari, Sandrine Codis, Garreth Martin, Yohan Dubois, Julien Devriendt, Christophe Pichon, Adrianne Slyz

Abstract: 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… ▽ More 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. △ Less

Submitted 11 November, 2019; originally announced November 2019.

Comments: 12 pages, 9 figures, accepted in MNRAS

Simulating MOS science on the ELT: Ly$α$ forest tomography

Authors: 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

Abstract: 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$α$, forest towards densely-sampled lines-of-sight of bright background sources, such as quasars and star forming galaxies. We investigate the scientific poten… ▽ More 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$α$, 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$α$ 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] △ Less

Submitted 6 November, 2019; v1 submitted 31 October, 2019; originally announced November 2019.

Comments: 28 pages, accepted to A&A; language edited

Journal ref: A&A 632, A94 (2019)

The impact of the connectivity of the cosmic web on the physical properties of galaxies at its nodes

Authors: 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

Abstract: 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 measurem… ▽ More 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). △ Less

Submitted 17 October, 2019; originally announced October 2019.

Comments: 16 pages, 15 figures, submitted to MNRAS

The impact of AGN feedback on galaxy intrinsic alignments in the Horizon simulations

Authors: Adam Soussana, Nora Elisa Chisari, Sandrine Codis, Ricarda S. Beckmann, Yohan Dubois, Julien Devriendt, Sebastien Peirani, Clotilde Laigle, Christophe Pichon, Adrianne Slyz

Abstract: 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 sensitiv… ▽ More 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. △ Less

Submitted 7 January, 2020; v1 submitted 30 August, 2019; originally announced August 2019.

Comments: 17 pages, 13 figures, accepted in MNRAS

How to Quench a Dwarf Galaxy: The Impact of Inhomogeneous Reionization on Dwarf Galaxies and Cosmic Filaments

Authors: Harley Katz, Marius Ramsoy, Joakim Rosdahl, Taysun Kimm, Jeremy Blaizot, Martin G. Haehnelt, Leo Michel-Dansac, Thibault Garel, Clotilde Laigle, Julien Devriendt, Adrianne Slyz

Abstract: We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below t… ▽ More We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below the atomic cooling threshold. This makes them an ideal tool for examining how reionization impacts star formation and the gas content of dwarf galaxies. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80% following reionization. The gas inflow rates into haloes with $M_{vir}<10^8M_{\odot}$ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses $M_{vir}<7\times10^7M_{\odot}$. Our simulations show that inflow suppression, rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by $Δz>3$ at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with $M_{vir}<10^8M_{\odot}$ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that local group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to $z=6$... △ Less

Submitted 26 February, 2020; v1 submitted 27 May, 2019; originally announced May 2019.

Comments: 21 pages, 11 figures, MNRAS Accepted

Modelling baryonic feedback for survey cosmology

Authors: Nora Elisa Chisari, Alexander J. Mead, Shahab Joudaki, Pedro Ferreira, Aurel Schneider, Joseph Mohr, Tilman Tröster, David Alonso, Ian G. McCarthy, Sergio Martin-Alvarez, Julien Devriendt, Adrianne Slyz, Marcel P. van Daalen

Abstract: Observational cosmology in the next decade will rely on probes of the distribution of matter in the redshift range between $0<z<3$ to elucidate the nature of dark matter and dark energy. In this redshift range, galaxy formation is known to have a significant impact on observables such as two-point correlations of galaxy shapes and positions, altering their amplitude and scale dependence beyond the… ▽ More Observational cosmology in the next decade will rely on probes of the distribution of matter in the redshift range between $0<z<3$ to elucidate the nature of dark matter and dark energy. In this redshift range, galaxy formation is known to have a significant impact on observables such as two-point correlations of galaxy shapes and positions, altering their amplitude and scale dependence beyond the expected statistical uncertainty of upcoming experiments at separations under 10 Mpc. Successful extraction of information in such a regime thus requires, at the very least, unbiased models for the impact of galaxy formation on the matter distribution, and can benefit from complementary observational priors. This work reviews the current state of the art in the modelling of baryons for cosmology, from numerical methods to approximate analytical prescriptions, and makes recommendations for studies in the next decade, including a discussion of potential probe combinations that can help constrain the role of baryons in cosmological studies. We focus, in particular, on the modelling of the matter power spectrum, $P(k,z)$, as a function of scale and redshift, and of the observables derived from this quantity. This work is the result of a workshop held at the University of Oxford in November of 2018. △ Less

Submitted 25 June, 2019; v1 submitted 15 May, 2019; originally announced May 2019.

Comments: 26 pages, 4 figures, a short review, accepted in The Open Journal of Astrophysics

Group connectivity in COSMOS: a tracer of mass assembly history

Authors: 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

Abstract: 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 us… ▽ More 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. △ Less

Submitted 26 April, 2019; originally announced April 2019.

Comments: 14 pages, 10 figures, 1 table; submitted to MNRAS; Comments are welcome!

Horizon-AGN virtual observatory - 1. SED-fitting performance and forecasts for future imaging surveys

Authors: C. Laigle, I. Davidzon, O. Ilbert, J. Devriendt, D. Kashino, C. Pichon, P. Capak, S. Arnouts, S. de la Torre, Y. Dubois, G. Gozaliasl, D. Le Borgne, S. Lilly, H. J. McCracken, M. Salvato, A. Slyz

Abstract: Using the ligthcone from the cosmological hydrodynamical simulation Horizon-AGN, we produced a photometric catalogue over $0<z<4$ with apparent magnitudes in COSMOS, DES, LSST-like, and Euclid-like filters at depths comparable to these surveys. The virtual photometry accounts for the complex star formation history and metal enrichment of Horizon-AGN galaxies, and consistently includes magnitude er… ▽ More Using the ligthcone from the cosmological hydrodynamical simulation Horizon-AGN, we produced a photometric catalogue over $0<z<4$ with apparent magnitudes in COSMOS, DES, LSST-like, and Euclid-like filters at depths comparable to these surveys. The virtual photometry accounts for the complex star formation history and metal enrichment of Horizon-AGN galaxies, and consistently includes magnitude errors, dust attenuation and absorption by inter-galactic medium. The COSMOS-like photometry is fitted in the same configuration as the COSMOS2015 catalogue. We then quantify random and systematic errors of photometric redshifts, stellar masses, and star-formation rates (SFR). Photometric redshifts and redshift errors capture the same dependencies on magnitude and redshift as found in COSMOS2015, excluding the impact of source extraction. COSMOS-like stellar masses are well recovered with a dispersion typically lower than 0.1 dex. The simple star formation histories and metallicities of the templates induce a systematic underestimation of stellar masses at $z<1.5$ by at most 0.12 dex. SFR estimates exhibit a dust-induced bimodality combined with a larger scatter (typically between 0.2 and 0.6 dex). We also use our mock catalogue to predict photometric redshifts and stellar masses in future imaging surveys. We stress that adding Euclid near-infrared photometry to the LSST-like baseline improves redshift accuracy especially at the faint end and decreases the outlier fraction by a factor $\sim$2. It also considerably improves stellar masses, reducing the scatter up to a factor 3. It would therefore be mutually beneficial for LSST and Euclid to work in synergy. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Comments: 29 pages, 23 figures, 6 tables; submitted to MNRAS; Comments are welcome

The formation and evolution of low-surface-brightness galaxies

Authors: G. Martin, S. Kaviraj, C. Laigle, J. E. G. Devriendt, R. A. Jackson, S. Peirani, Y. Dubois, C. Pichon, A. Slyz

Abstract: 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 cosmol… ▽ More 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. △ Less

Submitted 12 February, 2019; originally announced February 2019.

Comments: Accepted for publication in MNRAS

On the Observed Diversity of Star Formation Efficiencies in Giant Molecular Clouds

Authors: Kearn Grisdale, Oscar Agertz, Florent Renaud, Alessandro B. Romeo, Julien Devriendt, Adrianne Slyz

Abstract: Observations find a median star formation efficiency per free-fall time in Milky Way Giant Molecular Clouds (GMCs) on the order of $ε_{\rm ff}\sim 1\%$ with dispersions of $\sim0.5\,{\rm dex}$. The origin of this scatter in $ε_{\rm ff}$ is still debated and difficult to reproduce with analytical models. We track the formation, evolution and destruction of GMCs in a hydrodynamical simulation of a M… ▽ More Observations find a median star formation efficiency per free-fall time in Milky Way Giant Molecular Clouds (GMCs) on the order of $ε_{\rm ff}\sim 1\%$ with dispersions of $\sim0.5\,{\rm dex}$. The origin of this scatter in $ε_{\rm ff}$ is still debated and difficult to reproduce with analytical models. We track the formation, evolution and destruction of GMCs in a hydrodynamical simulation of a Milky Way-like galaxy and by deriving cloud properties in an observationally motivated way, measure the distribution of star formation efficiencies which are in excellent agreement with observations. We find no significant link between $ε_{\rm ff}$ and any measured global property of GMCs (e.g. gas mass, velocity dispersion). Instead, a wide range of efficiencies exist in the entire parameter space. From the cloud evolutionary tracks, we find that each cloud follow a \emph{unique} evolutionary path which gives rise to wide diversity in all properties. We argue that it is this diversity in cloud properties, above all else, that results in the dispersion of $ε_{\rm ff}$. △ Less

Submitted 30 April, 2019; v1 submitted 1 February, 2019; originally announced February 2019.

Comments: 10 pages, 6 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Probing Cosmic Dawn with Emission Lines: Predicting Infrared and Nebular Line Emission for ALMA and JWST

Authors: Harley Katz, Thomas P. Galligan, Taysun Kimm, Joakim Rosdahl, Martin G. Haehnelt, Jeremy Blaizot, Julien Devriendt, Adrianne Slyz, Nicolas Laporte, Richard Ellis

Abstract: Infrared and nebular lines provide some of our best probes of the physics regulating the properties of the interstellar medium (ISM) at high-redshift. However, interpreting the physical conditions of high-redshift galaxies directly from emission lines remains complicated due to inhomogeneities in temperature, density, metallicity, ionisation parameter, and spectral hardness. We present a new suite… ▽ More Infrared and nebular lines provide some of our best probes of the physics regulating the properties of the interstellar medium (ISM) at high-redshift. However, interpreting the physical conditions of high-redshift galaxies directly from emission lines remains complicated due to inhomogeneities in temperature, density, metallicity, ionisation parameter, and spectral hardness. We present a new suite of cosmological, radiation-hydrodynamics simulations, each centred on a massive Lyman-break galaxy that resolves such properties in an inhomogeneous ISM. Many of the simulated systems exhibit transient but well defined gaseous disks that appear as velocity gradients in [CII]~158.6$μ$m emission. Spatial and spectral offsets between [CII]~158.6$μ$m and [OIII]~88.33$μ$m are common, but not ubiquitous, as each line probes a different phase of the ISM. These systems fall on the local [CII]-SFR relation, consistent with newer observations that question previously observed [CII]~158.6$μ$m deficits. Our galaxies are consistent with the nebular line properties of observed $z\sim2-3$ galaxies and reproduce offsets on the BPT and mass-excitation diagrams compared to local galaxies due to higher star formation rate (SFR), excitation, and specific-SFR, as well as harder spectra from young, metal-poor binaries. We predict that local calibrations between H$α$ and [OII]~3727$Å$ luminosity and galaxy SFR apply up to $z>10$, as do the local relations between certain strong line diagnostics (R23 and [OIII]~5007$Å$/H$β$) and galaxy metallicity. Our new simulations are well suited to interpret the observations of line emission from current (ALMA and HST) and upcoming facilities (JWST and ngVLA). △ Less

Submitted 13 June, 2019; v1 submitted 4 January, 2019; originally announced January 2019.

Comments: 21 pages, 16 figures, MNRAS accepted