Julien Devriendt

Abstract:

We use the NEWHORIZON simulation to study the redshift evolution of bar properties and fractions within galaxies in the stellar masses range M⋆ = 107.25–1011.4M⊙ over the redshift range of z = 0.25–1.3. We select disc galaxies using stellar kinematics as a proxy for galaxy morphology. We employ two different automated bar detection methods, coupled with visual inspection, resulting in observable bar fractions of fbar = 0.070+0.018−0.012 at z ∼ 1.3, decreasing to fbar = 0.011+0.014−0.003 at z ∼ 0.25. Only one galaxy is visually confirmed as strongly barred in our sample. This bar is hosted by the most massive disc and only survives from z = 1.3 down to z = 0.7. Such a low bar fraction, in particular amongst Milky Way-like progenitors, highlights a missing bars problem, shared by literally all cosmological simulations with spatial resolution <100 pc to date. The analysis of linear growth rates, rotation curves, and derived summary statistics of the stellar, gas and dark matter components suggest that galaxies with stellar masses below 109.5−1010M⊙ in NEWHORIZON appear to be too dominated by dark matter relative to stellar content to form a bar, while more massive galaxies typically have formed large bulges that prevent bar persistence at low redshift. This investigation confirms that the evolution of the bar fraction puts stringent constraints on the assembly history of baryons and dark matter on to galaxies.

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 explore multiple ambient densities, magnetic fields, and fractions of supernova energy deposited as CRs (χCR), accounting for CR 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: pSNR=2.87×105(χCR+1)4.82(ncm−3)−0.196M⊙ km s−1. Assuming a 10 per cent fraction of SN energy in CRs leads to a 50 per cent 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 CR 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 polarization.

Abstract:

A recently proposed algorithm identifies voids in simulations as the regions associated with haloes 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 antihaloes’ with mass exceeding 4.38 × 1011 M⊙ (100 particles) at z = 0 in a full-sky sphere of radius 155 Mpc h−1 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.

Abstract:

ALMA observations of z > 6 galaxies hav e rev ealed abnormally high [O III ] 88 μm /[C II ] 158 μm ratios and [C II ] 158 μm deficits compared to local galaxies. The origin of this behaviour is unknown. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties, including ionization parameter, gas density, or photodissociation region (PDR) covering fraction. 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 which has been post-processed with CLOUDY to predict emission lines. We find that the observed z > 6 [O III ] 88 μm -SFR and [C II ] 158 μm -SFR relations can only be reproduced when the C/O abundance ratio is ∼8 ×lower than Solar and the total metal production is ∼4 ×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 [C II ] 158 μm -SFR and [C II ] 158 μm luminosity functions are predicted to converge to observed values at z ∼4.5. While we demonstrate that ionization parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [O III ] 88 μm /[C II ] 158 μm , observed values at z > 6 can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [C II ] 158 μm and [O III ] 88 μm can be used to predict the values of ionization 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 [O I ] 63 μm can be used as a replacement for [C II ] 158 μm in high-redshift galaxies where [C II ] 158 μm is unobserved and argue that more observation time should be used to target [O I ] 63 μm at z > 6. Future simulations will be needed to self-consistently address the numerous uncertainties surrounding a varying IMF at high redshift and the associated metal returns.