Julien Devriendt

Abstract:

Cosmic ray (CR) feedback has been proposed as a powerful mechanism for driving warm gas outflows in galaxies. We use cosmological magnetohydrodynamic simulations to investigate the impact of CR feedback on neutral hydrogen (H i) in a 1011 M⊙ dark matter halo at 2 < z < 4. To this end, we postprocess the simulations with ionizing radiative transfer and perform Monte Carlo Lyman-α (Lyα) transfer calculations. CR feedback reduces H i column densities around young stars, thereby allowing more Lyα photons to escape and consequently offering a better match to the Lyα luminosities of observed Lyα emitters. Although galaxies with CR-driven outflows have more extended H i in the circumgalactic medium, two Lyα line properties sensitive to optical depth and gas kinematic—the location of the red peak relative to the Lyα line center in velocity space (vred) and relative strength of the blue-to-red peaks (B/R)—cannot distinguish between the CR-driven and non-CR simulations. This is because Lyα photons propagate preferentially along low H i density channels created by the ionizing radiation, thereby limiting the scattering with volume-filling H i. In contrast, the observed low flux ratios between the valley and peak and the surface brightness profiles are better reproduced in the model with CR-driven outflows because the Lyα photons interact more before escaping, rather than being destroyed by dust as is the case in the non-CR simulation. We discuss the potential cause of the paucity of sightlines in simulations that exhibit prominent red peaks and large vred, which may require the presence of more volume-filling H i.

Abstract:

Supermassive black holes are prevalent at the centers of massive galaxies, and their masses scale with galaxy properties, increasing evidence suggesting that these trends continue to low stellar masses. Seeds are needed for supermassive black holes, especially at the highest redshifts explored by the James Webb Space Telescope. We study the hierarchical merging of galaxies via cosmological merger trees and argue that the seeds of supermassive black holes formed in nuclear star clusters via stellar black hole mergers at early epochs. Observable tracers include intermediate-mass black holes, nuclear star clusters, and early gas accretion in host dwarf galaxies, along with a potentially detectable stochastic gravitational-wave background, ejection of intermediate and supermassive black holes, and consequences of a significant population of early tidal disruption events and extreme mass ratio inspirals.

Abstract:

Recent observations by James Webb Space Telescope reveal an unexpectedly abundant population of rapidly growing supermassive black holes (SMBHs) in the early Universe, underscoring the need for improved models for their origin and growth. Employing new full radiative transfer hydrodynamical simulations of galaxy formation, we investigate the local and intergalactic feedback of SMBH progenitors for the Population III.1 (Pop III.1) scenario, i.e. efficient formation of supermassive stars from pristine, undisturbed dark matter minihaloes. Our cosmological simulations capture the R-type expansion phase of these Pop III.1 stars, with their H-ionizing photon luminosities of generating H ii regions that extend deep into the intergalactic medium, reaching comoving radii of . We vary both the Pop III.1 ionization flux and cosmological formation environments, finding the former regulates their final , whereas the latter is more important in setting their formation redshift. We use the results from our radiation-hydrodynamics simulations to estimate the cosmic number density of SMBHs, , expected from Pop III.1 progenitors. We find , consistent with the results inferred from recent observations of the local and high-redshift universe. Overall, this establishes Pop III.1 progenitors as viable candidates for the formation of the first SMBH, and emphasizes the importance of exploring heavy mass seed scenarios.

Abstract:

Dwarf galaxies serve as powerful laboratories for investigating the underlying physics of galaxy evolution including the impact of baryonic feedback processes and external environmental influences. We compare the visual and structural properties of dwarf galaxies in ultra-deep HSC–SSP imaging of the COSMOS field with those measured from realistic HSC-like synthetic observations of dwarfs generated by the Illustris TNG50 and NewHorizon simulations. Using Sérsic profile fitting and non-parametric morphological metrics (Gini, , asymmetry, and concentration), we evaluate the diversity of structural properties in observed and simulated galaxies. Our analysis shows that NewHorizon and TNG50 galaxies lie at opposite extremes of observed structural trends: NewHorizon produces diffuse, extended galaxies with shallow Sérsic indices, while TNG50 yields compact, concentrated systems with steep indices. Both simulations reproduce observed structural trends more closely at higher stellar masses () but fail to capture the full diversity of COSMOS dwarfs at lower masses. Non-parametric metrics further show that NewHorizon galaxies exhibit more uneven, clumpy light distributions while TNG50 galaxies have smoother but excessively concentrated profiles. These structural differences reflect underlying differences in their physical prescriptions and are likely driven by differing approaches to ISM gas physics, supernova feedback and star formation in addition to differences in numerical resolution. Our findings highlight the unique power of low-mass galaxies to constrain differences in simulation physics, especially star formation and feedback prescriptions. Upcoming surveys from facilities like the Vera C. Rubin Observatory and Euclid will enable more rigorous comparisons with simulations, offering deeper insights into the physical processes shaping galaxy evolution.

Abstract:

We use the cosmological hydrodynamical simulation NewHorizon to study the effects of the baryonic component on the inner mass profile of dark matter haloes of isolated galaxies (). Dark matter deficits (‘cores’) develop only in galaxies in a narrow range of stellar mass, . The lower stellar mass limit arises because a minimum amount of star formation is required to drive the baryonic outflows that redistribute dark matter and create a core. The upper limit roughly coincides with the total amount of dark matter initially contained within the innermost 2 kpc (), which roughly coincides with the stellar half-mass radius of these dwarfs. This enclosed mass is quite insensitive to the total virial mass of the system. The same upper limit applies to other simulations, like NIHAO and EAGLE-CHT10, despite their rather different galaxy formation efficiencies. This suggests that it is the galaxy total stellar mass that determines when a core is formed, and not the galaxy-to-dark halo mass ratio, as argued in earlier work. This is consistent with a back-of-the-envelope estimate for a SN-induced rate of orbital diffusion. Although NewHorizon dwarfs reproduce the observed diversity of rotation curves better than other simulations, there are significant differences in the gravitational importance of baryons in the inner regions of dwarfs compared to observations. These differences prevent us from concluding that cosmological simulations are currently fully able to account for the observed diversity of rotation curve shapes.