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.

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\sigma$ 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.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.

The Parker instability, a Rayleigh-Taylor like instability of thermal gas supported against gravity by magnetic fields and cosmic rays, is thought to be dynamically important for galaxy evolution, possibly promoting molecular cloud formation and the galactic dynamo. In previous work, we examined the effect of three different cosmic ray transport models on the Parker instability: decoupled ($\gamma_c = 0$), locked to the thermal gas ($\gamma_c = 4/3$) and coupled to the gas with streaming by self-confinement. We expand upon that work here by considering radiative cooling, a smooth gravitational potential, and simulations into the nonlinear regime. We determine that cosmic ray transport away from compression points, whether by diffusion or streaming, is the largest driver of the instability. Heating due to cosmic ray streaming is also destabilizing and especially affects the nonlinear regime. While cooling de-pressurizes the dense gas, streaming cosmic rays heat and inflate the diffuse extraplanar gas, greatly modifying the phase structure of the medium. In 3D, we find that the fastest growth favors short wavelength modes in the horizontal direction perpendicular to the background magnetic field; this is imprinted on Faraday rotation measure maps that may be used to detect the Parker instability. The modifications to the Parker instability that we observe in this work have large implications for the structure and evolution of galaxies, and they highlight the major role that cosmic rays play in shaping their environments.

We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $\mu$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $\mu$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um.

For the first time, we characterize the rigidity regimes of the diffusion coefficient $\kappa$ for arbitrary rigidities and guide fields, which we derive as a function of physical and numerical parameters. We show that at turbulence levels above 5\% of the total magnetic field, the approximation of an energy dependence $\kappa\propto E^{1/3}$ as predicted for a Kolmogorov spectrum within Quasi-Linear Theory does not hold. Consequently, a proper description of cosmic-ray propagation can only be achieved by using a turbulence-level dependent diffusion coefficient and can contribute to solving the Galactic cosmic-ray gradient problem.

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.

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.

Using hybrid simulations (kinetic ions--fluid electrons), we test the linear theory predictions of the cosmic ray (CR) streaming instability. We consider two types of CR distribution functions: a "hot" distribution where CRs are represented by a drifting power law in momentum and an anisotropic "beam" of monochromatic particles. Additionally, for each CR distribution we scan over different CR densities to transition from triggering the resonant to the non-resonant (Bell) streaming instability. We determine the growth rates of these instabilities in simulations by fitting an exponential curve during the linear stage, and we show that they agree well with the theoretical predictions as a function of wave number agree. We also examine the magnetic helicity as a function of time and wave number, finding a general good agreement with the predictions, as well as some unexpected non-linear features to the instability development.

This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

The Tayler instability (TI) is a non-axisymmetric linear instability of an axisymmetric toroidal magnetic field in magneto-hydrostatic equilibrium (MHSE). Spruit (1999, 2002) has proposed that in a differentially rotating radiative region of a star, the TI drives a dynamo which generates magnetic fields that can efficiently transport angular momentum; a parameterized version of this dynamo has been implemented in stellar structure and evolution codes and shown to be important for determining interior spin. Numerical simulations, however, have yet to definitively demonstrate the operation of the dynamo. A criterion for the MHSE to develop the TI was derived using fully-compressible magneto-hydrodynamics, while numerical simulations of dynamical processes in stars frequently use an anelastic approximation. This motivates us to derive a new anelastic Tayler instability (anTI) criterion. We find that some MHSE configurations unstable in the fully-compressible case, become stable in the anelastic case. We find and characterize the unstable modes of a simple family of cylindrical MHSE configurations using numerical calculations, and discuss the implications for fully non-linear anelastic simulations.

The nature of cosmic rays (CRs) and cosmic ray transport in galaxy clusters is probed by a number of observations. Radio observations reveal the synchrotron radiation of cosmic ray electrons (CRe) spiraling around cluster magnetic fields. $\gamma$-ray observations reveal hadronic reactions of cosmic ray protons (CRp) with ambient gas nuclei which produce pions. To date, no such cluster-wide $\gamma$-ray signal has been measured, putting an upper limit on the density of CRp present in clusters. But the presence of CRe implies some source of CRp, and consequently there must be some CRp loss mechanism. In this paper we quantify the observational constraints on this loss mechanism assuming that losses are dominated by CR transport, ultimately deriving a minimum diffusion coefficient of $\sim10^{31}$ cm$^2$ s$^{-1}$ in the Coma cluster. This lower limit on transport may help illuminate some unknown properties of the cluster field topology. Conversely, measurements of cluster field tangling scales can constrain other model parameters, such as the relative acceleration efficiency of protons to electrons. To be consistent with the Coma observations, protons cannot be accelerated more than 15 times more efficiently than electrons of the same energy.

Understanding the processes that drive galaxy formation and shape the observed properties of galaxies is one of the most interesting and challenging frontier problems of modern astrophysics. We now know that the evolution of galaxies is critically shaped by the energy injection from accreting supermassive black holes (SMBHs). However, it is unclear how exactly the physics of this feedback process affects galaxy formation and evolution. In particular, a major challenge is unraveling how the energy released near the SMBHs is distributed over nine orders of magnitude in distance throughout galaxies and their immediate environments. The best place to study the impact of SMBH feedback is in the hot atmospheres of massive galaxies, groups, and galaxy clusters, which host the most massive black holes in the Universe, and where we can directly image the impact of black holes on their surroundings. We identify critical questions and potential measurements that will likely transform our understanding of the physics of SMBH feedback and how it shapes galaxies, through detailed measurements of (i) the thermodynamic and velocity fluctuations in the intracluster medium (ICM) as well as (ii) the composition of the bubbles inflated by SMBHs in the centers of galaxy clusters, and their influence on the cluster gas and galaxy growth, using the next generation of high spectral and spatial resolution X-ray and microwave telescopes.

Galaxy clusters are the largest and most massive bound objects resulting from cosmic hierarchical structure formation. Baryons account for somewhat more than 10% of that mass, with roughly 90% of the baryonic matter distributed throughout the clusters as hot ($T>1$ keV), high-$\beta$, very weakly collisional plasma; the so-called "intracluster medium" (ICM). Cluster mergers, close gravitational encounters and accretion, along with violent feedback from galaxies and relativistic jets from active galactic nuclei, drive winds, gravity waves, turbulence and shocks within the ICM. Those dynamics, in turn, generate cluster-scale magnetic fields and accelerate and mediate the transport of high-energy charged particles. Kinetic-scale, collective plasma processes define the basic character and fundamental signatures of these ICM phenomena, which are observed primarily by X-ray and radio astronomers.

This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood: we lack a predictive theory for how it responds to perturbations, how it transports momentum and energy, and how it generates and amplifies magnetic fields. Such topics are foundational to the scientific study of plasmas, and are of intrinsic interest to those who regard plasma physics as a fundamental physics discipline. But these topics are also of extrinsic interest: addressing them directly informs upon our understanding of a wide variety of space and astrophysical systems, including accretion flows around supermassive black holes, the intracluster medium (ICM) between galaxies in clusters, and regions of the near-Earth solar wind. Specific recommendations to advance this field of study are discussed.

