Inhomogeneous Interstellar Medium Research Articles

CLASSY VII Lyα Profiles: The Structure and Kinematics of Neutral Gas and Implications for LyC Escape in Reionization-era Analogs* * Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the Data Archive at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555.

Lyα line profiles are a powerful probe of interstellar medium (ISM) structure, outflow speed, and Lyman-continuum escape fraction. In this paper, we present the Lyα line profiles of the Cosmic Origins Spectrograph (COS) Legacy Archive Spectroscopic SurveY, a sample rich in spectroscopic analogs of reionization-era galaxies. A large fraction of the spectra show a complex profile, consisting of a double-peaked Lyα emission profile in the bottom of a damped, Lyα absorption trough. Such profiles reveal an inhomogeneous ISM. We successfully fit the damped Lyα absorption and the Lyα emission profiles separately, but with complementary covering factors, a surprising result because this approach requires no Lyα exchange between high-N H i and low-N H i paths. The combined distribution of column densities is qualitatively similar to the bimodal distributions observed in numerical simulations. We find an inverse relation between Lyα peak separation and the [O iii]/[O ii] flux ratio, confirming that the covering fraction of Lyman-continuum-thin sightlines increases as the Lyα peak separation decreases. We combine measurements of Lyα peak separation and Lyα red peak asymmetry in a diagnostic diagram, which identifies six Lyman-continuum leakers in the COS Legacy Archive Spectrocopy SurveY (CLASSY) sample. We find a strong correlation between the Lyα trough velocity and the outflow velocity measured from interstellar absorption lines. We argue that greater vignetting of the blueshifted Lyα peak, relative to the redshifted peak, is the source of the well-known discrepancy between shell-model parameters and directly measured outflow properties. The CLASSY sample illustrates how scattering of Lyα photons outside the spectroscopic aperture reshapes Lyα profiles because the distances to these compact starbursts span a large range.

Constraining gas metal mixing strength in simulations using observations of the Milky Way’s disc

ABSTRACT This work explores the mixing rate of metals in the interstellar medium (ISM), comparing observational constraints from our solar neighbourhood to high resolution cosmological hydrodynamical simulations of Milky Way (MW)-like galaxies. The mixing rate, described by the coefficient C, is varied in simulations between 0 and 0.05, with resultant simulated galaxies compared to observations of metallicity dispersion in young star clusters, H ii regions and neutral gas in the disc of the MW. A value of C between 0.003125 and 0.0125 is found to self-consistently match a range of observables, with a best estimate of C = 0.0064 ± 0.0004. We demonstrate that the relationship between metal dispersion in young stars, H ii regions and neutral gas, versus the coefficient C, can be described by a power law. These constrained mixing rates infer a comparatively well-mixed ISM in the solar neighbourhood, at odds with some recent observations that have reported a highly inhomogeneous ISM. The degree of mixing suggested by this work is lower than what often employed in many hydrodynamical simulations. Our results have implications for studying the metallicity distribution of stars as well as of gas in the interstellar and circum-galactic media.

Momentum and energy injection by a wind-blown bubble into an inhomogeneous interstellar medium

Abstract We investigate the effect of mass-loading from embedded clouds on the evolution of wind-blown bubbles. We use 1D hydrodynamical calculations and assume that the clouds are numerous enough that they can be treated in the continuous limit, and that rapid mixing occurs so that the injected mass quickly merges with the global flow. The destruction of embedded clouds adds mass into the bubble, increasing its density. Mass-loading increases the temperature of the unshocked stellar wind due to the frictional drag, and reduces the temperature of the hot shocked gas as the available thermal energy is shared between more particles. Mass-loading may increase or decrease the volume-averaged bubble pressure. Mass-loaded bubbles are smaller, have less retained energy and lower radial momentum, but in all cases examined are still able to do significant PdV work on the swept-up gas. In this latter respect, the bubbles more closely resemble energy-conserving bubbles than the momentum-conserving-like behaviour of ‘quenched’ bubbles.

Parameterizing the Outflow from a Central Black Hole in Dwarf Spheroidal Galaxies: A 3D Hydrodynamic Simulation

