Search

Galaxies form and evolve via a multitude of complex physics. In this work, we investigate the role of cosmic ray (CR) feedback in galaxy evolution and reionisation, by examining its impact on the escape of ionising radiation from galaxies. For this purpose, we present two Sphinx cosmological radiation-magneto-hydrodynamics simulations, allowing for the first time a study of the impact of CR feedback on thousands of resolved galaxies during the Epoch of Reionisation (EoR). The simulations differ in their feedback prescriptions: one adopts a calibrated strong supernova (SN) feedback, while the other simulation reduces the strength of SN feedback and includes CR feedback instead. We show that both comparably regulate star formation, reasonably match observations of high-redshift UV luminosity functions, and produce a similar amount of hydrogen ionising photons. In contrast to the model with strong SN feedback, the model with CRs lead to incomplete reionisation, which is in strong disagreement with observational estimates of the reionisation history. This is due to CR feedback shaping the ISM differently, filling with gas the low density cavities otherwise carved by SN explosions. As a result, this reduces the escape of ionising photons, at any halo mass, and primarily in the close vicinity of the stars. Our study indicates that CR feedback regulates galaxy growth during the EoR, but negatively affects reionisation, a tension which paves the way for further exploration and refinement of existing galaxy formation and feedback models. Such improvements are crucial in order to capture and understand the process of reionisation and the underlying evolution of galaxies through cosmic time.

We present the first constraints on primordial magnetic fields from the Lyman-$\alpha$ forest using full cosmological hydrodynamic simulations. At the scales and redshifts probed by the data, the flux power spectrum is extremely sensitive to the extra power induced by primordial magnetic fields in the linear matter power spectrum, at a scale that we parametrize with $k_{\rm peak}$. We rely on a set of more than a quarter million flux models obtained by varying thermal, reionization histories and cosmological parameters. We find a hint of extra power that is well fitted by the PMF model with $B\sim 0.2$ nG, corresponding to $k_{\rm peak}\sim 20$ Mpc$^{-1}$. However, when applying very conservative assumptions on the modelling of the noise, we obtain a 3$\sigma$ C.L. lower limit $k_{\rm peak}> 30$ Mpc$^{-1}$ which translates into the tightest bounds on the strength of primordial intergalactic magnetic fields: $B < 0.30$ nG (for fixed, nearly scale-invariant $n_{\rm B}=-2.9$)., Comment: 5 pages + Supplemental material

The visibility of the Lyman-$\alpha$ (Ly$\alpha$) emission from reionization-epoch galaxies depends sensitively on the extent of the intrinsic \lya emission redwards of 1215.67~\AA. The prominent red peak resulting from resonant radiative transfer in the interstellar medium is often modelled as a single Gaussian. We use the \textsc{Azahar} simulation suite of a massive-reionization epoch galaxy to show that a significantly larger fraction of the \lya emission extends to $400$-$800$~km~s$^{-1}$, and thus significantly further to the red than predicted by a Gaussian line profile. A cycle of frequent galaxy mergers strongly modulates the \lya luminosity, the red peak velocity and its extended red wing emerging from the galaxy, which all also strongly vary with viewing angle. The \lya emission also depends sensitively on the implemented feedback, dust and star formation physics. Our simulations including cosmic rays reproduce the observed spectral properties of reionization epoch \lya emitters (LAEs) well if we assume that the \lya emission is affected by very little dust. The visibility of LAEs can be strongly underestimated if the extended red wings of the intrinsic \lya emission are not accounted for. We discuss implications for using the visibility of LAEs to constrain the evolution of the volume-averaged neutral fraction during reionization.