We investigate a mechanism for accelerating cool (10$^4$ K) clouds in the circumgalactic medium (CGM) with cosmic rays (CRs), possibly explaining some characteristics of observed high velocity clouds (HVCs). Enforcing CRs to stream down their pressure gradient into a region of slow streaming speed results in significant buildup of CR pressure which can accelerate the CGM. We present the results of the first two-dimensional magnetohydrodynamic (MHD) simulations of such `CR bottlenecks,' expanding on simpler simulations in 1D from \cite{wiener17a}. Although much more investigation is required, we find two main results. First, radiative cooling in the interfaces of these clouds is sufficient to keep the cloud intact to CR wave heating. Second, cloud acceleration depends almost linearly with the injected CR flux at low values (comparable to that expected from a Milky Way-like star formation rate), but scales sublinearly at higher CR fluxes in 1D simulations. 2D simulations show hints of sublinear dependence at high CR fluxes but are consistent with pure linear dependence up to the CR fluxes tested. It is therefore plausible to accelerate cool clouds in the CGM to speeds of hundreds of km s$^{-1}$.

We present Far-Infrared polarimetry observations of M82 at 53 and $154~\mu \rm{m}$ and NGC 253 at $89~\mu \rm{m}$, which were taken with HAWC+ in polarimetry mode on the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization of M82 at $53~\mu \rm{m}$ clearly shows a magnetic field geometry perpendicular to the disk in the hot dust emission. For M82 the polarization at $154~\mu \rm{m}$ shows a combination of field geometry perpendicular to the disk in the nuclear region, but closer to parallel to the disk away from the nucleus. The fractional polarization at $53~\mu \rm{m}$ $(154~\mu \rm{m})$ ranges from 7% (3%) off nucleus to 0.5% (0.3%) near the nucleus. A simple interpretation of the observations of M82 invokes a massive polar outflow, dragging the field along, from a region $\sim 700$~pc in diameter that has entrained some of the gas and dust, creating a vertical field geometry seen mostly in the hotter $(53~\mu \rm{m})$ dust emission. This outflow sits within a larger disk with a more typical planar geometry that more strongly contributes to the cooler $(154~\mu \rm{m})$ dust emission. For NGC 253, the polarization at $89~\mu \rm{m}$ is dominated by a planar geometry in the tilted disk, with weak indication of a vertical geometry above and below the plane from the nucleus. The polarization observations of NGC 253 at $53~\mu \rm{m}$ were of insufficient S/N for detailed analysis.

It has been known for half a century that the interstellar medium (ISM) of our Galaxy is structured on scales as small as a few hundred km, more than 10 orders of magnitude smaller than typical ISM structures and energy input scales. In this review we focus on neutral and ionized structures on spatial scales of a few to ~10^4 Astronomical Units (AU) which appear to be highly overpressured, as these have the most important role in the dynamics and energy balance of interstellar gas: the Tiny Scale Atomic Structure (TSAS) and Extreme Scattering Events (ESEs) as the most over-pressured example of the Tiny Scale Ionized Structures (TSIS). We review observational results and highlight key physical processes at AU scales. We present evidence for and against microstructures as part of a universal turbulent cascade and as discrete structures, and review their association with supernova remnants, the Local Bubble, and bright stars. We suggest a number of observational and theoretical programs that could clarify the nature of AU structures. TSAS and TSIS probe spatial scales in the range of what is expected for turbulent dissipation scales, therefore are of key importance for constraining exotic and not-well understood physical processes which have implications for many areas of astrophysics. The emerging picture is one in which a magnetized, turbulent cascade, driven hard by a local energy source and acting jointly with phenomena such as thermal instability, is the source of these microstructures.

A linearly unstable, sinusoidal $E \times B$ shear flow is examined in the gyrokinetic framework in both the linear and nonlinear regimes. In the linear regime, it is shown that the eigenmode spectrum is nearly identical to hydrodynamic shear flows, with a conjugate stable mode found at every unstable wavenumber. In the nonlinear regime, turbulent saturation of the instability is examined with and without the inclusion of a driving term that prevents nonlinear flattening of the mean flow, and a scale-independent radiative damping term that suppresses the excitation of conjugate stable modes. A simple fluid model for how momentum transport and partial flattening of the mean flow scale with the driving term is constructed, from which it is shown that, except at high radiative damping, stable modes play an important role in the turbulent state and yield significantly improved quantitative predictions when compared with corresponding models neglecting stable modes.

Galactic winds and fountains driven by supernova-heated gas play an integral role in re-distributing gas in galaxies, depositing metals in the circumgalactic medium (CGM), and quenching star formation. The interplay between these outflows and ram pressure stripping due to the galaxy's motion through an ambient medium may enhance these effects by converting fountain flows into expelled gas. In this paper, we present controlled, 3D simulations of ram pressure stripping combined with thermally driven, local outflows from clustered supernovae in an isolated disk galaxy modeled on the Large Magellanic Cloud (LMC), a dwarf satellite of the Milky Way on its first infall. Observational evidence of local outflows emanating from supergiant shells in the LMC and a trailing filament of HI gas originating from these regions - with no obvious Leading Arm counterpart - may represent a perfect example of this process. Our simulations present a proof-of-concept that ram pressure can convert fountain flows into expelled gas. We find that fountains launched near the peak star formation time of the LMC can comprise part of the LMC filament in the Trailing Stream, but with lower column densities than observed. Larger, more numerous outflows from the LMC may be possible and may contribute more mass, but higher inertia gas will lengthen the timescale for this gas to be swept away by ram pressure. Given the high resolution observations, increased knowledge of star formation histories, and growing evidence of multiphase, ionized outflows, the LMC is an ideal test-bed for future wind models.

Recent studies have found that cosmic ray transport plays an important role in feedback processes such as star formation and the launching of galactic winds. Although cosmic ray buoyancy is widely held to be a destabilizing force in galactic disks, the effect of cosmic ray transport on the stability of stratified systems has yet to be analyzed. We perform a stability analysis of a stratified layer for three different cosmic ray transport models: decoupled (Classic Parker), coupled with $\gamma_c=4/3$ but not streaming (Modified Parker), and finally coupled with streaming at the Alfv\'{e}n speed. When the compressibility of the cosmic rays is decreased the system becomes much more stable, but the addition of cosmic ray streaming to the Parker Instability severely destabilizes it. Through comparison of these three cases and analysis of the work contributions for the perturbed quantities of each system, we demonstrate that cosmic ray heating of the gas is responsible for the destabilization of the system. We find that a 3D system is unstable over a larger range of wavelengths than the 2D system. Therefore, the Parker Instability with cosmic ray streaming may play an important role in cosmic ray feedback.

Although several theories for the origin of cosmic rays in the region between the spectral `knee' and `ankle' exist, this problem is still unsolved. A variety of observations suggest that the transition from Galactic to extragalactic sources occurs in this energy range. In this work we examine whether a Galactic wind which eventually forms a termination shock far outside the Galactic plane can contribute as a possible source to the observed flux in the region of interest. Previous work by Bustard et al. (2017) estimated that particles can be accelerated up to energies above the `knee' up to $R_\mathrm{max} = 10^{16}$ eV for parameters drawn from a model of a Milky Way wind (Everett et al. 2017). A remaining question is whether the accelerated cosmic rays can propagate back into the Galaxy. To answer this crucial question, we simulate the propagation of the cosmic rays using the low energy extension of the CRPropa framework, based on the solution of the transport equation via stochastic differential equations. The setup includes all relevant processes, including three-dimensional anisotropic spatial diffusion, advection, and corresponding adiabatic cooling. We find that, assuming realistic parameters for the shock evolution, a possible Galactic termination shock can contribute significantly to the energy budget in the `knee' region and above. We estimate the resulting produced neutrino fluxes and find them to be below measurements from IceCube and limits by KM3NeT.

Galaxy cluster cores are pervaded by hot gas which radiates at far too high a rate to maintain any semblance of a steady state; this is referred to as the cooling flow problem. Of the many heating mechanisms that have been proposed to balance radiative cooling, one of the most attractive is dissipation of acoustic waves generated by Active Galactic Nuclei (AGN). Fabian (2005) showed that if the waves are nearly adiabatic, wave damping due to heat conduction and viscosity must be well below standard Coulomb rates in order to allow the waves to propagate throughout the core. Because of the importance of this result, we have revisited wave dissipation under galaxy cluster conditions in a way that accounts for the self limiting nature of dissipation by electron thermal conduction, allows the electron and ion temperature perturbations in the waves to evolve separately, and estimates kinetic effects by comparing to a semi-collisionless theory. While these effects considerably enlarge the toolkit for analyzing observations of wavelike structures and developing a quantitative theory for wave heating, the drastic reduction of transport coefficients proposed in Fabian (2005) remains the most viable path to acoustic wave heating of galaxy cluster cores.