Abstract Large galaxies harbor massive central black holes and their feedback exerts a substantial impact on their evolution. Recently, observations have suggested that dwarf galaxies might host black holes in their centers, but with lower masses (intermediate-mass black holes—IMBHs). The impact of IMBHs on the evolution of dwarf spheroidal galaxies (dSphs), however, has so far not been properly analyzed. In this work, we investigate the effects of an outflow from an IMBH on gas dynamics in dSphs by means of noncosmological, three-dimensional hydrodynamic simulations, letting the galactic gas distribution evolve over 3 Gyr under the influence of the IMBH’s outflow and supernova feedback. All simulations have a numerical resolution of 20.0 pc cell−1. Two scenarios are considered to infer differences in the propagation of the outflow, one with a homogeneous interstellar medium (ISM) and another one with inhomogeneities caused by supernova feedback. A minimal initial speed and a minimal initial density are required for the outflow to propagate, with the values depending on the conditions of the medium. In an unperturbed medium, the outflow propagates freely in both directions with the same velocity (lower than the initial one), removing a small fraction of gas from the galaxy (the exact fraction depends on the initial physical conditions of the outflow). However, in an inhomogeneous ISM, the impact of the outflow is substantially reduced, and its contribution to the removal of gas from the galaxy is almost negligible.

Discovery of O stars in the tidal Magellanic Bridge

The Magellanic Bridge, stretching between the Small and the Large Magellanic Cloud (SMC and LMC), is the nearest tidally stripped intergalactic environment. The Bridge has a significantly low average metallicity of Z ≲ 0.1 Z⊙. Here we report the first discovery of O-type stars in the Magellanic Bridge. Three massive O stars were identified thanks to the archival spectra obtained by the ESO’s Very Large Telescope FLAMES instrument. We analyze the spectra of each star using the Potsdam Wolf-Rayet (PoWR) non-local thermodynamic equilibrium model atmosphere code, which provides the physical parameters, ionizing photon fluxes, and surface abundances. The ages of the newly discovered O stars suggest that star formation in the Bridge is ongoing. Furthermore, the discovery of O stars in the Bridge implies that tidally stripped galactic tails containing low-density but highly dynamical gas are capable of producing massive O stars. The multi-epoch spectra indicate that all three O stars are binaries. Despite their spatial proximity to one another, these O stars are chemically distinct. One of them is a fast-rotating giant with nearly LMC-like abundances. The other two are main-sequence stars that rotate extremely slowly and are strongly metal depleted. We discover the most nitrogen-poor O star known to date. Taking into account the previous analyses of B stars in the Bridge, we interpret the various metal abundances as the signature of a chemically inhomogeneous interstellar medium (ISM), suggesting that the Bridge gas might have accreted during multiple episodes of tidal interaction between the Clouds. Attributing the lowest derived metal content to the primordial gas, the time of the initial formation of the Bridge may date back several billion years. Using the Gaia and Galex color-magnitude diagrams, we roughly estimate the total number of O stars in the Bridge and their total ionizing radiation. Comparing this with the energetics of the diffuse ISM, we find that the contribution of the hot stars to the ionizing radiation field in the Bridge is less than 10% and conclude that the main sources of ionizing photons are leaks from the LMC and SMC. We estimate a lower limit for the fraction of ionizing radiation that escapes from these two dwarf galaxies.

Numerical Simulations of Instability in the Shell of a Supernova Remnant Expanding in a Weakly in Inhomogeneous Interstellar Medium

Astronomical observations show that the supernova remnants, even with a close to spherical shape, usually have multiscale ripple-like distortions. For example 15 bends on the shock front are clearly visible in the remnant 0509-67.5. The global instability of the flow is considered as one of the possible mechanisms for generating such structures. In the frame of linear analysis [26] was shown that this instability has a resonance character. It means that the perturbations with a certain wavelength number should grow faster, therefore ripples in the remnant shell will manifest itself predominantly in a certain range of scales. In this paper we present the results of numerical simulations of the nonlinear stage of this instability, caused by small perturbations in the external environment, depending on their scale and intensity. The unpertubed gas is supposed to has a power-law spartial dependence ρ0(r) ~ r-ω, where ω is a constant. The blast wave generated by a supernova expolosion is descibed by a Sedov type similarity solution. We have developed two-dimensional numerical model of adiabatic flow with a blast wave in a comoving frame of reference based on parallel code AstroChemHydro [1]. It was shown that, according to the predictions of linear analysis, perturbations in the external flow amplify behind the front of the shock wave, which leads to the development of convective instability and the development of turbulence. The results of numerical simulations demonstrated that in shell-type flows (for omega = 2,7 and gamma = 4/3) external disturbances along with the characteristic rearrangement of the shock front and turbulization of the flow behind it, cause the formation of radially elongated filaments with a vortex structure behind the shock, the number of which is determined by the harmonic number of the perturbation l.

