We study the evolution of turbulent magnetic fields from a topological point of view, invoking commonplace mathematical tools from general topology and dynamical systems theory which connect magnetic field evolution to time reversal invariance, entropy increase and the second law of thermodynamics. We show that in fact magnetic topology is well-defined only in the phase space corresponding to a dynamical system governed by the induction equation. Hence the field's topology and stochasticity can be studied in terms of the corresponding phase space trajectories rather than the field lines in real Euclidean space. In fact, our results suggest that magnetic field lines should not be taken too literally because their existence and uniqueness and more importantly continuity in time require strong mathematical conditions, hardly satisfied in astrophysical systems. As for magnetic topology change, it is shown that the phase space topology is preserved in time for a magnetic field which, besides satisfying few continuity conditions, solves a time reversal invariant induction equation. What breaks the time symmetry in the induction equation is the presence of non-ideal plasma effects at small scales such as resistivity, which results from random collisions between diffusing electrons and other particles. The small scale, stochastic disturbances produced thereby are super-linearly amplified by Richardson diffusion in the turbulent cascade, which are eventually manifested as large scale reconnection events, somehow similar to stretching quantum fluctuations during inflation to seed large scale cosmological structures. This suggests that reconnection is rooted in the second law of thermodynamics that dictates entropy increase which in turn breaks the time symmetry.

Local magnetic reversals are an inseparable part of magnetohydrodynamic (MHD) turbulence whose collective outcome on an arbitrary scale in the inertial range may lead to a global stochastic reconnection event with a rate independent of small scale physics. We show that this picture is intimately related to the nanoflare theory of the solar corona. First, we argue that due to stochastic flux freezing, a generalized version of flux freezing in turbulence, the magnetic field follows the turbulent flow in a statistical sense. Bending and stretching an initially smooth field, therefore, the turbulence generally increases the magnetic spatial complexity. Strong magnetic shears associated with such a highly tangled field can trigger local reversals and field annihilations that convert magnetic energy into kinetic and thermal energy respectively. The former maintains the turbulence, which incidentally continues to entangle the field completing the cycle, while the latter enhances the heat generation in the dissipative range. We support this theoretical picture invoking recent analytical and numerical studies which suggest a correlation between magnetic complexity and magnetic energy dissipation. The amplification of multiple local, in-phase reversals by super-linear Richardson diffusion may initiate a global reconnection at larger scales, however, even in the absence of such a global stochastic reconnection, the small scale reversals will continue to interact with the turbulence. We employ conventional scaling laws of MHD turbulence to illustrate that these local events are indeed efficient in both enhancing the turbulence and generating heat. Finally, using an MHD numerical simulation, we show that the time evolution of the magnetic complexity is statistically correlated with the kinetic energy injection rate and/or magnetic-to-thermal energy conversion rate.

The level of spatial complexity associated with a given vector field on an arbitrary range of scales \iffalse ${\bf F(x}, t)$ can be quantified by a simple, time-dependent function $S(t)={1\over 2}(1-\hat{\bf F}_l.\hat{\bf F}_L)_{rms}$ with ${\bf F}_l$ (${\bf F}_L$) being the average field in a volume of scale $l$ ($L>l$) and the unit vector defined as $\hat{\bf F}={\bf F}/|{\bf F}|$. Thus,\fi can be quantified by a simple, scale-dependent function of time; $0\leq S(t)\leq 1$. Previous work has invoked kinetic and magnetic complexities, associated with velocity and magnetic fields ${\bf u(x}, t)$ and ${\bf B(x}, t)$, to study magnetic reconnection and diffusion in turbulent and magnetized fluids. In this paper, using the coarse-grained momentum equation, we argue that the fluid jets associated with magnetic reconnection events at an arbitrary scale $l$ in the turbulence inertial range are predominantly driven by the Lorentz force ${\bf{N}}_l=({\bf j\times B})_l-{\bf j}_l\times {\bf B}_l$. This force, is induced by the subscale currents and is analogous to the turbulent electromotive force ${\cal E}_l=({\bf u\times B})_l-{\bf u}_l\times {\bf B}_l$ in dynamo theories. Typically, high (low) magnetic complexities during reconnection imply large (small) spatial gradients for the magnetic field, i.e., strong (weak) Lorentz forces ${\bf N}_l$. Reconnection launches jets of fluid, hence the rate of change of kinetic complexity is expected to strongly correlate with the power injected by the Lorentz force ${\bf N}_l$. We test this prediction using an incompressible, homogeneous magnetohydrodynamic (MHD) simulation and associate it with previous results. It follows that the stronger (weaker) the turbulence, the more (less) complex the magnetic field and the stronger (weaker) the reconnection field and thus the ensuing reconnection.

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

We study the topology associated with physical vector and scalar fields. A mathematical object, e.g., a ball, can be continuously deformed, without tearing or gluing, to make other topologically equivalent objects, e.g., a cube or a solid disk. If tearing or gluing get involved, i.e., the deformation is not continuous anymore, the initial topology will consequently change giving rise to a topologically distinct object, e.g., a torus. This simple concept in general topology may be employed in the study of physical systems described by fields. Instead of continuously deforming objects, we can take a continuously evolving field, with an appropriately defined topology, such that the topology remains unchanged in time unless the system undergoes an important physical change, e.g., a transition to a different energy state. For instance, a sudden change in the magnetic topology in an energetically relaxing plasma, a process called reconnection, strongly affects the dynamics, e.g., it is involved in launching solar flares and generating large scale magnetic fields in astrophysical objects. In this topological formalism, the magnetic topology in a plasma can spontaneously change due to the presence of dissipative terms in the induction equation which break its time symmetry. We define a topology for the vector field $\bf F$ in the phase space $(\bf x, F)$. As for scalar fields represented by a perfect fluid, e.g., the inhomogeneous inflaton or Higgs fields, the fluid velocity $\bf u$ defines the corresponding topology. The vector field topology in its corresponding phase space $(\bf x, F)$ will be preserved in time if certain conditions including time reversal invariance are satisfied by the field and its governing differential equation.

Over the last few years it became clear that turbulent magnetic reconnection and magnetized turbulence are inseparable. It was not only shown that reconnection is responsible for violating the frozen-in condition in turbulence, but also that stochastic reconnection in 3D generates turbulence by itself. The actual mechanism responsible for this driving is still unknown. Processes such tearing mode or Kelvin-Helmholtz, among other plasma instabilities, could generate turbulence from irregular current sheets. We address the nature of driving mechanism for this process and consider a relative role of tearing and Kelvin-Helmholtz instabilities for the process of turbulence generation. In particular, we analyze the conditions for development of these two instabilities within three-dimensional reconnection regions. We show that both instabilities can excite turbulence fluctuations in reconnection regions. However, tearing mode has relatively slow growth rate, and at later times it becomes partially suppressed by transverse to the current sheet component of magnetic field, generated during the growth of turbulent fluctuations. On the contrary, the Kelvin-Helmholtz instability establishes quickly in the outflow region, and at later times it dominates the turbulence generation comparing to the contribution from tearing mode. Our results demonstrate that the tearing instability is subdominant compared to the the Kelvin-Helmholtz instability in terms of generation of turbulence in the 3D reconnection layers and therefore the self-driven reconnection is turbulent reconnection with tearing instability being important only at the initial stage of the reconnection.

Previous work has introduced scale-split energy density \psi for a given vector field B in order to quantify the field stochasticity S(t). Application to turbulent magnetic fields leads to the prediction that tangling magnetic field by turbulence increases magnetic stochasticity. An increasing stochasticity in turn leads to disalignments of the coarse-grained fields at smaller scales thus they average to weaker fields at larger scales upon coarse-graining. The field's resistance against tanglement by the turbulence may lead at some point to its sudden slippage through the fluid, decreasing the stochasticity and increasing the energy density. Thus the maxima (minima) of magnetic stochasticity are expected to approximately coincide with the minima (maxima) of energy density, occurrence of which corresponds to slippage of the magnetic field through the fluid. Field-fluid slippage, on the other hand, has been already found to be intimately related to magnetic reconnection. In this paper, we test these theoretical predictions numerically using a homogeneous, incompressible magnetohydrodynamic (MHD) simulation. Apart from expected small scale deviations, possibly due to e.g., intermittency and strong field annihilation, the theoretically predicted global relationship between stochasticity and magnetic energy is observed in different sub-volumes of the simulation box. This may indicate ubiquitous local field-fluid slippage and small scale reconnection events in MHD turbulence. We also show that the maximum magnetic stochasticity, i.e., \partial_t S(t)=0 & \partial^2_t S(t)<0, leads to sudden increases in kinetic stochasticity level which may correspond to fluid jets driven by the reconnecting field lines, i.e., reconnection. This suggests a new mathematical approach to the reconnection problem.