The Fermi bubbles, two giant structures above and below the Galactic center (GC), are among the most important discoveries of the Fermi Gamma-ray Space Telescope. Studying their physical origin has been providing valuable insights into cosmic-ray transport, the Galactic magnetic field, and past activity at the GC in the Milky Way galaxy. Despite their importance, the formation mechanism of the bubbles is still elusive. Over the past few years there have been numerous efforts, both observational and theoretical, to uncover the nature of the bubbles. In this article, we present an overview of the current status of our understanding of the bubbles' origin, and discuss possible future directions that will help to distinguish different scenarios of bubble formation.

The role of cosmic rays generated by supernovae and young stars has very recently begun to receive significant attention in studies of galaxy formation and evolution due to the realization that cosmic rays can efficiently accelerate galactic winds. Microscopic cosmic ray transport processes are fundamental for determining the efficiency of cosmic ray wind driving. Previous studies focused on modeling of cosmic ray transport either via constant diffusion coefficient or via streaming proportional to the Alfven speed. However, in predominantly cold, neutral gas, cosmic rays can propagate faster than in the ionized medium and the effective transport can be substantially larger; i.e., cosmic rays can decouple from the gas. We perform three-dimensional magnetohydrodynamical simulations of patches of galactic disks including the effects of cosmic rays. Our simulations include the decoupling of cosmic rays in the cold, neutral interstellar medium. We find that, compared to the ordinary diffusive cosmic ray transport case, accounting for the decoupling leads to significantly different wind properties such as the gas density and temperature, significantly broader spatial distribution of cosmic rays, and larger wind speed. These results have implications for X-ray, gamma-ray and radio emission, and for the magnetization and pollution of the circumgalactic medium by cosmic rays.

Motivated by recent interest in the phenomenon of waves transport in massive stars, we examine whether the heat-driven gravity (g) modes excited in slowly-pulsating B (SPB) stars can significantly modify the stars' internal rotation. We develop a formalism for the differential torque exerted by g modes, and implement this formalism using the GYRE oscillation code and the MESASTAR stellar evolution code. Focusing first on a $4.21$ $M_\odot$ model, we simulate 1,000 years of stellar evolution under the combined effects of the torque due to a single unstable prograde g mode (with an amplitude chosen on the basis of observational constraints), and diffusive angular momentum transport due to convection, overshooting, and rotational instabilities. We find that the g mode rapidly extracts angular momentum from the surface layers, depositing it deeper in the stellar interior. The angular momentum transport is so efficient that by the end of the simulation the initially non-rotating surface layers are spun in the retrograde direction to $\approx30\%$ of the critical rate. However, the additional inclusion of magnetic stresses in our simulations, almost completely inhibits this spin-up. Expanding our simulations to cover the whole instability strip, we show that the same general behavior is seen in all SPB stars. After providing some caveats to contextualize our results, we hypothesize that the observed slower surface rotation of SPB stars (as compared to other B-type stars) may be the direct consequence of the angular momentum transport that our simulations demonstrate.

Diffuse, extended radio emission in galaxy clusters, commonly referred to as radio halos, indicate the presence of high energy cosmic ray (CR) electrons and cluster-wide magnetic fields. We can predict from theory the expected surface brightness of a radio halo, given magnetic field and CR density profiles. Previous studies have shown that the nature of CR transport can radically effect the expected radio halo emission from clusters (Wiener et al. 2013). Reasonable levels of magnetohydrodynamic (MHD) wave damping can lead to significant CR streaming speeds. But a careful treatment of MHD waves in a high $\beta$ plasma, as expected in cluster environments, reveals damping rates may be enhanced by a factor of $\beta^{1/2}$. This leads to faster CR streaming and lower surface brightnesses than without this effect. In this work we re-examine the simplified, 1D Coma cluster simulations (with radial magnetic fields) of Wiener et al. (2013) and discuss observable consequences of this high $\beta$ damping. Future work is required to study this effect in more realistic simulations.

Magnetic fields play a pivotal role in the physics of interstellar medium in galaxies, but there are few observational constraints on how they evolve across cosmic time. Spatially resolved synchrotron polarization maps at radio wavelengths reveal well-ordered large-scale magnetic fields in nearby galaxies that are believed to grow from a seed field via a dynamo effect. To directly test and characterize this theory requires magnetic field strength and geometry measurements in cosmologically distant galaxies, which are challenging to obtain due to the limited sensitivity and angular resolution of current radio telescopes. Here, we report the cleanest measurements yet of magnetic fields in a galaxy beyond the local volume, free of the systematics traditional techniques would encounter. By exploiting the scenario where the polarized radio emission from a background source is gravitationally lensed by a foreground galaxy at z = 0.439 using broadband radio polarization data, we detected coherent $\mu$G magnetic fields in the lensing disk galaxy as seen 4.6 Gyrs ago, with similar strength and geometry to local volume galaxies. This is the highest redshift galaxy whose observed coherent magnetic field property is compatible with a mean-field dynamo origin.

The Cygnus X region contains giant molecular cloud complexes and populous associates of massive young stars. The discovery of spatially extended, hard gamma-ray emission in Cygnus X by both Milagro and Fermi indicates that Cygnus X is also a potential source of high-energy Galactic neutrinos. Here, we adapt our single-zone model for cosmic ray interactions in the central molecular zones of starburst galaxies for use in Cygnus X. We calculate the potential neutrino flux corresponding to the hard gamma-ray emission from the "Cygnus Cocoon" and to the soft, diffuse interstellar gamma-ray emission. We check our results by comparing the corresponding gamma-ray emission against the Fermi interstellar emission model and Milagro, ARGO-YBJ, and HAWC observations. In comparing our results against a recent IceCube analysis and the current sensitivity limits, we find that neutrino emission from the Cocoon has a large enough flux that it could plausibly be detected, provided hadronic interactions are occurring at sufficiently high energies. High-energy neutrinos from Cygnus X would provide direct evidence for the presence of as yet unidentified PeV energy accelerators in the Galactic disk.

We present the first kinematic study of extraplanar diffuse ionized gas (eDIG) in the nearby, face-on disk galaxy M83 using optical emission-line spectroscopy from the Robert Stobie Spectrograph on the Southern African Large Telescope. We use a Markov Chain Monte Carlo method to decompose the [NII]$\lambda\lambda$6548, 6583, H$\alpha$, and [SII]$\lambda\lambda$6717, 6731 emission lines into HII region and diffuse ionized gas emission. Extraplanar, diffuse gas is distinguished by its emission-line ratios ([NII]$\lambda$6583/H$\alpha \gtrsim 1.0$) and its rotational velocity lag with respect to the disk ($\Delta v = -24$ km/s in projection). With interesting implications for isotropy, the velocity dispersion of the diffuse gas, $\sigma = 96$ km/s, is a factor of a few higher in M83 than in the Milky Way and nearby, edge-on disk galaxies. The turbulent pressure gradient is sufficient to support the eDIG layer in dynamical equilibrium at an electron scale height of $h_{z} = 1$ kpc. However, this dynamical equilibrium model must be finely tuned to reproduce the rotational velocity lag. There is evidence of local bulk flows near star-forming regions in the disk, suggesting that the dynamical state of the gas may be intermediate between a dynamical equilibrium and a galactic fountain flow. As one of the first efforts to study eDIG kinematics in a face-on galaxy, this study demonstrates the feasibility of characterizing the radial distribution, bulk velocities, and vertical velocity dispersions in low-inclination systems.