Nonlinear diffusion of cosmic rays escaping from supernova remnants: Cold partially neutral atomic and molecular phases

Aims. We aim to elucidate cosmic ray (CR) propagation in the weakly ionized environments of supernova remnants (SNRs) basing our analysis on the cosmic ray cloud (CRC) model. Methods. We solved two transport equations simultaneously: one for the CR pressure and one for the Alfvén wave energy density where CRs are initially confined in the SNR shock. Cosmic rays trigger a streaming instability and produce slab-type resonant Alfvén modes. The self-generated turbulence is damped by ion-neutral collisions and by noncorrelated interaction with Alfvén modes generated at large scales. Results. We show that CRs leaking in cold dense phases such as those found in cold neutral medium (CNM) and diffuse molecular medium (DiM) can still be confined over distances of a few tens of parsecs from the CRC center for a few thousand years. At 10 TeV, CR diffusion can be suppressed by two or three orders of magnitude. This effect results from a reduced ion-neutral collision damping in the decoupled regime. We calculate the grammage of CRs in these environments. We find that in both single and multi-phase setups at 10 GeV, CNM and DiM media can produce grammage in the range 10–20 g cm−2 in the CNM and DiM phases. At 10 TeV, because of nonlinear propagation the grammage increases to values in the range 0.5–20 g cm−2 in these two phases. We also present preliminary calculations in inhomogeneous interstellar medium combining two or three different phases where we obtain the same trends.

Momentum and energy injection by a supernova remnant into an inhomogeneous medium

ABSTRACT We investigate the effect of mass-loading from embedded clouds on the evolution of supernova remnants and on the energy and momentum that they inject into an inhomogeneous interstellar medium. We use 1D hydrodynamical calculations and assume that the clouds are numerous enough that they can be treated in the continuous limit. The destruction of embedded clouds adds mass into the remnant, increasing its density and pressure, and decreasing its temperature. The remnant cools more quickly, is less able to do PdV work on the swept-up gas, and ultimately attains a lower final momentum (by up to a factor of two or more). We thus find that the injection of momentum is more sensitive to an inhomogeneous environment than previous work has suggested, and we provide fits to our results for the situation where the cloud mass is not limited. The behaviour of the remnant is more complex in situations where the cloud mass is finite and locally runs out. In the case of multiple supernovae in a clustered environment, later supernova explosions may encounter higher densities than previous explosions due to the prior liberation of mass from engulfed clouds. If the cloud mass is finite, later explosions may be able to create a sustained hot phase when earlier explosions have not been able to.

Probing cosmic dawn with emission lines: predicting infrared and nebular line emission for ALMA and JWST

ABSTRACT Infrared and nebular lines provide some of our best probes of the physics regulating the properties of the interstellar medium (ISM) at high redshift. However, interpreting the physical conditions of high-redshift galaxies directly from emission lines remains complicated due to inhomogeneities in temperature, density, metallicity, ionization parameter, and spectral hardness. We present a new suite of cosmological, radiation-hydrodynamics simulations, each centred on a massive Lyman-break galaxy that resolves such properties in an inhomogeneous ISM. Many of the simulated systems exhibit transient but well-defined gaseous discs that appear as velocity gradients in [C ii] 157.6 $\mu$m emission. Spatial and spectral offsets between [C ii] 157.6 $\mu$m and [O iii] 88.33 $\mu$m are common, but not ubiquitous, as each line probes a different phase of the ISM. These systems fall on the local [C ii]–SFR relation, consistent with newer observations that question previously observed [C ii] 157.6 $\mu$m deficits. Our galaxies are consistent with the nebular line properties of observed z ∼ 2–3 galaxies and reproduce offsets on the BPT and mass-excitation diagrams compared to local galaxies due to higher star formation rate (SFR), excitation, and specific-SFR, as well as harder spectra from young, metal-poor binaries. We predict that local calibrations between H α and [O ii] 3727$\, \mathring{\rm A}$ luminosity and galaxy SFR apply up to z > 10, as do the local relations between certain strong line diagnostics (R23 and [O iii] 5007$\, \mathring{\rm A}$/H β) and galaxy metallicity. Our new simulations are well suited to interpret the observations of line emission from current (ALMA and HST) and upcoming facilities (JWST and ngVLA).