Magnetic reconnection, topological change in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence magnetic lines constantly reconnect everywhere at all scales, making magnetic reconnection an intrinsic part of turbulent cascade. We note that this is inconsistent with the usual practice of regarding magnetic lines as persistent dynamical elements. A number of theoretical, numerical, and observational studies, starting with Lazarian & Vishniac (1999), demonstrated that 3D turbulence makes magnetic reconnection fast and that these two processes are intrinsically connected. We discuss the dramatic violation of the textbook concept of magnetic flux-freezing in the presence of turbulence and demonstrate that in the presence of turbulence the plasma effects are subdominant to turbulence as far as the magnetic reconnection is concerned. This justifies an MHD-like treatment of magnetic reconnection at scales much larger than the relevant plasma scales. We discuss numerical and observational evidences supporting the turbulent reconnection model. In particular, we show that tearing reconnection is suppressed in 3D and, unlike the 2D case, the 3D reconnection induces turbulence that makes reconnection independent of resistivity. We show that turbulent reconnection dramatically affects the key astrophysical processes, e.g., star formation, turbulent dynamo, acceleration of cosmic rays. We provide criticism of the concept of "reconnection-mediated turbulence" and explain why turbulent reconnection is very different from enhanced turbulent resistivity and hyper-resistivity, and why the latter has fatal conceptual flaws.

Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.

We develop a quantitative relationship between magnetic diffusion and the level of randomness, or stochasticity, of the diffusing magnetic field in a magnetized medium. A general mathematical formulation of magnetic stochasticity in turbulence has been developed in previous work in terms of the ${\cal L}_p$-norm $S_p(t)={1\over 2}|| 1-\hat{\bf B}_l.\hat{\bf B}_L||_p$, $p$th order magnetic stochasticity of the stochastic field ${\bf B}({\bf x}, t)$, based on the coarse-grained fields, ${\bf B}_l$ and ${\bf B}_L$, at different scales, $l\neq L$. For laminar flows, stochasticity level becomes the level of field self-entanglement or spatial complexity. In this paper, we establish a connection between magnetic stochasticity $S_p(t)$ and magnetic diffusion in magnetohydrodynamic (MHD) turbulence and use a homogeneous, incompressible MHD simulation to test this prediction. Our results agree with the well-known fact that magnetic diffusion in turbulent media follows the super-linear Richardson dispersion scheme. This is intimately related to stochastic magnetic reconnection in which super-linear Richardson diffusion broadens the matter outflow width and accelerates the reconnection process.

We review the basic concepts of magnetic reconnection and propose a general framework for the astrophysical reconnection at large scales. Magnetic reconnection is the rearrangement of magnetic field topology. The conventional Sweet-Parker scheme and some of its extensions presume a non-turbulent medium and generally produce very slow reconnection or an unstable configuration. However, the assumption of laminar flow is unrealistic in astrophysics since, even in an initially quiet environment, magnetic reconnection by itself can drive turbulence. The resulting turbulence has the potential to enhance the reconnection rate. This can lead to an unstable feedback loop as reconnection drives turbulence and turbulence drives reconnection. Stochastic reconnection was proposed, and subsequently tested by numerical simulations, for high $\beta$ plasmas with a magnetic Prandtl number of order unity, $Pr_m\sim 1$. This model predicts reconnection speeds comparable to the large scale turbulent eddy velocity. A recent study of stochastic reconnection for $Pr_m>1$ has shown that the width of the outflow layer and the ejection velocity of matter from the reconnection region seem to be unaffected by viscosity in typical astrophysical systems. However if $Pr_m>1$ viscosity can suppress small scale reconnection events near and below the Kolmogorov or viscous damping scale. This will produce a threshold for the suppression of large scale reconnection by viscosity when $Pr_m$ is larger than the root of the Reynolds number, $Pr_m>\sqrt{Re}$. For $Pr_m>1$ this leads to the spectral index $\sim-4/3$ for length scales between the viscous dissipation scale and eddies larger by roughly $Pr_m^{3/2}$.

Magnetic reconnection is a process that changes magnetic field topology in highly conducting fluids. Within the standard Sweet-Parker model, this process would be too slow to explain observations (e.g. solar flares). In reality, the process must be ubiquitous as astrophysical fluids are magnetized and motions of fluid elements necessarily entail crossing of magnetic frozen-in field lines and magnetic reconnection. In the presence of turbulence, the reconnection is independent of microscopic plasma properties, and may be much faster than previously thought, as proposed in Lazarian & Vishniac (1999) and tested in Kowal et al. (2009, 2012). However, the considered turbulence in the Lazarian-Vishniac model was imposed externally. In this work we consider reconnection-driven magnetized turbulence in realistic three-dimensional geometry initiated by stochastic noise. We demonstrate through numerical simulations that the stochastic reconnection is able to self-generate turbulence through interactions between its outflows. We analyze the statistical properties of velocity fluctuations using power spectra and anisotropy scaling, which demonstrates that the reconnection produces Kolmogorov-like turbulence, compatible with Goldreich-Sridhar (1995) model. Anisotropy statistics are, however, strongly affected by the dynamics of reconnection outflows. Once the broad turbulent region is formed, the typical anisotropy scaling $l_\parallel \propto l_\perp^{2/3}$ is formed, especially for high resolution models, were the broader range of scales is available. The decay of reconnection outflows to turbulent-like fluctuations, characterized by different anisotropy scalings, strongly depends on $\beta$ plasma. Moreover, the estimated reconnection rates are weakly dependent on the resolution, suggesting that no external processes are required to make reconnection fast.

In situ spacecraft data on the solar wind show events identified as magnetic reconnection with outflows and apparent "`$X$-lines" $10^{3-4}$ times ion scales. To understand the role of turbulence at these scales, we make a case study of an inertial-range reconnection event in a magnetohydrodynamic (MHD) simulation. We observe stochastic wandering of field-lines in space, breakdown of standard magnetic flux-freezing due to Richardson dispersion, and a broadened reconnection zone containing many current sheets. The coarse-grain magnetic geometry is like large-scale reconnection in the solar wind, however, with a hyperbolic flux-tube or "$X$-line" extending over integral length-scales.

Magnetic reconnection is a process of magnetic field topology change, which is one of the most fundamental processes in magnetized plasmas. In most astrophysical environments the Reynolds numbers are large and therefore the transition to turbulence is inevitable. This turbulence must be taken into account for any theory of magnetic reconnection, since the initially laminar configurations can transit to the turbulence state, what is demonstrated by 3D high resolution numerical simulations. We discuss ideas of how turbulence can modify reconnection with the focus on the Lazarian & Vishniac (1999) reconnection model and present numerical evidence supporting the model and demonstrate that it is closely connected to the concept of Richardson diffusion and compatible with the Lagrangian dynamics of magnetized fluids. We point out that the Generalized Ohm's Law, that accounts for turbulent motion, predicts the subdominance of the microphysical plasma effects for a realistically turbulent media. We show that on of the most dramatic consequences of turbulence is the violation of the generally accepted notion of magnetic flux freezing. This notion is a corner stone of most theories dealing with magnetized plasmas and therefore its change induces fundamental shifts in accepted paradigms like turbulent reconnection entailing the diffusion process that is essential for understanding star formation. We argue, that at sufficiently high Reynolds numbers the process of tearing reconnection should transfer to turbulent reconnection. We discuss flares predicted by turbulent reconnection and relate them to solar flares and gamma ray bursts. We analyze solar observations, measurements in the solar wind or heliospheric current sheet, and show their correspondence with turbulent reconnection predictions. Finally, we discuss 1st Order Fermi acceleration as a natural consequence of the turbulent reconnection.

Magnetic reconnection is a fundamental process of magnetic field topology change. We analyze the connection of this process with turbulence which is ubiquitous in astrophysical environments. We show how Lazarian & Vishniac (1999) model of turbulent reconnection is connected to the experimentally proven concept of Richardson diffusion and discuss how turbulence violates the generally accepted notion of magnetic flux freezing. We note that in environments that are laminar initially turbulence can develop as a result of magnetic reconnection and this can result in flares of magnetic reconnection in magnetically dominated media. In particular, magnetic reconnection can initially develop through tearing, but the transition to the turbulent state is expected for astrophysical systems. We show that turbulent reconnection predictions corresponds to the Solar and solar wind data.

