
Magnetohydrodynamic turbulence and magnetic reconnection are ubiquitous in
astrophysical environments. In most situations, these processes do not occur in
isolation, but interact with each other. This renders a comprehensive theory of
these processes highly challenging. Here, we propose a theory of
magnetohydrodynamic turbulence driven at large scale that selfconsistently
accounts for the mutual interplay with magnetic reconnection occurring at
smaller scales. Magnetic reconnection produces plasmoids that grow from
turbulencegenerated noise and eventually disrupt the sheetlike structures in
which they are born. The disruption of these structures leads to a modification
of the turbulent energy cascade, which, in turn, exerts a feedback effect on
the plasmoid formation via the turbulencegenerated noise. The energy spectrum
in this plasmoidmediated range steepens relative to the standard inertial
range and does not follow a simple power law. As a result of the complex
interplay between turbulence and reconnection, we also find that the length
scale which marks the beginning of the plasmoidmediated range and the
dissipation length scale do not obey true power laws. The transitional magnetic
Reynolds number above which the plasmoid formation becomes statistically
significant enough to affect the turbulent cascade is fairly modest, implying
that plasmoids are expected to modify the turbulent path to dissipation in many
astrophysical systems.

We present optical continuum lags for two Seyfert 1 galaxies, MCG+0811011
and NGC 2617, using monitoring data from a reverberation mapping campaign
carried out in 2014. Our light curves span the ugriz filters over four months,
with median cadences of 1.0 and 0.6 days for MCG+0811011 and NGC\,2617,
respectively, combined with roughly daily Xray and nearUV data from Swift for
NGC 2617. We find lags consistent with geometrically thin accretiondisk models
that predict a lagwavelength relation of $\tau \propto \lambda^{4/3}$.
However, the observed lags are larger than predictions based on standard
thindisk theory by factors of 3.3 for MCG+0811011 and 2.3 for NGC\,2617.
These differences can be explained if the mass accretion rates are larger than
inferred from the optical luminosity by a factor of 4.3 in MCG+0811011 and a
factor of 1.3 in NGC\,2617, although uncertainty in the SMBH masses determines
the significance of this result. While the Xray variability in NGC\,2617
precedes the UV/optical variability, the long 2.6 day lag is problematic for
coronal reprocessing models.

We revisit Parker's conjecture of current singularity formation in 3D
linetied plasmas using a recently developed numerical method, variational
integration for ideal magnetohydrodynamics in Lagrangian labeling. With the
frozenin equation builtin, the method is free of artificial reconnection, and
hence it is arguably an optimal tool for studying current singularity
formation. Using this method, the formation of current singularity has
previously been confirmed in the HahmKulsrudTaylor problem in 2D. In this
paper, we extend this problem to 3D linetied geometry. The linear solution,
which is singular in 2D, is found to be smooth for arbitrary system length.
However, with finite amplitude, the linear solution can become pathological
when the system is sufficiently long. The nonlinear solutions turn out to be
smooth for short systems. Nonetheless, the scaling of peak current density
versus system length suggests that the nonlinear solution may become singular
at finite length. With the results in hand, we can neither confirm nor rule out
this possibility conclusively, since we cannot obtain solutions with system
length near the extrapolated critical value.

The scaling of plasmoid instability maximum linear growth rate with respect
to Lundquist number $S$ in a SweetParker current sheet, $\gamma_{max}\sim
S^{1/4}$, indicates that at high $S$, the current sheet will break apart before
it approaches the SweetParker width. Therefore, a proper description for the
onset of the plasmoid instability must incorporate the evolving process of the
current sheet. We carry out a series of twodimensional simulations and develop
diagnostics to separate fluctuations from an evolving background. It is found
that the fluctuation amplitude starts to grow only when the linear growth rate
is sufficiently large ($\gamma_{max}\tau_{A}>O(1)$) to overcome convective
losses. The linear growth rate continues to rise until the sizes of plasmoids
become comparable to the inner layer width of the tearing mode. At this point
the current sheet is disrupted and the instability enters the early nonlinear
regime. The growth rate suddenly decreases, but the fluctuation amplitude
continues to grow until it reaches nonlinear saturation. We identify important
time scales of the instability development, as well as scalings for the linear
growth rate, current sheet width, and dominant wavenumber at current sheet
disruption. These scalings depend on not only the Lundquist number, but also
the initial noise amplitude. A phenomenological model that reproduces scalings
from simulation results is proposed. The model incorporates the effect of
reconnection outflow, which is crucial for yielding a critical Lundquist number
$S_{c}$ below which disruption does not occur. The critical Lundquist number
$S_{c}$ is not a constant value but has a weak dependence on the noise
amplitude.

