Simulations of radiative turbulent mixing layers Ji (季索清), Suoqing; Oh, S Peng; Masterson, Phillip
Monthly notices of the Royal Astronomical Society,
07/2019, Letnik:
487, Številka:
1
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
Radiative turbulent mixing layers should be ubiquitous in multi-phase gas with shear flow. They are a potentially attractive explanation for the high ions such as O vi seen in high-velocity ...clouds and the circumgalactic medium (CGM) of galaxies. We perform 3D magnetohydrohynamics (MHD) simulations with non-equilibrium (NEI) and photoionization modelling, with an eye towards testing simple analytic models. Even purely hydrodynamic collisional ionization equilibrium (CIE) calculations have column densities much lower than observations. Characteristic inflow and turbulent velocities are much less than the shear velocity, and the layer width $h \propto t_{\mathrm{cool}}^{1/2}$ rather than h ∝ tcool. Column densities are not independent of density or metallicity as analytic scalings predict, and show surprisingly weak dependence on shear velocity and density contrast. Radiative cooling, rather than Kelvin–Helmholtz instability, appears paramount in determining the saturated state. Low pressure due to fast cooling both seeds turbulence and sets the entrainment rate of hot gas, whose enthalpy flux, along with turbulent dissipation, energizes the layer. Regardless of initial geometry, magnetic fields are amplified and stabilize the mixing layer via magnetic tension, producing almost laminar flow and depressing column densities. NEI effects can boost column densities by factors of a few. Suppression of cooling by NEI or photoionization can, in principle, also increase O vi column densities, but, in practice, is unimportant for CGM conditions. To explain observations, sightlines must pierce hundreds or thousands of mixing layers, which may be plausible if the CGM exists as a ‘fog’ of tiny cloudlets.
This study is motivated by recent observations of ubiquitous interstellar density filaments and guided by modern theories of compressible magnetohydrodynamic (MHD) turbulence. The interstellar ...turbulence shapes the observed density structures. As the fundamental dynamics of compressible MHD turbulence, perpendicular turbulent mixing of density fluctuations entails elongated density structures aligned with the local magnetic field, accounting for low-density parallel filaments seen in diffuse atomic and molecular gas. The elongation of low-density parallel filaments depends on the turbulence anisotropy. When taking into account the partial ionization, we find that the minimum width of parallel filaments in the cold neutral medium and molecular clouds is determined by the neutral-ion decoupling scale perpendicular to magnetic field. In highly supersonic MHD turbulence in molecular clouds, both low-density parallel filaments due to anisotropic turbulent mixing and high-density filaments due to shock compression exist.
ABSTRACT
The microphysics of ∼ GeV cosmic ray (CR) transport on galactic scales remain deeply uncertain, with almost all studies adopting simple prescriptions (e.g. constant diffusivity). We explore ...different physically motivated, anisotropic, dynamical CR transport scalings in high-resolution cosmological Feedback In Realistic Environment (FIRE) simulations of dwarf and ∼L* galaxies where scattering rates vary with local plasma properties motivated by extrinsic turbulence (ET) or self-confinement (SC) scenarios, with varying assumptions about e.g. turbulent power spectra on un-resolved scales, Alfvén-wave damping, etc. We self-consistently predict observables including γ-rays (Lγ), grammage, residence times, and CR energy densities to constrain the models. We demonstrate many non-linear dynamical effects (not captured in simpler models) tend to enhance confinement. For example, in multiphase media, even allowing arbitrary fast transport in neutral gas does not substantially reduce CR residence times (or Lγ), as transport is rate-limited by the ionized WIM and ‘inner CGM’ gaseous halo (104–106 K gas within $\lesssim 10\!-\!30\,$ kpc), and Lγ can be dominated by trapping in small ‘patches’. Most physical ET models contribute negligible scattering of ∼1–10 GeV CRs, but it is crucial to account for anisotropy and damping (especially of fast modes) or else scattering rates would violate observations. We show that the most widely assumed scalings for SC models produce excessive confinement by factors ≳100 in the warm ionized medium (WIM) and inner CGM, where turbulent and Landau damping dominate. This suggests either a breakdown of quasi-linear theory used to derive the CR transport parameters in SC, or that other novel damping mechanisms dominate in intermediate-density ionized gas.
ABSTRACT
We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with ...explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$, $M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$ pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$, $M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses ($M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$) and z ≲ 1–2, CRs can suppress SF and stellar masses by factors ∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients $\kappa \gtrsim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around $\kappa \sim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 106 K) gas that would otherwise rain on to the galaxy. CR ‘heating’ (from collisional and streaming losses) is never dominant.
Turbulent dynamo field amplification has often been invoked to explain the strong field strengths in thin rims in supernova shocks (approx.100 micrograms) and in radio relics in galaxy clusters ...(approx. micrograms). We present high-resolution magnetohydrodynamic simulations of the interaction between pre-shock turbulence, clumping and shocks, to quantify the conditions under which turbulent dynamo amplification can be significant. We demonstrate numerically converged field amplification which scales with Alfven Mach number, B/B0 varies as MA, up to MA approx.150.This implies that the post-shock field strength is relatively independent of the seed field. Amplification is dominated by compression at low MA, and stretching (turbulent amplification) at high MA. For high MA, the B-field grows exponentially and saturates at equipartition with turbulence, while the vorticity jumps sharply at the shock and subsequently decays; the resulting field is orientated predominately along the shock normal (an effect only apparent in 3D and not 2D). This agrees with the radial field bias seen in supernova remnants. By contrast, for low MA, field amplification is mostly compressional, relatively modest, and results in a predominantly perpendicular field. The latter is consistent with the polarization seen in radio relics. Our results are relatively robust to the assumed level of gas clumping. Our results imply that the turbulent dynamo may be important for supernovae, but is only consistent with the field strength, and not geometry, for cluster radio relics. For the latter, this implies strong pre-existing B-fields in the ambient cluster outskirts.
ABSTRACT
We present the first simulations evolving resolved spectra of cosmic rays (CRs) from MeV–TeV energies (including electrons, positrons, (anti)protons, and heavier nuclei), in live ...kinetic-magnetohydrodynamics galaxy simulations with star formation and feedback. We utilize new numerical methods including terms often neglected in historical models, comparing Milky Way analogues with phenomenological scattering coefficients ν to Solar-neighbourhood Local interstellar medium (LISM) observations (spectra, B/C, e+/e−, $\mathrm{\bar{p}}/\mathrm{p}$, 10Be/9Be, ionization, and γ-rays). We show it is possible to reproduce observations with simple single-power-law injection and scattering coefficients (scaling with rigidity R), similar to previous (non-dynamical) calculations. We also find: (1) The circumgalactic medium in realistic galaxies necessarily imposes an $\sim 10\,$ kpc CR scattering halo, influencing the required ν(R). (2) Increasing the normalization of ν(R) re-normalizes CR secondary spectra but also changes primary spectral slopes, owing to source distribution and loss effects. (3) Diffusive/turbulent reacceleration is unimportant and generally sub-dominant to gyroresonant/streaming losses, which are sub-dominant to adiabatic/convective terms dominated by $\sim 0.1-1\,$ kpc turbulent/fountain motions. (4) CR spectra vary considerably across galaxies; certain features can arise from local structure rather than transport physics. (5) Systematic variation in CR ionization rates between LISM and molecular clouds (or Galactic position) arises naturally without invoking alternative sources. (6) Abundances of CNO nuclei require most CR acceleration occurs around when reverse shocks form in SNe, not in OB wind bubbles or later Sedov–Taylor stages of SNe remnants.
Type Ia supernovae (SNe Ia) play a crucial role as standardizable cosmological candles, though the nature of their progenitors is a subject of active investigation. Recent observational and ...theoretical work has pointed to merging white dwarf binaries, referred to as the double-degenerate channel, as the possible progenitor systems for some SNe Ia. Additionally, recent theoretical work suggests that mergers which fail to detonate may produce magnetized, rapidly rotating white dwarfs. In this paper, we present the first multidimensional simulations of the post-merger evolution of white dwarf binaries to include the effect of the magnetic field. In these systems, the two white dwarfs complete a final merger on a dynamical timescale, and are tidally disrupted, producing a rapidly rotating white dwarf merger surrounded by a hot corona and a thick, differentially rotating disk. The disk is strongly susceptible to the magnetorotational instability (MRI), and we demonstrate that this leads to the rapid growth of an initially dynamically weak magnetic field in the disk, the spin-down of the white dwarf merger, and to the subsequent central ignition of the white dwarf merger. Additionally, these magnetized models exhibit new features not present in prior hydrodynamic studies of white dwarf mergers, including the development of MRI turbulence in the hot disk, magnetized outflows carrying a significant fraction of the disk mass, and the magnetization of the white dwarf merger to field strengths ~2 x 10 super(8) G. We discuss the impact of our findings on the origins, circumstellar media, and observed properties of SNe Ia and magnetized white dwarfs.
ABSTRACT
We study the effects of cosmic rays (CRs) on outflows from star-forming galaxies in the circum and intergalactic medium (CGM/IGM), in high-resolution, fully cosmological FIRE-2 simulations ...(accounting for mechanical and radiative stellar feedback, magnetic fields, anisotropic conduction/viscosity/CR diffusion and streaming, and CR losses). We showed previously that massive ($M_{\rm halo}\gtrsim 10^{11}\, \mathrm{M}_{\odot }$), low-redshift (z ≲ 1–2) haloes can have CR pressure dominate over thermal CGM pressure and balance gravity, giving rise to a cooler CGM with an equilibrium density profile. This dramatically alters outflows. Absent CRs, high gas thermal pressure in massive haloes ‘traps’ galactic outflows near the disc, so they recycle. With CRs injected in supernovae as modelled here, the low-pressure halo allows ‘escape’ and CR pressure gradients continuously accelerate this material well into the IGM in ‘fast’ outflows, while lower-density gas at large radii is accelerated in situ into ‘slow’ outflows that extend to >Mpc scales. CGM/IGM outflow morphologies are radically altered: they become mostly volume-filling (with inflow in a thin mid-plane layer) and coherently biconical from the disc to >Mpc. The CR-driven outflows are primarily cool ($T\sim \! 10^{5}\,$ K) and low velocity. All of these effects weaken and eventually vanish at lower halo masses ($\lesssim 10^{11}\, \mathrm{M}_{\odot }$) or higher redshifts (z ≳ 1–2), reflecting the ratio of CR to thermal + gravitational pressure in the outer halo. We present a simple analytical model that explains all of the above phenomena. We caution that these predictions may depend on uncertain CR transport physics.
Abstract
Project AMIGA (Absorption Maps In the Gas of Andromeda) is a survey of the circumgalactic medium (CGM) of Andromeda (M31,
≃ 300 kpc) along 43 QSO sightlines at impact parameters 25 ≤
R
... ≤ 569 kpc (25 at
R
≲
). We use ultraviolet absorption measurements of Si
ii
, Si
iii
, Si
iv
, C
ii
, and C
iv
from the Hubble Space Telescope/Cosmic Origins Spectrograph and O
vi
from the Far Ultraviolet Spectroscopic Explorer to provide an unparalleled look at how the physical conditions and metals are distributed in the CGM of M31. We find that Si
iii
and O
vi
have a covering factor near unity for
R
≲ 1.2
and ≲1.9
, respectively, demonstrating that M31 has a very extended ∼10
4
–10
5.5
K ionized CGM. The metal and baryon masses of the 10
4
–10
5.5
K CGM gas within
are ≳10
8
and ≳4 × 10
10
(
Z
/0.3
Z
⊙
)
−1
M
⊙
, respectively. There is not much azimuthal variation in the column densities or kinematics, but there is with
R
. The CGM gas at
R
≲ 0.5
is more dynamic and has more complicated, multiphase structures than at larger radii, perhaps a result of more direct impact of galactic feedback in the inner regions of the CGM. Several absorbers are projected spatially and kinematically close to M31 dwarf satellites, but we show that those are unlikely to give rise to the observed absorption. Cosmological zoom simulations of ∼
L
* galaxies have O
vi
extending well beyond
as observed for M31 but do not reproduce well the radial column density profiles of the lower ions. However, some similar trends are also observed, such as the lower ions showing a larger dispersion in column density and stronger dependence on
R
than higher ions. Based on our findings, it is likely that the Milky Way has a ∼10
4
–10
5.5
K CGM as extended as for M31 and their CGM (especially the warm–hot gas probed by O
vi
) are overlapping.
Abstract We present a statistical study of the properties of diffuse H i in 10 nearby galaxies, comparing the H i detected by the single-dish telescope FAST (FEASTS program) and the interferometer ...Very Large Array (THINGS program), respectively. The THINGS observation missed H i with a median of 23% due to the short-spacing problem of interferometry and limited sensitivity. We extract the diffuse H i by subtracting the dense H i , which is obtained from the THINGS data with a uniform flux-density threshold, from the total H i detected by FAST. Among the sample, the median diffuse-H i fraction is 34%, and more diffuse H i is found in galaxies exhibiting more prominent tidal-interaction signatures. The diffuse H i we detected seems to be distributed in disk-like layers within a typical thickness of 1 kpc, different from the more halo-like diffuse H i detected around NGC 4631 in a previous study. Most of the diffuse H i is cospatial with the dense H i and has a typical column density of 10 17.7 –10 20.1 cm −2 . The diffuse and dense H i exhibit a similar rotational motion, but the former lags by a median of 25% in at least the inner disks, and its velocity dispersions are typically twice as high. Based on a simplified estimation of circumgalactic medium properties and assuming pressure equilibrium, the volume density of diffuse H i appears to be constant within each individual galaxy, implying its role as a cooling interface. Comparing with existing models, these results are consistent with a possible link between tidal interactions, the formation of diffuse H i , and gas accretion.