We study star formation variability, or burstiness, as a method to constrain and compare different galaxy formation models at high redshift using the Azahar simulation suite. The models range from magneto-hydrodynamics with a magneto-thermo-turbulent prescription for star formation (iMHD) to more sophisticated setups incorporating radiative transfer (RTiMHD) and cosmic ray physics (RTnsCRiMHD). Analysing a sample of galaxies at redshifts $z=4-10$, we find that the RTnsCRiMHD model exhibits more regular star formation periodicity compared to iMHD and RTiMHD, as revealed by the Lomb-Scargle periodogram. While the RTiMHD model captures a notable degree of stochasticity in star formation without cosmic rays, RTnsCRiMHD galaxies display even greater scatter in the burst intensity and in the scatter around the star-forming main sequence. To evaluate the burstiness in RTnsCRiMHD against observations, we generate a mock spectrum during a mini-quenching event at $z=7.5$. This spectrum aligns well with the low-mass quiescent galaxy JADES-GS-z7-01-QU observed at $z=7.3$, though some discrepancies attributed to stellar metallicity hint at a composite spectrum. Our findings highlight the importance of including complex physical processes like cosmic rays and radiative transfer in simulations to accurately capture the bursty nature of star formation in high-redshift galaxies. Future JWST observations, particularly regarding the scatter around the star-forming main sequence, have the potential to refine and guide the next generation of galaxy formation models., Comment: 10 pages, 4 figures, 1 table, comments welcome

Understanding the role that magnetic fields play on the stage of galaxy formation requires accurate methods for inferring the properties of extragalactic magnetic fields. Radio synchrotron emission has been the most promising avenue to infer magnetic field strengths across galaxies, with the application of a central assumption: that galactic cosmic rays are in energy equipartition with the magnetic field. In this work, we leverage flexible synthetic observations of a high-resolution magnetohydrodynamic simulation of a Milky Way-like galaxy to review whether true equipartition is capable of reproducing radio observations of galaxies, and investigate its impact on the inference of magnetic field strengths when varying the properties and density distribution of the cosmic rays. We find that imposing equipartition (regardless of scale length) results in cosmic ray electron densities that are unable to generate either the amplitude or the shape of the radio intensity profiles typically observed in spiral galaxies. Instead, observationally motivated smooth distributions of cosmic ray electrons across the galaxy provide a remarkable match to observations. We further demonstrate that assuming equipartition with those mock observations can lead to significant overestimation of the magnetic field strength. This misestimation varies with cosmic ray electron densities, cosmic ray spectrum power-law index, and galactic environment, aggravated in inter-arm regions and attenuated in star-forming regions. Our results promote caution when assuming equipartition in observations, and suggest that additional theoretical and numerical work is required to leverage the upcoming generation of radio observations poised to revolutionize our understanding of astrophysical magnetic fields., Comment: Updated to match current version after addressing referee's comments. Submitted to ApJ, 26 pages, 12 figures

Understanding the impact of baryonic physics on cosmic structure formation is crucial for accurate cosmological predictions, especially as we usher in the era of large galaxy surveys with the Rubin Observatory as well as the Euclid and Roman Space Telescopes. A key process that can redistribute matter across a large range of scales is feedback from accreting supermassive black holes. How exactly these active galactic nuclei (AGN) operate from sub-parsec to Mega-parsec scales however remains largely unknown. To understand this, we investigate how different AGN feedback models in the Fable simulation suite affect the cosmic evolution of the matter power spectrum (MPS). Our analysis reveals that AGN feedback significantly suppresses clustering at scales $k \sim 10\,h\,cMpc^{-1}$, with the strongest effect at redshift $z = 0$ causing a reduction of $\sim 10\%$ with respect to the dark matter-only simulation. This is due to the efficient feedback in both radio (low Eddington ratio) and quasar (high Eddington ratio) modes in our fiducial Fable model. We find that variations of the quasar and radio mode feedback with respect to the fiducial Fable model have distinct effects on the MPS redshift evolution, with the radio mode being more effective on larger scales and later epochs. Furthermore, MPS suppression is dominated by AGN feedback effects inside haloes at $z = 0$, while for $z \gtrsim 1$ the matter distribution both inside and outside of haloes shapes the MPS suppression. Hence, future observations probing earlier cosmic times beyond $z \sim 1$ will be instrumental in constraining the nature of AGN feedback., Comment: Submitted to MNRAS. 17 pages, 10 (+2) figures. Fable MPS raw data is available at: https://github.com/MartinAlvarezSergio/Fable_MPS

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 $\Lambda$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$.