The surface current on the plasmavacuum interface during a disruption event
involving kink instability can play an important role in driving current into
the vacuum vessel. However, there have been disagreements over the nature or
even the sign of the surface current in recent theoretical calculations based
on idealized stepfunction background plasma profiles. We revisit such
calculations by replacing stepfunction profiles with more realistic profiles
characterized by strong but finite gradient along the radial direction. It is
shown that the resulting surface current is no longer a deltafunction current
density, but a finite and smooth current density profile with internal
structure, concentrated within the region with strong plasma pressure gradient.
Moreover, this current density profile has peaks of both signs, unlike the
deltafunction case with a sign opposite to, or the same as the plasma current.
We show analytically and numerically that such current density can be separated
into two parts, with one of them, called the convective current density,
describing the transport of the background plasma density by the displacement,
and the other part that remains, called the residual current density. It is
argued that consideration of both types of current density is important and can
resolve past controversies.

A general theory of the onset and development of the plasmoid instability is
formulated by means of a principle of least time. The scaling relations for the
final aspect ratio, transition time to rapid onset, growth rate, and number of
plasmoids are derived, and shown to depend on the initial perturbation
amplitude $\left({\hat w}_0\right)$, the characteristic rate of current sheet
evolution $\left(1/\tau\right)$, and the Lundquist number $\left(S\right)$.
They are not simple power laws, and are proportional to $S^{\alpha}
\tau^{\beta} \left[\ln f(S,\tau,{\hat w}_0)\right]^\sigma$. The detailed
dynamics of the instability is also elucidated, and shown to comprise of a
period of quiescence followed by sudden growth over a short time scale.

A number of studies have considered how the rate of magnetic reconnection
scales in large and weakly collisional systems by the modelling of long
reconnecting current sheets. However, this setup neglects both the formation
of the current sheet and the coupling between the diffusion region and a larger
system that supplies the magnetic flux. Recent studies of magnetic island
merging, which naturally include these features, have found that ion kinetic
physics is crucial to describe the reconnection rate and global evolution of
such systems. In this paper, the effect of a guide field on reconnection during
island merging is considered. In contrast to the earlier current sheet studies,
we identify a limited range of guide fields for which the reconnection rate,
outflow velocity, and pileup magnetic field increase in magnitude as the guide
field increases. The HallMHD fluid model is found to reproduce kinetic
reconnection rates only for a sufficiently strong guide field, for which ion
inertia breaks the frozenin condition and the outflow becomes Alfvenic in the
kinetic system. The merging of large islands occurs on a longer timescale in
the zero guide field limit, which may in part be due to a mirrorlike
instability that occurs upstream of the reconnection region.

Magnetic reconnection is a fundamental plasma process involving an exchange
of magnetic energy to plasma kinetic energy through changes in the magnetic
field topology. In many astrophysical plasmas magnetic reconnection plays a key
role in the release of large amounts of energy \cite{hoshino1}, although making
direct measurements is challenging in the case of highenergy astrophysical
systems such as pulsar wind emissions \cite{lyubarsky1}, gammaray bursts
\cite{thompson1}, and jets from active galactic nuclei \cite{liu1}. Therefore,
laboratory studies of magnetic reconnection provide an important platform for
testing theories and characterising different regimes. Here we present
experimental measurements as well as numerical modeling of relativistic
magnetic reconnection driven by shortpulse, highintensity lasers that produce
relativistic plasma along with extremely strong magnetic fields. Evidence of
magnetic reconnection was identified by the plasma's Xray emission patterns,
changes to the electron energy spectrum, and by measuring the time over which
reconnection occurs. Accessing these relativistic conditions in the laboratory
allows for further investigation that may provide insight into unresolved areas
in space and astrophysics.

