Numerical studies of gas accretion onto supermassive black hole binaries (SMBHBs) have generally been limited to conditions where the circumbinary disk (CBD) is 10-100 times thicker than expected for disks in active galactic nuclei (AGN). This discrepancy arises from technical limitations, and also from publication bias toward replicating fiducial numerical models. Here we present the first systematic study of how the binary's orbital evolution varies with disk scale height. We report three key results: (1) Binary orbital evolution switches from outspiralling for warm disks (aspect ratio ~0.1), to inspiralling for more realistic cooler, thinner disks at a critical aspect ratio ~0.04, corresponding to orbital Mach number ~25. (2) The net torque on the binary arises from a competition between positive torque from gas orbiting close to the black holes, and negative torque from the inner edge of the CBD, which is denser for thinner disks. This leads to increasingly negative net torques on the binary for increasingly thin disks. (3) The accretion rate is modestly suppressed with increasing Mach number. We discuss how our results may influence modeling of the nano-Hz gravitational wave background, as well as estimates of the LISA merger event rate.

Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with previous studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary towards equal mass. Accretion variability depends sensitively on the numerical sink prescription, but the torque and relative accretion onto each component do not depend on the sink timescale. Positive torque and highly variable accretion occurs only for mass ratios greater than around $0.05$. This means that for mass ratios below $0.05$, the binary would migrate inward until the secondary accreted sufficient mass, after which it would execute a U-turn and migrate outward. We explore a range of viscosities, from $\alpha = 0.03$ to $\alpha = 0.15$, and find that this outward torque is proportional to the viscous torque, simply proportional to viscosity in this range. Dependence of accretion and torque on mass ratio is explored in detail, densely sampling mass ratios between $0.01$ and unity. For mass ratio $q > 0.6$, accretion variability is found to exhibit a distinct sawtooth pattern, typically with a five-orbit cycle that provides a "smoking gun" prediction for variable quasars observed over long periods, as a potential means to confirm the presence of a binary.

We present Chandra and VLA observations of GW170817 at ~521-743 days post merger, and a homogeneous analysis of the entire Chandra data set. We find that the late-time non-thermal emission follows the expected evolution from an off-axis relativistic jet, with a steep temporal decay $F_{\nu}\propto t^{-1.95\pm0.15}$ and a simple power-law spectrum $F_{\nu}\propto \nu^{-0.575\pm0.007}$. We present a new method to constrain the merger environment density based on diffuse X-ray emission from hot plasma in the host galaxy and we find $n\le 9.6 \times 10^{-3}\,\rm{cm^{-3}}$. This measurement is independent from inferences based on the jet afterglow modeling and allows us to partially solve for model degeneracies. The updated best-fitting model parameters with this density constraint are a fireball kinetic energy $E_0 = 1.5_{-1.1}^{+3.6}\times 10^{49}\,\rm{erg}$ ($E_{iso}= 2.1_{-1.5}^{+6.4}\times10^{52}\, \rm{erg}$), jet opening angle $\theta_{0}= 5.9^{+1.0}_{-0.7}\,\rm{deg}$ with characteristic Lorentz factor $\Gamma_j = 163_{-43}^{+23}$, expanding in a low-density medium with $n_0 = 2.5_{-1.9}^{+4.1} \times 10^{-3}\, \rm{cm^{-3}}$ and viewed $\theta_{obs} = 30.4^{+4.0}_{-3.4}\, \rm{deg}$ off-axis. The synchrotron emission originates from a power-law distribution of electrons with $p=2.15^{+0.01}_{-0.02}$. The shock microphysics parameters are constrained to $\epsilon_{\rm{e}} = 0.18_{-0.13}^{+0.30}$ and $\epsilon_{\rm{B}}=2.3_{-2.2}^{+16.0} \times 10^{-3}$. We investigate the presence of X-ray flares and find no statistically significant evidence of $\ge2.5\sigma$ of temporal variability at any time. Finally, we use our observations to constrain the properties of synchrotron emission from the deceleration of the fastest kilonova ejecta with energy $E_k^{KN}\propto (\Gamma\beta)^{-\alpha}$ into the environment, finding that shallow stratification indexes $\alpha\le6$ are disfavored.