Ly$\alpha$ emission is an exceptionally informative tracer of the life cycle of evolving galaxies and the escape of ionising photons. However, theoretical studies of Ly$\alpha$ emission are often limited by insufficient numerical resolution, incomplete sets of physical models, and poor line-of-sight (LOS) statistics. To overcome such limitations, we utilize here the novel PANDORA suite of high-resolution dwarf galaxy simulations that include a comprehensive set of state-of-the-art physical models for ionizing radiation, magnetic fields, supernova feedback and cosmic rays. We post-process the simulations with the radiative transfer code \textsc{RASCAS} to generate synthetic observations and compare to observed properties of Ly$\alpha$ emitters. Our simulated Ly$\alpha$ haloes are more extended than the spatial region from which the intrinsic emission emanates and our spatially resolved maps of spectral parameters of the Ly$\alpha$ emission are very sensitive to the underlying spatial distribution and kinematics of neutral hydrogen. Ly$\alpha$ and LyC emission display strongly varying signatures along different LOS depending on how each LOS intersects low-density channels generated by stellar feedback. Comparing galaxies simulated with different physics, we find the Ly$\alpha$ signatures to exhibit systematic offsets determined by the different levels of feedback strength and the clumpiness of the neutral gas. Despite this variance, and regardless of the different physics included in each model, we find universal correlations between Ly$\alpha$ observables and LyC escape fraction, demonstrating a robust connection between Ly$\alpha$ and LyC emission. Ly$\alpha$ observations from a large sample of dwarf galaxies should thus give strong constraints on their stellar feedback-regulated LyC escape and confirm their important role for the reionization of the Universe.

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., Comment: Accepted for publication in the ApJ. Replaced to match final version. 35 pages, 17 figures, 4 tables

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., Comment: 20 pages, 11 figures, submitted MNRAS

During the epoch of reionisation the first galaxies were enshrouded in pristine neutral gas, with one of the brightest emission lines in star-forming galaxies, Lyman-$\alpha$ (Ly$\alpha$), expected to remain undetected until the Universe became ionised. Providing an explanation for the surprising detection of Ly$\alpha$ in these early galaxies is a major challenge for extra-galactic studies. Recent JWST observations have reignited the debate on whether residence in an overdensity of galaxies is a it sufficient and necessary condition for Ly$\alpha$ to escape. Here, we take unique advantage of both high-resolution and high-sensitivity images from the JWST instrument NIRCam to reveal that all galaxies in a sample of z>7 Ly$\alpha$ emitters have close companions. We exploit novel on-the-fly radiative transfer magnetohydrodynamical simulations with cosmic ray feedback to show that galaxies with frequent mergers have very bursty star formation which drives episodes of high intrinsic Ly$\alpha$ emission and facilitates the escape of Ly$\alpha$ photons along channels cleared of neutral gas. We conclude that the rapid build up of stellar mass through mergers presents a compelling solution to the long-standing puzzle of the detection of Ly$\alpha$ emission deep into the epoch of reionisation., Comment: Published in Nature Astronomy. 45 pages, 9 figures, 2 tables

