Search

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.

We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($, Comment: 43 pages, 28 figures. Accepted for publication to ApJ

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.

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

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 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)

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

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

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

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

Observational cosmology in the next decade will rely on probes of the distribution of matter in the redshift range between $0, Comment: 26 pages, 4 figures, a short review, accepted in The Open Journal of Astrophysics

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)., 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., 24 pages, 20 figures, accepted for publication in MNRAS

Cosmological shock waves are tracers of the thermal history of the structures in the Universe. They play a crucial role in redistributing the energy within the cosmic structures and are also amongst the main ingredients of galaxy and galaxy cluster formation. Understanding this important function requires a proper description of the interplay between shocks and the different environments where they can be found. In this paper, an Adaptive Mesh Refinement (AMR) Eulerian cosmological simulation is analysed by means of a shock-finding algorithm that allows to generate shock wave maps. Based on the population of dark matter halos and on the distribution of density contrast in the simulation, we classify the shocks in five different environments. These range from galaxy clusters to voids. The shock distribution function and the shocks power spectrum are studied for these environments dynamics. We find that shock waves on different environments undergo different formation and evolution processes, showing as well different characteristics. We identify three different phases of formation, evolution and dissipation of these shock waves, and an intricate migration between distinct environments and scales. Shock waves initially form at external, low density regions and are merged and amplified through the collapse of structures. Shock waves and cosmic structures follow a parallel evolution. Later on, shocks start to detach from them and dissipate. We also find that most of the power that shock waves dissipate is found at scales of $k \sim0.5~\mbox{Mpc}^{-1}$ , with a secondary peak at $k \sim8~\mbox{Mpc}^{-1}$ . The evolution of the shocks power spectrum confirms that shock waves evolution is coupled and conditioned by their environment.