Joakim Rosdahl

In this paper, we explore the significant, non-linear impact that stellar winds have on H ii regions. We perform a parameter study using three-dimensional radiative magnetohydrodynamic simulations of wind and ultraviolet radiation feedback from a 35 M⊙ star formed self-consistently in a turbulent, self-gravitating cloud, similar to the Orion Nebula (M42) and its main ionizing source θ1 Ori C. Stellar winds suppress early radiative feedback by trapping ionizing radiation in the shell around the wind bubble. Rapid breakouts of warm photoionized gas (‘champagne flows’) still occur if the star forms close to the edge of the cloud. The impact of wind bubbles can be enhanced if we detect and remove numerical overcooling caused by shocks crossing grid cells. However, the majority of the energy in the wind bubble is still lost to turbulent mixing between the wind bubble and the gas around it. These results begin to converge if the spatial resolution at the wind bubble interface is increased by refining the grid on pressure gradients. Wind bubbles form a thin chimney close to the star, which then expands outwards as an extended plume once the wind bubble breaks out of the dense core the star formed in, allowing them to expand faster than a spherical wind bubble. We also find wind bubbles mixing completely with the photoionized gas when the H ii region breaks out of the cloud as a champagne flow, a process we term ‘hot champagne’.

Enshrouded in several well-known controversies, dwarf galaxies have been extensively studied to learn about the underlying cosmology, notwithstanding that physical processes regulating their properties are poorly understood. To shed light on these processes, we introduce the Pandora suite of 17 high-resolution (3.5 parsec half-cell side) dwarf galaxy formation cosmological simulations. Commencing with magneto-thermo-turbulent star formation and mechanical supernova (SN) feedback, we gradually increase the complexity of physics incorporated, ultimately leading to our full-physics models combining magnetism, on-the-fly radiative transfer and the corresponding stellar photoheating, and SN-accelerated cosmic rays. We investigate multiple combinations of these processes, comparing them with observations to constrain what are the main mechanisms determining dwarf galaxy properties. We find hydrodynamical ‘SN feedback-only’ simulations struggle to produce realistic dwarf galaxies, leading either to overquenched or too centrally concentrated, dispersion-dominated systems when compared to observed field dwarfs. Accounting for radiation with cosmic rays results in extended and rotationally supported systems. Spatially ‘distributed’ feedback leads to realistic stellar and H i masses, galaxy sizes, and integrated kinematics. Furthermore, resolved kinematic maps of our full-physics models predict kinematically distinct clumps and kinematic misalignments of stars, H i, and H ii after star formation events. Episodic star formation combined with its associated feedback induces more core-like dark matter central profiles, which our ‘SN feedback-only’ models struggle to achieve. Our results demonstrate the complexity of physical processes required to capture realistic dwarf galaxy properties, making tangible predictions for integral field unit surveys, radio synchrotron emission, and for galaxy and multiphase interstellar medium properties that JWST will probe.

The detection of Lyα nebulae around z ≳ 6 quasars reveals extended gas reservoirs around the first rapidly growing supermassive black holes. Observations of z > 6 quasars can be explained by cosmological models provided that the black holes by which they are powered evolve in rare, massive dark matter haloes. Whether these models also explain the observed extended Lyα emission remains an open question. We post-process a suite of cosmological, radiation-hydrodynamic simulations, targeting a quasar host halo at z > 6 with the Lyα radiative transfer code Rascas. A combination of recombination radiation from photo-ionized hydrogen and emission from collisionally excited gas powers Lyα nebulae with a surface brightness profile in a close agreement with observations. We also find that, even on its own, resonant scattering of the Lyα line associated to the quasar’s broad-line region can also generate Lyα emission on $\sim 100 \, \rm kpc$ scales, resulting in comparable agreement with observed surface brightness profiles. Even if powered by a broad quasar Lyα line, Lyα nebulae have narrow line-widths $\lesssim 1000 \, \rm km \, s^{-1}$, consistent with observational constraints. Even if there is no quasar, we find that halo gas cooling produces a faint, extended Lyα glow. However, to explain the brightest Lyα nebulae, our simulations unambiguously require quasar-powered outflows to clear out the galactic nucleus and allow the Lyα flux to escape and still remain resonant with halo gas. The close match between observations and simulations with quasar outflows suggest that AGN feedback already operates before $z \, = \, 6$ and confirms that high-z quasars reside in massive haloes tracing overdensities.

We measure escape fractions, fesc, of ionizing radiation from galaxies in the sphinx suite of cosmological radiation-hydrodynamical simulations of reionization, resolving haloes with $M_{\rm vir}\gtrsim 7.5 \times 10^7 \ {\rm {M}_{\odot }}$ with a minimum cell width of ≈10 pc. Our new and largest 20 co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with halo masses up to a few times 1011 M⊙. The simulated galaxies agree well with observational constraints of the ultraviolet (UV) luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over time-scales of a few Myr, due to its regulation by supernova and radiation feedback, and at any given time a tiny fraction of star-forming galaxies emits a large fraction of the ionizing radiation escaping into the intergalactic medium. Statistically, fesc peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies (M* ≈ 107 M⊙, M1500 ≈ −17, Z ≲ 5 × 10−3 Z⊙), dropping strongly for lower and higher masses, brighter and dimmer galaxies, and more metal-rich galaxies. The escape fraction correlates positively with both the short-term and long-term specific star formation rate. According to sphinx, galaxies too dim to be yet observed, with ${M_{1500}}\gtrsim -17$, provide about 55 per cent of the photons contributing to reionization. The global averaged fesc naturally decreases with decreasing redshift, as predicted by UV background models and low-redshift observations. This evolution is driven by decreasing specific star formation rates over cosmic time.

Cosmic rays (CRs) are thought to play an important role in galaxy evolution. We study their effect when coupled to other important sources of feedback, namely supernovae (SNe) and stellar radiation, by including CR anisotropic diffusion and radiative losses but neglecting CR streaming. Using the ramses-rt code, we perform the first radiation-magnetohydrodynamics simulations of isolated disc galaxies with and without CRs. We study galaxies embedded in dark matter haloes of 1010, 1011, and $10^{12}\, \rm M_{\odot }$ with a maximum resolution of $9 \, \rm pc$. We find that CRs reduce the star formation (SF) rate in our two dwarf galaxies by a factor of 2, with decreasing efficiency with increasing galaxy mass. They increase significantly the outflow mass loading factor in all our galaxies and make the outflows colder. We study the impact of the CR diffusion coefficient, exploring values from κ = 1027 to $\rm 3\times 10^{29}\, cm^2\, s^{-1}$. With a lower κ, CRs remain confined for longer on small scales and are consequently efficient in suppressing SF, whereas a higher diffusion coefficient reduces the effect on SF and increases the generation of cold outflows. Finally, we compare CR feedback to a calibrated ’strong’ SN feedback model known to sufficiently regulate SF in high-redshift cosmological simulations. We find that CR feedback is not sufficiently strong to replace this strong SN feedback. As they tend to smooth out the ISM and fill it with denser gas, CRs also lower the escape fraction of Lyman continuum photons from galaxies.