We present a revised and complete optical afterglow light curve of the binary neutron star merger GW170817, enabled by deep Hubble Space Telescope (HST) F606W observations at $\approx\!584$ days post-merger, which provide a robust optical template. The light curve spans $\approx 110-362$ days, and is fully consistent with emission from a relativistic structured jet viewed off-axis, as previously indicated by radio and X-ray data. Combined with contemporaneous radio and X-ray observations, we find no spectral evolution, with a weighted average spectral index of $\langle \beta \rangle = -0.583 \pm 0.013$, demonstrating that no synchrotron break frequencies evolve between the radio and X-ray bands over these timescales. We find that an extrapolation of the post-peak temporal slope of GW170817 to the luminosities of cosmological short GRBs matches their observed jet break times, suggesting that their explosion properties are similar, and that the primary difference in GW170817 is viewing angle. Additionally, we place a deep limit on the luminosity and mass of an underlying globular cluster of $L \lesssim 6.7 \times 10^{3}\,L_{\odot}$, or $M \lesssim 1.3 \times 10^{4}\,M_{\odot}$, at least 4 standard deviations below the peak of the globular cluster mass function of the host galaxy, NGC4993. This limit provides a direct and strong constraint that GW170817 did not form and merge in a globular cluster. As highlighted here, HST (and soon JWST) enables critical observations of the optical emission from neutron star merger jets and outflows.

We investigate the self-consistent particle acceleration physics associated with the development of the kink instability (KI) in nonrelativistic, electron-ion plasma jets. Using 3D fully kinetic particle-in-cell (PIC) simulations, we show that the KI efficiently converts the initial toroidal magnetic field energy into energetic ions. The accelerated ions form a nonthermal power-law tail in the energy spectrum, containing $\simeq10\%$ of the initial magnetic field energy, and with the maximum ion energy extending to the confinement energy of the jet. We find that the ions are efficiently accelerated by the concerted action of the motional electric field and highly tangled magnetic field that develop in the nonlinear phase of the KI: fast curvature drift motions of ions across magnetic field lines enable their acceleration along the electric field. We further investigate the role of Coulomb collisions on the ion acceleration efficiency, and identify the collisional threshold above which nonthermal ion acceleration is suppressed. Our results reveal how energetic ions may result from unstable nonrelativistic plasma jets in space and astrophysics, and provide constraints on the plasma conditions required to reproduce this acceleration mechanism in laboratory experiments.

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 field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

Magnetar flares excite strong Alfv\'{e}n waves in the magnetosphere of the neutron star. The wave energy can (1) dissipate in the magnetosphere, (2) convert to "fast modes" and possibly escape, and (3) penetrate the neutron star crust and dissipate there. We examine and compare the three options. Particularly challenging are nonlinear interactions between strong waves, which develop a cascade to small dissipative scales. This process can be studied in the framework of force-free electrodynamics (FFE). We perform three-dimensional FFE simulations to investigate Alfv\'{e}n wave dissipation, how long it takes, and how it depends on the initial wave amplitude on the driving scale. In the simulations, we launch two large Alfv\'{e}n wave packets that keep bouncing on closed magnetic field lines and collide repeatedly until the full turbulence spectrum develops. Besides dissipation due to the turbulent cascade, we find that in some simulations spurious energy losses occur immediately in the first collisions. This effect occurs in special cases where the FFE description breaks. It is explained with a simple one-dimensional model, which we examine in both FFE and full magnetohydrodynamic settings. We find that magnetospheric dissipation through nonlinear wave interactions is relatively slow and more energy is drained into the neutron star. The wave energy deposited into the star is promptly dissipated through plastic crustal flows induced at the bottom of the liquid ocean, and a fraction of the generated heat is radiated from the stellar surface.

Relativistic magnetized jets from active galaxies are among the most powerful cosmic accelerators, but their particle acceleration mechanisms remain a mystery. We present a new acceleration mechanism associated with the development of the helical kink instability in relativistic jets, which leads to the efficient conversion of the jet's magnetic energy into nonthermal particles. Large-scale three-dimensional ab initio simulations reveal that the formation of highly tangled magnetic fields and a large-scale inductive electric field throughout the kink-unstable region promotes rapid energization of the particles. The energy distribution of the accelerated particles develops a well-defined power-law tail extending to the radiation-reaction limited energy in the case of leptons, and to the confinement energy of the jet in the case of ions. When applied to the conditions of well-studied bright knots in jets from active galaxies, this mechanism can account for the spectrum of synchrotron and inverse Compton radiating particles, and offers a viable means of accelerating ultra-high-energy cosmic rays to $10^{20}$ eV.

