Source author record

Ellen G. Zweibel

Ellen G. Zweibel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2022arXiv

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- IV: Program overview and first results on the polarization fraction

We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation and sample selection of the program. Here, we focus on the general trends in the emissive polarization fraction. Far-infrared polarimetric observations trace the thermal polarized emission of magnetically aligned dust grains across the galaxy disks with polarization fractions of $P=0-15$% in the cold, $T_{\rm d} = [19,48]$ K, and dense, $\log_{10}(N_{\rm HI+H_{2}}) = [19.96,22.91]$, interstellar medium. The spiral galaxies show a median $\langle P_{154μm} \rangle = 3.3\pm0.9 $% across the disks. We report the first polarized spectrum of starburst galaxies showing a minimum within $89-154$ $μ$m. The falling $53-154$ $μ$m polarized spectrum may be due to a decrease in the dust grain alignment efficiency produced by variations in dust temperatures along the line-of-sight in the galactic outflow. We find that the starburst galaxies and the star-forming regions within normal galaxies have the lowest polarization fractions. We find that 50% (7 out of 14) of the galaxies require a broken power-law in the $P-N_{HI+H_{2}}$ and $P-T_{d}$ relations with three different trends. Group 1 has a relative increase of anisotropic random B-fields produced by compression or shear of B-fields in the galactic outflows, starburst rings, and inner-bar of galaxies; and Groups 2 and 3 have a relative increase of isotropic random B-fields driven by star-forming regions in the spiral arms, and/or an increase of dust grain alignment efficiency caused by shock-driven regions or evolutionary stages of a galaxy.

preprint2020arXiv

Cosmic Ray Driven Outflows from the Large Magellanic Cloud: Contributions to the LMC Filament

In this paper, we build from previous work (Bustard et al. 2018) and present simulations of recent (within the past Gyr), magnetized, cosmic ray driven outflows from the Large Magellanic Cloud (LMC), including our first attempts to explicitly use the derived star formation history of the LMC to seed outflow generation. We run a parameter set of simulations for different LMC gas masses and cosmic ray transport treatments, and we make preliminary comparisons to published outflow flux estimates, neutral and ionized hydrogen observations, and Faraday rotation measure maps. We additionally report on the gas mass that becomes unbound from the LMC disk and swept by ram pressure into the Trailing Magellanic Stream. We find that, even for our largest outburst, the mass contribution to the Stream is still quite small, as much of the outflow-turned-halo gas is shielded on the LMC's far-side due to the LMC's primarily face-on infall through the Milky Way halo over the past Gyr. On the LMC's near-side, past outflows have fought an uphill battle against ram pressure, with near-side halo mass being at least a factor of a few smaller than the far-side. Absorption line studies probing only the LMC foreground, then, may be severely underestimating the total mass of the LMC halo formed by outflows.

preprint2020arXiv

Small-scale structure traced by neutral hydrogen absorption in the direction of multiple-component radio continuum sources

We have studied the small scale distribution of atomic hydrogen (HI) using 21-cm absorption spectra against multiple-component background radio continuum sources from the 21-SPONGE survey and the Millennium Arecibo Absorption Line Survey. We have found $>5σ$ optical depth variations at a level of $\sim0.03-0.5$ between 13 out of 14 adjacent sightlines separated by a few arcseconds to a few arcminutes, suggesting the presence of neutral structures on spatial scales from a few to thousands of AU (which we refer to as tiny scale atomic structure, TSAS). The optical depth variations are strongest in directions where the HI column density and the fraction of HI in the cold neutral medium (CNM) are highest, which tend to be at low Galactic latitudes. By measuring changes in the properties of Gaussian components fitted to the absorption spectra, we find that changes in both the peak optical depth and the linewidth of TSAS absorption features contribute to the observed optical depth variations, while changes in the central velocity do not appear to strongly impact the observed variations. Both thermal and turbulent motions contribute appreciably to the linewidths, but the turbulence does not appear strong enough to confine overpressured TSAS. In a majority of cases, the TSAS column densities are sufficiently high that these structures can radiatively cool fast enough to maintain thermal equilibrium with their surroundings, even if they are overpressured. We also find that a majority of TSAS is associated with the CNM. For TSAS in the direction of the Taurus molecular cloud and the local Leo cold cloud, we estimate densities over an order of magnitude higher than typical CNM densities.