We investigate the effects of cosmic ray (CR) dynamics on cold, dense clouds embedded in a hot, tenuous galactic halo. If the magnetic field does not increase too much inside the cloud, the local reduction in Alfv\'en speed imposes a bottleneck on CRs streaming out from the star-forming galactic disk. The bottleneck flattens the upstream CR gradient in the hot gas, implying that multi-phase structure could have global effects on CR driven winds. A large CR pressure gradient can also develop on the outward-facing edge of the cloud. This pressure gradient has two independent effects. The CRs push the cloud upward, imparting it with momentum. On smaller scales, the CRs pressurize cold gas in the fronts, reducing its density, consistent with the low densities of cold gas inferred in recent COS observations of local $L_{*}$ galaxies. They also heat the material at the cloud edge, broadening the cloud-halo interface and causing an observable change in interface ionic abundances. Due to the much weaker temperature dependence of cosmic ray heating relative to thermal conductive heating, CR mediated fronts have a higher ratio of low to high ions compared to conduction fronts, in better agreement with observations. We investigate these effects separately using 1D simulations and analytic techniques.

There are two distinct breaks in the cosmic ray (CR) spectrum: the so-called "knee" around $3 \times 10^{15}$ eV and the so-called "ankle" around $10^{18}$ eV. Diffusive shock acceleration (DSA) at supernova remnant (SNR) shock fronts is thought to accelerate galactic CRs to energies below the knee, while an extragalactic origin is presumed for CRs with energies beyond the ankle. CRs with energies between $3 \times 10^{15}$ and $10^{18}$ eV, which we dub the "shin," have an unknown origin. It has been proposed that DSA at galactic wind termination shocks, rather than at SNR shocks, may accelerate CRs to these energies. This paper uses the galactic wind model of Bustard et al. (2016) to analyze whether galactic wind termination shocks may accelerate CRs to shin energies within a reasonable acceleration time and whether such CRs can subsequently diffuse back to the galaxy. We argue for acceleration times on the order of 100 Myrs rather than a few billion years, as assumed in some previous works, and we discuss prospects for magnetic field amplification at the shock front. Ultimately, we generously assume that the magnetic field is amplified to equipartition. This formalism allows us to obtain analytic formulae, applicable to any wind model, for CR acceleration. Even with generous assumptions, we find that very high wind velocities are required to set up the necessary conditions for acceleration beyond $10^{17}$ eV. We also estimate the luminosities of CRs accelerated by outflow termination shocks, including estimates for the Milky Way wind.

The observed scale heights of extraplanar diffuse ionized gas (eDIG) layers exceed their thermal scale heights by a factor of a few in the Milky Way and other nearby edge-on disk galaxies. Here, we test a dynamical equilibrium model of the extraplanar diffuse ionized gas layer in NGC 891, where we ask whether the thermal, turbulent, magnetic field, and cosmic ray pressure gradients are sufficient to support the layer. In optical emission line spectroscopy from the SparsePak integral field unit on the WIYN 3.5-meter telescope, the H-alpha emission in position-velocity space suggests that the eDIG is found in a ring between galactocentric radii of R_min <= R <= 8 kpc, where R_min >= 2 kpc. We find that the thermal (sigma_th = 11 km/s) and turbulent (sigma_turb = 25 km/s) velocity dispersions are insufficient to satisfy the hydrostatic equilibrium equation given an exponential electron scale height of h_z = 1.0 kpc. Using a literature analysis of radio continuum observations from the CHANG-ES survey, we demonstrate that the magnetic field and cosmic ray pressure gradients are sufficient to stably support the gas at R >= 8 kpc if the cosmic rays are sufficiently coupled to the system (gamma_cr = 1.45). Thus, a stable dynamical equilibrium model is viable only if the extraplanar diffuse ionized gas is found in a thin ring around R = 8 kpc, and non-equilibrium models such as a galactic fountain flow are of interest for further study.

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~\mu$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~\mu$G in strength, tidal interactions can produce large-scale regular field structure in the outskirts.

Narrow-band HST imaging has resolved the detailed internal structure of the 10 kpc diameter H alpha+[NII] emission line nebulosity in NGC4696, the central galaxy in the nearby Centaurus cluster, showing that the dusty, molecular, filaments have a width of about 60pc. Optical morphology and velocity measurements indicate that the filaments are dragged out by the bubbling action of the radio source as part of the AGN feedback cycle. Using the drag force we find that the magnetic field in the filaments is in approximate pressure equipartition with the hot gas. The filamentary nature of the cold gas continues inward, swirling around and within the Bondi accretion radius of the central black hole, revealing the magnetic nature of the gas flows in massive elliptical galaxies. HST imaging resolves the magnetic, dusty, molecular filaments at the centre of the Centaurus cluster to a swirl around and within the Bondi radius.

Galactic outflows play an important role in galactic evolution. Despite their importance, a detailed understanding of the physical mechanisms responsible for the driving of these winds is lacking. In an effort to gain more insight into the nature of these flows, we perform global three-dimensional magneto-hydrodynamical simulations of an isolated Milky Way-size starburst galaxy. We focus on the dynamical role of cosmic rays injected by supernovae, and specifically on the impact of the streaming and anisotropic diffusion of cosmic rays along the magnetic fields. We find that these microphysical effects can have a significant effect on the wind launching and mass loading factors depending on the details of the plasma physics. Due to the cosmic ray streaming instability, cosmic rays propagating in the interstellar medium scatter on self-excited Alfven waves and couple to the gas. When the wave growth due to the streaming instability is inhibited by some damping process, such as the turbulent damping, the cosmic ray coupling to the gas is weaker and their effective propagation speed faster than the Alfven speed. Alternatively, cosmic rays could scatter from "extrinsic turbulence" that is driven by another mechanism. We demonstrate that the presence of moderately super-Alfvenic cosmic ray streaming enhances the efficiency of galactic wind driving. Cosmic rays stream away from denser regions near the galactic disk along partially ordered magnetic fields and, in the process, accelerate more tenuous gas away from the galaxy. For cosmic ray acceleration efficiencies broadly consistent with the observational constraints, cosmic rays reduce the galactic star formation rates and significantly aid in launching galactic winds.

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.

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.

The Wisconsin Plasma Astrophysics Laboratory (WiPAL) is a flexible user facility designed to study a range of astrophysically relevant plasma processes as well as novel geometries that mimic astrophysical systems. A multi-cusp magnetic bucket constructed from strong samarium cobalt permanent magnets now confines a 10 m$^3$, fully ionized, magnetic-field free plasma in a spherical geometry. Plasma parameters of $ T_{e}\approx5$ to $20$ eV and $n_{e}\approx10^{11}$ to $5\times10^{12}$ cm$^{-3}$ provide an ideal testbed for a range of astrophysical experiments including self-exciting dynamos, collisionless magnetic reconnection, jet stability, stellar winds, and more. This article describes the capabilities of WiPAL along with several experiments, in both operating and planning stages, that illustrate the range of possibilities for future users.

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.

We present a study of the magneto-ionic medium in the Whirlpool galaxy (M51) using new wide-band multi-configuration polarization data at L band (1-2 GHz) obtained at the Karl G. Jansky Very Large Array. By fitting the observed diffuse complex polarization $Q$+$iU$ as a function of wavelength directly to various depolarization models, we find that polarized emission from M51 at 1-2 GHz originates from the top of the synchrotron disk and then experiences Faraday rotation in the near-side thermal halo of the galaxy. Thus, the scale height of the thermal gas must exceed that of the synchrotron emitting gas at L band. The observed Faraday depth distribution at L band is consistent with a halo field that comprises of a plane-parallel bisymmetric component and a vertical component which produces a Faraday rotation of $\sim$ $-$9 rad m$^{-2}$. The derived rotation measure structure functions indicate a characteristic scale of rotation measure fluctuations of less than 560 pc in the disk and approximately 1 kpc in the halo. The outer scale of turbulence of 1 kpc found in the halo of M51 is consistent with superbubbles and the Parker instability being the main energy injection mechanisms in galactic halos.

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.

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.

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.

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.