Magnetic reconnection is a process that changes magnetic field topology in highly conducting fluids. Traditionally, magnetic reconnection was associated mostly with solar flares. In reality, the process must be ubiquitous as astrophysical fluids are magnetized and motions of fluid elements necessarily entail crossing of magnetic frozen in field lines and magnetic reconnection. We consider magnetic reconnection in realistic 3D geometry in the presence of turbulence. This turbulence in most astrophysical settings is of pre-existing nature, but it also can be induced by magnetic reconnection itself. In this situation turbulent magnetic field wandering opens up reconnection outflow regions, making reconnection fast. We discuss Lazarian \& Vishniac (1999) model of turbulent reconnection, its numerical and observational testings, as well as its connection to the modern understanding of the Lagrangian properties of turbulent fluids. We show that the predicted dependences of the reconnection rates on the level of MHD turbulence make the generally accepted Goldreich \& Sridhar (1995) model of turbulence self-consistent. Similarly, we argue that the well-known Alfv\'en theorem on flux freezing is not valid for the turbulent fluids and therefore magnetic fields diffuse within turbulent volumes. This is an element of magnetic field dynamics that was not accounted by earlier theories. For instance, the theory of star formation that was developing assuming that it is only the drift of neutrals that can violate the otherwise perfect flux freezing, is affected and we discuss the consequences of the turbulent diffusion of magnetic fields mediated by reconnection.

Astrophysical media are turbulent and therefore reconnection should be treated in the presence of pre-existing turbulence. We consider the model of fast magnetic reconnection in Lazarian & Vishniac (1999) which predicts that the rate of reconnection is controlled by the intensity and the injection scale of turbulent motions. We provide new evidence of successful testing of the model and argue that the model presents a generic set up for astrophysical reconnection events. We study particle acceleration that takes place in volumes of driven turbulence as well turbulent volumes in the presence of large scale reconnection. We show that in the latter case the acceleration is of the first order Fermi type thus supporting the model of acceleration proposed in Gouveia dal Pino & Lazarian (2005).

We study a model of fast magnetic reconnection in the presence of weak turbulence proposed by Lazarian and Vishniac (1999) using three-dimensional direct numerical simulations. The model has been already successfully tested in Kowal et al. (2009) confirming the dependencies of the reconnection speed $V_{rec}$ on the turbulence injection power $P_{inj}$ and the injection scale $l_{inj}$ expressed by a constraint $V_{rec} \sim P_{inj}^{1/2} l_{inj}^{3/4}$ and no observed dependency on Ohmic resistivity. In Kowal et al. (2009), in order to drive turbulence, we injected velocity fluctuations in Fourier space with frequencies concentrated around $k_{inj}=1/l_{inj}$, as described in Alvelius (1999). In this paper we extend our previous studies by comparing fast magnetic reconnection under different mechanisms of turbulence injection by introducing a new way of turbulence driving. The new method injects velocity or magnetic eddies with a specified amplitude and scale in random locations directly in real space. We provide exact relations between the eddy parameters and turbulent power and injection scale. We performed simulations with new forcing in order to study turbulent power and injection scale dependencies. The results show no discrepancy between models with two different methods of turbulence driving exposing the same scalings in both cases. This is in agreement with the Lazarian and Vishniac (1999) predictions. In addition, we performed a series of models with varying viscosity $\nu$. Although Lazarian and Vishniac (1999) do not provide any prediction for this dependence, we report a weak relation between the reconnection speed with viscosity, $V_{rec}\sim\nu^{-1/4}$.

Astrophysical fluids are generically turbulent and this must be taken into account for most transport processes. We discuss how the preexisting turbulence modifies magnetic reconnection and how magnetic reconnection affects the MHD turbulent cascade. We show the intrinsic interdependence and interrelation of magnetic turbulence and magnetic reconnection, in particular, that strong magnetic turbulence in 3D requires reconnection and 3D magnetic turbulence entails fast reconnection. We follow the approach in Eyink, Lazarian & Vishniac 2011 to show that the expressions of fast magnetic reconnection in Lazarian & Vishniac 1999 can be recovered if Richardson diffusion of turbulent flows is used instead of ordinary Ohmic diffusion. This does not revive, however, the concept of magnetic turbulent diffusion which assumes that magnetic fields can be mixed up in a passive way down to a very small dissipation scales. On the contrary, we are dealing the reconnection of dynamically important magnetic field bundles which strongly resist bending and have well defined mean direction weakly perturbed by turbulence. We argue that in the presence of turbulence the very concept of flux-freezing requires modification. The diffusion that arises from magnetic turbulence can be called reconnection diffusion as it based on reconnection of magnetic field lines. The reconnection diffusion has important implications for the continuous transport processes in magnetized plasmas and for star formation. In addition, fast magnetic reconnection in turbulent media induces the First order Fermi acceleration of energetic particles, can explain solar flares and gamma ray bursts. However, the most dramatic consequence of these developments is the fact that the standard flux freezing concept must be radically modified in the presence of turbulence.

Magnetic field-lines in astrophysical plasmas are expected to be frozen-in at scales larger than the ion gyroradius. The rapid reconnection of magnetic flux structures with dimensions vastly larger than the gyroradius requires a breakdown in the standard Alfv\'en flux-freezing law. We attribute this breakdown to ubiquitous MHD plasma turbulence with power-law scaling ranges of velocity and magnetic energy spectra. Lagrangian particle trajectories in such environments become "spontaneously stochastic", so that infinitely-many magnetic field-lines are advected to each point and must be averaged to obtain the resultant magnetic field. The relative distance between initial magnetic field lines which arrive to the same final point depends upon the properties of two-particle turbulent dispersion. We develop predictions based on the phenomenological Goldreich & Sridhar theory of strong MHD turbulence and on weak MHD turbulence theory. We recover the predictions of the Lazarian & Vishniac theory for the reconnection rate of large-scale magnetic structures. Lazarian & Vishniac also invoked "spontaneous stochasticity", but of the field-lines rather than of the Lagrangian trajectories. More recent theories of fast magnetic reconnection appeal to microscopic plasma processes that lead to additional terms in the generalized Ohm's law, such as the collisionless Hall term. We estimate quantitatively the effect of such processes on the inertial-range turbulence dynamics and find them to be negligible in most astrophysical environments. For example, the predictions of the Lazarian-Vishniac theory are unchanged in Hall MHD turbulence with an extended inertial range, whenever the ion skin depth $\delta_i$ is much smaller than the turbulent integral length or injection-scale $L_i.$