We examine the possible role of turbulence in feeding the emission of gamma-ray bursts (GRBs). Turbulence may develop in a GRB jet as the result of hydrodynamic or current-driven instabilities. The jet carries dense radiation and the turbulence cascade can be damped by Compton drag, passing kinetic fluid energy to photons through scattering. We identify two regimes of turbulence dissipation: (1) "Viscous" - the turbulence cascade is Compton damped on a scale $\ell_{\rm damp}$ greater than the photon mean free path $\ell_\star$. Then turbulence energy is passed to photons via bulk Comptonization by smooth shear flows on scale $\ell_\star<\ell_{\rm damp}$. (2) "Collisionless" - the cascade avoids Compton damping and extends to microscopic plasma scales much smaller than $\ell_\star$. The collisionless dissipation energizes plasma particles, which radiate the received energy; how the dissipated power is partitioned between particles needs further investigation with kinetic simulations. We show that the dissipation regime switches from viscous to collisionless during the jet expansion, at a critical value of the jet optical depth which depends on the amplitude of turbulence. Turbulent GRB jets are expected to emit nonthermal photospheric radiation. Our analysis also suggests revisions of turbulent Comptonization in black hole accretion disks discussed in previous works.

We present synthetic radio images of the GRB170817A afterglow, computed from moving-mesh hydrodynamic simulations of binary neutron star merger outflows. Having expanded for nearly a year, the merger remnant is expected to subtend roughly 5 milli-arcseconds (mas) on the sky, potentially resolvable by very long baseline radio imaging techniques. Any observations revealing the radio centroid to be offset from the line-of-site to the merger would be the smoking gun of a jetted outflow. However, our results indicate that a measurement of the centroid position alone cannot independently determine whether that jet escaped successfully from the merger debris cloud, or was "choked," yielding a quasi-spherical explosion. We find that in both scenarios, the centroid exhibits superluminal proper motion away from the merger site at roughly 4 - 10 mas per day for at least the first 300 days. We argue that a successful strategy for differentiating among the explosion models will need to include multiple observations over the coming months - years. In particular, we find the time at which the centroid attains its maximum offset, and begins heading back toward the merger site, is considerably later if the jet was choked. Detecting a reversal of the centroid trajectory earlier than 600 days would uniquely identify a successful jet. Our results indicate the source might be resolved using VLBI radio observing techniques with roughly 1 mas resolution starting at roughly 400 days post-merger, and that the the angular extent of a successful jet is significantly smaller than that of a choked jet (4.5 versus 7 mas respectively).

We present new observations of the binary neutron star merger GW170817 at $\Delta t\approx 220-290$ days post-merger, at radio (Karl G. Jansky Very Large Array; VLA), X-ray (Chandra X-ray Observatory) and optical (Hubble Space Telescope; HST) wavelengths. These observations provide the first evidence for a turnover in the X-ray light curve, mirroring a decline in the radio emission at $\gtrsim5\sigma$ significance. The radio-to-X-ray spectral energy distribution exhibits no evolution into the declining phase. Our full multi-wavelength dataset is consistent with the predicted behavior of our previously published models of a successful structured jet expanding into a low-density circumbinary medium, but pure cocoon models with a choked jet cannot be ruled out. If future observations continue to track our predictions, we expect that the radio and X-ray emission will remain detectable until $\sim 1000$ days post-merger.

