We report on an analytical and numerical study of the dynamics of a three-dimensional array of identical magnetic flux tubes in the reduced-magnetohydrodynamic description of the plasma. We propose that the long-time evolution of this system is dictated by flux-tube mergers and that such mergers are dynamically constrained by the conservation of the pertinent (ideal) invariants, {\it viz.} the magnetic potential and axial fluxes of each tube. We also propose that in the direction perpendicular to the merging plane, flux tubes evolve in critically-balanced fashion. These notions allow us to construct an analytical model for how quantities such as the magnetic energy and the energy-containing scale evolve as functions of time. Of particular importance is the conclusion that, like its two-dimensional counterpart, this system exhibits an inverse transfer of magnetic energy that terminates only at the system scale. We perform direct numerical simulations that confirm these predictions and reveal other interesting aspects of the evolution of the system. We find, for example, that the early time evolution is characterized by a sharp decay of the initial magnetic energy, which we attribute to the ubiquitous formation of current sheets. We also show that a quantitatively similar inverse transfer of magnetic energy is observed when the initial condition is a random, small-scale magnetic seed field.

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

Nonthermal relativistic plasmas are ubiquitous in astrophysical systems like pulsar wind nebulae and active galactic nuclei, as inferred from their emission spectra. The underlying nonthermal particle acceleration (NTPA) processes have traditionally been modeled with a Fokker-Planck (FP) diffusion-advection equation in momentum space. In this paper, we directly test the FP framework in ab-initio kinetic simulations of driven magnetized turbulence in relativistic pair plasma. By statistically analyzing the motion of tracked particles, we demonstrate the diffusive nature of NTPA and measure the FP energy diffusion ($D$) and advection ($A$) coefficients as functions of particle energy $\gamma m_e c^2$. We find that $D(\gamma)$ scales as $\gamma^2$ in the high-energy nonthermal tail, in line with 2nd-order Fermi acceleration theory, but has a much weaker scaling at lower energies. We also find that $A$ is not negligible and reduces NTPA by tending to pull particles towards the peak of the particle energy distribution. This study provides strong support for the FP picture of turbulent NTPA, thereby enhancing our understanding of space and astrophysical plasmas.

Rapid gamma-ray flares pose an astrophysical puzzle, requiring mechanisms both to accelerate energetic particles and to produce fast observed variability. These dual requirements may be satisfied by collisionless relativistic magnetic reconnection. On the one hand, relativistic reconnection can energize gamma-ray emitting electrons. On the other, as previous kinetic simulations have shown, the reconnection acceleration mechanism preferentially focuses high-energy particles -- and their emitted photons -- into beams, which may create rapid blips in flux as they cross a telescope's line of sight. Using a series of 2D pair-plasma particle-in-cell simulations, we explicitly demonstrate the critical role played by radiative cooling in mediating the observable signatures of this `kinetic beaming' effect. Only in our efficiently cooled simulations do we measure kinetic beaming beyond one light crossing time of the reconnection layer. We find a correlation between the cooling strength and the photon energy range across which persistent kinetic beaming occurs: stronger cooling coincides with a wider range of beamed photon energies. We also apply our results to rapid gamma-ray flares in flat-spectrum radio quasars, suggesting that a paradigm of radiatively efficient kinetic beaming constrains relevant emission models. In particular, beaming-produced variability may be more easily realized in two-zone (e.g. spine-sheath) setups, with Compton seed photons originating in the jet itself, rather than in one-zone external Compton scenarios.

High-energy astrophysical systems frequently contain collisionless relativistic plasmas that are heated by turbulent cascades and cooled by emission of radiation. Understanding the nature of this radiative turbulence is a frontier of extreme plasma astrophysics. In this paper, we use particle-in-cell simulations to study the effects of external inverse Compton radiation on turbulence driven in an optically thin, relativistic pair plasma. We focus on the statistical steady state (where injected energy is balanced by radiated energy) and perform a parameter scan spanning from low magnetization to high magnetization ($0.04 \lesssim \sigma \lesssim 11$). We demonstrate that the global particle energy distributions are quasi-thermal in all simulations, with only a modest population of nonthermal energetic particles (extending the tail by a factor of $\sim 2$). This indicates that nonthermal particle acceleration (observed in similar non-radiative simulations) is quenched by strong radiative cooling. The quasi-thermal energy distributions are well fit by analytic models in which stochastic particle acceleration (due to, e.g., second-order Fermi mechanism or gyroresonant interactions) is balanced by the radiation reaction force. Despite the efficient thermalization of the plasma, nonthermal energetic particles do make a conspicuous appearance in the anisotropy of the global momentum distribution as highly variable, intermittent beams (for high magnetization cases). The beamed high-energy particles are spatially coincident with intermittent current sheets, suggesting that localized magnetic reconnection may be a mechanism for kinetic beaming. This beaming phenomenon may explain rapid flares observed in various astrophysical systems (such as blazar jets, the Crab nebula, and Sagittarius A*).

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.

Many relativistic plasma environments in high-energy astrophysics, including pulsar wind nebulae, hot accretion flows onto black holes, relativistic jets in active galactic nuclei and gamma-ray bursts, and giant radio lobes, are naturally turbulent. The plasma in these environments is often so hot that synchrotron and inverse-Compton (IC) radiative cooling becomes important. In this paper we investigate the general thermodynamic and radiative properties (and hence the observational appearance) of an optically thin relativistically hot plasma stirred by driven magnetohydrodynamic (MHD) turbulence and cooled by radiation. We find that if the system reaches a statistical equilibrium where turbulent heating is balanced by radiative cooling, the effective electron temperature tends to attain a universal value $\theta = kT_e/m_e c^2 \sim 1/\sqrt{\tau_T}$, where $\tau_T=n_e\sigma_T L \ll 1$ is the system's Thomson optical depth, essentially independent of the strength of turbulent driving or magnetic field. This is because both MHD turbulent dissipation and synchrotron cooling are proportional to the magnetic energy density. We also find that synchrotron self-Compton (SSC) cooling and perhaps a few higher-order IC components are automatically comparable to synchrotron in this regime. The overall broadband radiation spectrum then consists of several distinct components (synchrotron, SSC, etc.), well separated in photon energy (by a factor $\sim \tau_T^{-1}$) and roughly equal in power. The number of IC peaks is checked by Klein-Nishina effects and depends logarithmically on $\tau_T$ and the magnetic field. We also examine the limitations due to synchrotron self-absorption, explore applications to Crab PWN and blazar jets, and discuss links to radiative magnetic reconnection.

In this Letter we propose that coherent radio emission of Crab, other young energetic pulsars, and millisecond pulsars is produced in the magnetospheric current sheet beyond the light cylinder. We carry out global and local two-dimensional kinetic plasma simulations of reconnection to illustrate the coherent emission mechanism. Reconnection in the current sheet beyond the light cylinder proceeds in the very efficient plasmoid-dominated regime, and current layer gets fragmented into a dynamic chain of plasmoids which undergo successive coalescence. Mergers of sufficiently large plasmoids produce secondary perpendicular current sheets, which are also plasmoid-unstable. Collisions of plasmoids with each other and with the upstream magnetic field eject fast-magnetosonic waves, which propagate upstream across the background field and successfully escape from the plasma as electromagnetic waves that fall in the radio band. This model successfully explains many important features of the observed radio emission from Crab and other pulsars with high magnetic field at the light cylinder: phase coincidence with the high-energy emission, nano-second duration (nanoshots), and extreme instantaneous brightness of individual pulses.

Electron and ion energization (i.e., heating and nonthermal acceleration) is a fundamental, but poorly understood, outcome of plasma turbulence. In this work, we present new results on this topic from particle-in-cell simulations of driven turbulence in collisionless, relativistic electron-ion plasma. We focus on temperatures such that ions (protons) are sub-relativistic and electrons are ultra-relativistic, a regime relevant for high-energy astrophysical systems such as hot accretion flows onto black holes. We find that ions tend to be preferentially heated, gaining up to an order of magnitude more energy than electrons, and propose a simple empirical formula to describe the electron-ion energy partition as a function of the ratio of electron-to-ion gyroradii (which in turn is a function of initial temperatures and plasma beta). We also find that while efficient nonthermal particle acceleration occurs for both species in the ultra-relativistic regime, nonthermal electron populations are diminished with decreasing temperature whereas nonthermal ion populations are essentially unchanged. These results have implications for modeling and interpreting observations of hot accretion flows.

Magnetic energy transfer from small to large scales due to successive magnetic island coalescence is investigated. A solvable analytical model is introduced and shown to correctly capture the evolution of the main quantities of interest, as borne out by numerical simulations. Magnetic reconnection is identified as the key mechanism enabling the inverse transfer, and setting its properties: magnetic energy decays as $\tilde t^{-1}$, where $\tilde t$ is time normalized to the (appropriately defined) reconnection timescale; and the correlation length of the field grows as $\tilde t^{1/2}$. The magnetic energy spectrum is self-similar, and evolves as $\propto \tilde t ^{-3/2}k^{-2}$, where the $k$-dependence is imparted by the formation of thin current sheets.

