Source author record

Yi-Hsin Liu

Yi-Hsin Liu 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

20works
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

20 published item(s)

preprint2022arXiv

First-Principles Theory of the Rate of Magnetic Reconnection in Magnetospheric and Solar Plasmas

The rate of magnetic reconnection is of the utmost importance in a variety of processes because it controls, for example, the rate energy is released in solar flares, the speed of the Dungey convection cycle in Earth's magnetosphere, and the energy release rate in harmful geomagnetic substorms. It is known from numerical simulations and satellite observations that the rate is approximately 0.1 in normalized units, but despite years of effort, a full theoretical prediction has not been obtained. Here, we present a first-principles theory for the reconnection rate in non-relativistic electron-ion collisionless plasmas, and show that the same prediction explains why Sweet-Parker reconnection is considerably slower. The key consideration of this analysis is the pressure at the reconnection site (i.e., the x-line). We show that the Hall electromagnetic fields in antiparallel reconnection cause an energy void, equivalently a pressure depletion, at the x-line, so the reconnection exhaust opens out, enabling the fast rate of 0.1. If the energy can reach the x-line to replenish the pressure, the exhaust does not open out. In addition to heliospheric applications, these results are expected to impact reconnection studies in planetary magnetospheres, magnetically confined fusion devices, and astrophysical plasmas.

preprint2022arXiv

On the origin of "patchy" energy conversion in electron diffusion regions

During magnetic reconnection, field lines interconnect in electron diffusion regions (EDRs). In some EDRs the reconnection and energy conversion rates are controlled by a steady out-of-plane electric field. In other EDRs the energy conversion rate $\vec{J}\cdot\vec{E}'$ is "patchy", with electron-scale large-amplitude positive and negative peaks. We investigate 22 EDRs observed by NASA's Magnetospheric Multiscale (MMS) mission in a wide range of conditions to determine the cause of patchy $\vec{J}\cdot\vec{E}'$. The patchiness of the energy conversion is quantified and correlated with seven parameters describing various aspects of the asymptotic inflow regions that affect the structure, stability, and efficiency of reconnection. We find that (1) neither the guide field strength nor the asymmetries in the inflow ion pressure, electron pressure, reconnecting magnetic field strength, and number density are well correlated with the patchiness of the EDR energy conversion, (2) the out-of-plane axes of the 22 EDRs are typically fairly well aligned with the "preferred" axes, which bisect the time-averaged inflow magnetic fields and maximize the reconnection rate, and (3) the time-variability in the upstream magnetic field direction is best correlated with the patchiness of the EDR $\vec{J}\cdot\vec{E}'$. A 3-d fully-kinetic simulation of reconnection with a non-uniform inflow magnetic field is analyzed; the variation in the magnetic field generates secondary X-lines, which develop to maximize the reconnection rate for the time-varying inflow magnetic field. The results suggest that magnetopause reconnection, for which the inflow magnetic field direction is often highly variable, may commonly be patchy in space, at least at the electron scale.

preprint2021arXiv

Identification of active magnetic reconnection using magnetic flux transport in plasma turbulence

Magnetic reconnection has been suggested to play an important role in the dynamics and energetics of plasma turbulence by spacecraft observations, simulations and theory over the past two decades, and recently, by magnetosheath observations of MMS. A new method based on magnetic flux transport (MFT) has been developed to identify reconnection activity in turbulent plasmas. This method is applied to a gyrokinetic simulation of two-dimensional (2D) plasma turbulence. Results on the identification of three active reconnection X-points are reported. The first two X-points have developed bi-directional electron outflow jets. Beyond the category of electron-only reconnection, the third X-point does not have bi-directional electron outflow jets because the flow is modified by turbulence. In all cases, this method successfully identifies active reconnection through clear inward and outward flux transport around the X-points. This transport pattern defines reconnection and produces a new quadrupolar structure in the divergence of MFT. This method is expected to be applicable to spacecraft missions such as MMS, Parker Solar Probe, and Solar Orbiter.

preprint2020arXiv

A new Look at the Electron Diffusion Region in Asymmetric Magnetic Reconnection

