-
Galaxies with grains: unraveling dust evolution and extinction curves with hydrodynamical simulations
Authors:
Yohan Dubois,
Francisco Rodríguez Montero,
Corentin Guerra,
Maxime Trebitsch,
San Han,
Ricarda Beckmann,
Sukyoung K. Yi,
Joseph Lewis,
J. K. Jang
Abstract:
We introduce a model for dust evolution in the RAMSES code for simulations of galaxies with a resolved multiphase interstellar medium. Dust is modelled as a fluid transported with the gas component, and is decomposed into two sizes, 5 nm and 0.1 $μ\rm m$, and two chemical compositions for carbonaceous and silicate grains. Using a suite of isolated disc simulations with different masses and metalli…
▽ More
We introduce a model for dust evolution in the RAMSES code for simulations of galaxies with a resolved multiphase interstellar medium. Dust is modelled as a fluid transported with the gas component, and is decomposed into two sizes, 5 nm and 0.1 $μ\rm m$, and two chemical compositions for carbonaceous and silicate grains. Using a suite of isolated disc simulations with different masses and metallicities, the simulations can explore the role of these processes in shaping the key properties of dust in galaxies. The simulated Milky Way analogue reproduces the dust-to-metal mass ratio (DTM), depletion factors, size distribution and extinction curves of the Milky Way. Galaxies with lower metallicities reproduce the observed decrease in the DTM with metallicity at around a few 0.1 $\rm Z_\odot$. This break in the DTM corresponds to a galactic gas metallicity threshold that marks the transition from an ejecta-dominated to an accretion-dominated grain growth, and that is different for silicate and carbonaceous grains, with $\simeq$ 0.1 $\rm Z_\odot$ and $\simeq$ 0.5 $\rm Z_\odot$ respectively. This leads to more Magellanic Cloud-like extinction curves, i.e. with steeper slopes in the ultraviolet and a weaker bump feature at 217.5 nm, in galaxies with lower masses and lower metallicities. Steeper slopes in these galaxies are caused by the combination of the higher efficiency of gas accretion by silicate relative to carbonaceous grains and by the low rates of coagulation that preserves the amount of small silicate grains. Weak bumps are due to the overall inefficient accretion growth of carbonaceous dust at low metallicity, whose growth is mostly supported by the release of large grains in SN ejecta. We also show that the formation of CO molecules is a key component to limit the ability of carbonaceous dust to grow, in particular in low-metallicity gas-rich galaxies.
△ Less
Submitted 4 June, 2024; v1 submitted 28 February, 2024;
originally announced February 2024.
-
Extragalactic Magnetism with SOFIA (SALSA Legacy Program). VII. A Tomographic View of Far-infrared and Radio Polarimetric Observations through MHD Simulations of Galaxies
Authors:
Sergio Martin-Alvarez,
Enrique Lopez-Rodriguez,
Tara Dacunha,
Susan E. Clark,
Alejandro S. Borlaff,
Rainer Beck,
Francisco Rodríguez Montero,
S. Lyla Jung,
Julien Devriendt,
Adrianne Slyz,
Julia Roman-Duval,
Evangelia Ntormousi,
Mehrnoosh Tahani,
Kandaswamy Subramanian,
Daniel A. Dale,
Pamela M. Marcum,
Konstantinos Tassis,
Ignacio del Moral-Castro,
Le Ngoc Tram,
Matt J. Jarvis
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 24 March, 2024; v1 submitted 10 November, 2023;
originally announced November 2023.
-
The impact of cosmic rays on the interstellar medium and galactic outflows of Milky Way analogues
Authors:
Francisco Rodríguez Montero,
Sergio Martin-Alvarez,
Adrianne Slyz,
Julien Devriendt,
Yohan Dubois,
Debora Sijacki
Abstract:
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 ou…
▽ More
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.
△ Less
Submitted 25 July, 2023;
originally announced July 2023.
-
Momentum deposition of supernovae with cosmic rays
Authors:
Francisco Rodríguez Montero,
Sergio Martin-Alvarez,
Debora Sijacki,
Adrianne Slyz,
Julien Devriendt,
Yohan Dubois
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 ex…
▽ More
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 ($χ_{\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} (χ_{\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.
△ Less
Submitted 19 October, 2021;
originally announced October 2021.
-
Rapidly quenched galaxies in the Simba cosmological simulation and observations
Authors:
Yirui Zheng,
Romeel Dave,
Vivienne Wild,
Francisco Rodríguez Montero
Abstract:
A wide range of mechanisms have been put forward to explain the quenching of star formation in galaxies with cosmic time, however, the true balance of responsible mechanisms remains unknown. The identification and study of galaxies that have shut down their star formation on different timescales might elucidate which mechanisms dominate at different epochs and masses. Here we study the population…
▽ More
A wide range of mechanisms have been put forward to explain the quenching of star formation in galaxies with cosmic time, however, the true balance of responsible mechanisms remains unknown. The identification and study of galaxies that have shut down their star formation on different timescales might elucidate which mechanisms dominate at different epochs and masses. Here we study the population of rapidly quenched galaxies (RQGs) in the SIMBA cosmological hydrodynamic simulation at $0.5<z<2$, comparing directly to observational post-starburst galaxies in the UKIDSS Ultra Deep Survey via their colour distributions and mass functions. We find that the fraction of quiescent galaxies that are rapidly quenched in SIMBA is 59% (or 48% in terms of stellar mass), which is higher than observed. A similar "downsizing" of RQGs is observed in both SIMBA and the UDS, with RQGs at higher redshift having a higher average mass. However, SIMBA produces too many RQGs at $1<z_q<1.5$ and too few low mass RQGs at $0.5<z_q<1$. The precise colour distribution of SIMBA galaxies compared to the observations also indicates various inconsistencies in star formation and chemical enrichment histories, including an absence of short, intense starbursts. Our results will help inform the next generation of galaxy evolution models, particularly with respect to the quenching mechanisms employed.
△ Less
Submitted 30 March, 2022; v1 submitted 5 October, 2021;
originally announced October 2021.
-
Mergers, Starbursts, and Quenching in the Simba Simulation
Authors:
Francisco Rodríguez Montero,
Romeel Davé,
Vivienne Wild,
Daniel Anglés-Alcázar,
Desika Narayanan
Abstract:
We use the Simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\leq$ 4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass $M_*$ well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate sSFR$>t_{H}^{-1}$ to sSFR$<0.2t_{H}^{-1}$, where…
▽ More
We use the Simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\leq$ 4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass $M_*$ well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate sSFR$>t_{H}^{-1}$ to sSFR$<0.2t_{H}^{-1}$, where $t_{H}$ is the Hubble time. At $z\approx 0-3$, mergers show $\sim\times 2-3$ higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $\leq 1\%$ of the cosmic SF budget. At low masses, the increase in SFR in mergers is mostly attributed to an increase in the $H_2$ content, but for $M_*\geq 10^{10.5} M_{\odot}$ mergers also show an elevated star formation efficiency suggesting denser gas within merging galaxies. The merger rate for star-forming galaxies shows a rapid increase with redshift $\propto (1+z)^{3.5}$, but the quenching rate evolves much more slowly, $\propto (1+z)^{0.9}$; there are insufficient mergers to explain the quenching rate at $z\leq 1.5$. Simba first quenches galaxies at $z\geq 3$, with a number density in good agreement with observations. The quenching timescales $τ_q$ are strongly bimodal, with `slow' quenchings ($τ_q \sim 0.1t_{H}$) dominating overall, but `fast' quenchings ($τ_q\sim 0.01 t_H$) dominating in $M_*\sim 10^{10}-10^{10.5}M_{\odot}$ galaxies, likely induced by Simba's jet-mode black hole feedback. The delay time distribution between mergers and quenching events suggests no physical connection to either fast or slow quenching. Hence, Simba predicts that major mergers induce starbursts, but are unrelated to quenching in either fast or slow mode.
△ Less
Submitted 12 September, 2019; v1 submitted 29 July, 2019;
originally announced July 2019.