Particle-in-cell (PIC) simulations have shown that relativistic collisionless magnetic reconnection drives nonthermal particle acceleration (NTPA), potentially explaining high-energy (X-ray/$\gamma$-ray) synchrotron and/or inverse Compton (IC) radiation observed from various astrophysical sources. The radiation back-reaction force on radiating particles has been neglected in most of these simulations, even though radiative cooling considerably alters particle dynamics in many astrophysical environments where reconnection may be important. We present a radiative PIC study examining the effects of external IC cooling on the basic dynamics, NTPA, and radiative signatures of relativistic reconnection in pair plasmas. We find that, while the reconnection rate and overall dynamics are basically unchanged, IC cooling significantly influences NTPA: the particle spectra still show a hard power law (index $\geq -2$) as in nonradiative reconnection, but transition to a steeper power law that extends to a cooling-dependent cutoff. The steep power-law index fluctuates in time between roughly $-$3 and $-$5. The time-integrated photon spectra display corresponding power laws with indices $\approx -0.5$ and $\approx -1.1$, similar to those observed in hard X-ray spectra of accreting black holes.

Strong magnetic fields in magnetospheres of neutron stars (especially magnetars) and other astrophysical objects may release their energy in violent, intense episodes of magnetic reconnection. While reconnection has been studied extensively, the extreme field strength near neutron stars introduces new effects: synchrotron cooling and electron-positron pair production. Using massively parallel particle-in-cell simulations that self-consistently incorporate these new quantum-electrodynamic effects, we investigate relativistic magnetic reconnection in the strong-field regime. We show that reconnection in this regime can efficiently convert magnetic energy to X-ray and gamma-ray radiation and thus power bright high-energy astrophysical flares. Rapid radiative cooling causes strong plasma and magnetic field compression in compact plasmoids. In the most extreme cases, the field can approach the critical quantum limit, leading to copious pair production.

Black-hole driven relativistic astrophysical jets, such as blazars and gamma-ray bursts (GRBs), are powerful sources of electromagnetic radiation. Their emission is powered by some energy dissipation and particle acceleration mechanism operating in-situ at large distances from the black hole. We propose that the formation of the dissipative structures in the jet is controlled by the time variability of the accretion disc. We argue that the open magnetic field lines through the black hole, which drive a strongly magnetized jet, may have their polarity reversing over time-scales related to the growth of the magneto-rotational dynamo in the disc. Consequently, the jet is characterized by an alternating toroidal field polarity along its propagation axis, i.e., it is a "striped jet". Magnetic reconnection in the current sheets that form between the stripes dissipates the alternating-field energy and thus powers further jet acceleration. Here, we consider a jet with a broad distribution of stripe widths $l>l_{\rm min}$, above a dominant scale $l_{\rm min}$. We find that the bulk acceleration of the jet, driven by the annihilation of the stripes, is very gradual. The dissipation rate peaks at a distance $z_{\rm peak}\sim 10^6 R_g\, (\Gamma_\infty/30)^2\, (l_{\rm min}/1000R_g)$, where $R_{g}$ is the black-hole's gravititional radius and $\Gamma_{\infty}$ the jet's asymptotic Lorentz factor, and exhibits a very broad plateau extending by $\sim 4-5$ orders of magnitude in distance. The prolonged energy dissipation accounts for the flat-to-inverted long-wavelength spectra commonly observed in jets. The model can also account for the broad range of flaring timescales of blazars and the fact that their bulk acceleration appears to continue out to ~100 pc scales. In GRB jets, the model predicts comparable power for the photopsheric and Thomson-thin emission components.

We apply collisionless particle-in-cell simulations of relativistic pair plasmas to explore whether driven turbulence is a viable high-energy astrophysical particle accelerator. We characterize nonthermal particle distributions for varying system sizes up to $L/2\pi\rho_{e0} = 163$, where $L/2\pi$ is the driving scale and $\rho_{e0}$ is the initial characteristic Larmor radius. We show that turbulent particle acceleration produces power-law energy distributions that, when compared at a fixed number of large-scale dynamical times, slowly steepen with increasing system size. We demonstrate, however, that convergence is obtained by comparing the distributions at different times that increase with system size (approximately logarithmically). We suggest that the system-size dependence arises from the time required for particles to reach the highest accessible energies via Fermi acceleration. The converged power-law index of the energy distribution, $\alpha \approx 3.0$ for magnetization $\sigma = 3/8$, makes turbulence a possible explanation for nonthermal spectra observed in systems such as the Crab nebula.

Sweet-Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfv\'enic incompressible flows which narrow the sheet. Standard tearing theory (Furth et al. 1963; Coppi et al. 1976; Rutherford 1973) is not immediately applicable to such forming sheets for two reasons: first, because the flow introduces terms not present in the standard calculation; second, because the changing equilibrium introduces time dependence to terms which are constant in the standard calculation, complicating the formulation of an eigenvalue problem. This paper adapts standard tearing mode analysis to confront these challenges. In an initial phase when any perturbations are primarily governed by ideal MHD, a coordinate transformation reveals that the flow compresses and stretches perturbations. A multiple scale formulation describes how linear tearing mode theory (Furth et al. 1963; Coppi et al. 1976) can be applied to an equilibrium changing under flow, showing that the flow affects the separable exponential growth only implicitly, by making the standard scalings time-dependent. In the nonlinear Rutherford stage, the coordinate transformation shows that standard theory can be adapted by adding to the stationary rates time dependence and an additional term due to the strengthening equilibrium magnetic field. Overall, this understanding supports the use of flow-free scalings with slight modifications to study tearing in a forming sheet.

Magnetic reconnection in relativistic collisionless plasmas can accelerate particles and power high-energy emission in various astrophysical systems. Whereas most previous studies focused on relativistic reconnection in pair plasmas, less attention has been paid to electron-ion plasma reconnection, expected in black hole accretion flows and relativistic jets. We report a comprehensive particle-in-cell numerical investigation of reconnection in an electron-ion plasma, spanning a wide range of ambient ion magnetizations $\sigma_i$, from the semirelativistic regime (ultrarelativistic electrons but nonrelativistic ions, 0.001<<$\sigma_i$<<1) to the fully relativistic regime (both species are ultrarelativistic, $\sigma_i$>>1). We investigate how the reconnection rate, electron and ion plasma flows, electric and magnetic field structures, electron/ion energy partitioning, and nonthermal particle acceleration depend on $\sigma_i$. Our key findings are: (1) the reconnection rate is about 0.1 of the Alfvenic rate across all regimes; (2) electrons can form concentrated moderately relativistic outflows even in the semirelativistic, small-$\sigma_i$ regime; (3) while the released magnetic energy is partitioned equally between electrons and ions in the ultrarelativistic limit, the electron energy fraction declines gradually with decreased $\sigma_i$ and asymptotes to about 0.25 in the semirelativistic regime; (4) reconnection leads to efficient nonthermal electron acceleration with a $\sigma_i$-dependent power-law index, $p(\sigma_i) \simeq $const$+0.7 {\sigma_i}^{-1/2}$. These findings are important for understanding black hole systems and lend support to semirelativistic reconnection models for powering nonthermal emission in blazar jets, offering a natural explanation for the spectral indices observed in these systems.

We describe results from particle-in-cell simulations of driven turbulence in collisionless, magnetized, relativistic pair plasma. This physical regime provides a simple setting for investigating the basic properties of kinetic turbulence and is relevant for high-energy astrophysical systems such as pulsar wind nebulae and astrophysical jets. In this paper, we investigate the statistics of turbulent fluctuations in simulations on lattices of up to $1024^3$ cells and containing up to $2 \times 10^{11}$ particles. Due to the absence of a cooling mechanism in our simulations, turbulent energy dissipation reduces the magnetization parameter to order unity within a few dynamical times, causing turbulent motions to become sub-relativistic. In the developed stage, our results agree with predictions from magnetohydrodynamic turbulence phenomenology at inertial-range scales, including a power-law magnetic energy spectrum with index near $-5/3$, scale-dependent anisotropy of fluctuations described by critical balance, log-normal distributions for particle density and internal energy density (related by a $4/3$ adiabatic index, as predicted for an ultra-relativistic ideal gas), and the presence of intermittency. We also present possible signatures of a kinetic cascade by measuring power-law spectra for the magnetic, electric, and density fluctuations at sub-Larmor scales.

As a fundamental process converting magnetic to plasma energy in high-energy astrophysical plasmas, relativistic magnetic reconnection is a leading explanation for the acceleration of particles to the ultrarelativistic energies necessary to power nonthermal emission (especially X-rays and gamma-rays) in pulsar magnetospheres and pulsar wind nebulae, coronae and jets of accreting black holes, and gamma-ray bursts. An important objective of plasma astrophysics is therefore the characterization of nonthermal particle acceleration (NTPA) effected by reconnection. Reconnection-powered NTPA has been demonstrated over a wide range of physical conditions using large two-dimensional (2D) kinetic simulations. However, its robustness in realistic 3D reconnection -- in particular, whether the 3D relativistic drift-kink instability (RDKI) disrupts NTPA -- has not been systematically investigated, although pioneering 3D simulations have observed NTPA in isolated cases. Here we present the first comprehensive study of NTPA in 3D relativistic reconnection in collisionless electron-positron plasmas, characterizing NTPA as the strength of 3D effects is varied systematically via the length in the third dimension and the strength of the guide magnetic field. We find that, while the RDKI prominently perturbs 3D reconnecting current sheets, it does not suppress particle acceleration, even for zero guide field; fully 3D reconnection robustly and efficiently produces nonthermal power-law particle spectra closely resembling those obtained in 2D. This finding provides strong support for reconnection as the key mechanism powering high-energy flares in various astrophysical systems. We also show that strong guide fields significantly inhibit NTPA, slowing reconnection and limiting the energy available for plasma energization, yielding steeper and shorter power-law spectra.

