933 research outputs found
Cosmic ray feedback in the FIRE simulations: constraining cosmic ray propagation with GeV gamma ray emission
We present the implementation and the first results of cosmic ray (CR)
feedback in the Feedback In Realistic Environments (FIRE) simulations. We
investigate CR feedback in non-cosmological simulations of dwarf, sub-
starburst, and galaxies with different propagation models, including
advection, isotropic and anisotropic diffusion, and streaming along field lines
with different transport coefficients. We simulate CR diffusion and streaming
simultaneously in galaxies with high resolution, using a two moment method. We
forward-model and compare to observations of -ray emission from nearby
and starburst galaxies. We reproduce the -ray observations of dwarf and
galaxies with constant isotropic diffusion coefficient . Advection-only and streaming-only
models produce order-of-magnitude too large -ray luminosities in dwarf
and galaxies. We show that in models that match the -ray
observations, most CRs escape low-gas-density galaxies (e.g.\ dwarfs) before
significant collisional losses, while starburst galaxies are CR proton
calorimeters. While adiabatic losses can be significant, they occur only after
CRs escape galaxies, so they are only of secondary importance for -ray
emissivities. Models where CRs are ``trapped'' in the star-forming disk have
lower star formation efficiency, but these models are ruled out by -ray
observations. For models with constant that match the -ray
observations, CRs form extended halos with scale heights of several kpc to
several tens of kpc.Comment: 31 pages, 26 figures, accepted for publication in MNRA
The origin of ultra diffuse galaxies: stellar feedback and quenching
We test if the cosmological zoom-in simulations of isolated galaxies from the
FIRE project reproduce the properties of ultra diffuse galaxies. We show that
stellar feedback-generated outflows that dynamically heat galactic stars,
together with a passively aging stellar population after imposed quenching
(from e.g. infall into a galaxy cluster), naturally reproduce the observed
population of red UDGs, without the need for high spin halos or dynamical
influence from their host cluster. We reproduce the range of surface
brightness, radius and absolute magnitude of the observed z=0 red UDGs by
quenching simulated galaxies at a range of different times. They represent a
mostly uniform population of dark matter-dominated galaxies with M_star ~1e8
Msun, low metallicity and a broad range of ages. The most massive simulated
UDGs require earliest quenching and are therefore the oldest. Our simulations
provide a good match to the central enclosed masses and the velocity
dispersions of the observed UDGs (20-50 km/s). The enclosed masses of the
simulated UDGs remain largely fixed across a broad range of quenching times
because the central regions of their dark matter halos complete their growth
early. A typical UDG forms in a dwarf halo mass range of Mh~4e10-1e11 Msun. The
most massive red UDG in our sample requires quenching at z~3 when its halo
reached Mh ~ 1e11 Msun. If it, instead, continues growing in the field, by z=0
its halo mass reaches > 5e11 Msun, comparable to the halo of an L* galaxy. If
our simulated dwarfs are not quenched, they evolve into bluer low-surface
brightness galaxies with mass-to-light ratios similar to observed field dwarfs.
While our simulation sample covers a limited range of formation histories and
halo masses, we predict that UDG is a common, and perhaps even dominant, galaxy
type around Ms~1e8 Msun, both in the field and in clusters.Comment: 20 pages, 13 figures; match the MNRAS accepted versio
Breathing FIRE: How Stellar Feedback Drives Radial Migration, Rapid Size Fluctuations, and Population Gradients in Low-Mass Galaxies
We examine the effects of stellar feedback and bursty star formation on
low-mass galaxies ()
using the FIRE (Feedback in Realistic Environments) simulations. While previous
studies emphasized the impact of feedback on dark matter profiles, we
investigate the impact on the stellar component: kinematics, radial migration,
size evolution, and population gradients. Feedback-driven outflows/inflows
drive significant radial stellar migration over both short and long timescales
via two processes: (1) outflowing/infalling gas can remain star-forming,
producing young stars that migrate within their first , and (2) gas outflows/inflows drive strong fluctuations in the
global potential, transferring energy to all stars. These processes produce
several dramatic effects. First, galaxies' effective radii can fluctuate by
factors of over , and these rapid size fluctuations
can account for much of the observed scatter in radius at fixed
Second, the cumulative effects of many outflow/infall episodes steadily heat
stellar orbits, causing old stars to migrate outward most strongly. This
age-dependent radial migration mixes---and even inverts---intrinsic age and
metallicity gradients. Thus, the galactic-archaeology approach of calculating
radial star-formation histories from stellar populations at can be
severely biased. These effects are strongest at , the same regime where feedback most
efficiently cores galaxies. Thus, detailed measurements of stellar kinematics
in low-mass galaxies can strongly constrain feedback models and test baryonic
solutions to small-scale problems in CDM.Comment: Accepted to ApJ (820, 131) with minor revisions from v1. Figure 4 now
includes dark matter. Main results in Figures 7 and 1
On the deuterium abundance and the importance of stellar mass loss in the interstellar and intergalactic medium
We quantify the gas-phase abundance of deuterium and fractional contribution
of stellar mass loss to the gas in cosmological zoom-in simulations from the
Feedback In Realistic Environments project. At low metallicity, our simulations
confirm that the deuterium abundance is very close to the primordial value. The
chemical evolution of the deuterium abundance that we derive here agrees
quantitatively with analytical chemical evolution models. We furthermore find
that the relation between the deuterium and oxygen abundance exhibits very
little scatter. We compare our simulations to existing high-redshift
observations in order to determine a primordial deuterium fraction of 2.549 +/-
0.033 x 10^-5 and stress that future observations at higher metallicity can
also be used to constrain this value. At fixed metallicity, the deuterium
fraction decreases slightly with decreasing redshift, due to the increased
importance of mass loss from intermediate-mass stars. We find that the
evolution of the average deuterium fraction in a galaxy correlates with its
star formation history. Our simulations are consistent with observations of the
Milky Way's interstellar medium: the deuterium fraction at the solar circle is
85-92 per cent of the primordial deuterium fraction. We use our simulations to
make predictions for future observations. In particular, the deuterium
abundance is lower at smaller galactocentric radii and in higher mass galaxies,
showing that stellar mass loss is more important for fuelling star formation in
these regimes (and can even dominate). Gas accreting onto galaxies has a
deuterium fraction above that of the galaxies' interstellar medium, but below
the primordial fraction, because it is a mix of gas accreting from the
intergalactic medium and gas previously ejected or stripped from galaxies.Comment: Accepted for publication in MNRAS. Revised version: expanded
discussion and added Figure 2 (residual dependence on iron abundance
Properties of the circumgalactic medium in cosmic ray-dominated galaxy haloes
We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass haloes hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, and collisional and streaming losses, with constant parallel diffusivity κ∼3×10²⁹ cm² s⁻¹ chosen to match γ-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass haloes at z ≲ 1–2. The gas in these ‘CR-dominated’ haloes differs significantly from runs without CRs: the gas is primarily cool (a few ∼10⁴), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the ‘low’ and ‘mid’ ions in this diffuse cool gas is dominated by photoionization, with O VI columns ≳10^(14.5) cm⁻² at distances ≳150kpc. CR and thermal gas pressure are locally anticorrelated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same haloes are primarily warm/hot (T≳10⁵) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase
But What About... Cosmic Rays, Magnetic Fields, Conduction, & Viscosity in Galaxy Formation
We present a suite of high-resolution cosmological simulations, using the
FIRE-2 feedback physics together with explicit treatment of magnetic fields,
anisotropic conduction and viscosity, and cosmic rays (CRs) injected by
supernovae (including anisotropic diffusion, streaming, adiabatic, hadronic and
Coulomb losses). We survey systems from ultra-faint dwarf (, ) through Milky Way
masses, systematically vary CR parameters (e.g. the diffusion coefficient
and streaming velocity), and study an ensemble of galaxy properties
(masses, star formation histories, mass profiles, phase structure,
morphologies). We confirm previous conclusions that magnetic fields,
conduction, and viscosity on resolved (pc) scales have small
effects on bulk galaxy properties. CRs have relatively weak effects on all
galaxy properties studied in dwarfs (, ), or at high redshifts (), for
any physically-reasonable parameters. However at higher masses () and , CRs can suppress star
formation by factors , given relatively high effective diffusion
coefficients . At lower
, CRs take too long to escape dense star-forming gas and lose energy to
hadronic collisions, producing negligible effects on galaxies and violating
empirical constraints from -ray emission. But around , CRs escape the galaxy and build up a
CR-pressure-dominated halo which supports dense, cool ( K) gas
that would otherwise rain onto the galaxy. CR heating (from collisional and
streaming losses) is never dominant.Comment: 35 pages, 23 figures. Updated to match published (MNRAS) versio
The failure of stellar feedback, magnetic fields, conduction, and morphological quenching in maintaining red galaxies
The quenching "maintenance'" and related "cooling flow" problems are
important in galaxies from Milky Way mass through clusters. We investigate this
in halos with masses , using
non-cosmological high-resolution hydrodynamic simulations with the FIRE-2
(Feedback In Realistic Environments) stellar feedback model. We specifically
focus on physics present without AGN, and show that various proposed "non-AGN"
solution mechanisms in the literature, including Type Ia supernovae, shocked
AGB winds, other forms of stellar feedback (e.g. cosmic rays), magnetic fields,
Spitzer-Braginskii conduction, or "morphological quenching" do not halt or
substantially reduce cooling flows nor maintain "quenched" galaxies in this
mass range. We show that stellar feedback (including cosmic rays from SNe)
alters the balance of cold/warm gas and the rate at which the cooled gas within
the galaxy turns into stars, but not the net baryonic inflow. If anything,
outflowing metals and dense gas promote additional cooling. Conduction is
important only in the most massive halos, as expected, but even at reduces inflow only by a factor (owing to
saturation effects and anisotropic suppression). Changing the morphology of the
galaxies only slightly alters their Toomre- parameter, and has no effect on
cooling (as expected), so has essentially no effect on cooling flows or
maintaining quenching. This all supports the idea that additional physics,
e.g., AGN feedback, must be important in massive galaxies.Comment: 16 pages, 12 figure
CHIANTI - an Atomic Database for Emission Lines. Paper VI: Proton Rates and Other Improvements
The CHIANTI atomic database contains atomic energy levels, wavelengths,
radiative transition probabilities and electron excitation data for a large
number of ions of astrophysical interest. Version 4 has been released, and
proton excitation data is now included, principally for ground configuration
levels that are close in energy. The fitting procedure for excitation data,
both electrons and protons, has been extended to allow 9 point spline fits in
addition to the previous 5 point spline fits. This allows higher quality fits
to data from close-coupling calculations where resonances can lead to
significant structure in the Maxwellian-averaged collision strengths. The
effects of photoexcitation and stimulated emission by a blackbody radiation
field in a spherical geometry on the level balance equations of the CHIANTI
ions can now be studied following modifications to the CHIANTI software. With
the addition of H I, He I and N I, the first neutral species have been added to
CHIANTI. Many updates to existing ion data-sets are described, while several
new ions have been added to the database, including Ar IV, Fe VI and Ni XXI.
The two-photon continuum is now included in the spectral synthesis routines,
and a new code for calculating the relativistic free-free continuum has been
added. The treatment of the free-bound continuum has also been updated.Comment: CHIANTI is available at http://wwwsolar.nrl.navy.mil/chianti.htm
- …