Modelling the bow shock Pulsar Wind Nebulae propagating through a non-uniform ISM

Many pulsars propagate through the interstellar medium (ISM) with supersonic velocities, and their pulsar winds interact with the interstellar medium (ISM), forming bow shocks and magnetotails (PWN). We model the propagation of pulsars through the inhomogeneous ISM using non-relativistic axisymmetric magneto-hydrodynamic (MHD) simulations. We take into account the wind from the star, and the azimuthal and poloidal components of the magnetic field, and investigate the PWN at different levels of magnetization (the ratio of magnetic to matter energy-densities) in the wind. We consider the interaction of PWN with small-scale and large-scale imhomogeneities in the ISM at different values of magnetization. We conclude that the inhomogeneities in the ISM can change the shapes of the bow shocks and magnetotails at different values of the magnetization.We compare the results of our simulations with the images of the Guitar Nebula and other PWN that show irregularities in the shapes of their bow shocks and magnetotails. We conclude that these irregularities may be caused by the interaction of PWN with the inhomogeneities in the ISM.

Enhancement of Feedback Efficiency by Active Galactic Nucleus Outflows via the Magnetic Tension Force in the Inhomogeneous Interstellar Medium

Abstract By performing three-dimensional magnetohydrodynamics simulations of subrelativistic jets and disk winds propagating into the magnetized inhomogeneous interstellar medium (ISM), we investigate the magnetic effects on the active galactic nucleus feedback. Our simulations reveal that the magnetic tension force promotes the acceleration of the dense gas clouds, since the magnetic field lines, which are initially straight, bend around the gas clouds. In the jet models, the velocity dispersion of the clouds increases with an increase in the initial magnetic fields. The increment of the kinetic energy of the clouds is proportional to the initial magnetic fields, implying that the magnetic tension force increases the energy conversion efficiency from the jet to the gas clouds. Through simulations of the mildly collimated disk wind and the funnel-shaped disk wind, we confirm that such an enhancement of the energy conversion efficiency via the magnetic fields appears even if the energy is injected via the disk winds. The enhancement of the acceleration of the dense part of the magnetized ISM via the magnetic tension force will occur wherever the magnetized inhomogeneous matter is blown away.

MODELING DUST EVOLUTION IN GALAXIES WITH A MULTIPHASE, INHOMOGENEOUS ISM

ABSTRACT We develop a model of dust evolution in a multiphase, inhomogeneous interstellar medium (ISM) using hydrodynamical simulations of giant molecular clouds in a Milky Way–like spiral galaxy. We improve the treatment of dust growth by accretion in the ISM to investigate the role of the temperature-dependent sticking coefficient and ion–grain interactions. From detailed observational data on the gas-phase Si abundances measured in the local Galaxy, we derive a relation between the average and the local gas density that we use as a critical constraint for the models. This relation requires a sticking coefficient that decreases with the gas temperature. The relation predicted by the models reproduces the slope of −0.5 for the observed relation in cold clouds, which is steeper than that for the warm medium and is explained by dust growth. We find that growth occurs in the cold medium for all adopted values of the minimum grain size a min from 1 to 5 nm. For the classical cutoff of , the Coulomb repulsion results in slower accretion and higher than the observed values. For , the Coulomb interactions enhance the growth rate, steepen the slope of the – relation, and provide a better match to observations. The rates of dust re-formation in the ISM by far exceed the rates of dust production by stellar sources. After the initial 140 Myr, the cycle of matter in and out of dust reaches a steady state, in which the dust growth balances the destruction on a similar timescale of 350 Myr.

MOMENTUM INJECTION BY SUPERNOVAE IN THE INTERSTELLAR MEDIUM