Magnetic field embedded in a perfectly conducting fluid preserves its topology for all time. Although ionized astrophysical objects, like stars and galactic disks, are almost perfectly conducting, they show indications of changes in topology, `magnetic reconnection', on dynamical time scales. Reconnection can be observed directly in the solar corona, but can also be inferred from the existence of large scale dynamo activity inside stellar interiors. Solar flares and gamma ray busts are usually associated with magnetic reconnection. Previous work has concentrated on showing how reconnection can be rapid in plasmas with very small collision rates. Here we present numerical evidence, based on three dimensional simulations, that reconnection in a turbulent fluid occurs at a speed comparable to the rms velocity of the turbulence, regardless of the value of the resistivity. In particular, this is true for turbulent pressures much weaker than the magnetic field pressure so that the magnetic field lines are only slightly bent by the turbulence. These results are consistent with the proposal by Lazarian and Vishniac (1999) that reconnection is controlled by the stochastic diffusion of magnetic field lines, which produces a broad outflow of plasma from the reconnection zone. This work implies that reconnection in a turbulent fluid typically takes place in approximately a single eddy turnover time, with broad implications for dynamo activity and particle acceleration throughout the universe. In contrast, the reconnection in 2D configurations in the presence of turbulence depends on resistivity, i.e. is slow.

Our numerical simulations show that the reconnection of magnetic field becomes fast in the presence of weak turbulence in the way consistent with the Lazarian and Vishniac (1999) model of fast reconnection. We trace particles within our numerical simulations and show that the particles can be efficiently accelerated via the first order Fermi acceleration. We discuss the acceleration arising from reconnection as a possible origin of the anomalous cosmic rays measured by Voyagers.

We present a constructive numerical example of fast magnetic reconnection in a three dimensional periodic box. Reconnection is initiated by a strong, localized perturbation to the field lines. The solution is intrinsically three dimensional, and its gross properties do not depend on the details of the simulations. $\sim 50%$ of the magnetic energy is released in an event which lasts about one Alfven time, but only after a delay during which the field lines evolve into a critical configuration. We present a physical picture of the process. The reconnection regions are dynamical and mutually interacting. In the comoving frame of these regions, reconnection occurs through an X-like point, analogous to Petschek reconnection. The dynamics appear to be driven by global flows, not local processes.

We study two dimensional turbulent magnetic reconnection in a compressible fluid in the gas pressure dominated limit. We use open boundary conditions and start from a Harris current sheet configuration with a uniform total pressure. A small perturbation to the vector potential initiates laminar reconnection at the Sweet-Parker rate, which is allowed to evolve for several dynamical times. Subsequently sub-Alfvenic turbulence is produced through random forcing at small wave numbers. The magnetic field topology near the current sheet is strongly affected by the turbulence. However, we find that the resulting reconnection speed depends on the resistivity. In contrast to previous results in three dimensions, we find no evidence for fast reconnection. The reconnection speed exhibits large variations but the time averages increase smoothly with the strength of the turbulence.

We study the effects of turbulence on magnetic reconnection using 3D numerical simulations. This is the first attempt to test a model of fast magnetic reconnection in the presence of weak turbulence proposed by Lazarian & Vishniac (1999). This model predicts that weak turbulence, generically present in most of astrophysical systems, enhances the rate of reconnection by reducing the transverse scale for reconnection events and by allowing many independent flux reconnection events to occur simultaneously. As a result the reconnection speed becomes independent of Ohmic resistivity and is determined by the magnetic field wandering induced by turbulence. To quantify the reconnection speed we use both an intuitive definition, i.e. the speed of the reconnected flux inflow, as well as a more sophisticated definition based on a formally derived analytical expression. Our results confirm the predictions of the Lazarian & Vishniac model. In particular, we find that Vrec Pinj^(1/2), as predicted by the model. The dependence on the injection scale for some of our models is a bit weaker than expected, i.e. l^(3/4), compared to the predicted linear dependence on the injection scale, which may require some refinement of the model or may be due to the effects like finite size of the excitation region. The reconnection speed was found to depend on the expected rate of magnetic field wandering and not on the magnitude of the guide field. In our models, we see no dependence on the guide field when its strength is comparable to the reconnected component. More importantly, while in the absence of turbulence we successfully reproduce the Sweet-Parker scaling of reconnection, in the presence of turbulence we do not observe any dependence on Ohmic resistivity, confirming that our reconnection is fast.

Simulations of the magnetorotational instability (MRI) in a homogeneous shearing box have shown that the asymptotic strength of the magnetic field declines steeply with increasing resolution. Here I model the MRI driven dynamo as a large scale dynamo driven by the vertical magnetic helicity flux. This growth is balanced by large scale mixing driven by a secondary instability. The saturated magnetic energy density depends almost linearly on the vertical height of the typical eddies. The MRI can drive eddies with arbitrarily large vertical wavenumber, so the eddy thickness is either set by diffusive effects, by the magnetic tension of a large scale vertical field component, or by magnetic buoyancy effects. In homogeneous, zero magnetic flux, simulations only the first effect applies and the saturated limit of the dynamo is determined by explicit or numerical diffusion. The exact result depends on the numerical details, but is consistent with previous work, including the claim that the saturated field energy scales as the gas pressure to the one quarter power (which we interpret as an artifact of numerical dissipation). The magnetic energy density in a homogeneous shearing box will tend to zero as the resolution of the simulation increases, but this has no consequences for the dynamo or for angular momentum transport in real accretion disks. The claim that the saturated state depends on the magnetic Prandtl number may also be an artifact of simulations in which microphysical transport coefficients set the MRI eddy thickness. Finally, the efficiency of the MRI dynamo is a function of the ratio of the Alfv\'en velocity to the product of the pressure scale height and the local shear. As this approaches unity from below the dynamo reaches maximum efficiency.

Astrophysical fluids are generically turbulent, which means that frozen-in magnetic fields are, at least, weakly stochastic. Therefore realistic studies of astrophysical magnetic reconnection should include the effects of stochastic magnetic field. In the paper we discuss and test numerically the Lazarian & Vishniac (1999) model of magnetic field reconnection of weakly stochastic fields. The turbulence in the model is assumed to be subAlfvenic, with the magnetic field only slightly perturbed. The model predicts that the degree of magnetic field stochasticity controls the reconnection rate and that the reconnection can be fast independently on the presence or absence of anomalous plasma effects. For testing of the model we use 3D MHD simulations. To measure the reconnection rate we employ both the inflow of magnetic flux and a more sophisticated measure that we introduce in the paper. Both measures of reconnection provide consistent results. Our testing successfully reproduces the dependences predicted by the model, including the variations of the reconnection speed with the variations of the injection scale of turbulence driving as well as the intensity of driving. We conclude that, while anomalous and Hall-MHD effects in particular circumstances may be important for the initiation of reconnection, the generic astrophysical reconnection is fast due to turbulence, irrespectively of the microphysical plasma effects involved. This conclusion justifies numerical modeling of many astrophysical environments, e.g. interstellar medium, for which plasma-effect-based collisionless reconnection is not applicable.

We discuss the model of magnetic field reconnection in the presence of turbulence introduced by us approximately ten years ago. The model does not require any plasma effects to be involved in order to make the reconnection fast. In fact, it shows that the degree of magnetic field stochasticity controls the reconnection. The turbulence in the model is assumed to be subAlfvenic, with the magnetic field only slightly perturbed. This ensures that the reconnection happens in generic astrophysical environments and the model does not appeal to any unphysical concepts, similar to the turbulent magnetic diffusivity concept, which is employed in the kinematic magnetic dynamo. The interest to that model has recently increased due to successful numerical testings of the model predictions. In view of this, we discuss implications of the model, including the first-order Fermi acceleration of cosmic rays, that the model naturally entails, bursts of reconnection, that can be associated with Solar flares, as well as, removal of magnetic flux during star-formation.

We study the effects of turbulence on magnetic reconnection using three-dimensional numerical simulations. This is the first attempt to test a model of fast magnetic reconnection proposed by Lazarian & Vishniac (1999), which assumes the presence of weak, small-scale magnetic field structure near the current sheet. This affects the rate of reconnection by reducing the transverse scale for reconnection flows and by allowing many independent flux reconnection events to occur simultaneously. We performed a number of simulations to test the dependencies of the reconnection speed, defined as the ratio of the inflow velocity to the Alfven speed, on the turbulence power, the injection scale and resistivity. Our results show that turbulence significantly affects the topology of magnetic field near the diffusion region and increases the thickness of the outflow region. We confirm the predictions of the Lazarian & Vishniac model. In particular, we report the growth of the reconnection speed proportional to ~ V^2, where V is the amplitude of velocity at the injection scale. It depends on the injection scale l as ~ (l/L)^(2/3), where L is the size of the system, which is somewhat faster but still roughly consistent with the theoretical expectations. We also show that for 3D reconnection the Ohmic resistivity is important in the local reconnection events only, and the global reconnection rate in the presence of turbulence does not depend on it.

We present numerical simulations of driven magnetohydrodynamic (MHD) turbulence with weak/moderate imposed magnetic fields. The main goal is to clarify dynamics of magnetic field growth. We also investigate the effects of the imposed magnetic fields on the MHD turbulence, including, as a limit, the case of zero external field. Our findings are as follows. First, when we start off simulations with weak mean magnetic field only (or with small scale random field with zero imposed field), we observe that there is a stage at which magnetic energy density grows linearly with time. Runs with different numerical resolutions and/or different simulation parameters show consistent results for the growth rate at the linear stage. Second, we find that, when the strength of the external field increases, the equilibrium kinetic energy density drops by roughly the product of the rms velocity and the strength of the external field. The equilibrium magnetic energy density rises by roughly the same amount. Third, when the external magnetic field is not very strong (say, less than ~0.2 times the rms velocity when measured in the units of Alfven speed), the turbulence at large scales remains statistically isotropic, i.e. there is no apparent global anisotropy of order B_0/v. We discuss implications of our results on astrophysical fluids.

We present the results of standard one-dimensional test problems in relativistic hydrodynamics using Glimm's (random choice) method, and compare them to results obtained using finite differencing methods. For problems containing profiles with sharp edges, such as shocks, we find Glimm's method yields global errors ~1-3 orders of magnitude smaller than the traditional techniques. The strongest differences are seen for problems in which a shear field is superposed. For smooth flows, Glimm's method is inferior to standard methods. The location of specific features can be off by up to two grid points with respect to an exact solution in Glimm's method, and furthermore curved states are not modeled optimally since the method idealizes solutions as being composed of piecewise constant states. Thus although Glimm's method is superior at correctly resolving sharp features, especially in the presence of shear, for realistic applications in which one typically finds smooth flows plus strong gradients or discontinuities, standard FD methods yield smaller global errors. Glimm's method may prove useful in certain applications such as GRB afterglow shock propagation into a uniform medium.

We report results of a nine-year monitoring effort on the unusual classical Be with a robotic ground-based (APT) B,V-filtered telescope as well as simultaneous observations in 2004 November with this instrument and the RXTE (X-ray) telescope. Our observations disclosed no correlated optical response to the rapid X-ray flares in this star, nor did the star show any sustained flux changes during the course of either of the two monitored nights in either wavelength regime. Our optical light curves reveal that gamma Cas undergoes \~3%-amplitude cycles with lengths of 60--90 days. Over the nine days we monitored the star with the RXTE, the X-ray flux varied in phase with its optical cycle and with an amplitude predicted from correlated optical/X-ray data from an earlier paper. The amplitudes of the V magnitude cycles are 30--40% larger than the B amplitudes, suggesting the seat of the cycles is circumstellar. The cycle lengths constantly change and can damp or grow on timescales as short as 13 days. We have also discovered a coherent period of 1.21581 +/-0.00002 days in all our data, which is consistent only with rotation. The full amplitude of this variation is 0.0060 in both filters. The derived waveform, somewhat surprisingly, is almost sawtooth in shape. This variation probably originates on the star's surface. This circumstance hints at the existence of a strong magnetic field with a complex topology and an associated heterogeneous surface composition.

A photon or neutrino gas--semi-contained by a baryonic species through scattering--comprises a rather peculiar MHD fluid where the magnetic field is truly frozen only to the co-moving volume associated with the mass density. Although radiative diffusion precludes an adiabatic treatment of compressive perturbations, we show that the energy equation may be cast in "quasi-adiabatic" form for exponentially growing non-propagating wave modes. Defining a generalized quasi-adiabatic index leads to a relatively straightforward dispersion relation for non-axisymmetric magnetorotational modes in the horizontal regime when an accretion disk has comparable stress contributions from diffusive and non-diffusive particle species. This analysis is generally applicable to optically thick, neutrino-cooled disks since the pressure contributions from photons, pairs and neutrinos, all have the same temperature dependence whereas only the neutrino component has radiative heat conduction properties on the time and length scales of the instability. We discuss the energy deposition process and the temporal and spatial properties of the ensuing turbulent disk structure on the basis of the derived dispersion relation.

We consider stochastic reconnection in a magnetized, partially ionized medium. Stochastic reconnection is a generic effect, due to field line wandering, in which the speed of reconnection is determined by the ability of ejected plasma to diffuse away from the current sheet along magnetic field lines, rather than by the details of current sheet structure. We consider the limit of weak stochasticity, so that the mean magnetic field energy density is greater than either the turbulent kinetic energy density or the energy density associated with the fluctuating component of the field. We consider field line stochasticity generated through a turbulent cascade, which leads us to consider the effect of neutral drag on the turbulent cascade of energy. In a collisionless plasma, neutral particle viscosity and ion-neutral drag will damp mid-scale turbulent motions, but the power spectrum of the magnetic perturbations extends below the viscous cutoff scale. We give a simple physical picture of the magnetic field structure below this cutoff, consistent with numerical experiments. We provide arguments for the reemergence of the turbulent cascade well below the viscous cut-off scale and derive estimates for field line diffusion on all scales. We note that this explains the persistence of a single power law form for the turbulent power spectrum of the interstellar medium, from scales of tens of parsecs down to thousands of kilometers. We find that under typical conditions in the ISM stochastic reconnection speeds are reduced by the presence of neutrals, but by no more than an order of magnitude.

We present the first unrestricted, three-dimensional relativistic hydrodynamical calculations of the blob of gas associated with the jet producing a gamma-ray burst. We investigate the deceleration phase of the blob corresponding to the time when afterglow radiation is produced, concentrating on the transition in which the relativistic beaming 1/gamma goes from being less than theta, where gamma is the bulk Lorentz factor and theta is the angular width of the jet, to 1/gamma > theta. We study the time dependent evolution of the physical parameters associated with the jet, both parallel to the direction of motion and perpendicular to it. We calculate light curves for observers at varying angles with respect to the velocity vector of the blob, assuming optically thin emission that scales with the local pressure. Our main findings are that (i) gas ahead of the advancing blob does not accrete onto and merge with the blob material but rather flows around the blob, (ii) the decay light curve steepens at a time corresponding roughly to 1/gamma equals theta (in accord with earlier studies), and (iii) the rate of decrease of the forward component of momentum in the blob is well-fit by a simple model in which the gas in front of the blob exerts a drag force on the blob, and the cross sectional area of the blob increases quadratically with laboratory time (or distance).

We compare the properties of ordinary strong magnetohydrodynamic (MHD) turbulence in a strongly magnetized medium with the recently discovered viscosity-damped regime. We focus on energy spectra, anisotropy, and intermittency. Our most surprising conclusion is that in ordinary strong MHD turbulence the velocity and magnetic fields show different high-order structure function scalings. Moreover this scaling depends on whether the intermittency is viewed in a global or local system of reference. This reconciles seemingly contradictory earlier results. On the other hand, the intermittency scaling for viscosity-damped turbulence is very different, and difficult to understand in terms of the usual phenomenological models for intermittency in turbulence. Our remaining results are in reasonable agreement with expectations. First, we find that our high resolution simulations for ordinary MHD turbulence show that the energy spectra are {\it compatible} with a Kolmogorov spectrum, while viscosity-damped turbulence shows a shallow $k^{-1}$ spectrum for the magnetic fluctuations. Second, a new numerical technique confirms that ordinary MHD turbulence exhibits Goldreich-Sridhar type anisotropy, while viscosity-damped MHD turbulence shows extremely anisotropic eddy structures. Finally, we show that many properties of incompressible turbulence for both the ordinary and viscosity-damped regimes carry over to the case of compressible turbulence.

Astrophysical objects with negligible resistivity are often threaded by large scale magnetic fields. The generation of these fields is somewhat mysterious, since a magnetic field in a perfectly conducting fluid cannot change the flux threading a fluid element, or the field topology. Classical dynamo theory evades this limit by assuming that magnetic reconnection is fast, even for vanishing resistivity, and that the large scale field can be generated by the action of kinetic helicity. Both these claims have been severely criticized, and the latter appears to conflict with strong theoretical arguments based on magnetic helicity conservation and a series of numerical simulations. Here we discuss recent efforts to explain fast magnetic reconnection through the topological effects of a weak stochastic magnetic field component. We also show how mean-field dynamo theory can be recast in a form which respects magnetic helicity conservation, and how this changes our understanding of astrophysical dynamos. Finally, we comment briefly on why an asymmetry between small scale magnetic and velocity fields is necessary for dynamo action, and how it can arise naturally.

Turbulence is the most common state of astrophysical flows. In typical astrophysical fluids, turbulence is accompanied by strong magnetic fields, which has a large impact on the dynamics of the turbulent cascade. Recently, there has been a significant breakthrough on the theory of magnetohydrodynamic (MHD) turbulence. For the first time we have a scaling model that is supported by both observations and numerical simulations. We review recent progress in studies of both incompressible and compressible turbulence. We compare Iroshnikov-Kraichnan and Goldreich-Sridhar models, and discuss scalings of Alfv\'en, slow, and fast waves. We also discuss the completely new regime of MHD turbulence that happens below the scale at which hydrodynamic turbulent motions are damped by viscosity. In the case of the partially ionized diffuse interstellar gas the viscosity is due to neutrals and truncates the turbulent cascade at $\sim$parsec scales. We show that below this scale magnetic fluctuations with a shallow spectrum persist and discuss the possibility of a resumption of the MHD cascade after ions and neutrals decouple. We discuss the implications of this new insight into MHD turbulence for cosmic ray transport, grain dynamics, etc., and how to test theoretical predictions against observations.

In astrophysical situations, e.g. in the interstellar medium (ISM), neutrals can provide viscous damping on scales much larger than the magnetic diffusion scale. Through numerical simulations, we have found that the magnetic field can have a rich structure below the dissipation cutoff scale. This implies that magnetic fields in the ISM can have structures on scales much smaller than parsec scales. Our results show that the magnetic energy contained in a wavenumber band is independent of the wavenumber and magnetic structures are intermittent and extremely anisotropic. We discuss the relation between our results and the formation of the tiny-scale atomic structure (TSAS).

We analyze 3D numerical simulations of driven incompressible magnetohydrodynamic (MHD) turbulence in a periodic box threaded by a moderately strong external magnetic field. We sum over nonlinear interactions within Fourier wavebands and find that the time scale for the energy cascade is consistent with the Goldreich-Sridhar model of strong MHD turbulence. Using higher order longitudinal structure functions we show that the turbulent motions in the plane perpendicular to the local mean magnetic field are similar to ordinary hydrodynamic turbulence while motions parallel to the field are consistent with a scaling correction which arises from the eddy anisotropy. We present the structure tensor describing velocity statistics of Alfvenic and pseudo-Alfvenic turbulence. Finally, we confirm that an imbalance of energy moving up and down magnetic field lines leads to a slow decay of turbulent motions and speculate that this imbalance is common in the interstellar medium where injection of energy is intermittent both in time and space.

We construct a magnetic helicity conserving dynamo theory which incorporates a calculated magnetic helicity current. In this model the fluid helicity plays a small role in large scale magnetic field generation. Instead, the dynamo process is dominated by a new quantity, derived from asymmetries in the second derivative of the velocity correlation function, closely related to the `twist and fold' dynamo model. The turbulent damping term is, as expected, almost unchanged. Numerical simulations with a spatially constant fluid helicity and vanishing resistivity are not expected to generate large scale fields in equipartition with the turbulent energy density. The prospects for driving a fast dynamo under these circumstances are uncertain, but if it is possible, then the field must be largely force-free. On the other hand, there is an efficient analog to the $\alpha-\Omega$ dynamo. Systems whose turbulence is driven by some anisotropic local instability in a shearing flow, like real stars and accretion disks, and some computer simulations, may successfully drive the generation of strong large scale magnetic fields, provided that $\partial_r\Omega< \partial_\theta v_z\omega_\theta>>0$. We show that this criterion is usually satisfied. Such dynamos will include a persistent, spatially coherent vertical magnetic helicity current with the same sign as $-\partial_r\Omega$, that is, positive for an accretion disk and negative for the Sun. We comment on the role of random magnetic helicity currents in storing turbulent energy in a disordered magnetic field, which will generate an equipartition, disordered field in a turbulent medium, and also a declining long wavelength tail to the power spectrum. As a result, calculations of the galactic `seed' field are largely irrelevant.