We present the analysis of the magnetic field ($B$-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby ($<20$ Mpc) galaxies with resolved (5 arcsec-18 arcsec; $90$ pc--$1$ kpc) imaging polarimetric observations using HAWC+/SOFIA from $53$ to $214$ \um. We compute the magnetic pitch angle ($\Psi_{B}$) profiles as a function of the galactrocentric radius. We introduce a new magnetic alignment parameter ($\zeta$) to estimate the disordered-to-ordered $B$-field ratio in spiral $B$-fields. We find FIR and radio wavelengths to not generally trace the same $B$-field morphology in galaxies. The $\Psi_{B}$ profiles tend to be more ordered with galactocentric radius in radio ($\zeta_{\rm{6cm}} = 0.93\pm0.03$) than in FIR ($\zeta_{\rm{154\mu m}} = 0.84\pm0.14$). For spiral galaxies, FIR $B$-fields are $2-75$\% more turbulent than the radio $B$-fields. For starburst galaxies, we find that FIR polarization is a better tracer of the $B$-fields along the galactic outflows than radio polarization. Our results suggest that the $B$-fields associated with dense, dusty, turbulent star-forming regions, those traced at FIR, are less ordered than warmer, less-dense regions, those traced at radio, of the interstellar medium. The FIR $B$-fields seem to be more sensitive to the activity of the star-forming regions and the morphology of the molecular clouds within a vertical height of few hundred pc in the disk of spiral galaxies than the radio $B$-fields., Comment: 28 pages, 14 figures. Accepted for publication in The Astrophysical Journal

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 thermo-turbulent star formation and mechanical supernova feedback, we gradually increase the complexity of physics incorporated leading to full-physics models combining magnetism, on-the-fly radiative transfer and the corresponding stellar photoheating, and SN-accelerated cosmic rays. We investigate 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 HI masses as well as kinematics. Furthermore, resolved kinematic maps of our full-physics models predict kinematically distinct clumps and kinematic misalignments of stars, HI and HII 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 multi-phase interstellar medium properties that JWST will probe., Comment: Replaced to match the MNRAS published version

Mergers are thought to be a fundamental channel for galaxy growth, perturbing the gas dynamics and the magnetic fields (B-fields) in the interstellar medium (ISM). However, the mechanisms that amplify and dissipate B-fields during a merger remain unclear. We characterize the morphology of the ordered B-fields in the multi-phase ISM of the closest merger of two spiral galaxies, the Antennae galaxies. We compare the inferred B-fields using $154~\mu$m thermal dust and $11$ cm radio synchrotron emission polarimetric observations. We find that the $154~\mu$m B-fields are more ordered across the Antennae galaxies than the $11$ cm B-fields. The turbulent-to-ordered $154~\mu$m B-field increases at the galaxy cores and star-forming regions. The relic spiral arm has an ordered spiral $154~\mu$m B-field, while the $11$ cm B-field is radial. The $154~\mu$m B-field may be dominated by turbulent dynamos with high $^{12}$CO(1-0) velocity dispersion driven by star-forming regions, while the $11$ cm B-field is cospatial with high HI velocity dispersion driven by galaxy interaction. This result shows the dissociation between the warm gas mainly disturbed by the merger, and the dense gas still following the dynamics of the relic spiral arm. We find a $\sim8.9$ kpc scale ordered B-field connecting the two galaxies. The base of the tidal tail is cospatial with the HI and $^{12}$CO(1-0) emission and has compressed and/or sheared $154~\mu$m and $11$ cm B-fields driven by the merger. We suggest that amplify B-fields, with respect to the rest of the system and other spiral galaxies, may be supporting the gas flow between both galaxies and the tidal tail., Comment: 11 pages, 5 figures, Accepted for publication in ApJ Letters

We measure escape fractions, $f_{\rm esc}$, of ionizing radiation from galaxies in the SPHINX suite of cosmological radiation-hydrodynamical simulations of reionization, resolving halos with $M_{\rm vir} \gtrapprox 7.5 \times 10^7 \ M_{\odot}$ with a minimum cell width of $\approx 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 $10^{11} \ M_{\odot}$. The simulated galaxies agree well with observational constraints of the UV luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over timescales of a few Myrs, 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 inter-galactic medium. Statistically, $f_{\rm esc}$ peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies ($M_{*} \approx 10^7 \ M_{\odot}$, $M_{1500} \approx -17$, $Z\lesssim 5 \times 10^{-3} \ Z_{\odot}$), 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} \gtrapprox -17$, provide about $55$ percent of the photons contributing to reionization. The global averaged $f_{\rm esc}$ 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., Comment: 20 figures, 21 pages plus appendices. Accepted for publication in MNRAS

We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $\mu$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation and sample selection of the program. Here, we focus on the general trends in the emissive polarization fraction. Far-infrared polarimetric observations trace the thermal polarized emission of magnetically aligned dust grains across the galaxy disks with polarization fractions of $P=0-15$% in the cold, $T_{\rm d} = [19,48]$ K, and dense, $\log_{10}(N_{\rm HI+H_{2}}) = [19.96,22.91]$, interstellar medium. The spiral galaxies show a median $\langle P_{154\mu m} \rangle = 3.3\pm0.9 $% across the disks. We report the first polarized spectrum of starburst galaxies showing a minimum within $89-154$ $\mu$m. The falling $53-154$ $\mu$m polarized spectrum may be due to a decrease in the dust grain alignment efficiency produced by variations in dust temperatures along the line-of-sight in the galactic outflow. We find that the starburst galaxies and the star-forming regions within normal galaxies have the lowest polarization fractions. We find that 50% (7 out of 14) of the galaxies require a broken power-law in the $P-N_{HI+H_{2}}$ and $P-T_{d}$ relations with three different trends. Group 1 has a relative increase of anisotropic random B-fields produced by compression or shear of B-fields in the galactic outflows, starburst rings, and inner-bar of galaxies; and Groups 2 and 3 have a relative increase of isotropic random B-fields driven by star-forming regions in the spiral arms, and/or an increase of dust grain alignment efficiency caused by shock-driven regions or evolutionary stages of a galaxy., Comment: 43 pages, 28 figures. Accepted for publication to ApJ