Supernova (SN) explosions deposit prodigious energy and momentum in their environments, with the former regulating multiphase thermal structure and the latter regulating turbulence and star formation rates in the interstellar medium (ISM). However, systematic studies quantifying the impact of SNe in realistic inhomogeneous ISM conditions have been lacking. Using three-dimensional hydrodynamic simulations, we investigate the dependence of radial momentum injection on both physical conditions (considering a range of mean density n=0.1-100) and numerical parameters. Our inhomogeneous simulations adopt two-phase background states that result from thermal instability in atomic gas. Although the SNR morphology becomes highly complex for inhomogeneous backgrounds, the radial momentum injection is remarkably insensitive to environmental details. For our two-phase simulations, the final momentum produced by a single SN is given by 2.8*10^5 M_sun*km/s n^{-0.17}. This is only 5% less than the momentum injection for a homogeneous environment with the same mean density, and only 30% greater than the momentum at the time of shell formation. The maximum mass in hot gas is quite insensitive to environmental inhomogeneity. Strong magnetic fields alter the hot gas mass at very late times, but the momentum injection remains the same. Initial experiments with multiple spatially-correlated SNe show a momentum per event nearly as large as single-SN cases. We also present a full numerical parameter study to assess convergence requirements. For convergence in the momentum and other quantities, we find that the numerical resolution dx and the initial size of the SNR r_init must satisfy dx, r_init<r_sf/3, where the shell formation radius is given by r_sf = 30 pc n^{-0.46} for two-phase models (or 30% smaller for a homogeneous medium).

AGN feedback by relativistic jets

AbstractFeedback provided by relativistic jets may be effective in shaping the galaxy luminosity function. The quenching mode (quasar mode) at redshifts ~2-3 potentially disperses gas in star-forming galaxies. The maintenance mode (radio mode) heats the gas in galaxy clusters counteracting cooling flows. A number of authors have examined the effect of relativistic jets in dispersing clouds in the kpc-scale inhomogeneous interstellar medium of evolving galaxies. We have also investigated a particular case of maintenance-mode feedback in our simulation of the iconic radio galaxy / cooling flow cluster Hydra A. Modelling of the knots produced by the jets in the inner 10 kpc provides an estimate of 0.8 – 0.9 c for the velocities of the jets in agreement with other velocity estimates for FR1 jets. The addition of jet precession provides realistic simulations of the morphology of the Hydra A radio source and raises interesting questions as to the role of black hole and disk precession, in general, in galaxy formation.

XMM-Newtonobservation of the Galactic supernova remnant W51C (G49.1–0.1)

The supernova remnant (SNR) W51C is a Galactic object located in a strongly inhomogeneous interstellar medium with signs of an interaction of the SNR blast wave with dense molecular gas. Diffuse X-ray emission from the interior of the SNR can reveal element abundances in the different emission regions and shed light on the type of supernova (SN) explosion and its progenitor. The hard X-ray emission helps to identify possible candidates for a pulsar formed in the SN explosion and for its pulsar wind nebula (PWN). We have analysed X-ray data obtained with XMM-Newton. Spectral analyses in selected regions were performed. Ejecta emission in the bright western part of the SNR, located next to a complex of dense molecular gas, was confirmed. The Ne and Mg abundances suggest a massive progenitor with a mass of > 20 M_sun. Two extended regions emitting hard X-rays were identified (corresponding to the known sources [KLS2002] HX3 west and CXO J192318.5+140305 discovered with ASCA and Chandra, respectively), each of which has an additional point source inside and shows a power-law spectrum with Gamma ~ 1.8. Based on their X-ray emission, both sources can be classified as PWN candidates.

Lymanαline and continuum radiative transfer in a clumpy interstellar medium

Aims: We aim to study the effects of an inhomogeneous interstellar medium (ISM) on the strength and the shape of the Lyman alpha (Lya) line in starburst galaxies. Methods: Using our 3D Monte Carlo Lya radiation transfer code, we study the radiative transfer of Lya, UV and optical continuum photons in homogeneous and clumpy shells of neutral hydrogen and dust surrounding a central source. Our simulations predict the Lya and continuum escape fraction, the Lya equivalent width EW(Lya), the Lya line profile and their dependence on the gas geometry and the main input physical parameters. Results: The ISM clumpiness is found to have a strong impact on the Lya line radiative transfer, entailing a strong dependence of the emergent features of the Lya line (escape fraction, EW(Lya)) on the ISM morphology. Although a clumpy and dusty ISM appears more transparent to radiation (both line and continuum) compared to an equivalent homogeneous ISM of equal dust optical depth, we find that the Lya photons are, in general, still more attenuated than UV continuum radiation. As a consequence, the observed Lya equivalent width (EWobs(Lya)) is lower than the intrinsic one (EWint(Lya)) for nearly all clumpy ISM configurations considered. There are, however, special conditions under which Lya photons escape more easily than the continuum, resulting in an enhanced EWobs(Lya). The requirement for this to happen is that the ISM is almost static (galactic outflows < 200 km/s), extremely clumpy (with density contrasts >10^7 in HI between clumps and the interclump medium), and very dusty (E(B-V) > 0.30). When these conditions are fulfilled the emergent Lya line profile shows no velocity shift and little asymmetry. Given the asymmetry and velocity shifts generally observed in starburst galaxies with Lya emission, we conclude that clumping is unlikely to significantly enhance their relative Lya/UV transmission.