preprint2019arXiv

A Dynamical Study of Extraplanar Diffuse Ionized Gas in NGC 5775

The structure and kinematics of gaseous, disk-halo interfaces are imprinted with the processes that transfer mass, metals, and energy between galactic disks and their environments. We study the extraplanar diffuse ionized gas (eDIG) layer in the interacting, star-forming galaxy NGC 5775 to better understand the consequences of star-formation feedback on the dynamical state of the thick-disk interstellar medium (ISM). Combining emission-line spectroscopy from the Robert Stobie Spectrograph on the Southern African Large Telescope with radio continuum observations from Continuum Halos in Nearby Galaxies - an EVLA Survey, we ask whether thermal, turbulent, magnetic field, and cosmic-ray pressure gradients can stably support the eDIG layer in dynamical equilibrium. This model fails to reproduce the observed exponential electron scale heights of the eDIG thick disk and halo on the northeast ($h_{z,e} = 0.6, 7.5$ kpc) and southwest ($h_{z,e} = 0.8, 3.6$ kpc) sides of the galaxy at $R < 11$ kpc. We report the first definitive detection of an increasing eDIG velocity dispersion as a function of height above the disk. Blueshifted gas along the minor axis at large distances from the midplane hints at a disk-halo circulation and/or ram pressure effects caused by the ongoing interaction with NGC 5774. This work motivates further integral field unit and/or Fabry-Perot spectroscopy of galaxies with a range of star-formation rates to develop a spatially-resolved understanding of the role of star-formation feedback in shaping the kinematics of the disk-halo interface.

preprint2019arXiv

The Role of Pressure Anisotropy in Cosmic Ray Hydrodynamics

Cosmic ray propagation in the Milky Way and other galaxies is largely diffusive, with mean free path determined primarily by pitch angle scattering from hydromagnetic waves with wavelength of order the cosmic ray gyroradius. In the theory of cosmic ray self confinement, the waves are generated by instabilities driven by the cosmic rays themselves. The dominant instability is due to bulk motion, or streaming, of the cosmic rays, parallel to the background magnetic field B, and transfers cosmic ray momentum and energy to the thermal gas as well as confining the cosmic rays. Classical arguments and recent numerical simulations show that self confinement due to the streaming instability breaks down unless the cosmic ray pressure and thermal gas density gradients parallel to B are aligned, a condition which is unlikely to always be satisfied We investigate an alternative mechanism for cosmic ray self confinement and heating of thermal gas based on pressure anisotropy instability. Although pressure anisotropy is demonstrably less effective than streaming instability as a self confinement and heating mechanism on global scales, it may be important on mesoscales, particularly near sites of cosmic ray injection.

preprint2016arXiv

A Versatile Family of Galactic Wind Models

We present a versatile family of model galactic outflows including non-uniform mass and energy source distributions, a gravitational potential from an extended mass source, and radiative losses. The model easily produces steady-state wind solutions for a range of mass-loading factors, energy-loading factors, galaxy mass and galaxy radius. We find that, with radiative losses included, highly mass-loaded winds must be driven at high central temperatures, whereas low mass-loaded winds can be driven at low temperatures just above the peak of the cooling curve, meaning radiative losses can drastically affect the wind solution even for low mass-loading factors. By including radiative losses, we are able to show that subsonic flows can be ignored as a possible mechanism for expelling mass and energy from a galaxy compared to the more efficient transonic solutions. Specifically, the transonic solutions with low mass-loading and high energy-loading are the most efficient. Our model also produces low-temperature, high-velocity winds that could explain the prevalence of low-temperature material in observed outflows. Finally, we show that our model, unlike the well-known Chevalier & Clegg (1985) model, can reproduce the observed linear relationship between wind X-ray luminosity and star formation rate (SFR) over a large range of SFR from $1-1000$ M$_{\odot}$/yr assuming the wind mass-loading factor is higher for low-mass, and hence, low-SFR galaxies. We also constrain the allowed mass-loading factors that can fit the observed X-ray luminosity vs. SFR trend, further suggesting an inverse relationship between mass-loading and SFR as explored in advanced numerical simulations.

preprint2016arXiv