We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipeline steps to obtain homogeneously reduced high-level data products of polarimetric maps of galaxies for use in scientific analysis. Our approach has a general design and can be applied to sources smaller than the field-of-view of the HAWC+ array in any given band. We estimate that the OTFMAP polarimetric mode offers a reduction of observing overheads by a factor 2.34, and an improvement in sensitivity by a factor 1.80 when compared to previously obtained polarimetric observations using the chopping and nodding mode. The OTFMAP is a significant optimization of the polarimetric mode of HAWC+ as it ultimately reduces the cost of operations of SOFIA/HAWC+ by increasing the science collected per hour of observation up to an overall factor of 2.49. The OTFMAP polarimetric mode is the standard observing strategy of SALSA. The results and quantitative analysis of this first data release are presented in Papers IV and V of the series., Comment: 26 pages, 16 figures, 5 tables. Submitted to ApJ (Comments welcome)

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 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 $10^{10}$, $10^{11}$ and $10^{12}\, \rm M_{\odot}$ with a maximum resolution of $9 \,\rm pc$. We find that CRs reduce star formation rate in our two dwarf galaxies by a factor 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 $\kappa = 10^{27}$ to $\rm 3\times 10^{29}\, cm^2\, s^{-1}$. With lower $\kappa$, CRs remain confined for longer on small scales and are consequently efficient in suppressing star formation, whereas a higher diffusion coefficient reduces the effect on star formation and increases the generation of cold outflows. Finally, we compare CR feedback to a calibrated 'strong' supernova feedback model known to sufficiently regulate star formation in high-redshift cosmological simulations. We find that CR feedback is not sufficiently strong to replace this strong supernova 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.

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}\mu$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., Comment: Replaced to match version accepted for publication in MNRAS. 20 pages, 17 figures and 2 appendices

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 ($\chi_{\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} (\chi_{\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., Comment: 17 pages, 11 figures, submitted to MNRAS

ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88\mu m}$/[CII]$_{\rm 158\mu m}$ ratios and [CII]$_{\rm 158\mu 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\mu 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), Comment: 20 pages, 18 figures, Accepted by MNRAS

We investigate the effect of the Biermann battery during the Epoch of Reionization (EoR) using cosmological Adaptive Mesh Refinement simulations within the framework of the SPHINX project. We develop a novel numerical technique to solve for the Biermann battery term in the Constrained Transport method, preserving both the zero divergence of the magnetic field and the absence of Biermann battery for isothermal flows. The structure-preserving nature of our numerical method turns out to be very important to minimise numerical errors during validation tests of the propagation of a Str\"omgren sphere and of a Sedov blast wave. We then use this new method to model the evolution of a 2.5 and 5 co-moving Mpc cosmological box with a state-of-the-art galaxy formation model within the RAMSES code. Contrary to previous findings, we show that three different Biermann battery channels emerge: the first one is associated with linear perturbations before the EoR, the second one is the classical Biermann battery associated with reionization fronts during the EoR, and the third one is associated with strong, supernova-driven outflows. While the two former channels generate spontaneously volume-filling magnetic fields with a strength on the order or below $10^{-20}$ G, the latter, owing to the higher plasma temperature and a marginally-resolved turbulent dynamo, reaches a field strength as high as $10^{-18}$ G in the intergalactic medium around massive haloes., Comment: 14 pages, 18 figures, accepted in MNRAS

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 ($\tau_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 $\tau_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$., Comment: 29 pages, 27 figures, accepted MNRAS

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

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

Optically-compact star-forming galaxies (SFGs) have been proposed as immediate progenitors of quiescent galaxies, although their origin and nature are debated. Were they formed in slow secular processes or in rapid merger-driven starbursts? Addressing this question would provide fundamental insight into how quenching occurs. We explore the location of the general population of galaxies with respect to fundamental star-forming and structural relations, identify compact SFGs based on their stellar core densities, and study three diagnostics of the burstiness of star formation: 1) Star formation efficiency, 2) interstellar medium (ISM), and 3) radio emission. The overall distribution of galaxies in the fundamental relations points towards a smooth transition towards quiescence while galaxies grow their stellar cores, although some galaxies suddenly increase their specific star-formation rate when they become compact. From their star formation efficiencies compact and extended SFGs appear similar. In relation to the ISM diagnostic, by studying the CO excitation, the density of the neutral gas, and the strength of the ultraviolet radiation field, compact SFGs resemble galaxies located in the upper envelope of the SFGs main sequence, although yet based on a small sample size. Regarding the radio emission diagnostic we find that galaxies become increasingly compact as the starburst ages, implying that at least some compact SFGs are old starbursts. We suggest that compact SFGs could be starburts winding down and eventually crossing the main sequence towards quiescence., Comment: Accepted by ApJ. 18 pages, 8 figures

