Source author record

William Daughton

William Daughton appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

Magnetic reconnection in the era of exascale computing and multiscale experiments

Astrophysical plasmas have the remarkable ability to preserve magnetic topology, which inevitably gives rise to the accumulation of magnetic energy within stressed regions including current sheets. This stored energy is often released explosively through the process of magnetic reconnection, which produces a reconfiguration of the magnetic field, along with high-speed flows, thermal heating, and nonthermal particle acceleration. Either collisional or kinetic dissipation mechanisms are required to overcome the topological constraints, both of which have been predicted by theory and validated with in situ spacecraft observations or laboratory experiments. However, major challenges remain in understanding magnetic reconnection in large systems, such as the solar corona, where the collisionality is weak and the kinetic scales are vanishingly small in comparison to macroscopic scales. The plasmoid instability or formation of multiple plasmoids in long reconnecting current sheets is one possible multiscale solution for bridging this vast range of scales, and new laboratory experiments are poised to study these regimes. In conjunction with these efforts, we anticipate that the coming era of exascale computing, together with the next generation of observational capabilities, will enable new progress on a range of challenging problems, including the energy build-up and onset of reconnection, partially ionized regimes, the influence of magnetic turbulence, and particle acceleration.

preprint2020arXiv

Mode Converting Alfvén Waves from Magnetic Reconnection Enhancing the Energy Source for the Aurora Borealis

Previous studies have concluded that the Hall magnetic field structures generated during magnetic reconnection are carried away by kinetic Alfvén waves (KAW). Here we apply a kinetic simulation with an ion/electron mass ratio closer to its natural value and find that much-reduced damping rates permit the KAW to convert into shear Alfvén waves (SAW). For magnetotail reconnection these SAW provides efficient transport of wave energy, enhancing the energy input for the Aurora Borealis by orders of magnitude above previous estimates.

preprint2020arXiv

Recent Progress on Particle Acceleration and Reconnection Physics during Magnetic Reconnectionin the Magnetically-dominated Relativistic Regime

Magnetic reconnection in strongly magnetized astrophysical plasma environments is believed to be the primary process for fast energy release and particle energization. Currently there is strong interest in relativistic magnetic reconnection, in that it may provide a new explanation for high-energy particle acceleration and radiation in strongly magnetized astrophysical systems. We review recent advances in particle acceleration and reconnection physics in the magnetically-dominated regime. More discussion is focused on the physics of particle acceleration, power-law formation as well as the reconnection rate problem. In addition, we provide an outlook for studying reconnection acceleration mechanisms and kinetic physics in the next step.

preprint2016arXiv

A model of global magnetic reconnection rate in relativistic collisionless plasmas

A model of global magnetic reconnection rate in relativistic collisionless plasmas is developed and validated by the fully kinetic simulation. Through considering the force balance at the upstream and downstream of the diffusion region, we show that the global rate is bounded by a value $\sim 0.3$ even when the local rate goes up to $\sim O(1)$ and the local inflow speed approaches the speed of light in strongly magnetized plasmas. The derived model is general and can be applied to magnetic reconnection under widely different circumstances.

preprint2016arXiv

A Two-Fluid Study of Oblique Tearing Modes in a Force-Free Current Sheet

Kinetic simulations have demonstrated that three-dimensional reconnection in collisionless regimes proceeds through the formation and interaction of magnetic flux ropes, which are generated due to the growth of tearing instabilities at multiple resonance surfaces. Since kinetic simulations are intrinsically expensive, it is desirable to explore the feasibility of reduced two-fluid models to capture this complex evolution, particularly, in the strong guide field regime, where two-fluid models are better justified. With this goal in mind, this paper compares the evolution of the collisionless tearing instability in a force-free current sheet with a two-fluid model and fully kinetic simulations. Our results indicate that the most unstable modes are oblique for guide fields larger than the reconnecting field, in agreement with the kinetic results. The standard two-fluid tearing theory is extended to address the tearing instability at oblique angles. The resulting theory yields a flat oblique spectrum and underestimates the growth of oblique modes in a similar manner to kinetic theory relative to kinetic simulations.