We present numerical simulations of energetic flows propagating through the debris cloud of a binary neutron star (BNS) merger. Starting from the scale of the central engine, we use a moving-mesh hydrodynamics code to simulate the complete dynamical evolution of the produced relativistic jets. We compute synchrotron emission directly from the simulations and present multi-band light curves of the early (sub-day) through late (weeks to years) afterglow stages. Our work systematically compares two distinct models for the central engine, referred to as the narrow and wide engine scenario, which is associated with a successful structured jet and a quasi-isotropic explosion respectively. Both engine models naturally evolve angular and radial structure through hydrodynamical interaction with the merger debris cloud. They both also result in a relativistic blast wave capable of producing the observed multi-band afterglow data. However, we find that the narrow and wide engine scenario might be differentiated by a new emission component that we refer to as a merger flash. This component is a consequence of applying the synchrotron radiation model to the shocked optically thin merger cloud. Such modeling is appropriate if injection of non-thermal electrons is sustained in the breakout relativistic shell, for example by internal shocks or magnetic reconnection. The rapidly declining signature may be detectable for future BNS mergers during the first minutes to day following the GW chirp. Furthermore, its non-detection for the GRB170817A event may disfavor the wide, quasi-isotropic explosion model.

We report deep Chandra, HST and VLA observations of the binary neutron star event GW170817 at $t<160$ d after merger. These observations show that GW170817 has been steadily brightening with time and might have now reached its peak, and constrain the emission process as non-thermal synchrotron emission where the cooling frequency $\nu_c$ is above the X-ray band and the synchrotron frequency $\nu_m$ is below the radio band. The very simple power-law spectrum extending for eight orders of magnitude in frequency enables the most precise measurement of the index $p$ of the distribution of non-thermal relativistic electrons $N(\gamma)\propto \gamma^{-p}$ accelerated by a shock launched by a NS-NS merger to date. We find $p=2.17\pm0.01$, which indicates that radiation from ejecta with $\Gamma\sim3-10$ dominates the observed emission. While constraining the nature of the emission process, these observations do \emph{not} constrain the nature of the relativistic ejecta. We employ simulations of explosive outflows launched in NS ejecta clouds to show that the spectral and temporal evolution of the non-thermal emission from GW170817 is consistent with both emission from radially stratified quasi-spherical ejecta traveling at mildly relativistic speeds, \emph{and} emission from off-axis collimated ejecta characterized by a narrow cone of ultra-relativistic material with slower wings extending to larger angles. In the latter scenario, GW170817 harbored a normal SGRB directed away from our line of sight. Observations at $t\le 200$ days are unlikely to settle the debate as in both scenarios the observed emission is effectively dominated by radiation from mildly relativistic material.

We study the magnetic field evolution of an isolated spiral galaxy, using isolated Milky Way-mass galaxy formation simulations and a novel prescription for magnetohydrodynamic (MHD) supernova feedback. Our main result is that a galactic dynamo can be seeded and driven by supernova explosions, resulting in magnetic fields whose strength and morphology is consistent with observations. In our model, supernovae supply thermal energy, and a low level magnetic field along with their ejecta. The thermal expansion drives turbulence, which serves a dual role by efficiently mixing the magnetic field into the interstellar medium, and amplifying it by means of turbulent dynamo. The computational prescription for MHD supernova feedback has been implemented within the publicly available ENZO code, and is fully described in this paper. This improves upon ENZO's existing modules for hydrodynamic feedback from stars and active galaxies. We find that the field attains $\mu G$-levels over Gyr-time scales throughout the disk. The field also develops large-scale structure, which appears to be correlated with the disk's spiral arm density structure. We find that seeding of the galactic dynamo by supernova ejecta predicts a persistent correlation between gas metallicity and magnetic field strength. We also generate all-sky maps of the Faraday rotation measure from the simulation-predicted magnetic field, and present a direct comparison with observations.

We present a model for magnetic energy dissipation in a pulsar wind nebula. Better understanding of this process is required to assess the likelihood that certain astrophysical transients may be powered by the spin-down of a "millisecond magnetar." Examples include superluminous supernovae, gamma-ray bursts, and anticipated electromagnetic counterparts to gravitational wave detections of binary neutron star coalescence. Our model leverages recent progress in the theory of turbulent magnetic relaxation to specify a dissipative closure of the stationary magnetohydrodynamic (MHD) wind equations, yielding predictions of the magnetic energy dissipation rate throughout the nebula. Synchrotron losses are treated self-consistently. To demonstrate the model's efficacy, we show that it can reproduce many features of the Crab Nebula, including its expansion speed, radiative efficiency, peak photon energy, and mean magnetic field strength. Unlike ideal MHD models of the Crab (which lead to the so-called sigma problem) our model accounts for the transition from ultra to weakly magnetized plasma flow, and for the associated heating of relativistic electrons. We discuss how the predicted heating rates may be utilized to improve upon models of particle transport and acceleration in pulsar wind nebulae. We also discuss implications for the Crab Nebula's gamma-ray flares, and point out potential modifications to models of astrophysical transients invoking the spin-down of a millisecond magnetar.