We present results from particle-in-cell simulations of driven turbulence in magnetized, collisionless, and relativistic pair plasma. We find that fluctuations are consistent with the classical $k_\perp^{-5/3}$ magnetic energy spectrum at fluid scales and a steeper $k_\perp^{-4}$ spectrum at sub-Larmor scales, where $k_\perp$ is the wavevector perpendicular to the mean field. We demonstrate the development of a non-thermal, power-law particle energy distribution, $f(E) \sim E^{-\alpha}$, with index $\alpha$ that decreases with increasing magnetization and increases with increasing system size (relative to the characteristic Larmor radius). Our simulations indicate that turbulence can be a viable source of energetic particles in high-energy astrophysical systems, such as pulsar wind nebulae, if scalings asymptotically become insensitive to the system size.

We carry out a systematic study of the dispersion relation for linear electrostatic waves in an arbitrarily degenerate quantum electron plasma. We solve for the complex frequency spectrum for arbitrary values of wavenumber $k$ and level of degeneracy $\mu$. Our finding is that for large $k$ and high $\mu$ the real part of the frequency $\omega_{r}$ grows linearly with $k$ and scales with $\mu$ only because of the scaling of the Fermi energy. In this regime the relative Landau damping rate $\gamma/\omega_{r}$ becomes independent of $k$ and varies inversly with $\mu$. Thus, damping is weak but finite at moderate levels of degeneracy for short wavelengths.

The recent realization that Sweet-Parker current sheets are violently unstable to the secondary tearing (plasmoid) instability implies that such current sheets cannot occur in real systems. This suggests that, in order to understand the onset of magnetic reconnection, one needs to consider the growth of the tearing instability in a current layer as it is being formed. Such an analysis is performed here in the context of nonlinear resistive MHD for a generic time-dependent equilibrium representing a gradually forming current sheet. It is shown that two onset regimes, single-island and multi-island, are possible, depending on the rate of current sheet formation. A simple model is used to compute the criterion for transition between these two regimes, as well as the reconnection onset time and the current sheet parameters at that moment. For typical solar corona parameters this model yields results consistent with observations.

Turbulence is ubiquitous in plasmas, leading to rich dynamics characterized by irregularity, irreversibility, energy fluctuations across many scales, and energy transfer across many scales. Another fundamental and generic feature of turbulence, although sometimes overlooked, is the inhomogeneous dissipation of energy in space and in time. This is a consequence of intermittency, the scale-dependent inhomogeneity of dynamics caused by fluctuations in the turbulent cascade. Intermittency causes turbulent plasmas to self-organize into coherent dissipative structures, which may govern heating, diffusion, particle acceleration, and radiation emissions. In this paper, we present recent progress on understanding intermittency in incompressible magnetohydrodynamic turbulence with a strong guide field. We focus on the statistical analysis of intermittent dissipative structures, which occupy a small fraction of the volume but arguably account for the majority of energy dissipation. We show that, in our numerical simulations, intermittent structures in the current density, vorticity, and Els\"{a}sser vorticities all have nearly identical statistical properties. We propose phenomenological explanations for the scalings based on general considerations of Els\"{a}sser vorticity structures. Finally, we examine the broader implications of intermittency for astrophysical systems.

Using two-dimensional particle-in-cell simulations, we characterize the energy spectra of particles accelerated by relativistic magnetic reconnection (without guide field) in collisionless electron-positron plasmas, for a wide range of upstream magnetizations $\sigma$ and system sizes $L$. The particle spectra are well-represented by a power law $\gamma^{-\alpha}$, with a combination of exponential and super-exponential high-energy cutoffs, proportional to $\sigma$ and $L$, respectively. For large $L$ and $\sigma$, the power-law index $\alpha$ approaches about 1.2.

We investigate the distribution of particle acceleration sites, independently of the actual acceleration mechanism, during plasmoid-dominated, relativistic collisionless magnetic reconnection by analyzing the results of a particle-in-cell numerical simulation. The simulation is initiated with Harris-type current layers in pair plasma with no guide magnetic field, negligible radiative losses, no initial perturbation, and using periodic boundary conditions. We find that the plasmoids develop a robust internal structure, with colder dense cores and hotter outer shells, that is recovered after each plasmoid merger on a dynamical time scale. We use spacetime diagrams of the reconnection layers to probe the evolution of plasmoids, and in this context we investigate the individual particle histories for a representative sample of energetic electrons. We distinguish three classes of particle acceleration sites associated with (1) magnetic X-points, (2) regions between merging plasmoids, and (3) the trailing edges of accelerating plasmoids. We evaluate the contribution of each class of acceleration sites to the final energy distribution of energetic electrons -- magnetic X-points dominate at moderate energies, and the regions between merging plasmoids dominate at higher energies. We also identify the dominant acceleration scenarios, in order of decreasing importance -- (1) single acceleration between merging plasmoids, (2) single acceleration at a magnetic X-point, and (3) acceleration at a magnetic X-point followed by acceleration in a plasmoid. Particle acceleration is absent only in the vicinity of stationary plasmoids. The effect of magnetic mirrors due to plasmoid contraction does not appear to be significant in relativistic reconnection.

I review a new rapidly growing area of high-energy plasma astrophysics --- radiative magnetic reconnection, i.e., a reconnection regime where radiation reaction influences reconnection dynamics, energetics, and nonthermal particle acceleration. This influence be may be manifested via a number of astrophysically important radiative effects, such as radiation-reaction limits on particle acceleration, radiative cooling, radiative resistivity, braking of reconnection outflows by radiation drag, radiation pressure, viscosity, and even pair creation at highest energy densities. Self-consistent inclusion of these effects in magnetic reconnection theory and modeling calls for serious modifications to our overall theoretical approach to the problem. In addition, prompt reconnection-powered radiation often represents our only observational diagnostic tool for studying remote astrophysical systems; this underscores the importance of developing predictive modeling capabilities to connect the underlying physical conditions in a reconnecting system to observable radiative signatures. This Chapter gives an overview of recent theoretical progress in developing basic physical understanding of radiative reconnection, with a special emphasis on astrophysically important radiation mechanisms like synchrotron, curvature, and inverse-Compton. It also offers a broad review of key high-energy astrophysical applications of radiative reconnection, such as: pulsar wind nebulae and magnetospheres, accreting black-hole coronae and jets in XRBs and AGN, magnetospheres of magnetars, and Gamma-Ray Bursts. Finally, this Chapter discusses the most critical open questions and outlines the directions for future research in this exciting new frontier of plasma astrophysics.

(abridged) Magnetic reconnection is the topological reconfiguration of the magnetic field in a plasma, accompanied by the violent release of energy and particle acceleration. Reconnection is as ubiquitous as plasmas themselves, with solar flares perhaps the most popular example. Over the last few years, the theoretical understanding of magnetic reconnection in large-scale fluid systems has undergone a major paradigm shift. The steady-state model of reconnection described by the famous Sweet-Parker (SP) theory, which dominated the field for ~50 years, has been replaced with an essentially time-dependent, bursty picture of the reconnection layer, dominated by the continuous formation and ejection of multiple secondary islands (plasmoids). Whereas in the SP model reconnection was predicted to be slow, a major implication of this new paradigm is that reconnection in fluid systems is fast (i.e., independent of the Lundquist number), provided that the system is large enough. This conceptual shift hinges on the realization that SP-like current layers are violently unstable to the plasmoid instability - implying, therefore, that such current sheets are super-critically unstable and thus can never form in the first place. This suggests that the formation of a current sheet and the subsequent reconnection process cannot be decoupled, as is commonly assumed. This paper provides an introductory-level overview of the recent developments in reconnection theory and simulations that led to this essentially new framework. We briefly discuss the role played by the plasmoid instability in selected applications, and describe some of the outstanding challenges that remain at the frontier of this subject. Amongst these are the analytical and numerical extension of the plasmoid instability to (i) 3D and (ii) non-MHD regimes. New results are reported in both cases.

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

Energy dissipation is highly intermittent in turbulent plasmas, being localized in coherent structures such as current sheets. The statistical analysis of spatial dissipative structures is an effective approach to studying turbulence. In this paper, we generalize this methodology to investigate four-dimensional spatiotemporal structures, i.e., dissipative processes representing sets of interacting coherent structures, which correspond to flares in astrophysical systems. We develop methods for identifying and characterizing these processes, and then perform a statistical analysis of dissipative processes in numerical simulations of driven magnetohydrodynamic turbulence. We find that processes are often highly complex, long-lived, and weakly asymmetric in time. They exhibit robust power-law probability distributions and scaling relations, including a distribution of dissipated energy with power-law index near -1.75, indicating that intense dissipative events dominate the overall energy dissipation. We compare our results with the previously observed statistical properties of solar flares.

Energy dissipation in magnetohydrodynamic (MHD) turbulence is known to be highly intermittent in space, being concentrated in sheet-like coherent structures. Much less is known about intermittency in time, another fundamental aspect of turbulence which has great importance for observations of solar flares and other space/astrophysical phenomena. In this Letter, we investigate the temporal intermittency of energy dissipation in numerical simulations of MHD turbulence. We consider four-dimensional spatiotemporal structures, "flare events", responsible for a large fraction of the energy dissipation. We find that although the flare events are often highly complex, they exhibit robust power-law distributions and scaling relations. We find that the probability distribution of dissipated energy has a power law index close to -1.75, similar to observations of solar flares, indicating that intense dissipative events dominate the heating of the system. We also discuss the temporal asymmetry of flare events as a signature of the turbulent cascade.