preprint2016arXiv

Efficient Production of High-energy Nonthermal Particles during Magnetic Reconnection in a Magnetically-dominated Ion-Electron Plasma

Magnetic reconnection is a leading mechanism for dissipating magnetic energy and accelerating nonthermal particles in Poynting-flux dominated flows. In this letter, we investigate nonthermal particle acceleration during magnetic reconnection in a magnetically-dominated ion-electron plasma using fully kinetic simulations. For an ion-electron plasma with the total magnetization $σ_0=B^2/(4πn(m_i+m_e)c^2)$, the magnetization for each species is $σ_i \sim σ_0$ and $σ_e \sim (m_i/m_e) σ_0$, respectively. We have studied the magnetically dominated regime by varying $σ_{e} = 10^3 - 10^5$ with initial ion and electron temperatures $T_i = T_e = 5 - 20 m_ec^2$ and mass ratio $m_i/m_e = 1 - 1836$. The results demonstrate that reconnection quickly establishes power-law energy distributions for both electrons and ions within several ($2-3$) light-crossing times. For the cases with periodic boundary conditions, the power-law index is $1<s<2$ for both electrons and ions. The hard spectra limit the power-law energies for electrons and ions to be $γ_{be} \sim σ_e$ and $γ_{bi} \sim σ_i$, respectively. The main acceleration mechanism is a Fermi-like acceleration through the drift motions of charged particles. When comparing the spectra for electrons and ions in momentum space, the spectral indices $s_p$ are identical as predicted in Fermi acceleration. We also find that the bulk flow can carry a significant amount of energy during the simulations. We discuss the implication of this study in the context of Poynting-flux dominated jets and pulsar winds especially the applications for explaining the nonthermal high-energy emissions.

preprint2016arXiv

Particle Acceleration during Magnetic Reconnection in a Low-beta Pair Plasma

Plasma energization through magnetic reconnection in the magnetically-dominated regime featured by low plasma beta ($β= 8 πnkT_0/B^2 \ll 1$) and/or high magnetization ($σ= B^2/(4 πnmc^2) \gg 1$) is important in a series of astrophysical systems such as solar flares, pulsar wind nebula, and relativistic jets from black holes, etc. In this paper, we review the recent progress on kinetic simulations of this process and further discuss plasma dynamics and particle acceleration in a low-$β$ reconnection layer that consists of electron-positron pairs. We also examine the effect of different initial thermal temperatures on the resulting particle energy spectra. While earlier papers have concluded that the spectral index is smaller for higher $σ$, our simulations show that the spectral index approaches $p=1$ for sufficiently low plasma $β$, even if $σ\sim 1$. Since this predicted spectral index in the idealized limit is harder than most observations, it is important to consider effects that can lead to a softer spectrum such as open boundary simulations. We also remark that the effects of 3D reconnection physics and turbulence on reconnection need to be addressed in the future.

preprint2016arXiv

Two-Stage Bulk Electron Heating in the Diffusion Region of Anti-Parallel Symmetric Reconnection

Electron bulk energization in the diffusion region during anti-parallel symmetric reconnection entails two stages. First, the inflowing electrons are adiabatically trapped and energized by an ambipolar parallel electric field. Next, the electrons gain energy from the reconnection electric field as they undergo meandering motion. These collisionless mechanisms have been decribed previously, and they lead to highly-structured electron velocity distributions. Nevertheless, a simplified control-volume analysis gives estimates for how the net effective heating scales with the upstream plasma conditions in agreement with fully kinetic simulations and spacecraft observations.

preprint2015arXiv

Double layer electric fields aiding the production of energetic flat-top distributions and superthermal electrons within the exhausts from magnetic reconnection