Many powerful and variable gamma-ray sources, including pulsar wind nebulae, active galactic nuclei and gamma-ray bursts, seem capable of accelerating particles to gamma-ray emitting energies efficiently over very short time scales. These are likely due to rapid dissipation of electromagnetic energy in a highly magnetized, relativistic plasma. In order to understand the generic features of such processes, we have investigated simple models based on relaxation of unstable force-free magnetostatic equilibria. In this work, we make the connection between the corresponding plasma dynamics and the expected radiation signal, using 2D particle-in-cell simulations that self-consistently include synchrotron radiation reaction. We focus on the lowest order unstable force-free equilibrium in a 2D periodic box. We find that rapid variability, with modest apparent radiation efficiency as perceived by a fixed observer, can be produced during the evolution of the instability. The "flares" are accompanied by an increased polarization degree in the high energy band, with rapid variation in the polarization angle. Furthermore, the separation between the acceleration sites and the synchrotron radiation sites for the highest energy particles facilitates acceleration beyond the synchrotron radiation reaction limit. We also discuss the dynamical consequences of radiation reaction, and some astrophysical applications of this model. Our current simulations with numerically tractable parameters are not yet able to reproduce the most dramatic gamma-ray flares, e.g., from Crab Nebula. Higher magnetization studies are promising and will be carried out in the future.

We present the results of particle-in-cell numerical pair plasma simulations of relativistic 2D magnetostatic equilibria known as the 'ABC' fields. In particular, we focus on the lowest-order unstable configuration consisting of two minima and two maxima of the magnetic vector potential. Breaking of the initial symmetry leads to exponential growth of the electric energy and to the formation of two current layers, which is consistent with the picture of 'X-point collapse' first described by Syrovatskii. Magnetic reconnection within the layers heats a fraction of particles to very high energies. After the saturation of the linear instability, the current layers are disrupted and the system evolves chaotically, diffusing the particle energies in a stochastic second-order Fermi process leading to the formation of power-law energy distributions. The power-law slopes harden with the increasing mean magnetization, but they are significantly softer than those produced in simulations initiated from Harris-type layers. The maximum particle energy is proportional to the mean magnetization, which is attributed partly to the increase of the effective electric field and partly to the increase of the acceleration time scale. We describe in detail the evolving structure of the dynamical current layers, and report on the conservation of magnetic helicity. These results can be applied to highly magnetized astrophysical environments, where ideal plasma instabilities trigger rapid magnetic dissipation with efficient particle acceleration and flares of high-energy radiation.

We interpret $\gamma$-ray flares from the Crab Nebula as the signature of turbulence in the pulsar's electromagnetic outflow. Turbulence is triggered upstream by dynamical instability of the wind's oscillating magnetic field, and accelerates non-thermal particles. On impacting the wind termination shock, those particles emit a distinct synchrotron component $F_{\nu,\rm flare}$, which is constantly modulated by intermittency of the upstream plasma flow. Flares are observed when the high-energy cutoff of $F_{\nu,\rm flare}$ emerges above the fast-declining nebular emission around 0.1 - 1 GeV. Simulations carried out in the force-free electrodynamics approximation predict the striped wind to become fully turbulent well ahead of the wind termination shock, provided its terminal Lorentz factor is $\lesssim 10^4$.

Freely decaying relativistic force-free turbulence is studied for the first time. We initiate the magnetic field at a short wavelength and simulate its relaxation toward equilibrium on two and three dimensional periodic domains, in both helical and non-helical settings. Force-free turbulent relaxation is found to exhibit an inverse cascade in all settings, and in 3D to have a magnetic energy spectrum consistent with the Kolmogorov $5/3$ power law. 3D relaxations also obey the Taylor hypothesis; they settle promptly into the lowest energy configuration allowed by conservation of the total magnetic helicity. But in 2D, the relaxed state is a force-free equilibrium whose energy greatly exceeds the Taylor minimum, and which contains persistent force-free current layers and isolated flux tubes. We explain this behavior in terms of additional topological invariants that exist only in two dimensions, namely the helicity enclosed within each level surface of the magnetic potential function. The speed and completeness of turbulent magnetic free energy discharge could help account for rapidly variable gamma-ray emission from the Crab Nebula, gamma-ray bursts, blazars, and radio galaxies.