The $Fermi$ gamma-ray telescope discovered a pair of bubbles at the Galactic center. These structures are spatially-correlated with the microwave emission detected by the WMAP and Planck satellites. These bubbles were likely inflated by a jet launched from the vicinity of a supermassive black hole in the Galactic center. Using MHD simulations, which self-consistently include interactions between cosmic rays and magnetic fields, we build models of the supersonic jet propagation, cosmic ray transport, and the magnetic field amplification within the $Fermi$ bubbles. Our key findings are that: (1) the synthetic $Fermi$ gamma-ray and WMAP microwave spectra based on our simulations are consistent with the observations, suggesting that a single population of cosmic ray leptons may simultaneously explain the emission across a range of photon energies; (2) the model fits the observed centrally-peaked microwave emission if a second, more recent, pair of jets embedded in the $Fermi$ bubbles is included in the model. This is consistent with the observationally-based suggestion made by Su & Finkbeiner (2012); (3) the radio emission from the bubbles is expected to be strongly polarized due to the relatively high level of field ordering caused by elongated turbulent vortices. This effect is caused by the interaction of the shocks driven by the jets with the preexisting interstellar medium turbulence; (4) a layer of enhanced rotation measure in the shock-compressed region could exist in the bubble vicinity but the level of this enhancement depends on the details of the magnetic topology.

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.

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.

The origin of the Fermi bubbles and the microwave haze is yet to be determined. To disentangle different models requires detailed comparisons between theoretical predictions and multi-wavelength observations. Our previous simulations have demonstrated that the primary features of the Fermi bubbles could be successfully reproduced by recent jet activity from the central active galactic nucleus (AGN). In this work, we generate gamma-ray and microwave maps and spectra based on the simulated properties of cosmic rays (CRs) and magnetic fields in order to examine whether the observed bubble and haze emission could be explained by leptons contained in the AGN jets. We also investigate the model predictions of the polarization properties of the Fermi bubbles. We find that: (1) The same population of leptons can simultaneously explain the bubble and haze emission given that the magnetic fields within the bubbles are very close to the exponentially distributed ambient field, which can be explained by mixing in of the ambient field followed by turbulent field amplification; (2) The centrally peaked microwave profile suggests CR replenishment, which is consistent with the presence of a more recent second jet event; (3) The bubble interior exhibits a high degree of polarization because of ordered radial magnetic field lines stretched by elongated vortices behind the shocks; highly-polarized signals could also be observed inside the draping layer; (4) Enhancement of rotation measures could exist within the shock-compressed layer because of increased gas density and more amplified and ordered magnetic fields. We discuss the possibility that the deficient haze emission at b<-35 degrees is due to the suppression of magnetic fields, which is consistent with the existence of lower-energy CRs causing the polarized emission at 2.3 GHz. Possible AGN jet composition in the leptonic scenario is also discussed.

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 \gamma-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 \gamma-ray spectra of M82. Chi-squared tests are used with radio and \gamma-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 \mu 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.

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.

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, $\Lambda 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\'en 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.

The Fermi Gamma-ray Space Telescope reveals two large bubbles in the Galaxy, which extend nearly symmetrically ~50 degrees above and below the Galactic center (GC). Using three-dimensional (3D) magnetohydrodynamic (MHD) simulations that self-consistently include the dynamical interaction between cosmic rays (CR) and thermal gas, and anisotropic CR diffusion along the magnetic field lines, we show that the key characteristics of the observed gamma-ray bubbles and the spatially-correlated X-ray features in ROSAT 1.5 keV map can be successfully reproduced by a recent jet activity from the central active galactic nucleus (AGN). We find that after taking into account the projection of the 3D bubbles onto the sky, the physical heights of the bubbles can be much smaller than previously thought, greatly reducing the formation time of the bubbles to about a Myr. This relatively small bubble age is needed to reconcile the simulations with the upper limit of bubbles ages estimated from the cooling time of high-energy electrons. No additional physical mechanisms are required to suppress large-scale hydrodynamic instabilities because the evolution time is too short for them to develop. The simulated CR bubbles are edge-brightened, which is consistent with the observed projected flat surface brightness distribution. Furthermore, we demonstrate that the sharp edges of the observed bubbles can be due to anisotropic CR diffusion along magnetic field lines that drape around the bubbles during their supersonic expansion, with suppressed perpendicular diffusion across the bubble surface. Possible causes of the slight bends of the Fermi bubbles to the west are also discussed.

Typical flows in stellar interiors are much slower than the speed of sound. To follow the slow evolution of subsonic motions, various sound-proof equations are in wide use, particularly in stellar astrophysical fluid dynamics. These low-Mach number equations include the anelastic equations. Generally, these equations are valid in nearly adiabatically stratified regions like stellar convection zones, but may not be valid in the sub-adiabatic, stably stratified stellar radiative interiors. Understanding the coupling between the convection zone and the radiative interior is a problem of crucial interest and may have strong implications for solar and stellar dynamo theories as the interface between the two, called the tachocline in the Sun, plays a crucial role in many solar dynamo theories. Here we study the properties of gravity waves in stably-stratified atmospheres. In particular, we explore how gravity waves are handled in various sound-proof equations. We find that some anelastic treatments fail to conserve energy in stably-stratified atmospheres, instead conserving pseudo-energies that depend on the stratification, and we demonstrate this numerically. One anelastic equation set does conserve energy in all atmospheres and we provide recommendations for converting low-Mach number anelastic codes to this set of equations.

We explore the suggestions by Uzdensky (2007) and Cassak et al. (2008) that coronal loops heated by magnetic reconnection should self-organize to a state of marginal collisionality. We discuss their model of coronal loop dynamics with a one-dimensional hydrodynamic calculation. We assume that many current sheets are present, with a distribution of thicknesses, but that only current sheets thinner than the ion skin depth can rapidly reconnect. This assumption naturally causes a density dependent heating rate which is actively regulated by the plasma. We report 9 numerical simulation results of coronal loop hydrodynamics in which the absolute values of the heating rates are different but their density dependences are the same. We find two regimes of behavior, depending on the amplitude of the heating rate. In the case that the amplitude of heating is below a threshold value, the loop is in stable equilibrium. Typically the upper and less dense part of coronal loop is collisionlessly heated and conductively cooled. When the amplitude of heating is above the threshold, the conductive flux to the lower atmosphere required to balance collisionless heating drives an evaporative flow which quenches fast reconnection, ultimately cooling and draining the loop until the cycle begins again. The key elements of this cycle are gravity and the density dependence of the heating function. Some additional factors are present, including pressure driven flows from the loop top, which carry a large enthalpy flux and play an important role in reducing the density. We find that on average the density of the system is close to the marginally collisionless value.

An emerging theme in modern astrophysics is the connection between astronomical observations and the underlying physical phenomena that drive our cosmos. Both the mechanisms responsible for the observed astrophysical phenomena and the tools used to probe such phenomena - the radiation and particle spectra we observe - have their roots in atomic, molecular, condensed matter, plasma, nuclear and particle physics. Chemistry is implicitly included in both molecular and condensed matter physics. This connection is the theme of the present report, which provides a broad, though non-exhaustive, overview of progress in our understanding of the cosmos resulting from recent theoretical and experimental advances in what is commonly called laboratory astrophysics. This work, carried out by a diverse community of laboratory astrophysicists, is increasingly important as astrophysics transitions into an era of precise measurement and high fidelity modeling.