Using a kinetic simulation of magnetic reconnection it was recently shown that magnetic-field-aligned electric fields (E||) can be present over large spatial scales in reconnection exhausts. The largest values of E|| are observed within double layers. The existence of double layers in the Earth's magnetosphere is well documented. In our simulation their formation is triggered by large parallel streaming of electrons into the reconnection region. These parallel electron fluxes are required for maintaining quasi-neutrality of the reconnection region and increase with decreasing values of the normalized electron pressure upstream of the reconnection region. A threshold normalized pressure is derived for strong double layers to develop. We also document how the electron confinement, provided in part by the structure in E||, allows sustained energization by perpendicular electric fields. The energization is a consequence of the confined electrons' chaotic orbital motion that includes drifts aligned with the reconnection electric field. The level of energization is proportional to the initial particle energy and therefore is enhanced by the initial energy boost of the acceleration potential, acquired by electrons entering the region. The mechanism is effective in an extended region of the reconnection exhaust allowing for the generation of superthermal electrons in reconnection scenarios, including those with only a single x-line. An expression for the phase-space distribution of the superthermal electrons is derived, providing an accurate match to the kinetic simulation results. The numerical and analytical results agree with detailed spacecraft observations recorded during reconnection events in the Earth's magnetotail.

preprint2015arXiv

Particle Acceleration and Plasma Dynamics during Magnetic Reconnection in the Magnetically-dominated Regime

Magnetic reconnection is thought to be the driver for many explosive phenomena in the universe. The energy release and particle acceleration during reconnection have been proposed as a mechanism for producing high-energy emissions and cosmic rays. We carry out two- and three-dimensional kinetic simulations to investigate relativistic magnetic reconnection and the associated particle acceleration. The simulations focus on electron-positron plasmas starting with a magnetically dominated, force-free current sheet ($σ\equiv B^2/(4πn_e m_e c^2) \gg 1$). For this limit, we demonstrate that relativistic reconnection is highly efficient at accelerating particles through a first-order Fermi process accomplished by the curvature drift of particles along the electric field induced by the relativistic flows. This mechanism gives rise to the formation of hard power-law spectra $f \propto (γ-1)^{-p}$ and approaches $p = 1$ for sufficiently large $σ$ and system size. Eventually most of the available magnetic free energy is converted into nonthermal particle kinetic energy. An analytic model is presented to explain the key results and predict a general condition for the formation of power-law distributions. The development of reconnection in these regimes leads to relativistic inflow and outflow speeds and enhanced reconnection rates relative to non-relativistic regimes. In the three-dimensional simulation, the interplay between secondary kink and tearing instabilities leads to strong magnetic turbulence, but does not significantly change the energy conversion, reconnection rate, or particle acceleration. This study suggests that relativistic reconnection sites are strong sources of nonthermal particles, which may have important implications to a variety of high-energy astrophysical problems.

preprint2015arXiv

The island coalescence problem: scaling of reconnection in extended fluid models including higher-order moments

As modeling of collisionless magnetic reconnection in most space plasmas with realistic parameters is beyond the capability of today's simulations, due to the separation between global and kinetic length scales, it is important to establish scaling relations in model problems so as to extrapolate to realistic scales. Recently, large scale particle-in-cell (PIC) simulations of island coalescence have shown that the time averaged reconnection rate decreases with system size, while fluid systems at such large scales in the Hall regime have not been studied. Here we perform the complementary resistive MHD, Hall MHD and two fluid simulations using a ten-moment model with the same geometry. In contrast to the standard Harris sheet reconnection problem, Hall MHD is insufficient to capture the physics of the reconnection region. Additionally, motivated by the results of a recent set of hybrid simulations which show the importance of ion kinetics in this geometry, we evaluate the efficacy of the ten-moment model in reproducing such results.

preprint2014arXiv

Formation of Hard Power-laws in the Energetic Particle Spectra Resulting from Relativistic Magnetic Reconnection