Jets associated with Active Galactic Nuclei (AGN) have been observed for almost a century, initially at optical and radio wavelengths. They are now widely accepted as "exhausts" produced electromagnetically by the central, spinning, massive black hole and its orbiting, accreting gas. Observations at X-ray and, especially, gamma-ray energies have transformed our understanding of how these jets evolve dynamically, accelerate electrons (and positrons) and radiate throughout the entire electromagnetic spectrum. Some new approaches to modeling the powerful and rapidly variable TeV emission observed from many blazars are sketched. Observations at the highest TeV energies, to which the High Altitude Water Cherenkov Gamma-Ray Observatory (HAWC) will contribute, promise crucial discrimination between rival models of AGN jets.

The merger of binary neutron stars (BNSs) can lead to large amplifications of the magnetic field due to the development of turbulence and instabilities in the fluid, such as the Kelvin-Helmholtz shear instability, which drive small-scale dynamo activity. In order to properly resolve such instabilities and obtain the correct magnetic field amplification, one would need to employ resolutions that are currently unfeasible in global general relativistic magnetohydrodynamic (GRMHD) simulations of BNS mergers. Here, we present a subgrid model that allows global simulations to take into account the small-scale amplification of the magnetic field which is caused by the development of turbulence during BNS mergers. Assuming dynamo saturation, we show that magnetar-level fields ($\sim 10^{16}\,{\rm G}$) can be easily reached, and should therefore be expected from the merger of magnetized BNSs. The total magnetic energy can reach values up to $\sim 10^{51}\,{\rm erg}$ and the post-merger remnant can therefore emit strong electromagnetic signals and possibly produce short gamma-ray bursts.

In order to understand the conditions which lead a highly magnetized, relativistic plasma to become unstable, and in such cases how the plasma evolves, we study a prototypical class of magnetostatic equilibria where the magnetic field satisfies $\nabla \times\mathbf B = \alpha \mathbf B$, where $\alpha$ is spatially uniform, on a periodic domain. Using numerical solutions we show that generic examples of such equilibria are unstable to ideal modes (including incompressible ones) which are marked by exponential growth in the linear phase. We characterize the unstable mode, showing how it can be understood in terms of merging magnetic and current structures, and explicitly demonstrate its instability using the energy principle. Following the nonlinear evolution of these solutions, we find that they rapidly develop regions with relativistic velocities and electric fields of comparable magnitude to the magnetic field, liberating magnetic energy on dynamical timescales and eventually settling into a configuration with the largest allowable wavelength. These properties make such solutions a promising setting for exploring the mechanisms behind extreme cosmic sources of gamma rays.

Enormous progress has been made on observing stellar magnetism in stars from the main sequence through to compact objects. Recent data have thrown into sharper relief the vexed question of the origin of stellar magnetic fields, which remains one of the main unanswered questions in astrophysics. In this chapter we review recent work in this area of research. In particular, we look at the fossil field hypothesis which links magnetism in compact stars to magnetism in main sequence and pre-main sequence stars and we consider why its feasibility has now been questioned particularly in the context of highly magnetic white dwarfs. We also review the fossil versus dynamo debate in the context of neutron stars and the roles played by key physical processes such as buoyancy, helicity, and superfluid turbulence,in the generation and stability of neutron star fields. Independent information on the internal magnetic field of neutron stars will come from future gravitational wave detections. Thus we maybe at the dawn of a new era of exciting discoveries in compact star magnetism driven by the opening of a new, non-electromagnetic observational window. We also review recent advances in the theory and computation of magnetohydrodynamic turbulence as it applies to stellar magnetism and dynamo theory. These advances offer insight into the action of stellar dynamos as well as processes whichcontrol the diffusive magnetic flux transport in stars.