We have previously found a temperature-dependent upflow in the dimming region following a coronal mass ejection (CME) observed by the {\it Hinode} EUV Imaging Spectrometer (EIS). In this paper, we reanalyzed the observations along with previous work on this event, and provided boundary conditions for modeling. We found that the intensity in the dimming region dramatically drops within 30 minutes from the flare onset, and the dimming region reaches the equilibrium stage after $\sim$1 hour later. The temperature-dependent upflows were observed during the equilibrium stage by EIS. The cross sectional area of the fluxtube in the dimming region does not appear to expand significantly. From the observational constraints, we reconstructed the temperature-dependent upflow by using a new method which considers the mass and momentum conservation law, and demonstrated the height variation of plasma conditions in the dimming region. We found that a super radial expansion of the cross sectional area is required to satisfy the mass conservation and momentum equations. There is a steep temperature and velocity gradient of around 7 Mm from the solar surface. This result may suggest that the strong heating occurred above 7 Mm from the solar surface in the dimming region. We also showed that the ionization equilibrium assumption in the dimming region is violated especially in the higher temperature range.

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}).

We examine the magnetic field in NGC 4214, a nearby irregular galaxy, using multi-wavelength radio continuum polarization data from the Very Large Array. We find that the global radio continuum spectrum shows signs that free-free absorption and/or synchrotron losses may be important. The 3cm radio continuum morphology is similar to that of the Halpha, while the 20cm emission is more diffuse. We estimate that 50% of the radio continuum emission in the center of the galaxy is thermal. Our estimate of the magnetic field strength is $30\pm 9.5$ \uG\ in the center and $10\pm3$ \uG\ at the edges. We find that the hot gas, magnetic, and the gravitational pressures are all the same order of magnitude. Inside the central star forming regions, we find that the thermal and turbulent pressures of the HII regions dominate the pressure balance. We do not detect any significant polarization on size scales greater than 200 pc. We place an upper limit of 8 \uG\ on the uniform field strength in this galaxy. We suggest that the diffuse synchrotron region, seen to the north of the main body of emission at 20cm, is elongated due to a uniform magnetic field with a maximum field strength of 7.6 \uG. We find that, while the shear in NGC 4214 is comparable to that of the Milky Way, the supernova rate is half that of the Milky Way and suggest that the star formation episode in NGC 4214 needs additional time to build up enough turbulence to drive an $\alpha-\omega$ dynamo.

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.

Plasma experiments in laboratory settings offer unique opportunities to address fundamental aspects of the solar dynamo and magnetism in the solar atmosphere. We argue here that ground-based laboratory experiments have direct connections to NASA based missions and NSF programs, and that a small investment in laboratory heliophysics may have a high payoff. We advocate for broad involvement in community-scale plasma experiments.

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.

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 \AA. 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.

We present a study of the vertical magnetic field of the Milky Way towards the Galactic poles, determined from observations of Faraday rotation toward more than 1000 polarized extragalactic radio sources at Galactic latitudes |b| > 77 degs, using the Westerbork Radio Synthesis Telescope and the Australia Telescope Compact Array. We find median rotation measures (RMs) of 0.0 +/- 0.5 rad/m^2 and +6.3 +/- 0.7 rad/m^2 toward the north and south Galactic poles, respectively, demonstrating that there is no coherent vertical magnetic field in the Milky Way at the Sun's position. If this is a global property of the Milky Way's magnetism, then the lack of symmetry across the disk rules out pure dipole or quadrupole geometries for the Galactic magnetic field. The angular fluctuations in RM seen in our data show no preferred scale within the range ~ 0.1 to 25 degs. The observed standard deviation in RM of ~ 9 rad/m^2 then implies an upper limit of ~1microGauss on the strength of the random magnetic field in the warm ionized medium at high Galactic latitudes.

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.

(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.

We present photometry of the outer star clusters in NGC 1275, the brightest galaxy in the Perseus cluster. The observations were taken using the Hubble Space Telescope Advanced Camera for Surveys. We focus on two stellar regions in the south and south-east, far from the nucleus of the low velocity system (~22 kpc). These regions of extended star formation trace the H alpha filaments, drawn out by rising radio bubbles. In both regions bimodal distributions of colour (B-R)_0 against magnitude are apparent, suggesting two populations of star clusters with different ages; most of the H alpha filaments show no detectable star formation. The younger, bluer population is found to be concentrated along the filaments while the older population is dispersed evenly about the galaxy. We construct colour-magnitude diagrams and derive ages of at most 10^8 years for the younger population, a factor of 10 younger than the young population of star clusters in the inner regions of NGC 1275. We conclude that a formation mechanism or event different to that for the young inner population is needed to explain the outer star clusters and suggest that formation from the filaments, triggered by a buoyant radio bubble either rising above or below these filaments, is the most likely mechanism.

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.

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.

The recent demonstration of current singularity formation by Low et al. assumes that potential fields will remain potential under simple expansion or compression (Low 2006, 2007; Janse & Low 2009). An explicit counterexample to their key assumption is constructed. Our findings suggest that their results may need to be reconsidered.

We discuss the outstanding issues of the interstellar medium which will depend on the application of knowledge from plasma physics. We particularly advocate attention to recent developments in experimental plasma physics, and urge that the astronomical community consider support of these experiments in the next decade.

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.

The major questions relevant to star and planet formation are: What controls the rate, efficiency, spatial clustering, multiplicity, and initial mass function of star formation, now and in the past? What are the major feedback mechanisms through which star formation affects its environment? What controls the formation and orbital parameters of planets, especially terrestrial planets? These questions cannot be fully addressed without understanding key magnetohydrodynamics (MHD) and plasma physics processes. Although some of these basic problems have long been considered intractable, attacking them through a combination of laboratory experiment, theory, and numerical simulation is now feasible, and would be fruitful. Achieving a better understanding of these processes is critical to interpreting observations, and will form an important component of astrophysical models. These models in turn will serve as inputs to other areas of astrophysics, e.g. cosmology and galaxy formation.

In this White Paper, we emphasize the need for and the important role of plasma astrophysics in the studies of formation, evolution of, and feedback by Active Galaxies. We make three specific recommendations: 1) We need to significantly increase the resolution of VLA, perhaps by building an EVLA-II at a modest cost. This will provide the angular resolution to study jets at kpc scales, where, for example, detailed Faraday rotation diagnosis can be done at 1GHz transverse to jets; 2) We need to build coordinated programs among NSF, NASA, and DOE to support laboratory plasma experiments (including liquid metal) that are designed to study key astrophysical processes, such as magneto-rotational instability (origin of angular momentum transport), dynamo (origin of magnetic fields), jet launching and stability. Experiments allowing access to relativistic plasma regime (perhaps by intense lasers and magnetic fields) will be very helpful for understanding the stability and dissipation physics of jets from Supermassive Black Holes; 3) Again through the coordinated support among the three Agencies, we need to invest in developing comprehensive theory and advanced simulation tools to study the accretion disks and jets in relativistic plasma physics regime, especially in connecting large scale fluid scale phenomena with relativistic kinetic dissipation physics through which multi-wavelength radiation is produced.

The giant elliptical galaxy NGC 1275, at the centre of the Perseus cluster, is surrounded by a well-known giant nebulosity of emission-line filaments, which are plausibly about >10^8 yr old. The filaments are dragged out from the centre of the galaxy by the radio bubbles rising buoyantly in the hot intracluster gas before later falling back. They act as dramatic markers of the feedback process by which energy is transferred from the central massive black hole to the surrounding gas. The mechanism by which the filaments are stabilized against tidal shear and dissipation into the surrounding 4x10^7 K gas has been unclear. Here we report new observations that resolve thread-like structures in the filaments. Some threads extend over 6 kpc, yet are only 70 pc wide. We conclude that magnetic fields in the threads, in pressure balance with the surrounding gas, stabilize the filaments, so allowing a large mass of cold gas to accumulate and delay star formation.