Observational cosmology in the next decade will rely on probes of the distribution of matter in the redshift range between $0

Despite their ubiquity, the origin of cosmic magnetic fields remains unknown. Various mechanisms have been proposed for their existence including primordial fields generated by inflation, or amplification and injection by compact astrophysical objects. Separating the potential impact of each magnetogenesis scenario on the magnitude and orientation of the magnetic field and their impact on gas dynamics may give insight into the physics that magnetised our Universe. In this work, we demonstrate that because the induction equation and solenoidal constraint are linear with $B$, the contribution from different sources of magnetic field can be separated in cosmological magnetohydrodynamics simulations and their evolution and influence on the gas dynamics can be tracked. Exploiting this property, we develop a magnetic field tracer algorithm for cosmological simulations that can track the origin and evolution of different components of the magnetic field. We present a suite of cosmological magnetohydrodynamical RAMSES simulations that employ this algorithm where the primordial field strength is varied to determine the contributions of the primordial and supernovae-injected magnetic fields to the total magnetic energy as a function of time and spatial location. We find that, for our specific model, the supernova-injected fields rarely penetrate far from haloes, despite often dominating the total magnetic energy in the simulations. The magnetic energy density from the supernova-injected field scales with density with a power-law slope steeper than 4/3 and often dominates the total magnetic energy inside of haloes. However, the star formation rates in our simulations are not affected by the presence of magnetic fields, for the ranges of primordial field strengths examined. These simulations represent a first demonstration of the magnetic field tracer algorithm (abridged)., Comment: 12 pages, 9 figures, MNRAS Accepted

Arguably the main challenge of galactic magnetism studies is to explain how the interstellar medium of galaxies reaches energetic equipartition despite the extremely weak cosmic primordial magnetic fields that are originally predicted to thread the inter-galactic medium. Previous numerical studies of isolated galaxies suggest that a fast dynamo amplification might suffice to bridge the gap spanning many orders of magnitude in strength between the weak early Universe magnetic fields and the ones observed in high redshift galaxies. To better understand their evolution in the cosmological context of hierarchical galaxy growth, we probe the amplification process undergone by the cosmic magnetic field within a spiral galaxy to unprecedented accuracy by means of a suite of constrained transport magnetohydrodynamical adaptive mesh refinement cosmological zoom simulations with different stellar feedback prescriptions. A galactic turbulent dynamo is found to be naturally excited in this cosmological environment, being responsible for most of the amplification of the magnetic energy. Indeed, we find that the magnetic energy spectra of simulated galaxies display telltale inverse cascades. Overall, the amplification process can be divided in three main phases, which are related to different physical mechanisms driving galaxy evolution: an initial collapse phase, an accretion-driven phase, and a feedback-driven phase. While different feedback models affect the magnetic field amplification differently, all tested models prove to be subdominant at early epochs, before the feedback-driven phase is reached. Thus the three-phase evolution paradigm is found to be quite robust vis-a-vis feedback prescriptions., Comment: 24 pages, 20 figures, accepted for publication in MNRAS