A new look at the structure of the electron diffusion region in collisionless magnetic reconnection is presented. The research is based on a particle-in-cell simulation of asymmetric magnetic reconnection, which include a temperature gradient across the current layer in addition to density and magnetic field gradient. We find that none of X-point, flow stagnation point, and local current density peak coincide. Current and energy balance analyses around the flow stagnation point and current density peak show consistently that current dissipation is associated with the divergence of nongyrotropic electron pressure. Furthermore, the same pressure terms, when combined with shear-type gradients of the electron flow velocity, also serve to maintain local thermal energy against convective losses. These effects are similar to those found also in symmetric magnetic reconnection. In addition, we find here significant effects related to the convection of current, which we can relate to a generalized diamagnetic drift by the nongyrotropic pressure divergence. Therefore, only part of the pressure force serves to dissipate the current density. However, the prior conclusion that the role of the reconnection electric field is to maintain the current density, which was obtained for a symmetric system, applies here as well. Finally, we discuss related features of electron distribution function in the EDR.

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.

preprint2019arXiv

The Critical Role of Collisionless Plasma Energization on the Structure of Relativistic Magnetic Reconnection

During magnetically dominated relativistic reconnection, inflowing plasma depletes the initial relativistic pressure at the x-line and collisionless plasma heating inside the diffusion region is insufficient to overcome this loss. The resulting pressure drop causes a collapse at the x-line, essentially a localization mechanism of the diffusion region necessary for fast reconnection. The extension of this low-pressure region further explains the bursty nature of anti-parallel reconnection because a once opened outflow exhaust can also collapse, which repeatedly triggers secondary tearing islands. However, a stable single x-line reconnection can be achieved when an external guide field exists, since the reconnecting magnetic field component rotates out of the reconnection plane at outflows, providing additional magnetic pressure to sustain the opened exhausts.

preprint2019arXiv

Three-Dimensional X-line Spreading in Asymmetric Magnetic Reconnection

The spreading of the X-line out of the reconnection plane under a strong guide field is investigated using large-scale three-dimensional (3D) particle-in-cell (PIC) simulations in asymmetric magnetic reconnection. A simulation with a thick, ion-scale equilibrium current sheet (CS) reveals that the X-line spreads at the ambient ion/electron drift speeds, significantly slower than the Alfvén speed based on the guide field $V_{Ag}$. Additional simulations with a thinner, sub-ion-scale CS show that the X-line spreads at $V_{Ag}$ (Alfvénic spreading), much higher than the ambient species drifts. An Alfvénic signal consistent with kinetic Alfvén waves develops and propagates, leading to CS thinning and extending, which then results in reconnection onset. The continuous onset of reconnection in the signal propagation direction manifests as Alfvénic X-line spreading. The strong dependence on the CS thickness of the spreading speeds, and the X-line orientation are consistent with the collisionless tearing instability. Our simulations indicate that when the collisionless tearing growth is sufficiently strong in a thinner CS such that $γ/Ω_{ci}\gtrsim\mathcal{O}(1)$, Alfvénic X-line spreading can take place. Our results compare favorably with a number of numerical simulations and recent magnetopause observations. A key implications is that the magnetopause CS is typically too thick for Alfvénic X-line spreading to effectively take place.

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

Suppression of collisionless magnetic reconnection in asymmetric current sheets

Using fully kinetic simulations, we study the suppression of asymmetric reconnection in the limit where the diamagnetic drift speed >> Alfven speed and the magnetic shear angle is moderate. We demonstrate that the slippage between electrons and the magnetic flux facilitates reconnection, and can even result in fast reconnection that lacks one of the outflow jets. Through comparing a case where the diamagnetic drift is supported by the temperature gradient with a companion case that has a density gradient instead, we identify a robust suppression mechanism. The drift of the x-line is slowed down locally by the asymmetric nature of the current sheet and the resulting tearing modes, then the x-line is run over and swallowed by the faster-moving following flux.

preprint2015arXiv

Orientation of x-lines in asymmetric magnetic reconnection - mass ratio dependency

Using fully kinetic simulations, we study the x-line orientation of magnetic reconnection in an asymmetric configuration. A spatially localized perturbation is employed to induce a single x-line, that has sufficient freedom to choose its orientation in three-dimensional systems. The effect of ion to electron mass ratio is investigated, and the x-line appears to bisect the magnetic shear angle across the current sheet in the large mass ratio limit. The orientation can generally be deduced by scanning through corresponding 2D simulations to find the reconnection plane that maximizes the peak reconnection electric field. The deviation from the bisection angle in the lower mass ratio limit can be explained by the physics of tearing instability.

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.

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