The free decay of non-helical relativistic magnetohydrodynamic turbulence is studied numerically, and found to exhibit cascading of magnetic energy toward large scales. Evolution of the magnetic energy spectrum $P_M(k,t)$ is self-similar in time and well modeled by a broken power law with sub-inertial and inertial range indices very close to $7/2$ and $-2$ respectively. The magnetic coherence scale is found to grow in time as $t^{2/5}$, much too slow to account for optical polarization of gamma-ray burst afterglow emission if magnetic energy is to be supplied only at microphysical length scales. No bursty or explosive energy loss is observed in relativistic MHD turbulence having modest magnetization, which constrains magnetic reconnection models for rapid time variability of GRB prompt emission, blazars and the Crab nebula.

The simultaneous detection of electromagnetic and gravitational wave emission from merging neutron star binaries would aid greatly in their discovery and interpretation. By studying turbulent amplification of magnetic fields in local high-resolution simulations of neutron star merger conditions, we demonstrate that magnetar-level (~10^16) G fields are present throughout the merger duration. We find that the small-scale turbulent dynamo converts 60% of the randomized kinetic energy into magnetic fields on a merger time scale. Since turbulent magnetic energy dissipates through reconnection events which accelerate relativistic electrons, turbulence may facilitate the conversion of orbital kinetic energy into radiation. If 10^-4 of the ~ 10^53 erg of orbital kinetic available gets processed through reconnection, and creates radiation in the 15-150 keV band, then the fluence at 200 Mpc would be 10^-7 erg/cm^2, potentially rendering most merging neutron stars in the advanced LIGO and Virgo detection volumes detectable by Swift BAT.

High-resolution numerical simulations are utilized to examine isotropic turbulence in a compressible fluid when long wavelength velocity fluctuations approach light speed. Spectral analysis reveals an inertial sub-range of relativistic motions with a broadly 5/3 index. The use of generalized Lorentz-covariant structure functions based on the four-velocity is proposed. These structure functions extend the She-Leveque model for intermittency into the relativistic regime.

Models for GRB outflows invoke turbulence in relativistically hot magnetized fluids. In order to investigate these conditions we have performed high-resolution three-dimensional numerical simulations of relativistic magneto-hydrodynamical (RMHD) turbulence. We find that magnetic energy is amplified to several percent of the total energy density by turbulent twisting and folding of magnetic field lines. Values of epsilon_B near 1% are thus naturally expected. We study the dependence of saturated magnetic field energy fraction as a function of Mach number and relativistic temperature. We then present power spectra of the turbulent kinetic and magnetic energies. We also present solenoidal (curl-like) and dilatational (divergence-like) power spectra of kinetic energy. We propose that relativistic effects introduce novel couplings between these spectral components. The case we explore in most detail is for equal amounts of thermal and rest mass energy, corresponding to conditions after collisions of shells with relative Lorentz factors of several. These conditions are relevant in models for internal shocks, for the late afterglow phase, for cocoon material along the edge of a relativistic jet as it propagates through a star, as well neutron stars merging with each other and with black hole companions. We find that relativistic turbulence decays extremely quickly, on a sound crossing time of an eddy. Models invoking sustained relativistic turbulence to explain variability in GRB prompt emission are thus strongly disfavored unless a persistant driving of the turbulence is maintained for the duration of the prompt emission.

A wide variety of astrophysical phenomena involve the flow of turbulent magnetized gas with relativistic velocity or energy density. Examples include gamma-ray bursts, active galactic nuclei, pulsars, magnetars, micro-quasars, merging neutron stars, X-ray binaries, some supernovae, and the early universe. In order to elucidate the basic properties of the relativistic magnetohydrodynamical (RMHD) turbulence present in these systems, we present results from numerical simulations of fully developed driven turbulence in a relativistically warm, weakly magnetized and mildly compressible ideal fluid. We have evolved the RMHD equations for many dynamical times on a uniform grid with 1024^3 zones using a high order Godunov code. We observe the growth of magnetic energy from a seed field through saturation at about 1% of the total fluid energy. We compute the power spectrum of velocity and density-weighted velocity and conclude that the inertial scaling is consistent with a slope of -5/3. We compute the longitudinal and transverse velocity structure functions of order p up to 11, and discuss their possible deviation from the expected scaling for non-relativistic media. We also compute the scale-dependent distortion of coherent velocity structures with respect to the local magnetic field, finding a weaker scale dependence than is expected for incompressible non-relativistic flows with a strong mean field.