We perform direct 3-dimensional numerical simulations for magnetohydrodynamic (MHD) turbulence in a periodic box of size $2\pi$ threaded by strong uniform magnetic fields. We use a pseudo-spectral code with hyperviscosity and hyperdiffusivity to solve the incompressible MHD equations. We analyze the structure of the eddies as a function of scale. A straightforward calculation of anisotropy in wavevector space shows that the anisotropy is scale-{\it independent}. We discuss why this is {\it not} the true scaling law and how the curvature of large-scale magnetic fields affects the power spectrum and leads to the wrong conclusion. When we correct for this effect, we find that the anisotropy of eddies depends on their size: smaller eddies are more elongated than larger ones along {\it local} magnetic field lines. The results are consistent with the scaling law $\tilde{k}_{\parallel} \sim \tilde{k}_{\perp}^{2/3}$ proposed by Goldreich and Sridhar (1995, 1997). Here $\tilde{k}_{\|}$ (and $\tilde{k}_{\perp}$) are wavenumbers measured relative to the local magnetic field direction. However, we see some systematic deviations which may be a sign of limitations to the model, or our inability to fully resolve the inertial range of turbulence in our simulations.

We have tested the ability of driven turbulence to generate magnetic field structure from a weak uniform field using three dimensional numerical simulations of incompressible turbulence. We used a pseudo-spectral code with a numerical resolution of up to $144^3$ collocation points. We find that the magnetic fields are amplified through field line stretching at a rate proportional to the difference between the velocity and the magnetic field strength times a constant. Equipartition between the kinetic and magnetic energy densities occurs at a scale somewhat smaller than the kinetic energy peak. Above the equipartition scale the velocity structure is, as expected, nearly isotropic. The magnetic field structure at these scales is uncertain, but the field correlation function is very weak. At the equipartition scale the magnetic fields show only a moderate degree of anisotropy, so that the typical radius of curvature of field lines is comparable to the typical perpendicular scale for field reversal. In other words, there are few field reversals within eddies at the equipartition scale, and no fine-grained series of reversals at smaller scales. At scales below the equipartition scale, both velocity and magnetic structures are anisotropic; the eddies are stretched along the local magnetic field lines, and the magnetic energy dominates the kinetic energy on the same scale by a factor which increases at higher wavenumbers. We do not show a scale-free inertial range, but the power spectra are a function of resolution and/or the imposed viscosity and resistivity. Our results are consistent with the emergence of a scale-free inertial range at higher Reynolds numbers.