We present the first results from an optical reverberation mapping campaign
executed in 2014, targeting the active galactic nuclei (AGN) MCG+0811011, NGC
2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting
observational properties, including a "changing look" AGN and a broadline
radio galaxy. Based on continuumH$\beta$ lags, we measure black hole masses
for all five targets. We also obtain H$\gamma$ and He{\sc ii}\,$\lambda 4686$
lags for all objects except 3C 382. The He{\sc ii}\,$\lambda 4686$ lags
indicate radial stratification of the BLR, and the masses derived from
different emission lines are in general agreement. The relative responsivities
of these lines are also in qualitative agreement with photoionization models.
These spectra have extremely high signaltonoise ratios (100300 per pixel)
and there are excellent prospects for obtaining velocityresolved reverberation
signatures.

Although the ShakuraSunyaev $\alpha$ viscosity prescription has been highly
successful in characterizing myriad astrophysical environments, it has proven
to be partly inadequate in modelling turbulent stresses driven by the MRI.
Hence, we adopt the approach employed by \citet{GIO03}, but in the context of
Hall magnetohydrodynamics (MHD), to study MRI turbulence. We utilize the exact
evolution equations for the stresses, and the nonlinear terms are closed
through the invocation of dimensional analysis and physical considerations. We
demonstrate that the inclusion of the Hall term leads to nontrivial results,
including the modification of the Reynolds and Maxwell stresses, as well as the
(asymptotic) nonequipartition between the kinetic and magnetic energies; the
latter issue is also addressed via the analysis of nonlinear waves. The
asymptotic ratio of the kinetic and magnetic energies is shown to be
\emph{independent} of the choice of initial conditions, but it is governed by
the \emph{Hall parameter}. We contrast our model with the Kazantsev
prescription from smallscale dynamo theory, and the Hall term does not
contribute in the latter approach, illustrating the limitations of the
Kazantsev formalism. We indicate potential astrophysical applications of our
model, including the solar wind where a lack of equipartition has been
observed.

Recently a variational integrator for ideal magnetohydrodynamics in
Lagrangian labeling has been developed. Its builtin frozenin equation makes
it optimal for studying current sheet formation. We use this scheme to study
the HahmKulsrudTaylor problem, which considers the response of a 2D plasma
magnetized by a sheared field under sinusoidal boundary forcing. We obtain an
equilibrium solution that preserves the magnetic topology of the initial field
exactly, with a fluid mapping that is nondifferentiable. Unlike previous
studies that examine the current density output, we identify a singular current
sheet from the fluid mapping. These results are benchmarked with a constrained
GradShafranov solver. The same signature of current singularity can be found
in other cases with more complex magnetic topologies.

It has been established that the SweetParker current layer in high Lundquist
number reconnection is unstable to the superAlfv\'enic plasmoid instability.
Past twodimensional magnetohydrodynamic simulations have demonstrated that the
plasmoid instability leads to a new regime where the SweetParker current layer
changes into a chain of plasmoids connected by secondary current sheets, and
the averaged reconnection rate becomes nearly independent of the Lundquist
number. In this work, threedimensional simulation with a guide field shows
that the additional degree of freedom allows plasmoid instabilities to grow at
oblique angles, which interact and lead to selfgenerated turbulent
reconnection. The averaged reconnection rate in the selfgenerated turbulent
state is of the order of a hundredth of the characteristic Alfv\'en speed,
which is similar to the twodimensional result but is an order of magnitude
lower than the fastest reconnection rate reported in recent studies of
externally driven threedimensional turbulent reconnection. Kinematic and
magnetic energy fluctuations both form elongated eddies along the direction of
local magnetic field, which is a signature of anisotropic magnetohydrodynamic
turbulence. Both energy fluctuations satisfy powerlaw spectra in the inertial
range, where the magnetic energy spectral index is in the range from $2.3$ to
$2.1$, while the kinetic energy spectral index is slightly steeper, in the
range from $2.5$ to $2.3$. The anisotropy of turbulence eddies is found to be
nearly scaleindependent, in contrast with the prediction of the
GoldreichSridhar theory for anisotropic turbulence in a homogeneous plasma
permeated by a uniform magnetic field.