Certain classes of astrophysical objects, namely magnetars and central engines of supernovae and gamma-ray bursts (GRBs), are characterized by extreme physical conditions not encountered elsewhere in the Universe. In particular, they possess magnetic fields that exceed the critical quantum field of 44 teragauss. Figuring out how these complex ultra-magnetized systems work requires understanding various plasma processes, both small-scale kinetic and large-scale magnetohydrodynamic (MHD). However, an ultra-strong magnetic field modifies the underlying physics to such an extent that many relevant plasma-physical problems call for building QED-based relativistic quantum plasma physics. In this review, after describing the extreme astrophysical systems of interest and identifying the key relevant plasma-physical problems, we survey the recent progress in the development of such a theory. We discuss how a super-critical field modifies the properties of vacuum and matter and outline the basic theoretical framework for describing both non-relativistic and relativistic quantum plasmas. We then turn to astrophysical applications of relativistic QED plasma physics relevant to magnetar magnetospheres and central engines of supernovae and long GRBs. Specifically, we discuss propagation of light through a magnetar magnetosphere; large-scale MHD processes driving magnetar activity and GRB jet launching and propagation; energy-transport processes governing the thermodynamics of extreme plasma environments; micro-scale kinetic plasma processes important in the interaction of intense magnetospheric electric currents with a magnetar's surface; and magnetic reconnection of ultra-strong magnetic fields. Finally, we point out that future progress will require the development of numerical modeling capabilities.

The discovery of rapid synchrotron gamma-ray flares above 100 MeV from the Crab Nebula has attracted new interest in alternative particle acceleration mechanisms in pulsar wind nebulae. Diffuse shock-acceleration fails to explain the flares because particle acceleration and emission occur during a single or even sub-Larmor timescale. In this regime, the synchrotron energy losses induce a drag force on the particle motion that balances the electric acceleration and prevents the emission of synchrotron radiation above 160 MeV. Previous analytical studies and 2D particle-in-cell (PIC) simulations indicate that relativistic reconnection is a viable mechanism to circumvent the above difficulties. The reconnection electric field localized at X-points linearly accelerates particles with little radiative energy losses. In this paper, we check whether this mechanism survives in 3D, using a set of large PIC simulations with radiation reaction force and with a guide field. In agreement with earlier works, we find that the relativistic drift kink instability deforms and then disrupts the layer, resulting in significant plasma heating but few non-thermal particles. A moderate guide field stabilizes the layer and enables particle acceleration. We report that 3D magnetic reconnection can accelerate particles above the standard radiation reaction limit, although the effect is less pronounced than in 2D with no guide field. We confirm that the highest energy particles form compact bunches within magnetic flux ropes, and a beam tightly confined within the reconnection layer, which could result in the observed Crab flares when, by chance, the beam crosses our line of sight.

The Crab Nebula was formed after the collapse of a massive star about a thousand years ago, leaving behind a pulsar that inflates a bubble of ultra-relativistic electron-positron pairs permeated with magnetic field. The observation of brief but bright flares of energetic gamma rays suggests that pairs are accelerated to PeV energies within a few days; such rapid acceleration cannot be driven by shocks. Here, it is argued that the flares may be the smoking gun of magnetic dissipation in the Nebula. Using 2D and 3D particle-in-cell simulations, it is shown that the observations are consistent with relativistic magnetic reconnection, where pairs are subject to strong radiative cooling. The Crab flares may highlight the importance of relativistic magnetic reconnection in astrophysical sources.

It is generally accepted that astrophysical sources cannot emit synchrotron radiation above 160 MeV in their rest frame. This limit is given by the balance between the accelerating electric force and the radiation reaction force acting on the electrons. The discovery of synchrotron gamma-ray flares in the Crab Nebula, well above this limit, challenges this classical picture of particle acceleration. To overcome this limit, particles must accelerate in a region of high electric field and low magnetic field. This is possible only with a non-ideal magnetohydrodynamic process, like magnetic reconnection. We present the first numerical evidence of particle acceleration beyond the synchrotron burnoff limit, using a set of 2D particle-in-cell simulations of ultra-relativistic pair plasma reconnection. We use a new code, Zeltron, that includes self-consistently the radiation reaction force in the equation of motion of the particles. We demonstrate that the most energetic particles move back and forth across the reconnection layer, following relativistic Speiser orbits. These particles then radiate >160 MeV synchrotron radiation rapidly, within a fraction of a full gyration, after they exit the layer. Our analysis shows that the high-energy synchrotron flux is highly variable in time because of the strong anisotropy and inhomogeneity of the energetic particles. We discover a robust positive correlation between the flux and the cut-off energy of the emitted radiation, mimicking the effect of relativistic Doppler amplification. A strong guide field quenches the emission of >160 MeV synchrotron radiation. Our results are consistent with the observed properties of the Crab flares, supporting the reconnection scenario.

We develop a framework for studying the statistical properties of current sheets in numerical simulations of 3D magnetohydrodynamic (MHD) turbulence. We describe an algorithm that identifies current sheets in a simulation snapshot and then determines their geometrical properties (including length, width, and thickness) and intensities (peak current density and total energy dissipation rate). We then apply this procedure to simulations of reduced MHD turbulence and perform a statistical analysis on the obtained population of current sheets. We evaluate the role of reconnection by separately studying the populations of current sheets which contain magnetic X-points and those which do not. We find that the statistical properties of the two populations are different in general. We compare the scaling of these properties to phenomenological predictions obtained for the inertial range of MHD turbulence. Finally, we test whether the reconnecting current sheets are consistent with the Sweet-Parker model.

The magnetosphere of a rotating pulsar naturally develops a current sheet beyond the light cylinder (LC). Magnetic reconnection in this current sheet inevitably dissipates a nontrivial fraction of the pulsar spin-down power within a few LC radii. We develop a basic physical picture of reconnection in this environment and discuss its implications for the observed pulsed gamma-ray emission. We argue that reconnection proceeds in the plasmoid-dominated regime, via an hierarchical chain of multiple secondary islands/flux ropes. The inter-plasmoid reconnection layers are subject to strong synchrotron cooling, leading to significant plasma compression. Using the conditions of pressure balance across these current layers, the balance between the heating by magnetic energy dissipation and synchrotron cooling, and Ampere's law, we obtain simple estimates for key parameters of the layers --- temperature, density, and layer thickness. In the comoving frame of the relativistic pulsar wind just outside of the equatorial current sheet, these basic parameters are uniquely determined by the strength of the reconnecting upstream magnetic field. For the case of the Crab pulsar, we find them to be of order 10 GeV, $10^{13} cm^{-3}$, and 10 cm, respectively. After accounting for the bulk Doppler boosting due to the pulsar wind, the synchrotron and inverse-Compton emission from the reconnecting current sheet can explain the observed pulsed high-energy (GeV) and VHE (~100 GeV) radiation, respectively. Also, we suggest that the rapid relative motions of the secondary plasmoids in the hierarchical chain may contribute to the production of the pulsar radio emission.