ASYMMETRIES OF AGN JETS IN INHOMOGENEOUS MEDIA

Considerable asymmetries in jets from active galactic nuclei (AGN) and associated double radio sources can be caused by an inhomogeneous interstellar medium of the host galaxy. These asymmetries can easily be estimated by 1D propagation models, but hydrodynamical simulations have shown that the actual asymmetries can be considerably larger. With a set of smaller-scale hydrodynamical simulations we examine these asymmetries, and find they are typically a factor of ~ 3 larger than in 1D models. We conclude that, at high redshift, large asymmetries in radio sources are expected in gas-rich galaxies with a clumpy interstellar medium.

Three-dimensional modeling of ionized gas

The first generation of stars that formed directly from the primordial gas played a crucial role in the early phase of the reionization of the universe. Because of the short lifetimes of these stars the metals produced in their cores were quickly returned to the environment, from which early PopII stars with a different initial mass function and different SEDs were formed, already much earlier than the time at which the universe became completely reionized at a redshift of z~6. Using a state-of-the-art model atmosphere code we calculate realistic SEDs of very massive stars (VMSs) of different metallicities to serve as input for the 3-dimensional multi-frequency radiative transfer code we have developed to simulate the temporal evolution of the ionization of the inhomogeneous interstellar and intergalactic medium, using multiple stellar clusters as sources of ionizing radiation. Our tool handles distributions of numerous radiative sources characterized by high resolution synthetic SEDs, and yields occupation numbers of the required energy levels of the most important elements which are treated in NLTE consistently with the 3d radiative transfer. We further demonstrate that the increasing metallicity of the radiative sources in the transition from PopIII stars to PopII stars has a strong impact on the hardness of the emitted spectrum, and hence on the reionization history of helium. A top-heavy stellar mass distribution characterized by VMSs forming in chemically evolved clusters of high core mass density may not only provide the progenitors of intermediate-mass and supermassive black holes (SMBHs), but also play an important role for the reionization of HeII. The number of VMSs required to reionize HeII by a redshift of z~2.5 is astonishingly close to the number of VMSs required to explain galactic SMBHs if one assumes that these have been formed by mergers of smaller black holes.

Active galactic nucleus feedback and triggering of star formation in galaxies

Feedback from the central black hole in active galactic nuclei (AGN) may be responsible for establishing the observed MBH-sigma relation and limiting the bulge stellar mass of the host galaxy. Here we explore the possibility of AGN feedback triggering star formation in the host galaxy. We consider a shell of dusty gas, driven outwards by radiation pressure, and analyse its escape/trapping condition in the galactic halo for different underlying dark matter potentials. In the isothermal potential, we obtain that the standard condition setting the observed MBH-sigma relation is not sufficient to clear gas out of the entire galaxy; whereas the same condition is formally sufficient in the case of the Hernquist and Navarro-Frenk-White profiles. The squeezing and compression of the inhomogeneous interstellar medium during the ejection process can trigger star formation within the feedback-driven shell. We estimate the resulting star formation rate and total additional stellar mass. In this picture, new stars are formed at increasingly larger radii and successively populate the outer regions of the host galaxy. This characteristic pattern may be compared with the observed 'inside-out' growth of massive galaxies.

Turbulent magnetic field amplification driven by cosmic ray pressure gradients

Observations of non-thermal emission from several supernova remnants suggest that magnetic fields close to the blastwave are much stronger than would be naively expected from simple shock compression of the field permeating the interstellar medium (ISM). We present a simple model which is capable of achieving sufficient magnetic field amplification to explain the observations. We propose that the cosmic-ray pressure gradient acting on the inhomogeneous ISM upstream of the supernova blastwave induces strong turbulence upstream of the supernova blastwave. The turbulence is generated through the differential acceleration of the upstream ISM which occurs as a result of density inhomogeneities in the ISM. This turbulence then amplifies the pre-existing magnetic field. Numerical simulations are presented which demonstrate that amplification factors of 20 or more are easily achievable by this mechanism when reasonable parameters for the ISM and supernova blastwave are assumed. The length scale over which this amplification occurs is that of the diffusion length of the highest energy non-thermal particles.