Reconnection is the process by which magnetic fields in a conducting fluid change their topology. This process is essential for understanding a wide variety of astrophysical processes, including stellar and galactic dynamos and astrophysical turbulence. To account for solar flares, solar cycles and the structure of the galactic magnetic field reconnection must be fast, propagating with a speed close to the Alfven speed. We show that the presence of a random magnetic field component substantially enhances the reconnection rate and enables fast reconnection, i.e. reconnection that does not depend on fluid resistivity. The enhancement of the reconnection rate is achieved via a combination of two effects. First of all, only small segments of magnetic field lines are subject to direct Ohmic annihilation. Thus the fraction of magnetic energy that goes directly into fluid heating goes to zero as fluid resistivity vanishes. However, the most important enhancement comes from the fact that unlike the laminar fluid case where reconnection is constrained to proceed line by line, the presence of turbulence enables many magnetic field lines to enter the reconnection zone simultaneously. A significant fraction of magnetic energy goes into MHD turbulence and this enhances reconnection rates through an increase in the field stochasticity. In this way magnetic reconnection becomes fast when field stochasticity is accounted for. As a consequence solar and galactic dynamos are also fast, i.e. do not depend on fluid resistivity.

We consider vertical heat transport in Keplerian accretion disks, including the effects of radiation, convection, and turbulent mixing driven by the Balbus-Hawley instability, in astronomical systems ranging from dwarf novae (DNe), and soft X-ray transients (SXTs), to active galactic nuclei (AGN). We propose a modified, anisotropic form of mixing-length theory, which includes radiative and turbulent damping. We also include turbulent heat transport, which acts everywhere within disks, regardless of whether or not they are stably stratified, and can move entropy in either direction. We have generated a series of vertical structure models and thermal equilibrium curves using the scaling law for the viscosity parameter $\alpha$ suggested by the exponential decay of the X-ray luminosity in SXTs. We have also included equilibrium curves for DNe using an $\alpha$ which is constant down to a small magnetic Reynolds number ($\sim 10^4$). Our models indicate that weak convection is usually eliminated by turbulent radial mixing. The substitution of turbulent heat transport for convection is more important on the unstable branches of thermal equilibrium S-curves when $\alpha$ is larger. The low temperature turnover points $\Sigma_{max}$ on the equilibrium S-curves are significantly reduced by turbulent mixing in DNe and SXT disks. However, in AGN disks the standard mixing-length theory for convection is still a useful approximation when we use the scaling law for $\alpha$, since these disks are very thin at the relevant radii. In accordance with previous work, we find that constant $\alpha$ models give almost vertical S-curves in the $\Sigma-T$ plane and consequently imply very slow, possibly oscillating, cooling waves.

The accretion-powered, X-ray pulsar 4U 1626-67 has recently shown an abrupt torque reversal accompanied by a dramatic spectral transition and a relatively small luminosity change. The time-averaged X-ray spectrum during spin-down is considerably harder than during spin-up. The observed torque reversal can be explained by an accretion flow transition triggered by a gradual change in the mass accretion rate. The sudden transition to spin-down is caused by a change in the accretion flow rotation from Keplerian to sub-Keplerian. 4U 1626-67 is estimated to be near spin equilibrium with a mass accretion rate Mdot~2x10**16 g/s, Mdot decreasing at a rate ~6x10**14 g/s/yr, and a polar surface magnetic field of ~2b_p**{-1/2} 10^**12G where b_p is the magnetic pitch. During spin-up, the Keplerian flow remains geometrically thin and cool. During spin-down, the sub-Keplerian flow becomes geometrically thick and hot. Soft photons from near the stellar surface are Compton up-scattered by the hot accretion flow during spin-down while during spin-up such scattering is unlikely due to the small scale-height and low temperature of the flow. This mechanism accounts for the observed spectral hardening and small luminosity change. The scattering occurs in a hot radially falling column of material with a scattering depth ~0.3 and a temperature ~10^9K. The X-ray luminosity at energies >5keV could be a poor indicator of the mass accretion rate. We briefly discuss the possible application of this mechanism to GX 1+4, although there are indications that this system is significantly different from other torque-reversal systems.

We examine the effect of weak, small scale magnetic field structure on the rate of reconnection in a strongly magnetized plasma. This affects the rate of reconnection by reducing the transverse scale for reconnection flows, and by allowing many independent flux reconnection events to occur simultaneously. Allowing only for the first effect and using Goldreich and Sridhar's model of strong turbulence in a magnetized plasma with negligible intermittency, we find that the lower limit for the reconnection speed is the Alfven speed times the Lundquist number to the power (-3/16). The upper limit on the reconnection speed is typically a large fraction of Alfven speed. We argue that generic reconnection in turbulent plasmas will normally occur at close to this upper limit. The fraction of magnetic energy that goes directly into electron heating scales as Lundquist number to the power (-2/5) and the thickness of the current sheet scales as the Lundquist number to the power (-3/5). A significant fraction of the magnetic energy goes into high frequency Alfven waves. We claim that the qualitative sense of these conclusions, that reconnection is fast even though current sheets are narrow, is almost independent of the local physics of reconnection and the nature of the turbulent cascade. As the consequence of this the Galactic and Solar dynamos are generically fast, i.e. do not depend on the plasma resistivity.