We review the literature concerning how the cosmic magnetic fields pervading nearly all galaxies actually got started. some observational evidence involves the chemical abundance of the light elements Be and B, while another one is based on strong magnetic fields seen in high red shift galaxies. Seed fields, whose strength is of order 10^{-20} gauss, easily sprung up in the era preceding galaxy formation. Several mechanisms are proposed to amplify these seed fields to microgauss strengths. The standard mechanism is the Alpha-Omega dynamo theory. It has a major difficulty that makes unlikely to provide the sole origin. The difficulty is rooted in the fact that the total flux is constant. This implies that flux must be removed from the galactic discs. This requires that the field and flux be separated, for otherwise interstellar mass must be removed from the deep galactic gravitational and then their strength increased by the alpha omega theory.

We present numerical simulations of the growth and saturation of the Kelvin-Helmholtz instability in a compressible fluid layer with and without a weak magnetic field. In the absence of a magnetic field, the instability generates a single eddy which flattens the velocity profile, stabilizing it against further perturbations. Adding a weak magnetic field - weak in the sense that it has almost no effect on the linear instability - leads to a complex flow morphology driven by MHD forces and to enhanced broadening of the layer, due to Maxwell stresses. We corroborate earlier studies which showed that magnetic fields destroy the large scale eddy structure through periodic cycles of windup and resistive decay, but we show that the rate of decay decreases with decreasing plasma resistivity, at least within the range of resistivity accessible to our simulations. Magnetization increases the efficiency of momentum transport, and the transport increases with decreasing resistivity.

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.

Magnetic fields are an important component of the interstellar medium, especially in low-mass galaxies like irregulars where the magnetic pressure may be significant. However, few irregular galaxies have observed magnetic field structures. Using the VLA, the GBT, and the ATCA, we have observed several irregular galaxies in the radio continuum to determine their magnetic field structures. Here we report on our results for the galaxies NGC 4214 and NGC 1569.

It is well known that the radial displacement of the $m=1$ internal kink mode in a periodic screw pinch has a steep jump at the resonant surface where $\mathbf{k}\cdot\mathbf{B}=0$. In a line-tied system, relevant to solar and astrophysical plasmas, the resonant surface is no longer a valid concept. It is then of interest to see how line-tying alters the aforementioned result for a periodic system. If the line-tied kink also produces a steep gradient, corresponding to a thin current layer, it may lead to strong resistive effects even with weak dissipation. Numerical solution of the eigenmode equations shows that the fastest growing kink mode in a line-tied system still possesses a jump in the radial displacement at the location coincident with the resonant surface of the fastest growing mode in the periodic counterpart. However, line-tying thickens the inner layer and slows down the growth rate. As the system length $L$ approaches infinity, both the inner layer thickness and the growth rate approach the periodic values. In the limit of small $\epsilon\sim B_{\phi}/B_{z}$, the critical length for instability $L_{c}\sim\epsilon^{-3}$. The relative increase in the inner layer thickness due to line-tying scales as $\epsilon^{-1}(L_{c}/L)^{2.5}$.

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.

We have calculated the radiation field, dust grain temperatures, and far infrared emissivity of numerical models of turbulent molecular clouds. When compared to a uniform cloud of the same mean optical depth, most of the volume inside the turbulent cloud is brighter, but most of the mass is darker. There is little mean attenuation from center to edge, and clumping causes the radiation field to be somewhat bluer. There is also a large dispersion, typically by a few orders of magnitude, of all quantities relative to their means. However, despite the scatter, the 850 micron emission maps are well correlated with surface density. The fraction of mass as a function of intensity can be reproduced by a simple hierarchical model of density structure.

We consider the stability of a configuration consisting of a vertical magnetic field in a planar flow on elliptical streamlines in ideal hydromagnetics. In the absence of a magnetic field the elliptical flow is universally unstable (the ``elliptical instability''). We find this universal instability persists in the presence of magnetic fields of arbitrary strength, although the growthrate decreases somewhat. We also find further instabilities due to the presence of the magnetic field. One of these, a destabilization of Alfven waves, requires the magnetic parameter to exceed a certain critical value. A second, involving a mixing of hydrodynamic and magnetic modes, occurs for all magnetic-field strengths. These instabilities may be important in tidally distorted or otherwise elliptical disks. A disk of finite thickness is stable if the magnetic fieldstrength exceeds a critical value, similar to the fieldstrength which suppresses the magnetorotational instability.

We present calculations of magnetic field evolution by the Hall effect and Ohmic decay in neutron star crusts. We first discuss when and for which neutron stars the Hall effect is important relative to Ohmic decay. We compute the evolution of an initial field distribution by Ohmic decay, and give approximate analytic formulas for both the surface and interior fields as a function of time. We then discuss the subsequent evolution by the Hall effect, giving several examples which illustrate how the field configuration evolves. Hall wave eigenfunctions are computed, including the effect of the large density change across the crust. We estimate the response of the crust to the magnetic stresses induced by Hall waves, and give a detailed discussion of the boundary conditions at the solid-liquid interface. Finally, we discuss the implications for the Hall cascade proposed by Goldreich & Reisenegger.

Ambient interstellar material may become entrained in outflows from massive stars as a result of shear flow instabilities. We study the linear theory of the Kelvin - Helmholtz instability, the simplest example of shear flow instability, in a partially ionized medium. We model the interaction as a two fluid system (charged and neutral) in a planar geometry. Our principal result is that for much of the relevant parameter space, neutrals and ions are sufficiently decoupled that the neutrals are unstable while the ions are held in place by the magnetic field. Thus, we predict that there should be a detectably narrower line profile in ionized species tracing the outflow compared with neutral species since ionized species are not participating in the turbulent interface with the ambient ISM. Since the magnetic field is frozen to the plasma, it is not tangled by the turbulence in the boundary layer.

Increasingly sophisticated observational tools and techniques are now being developed for probing the nature of interstellar turbulence. At the same time, theoretical advances in understanding the nature of turbulence and its effects on the structure of the ISM and on star formation are occurring at a rapid pace, aided in part by numerical simulations. These increased capabilities on both fronts open new opportunities for strengthening the links between observation and theory,and for meaningful comparisons between the two.

Under ideal MHD conditions the magnetic field strength should be correlated with density in the interstellar medium (ISM). However, observations indicate that this correlation is weak. Ambipolar diffusion can decrease the flux-to-mass ratio in weakly ionized media; however, it is generally thought to be too slow to play a significant role in the ISM except in the densest molecular clouds. Turbulence is often invoked in astrophysical problems to increase transport rates above the (very slow) laminar values predicted by kinetic theory. We describe a series of numerical experiments addressing the problem of turbulent transport of magnetic fields in weakly ionized gases. We show, subject to various geometrical and physical restrictions, that turbulence in a weakly ionized medium rapidly diffuses the magnetic flux to mass ratio through the buildup of appreciable ion-neutral drifts on small scales. These results are applicable to the fieldstrength - density correlation in the ISM, as well as the merging of flux systems such as protostar and accretion disk fields or protostellar jets with ambient matter, and the vertical transport of galactic magnetic fields.

There is evidence that cosmic rays were present in galaxies at moderately high redshift. This suggests that magnetic fields were also present. If cosmic rays and magnetic fields must always be close to equipartition, as they are to an order of magnitude within the local universe, this would provide a powerful constraint on theories of the origin and evolution of magnetic fields in galaxies. We evaluate current theories of cosmic ray acceleration and confinement at very low magnetic fieldstrength, and find that the properties of small scale turbulence are fundamentally changed in the presence of cosmic rays. As a result, magnetic fields several orders of magnitude weaker than present galactic fields can accelerate and retain a population of relativistic cosmic rays, provided that the fields are coherent on lengthscales greater than a cosmic ray gyroradius.