A 2D linear theory of the instability of Sweet-Parker (SP) current sheets is developed in the framework of Reduced MHD. A local analysis is performed taking into account the dependence of a generic equilibrium profile on the outflow coordinate. The plasmoid instability [Loureiro et al, Phys. Plasmas {\bf 14}, 100703 (2007)] is recovered, i.e., current sheets are unstable to the formation of a large-wave-number chain of plasmoids ($k_{\rm max}\Lsheet \sim S^{3/8}$, where $k_{\rm max}$ is the wave-number of fastest growing mode, $S=\Lsheet V_A/\eta$ is the Lundquist number, $\Lsheet$ is the length of the sheet, $V_A$ is the Alfv\'en speed and $\eta$ is the plasma resistivity), which grows super-Alfv\'enically fast ($\gmax\tau_A\sim S^{1/4}$, where $\gmax$ is the maximum growth rate, and $\tau_A=\Lsheet/V_A$). For typical background profiles, the growth rate and the wave-number are found to {\it increase} in the outflow direction. This is due to the presence of another mode, the Kelvin-Helmholtz (KH) instability, which is triggered at the periphery of the layer, where the outflow velocity exceeds the Alfv\'en speed associated with the upstream magnetic field. The KH instability grows even faster than the plasmoid instability, $\gmax \tau_A \sim k_{\rm max} \Lsheet\sim S^{1/2}$. The effect of viscosity ($\nu$) on the plasmoid instability is also addressed. In the limit of large magnetic Prandtl numbers, $Pm=\nu/\eta$, it is found that $\gmax\sim S^{1/4}Pm^{-5/8}$ and $k_{\rm max} \Lsheet\sim S^{3/8}Pm^{-3/16}$, leading to the prediction that the critical Lundquist number for plasmoid instability in the $Pm\gg1$ regime is $\Scrit\sim 10^4Pm^{1/2}$. These results are verified via direct numerical simulation of the linearized equations, using a new, analytical 2D SP equilibrium solution.

We study theoretical implications of a rapid Very-High-Energy (VHE) flare detected by MAGIC in the Flat-Spectrum Radio Quasar PKS 1222+216. The minimum distance from the jet origin at which this flare could be produced is 0.5 pc. A moderate Doppler factor of the VHE source, D_{VHE} ~ 20, is allowed by all opacity constraints. The concurrent High-Energy (HE) emission observed by Fermi provides estimates of the total jet power and the jet magnetic field strength. Energetic constraints for the VHE flare are extremely tight: for an isotropic particle distribution they require a huge co-moving energy density in the emitting region and a very efficient radiative process. We disfavor hadronic processes due to their low radiative efficiency, as well as the synchrotron scenario recently proposed for the case of HE flares in the Crab Nebula, since the parameters needed to overcome the radiative losses are quite extreme. The VHE emission can be explained by the Synchrotron Self-Compton (SSC) mechanism for D_{VHE} ~ 20 or by the External Radiation Compton (ERC) mechanism involving the infrared radiation of the dusty torus for D_{VHE} ~ 50. After discussing several alternative scenarios, we propose that the extreme energy density constraint can be satisfied when the emission comes from highly anisotropic short-lived bunches of particles formed by the kinetic beaming mechanism in magnetic reconnection sites. By focusing the emitting particles into very narrow beams, this mechanism allows one to relax the causality constraint on the source size, decreasing the required energy density by 4 orders of magnitude.

We report on the first study of the angular distribution of energetic particles and radiation generated in relativistic collisionless electron-positron pair plasma reconnection, using two-dimensional particle-in-cell simulations. We discover a strong anisotropy of the particles accelerated by reconnection and the associated strong beaming of their radiation. The focusing of particles and radiation increases with their energy; in this sense, this "kinetic beaming" effect differs fundamentally from the relativistic Doppler beaming usually invoked in high-energy astrophysics, in which all photons are focused and boosted achromatically. We also present, for the first time, the modeling of the synchrotron emission as seen by an external observer during the reconnection process. The expected lightcurves comprise several bright symmetric sub-flares emitted by the energetic beam of particles sweeping across the line of sight intermittently, and exhibit super-fast time variability as short as about one tenth of the system light-crossing time. The concentration of the energetic particles into compact regions inside magnetic islands and particle anisotropy explain the rapid variability. This radiative signature of reconnection can account for the brightness and variability of the gamma-ray flares in the Crab Nebula and in blazars.

In this paper we consider two outstanding intertwined problems in modern high-energy astrophysics: (1) the vertical thermal structure of an optically thick accretion disk heated by the dissipation of magnetohydrodynamic (MHD) turbulence driven by the magneto-rotational instability (MRI), and (2) determining the fraction of the accretion power released in the corona above the disk. For simplicity, we consider a gas-pressure-dominated disk and assume a constant opacity. We argue that the local turbulent dissipation rate due to the disruption of MRI channel flows by secondary parasitic instabilities should be uniform across most of the disk, almost up to the disk photosphere. We then obtain a self-consistent analytical solution for the vertical thermal structure of the disk, governed by the balance between the heating by MRI turbulence and the cooling by radiative diffusion. Next, we argue that the coronal power fraction is determined by the competition between the Parker instability, viewed as a parasitic instability feeding off of MRI channel flows, and other parasitic instabilities. We show that the Parker instability inevitably becomes important near the disk surface, leading to a certain lower limit on the coronal power. While most of the analysis in this paper focuses on the case of a disk threaded by an externally imposed vertical magnetic field, we also discuss the zero-net-flux case, in which the magnetic field is produced by the MRI dynamo itself, and show that most of our arguments and conclusions should be valid in this case as well.

The gamma-ray space telescopes AGILE and Fermi detected short and bright synchrotron gamma-ray flares at photon energies above 100 MeV in the Crab Nebula. This discovery suggests that electron-positron pairs in the nebula are accelerated to PeV energies in a milliGauss magnetic field, which is difficult to explain with classical models of particle acceleration and pulsar wind nebulae. We investigate whether particle acceleration in a magnetic reconnection layer can account for the puzzling properties of the flares. We numerically integrate relativistic test-particle orbits in the vicinity of the layer, including the radiation reaction force, and using analytical expressions for the large-scale electromagnetic fields. As they get accelerated by the reconnection electric field, the particles are focused deep inside the current layer where the magnetic field is small. The electrons suffer less from synchrotron losses and are accelerated to extremely high energies. Population studies show that, at the end of the layer, the particle distribution piles up at the maximum energy given by the electric potential drop and is focused into a thin fan beam. Applying this model to the Crab Nebula, we find that the emerging synchrotron emission spectrum peaks above 100 MeV and is close to the spectral shape of a single electron. The flare inverse Compton emission is negligible and no detectable emission is expected at other wavelengths. This mechanism provides a plausible explanation for the gamma-ray flares in the Crab Nebula and could be at work in other astrophysical objects such as relativistic jets in active galactic nuclei.

Prompt gamma-ray burst (GRB) emission requires some mechanism to dissipate an ultrarelativistic jet. Internal shocks or some form of electromagnetic dissipation are candidate mechanisms. Any mechanism needs to answer basic questions, such as what is the origin of variability, what radius does dissipation occur at, and how does efficient prompt emission occur. These mechanisms also need to be consistent with how ultrarelativistic jets form and stay baryon pure despite turbulence and electromagnetic reconnection near the compact object and despite stellar entrainment within the collapsar model. We use the latest magnetohydrodynamical models of ultrarelativistic jets to explore some of these questions in the context of electromagnetic dissipation due to the slow collisional and fast collisionless reconnection mechanisms, as often associated with Sweet-Parker and Petschek reconnection, respectively. For a highly magnetized ultrarelativistic jet and typical collapsar parameters, we find that significant electromagnetic dissipation may be avoided until it proceeds catastrophically near the jet photosphere at large radii ($r\sim 10^{13}--10^{14}{\rm cm}$), by which the jet obtains a high Lorentz factor ($\gamma\sim 100--1000), has a luminosity of $L_j \sim 10^{50}--10^{51}\ergs$, has observer variability timescales of order 1s (ranging from 0.001-10s), achieves $\gamma\theta_j\sim 10--20 (for opening half-angle $\theta_j$) and so is able to produce jet breaks, and has comparable energy available for both prompt and afterglow emission. This reconnection switch mechanism allows for highly efficient conversion of electromagnetic energy into prompt emission and associates the observed prompt GRB pulse temporal structure with dissipation timescales of some number of reconnecting current sheets embedded in the jet.[abridged]

A numerical study of magnetic reconnection in the large-Lundquist-number ($S$), plasmoid-dominated regime is carried out for $S$ up to $10^7$. The theoretical model of Uzdensky {\it et al.} [Phys. Rev. Lett. {\bf 105}, 235002 (2010)] is confirmed and partially amended. The normalized reconnection rate is $\normEeff\sim 0.02$ independently of $S$ for $S\gg10^4$. The plasmoid flux ($\Psi$) and half-width ($w_x$) distribution functions scale as $f(\Psi)\sim \Psi^{-2}$ and $f(w_x)\sim w_x^{-2}$. The joint distribution of $\Psi$ and $w_x$ shows that plasmoids populate a triangular region $w_x\gtrsim\Psi/B_0$, where $B_0$ is the reconnecting field. It is argued that this feature is due to plasmoid coalescence. Macroscopic "monster" plasmoids with $w_x\sim 10$% of the system size are shown to emerge in just a few Alfv\'en times, independently of $S$, suggesting that large disruptive events are an inevitable feature of large-$S$ reconnection.

The recent discovery of day-long gamma-ray flares in the Crab Nebula, presumed to be synchrotron emission by PeV (10^{15} eV) electrons in milligauss magnetic fields, presents a strong challenge to particle acceleration models. The observed photon energies exceed the upper limit (~100 MeV) obtained by balancing the acceleration rate and synchrotron radiation losses under standard conditions where the electric field is smaller than the magnetic field. We argue that a linear electric accelerator, operating at magnetic reconnection sites, is able to circumvent this difficulty. Sufficiently energetic electrons have gyroradii so large that their motion is insensitive to small-scale turbulent structures in the reconnection layer and is controlled only by large-scale fields. We show that such particles are guided into the reconnection layer by the reversing magnetic field as they are accelerated by the reconnection electric field. As these electrons become confined within the current sheet, they experience a decreasing perpendicular magnetic field that may drop below the accelerating electric field. This enables them to reach higher energies before suffering radiation losses and hence to emit synchrotron radiation in excess of the 100 MeV limit, providing a natural resolution to the Crab gamma-ray flare paradox.

Magnetic reconnection is a basic plasma process of dramatic rearrangement of magnetic topology, often leading to a violent release of magnetic energy. It is important in magnetic fusion and in space and solar physics --- areas that have so far provided the context for most of reconnection research. Importantly, these environments consist just of electrons and ions and the dissipated energy always stays with the plasma. In contrast, in this paper I introduce a new direction of research, motivated by several important problems in high-energy astrophysics --- reconnection in high energy density (HED) radiative plasmas, where radiation pressure and radiative cooling become dominant factors in the pressure and energy balance. I identify the key processes distinguishing HED reconnection: special-relativistic effects; radiative effects (radiative cooling, radiation pressure, and Compton resistivity); and, at the most extreme end, QED effects, including pair creation. I then discuss the main astrophysical applications --- situations with magnetar-strength fields (exceeding the quantum critical field of about 4 x 10^13 G): giant SGR flares and magnetically-powered central engines and jets of GRBs. Here, magnetic energy density is so high that its dissipation heats the plasma to MeV temperatures. Electron-positron pairs are then copiously produced, making the reconnection layer highly collisional and dressing it in a thick pair coat that traps radiation. The pressure is dominated by radiation and pairs. Yet, radiation diffusion across the layer may be faster than the global Alfv\'en transit time; then, radiative cooling governs the thermodynamics and reconnection becomes a radiative transfer problem, greatly affected by the ultra-strong magnetic field. This overall picture is very different from our traditional picture of reconnection and thus represents a new frontier in reconnection research.

A conceptual model of resistive magnetic reconnection via a stochastic plasmoid chain is proposed. The global reconnection rate is shown to be independent of the Lundquist number. The distribution of fluxes in the plasmoids is shown to be an inverse square law. It is argued that there is a finite probability of emergence of abnormally large plasmoids, which can disrupt the chain (and may be responsible for observable large abrupt events in solar flares and sawtooth crashes). A criterion for the transition from magnetohydrodynamic to collisionless regime is provided.

We construct a numerical model of emission from minijets, localized flows driven by magnetic reconnection inside Poynting-flux-dominated jets proposed to explain the ultrafast variability of blazars. The geometrical structure of the model consists of two wedge-like regions of relativistically flowing gas, separated by a stationary shock. The dynamics is based on solutions of relativistic magnetic reconnection with a guide field from Lyubarsky (2005). Electron distributions in each region are chosen to the match the pressure and density of the local plasma. Synchrotron emission from both regions is used to calculate Compton scattering, Compton drag and photon-photon opacity effects, with exact treatment of anisotropy and the Klein-Nishina regime. Radiative effects on plasma are taken into account, including the dependence of pressure on electron radiative losses and adiabatic heating of the flow decelerating under Compton drag. The results are applied to the July 2006 flare in the BL Lac object PKS 2155-304, with the aim of matching TeV flux measurements by H.E.S.S. with models that satisfy the variability constraints, while keeping X-ray emission below simultaneous Chandra observations. We find that models of isolated minijets with a significant guide field overproduce X-ray emission, and that we must take into account the radiative interaction of oppositely-oriented minijets in order to achieve a high enough dominance by Comptonized TeV radiation. We argue that such interactions are likely to occur in a jet where there is substantial internal reconnection, producing a large number of misaligned minijets. Finally, we show that large jet magnetizations are indeed required to satisfy all observational constraints and that the effective Lorentz factor of the minijet plasma has to be larger than 50, in agreement with earlier one-zone estimates.

Magnetic reconnection, a fundamental plasma process associated with a rapid dissipation of magnetic energy, is believed to power many disruptive phenomena in laboratory plasma devices, the Earth magnetosphere, and the solar corona. Traditional reconnection research, geared towards these rather tenuous environments, has justifiably ignored the effects of radiation on the reconnection process. However, in many reconnecting systems in high-energy astrophysics (e.g., accretion-disk coronae, relativistic jets, magnetar flares) and, potentially, in powerful laser plasma and z-pinch experiments, the energy density is so high that radiation, in particular radiative cooling, may start to play an important role. This observation motivates the development of a theory of high-energy-density radiative magnetic reconnection. As a first step towards this goal, we present in this paper a simple Sweet--Parker-like theory of non-relativistic resistive-MHD reconnection with strong radiative cooling. First, we show how, in the absence of a guide magnetic field, intense cooling leads to a strong compression of the plasma in the reconnection layer, resulting in a higher reconnection rate. The compression ratio and the layer temperature are determined by the balance between ohmic heating and radiative cooling. The lower temperature in the radiatively-cooled layer leads to a higher Spitzer resistivity and hence to an extra enhancement of the reconnection rate. We then apply our general theory to several specific astrophysically important radiative processes (bremsstrahlung, cyclotron, and inverse-Compton) in the optically thin regime, for both the zero- and strong-guide-field cases. We derive specific expressions for key reconnection parameters, including the reconnection rate. We also discuss the limitations and conditions for applicability of our theory.

The jet of the radio galaxy M87 is misaligned, resulting in a Doppler factor delta~1 for emission of plasma moving parallel to the jet. This makes the observed fast TeV flares on timescales of t_v~5R_g/c harder to understand as emission from the jet. In previous work, we have proposed a jets-in-a-jet model for the ultra-fast TeV flares with t_v<<R_g/c seen in Mrk 501 and PKS 2155-304. Here, we show that about half of the minijets beam their emission outside the jet cone. Minijets emitting off the jet axis result in rapidly evolving TeV (and maybe lower energy) flares that can be observed in nearby radio galaxies. The TeV flaring from M87 fits well into this picture, if M87 is a misaligned blazar.

Two-dimensional numerical simulations of the effect of background turbulence on 2D resistive magnetic reconnection are presented. For sufficiently small values of the resistivity ($\eta$) and moderate values of the turbulent power ($\epsilon$), the reconnection rate is found to have a much weaker dependence on $\eta$ than the Sweet-Parker scaling of $\eta^{1/2}$ and is even consistent with an $\eta-$independent value. For a given value of $\eta$, the dependence of the reconnection rate on the turbulent power exhibits a critical threshold in $\epsilon$ above which the reconnection rate is significantly enhanced.

The fast TeV variability of the blazars Mrk 501 and PKS 2155--304 implies a compact emitting region that moves with a bulk Lorentz factor of Gamma_{em}~100 toward the observer. The Lorentz factor is clearly in excess of the jet Lorentz factors Gamma_j\simless 10 measured on sub-pc scales in these sources. We propose that the TeV emission originates from compact emitting regions that move relativistically {\it within} a jet of bulk Gamma_j~10. This can be physically realized in a Poynting flux-dominated jet. We show that if a large fraction of the luminosity of the jet is prone to magnetic dissipation through reconnection, then material outflowing from the reconnection regions can efficiently power the observed TeV flares through synchrotron-self-Compton emission. The model predicts simultaneous far UV/soft X-ray flares.

A detailed numerical study of magnetic reconnection in resistive MHD for very large, previously inaccessible, Lundquist numbers ($10^4\le S\le 10^8$) is reported. Large-aspect-ratio Sweet-Parker current sheets are shown to be unstable to super-Alfv\'enically fast formation of plasmoid (magnetic-island) chains. The plasmoid number scales as $S^{3/8}$ and the instability growth rate in the linear stage as $S^{1/4}$, in agreement with the theory by Loureiro et al. [Phys. Plasmas {\bf 14}, 100703 (2007)]. In the nonlinear regime, plasmoids continue to grow faster than they are ejected and completely disrupt the reconnection layer. These results suggest that high-Lundquist-number reconnection is inherently time-dependent and hence call for a substantial revision of the standard Sweet-Parker quasi-stationary picture for $S>10^4$.

This is a white paper submitted to the Stars and Stellar Evolution (SSE) Science Frontier Panel (SFP) of the NRC's Astronomy and Astrophysics 2010 Decadal Survey. The white paper is endorsed by the American Physical Society's (APS) Topical Group on Plasma Astrophysics (GPAP).

This is a white paper submitted to the Stars and Stellar Evolution (SSE) Science Frontier Panel (SFP) of the NRC's 2010 Astronomy and Astrophysics Decadal Survey. The white paper is endorsed by the NSF Physics Frontier Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas (CMSO).

A simple Sweet--Parker-like model for the electron current layer in resistive Hall magnetohydrodynamic (MHD) reconnection is presented, with the focus on the collisionless limit. The derivation readily recovers the main results obtained recently by Malyshkin [PRL, 101, 225001 (2008)] and others, but is much quicker and more physically transparent. In particular, it highlights the role of resistive drag in determining the electron outflow velocity. The principal limitations of any such approach are discussed.

We point out that a conventional construction placed upon observations of accreting black holes, in which their nonthermal X-ray spectra are produced by inverse comptonization in a coronal plasma, suggests that the plasma is marginally collisionless. Recent developments in plasma physics indicate that fast reconnection takes place only in collisionless plasmas. As has recently been suggested for the Sun's corona, such marginal states may result from a combination of energy balance and the requirements of fast magnetic reconnection.

This paper provides a concise summary of the current status of the research and future perspectives discussed in the Mini-Conference on Angular Momentum Transport in Laboratory and Nature. This Mini-conference, sponsored by the Topical Group on Plasma Astrophysics, was held as part of the American Physical Society's Division of Plasma Physics 2007 Annual Meeting (November 12--16, 2007). This Mini-conference covers a wide range of phenomena happening in fluids and plasmas, either in laboratory or in nature. The purpose of this paper is not to comprehensively review these phenomena, but to provide a starting point for interested readers to refer to related research in areas other than their own.

We present a physics-based statistical theory of a force-free magnetic field in the corona above a turbulent accretion disk. The field is represented by a statistical ensemble of loops tied to the disk. Each loop evolves under several physical processes: Keplerian shear, turbulent random walk of the disk footpoints, and reconnection with other loops. To build a statistical description, we introduce the distribution function of loops over their sizes and construct a kinetic equation that governs its evolution. This loop kinetic equation is formally analogous to Boltzmann's kinetic equation, with loop-loop reconnection described by a binary collision integral. A dimensionless parameter is introduced to scale the (unknown) overall rate of reconnection relative to Keplerian shear. After solving for the loop distribution function numerically, we calculate self-consistently the distribution of the mean magnetic pressure and dissipation rate with height, and the equilibrium shapes of loops of different sizes. We also compute the energy and torque associated with a given loop, as well as the total magnetic energy and torque in the corona. We explore the dependence of these quantities on the reconnection parameter and find that they can be greatly enhanced if reconnection between loops is suppressed.

I propose a new paradigm for solar coronal heating viewed as a self-regulating process keeping the plasma marginally collisionless. The mechanism is based on the coupling between two effects. First, coronal density controls the plasma collisionality and hence the transition between the slow collisional Sweet-Parker and the fast collisionless reconnection regimes. In turn, coronal energy release leads to chromospheric evaporation, increasing the density and thus inhibiting subsequent reconnection of the newly-reconnected loops. As a result, statistically, the density fluctuates around some critical level, comparable to that observed in the corona. In the long run, coronal heating can be represented by repeating cycles of fast reconnection events (nano-flares), evaporation episodes, and long periods of slow magnetic stress build-up and radiative cooling of the coronal plasma.

I propose that solar coronal heating is a self-regulating process that keeps the coronal plasma roughly marginally collisionless. The self-regulating mechanism is based on the interplay of two effects. First, plasma density controls coronal energy release via the transition between the slow collisional Sweet-Parker regime and the fast collisionless reconnection regime. This transition takes place when the Sweet--Parker layer becomes thinner than the characteristic collisionless reconnection scale. I present a simple criterion for this transition in terms of the upstream plasma density (n_e), the reconnecting (B_0) and guide (B_z) magnetic field components, and the global length (L) of the reconnection layer: L < 6.10^9 cm [n_e/(10^{10}/cm^3)]^(-3) (B_0/30G)^4 (B_0/B_z)^2. Next, coronal energy release by reconnection raises the ambient plasma density via chromospheric evaporation and this, in turn, temporarily inhibits subsequent reconnection involving the newly-reconnected loops. Over time, however, radiative cooling gradually lowers the density again below the critical value and fast reconnection again becomes possible. As a result, the density is highly inhomogeneous and intermittent but, statistically, does not deviate strongly from the critical value which is comparable with the observed coronal density. Thus, in the long run, the coronal heating process can be represented by repeating cycles that consist of fast reconnection events (i.e., nanoflares), followed by rapid evaporation episodes, followed by relatively long periods (1-hour) during which magnetic stresses build up and simultaneously the plasma cools down and precipitates.

It has been suggested that magnetic fields play a dynamically-important role in core-collapse explosions of massive stars. In particular, they may be important in the collapsar scenario for gamma-ray bursts (GRB), where the central engine is a hyper-accreting black hole or a millisecond magnetar. The present paper is focussed on the magnetar scenario, with a specific emphasis on the interaction of the magnetar magnetosphere with the infalling stellar envelope. First, the ``Pulsar-in-a-Cavity'' problem is introduced as a paradigm for a magnetar inside a collapsing star. The basic set-up of this fundamental plasma-physics problem is described, outlining its main features, and simple estimates are derived for the evolution of the magnetic field. In the context of a collapsing star, it is proposed that, at first, the ram pressure of the infalling plasma acts to confine the magnetosphere, enabling a gradual build-up of the magnetic pressure. At some point, the growing magnetic pressure overtakes the (decreasing) ram pressure of the gas, resulting in a magnetically-driven explosion. The explosion should be highly anisotropic, as the hoop-stress of the toroidal field, confined by the surrounding stellar matter, collimates the magnetically-dominated outflow into two beamed magnetic-tower jets. This creates a clean narrow channel for the escape of energy from the central engine through the star, as required for GRBs. In addition, the delayed onset of the collimated-explosion phase can explain the production of large quantities of Nickel-56, as suggested by the GRB-Supernova connection. Finally, the prospects for numerical simulations of this scenario are discussed.

We consider a newly-born millisecond magnetar, focusing on its interaction with the dense stellar plasma in which it is initially embedded. We argue that the confining pressure and inertia of the surrounding plasma acts to collimate the magnetar's Poynting-flux-dominated outflow into tightly beamed jets and increases its magnetic luminosity. We propose this process as an essential ingredient in the magnetar model for gamma-ray burst and asymmetric supernova central engines. We introduce the ``pulsar-in-a-cavity'' as an important model problem representing a magnetized rotating neutron star inside a collapsing star. We describe its essential properties and derive simple estimates for the evolution of the magnetic field and the resulting spin-down power. We find that the infalling stellar mantle confines the magnetosphere, enabling a gradual build-up of the toroidal magnetic field due to continuous twisting. The growing magnetic pressure eventually becomes dominant, resulting in a magnetically-driven explosion. The initial phase of the explosion is quasi-isotropic, potentially exposing a sufficient amount of material to $^{56}$Ni-producing temperatures to result in a bright supernova. However, if significant expansion of the star occurs prior to the explosion, then very little $^{56}$Ni is produced and no supernova is expected. In either case, hoop stress subsequently collimates the magnetically-dominated outflow, leading to the formation of a magnetic tower. After the star explodes, the decrease in bounding pressure causes the magnetic outflow to become less beamed. However, episodes of late fallback can reform the beamed outflow, which may be responsible for late X-ray flares.

I present a novel view on the problem of solar coronal heating. In my picture, coronal heating should be viewed as a self-regulating process that works to keep the coronal plasma marginally collisionless. The self-regulating mechanism is based on the interplay between two effects: (1) Plasma density controls coronal energy release via the transition between the slow collisional Sweet-Parker regime and the fast collisionless reconnection regime; (2) In turn, coronal energy release through reconnection leads to an increase in the ambient plasma density via chromospheric evaporation, which temporarily shuts off any subsequent reconnection involving the newly-reconnected loops.

The main subject of my talk is the question: in what kind of astrophysical systems magnetic reconnection is interesting and/or important? To address this question, I first put forward three general criteria for selecting the relevant astrophysical environments. Namely, reconnection should be: fast; energetically important; and observable. From this, I deduce that the gas density should be low, so that the plasma is: collisionless; force-free; and optically thin. Thus, for example, the requirement that reconnection is fast implies that Petschek's reconnection mechanism must be operating, which is possible, apparently, only in the collisionless regime. Next, I argue that the force-free condition implies that the magnetic field be produced in, and anchored by, a nearby dense massive object, e.g., a star or a disk, strongly stratified by gravity. I then stress the importance of field-line opening (e.g., by differential rotation) as a means to form a reconnecting current sheet. Correspondingly, I suggest the Y-point helmet streamer as a generic prototypical magnetic configuration relevant to large-scale reconnection in astrophysics. Finally, I discuss several specific astrophysical systems where the above criteria are met: stellar coronae, magnetically-interacting star--disk systems, and magnetized coronae above turbulent accretion disks. In the Appendix, I apply the ideas put forward in this talk to the solar coronal heating problem.

A quadrupole pattern of the out-of-plane component of the magnetic field inside a reconnection region is seen as an important signature of the Hall-magnetohydrodynamic (Hall-MHD) regime of reconnection. It has been first observed in numerical simulations and just recently confirmed in the MRX (Magnetic Reconnection Experiment) [M. Yamada, H. Ji, S. Hsu, T. Carter, R. Kulsrud, N. Bertz, F. Jobes, Y. Ono, and F. Perkins, Phys. Plasmas 4, 1936 (1997)] and also seen in spacecraft observations of Earth's magnetosphere. In this study, the physical origin of the quadrupole field is analyzed and traced to a current of electrons that flows along the lines in and out of the inner reconnection region to maintain charge neutrality. The role of the quadrupole magnetic field in the overall dynamics of the reconnection process is discussed. In addition, the bipolar poloidal electric field is estimated and its effect on ion motions is emphasized.

We propose a magnetic mechanism for the collimated explosion of a massive star relevant for GRBs, XRFs and asymmetric supernovae. We apply Lynden-Bell's magnetic tower scenario to the interior of a massive rotating star after the core has collapsed to form a black hole with an accretion disk or a millisecond magnetar acting as a central engine. We solve the force-free Grad-Shafranov equation to calculate the magnetic structure and growth of a tower embedded in a stellar environment. The pressure of the toroidal magnetic field, continuously generated by differential rotation of the central engine, drives a rapid expansion which becomes vertically collimated after lateral force balance with the surrounding gas pressure is reached. The collimation naturally occurs because hoop stress concentrates magnetic field toward the rotation axis and inhibits lateral expansion. This leads to the growth of a self-collimated magnetic tower. When embedded in a massive star, the supersonic expansion of the tower drives a strong bow shock behind which an over-pressured cocoon forms. The cocoon confines the tower by supplying collimating pressure and provides stabilization against disruption due to MHD instabilities. Because the tower consists of closed field lines starting and ending on the central engine, mixing of baryons from the cocoon into the tower is suppressed. The channel cleared by the growing tower is thus plausibly free of baryons and allows the escape of magnetic energy from the central engine through the star. While propagating down the stellar density gradient, the tower accelerates and becomes relativistic. During the expansion, fast collisionless reconnection becomes possible resulting in dissipation of magnetic energy which may be responsible for GRB prompt emission.

We suggest a possible mechanism of efficient transport of the large-scale external magnetic flux inward through a turbulent accretion disk. The outward drift by turbulent diffusion which limits the effectiveness of this process is reduced if the external flux passes through the disk in concentrated patches. Angular momentum loss by a magnetocentrifugal wind from these patches of strong field adds to the inward drift. We propose that magnetic flux accumulating in this way at the center of the disk provides the `second parameter' determining X-ray spectral states of black hole binaries, and the presence of relativistic outflows.

I consider a stationary axisymmetric force-free degenerate magnetosphere of a rotating Kerr black hole surrounded by a thin Keplerian infinitely-conducting accretion disk. I focus on the closed-field geometry with a direct magnetic coupling between the disk and the event horizon. I first present a simple physical argument that shows how the black hole's rotation limits the radial extent of the force-free link. I then confirm this result by solving numerically the general-relativistic force-free Grad--Shafranov equation in the magnetosphere, using the regularity condition at the inner light cylinder to determine the poloidal current. I indeed find that force-free solutions exist only when the magnetic link between the hole and the disk has a limited extent on the disk surface. I chart out the maximum allowable size of this magnetically-connected part of the disk as a function of the black hole spin. I also compute the angular momentum and energy transfer between the hole and the disk that takes place via the direct magnetic link. I find that both of these quantities grow rapidly and that their deposition becomes highly concentrated near the inner edge of the disk as the black hole spin is increased.

We consider a stationary axisymmetric force-free degenerate magnetosphere of a rotating Kerr black hole surrounded by a thin Keplerian infinitely-conducting accretion disk. We focus on the closed-field geometry characterize by a direct magnetic coupling between the disk and the hole's event horizon. We first argue that the hole's rotation necessarily limits the radial extent of the force-free link on the disk surface: the faster the hole rotates, the smaller the magnetically-connected inner part of the disk has to be. We then show that this is indeed the case by solving numerically the Grad--Shafranov equation--the main differential equation describing the structure of the magnetosphere. An important element in our approach is the use of the regularity condition at the inner light cylinder to fix the poloidal current as a function of the poloidal magnetic flux. As an outcome of our computations, we are able to chart out the maximum allowable size of the portion of the disk that is magnetically connected to the hole as a function of the black hole spin. We also calculate the angular momentum and energy transfer between the hole and the disk that takes place via the direct magnetic link. We find that both of these quantities grow rapidly and that their deposition becomes highly concentrated near the inner edge of the disk as the black hole spin is increased.

In this paper I study the magnetosphere of a black hole that is connected by the magnetic field to a thin conducting Keplerian disk. I consider the case of a Schwarzschild black hole only, leaving the more interesting but difficult case of a Kerr black hole to a future study. I assume that the magnetosphere is ideal, stationary, axisymmetric, and force-free. I pay a special attention to the two singular surfaces present in the system, i.e., the event horizon and the inner light cylinder; I use the regularity condition at the light cylinder to determine the poloidal electric current as a function of poloidal magnetic flux. I solve numerically the Grad--Shafranov equation, which governs the structure of the magnetosphere, for two cases: the case of a nonrotating disk and the case of a Keplerian disk. I find that, in both cases, the poloidal flux function on the horizon matches a simple analytical expression corresponding to a radial magnetic field that is uniform on the horizon. Using this result, I express the poloidal current as an explicit function of the flux and find a perfect agreement between this analytical expression and my numerical results.

In this review I consider modern theoretical models of coupled star--disk magnetospheres. I discuss a number of models, both stationary and time-dependent, and examine what physical conditions govern the selection of a preferred model.

We investigate the axisymmetric magnetosphere of an aligned rotating magnetic dipole surrounded by an ideal force-free plasma. We concentrate on the magnetic field structure around the point of intersection of the separatrix between the open and closed field-line regions and the equatorial plane. We first study the case where this intersection point is located at the Light Cylinder. We find that in this case the separatrix equilibrium condition implies that all the poloidal current must return to the pulsar in the open-field region, i.e., that there should be no finite current carried by the separatrix/equator current sheet. We then perform an asymptotic analysis of the pulsar equation near the intersection point and find a unique self-similar solution; however, a Light Surface inevitably emerges right outside the Light Cylinder. We then perform a similar analysis for the situation where the intersection point lies somewhere inside the Light Cylinder, in which case a finite current flowing along the separatrix and the equator is allowed. We find a very simple behavior in this case, characterized by a 90-degree angle between the separatrix and the equator and by finite vertical field in the closed-field region. Finally, we discuss the implications of our results for global numerical studies of pulsar magnetospheres.

Recent numerical simulations of magnetic reconnection in two dimensions have shown that, when the resistivity is strongly localized, the reconnection region develops a Petschek-like structure, with the width of the inner diffusion region being of the order of the resistivity localization scale. In this paper, we combine this fact with a realistic model for locally-enhanced anomalous resistivity generated by current-driven microturbulence. The result is a qualitative model of the reconnection layer where the size of Petschek's diffusion region and, therefore, the final reconnection rate are determined self-consistently in terms of the main parameters of the functional dependence of anomalous resistivity on the electric current density. We then consider anomalous resistivity due to ion-acoustic turbulence as a particular case. This enables us to express the reconnection region's parameters directly in terms of the basic parameters of the plasma. Finally, we apply this reconnection model to solar flares and obtain specific predictions for typical reconnection times, which are very consistent with observations.

The dynamics of the steady-state Sweet--Parker-type reconnection is analyzed in relativistic regime when energy density in the inflowing region is dominated by magnetic field. The structure of reconnection layer (its thickness, inflow and outflow velocities) depends on the ratio of two large dimensionless parameters of the problem - magnetization parameter $\sigma \gg 1$ (the ratio of the magnetic to particle energy-densities in the inflowing region) and the Lundquist number $S$. The inflow velocity may be relativistic (for $ S< \sigma$) or non-relativistic (for $ S > \sigma$), while the outflowing plasma is moving always relativisticly. For extremely magnetized plasmas with $\sigma \geq S^2$, the inflow four-velocity becomes of the order of the \Alfven four-velocity.

This paper discusses a quasi-static evolution of a force-free magnetic field under slow sheared footpoint motions on the plasma's boundary, an important problem with applications to the solar and accretion disk coronae. The main qualitative features of the evolution (such as field-line expansion and opening) are considered and a comparison is made between two different geometrical settings: the Cartesian case with translational symmetry along a straight line, and the axisymmetric case with axial symmetry around the rotation axis. The main question addressed in the paper is whether a continuous sequence of force-free equilibria describes the evolution at arbitrarily large values of the footpoint displacement or the sequence ends abruptly and the system exhibits a loss of equilibrium at a finite footpoint displacement. After a formal description of the problem, a review/discussion of the extensive previous work on the subject is given. After that, a series of simple scaling-type arguments, explaining the key essential reason for the main qualitative difference between the two geometry types, is presented. It is found that, in the Cartesian case, force-free equilibria exist at arbitrarily large values of shear and the field approaches the open state only at infinite shear, whereas in the axisymmetric case the field opens up already at a finite shear.

In this paper I analyze the process of formation of thin current structures in the magnetosphere of a conducting accretion disk in response to the field-line twisting brought about by the rotation of the disk relative to the central star. I consider an axisymmetric force-free magnetically-linked star--disk configuration and investigate the expansion of the poloidal field lines and partial field-line opening caused by the differential rotation between the star and a nonuniformly-rotating disk. I present a simple analytical model that describes the asymptotic behavior of the field in the strong-expansion limit. I demonstrate the existence of a finite (of order one radian) critical twist angle, beyond which the poloidal field starts inflating very rapidly. If the relative star--disk twist is enhanced locally, in some finite part of the disk (which may be the case for a Keplerian disk that extends inward significantly closer to the central star than the corotation radius), then, as the twist is increased by a finite amount, the field approaches a partially-open configuration, with some field lines going out to infinity. Simultaneous with this partial field opening, a very thin, radially extended current layer forms, thus laying out a way towards reconnection in the disk magnetosphere. Reconnection, in turn, leads to a very interesting scenario for a quasi-periodic behavior of magnetically-linked star--disk systems with successive cycles of field inflation, opening, and reconnection.

Disk accretion onto a magnetized star occurs in a variety of astrophysical contexts, from young stars to X-ray pulsars. The magnetohydrodynamic interaction between the stellar field and the accreting matter can have a strong effect on the disk structure, the transfer of mass and angular momentum between the disk and the star, and the production of bipolar outflows, e.g., plasma jets. We study a key element of this interaction - the time evolution of the magnetic field configuration brought about by the relative rotation between the disk and the star - using simplified, largely semianalytic, models. We first discuss the rapid inflation and opening up of the magnetic field lines in the corona above the accretion disk, which is caused by the differential rotation twisting. Then we consider additional physical effects that tend to limit this expansion, such as the effect of plasma inertia and the possibility of reconnection in the disk's corona, the latter possibly leading to repeated cycles in the evolution. We also derive the condition for the existence of a steady state for a resistive disk and conclude that a steady state configuration is not realistically possible. Finally, we generalize our analysis of the opening of magnetic field lines by using a non-self-similar numerical model that applies to an arbitrarily rotating (e.g. keplerian) disk.

In this paper we present a two-dimensional, time dependent, numerical simulation of a reconnection current layer in incompressible resistive magnetohydrodynamics with uniform resistivity in the limit of very large Lundquist numbers. We use realistic boundary conditions derived consistently from the outside magnetic field, and we also take into account the effect of the back pressure from flow into the the separatrix region. We find that within a few Alfven times the system evolves from an arbitrary initial state to a steady state consistent with the Sweet--Parker model, even if the initial state is Petschek-like.

In this paper we present a two-dimensional numerical simulation of a reconnection current layer in incompressible resistive magnetohydrodynamics with uniform resistivity in the limit of very large Lundquist numbers. We use realistic boundary conditions derived consistently from the outside magnetic field, and we also take into account the effect of the backpressure from the flow into the separatrix region. We find that within a few Alfven times the system reaches a steady state consistent with the Sweet--Parker model, even if the initial state is Petschek-like.