Dynamos in astrophysical disks are usually explained in terms of the standard alpha-omega mean field dynamo model where the local helicity generates a radial field component from an azimuthal field. The subsequent shearing of the radial field gives rise to exponentially growing dynamo modes. There are several problems with this model. The exponentiation time for the galactic dynamo is hard to calculate, but is probably uncomfortably long. Moreover, numerical simulations of magnetic fields in shearing flows indicate that the presence of a dynamo does not depend on a non-zero average helicity. However, these difficulties can be overcome by including a fluctuating helicity driven by hydrodynamic or magnetic instabilities. Unlike traditional disk dynamo models, this `incoherent' dynamo does not depend on the presence of systematic fluid helicity or any kind of vertical symmetry breaking. It will depend on geometry, in the sense that the dynamo growth rate becomes smaller for very thin disks, in agreement with constraints taken from the study of X-ray novae. In this picture the galactic dynamo will operate efficiently, but the resulting field will have a radial coherence length which is a fraction of the galactic radius.

We discuss the role of tearing instabilities in magnetic reconnection. In three dimensions this instability leads to the formation of strong Alfvenic waves that remove plasma efficiently from the reconnection layer. As a result the instability proceeds at high rates while staying close to the linear regime. Our calculations show that for a resistive fluid the reconnection speed scales as the product of the Alfven speed V_A over the magnetic Reynolds number to the power -0.3. In the limit of vanishing resistivity, tearing modes proceed at a non-zero rate, driven by the electron inertia term, giving rise to a reconnection speed V_A (c/\omega_p L_x)^{3/5}, where \omega_p is the plasma frequency and L_x is the transverse scale of the reconnection layer. Formally this solves the problem of fast reconnection, but in practice this reconnection speed is small.

Accretion disk outbursts, and their subsequent decline, offer a unique opportunity to constrain the physics of angular momentum transport in hot accretion disks. Recent work has centered on the claim by Cannizzo et al. that the exponential decay of luminosity following an outburst in black hole accretion disk systems is only consistent with a particular form for the dimensionless viscosity, $\alpha=35(c_s/r\Omega)^{3/2}$. This result can be understood in terms of a simple model of the evolution of cooling fronts in accretion disks. In particular, the cooling front speed during decline is $\sim \alpha_F c_{s,F}(c_{s,F}/r\Omega)^{n}$, where $F$ denotes the position of the cooling front, and the exact value of $n$ depends on the hot state opacity, (although generally $n\approx 1/2$). Setting this speed proportional to $r$ constrains the functional form of $\alpha$ in the hot phase of the disk, which sets it apart from previous arguments based on the relative durations of outburst and quiescence. However, it remains uncertain how well we know the exponent $n$. In addition, more work is needed to clarify the role of irradiation in these systems and its effect on the cooling front evolution.

We discuss the role of ambipolar diffusion for simple reconnection in a partially ionized gas, following the reconnection geometry of Parker and Sweet. When the recombination time is short the mobility and reconnection of the magnetic field is substantially enhanced as matter escapes from the reconnection region via ambipolar diffusion. Our analysis shows that in the interstellar medium it is the recombination rate that usually limits the rate of reconnection. We show that heating effects can reduce this rate by increasing the recombination time and raising the local ion pressure. In the colder parts of the ISM, when temperatures are or the order of 100K or less, we obtain a significant enhancement over the usual Sweet-Parker rate, but only in dense molecular clouds will the reconnection velocity exceed 0.001 times the Alfv\'en speed. The ratio of the ion orbital radius to the reconnection layer thickness is typically a few percent, except in dense molecular clouds where it can approach unity. We briefly discuss prospects for obtaining much faster reconnection speeds in astrophysical plasmas.

Accretion-powered X-ray pulsars 4U 1626-67, GX 1+4, and OAO 1657-415 have recently shown puzzling torque reversals. These reversals are characterized by short time scales, on the order of days, nearly identical spin-up and spin-down rates, and very small changes in X-ray luminosity. We propose that this phenomenon is the result of sudden dynamical changes in the accretion disks triggered by a gradual variation of mass accretion rates. These sudden torque reversals may occur at a critical accretion rate $\sim$ $10^{15}-10^{16}g~s^{-1}$ when the system makes a transition from (to) a primarily Keplerian flow to (from) a substantially sub-Keplerian, radial advective flow in the inner disk. For systems near spin equilibrium, the spin-up torques in the Keplerian state are slightly larger than the spin-down torques in the advective state, in agreement with observation. The abrupt reversals could be a signature of pulsar systems near spin equilibrium with the mass accretion rates modulated on a time scale of a year or longer near the critical accretion rate. It is interesting that cataclysmic variables and black hole soft X-ray transients change their X-ray emission properties at accretion rates similar to the pulsars' critical rate. We speculate that the dynamical change in pulsar systems shares a common physical origin with white dwarf and black hole accretion disk systems.

Recent work has shown that the speed of the cooling front in soft X-ray transients may be an important clue in understanding the nature of accretion disk viscosity. In a previous paper (Vishniac and Wheeler 1996) we derived the scaling law for the cooling front speed. Here we derive a similarity solution for the hot inner part of disks undergoing cooling. This solution is exact in the limit of a thin disk, power law opacities, and a minimum hot state column density which is an infinitesimal fraction of the maximum cold state density. For a disk of finite thickness the largest error is in the ratio of the mass flow across the cooling front to the mass flow at small radii. Comparison to the numerical simulations of Cannizzo et al. (1995) inidcates that the errors in the other parameters do not exceed $(c_{sF}/r_F\Omega_F)^q$, that is, the ratio of the sound speed at the disk midplane to its orbital velocity, evaluated at the cooling front, to the qth power. Here $q\approx 1/2$. Its precise value is determined by the relevant hot state opacity law and the functional form of the dimensionless viscosity.

The gravitational instability of cosmological pancakes composed of collisionless dark matter in an Einstein-de Sitter universe is investigated numerically to demonstrate that pancakes are unstable with respect to fragmentation and the formation of filaments. A ``pancake'' is defined here as the nonlinear outcome of the growth of a 1D, sinusoidal, plane-wave, adiabatic density perturbation. We have used high resolution, 2D, N-body simulations by the Particle-Mesh (PM) method to study the response of pancakes to perturbation by either symmetric (density) or antisymmetric (bending or rippling) modes, with corresponding wavevectors k_s and k_a transverse to the wavevector k_p of the unperturbed pancake plane-wave. We consider dark matter which is initially ``cold'' (i.e. with no random thermal velocity in the initial conditions). We also investigate the effect of a finite, random, isotropic, initial velocity dispersion (i.e. initial thermal velocity) on the fate of pancake collapse and instability. Pancakes are shown to be gravitationally unstable with respect to all perturbations of wavelength l<l_p (where l_p= 2pi/k_p). These results are in contradiction with the expectations of an approximate, thin-sheet energy argument.

We consider the radial buoyancy of vertical magnetic field lines in radiation and gas pressure dominated accretion disks. We find that in addition to radial drift driven by turbulent diffusion and biased by the global field geometry, there are buoyancy effects which tend to move magnetic flux outward. In gas pressure dominated disks the poloidal magnetic field will move outward at a rate comparable to its inward advection. On the other hand, in a radiation pressure dominated disk the poloidal magnetic field will usually move outward faster than it is advected inward. This implies that the fields in disks in active galactic nuclei are generated at small radii by an internal disk dynamo. This conclusion can be avoided if the external field imposes a supersonic Alfv\'en speed within the disk without giving rise to interchange instabilities. In any case we note that variations in the mass transfer rate will lead directly to a modulation of the nonthermal emission from the disk system.

We examine the speed of inward traveling cooling fronts in accretion disks. We show that their speed is determined by the rarefaction wave that precedes them and is approximately $\alpha_F c_{F} (H/r)^q$, where $\alpha_F$ is the dimensionless viscosity, $c_{F}$ is the sound speed, $r$ is the radial coordinate, $H$ is the disk thickness, and all quantities are evaluated at the cooling front. The scaling exponent $q$ lies in the interval $[0,1]$, depending on the slope of the $(T,\Sigma)$ relation in the hot state. For a Kramer's law opacity and $\alpha\propto (H/r)^n$, where $n$ is of order unity, we find that $q\sim 1/2$. This supports the numerical work of Cannizzo, Chen and Livio (1995) and their conclusion that $n\approx3/2$ is necessary to reproduce the exponential decay of luminosity in black hole X-ray binary systems. Our results are insensitive to the structure of the disk outside of the radius where rapid cooling sets in. In particular, the width of the rapid cooling zone is a consequence of the cooling front speed rather than its cause. We conclude that the exponential luminosity decay of cooling disks is probably compatible with the wave-driven dynamo model. It is not compatible with models with separate, constant values of $\alpha$ for the hot and cold states.

We use the mean-field dynamo equations to show that an incoherent alpha effect in mirror-symmetric turbulence in a shearing flow can generate a large scale, coherent magnetic field. We illustrate this effect with simulations of a few simple systems. In accretion disks, this process can lead to axisymmetric magnetic domains whose radial and vertical dimensions will be comparable to the disk height. This process may be responsible for observations of dynamo activity seen in simulations of dynamo-generated turbulence involving, for example, the Balbus-Hawley instability. In this case the magnetic field strength will saturate at $\sim (h/r)^2$ times the ambient pressure in real accretion disks. The resultant dimensionless viscosity will be of the same order. In numerical simulations the azimuthal extent of the simulated annulus should be substituted for $r$. We compare the predictions of this model to numerical simulations previously reported by Brandenburg et al. (1995). In a radiation pressure dominated environment this estimate for viscosity should be reduced by a factor of $(P_{gas}/P_{radiation})^6$ due to magnetic buoyancy.

Tidally distorted accretion disks in binary star systems are subject to a local hydrodynamic instability which excites $m=1$ internal waves. This instability is three dimensional and approximately incompressible. We study the global aspects of this local instability using equations derived under the shearing sheet approximation, where the effects of the azimuthal variation along distorted orbital trajectories are included in source terms which oscillate with local orbital phase. Linear analyses show that the excitation of the instability is essentially local, i.e. insensitive to radial boundary conditions. The region of rapid growth feeds waves into the region of slow or negligible growth, allowing the instability to become global. The global growth rate depends the maximum local growth rate, the size of the rapid growth region, and the local group velocity. We present an empirical expression for the global growth rate. We note that the local nature of the instability allows the excitation of waves with $m\ne 1$ when the local growth rate is large.

Recent work on the structure of magnetic fields in a turbulent medium gives predictions for the properties of the magnetic flux tubes as a function of the Mach number and scale of the turbulence, and the resistivity and viscosity of the fluid. Here I discuss the implications of this work for accretion disks. I show that although accretion disk flux tubes are usually almost completely evacuated, they are nevertheless less buoyant than previous estimates have suggested. I also note that vertical magnetic flux tends to be ejected from the outer edge of accretion disks, so evidence for continued magnetic activity in such systems should be interpreted as supporting the existence of dynamo activity.

Accretion disks in close binary systems are subject to a tidally driven parametric instability which leads to the growth of internal waves near the outer edges of such disks (Goodman 1993). These waves are important in understanding the torque exerted on a disk by tidal forces and may play a role in the structure of the disk at small radii. Here we calculate the growth rate of this instability, including the effects of vertical structure and fluid compressibility. We find growth rates which are only slightly different from Goodman's original results, except that near the vertical resonance radius the growth rate can have an extremely broad and strong peak when the disk is stably stratified in the vertical direction. Higher order modes, in the sense of increasing number of vertical nodes, have similar growth rates. Our results differ from a previous calculation along these lines by Lubow et al. (1993). The difference is mostly due to their neglect of radial gradients in the tidally distorted streamlines.

We suggest a new model for the structure of a magnetic field embedded high $\beta$ turbulent plasma, based on the popular notion that the magnetic field will tend to separate into individual flux tubes. We point out that interactions between the flux tubes will be dominated by coherent effects stemming from the turbulent wakes created as the fluid streams by the flux tubes. Balancing the attraction caused by shielding effects with turbulent diffusion we find that flux tubes have typical radii comparable to the local Mach number squared times the large scale eddy length, are arranged in a one dimensional fractal pattern, have a radius of curvature comparable to the largest scale eddies in the turbulence, and have an internal magnetic pressure comparable to the ambient pressure. When the average magnetic energy density is much less than the turbulent energy density the radius, internal magnetic field and curvature scale of the flux tubes will be smaller than these estimates. Realistic resistivity does not alter the macroscopic properties of the fluid or the large scale magnetic field. In either case we show that the Sweet-Parker reconnection rate is much faster than an eddy turnover time. Realistic stellar plasmas are expected to either be in the ideal limit (e.g. the solar photosphere) or the resistive limit (most of the solar convection zone). All current numerical simulations of three dimensional MHD turbulence are in the viscous regime and are inapplicable to stars or accretion disks.

I review recent work on the radial transport of angular momentum in ionized, Keplerian accretion disks. Proposed mechanisms include hydrodynamic and MHD local instabilities and long range effects mediated by wave transport. The most promising models incorporate the Velikhov-Chandrasekhar instability, caused by an instability of the magnetic field embedded in a differentially rotating disk. This has the important feature that the induced turbulent motions necessarily transport angular momentum outward. By contrast, convective modes may transport angular momentum in either direction. Combining the magnetic field instability with an $\alpha-\Omega$ dynamo driven by internal waves leads to a model in which the dimensionless viscosity scales as $(H/r)^{4/3}$. However, this model has a phenomenology which is quite different from the $\alpha$ disk model. For example, an active disk implies some source of excitation for the internal waves. In binary systems with a mass ratio of order unity the most likely exciting mechanism is a parametric instability due to tidal forces. This implies that in systems where the accretion stream is intermittent, like MV Lyrae or TT Ari, epochs when the mass flow is absent or very small will be epochs in which the disk shrinks and becomes relatively inactive and dark. This model also implies that forced vertical mixing is important, even in convectively stable disks. I discuss various observational tests of this model and the focus of current theoretical work.

Recent work on the transport of angular momentum in accretion disks suggests that the Velikhov-Chandrasekhar instability, in which a large scale magnetic field generates small scale eddys in a shearing environment, may be ultimately responsible for this process. Although there is considerable controversy about the origin and maintenance of this field in accretion disks, it turns out that it is possible to argue, quite generally, using scaling arguments, that this process is sensitive to the total pressure in an AGN disk, rather than the pressure contributed by gas alone. We conclude that the resolution of the conceptual difficulties implied by the presence of strong thermal and viscous instabilities in radiation pressure and electron scattering dominated does not lie in models that couple the total dissipation rate to the gas pressure alone, or to some weighted mean of the gas and radiation pressures.

(substantial changes to section 3.2, otherwise minor) We present an analysis of the hydrodynamic stability of a cold slab bounded by two accretion shocks. Previous numerical work has shown that when the Mach number of the shock is large the slab is unstable. Here we show that to linear order both the bending and breathing modes of such a slab are stable. However, nonlinear effects will tend to soften the restoring forces for bending modes, and when the slab displacement is comparable to its thickness this gives rise to a nonlinear instability. The growth rate of the instability, above this threshold but for small bending angles, is $\sim c_sk (k\eta)^{1/2}$, where $\eta$ is the slab displacement. When the bending angle is large the slab will contain a local vorticity comparable to $c_s/L$, where $L$ is the slab thickness. We discuss the implications of this work for gravitational instabilities of slabs. Finally, we examine the cases of a decelerating slab bounded by a single shock and a stationary slab bounded on one side by thermal pressure. The latter case is stable, but appears to be a special case. The former case is subject to a nonlinear overstability driven by deceleration effects. We conclude that shock bounded slabs with a high density compression ratio generically produce substructure with a strong local shear, a bulk velocity dispersion like the sound speed in the cold layer and a characteristic scale comparable to the slab thickness.

We review mechanisms for the transport of angular momentum in accretion disks that are low mass, in the sense that the gravitational forces produced by the material in these disks has a negligible effect on disk dynamics. There is no established consensus on how this transport takes place. We note that for phenomenological reasons the traditional $\alpha$ model is probably not a good description of real disks. We discuss models in which angular momentum transport is driven by shocks, and by magnetic field instabilities. The latter is more promising, but requires a dynamo. We note that the direction of angular momentum transport due to convection in a conducting disk is not known as competing mechanisms are at work. We briefly discuss a number of possible dynamo mechanisms, and their problems. We then give a detailed exposition of the internal wave driven dynamo model, in which internal waves excited at large radii drive an $\alpha-\Omega$ dynamo.

The galactic magnetic field is commonly supposed to be due to a dynamo acting on some large scale seed field. A major difficulty with this idea is that estimates of reasonable seed field strengths tend to be quite low, on the order of $\sim10^{-20}$ gauss. Here we examine the contribution due to the flux entrained in winds from protostars formed in the first dynamo e-folding time of a galaxy's existence. Using a minimal estimate of a protostellar magnetic field we find that if each protostar ejects a single current ring, sufficient to maintain flux freezing in the wind, than the large scale average dipole field from all such current rings will be at least 5 orders of magnitude larger than previous seed field estimates. Allowing for a reasonable amount of magnetic activity in protostars during an extended period of mass loss increases this to a dipole seed field of $\sim10^{-12}$ gauss. For the purposes of producing a seed field it is irrelevant whether or not this initial injection of flux takes place in a newly formed galactic disk, or in star forming proto-galactic clouds. The compression of this dipole field into a thin disk will lead to a large scale $B_r\sim 10^{-10.5}$ gauss. Initially, field strengths on smaller scales will be larger, but nowhere near current levels.