Magnetic reconnection plays an essential role in the generation and evolution of astrophysical magnetic fields. The best tested and most robust reconnection theory is that of Parker and Sweet. According to this theory, the reconnection rate scales with magnetic diffusivity lambda as lambda^0.5. In the interstellar medium, the Parker-Sweet reconnection rate is far too slow to be of interest. Thus, a mechanism for fast reconnection seems to be required. We have studied the magnetic merging of two oppositely directed flux systems in weakly ionized, but highly conducting, compressible gas. In such systems, ambipolar diffusion steepens the magnetic profile, leading to a thin current sheet. If the ion pressure is small enough, and the recombination of ions is fast enough, the resulting rate of magnetic merging is fast, and independent of lambda. Slow recombination or sufficiently large ion pressure leads to slower merging which scales with lambda as lambda^0.5. We derive a criterion for distinguishing these two regimes, and discuss applications to the weakly ionized ISM and to protoplanetary accretion disks.

The generation and evolution of astrophysical magnetic fields occurs largely through the action of turbulence. In many situations, the magnetic field is strong enough to influence many important properties of turbulence itself. Numerical simulation of magnetized turbulence is especially challenging in the astrophysical regime because of the high magnetic Reynolds numbers involved, but some aspects of this difficulty can be avoided in weakly ionized systems.

The interstellar magnetic field strength and gas density are observed to be correlated, but there is a large dispersion in this relation. In particular, the magnetic field is often observed to be weaker than expected. It is usually assumed that the field is frozen into the gas, except at very low ionization fraction, in which case ion-neutral drift, or ambipolar diffusion, leads to appreciable slip of the field and tends to make it more uniform. We show that turbulence enhances the rate of ambipolar drift, and that this may explain the surprisingly low fieldstrengths sometimes seen in dense interstellar gas.

The dynamical state of star-forming molecular clouds cannot be understood without determining the structure and strength of their magnetic fields. Measurements of polarized far-infrared radiation from thermally aligned dust grains are used to map the orientation of the field and estimate its strength, but the accuracy of the results has remained in doubt. In order to assess the reliability of this method, we apply it to simulated far-infrared polarization maps derived from three-dimensional simulations of supersonic magnetohydrodynamical turbulence, and compare the estimated values to the known magnetic field strengths in the simulations. We investigate the effects of limited telescope resolution and self-gravity on the structure of the maps. Limited observational resolution affects the field structure such that small scale variations can be completely suppressed, thus giving the impression of a very homogeneous field. The Chandrasekhar-Fermi method of estimating the mean magnetic field in a turbulent medium is tested, and we suggest an extension to measure the rms field. Both methods yield results within a factor of 2 for field strengths typical of molecular clouds, with the modified version returning more reliable estimates for slightly weaker fields. However, neither method alone works well for very weak fields, missing them by a factor of up to 150. Taking the geometric mean of both methods estimates even the weakest fields accurately within a factor of 2.5. Limited telescope resolution leads to a systematic overestimation of the field strengths for all methods. We discuss the effects responsible for this overestimation and show how to extract information on the underlying (turbulent) power spectrum.

We investigate whether the magnetic field of an accreting neutron star may be diamagnetically screened by the accreted matter. We assume the freshly accumulated material is unmagnetized, and calculate the rate at which the intrinsic stellar magnetic flux is transported into it by Ohmic diffusion. We calculate the one-dimensional steady-state magnetic field profiles, and show that the magnetic field strength decreases as one moves up through the outer crust and ocean by roughly (Mdot/0.02 Mdot_Edd) orders of magnitude, where Mdot is the accretion rate and Mdot_Edd the Eddington accretion rate. We show that buoyancy instabilities set a limit to the strength of any buried field of roughly 10^10-10^11 G. Our results show that magnetic screening is ineffective for Mdot<0.01 Mdot_Edd, so that, no matter how the accreted material joins onto the star, the underlying stellar field should always be evident. In this respect, we point out the only known persistently-pulsing accreting X-ray millisecond pulsar, SAX J1808.4-3658, has an accretion rate of 10^-3 Mdot_Edd, far below the regime where magnetic screening can play a role. Most steadily accreting neutron stars in low-mass X-ray binaries in our Galaxy accrete at rates where screening would be effective if the simplified magnetic and accretion geometry we adopt were correct. If screened, then the underlying field will emerge after accretion halts, on a timescale of only 100--1000 years, set by the Ohmic diffusion time across the outer crust. It thus seems unlikely that screening alone can explain the low magnetic fields of the millisecond radio pulsars.

We investigate the generation of magnetic field by the Biermann battery in cosmological ionization fronts, using new simulations of the reionization of the universe by stars in protogalaxies. Two mechanisms are primarily responsible for magnetogenesis: i) the breakout of I-fronts from protogalaxies, and ii) the propagation of I-fronts through the high density neutral filaments which are part of the cosmic web. The first mechanism is dominant prior to overlapping of ionized regions (z ~ 7), whereas the second continues to operate even after that epoch. However, after overlap the field strength increase is largely due to the gas compression occurring as cosmic structures form. As a consequence, the magnetic field at z ~ 5 closely traces the gas density, and it is highly ordered on megaparsec scales. The mean mass-weighted field strength is B_0 ~ 10^{-19} G in the simulation box. There is a relatively well-defined, nearly linear correlation between B_0 and the baryonic mass of virialized objects, with B_0 ~ 10^{-18} G in the most massive objects (M ~ 10^9 M_sun) in our simulations. This is a lower limit, as lack of numerical resolution prevents us from following small scale dynamical processes which could amplify the field in protogalaxies. Although the field strengths we compute are probably adequate as seed fields for a galactic dynamo, the field is too small to have had significant effects on galaxy formation, on thermal conduction, or on cosmic ray transport in the intergalactic medium. It could, however, be observed in the intergalactic medium through innovative methods based on analysis of gamma-ray burst photon arrival times.

We simulate fragmentation and gravitational collapse of cold, magnetized molecular clouds. We explore the nonlinear development of an instability mediated by ambipolar diffusion, in which the collapse rate is intermediate to fast gravitational collapse and slow quasistatic collapse. Initially uniform stable clouds fragment into elongated clumps with masses largely determined by the cloud temperature, but substantially larger than the thermal Jeans mass. The clumps are asymmetric, with significant rotation and vorticity, and lose magnetic flux as they collapse. The clump shapes, intermediate collapse rates, and infall profiles may help explain observations not easily fit by contemporary slow or rapid collapse models.

We describe an approach for finding the eigenfrequencies of solar acoustic modes (p modes) in a convective envelope in the WKB limit. This approximation restricts us to examining the effects of fluid motions which are large compared to the mode wavelength, but allows us to treat the three-dimensional mode as a localized ray. The method of adiabatic switching is then used to investigate the frequency shifts resulting from simple perturbations to a polytropic model of the convection zone as well as from two basic models of a convective cell. We find that although solely depth-dependent perturbations can give frequency shifts which are first order in the strength of the perturbation, models of convective cells generate downward frequency shifts which are second order in the perturbation strength. These results may have implications for resolving the differences between eigenfrequencies derived from solar models and those found from helioseismic observations.

We develop an analytical theory of reconnection between colliding, twisted magnetic flux tubes. Our analysis is restricted to direct collisions between parallel tubes and is based on the collision dynamics worked out by Bogdan (1984). We show that there is a range of collision velocities for which neutral point reconnection of the Parker-Sweet type can occur, and a smaller range for which reconnection leads to coalescence. Mean velocities within the solar convection zone are probably significantly greater than the upper limit for coalescence. This suggests that the majority of flux tube collisions do not result in merging, unless the frictional coupling of the tubes to the background flow is extremely strong.