In the interaction of highpower laser beams with solid density plasma there
are a number of mechanisms that generate strong magnetic fields. Such fields
subsequently inhibit or redirect electron flows, but can themselves be advected
by heat fluxes, resulting in complex interplay between thermal transport and
magnetic fields.We show that for heating by multiple laser spots reconnection
of magnetic field lines can occur, mediated by these heat fluxes, using a fully
implicit 2D VlasovFokkerPlanck code. Under such conditions, the reconnection
rate is dictated by heat flows rather than Alfv\`enic flows. We find that this
mechanism is only relevant in a high $\beta$ plasma. However, the Hall
parameter $\omega_c \tau_{ei}$ can be large so that thermal transport is
strongly modified by these magnetic fields, which can impact longer time scale
temperature homogeneity and ion dynamics in the system.

Newcomb's Lagrangian for ideal magnetohydrodynamics (MHD) in Lagrangian
labeling is discretized using discrete exterior calculus. Variational
integrators for ideal MHD are derived thereafter. Besides being symplectic and
momentumpreserving, the schemes inherit builtin advection equations from
Newcomb's formulation, and therefore avoid solving them and the accompanying
error and dissipation. We implement the method in 2D and show that numerical
reconnection does not take place when singular current sheets are present. We
then apply it to studying the dynamics of the ideal coalescence instability
with multiple islands. The relaxed equilibrium state with embedded current
sheets is obtained numerically.

We introduce an extensible multifluid moment model in the context of
collisionless magnetic reconnection. This model evolves full Maxwell equations,
and simultaneously moments of the VlasovMaxwell equation for each species in
the plasma. Effects like electron inertia and pressure gradient are
selfconsistently embedded in the resulting multifluid moment equations,
without the need to explicitly solving a generalized Ohms's law. Two limits of
the multifluid moment model are discussed, namely, the fivemoment limit that
evolves a scalar pressures for each species, and the tenmoment limit that
evolves the full anisotropic, nongyrotropic pressure tensor for each species.
We first demonstrate, analytically and numerically, that the fivemoment model
reduces to the widely used Hall Magnetohydrodynamics (Hall MHD) model under the
assumptions of vanishing electron inertia, infinite speed of light, and
quasineutrality. Then, we compare tenmoment and fully kinetic
ParticleInCell (PIC) simulations of a large scale Harris sheet reconnection
problem, where the tenmoment equations are closed with a local linear
collisionless approximation for the heat flux. The tenmoment simulation gives
reasonable agreement with the PIC results regarding the structures and
magnitudes of the electron flows, the polarities and magnitudes of elements of
the electron pressure tensor, and the decomposition of the generalized Ohm's
law. Possible ways to improve the simple local closure towards a nonlocal fully
threedimensional closure are also discussed.

Magnetic fields without a direction of continuous symmetry have the generic
feature that neighboring field lines exponentiate away from each other and
become stochastic, hence the ideal constraint of preserving magnetic field line
connectivity becomes exponentially sensitive to small deviations from ideal
Ohm's law. The idea of breaking field line connectivity by stochasticity as a
mechanism for fast reconnection is tested with numerical simulations based on
reduced magnetohydrodynamics equations with a strong guide field linetied to
two perfectly conducting end plates. Starting from an ideally stable forcefree
equilibrium, the system is allowed to undergo resistive relaxation. Two
distinct phases are found in the process of resistive relaxation. During the
quasistatic phase, rapid change of field line connectivity and strong induced
flow are found in regions of high field line exponentiation. However, although
the field line connectivity of individual field lines can change rapidly, the
overall pattern of field line mapping appears to deform gradually. From this
perspective, field line exponentiation appears to cause enhanced diffusion
rather than reconnection. In some cases, resistive quasistatic evolution can
cause the ideally stable initial equilibrium to cross a stability threshold,
leading to formation of intense current filaments and rapid change of field
line mapping into a qualitatively different pattern. It is in this onset phase
that the change of field line connectivity is more appropriately designated as
magnetic reconnection. Our results show that rapid change of field line
connectivity appears to be a necessary, but not a sufficient condition for fast
reconnection.

Supraarcade downflows (hereafter referred to as SADs) are lowemission,
elongated, fingerlike features usually observed in activeregion coronae above
posteruption flare arcades. Observations exhibit downward moving SADs
intertwined with bright upward moving spikes. Whereas SADs are dark voids,
spikes are brighter, denser structures. Although SADs have been observed for
decades, the mechanism of formation of SADs remains an open issue. In our
threedimensional resistive magnetohydrodynamic simulations, we demonstrate
that secondary RayleighTaylor type instabilities develop in the downstream
region of a reconnecting current sheet. The instability results in the
formation of lowdensity coherent structures that resemble SADs, and
highdensity structures that appear to be spikelike. Comparison between the
simulation results and observations suggests that secondary RayleighTaylor
type instabilities in the exhaust of reconnecting current sheets provide a
plausible mechanism for observed SADs and spikes.

The constraint imposed by magnetic helicity conservation on the alpha effect
is considered for both magnetically and flow dominated selforganizing plasmas.
Direct numerical simulations are presented for a dominant contribution to the
alpha effect, which can be cast in the functional form of a total divergence of
an averaged helicity flux, called the helicityfluxdriven alpha ( H$\alpha$)
effect. Direct numerical simulations of the H$\alpha$ effect are prese nted for
two examplesthe magnetically dominated toroidal plasma unstable to tearing
modes, and the flowdominated accretion disk.

Filamentation due to the growth of a Weibeltype instability was observed in
the interaction of a pair of counterstreaming, ablativelydriven plasma flows,
in a supersonic, collisionless regime relevant to astrophysical collisionless
shocks. The flows were created by irradiating a pair of opposing plastic (CH)
foils with 1.8 kJ, 2ns laser pulses on the OMEGA EP laser system. Ultrafast
laserdriven proton radiography was used to image the Weibelgenerated
electromagnetic fields. The experimental observations are in good agreement
with the analytical theory of the Weibel instability and with particleincell
simulations.

Magnetic reconnection mediated by the hyperresistive plasmoid instability is
studied with both linear analysis and nonlinear simulations. The linear growth
rate is found to scale as $S_{H}^{1/6}$ with respect to the hyperresistive
Lundquist number $S_{H}\equiv L^{3}V_{A}/\eta_{H}$, where $L$ is the system
size, $V_{A}$ is the Alfv\'en velocity, and $\eta_{H}$ is the
hyperresistivity. In the nonlinear regime, reconnection rate becomes nearly
independent of $S_{H}$, the number of plasmoids scales as $S_{H}^{1/2}$, and
the secondary current sheet length and width both scale as $S_{H}^{1/2}$.
These scalings are consistent with a heuristic argument assuming secondary
current sheets are close to marginal stability. The distribution of plasmoids
as a function of the enclosed flux $\psi$ is found to obey a $\psi^{1}$ power
law over an extended range, followed by a rapid fall off for large plasmoids.
These results are compared with those from resistive magnetohydrodynamic
studies.

The HahmKulsrud (HK) [T. S. Hahm and R. M. Kulsrud, Phys. Fluids {\bf 28},
2412 (1985)] solutions for a magnetically sheared plasma slab driven by a
resonant periodic boundary perturbation illustrate fully shielded (current
sheet) and fully reconnected (magnetic island) responses. On the global scale,
reconnection involves solving a magnetohydrodynamic (MHD) equilibrium problem.
In systems with a continuous symmetry such MHD equilibria are typically found
by solving the GradShafranov equation, and in slab geometry the elliptic
operator in this equation is the 2D Laplacian. Thus, assuming appropriate
pressure and poloidal current profiles, a conformal mapping method can be used
to transform one solution into another with different boundary conditions,
giving a continuous sequence of solutions in the form of partially reconnected
magnetic islands (plasmoids) separated by Syrovatsky current sheets. The two HK
solutions appear as special cases.

Our understanding of magnetic reconnection in resistive magnetohydrodynamics
has gone through a fundamental change in recent years. The conventional wisdom
is that magnetic reconnection mediated by resistivity is slow in laminar high
Lundquist ($S$) plasmas, constrained by the scaling of the reconnection rate
predicted by SweetParker theory. However, recent studies have shown that when
$S$ exceeds a critical value $\sim10^{4}$, the SweetParker current sheet is
unstable to a superAlfv\'enic plasmoid instability, with a linear growth rate
that scales as $S^{1/4}$. In the fully developed statistical steady state of
twodimensional resistive magnetohydrodynamic simulations, the normalized
average reconnection rate is approximately 0.01, nearly independent of $S$, and
the distribution function $f(\psi)$ of plasmoid magnetic flux $\psi$ follows a
power law $f(\psi)\sim\psi^{1}$. When Hall effects are included, the plasmoid
instability may trigger onset of Hall reconnection even when the conventional
criterion for onset is not satisfied. The rich variety of possible reconnection
dynamics is organized in the framework of a phase diagram.

In a recent numerical study [Ng et al., Astrophys. J. {\bf 747}, 109, 2012],
with a threedimensional model of coronal heating using reduced
magnetohydrodynamics (RMHD), we have obtained scaling results of heating rate
versus Lundquist number based on a series of runs in which random photospheric
motions are imposed for hundreds to thousands of \al time in order to obtain
converged statistical values. The heating rate found in these simulations
saturate to a level that is independent of the Lundquist number. This scaling
result was also supported by an analysis with the assumption of the
SweetParker scaling of the current sheets, as well as how the width, length
and number of current sheets scale with Lundquist number. In order to test
these assumptions, we have implemented an automated routine to analyze
thousands of current sheets in these simulations and return statistical
scalings for these quantities. It is found that the SweetParker scaling is
justified. However, some discrepancies are also found and require further
study.

The distribution function $f(\psi)$ of magnetic flux $\psi$ in plasmoids
formed in highLundquistnumber current sheets is studied by means of an
analytic phenomenological model and direct numerical simulations. The
distribution function is shown to follow a power law $f(\psi)\sim\psi^{1}$,
which differs from other recent theoretical predictions. Physical explanations
are given for the discrepant predictions of other theoretical models.

Parker's model is one of the most discussed mechanisms for coronal heating
and has generated much debate. We have recently obtained new scaling results in
a twodimensional (2D) version of this problem suggesting that the heating rate
becomes independent of resistivity in a statistical steady state [Ng and
Bhattacharjee, Astrophys. J., 675, 899 (2008)]. Our numerical work has now been
extended to 3D by means of largescale numerical simulations. Random
photospheric footpoint motion is applied for a time much longer than the
correlation time of the motion to obtain converged average coronal heating
rates. Simulations are done for different values of the Lundquist number to
determine scaling. In the highLundquist number limit, the coronal heating rate
obtained so far is consistent with a trend that is independent of the Lundquist
number, as predicted by previous analysis as well as 2D simulations. In the
same limit the average magnetic energy built up by the random footpoint motion
tends to have a much weaker dependence on the Lundquist number than that in the
2D simulations, due to the formation of strong current layers and subsequent
disruption when the equilibrium becomes unstable. We will present scaling
analysis showing that when the dissipation time is comparable or larger than
the correlation time of the random footpoint motion, the heating rate tends to
become independent of Lundquist number, and that the magnetic energy production
is also reduced significantly.