The effects of strong temperature anisotropy on the kinetic structure of collisionless slow shocks and reconnection exhausts. Part I: PIC simulations

A 2-D Riemann problem is designed to study the development and dynamics of the slow shocks that are thought to form at the boundaries of reconnection exhausts. Simulations are carried out for varying ratios of normal magnetic field to the transverse upstream magnetic field (i.e., propagation angle with respect to the upstream magnetic field). When the angle is sufficiently oblique, the simulations reveal a large firehose-sense (P_parallel>P_perpendicular) temperature anisotropy in the downstream region, accompanied by a transition from a coplanar slow shock to a non-coplanar rotational mode. In the downstream region the firehose stability parameter epsilon=1-mu_0(P_parallel-P_perpendicular)/ B^2 tends to lock in to 0.25. This balance arises from the competition between counterstreaming ions, which drives epsilon down, and the scattering due to ion inertial scale waves, which are driven unstable by the downstream rotational wave. At very oblique propagating angles, 2-D turbulence also develops in the downstream region.

preprint2011arXiv

The effects of strong temperature anisotropy on the kinetic structure of collisionless slow shocks and reconnection exhausts. Part II: Theory

Simulations of collisionless oblique propagating slow shocks have revealed the existence of a transition associated with a critical temperature anisotropy epsilon=1-mu_0(P_parallel-P_perpendicular)/ B^2 = 0.25 (Liu, Drake and Swisdak (2011)). An explanation for this phenomenon is proposed here based on anisotropic fluid theory, in particular the Anisotropic Derivative Nonlinear-Schrodinger-Burgers equation, with an intuitive model of the energy closure for the downstream counter-streaming ions. The anisotropy value of 0.25 is significant because it is closely related to the degeneracy point of the slow and intermediate modes, and corresponds to the lower bound of the coplanar to non-coplanar transition that occurs inside a compound slow shock (SS)/rotational discontinuity (RD) wave. This work implies that it is a pair of compound SS/RD waves that bound the outflows in magnetic reconnection, instead of a pair of switch-off slow shocks as in Petschek's model. This fact might explain the rareness of in-situ observations of Petschek-reconnection-associated switch-off slow shocks.

preprint2011arXiv

The structure of the magnetic reconnection exhaust boundary

The structure of shocks that form at the exhaust boundaries during collisionless reconnection of anti-parallel fields is studied using particle-in-cell (PIC) simulations and modeling based on the anisotropic magnetohydrodynamic equations. Large-scale PIC simulations of reconnection and companion Riemann simulations of shock development demonstrate that the pressure anisotropy produced by counterstreaming ions within the exhaust prevents the development of classical Petschek switch-off-slow shocks (SSS). The shock structure that does develop is controlled by the firehose stability parameter epsilon=1-mu_0(P_parallel-P_perpendicular)/ B^2 through its influence on the speed order of the intermediate and slow waves. Here P_parallel and P_perpendicular are the pressure parallel and perpendicular to the local magnetic field. The exhaust boundary is made up of a series of two shocks and a rotational wave. The first shock takes epsilon from unity upstream to a plateau of 0.25 downstream. The condition epsilon =0.25 is special because at this value the speeds of nonlinear slow and intermediate waves are degenerate. The second slow shock leaves epsilon=0.25 unchanged but further reduces the amplitude of the reconnecting magnetic field. Finally, in the core of the exhaust epsilon drops further and the transition is completed by a rotation of the reconnecting field into the out-of-plane direction. The acceleration of the exhaust takes place across the two slow shocks but not during the final rotation. The result is that the outflow speed falls below that expected from the Walen condition based on the asymptotic magnetic field. A simple analytic expression is given for the critical value of epsilon within the exhaust below which SSSs no longer bound the reconnection outflow.

preprint2009arXiv

The Weibel Instability inside the Electron-Positron Harris Sheet

Recent full-particle simulations of electron-positron reconnection have revealed that the Weibel instability plays an active role in controlling the dynamics of the current layer and maintaining fast reconnection. A four-beam model is developed to explore the development of the instability within a narrow current layer characteristic of reconnection. The problem is reduced to two coupled second-order differential equations, whose growing eigenmodes are obtained via both asymptotic approximations and finite difference methods. Full particle simulations confirm the linear theory and help probe the nonlinear development of the instability. The current layer broadening in the reconnection outflow jet is linked to the scattering of high-velocity streaming particles in the Weibel-generated, out-of-plane magnetic field.