Detection of a ~20 kpc coherent magnetic field in the outskirt of merging spirals: the Antennae galaxies

We present a study of the magnetic field properties of NGC 4038/9 (the `Antennae' galaxies), the closest example of a late stage merger of two spiral galaxies. Wideband polarimetric observations were performed using the Karl G. Jansky Very Large Array between 2 and 4 GHz. Rotation measure synthesis and Faraday depolarization analysis was performed to probe the magnetic field strength and structure at spatial resolution of $\sim1$ kpc. Highly polarized emission from the southern tidal tail is detected with intrinsic fractional polarization close to the theoretical maximum ($0.62\pm0.18$), estimated by fitting the Faraday depolarization with a volume that is both synchrotron emitting and Faraday rotating containing random magnetic fields. Magnetic fields are well aligned along the tidal tail and the Faraday depths shows large-scale smooth variations preserving its sign. This suggests the field in the plane of the sky to be regular up to $\sim20$ kpc, which is the largest detected regular field structure on galactic scales. The equipartition field strength of $\sim8.5~μ$G of the regular field in the tidal tail is reached within a few 100 Myr, likely generated by stretching of the galactic disc field by a factor of 4--9 during the tidal interaction. The regular field strength is greater than the turbulent fields in the tidal tail. Our study comprehensively demonstrates, although the magnetic fields within the merging bodies are dominated by strong turbulent magnetic fields of $\sim20~μ$G in strength, tidal interactions can produce large-scale regular field structure in the outskirts.

preprint2015arXiv

Cosmic Rays, Gamma-Rays, and Neutrinos in the Starburst Nuclei of Arp 220

The cores of Arp 220, the closest ultra-luminous infrared starburst galaxy, provide an opportunity to study interactions of cosmic rays under extreme conditions. In this paper, we model the populations of cosmic rays produced by supernovae in the central molecular zones of both starburst nuclei. We find that ~65 - 100% of cosmic rays are absorbed in these regions due to their huge molecular gas contents, and thus, the nuclei of Arp 220 nearly complete proton calorimeters. As the cosmic ray protons collide with the interstellar medium, they produce secondary electrons that are also contained within the system and radiate synchrotron emission. Using results from chi-squared tests between the model and the observed radio spectral energy distribution, we predict the emergent gamma-ray and high-energy neutrino spectra and find the magnetic field to be at milligauss levels. Because of the extremely intense far-infrared radiation fields, the gamma-ray spectrum steepens significantly at TeV energies due to gamma-gamma absorption.

preprint2015arXiv

Effects of Line-tying on Resistive Tearing Instability in Slab Geometry

The effects of line-tying on resistive tearing instability in slab geometry is studied within the framework of reduced magnetohydrodynamics (RMHD).\citep{KadomtsevP1974,Strauss1976} It is found that line-tying has a stabilizing effect. The tearing mode is stabilized when the system length $L$ is shorter than a critical length $L_{c}$, which is independent of the resistivity $η$. When $L$ is not too much longer than $L_{c}$, the growthrate $γ$ is proportional to $η$ . When $L$ is sufficiently long, the tearing mode scaling $γ\simη^{3/5}$ is recovered. The transition from $γ\simη$ to $γ\simη^{3/5}$ occurs at a transition length $L_{t}\simη^{-2/5}$.

preprint2015arXiv

Equipartition and Cosmic Ray Energy Densities in Central Molecular Zones of Starbursts

The energy densities in magnetic fields and cosmic rays (CRs) in galaxies are often assumed to be in equipartition, allowing for an indirect estimate of the magnetic field strength from the observed radio synchrotron spectrum. However, both primary and secondary CRs contribute to the synchrotron spectrum, and the CR electrons also loose energy via bremsstrahlung and inverse Compton. While classical equipartition formulae avoid these intricacies, there have been recent revisions that account for the extreme conditions in starbursts. Yet, the application of the equipartition formula to starburst environments also presupposes that timescales are long enough to reach equilibrium. Here, we test equipartition in the central molecular zones (CMZs) of nearby starburst galaxies by modeling the observed gamma-ray spectra, which provide a direct measure of the CR energy density, and the radio spectra, which provide a probe of the magnetic field strength. We find that in starbursts, the magnetic field energy density is significantly larger than the CR energy density, demonstrating that the equipartition argument is frequently invalid for CMZs.

preprint2014arXiv

Conduction in low Mach number flows: Part I Linear & weakly nonlinear regimes

Thermal conduction is an important energy transfer and damping mechanism in astrophysical flows. Fourier's law - the heat flux is proportional to the negative temperature gradient, leading to temperature diffusion - is a well-known empirical model of thermal conduction. However, entropy diffusion has emerged as an alternative thermal conduction model, despite not ensuring the monotonicity of entropy. This paper investigates the differences between temperature and entropy diffusion for both linear internal gravity waves and weakly nonlinear convection. In addition to simulating the two thermal conduction models with the fully compressible Navier-Stokes equations, we also study their effects in the reduced, "sound-proof" anelastic and pseudo-incompressible equations. We find that in the linear and weakly nonlinear regimes, temperature and entropy diffusion give quantitatively similar results, although there are some larger errors in the pseudo-incompressible equations with temperature diffusion due to inaccuracies in the equation of state. Extrapolating our weakly nonlinear results, we speculate that differences between temperature and entropy diffusion might become more important for strongly turbulent convection.

preprint2014arXiv

The Cosmic Ray Population of the Galactic Central Molecular Zone

The conditions in the Galactic Center are often compared with those in starburst systems, which contain higher supernova rates, stronger magnetic fields, more intense radiation fields, and larger amounts of dense molecular gas than in our own Galactic disk. Interactions between such an augmented interstellar medium and cosmic rays result in brighter radio and gamma-ray emission. Here, we test how well the comparisons between the Galactic Center and starburst galaxies hold by applying a model for cosmic ray interactions to the Galactic Center to predict the resulting gamma-ray emission. The model only partially explains the observed gamma-ray and radio emission. The model for the gamma-ray spectrum agrees with the data at TeV energies but not at GeV energies. Additionally, as the fits of the model to the radio and gamma-ray spectra require significant differences in the optimal wind speed and magnetic field strength, we find that the single-zone model alone cannot account for the observed emission from the Galactic Center. Our model is improved by including a soft, additional cosmic-ray population. We assess such a cosmic ray population and its potential sources and find that a cosmic-ray electron spectrum is energetically favored over a cosmic-ray proton spectrum.

preprint2014arXiv

The Transport of Cosmic Rays Across Magnetic Fieldlines

The long residence times and small anisotropies of cosmic rays suggest that they are well confined and well scattered by the Galactic magnetic field. Due to the disklike shape of the confinement volume, transport in the vertical direction, perpendicular to the mean Galactic magnetic field, is key to cosmic ray escape. It has long been recognized that this vertical transport depends both on the vertical component of the fieldlines themselves and on the extent to which the cosmic rays are tied to the fieldlines. In this paper we use magnetic fields with very simple spatial and temporal structure to isolate some important features of cross fieldline transport. We show that even simple magnetic nonuniformities combined with pitch angle scattering can enhance cross fieldline transport by several orders of magnitude, while pitch angle scattering is unnecessary for enhanced transport if the field is chaotic. Nevertheless, perpendicular transport is much less than parallel transport in all the cases we study. We apply the results to confinement of cosmic rays in the Fermi Bubbles.

preprint2013arXiv

Active Galactic Nuclei, Neutrinos, and Interacting Cosmic Rays in NGC 253 & NGC 1068

The galaxies M82, NGC 253, NGC 1068, and NGC 4945 have been detected in gamma-rays by Fermi. Previously, we developed and tested a model for cosmic ray interactions in the starburst galaxy M82. Now, we aim to explore the differences between starburst and active galactic nuclei (AGN) environments by applying our self-consistent model to the starburst galaxy NGC 253 and the Seyfert galaxy NGC 1068. Assuming constant cosmic-ray acceleration efficiency by supernova remnants with Milky-Way parameters, we calculate the cosmic-ray proton and primary and secondary electron/positron populations, predict the radio and gamma-ray spectra, and compare with published measurements. We find that our models easily fits the observed gamma-ray spectrum for NGC 253 while constraining the cosmic ray source spectral index and acceleration efficiency. However, we encountered difficultly modeling the observed radio data and constraining the speed of the galactic wind and the magnetic field strength, unless the gas mass is less than currently preferred values. Additionally, our starburst model consistently underestimates the observed gamma-ray flux and overestimates the radio flux for NGC 1068; these issues would be resolved if the AGN is the primary source of gamma-rays. We discuss the implications of these results and make predictions for the neutrino fluxes for both galaxies.

preprint2013arXiv

Cosmic Ray Heating of the Warm Ionized Medium

Observations of line ratios in the Milky Way's warm ionized medium (WIM) suggest that photoionization is not the only heating mechanism present. For the additional heating to explain the discrepancy it would have to have a weaker dependence on the gas density than the cooling rate, $Λn_e^2$. \cite{reynolds99} suggested turbulent dissipation or magnetic field reconnection as possible heating sources. We investigate here the viability of MHD-wave mediated cosmic ray heating as a supplemental heating source. This heating rate depends on the gas density only through its linear dependence on the Alfvén speed, which goes as $n_e^{-1/2}$. We show that, scaled to appropriate values of cosmic ray energy density, cosmic ray heating can be significant. Furthermore, this heating is stable to perturbations. These results should also apply to warm ionized gas in other galaxies.

preprint2013arXiv

Cosmic Ray Sampling of a Clumpy Interstellar Medium

How cosmic rays sample the multi-phase interstellar medium (ISM) in starburst galaxies has important implications for many science goals, including evaluating the cosmic ray calorimeter model for these systems, predicting their neutrino fluxes, and modeling their winds. Here, we use Monte Carlo simulations to study cosmic ray sampling of a simple, two-phase ISM under conditions similar to those of the prototypical starburst galaxy M82. The assumption that cosmic rays sample the mean density of the ISM in the starburst region is assessed over a multi-dimensional parameter space where we vary the number of molecular clouds, the galactic wind speed, the extent to which the magnetic field is tangled, and the cosmic ray injection mechanism. We evaluate the ratio of the emissivity from pion production in molecular clouds to the emissivity that would be observed if the cosmic rays sampled the mean density, and seek areas of parameter space where this ratio differs significantly from unity. The assumption that cosmic rays sample the mean density holds over much of parameter space; however, this assumption begins to break down for high cloud density, injection close to the clouds, and a very tangled magnetic field. We conclude by evaluating the extent to which our simulated starburst region behaves as a proton calorimeter and constructing the time-dependent spectrum of a burst of cosmic rays.

preprint2013arXiv

Energy conservation and gravity waves in sound-proof treatments of stellar interiors: Part II Lagrangian constrained analysis

The speed of sound greatly exceeds typical flow velocities in many stellar and planetary interiors. To follow the slow evolution of subsonic motions, various sound-proof models attempt to remove fast acoustic waves whilst retaining stratified convection and buoyancy dynamics. In astrophysics, anelastic models typically receive the most attention in the class of sound-filtered stratified models. Generally, anelastic models remain valid in nearly adiabatically stratified regions like stellar convection zones, but may break down in strongly sub-adiabatic, stably stratified layers common in stellar radiative zones. However, studying stellar rotation, circulation, and dynamos requires understanding the complex coupling between convection and radiative zones, and this requires robust equations valid in both regimes. Here we extend the analysis of equation sets begun in Brown Vasil & Zweibel 2012, which studied anelastic models, to two types of pseudo-incompressible models. This class of models has received attention in atmospheric applications, and more recently in studies of white-dwarf supernovae progenitors. We demonstrate that one model conserves energy but the other does not. We use Lagrangian variational methods to extend the energy conserving model to a general equation of state, and dub the resulting equation set the Generalized Pseudo-Incompressible (GPI) model. We show that the GPI equations suitably capture low frequency phenomena in both convection and radiative zones in stars and other stratified systems, and we provide recommendations for converting low-Mach number codes to this equation set.

preprint2013arXiv

Lessons from Comparisons between the Nuclear Region of the Milky Way & Those in Nearby Spirals

The Milky Way appears is a typical barred spiral, and comparisons can be made between its nuclear region and those of structurally similar nearby spirals. Maffei 2, M83, IC 342 and NGC 253 are nearby systems whose nuclear region properties contrast with those of the Milky Way. Stellar masses derived from NIR photometery, molecular gas masses and star formation rates allow us to assess the evolutionary states of this set of nuclear regions. These data suggest similarities between nuclear regions in terms of their stellar content while highlighting significant differences in current star formation rates. In particular current star formation rates appear to cover a larger range than expected based on the molecular gas masses. This behavior is consistent with nuclear region star formation experiencing episodic variations. Under this hypothesis the Milky Way's nuclear region currently may be in a low star formation rate phase.

preprint2013arXiv

Winds, Clumps, and Interacting Cosmic Rays in M82

We construct a family of models for the evolution of energetic particles in the starburst galaxy M82 and compare them to observations to test the calorimeter assumption that all cosmic ray energy is radiated in the starburst region. Assuming constant cosmic ray acceleration efficiency with Milky Way parameters, we calculate the cosmic-ray proton and primary and secondary electron/positron populations as a function of energy. Cosmic rays are injected with Galactic energy distributions and electron-to-proton ratio via type II supernovae at the observed rate of 0.07/yr. From the cosmic ray spectra, we predict the radio synchrotron and γ-ray spectra. To more accurately model the radio spectrum, we incorporate a multiphase interstellar medium in the starburst region of M82. Our model interstellar medium is highly fragmented with compact dense molecular clouds and dense photoionized gas, both embedded in a hot, low density medium in overall pressure equilibrium. The spectra predicted by this one-zone model are compared to the observed radio and γ-ray spectra of M82. Chi-squared tests are used with radio and γ-ray observations and a range of model predictions to find the best-fit parameters. The best-fit model yields constraints on key parameters in the starburst zone of M82, including a magnetic field strength of ~250 μG and a wind advection speed in the range of 300-700 km/s. We find that M82 is a good electron calorimeter but not an ideal cosmic-ray proton calorimeter and discuss the implications of our results for the astrophysics of the far infrared-radio correlation in starburst galaxies.

preprint2011arXiv

Onset of Fast Magnetic Reconnection in Partially Ionized Gases

We consider quasi-stationary two-dimensional magnetic reconnection in a partially ionized incompressible plasma. We find that when the plasma is weakly ionized and the collisions between the ions and the neutral particles are significant, the transition to fast collisionless reconnection due to the Hall effect in the generalized Ohm's law is expected to occur at much lower values of the Lundquist number, as compared to a fully ionized plasma case. We estimate that these conditions for fast reconnection are satisfied in molecular clouds and in protostellar disks.

preprint2011arXiv

The Interaction of Cosmic Rays with Diffuse Clouds

We study the change in cosmic-ray pressure, the change in cosmic-ray density, and the level of cosmic-ray induced heating via Alfven-wave damping when cosmic rays move from a hot ionized plasma to a cool cloud embedded in that plasma. The general analysis method outlined here can apply to diffuse clouds in either the ionized interstellar medium or in galactic winds. We introduce a general-purpose model of cosmic-ray diffusion building upon the hydrodynamic approximation for cosmic rays (from McKenzie & Voelk and Breitschwerdt and collaborators). Our improved method self-consistently derives the cosmic-ray flux and diffusivity under the assumption that the streaming instability is the dominant mechanism for setting the cosmic-ray flux and diffusion. We find that, as expected, cosmic rays do not couple to gas within cool clouds (cosmic rays exert no forces inside of cool clouds), that the cosmic-ray density does not increase within clouds (it may slightly decrease in general, and decrease by an order of magnitude in some cases), and that cosmic-ray heating (via Alfven-wave damping and not collisional effects as for ~10 MeV cosmic rays) is only important under the conditions of relatively strong (10 micro-Gauss) magnetic fields or high cosmic-ray pressure (~10^{-11} ergs cm^{-3}).

preprint2010arXiv

Ambipolar Diffusion-Mediated Thermal Fronts in the Neutral ISM

In a thermally bistable medium, cold, dense gas is separated from warm, rareified gas by thin phase transition layers, or fronts, in which heating, radiative cooling, thermal conduction, and convection of material are balanced. We calculate the steady-state structure of such fronts in the presence of magnetic fields, including the processes of ion-neutral drift and ion-neutral frictional heating. We find that ambipolar diffusion efficiently transports the magnetic field across the fronts, leading to a flat magnetic field strength profile. The thermal profiles of such fronts are not significantly different from those of unmagnetized fronts. The near uniformity of the magnetic field strength across a front is consistent with the flat field strength-gas density relation that is observed in diffuse interstellar gas.

preprint2010arXiv

Effects of Line-tying on Magnetohydrodynamic Instabilities and Current Sheet Formation

An overview of some recent progress on magnetohydrodynamic stability and current sheet formation in a line-tied system is given. Key results on the linear stability of the ideal internal kink mode and resistive tearing mode are summarized. For nonlinear problems, a counterexample to the recent demonstration of current sheet formation by Low \emph{et al}. [B. C. Low and Å. M. Janse, Astrophys. J. \textbf{696}, 821 (2009)] is presented, and the governing equations for quasi-static evolution of a boundary driven, line-tied magnetic field are derived. Some open questions and possible strategies to resolve them are discussed.

preprint2010arXiv

Synchrotron Constraints on a Hybrid Cosmic-Ray and Thermally-Driven Galactic Wind

Cosmic rays and magnetic fields can substantially impact the launching of large-scale galactic winds. Many researchers have investigated the role of cosmic rays; our group previously showed that a cosmic-ray and thermally-driven wind could explain soft X-ray emission towards the center of the Galaxy. In this paper, we calculate the synchrotron emission from our original wind model and compare it to observations; the synchrotron data shows that earlier assumptions about the launching conditions of the wind must be changed: we are required to improve that earlier model by restricting the launching region to the domain of the inner "Molecular Ring", and by decreasing the magnetic field strength from the previously assumed maximum strength. With these physically-motived modifications, we find that a wind model can fit both the radio synchrotron and the X-ray emission, although that model is required to have a higher gas pressure and density than the previous model in order to reproduce the observed X-ray emission measure within the smaller `footprint'. The drop in magnetic field also decreases the effect of cosmic-ray heating, requiring a higher temperature at the base of the wind than the previous model.

preprint2010arXiv

The Role of the Magnetic Field in the Interstellar Medium of the Post-Starburst Dwarf Irregular Galaxy NGC 1569

(abridged) NGC 1569 is a nearby dwarf irregular galaxy which underwent an intense burst of star formation 10 to 40 Myr ago. We present observations that reach surface brightnesses two to eighty times fainter than previous radio continuum observations and the first radio continuum polarization observations. These observations allow us to probe the relationship of the magnetic field of NGC 1569 to the rest of its interstellar medium. We confirm the presence of an extended radio continuum halo at 20 cm and see for the first time the radio continuum feature associated with the western Halpha arm at wavelengths shorter than 20cm. The spectral index trends in this galaxy support the theory that there is a convective wind at work in this galaxy. We derive a total magnetic field strength of 38 microG in the central regions and 10-15 microG in the halo. The magnetic field is largely random in the center of the galaxy; the uniform field is ~3-9 microG and is strongest in the halo. We find that the magnetic pressure is the same order of magnitude but, in general, a factor of a few less than the other components of the interstellar medium in this galaxy. The uniform magnetic field in NGC 1569 is closely associated with the Halpha bubbles and filaments. We suggest that a supernova-driven dynamo may be operating in this galaxy. The outflow of hot gas from NGC 1569 is clearly shaping the magnetic field, but the magnetic field in turn may be aiding the outflow by channeling gas out of the disk of the galaxy. Dwarf galaxies with extended radio continuum halos like that of NGC 1569 may play an important role in magnetizing the intergalactic medium.

preprint2010arXiv

Violation of Richardson's Criterion via Introduction of a Magnetic Field

Shear flow instabilities can profoundly affect the diffusion of momentum in jets, stars, and disks. The Richardson criterion gives a sufficient condition for instability of a shear flow in a stratified medium. The velocity gradient $V'$ can only destabilize a stably stratified medium with squared Brunt-Vaisala frequency $N^2$ if $V'^2/4>N^2$. We find this is no longer true when the medium is a magnetized plasma. We investigate the effect of stable stratification on magnetic field and velocity profiles unstable to magneto-shear instabilities, i.e., instabilities which require the presence of both magnetic field and shear flow. We show that a family of profiles originally studied by Tatsuno & Dorland (2006) remain unstable even when $V'^2/4<N^2$, violating the Richardson criterion. However, not all magnetic fields can result in a violation of the Richardson criterion. We consider a class of flows originally considered by Kent (1968), which are destabilized by a constant magnetic field, and show that they become stable when $V'^2/4<N^2$, as predicted by the Richardson criterion. This suggests that magnetic free energy is required to violate the Richardson criterion. This work implies that the Richardson criterion cannot be used when evaluating the ideal stability of a sheared, stably stratified, and magnetized plasma. We briefly discuss the implications for astrophysical systems.

preprint2009arXiv

Environments for Magnetic Field Amplification by Cosmic Rays

We consider a recently discovered class of instabilities, driven by cosmic ray streaming, in a variety of environments. We show that although these instabilities have been discussed primarily in the context of supernova driven interstellar shocks, they can also operate in the intergalactic medium and in galaxies with weak magnetic fields, where, as a strong source of helical magnetic fluctuations, they could contribute to the overall evolution of the magnetic field. Within the Milky Way, these instabilities are strongest in warm ionized gas, and appear to be weak in hot, low density gas unless the injection efficiency of cosmic rays is very high.

preprint2009arXiv

MHD Stability of ISM Phase Transition Layers I: Magnetic Field Orthogonal to Front

We consider the scenario of a magnetic field orthogonal to a front separating two media of different temperatures and densities, such as cold and warm interstellar gas, in a 2-D plane-parallel geometry. A linear stability analysis is performed to assess the behavior of both evaporation and condensation fronts when subject to incompressible, corrugational perturbations with wavelengths larger than the thickness of the front. We discuss the behavior of fronts in both super-Alfvenic and sub-Alfvenic flows. Since the propagation speed of fronts is slow in the ISM, it is the sub-Alfvenic regime that is relevant, and magnetic fields are a significant influence on front dynamics. In this case we find that evaporation fronts, which are unstable in the hydrodynamic regime, are stabilized. Condensation fronts are unstable, but for parameters typical of the neutral ISM the growth rates are so slow that steady state fronts are effectively stable. However, the instability may become important if condensation proceeds at a sufficiently fast rate. This paper is the first in a series exploring the linear and nonlinear effects of magnetic field strength and orientation on the corrugational instability, with the ultimate goal of addressing outstanding questions about small-scale ISM structure.

preprint2007arXiv

The Milky Way's Kiloparsec Scale Wind: A Hybrid Cosmic-Ray and Thermally Driven Outflow

We apply a wind model, driven by combined cosmic-ray and thermal-gas pressure, to the Milky Way, and show that the observed Galactic diffuse soft X-ray emission can be better explained by a wind than by previous static gas models. We find that cosmic-ray pressure is essential to driving the observed wind. Having thus defined a "best-fit" model for a Galactic wind, we explore variations in the base parameters and show how the wind's properties vary with changes in gas pressure, cosmic-ray pressure and density. We demonstrate the importance of cosmic rays in launching winds, and the effect cosmic rays have on wind dynamics. In addition, this model adds support to the hypothesis of Breitschwerdt and collaborators that such a wind may help explain the relatively small gradient observed in gamma-ray emission as a function of galactocentric radius.

preprint2006arXiv

The Weak Field Limit of the Magnetorotational Instability

We investigate the behavior of the magneto-rotational instability in the limit of extremely weak magnetic field, i.e., as the ratio of ion cyclotron frequency to orbital frequency (X) becomes small. Considered only in terms of cold two-fluid theory, instability persists to arbitrarily small values of X, and the maximum growth rate is of order the orbital frequency except for the range m_e/m_i < |X| < 1, where it can be rather smaller. In this range, field aligned with rotation (X > 0) produces slower growth than anti-aligned field (X < 0). The maximum growth rate is generally achieved at smaller and smaller wavelengths as |X| diminishes. When |X| < m_e/m_i, new unstable "electromagnetic-rotational" modes appear that do not depend on the equilibrium magnetic field. Because the most rapidly-growing modes have extremely short wavelengths when |X| is small, they are often subject to viscous or resistive damping, which can result in suppressing all but the longest wavelengths, for which growth is much slower. We find that this sort of damping is likely to curtail severely the frequently-invoked mechanism for cosmological magnetic field growth in which a magnetic field seeded by the Biermann battery is then amplified by the magneto-rotational instability. On the other hand, the small |X| case may introduce interesting effects in weakly-ionized disks in which dust grains carry most of the electric charge.