Using fully kinetic simulations, we demonstrate that magnetic reconnection in relativistic plasmas is highly efficient at accelerating particles through a first-order Fermi process resulting from the curvature drift of particles in the direction of the electric field induced by the relativistic flows. This mechanism gives rise to the formation of hard power-law spectra in parameter regimes where the energy density in the reconnecting field exceeds the rest mass energy density $σ\equiv B^2/(4 πn m_ec^2) > 1$ and when the system size is sufficiently large. In the limit $σ\gg 1$, the spectral index approaches $p=1$ and most of the available energy is converted into non-thermal particles. A simple analytic model is proposed which explains these key features and predicts a general condition under which hard power-law spectra will be generated from magnetic reconnection.

preprint2014arXiv

Scaling of magnetic reconnection in relativistic collisionless plasmas

Using fully kinetic simulations, we study the scaling of the inflow speed of collisionless magnetic reconnection from the non-relativistic to ultra-relativistic limit. In the anti-parallel configuration, the inflow speed increases with the upstream magnetization parameter $σ$ and approaches the light speed when $σ> O(100)$, leading to an enhanced reconnection rate. In all regimes, the divergence of pressure tensor is the dominant term responsible for breaking the frozen-in condition at the x-line. The observed scaling agrees well with a simple model that accounts for the Lorentz contraction of the plasma passing through the diffusion region. The results demonstrate that the aspect ratio of the diffusion region remains $\sim 0.1$ in both the non-relativistic and relativistic limits.

preprint2011arXiv

Phase Diagram for Magnetic Reconnection in Heliophysical, Astrophysical and Laboratory Plasmas

Recent progress in understanding the physics of magnetic reconnection is conveniently summarized in terms of a phase diagram which organizes the essential dynamics for a wide variety of applications in heliophysics, laboratory and astrophysics. The two key dimensionless parameters are the Lundquist number and the macrosopic system size in units of the ion sound gyroradius. In addition to the conventional single X-line collisional and collisionless phases, multiple X-line reconnection phases arise due to the presence of the plasmoid instability either in collisional and collisionless current sheets. In particular, there exists a unique phase termed "multiple X-line hybrid phase" where a hierarchy of collisional islands or plasmoids is terminated by a collisionless current sheet, resulting in a rapid coupling between the macroscopic and kinetic scales and a mixture of collisional and collisionless dynamics. The new phases involving multiple X-lines and collisionless physics may be important for the emerging applications of magnetic reconnection to accelerate charged particles beyond their thermal speeds. A large number of heliophysical and astrophysical plasmas are surveyed and grouped in the phase diagram: Earth's magnetosphere, solar plasmas (chromosphere, corona, wind and tachocline), galactic plasmas (molecular clouds, interstellar media, accretion disks and their coronae, Crab nebula, Sgr A*, gamma ray bursts, magnetars), extragalactic plasmas (Active Galactic Nuclei disks and their coronae, galaxy clusters, radio lobes, and extragalactic jets). Significance of laboratory experiments, including a next generation reconnection experiment, is also discussed.

preprint2011arXiv

Vlasov simulation in multiple spatial dimensions

A long-standing challenge encountered in modeling plasma dynamics is achieving practical Vlasov equation simulation in multiple spatial dimensions over large length and time scales. While direct multi-dimension Vlasov simulation methods using adaptive mesh methods [J. W. Banks et al., Physics of Plasmas 18, no. 5 (2011): 052102; B. I. Cohen et al., November 10, 2010, http://meetings.aps.org/link/BAPS.2010.DPP.NP9.142] have recently shown promising results, in this paper we present an alternative, the Vlasov Multi Dimensional (VMD) model, that is specifically designed to take advantage of solution properties in regimes when plasma waves are confined to a narrow cone, as may be the case for stimulated Raman scatter in large optic f# laser beams. Perpendicular grid spacing large compared to a Debye length is then possible without instability, enabling an order 10 decrease in required computational resources compared to standard particle in cell (PIC) methods in 2D, with another reduction of that order in 3D. Further advantage compared to PIC methods accrues in regimes where particle noise is an issue. VMD and PIC results in a 2D model of localized Langmuir waves are in qualitative agreement.