Source author record

Fan Guo

Fan Guo 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

47works
10topics
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

47 published item(s)

preprint2026arXiv

Heterogeneous Multi-Expert Reinforcement Learning for Long-Horizon Multi-Goal Tasks in Autonomous Forklifts

Autonomous mobile manipulation in unstructured warehouses requires a balance between efficient large-scale navigation and high-precision object interaction. Traditional end-to-end learning approaches often struggle to handle the conflicting demands of these distinct phases. Navigation relies on robust decision-making over large spaces, while manipulation needs high sensitivity to fine local details. Forcing a single network to learn these different objectives simultaneously often causes optimization interference, where improving one task degrades the other. To address these limitations, we propose a Heterogeneous Multi-Expert Reinforcement Learning (HMER) framework tailored for autonomous forklifts. HMER decomposes long-horizon tasks into specialized sub-policies controlled by a Semantic Task Planner. This structure separates macro-level navigation from micro-level manipulation, allowing each expert to focus on its specific action space without interference. The planner coordinates the sequential execution of these experts, bridging the gap between task planning and continuous control. Furthermore, to solve the problem of sparse exploration, we introduce a Hybrid Imitation-Reinforcement Training Strategy. This method uses expert demonstrations to initialize the policy and Reinforcement Learning for fine-tuning. Experiments in Gazebo simulations show that HMER significantly outperforms sequential and end-to-end baselines. Our method achieves a task success rate of 94.2\% (compared to 62.5\% for baselines), reduces operation time by 21.4\%, and maintains placement error within 1.5 cm, validating its efficacy for precise material handling.

preprint2022arXiv

A model of double coronal hard X-ray sources in solar flares

A number of double coronal X-ray sources have been observed during solar flares by RHESSI, where the two sources reside at different sides of the inferred reconnection site. However, where and how are these X-ray-emitting electrons accelerated remains unclear. Here we present the first model of the double coronal hard X-ray (HXR) sources, where electrons are accelerated by a pair of termination shocks driven by bi-directional fast reconnection outflows. We model the acceleration and transport of electrons in the flare region by numerically solving the Parker transport equation using velocity and magnetic fields from the macroscopic magnetohydrodynamic simulation of a flux rope eruption. We show that electrons can be efficiently accelerated by the termination shocks and high-energy electrons mainly concentrate around the two shocks. The synthetic HXR emission images display two distinct sources extending to $>$100 keV below and above the reconnection region, with the upper source much fainter than the lower one. The HXR energy spectra of the two coronal sources show similar spectral slopes, consistent with the observations. Our simulation results suggest that the flare termination shock can be a promising particle acceleration mechanism in explaining the double-source nonthermal emissions in solar flares.

preprint2022arXiv

Advancing Theory and Modeling Efforts in Heliophysics

Heliophysics theory and modeling build understanding from fundamental principles to motivate, interpret, and predict observations. Together with observational analysis, they constitute a comprehensive scientific program in heliophysics. As observations and data analysis become increasingly detailed, it is critical that theory and modeling develop more quantitative predictions and iterate with observations. Advanced theory and modeling can inspire and greatly improve the design of new instruments and increase their chance of success. In addition, in order to build physics-based space weather forecast models, it is important to keep developing and testing new theories, and maintaining constant communications with theory and modeling. Maintaining a sustainable effort in theory and modeling is critically important to heliophysics. We recommend that all funding agencies join forces and consider expanding current and creating new theory and modeling programs--especially, 1. NASA should restore the HTMS program to its original support level to meet the critical needs of heliophysics science; 2. a Strategic Research Model program needs to be created to support model development for next-generation basic research codes; 3. new programs must be created for addressing mission-critical theory and modeling needs; and 4. enhanced programs are urgently required for training the next generation of theorists and modelers.

preprint2022arXiv

Double-power-law feature of energetic particles accelerated at coronal shocks

Recent observations have shown that in many large solar energetic particle (SEP) events the event-integrated differential spectra resemble double power laws. We perform numerical modeling of particle acceleration at coronal shocks propagating through a streamer-like magnetic field by solving the Parker transport equation, including protons and heavier ions. We find that for all ion species the energy spectra integrated over the simulation domain can be described by a double power law, and the break energy depends on the ion charge-to-mass ratio as $E_B \sim (Q/A)^α$, with $α$ varying from 0.16 to 1.2 by considering different turbulence spectral indices. We suggest that the double power law distribution may emerge as a result of the superposition of energetic particles from different source regions where the acceleration rates differ significantly due to particle diffusion. The diffusion and mixing of energetic particles could also provide an explanation for the increase of Fe/O at high energies as observed in some SEP events. Although further mixing processes may occur, our simulations indicate that either power-law break or rollover can occur near the Sun and predict that the spectral forms vary significantly along the shock front, which may be examined by upcoming near-Sun SEP measurements from Parker Solar Probe and Solar Orbiter.

preprint2022arXiv

Modeling Electron Acceleration and Transport in the Early Impulsive Phase of the 2017 September 10 Solar Flare

The X8.2-class limb flare on September 10, 2017 is among the best studied solar flare events owing to its great similarity to the standard flare model and the broad coverage by multiple spacecraft and ground-based observations. These multiwavelength observations indicate that electron acceleration and transport are efficient in the reconnection and flare looptop regions. However, there lacks a comprehensive model for explaining and interpreting the multi-faceted observations. In this work, we model the electron acceleration and transport in the early impulsive phase of this flare. We solve the Parker transport equation that includes the primary acceleration mechanism during magnetic reconnection in the large-scale flare region modeled by MHD simulations. We find that electrons are accelerated up to several MeV and fill a large volume of the reconnection region, similar to the observations shown in microwaves. The electron spatial distribution and spectral shape in the looptop region agree well with those derived from the microwave and hard X-ray emissions before magnetic islands grow large and dominate the acceleration. Future emission modelings using the electron maps will enable direct comparison with microwave and hard X-ray observations. These results shed new light on the electron acceleration and transport in a broad region of solar flares within a data-constrained realistic flare geometry.

preprint2022arXiv

Particle Acceleration in Magnetic Reconnection with Ad hoc Pitch-angle Scattering

Particle acceleration during magnetic reconnection is a long-standing topic in space, solar and astrophysical plasmas. Recent 3D particle-in-cell simulations of magnetic reconnection show that particles can leave flux ropes due to 3D field-line chaos, allowing particles to access additional acceleration sites, gain more energy through Fermi acceleration, and develop a power-law energy distribution. This 3D effect does not exist in traditional 2D simulations, where particles are artificially confined to magnetic islands due to their restricted motions across field lines. Full 3D simulations, however, are prohibitively expensive for most studies. Here, we attempt to reproduce 3D results in 2D simulations by introducing ad hoc pitch-angle scattering to a small fraction of the particles. We show that scattered particles are able to transport out of 2D islands and achieve more efficient Fermi acceleration, leading to a significant increase of energetic particle flux. We also study how the scattering frequency influences the nonthermal particle spectra. This study helps achieve a complete picture of particle acceleration in magnetic reconnection.

preprint2021arXiv

Hybrid Paricle-in-Cell Simulations of Electromagnetic Coupling and Waves From Streaming Burst Debris

Various systems can be modeled as a point-like explosion of ionized debris into a magnetized, collisionless background plasma -- including astrophysical examples, active experiments in space, and laser-driven laboratory experiments. Debris streaming from the explosion parallel to the magnetic field may drive multiple resonant and non-resonant ion-ion beam instabilities, some of which can efficiently couple the debris energy to the background and may even support the formation of shocks. We present a large-scale hybrid (kinetic ions + fluid electrons) particle-in-cell (PIC) simulation, extending hundreds of ion inertial lengths from a 3-D explosion, that resolves these instabilities. We show that the character of these instabilities differs notably from the 1-D equivalent by the presence of unique transverse structure. Additional 2-D simulations explore how the debris beam length, width, density, and speed affect debris-background coupling, with implications for the generation of quasi-parallel shocks.

preprint2021arXiv

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets--II. Connection with $γ$-rays

It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multi-wavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multi-wavelength flares. Nevertheless, previous works have not directly evaluated $γ$-ray signatures from first-principle simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multi-wavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and {\it Fermi-LAT} $γ$-ray bands as well as closely correlated optical and $γ$-ray flares. The optical polarization angle swings are also accompanied by $γ$-ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fast $γ$-ray flares or orphan $γ$-ray flares under the leptonic scenario, complementary to the frequently considered mini-jet scenario. It may also infer neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very high energy.

preprint2021arXiv

Role of Parallel Solenoidal Electric Field on Energy Conversion in 2.5D Decaying Turbulence with a Guide Magnetic Field

We perform 2.5D particle-in-cell simulations of decaying turbulence in the presence of a guide (out-of-plane) background magnetic field. The fluctuating magnetic field initially consists of Fourier modes at low wavenumbers (long wavelengths). With time, the electromagnetic energy is converted to plasma kinetic energy (bulk flow+thermal energy) at the rate per unit volume of ${\pp J}\cdot{\pp E}$ for current density ${\pp J}$ and electric field ${\pp E}$. Such decaying turbulence is well known to evolve toward a state with strongly intermittent plasma current. Here we decompose the electric field into components that are irrotational, ${\pp E}_{\rm ir}$, and solenoidal (divergence-free), ${\pp E}_{\rm so}$. ${\pp E}_{\rm ir}$ is associated with charge separation, and ${\pp J}\cdot{\pp E}_{\rm ir}$ is a rate of energy transfer between ions and electrons with little net change in plasma kinetic energy. Therefore, the net rate of conversion of electromagnetic energy to plasma kinetic energy is strongly dominated by ${\pp J}\cdot{\pp E}_{\rm so}$, and for a strong guide magnetic field, this mainly involves the component ${\pp E}_{\rm so,\parallel}$ parallel to the total magnetic field ${\pp B}$. We examine various indicators of the spatial distribution of the energy transfer rate {\bf J$_\parallel\cdot$E$_{so,\parallel}$}, which relates to magnetic reconnection, the best of which are 1) the ratio of the out-of-plane electric field to the in-plane magnetic field, 2) the out-of-plane component of the non-ideal electric field, and 3) the magnitude of the estimate of current helicity.

preprint2021arXiv

Shock Propagation and Associated Particle Acceleration in the Presence of Ambient Solar-Wind Turbulence

The topic of this review paper is on the influence of solar wind turbulence on shock propagation and its consequence on the acceleration and transport of energetic particles at shocks. As the interplanetary shocks sweep through the turbulent solar wind, the shock surfaces fluctuate and ripple in a range of different scales. We discuss particle acceleration at rippled shocks in the presence of ambient solar-wind turbulence. This strongly affects particle acceleration and transport of energetic particles (both ions and electrons) at shock fronts. In particular, we point out that the effects of upstream turbulence is critical for understanding the variability of energetic particles at shocks. Moreover, the presence of pre-existing upstream turbulence significantly enhances the trapping near the shock of low-energy charged particles, including near the thermal energy of the incident plasma, even when the shock propagates normal to the average magnetic field. Pre-existing turbulence, always present in space plasmas, provides a means for the efficient acceleration of low-energy particles and overcoming the well known injection problem at shocks.

preprint2021arXiv

Slow Shock Formation Upstream of Reconnecting Current Sheets

The formation, development and impact of slow shocks in the upstream region of reconnecting current layers are explored. Slow shocks have been documented in the upstream region of magnetohydrodynamic (MHD) simulations of magnetic reconnection as well as in similar simulations with the {\it kglobal} kinetic macroscale simulation model. They are therefore a candidate mechanism for preheating the plasma that is injected into the current layers that facilitate magnetic energy release in solar flares. Of particular interest is their potential role in producing the hot thermal component of electrons in flares. During multi-island reconnection, the formation and merging of flux ropes in the reconnecting current layer drives plasma flows and pressure disturbances in the upstream region. These pressure disturbances steepen into slow shocks that propagate along the reconnecting component of the magnetic field and satisfy the expected Rankine-Hugoniot jump conditions. Plasma heating arises from both compression across the shock and the parallel electric field that develops to maintain charge neutrality in a kinetic system. Shocks are weaker at lower plasma $β$, where shock steepening is slow. While these upstream slow shocks are intrinsic to the dynamics of multi-island reconnection, their contribution to electron heating remains relatively minor compared with that from Fermi reflection and the parallel electric fields that bound the reconnection outflow.

preprint2021arXiv

Variable Ion Compositions of Solar Energetic Particle Events in the Inner Heliosphere: A Field-line Braiding Model with Compound Injections

We propose a model for interpreting highly variable ion composition ratios in solar energetic particles (SEP) events recently observed by Parker Solar Probe (PSP) at $0.3 - 0.45$ astronomical unit. We use numerical simulations to calculate SEP propagation in a turbulent interplanetary magnetic field with a Kolmogorov power spectrum from large scale down to the gyration scale of energetic particles. We show that when the source regions of different species are offset by a distance comparable to the size of the source regions, the observed energetic particle composition He/H can be strongly variable over more than two orders of magnitude, even if the source ratio is at the nominal value. Assuming a $^3$He/$^4$He source ratio of $10 \%$ in impulsive $^3$He-rich events and the same spatial offset of the source regions, the $^3$He/$^4$He ratio at observation sites also vary considerably. The variability of the ion composition ratios depends on the radial distance, which can be tested by observations made at different radial locations. We discuss the implication of these results on the variability of ion composition of impulsive events and on further PSP and Solar Orbiter observations close to the Sun.

preprint2020arXiv

Dynamical modulation of solar flare electron acceleration due to plasmoid-shock interactions in the looptop region

A fast-mode shock can form in the front of reconnection outflows and has been suggested as a promising site for particle acceleration in solar flares. Recent development of magnetic reconnection has shown that numerous plasmoids can be produced in a large-scale current layer. Here we investigate the dynamical modulation of electron acceleration in the looptop region when plasmoids intermittently arrive at the shock by combining magnetohydrodynamics simulations with a particle kinetic model. As plasmoids interact with the shock, the looptop region exhibits various compressible structures that modulate the production of energetic electrons. The energetic electron population varies rapidly in both time and space. The number of 5$-$10 keV electrons correlates well with the area with compression, while that of $>$50 keV electrons shows good correlation with strong compression area but only moderate correlation with shock parameters. We further examine the impacts of the first plasmoid, which marks the transition from a quasi-steady shock front to a distorted and dynamical shock. The number of energetic electrons is reduced by $\sim 20\%$ at 15$-$25 keV and nearly 40\% for 25$-$50 keV, while the number of 5$-$10 keV electrons increases. In addition, the electron energy spectrum above 10 keV evolves softer with time. We also find double or even multiple distinct sources can develop in the looptop region when the plasmoids move across the shock. Our simulations have strong implications to the interpretation of nonthermal looptop sources, as well as the commonly observed fast temporal variations in flare emissions, including the quasi-periodic pulsations.

preprint2020arXiv

Exploring the acceleration mechanisms for particle injection and power-law formation during trans-relativistic magnetic reconnection

Magnetic reconnection in the relativistic and trans-relativistic regimes is able to accelerate particles to hard power law energy spectra $f \propto γ^{-p}$ (approaching $p=1$). The underlying acceleration mechanism that determines the spectral shape is currently a topic of intense investigation. By means of fully kinetic plasma simulations, we carry out a study of particle acceleration during magnetic reconnection in the trans-relativistic regime of a proton-electron plasma. While earlier work in this parameter regime has focused on the effects of electric field parallel to the local magnetic field on the particle injection (from thermal energy to the lower energy bound of the power-law spectrum), here we examine the roles of both parallel and perpendicular electric fields to gain a more complete understanding on the injection process and further development of a power-law spectrum. We show that the parallel electric field does contribute significantly to particle injection, and is more important in the initial phase of magnetic reconnection. However, as the simulation proceeds, the acceleration by the perpendicular electric field becomes more important for particle injection and completely dominates the acceleration responsible for the high-energy power-law spectrum. This holds robustly, in particular for longer reconnection times and larger systems, i.e. in simulations that are more indicative of the processes in astrophysical sources.

preprint2020arXiv

Fermi-type particle acceleration from magnetic reconnection at the termination shock of a relativistic striped wind

An oblique-rotating pulsar generates a relativistic striped wind in a pulsar wind nebula (PWN). The termination shock of the PWN compresses the Poynting-flux-dominated flow and drives magnetic reconnection. By carrying out particle-in-cell (PIC) simulations of the termination shock of the PWN, we study the shock structure as well as the energy conversion processes and particle acceleration mechanisms. With the recent advances in the numerical methods, we extend the simulations to the ultra-relativistic regime with bulk Lorentz factor up to γ_{0}=10^{6}. Magnetic reconnection at the termination shock is highly efficient at converting magnetic energy to particle kinetic energy and accelerating particles to high energies. We find that the resulting energy spectra crucially depend on λ/d_{e}. When λ/d_{e} is large (λ\gtrsim40d_{e}) , the downstream particle spectra form a power-law distribution in the magnetically dominated relativistic wind regime with upstream magnetization parameter σ_{0}=10. By analyzing particle trajectories and statistical quantities relevant to particle energization, we find that Fermi-type mechanism dominates the particle acceleration and power-law formation. We find that the results for particle acceleration are scalable as γ_{0} and σ_{0} increase to large values. The maximum energy for electrons and positrons can reach hundreds of TeV if the wind has a bulk Lorentz factor γ_{0}\approx10^{6} and magnetization parameter σ_{0}=10, which can explain the recent observations of high-energy gamma-rays from pulsar wind nebulae (PWNe).

preprint2020arXiv

Measurement of magnetic field and relativistic electrons along a solar flare current sheet

In the standard model of solar flares, a large-scale reconnection current sheet is postulated as the central engine for powering the flare energy release and accelerating particles. However, where and how the energy release and particle acceleration occur remain unclear due to the lack of measurements for the magnetic properties of the current sheet. Here we report the measurement of spatially-resolved magnetic field and flare-accelerated relativistic electrons along a current-sheet feature in a solar flare. The measured magnetic field profile shows a local maximum where the reconnecting field lines of opposite polarities closely approach each other, known as the reconnection X point. The measurements also reveal a local minimum near the bottom of the current sheet above the flare loop-top, referred to as a "magnetic bottle". This spatial structure agrees with theoretical predictions and numerical modeling results. A strong reconnection electric field of ~4000 V/m is inferred near the X point. This location, however, shows a local depletion of microwave-emitting relativistic electrons. These electrons concentrate instead at or near the magnetic bottle structure, where more than 99% of them reside at each instant. Our observations suggest that the loop-top magnetic bottle is likely the primary site for accelerating and/or confining the relativistic electrons.

preprint2020arXiv

PIC simulations of microinstabilities and waves at near-Sun solar wind perpendicular shocks: Predictions for Parker Solar Probe and Solar Orbiter

Microinstabilities and waves excited at moderate-Mach-number perpendicular shocks in the near-Sun solar wind are investigated by full particle-in-cell (PIC) simulations. By analyzing the dispersion relation of fluctuating field components directly issued from the shock simulation, we obtain key findings concerning wave excitations at the shock front: (1) at the leading edge of the foot, two types of electrostatic (ES) waves are observed. The relative drift of the reflected ions versus the electrons triggers an electron cyclotron drift instability (ECDI) which excites the first ES wave. Because the bulk velocity of gyro-reflected ions shifts to the direction of the shock front, the resulting ES wave propagates oblique to the shock normal. Immediately, a fraction of incident electrons are accelerated by this ES wave and a ring-like velocity distribution is generated. They can couple with the hot Maxwellian core and excite the second ES wave around the upper hybrid frequency. (2) from the middle of the foot all the way to the ramp, electrons can couple with both incident and reflected ions. ES waves excited by ECDI in different directions propagate across each other. Electromagnetic (EM) waves (X mode) emitted toward upstream are observed in both regions. They are probably induced by a small fraction of relativistic electrons. Results shed new insight on the mechanism for the occurrence of ES wave excitations and possible EM wave emissions at young CME-driven shocks in the near-Sun solar wind.

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.

preprint2020arXiv

Studying particle acceleration from driven magnetic reconnection at the termination shock of a relativistic striped wind using particle-in-cell simulations

A rotating pulsar creates a surrounding pulsar wind nebula (PWN) by steadily releasing an energetic wind into the interior of the expanding shockwave of supernova remnant or interstellar medium. At the termination shock of a PWN, the Poynting-flux-dominated relativistic striped wind is compressed. Magnetic reconnection is driven by the compression and converts magnetic energy into particle kinetic energy and accelerating particles to high energies. We carrying out particle-in-cell (PIC) simulations to study the shock structure as well as the energy conversion and particle acceleration mechanism. By analyzing particle trajectories, we find that many particles are accelerated by Fermi-type mechanism. The maximum energy for electrons and positrons can reach hundreds of TeV.

preprint2020arXiv

Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations

The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov instability (NCI) growth rate for arbitrary time steps allowed by the Courant condition. While NCI from higher order resonances is still present, the numerical tests show that for smaller time steps, the numerical instability grows much slower than using the optimal time step found in previous studies. The WT scheme is efficient for improving the quality and flexibility of relativistic PIC simulations.

preprint2019arXiv

Energy Dissipation and Entropy in Collisionless Plasma

It is well known that collisionless systems are dissipation free from the perspective of particle collision and thus conserve entropy. On the other hand, processes such as magnetic reconnection and turbulence appear to convert large-scale magnetic energy into heat. In this paper, we investigate the energization and heating of collisionless plasma. The dissipation process is discussed in terms of fluid entropy in both isotropic and gyrotropic forms. Evolution equations for the entropy are derived and they reveal mechanisms that lead to changes in fluid entropy. These equations are verified by a collisionless particle-in-cell simulation of multiple reconnecting current sheets. In addition to previous findings regarding the pressure tensor, we emphasize the role of heat flux in the dissipation process.

preprint2019arXiv

Heating of Heavy Ions in Low-beta Compressible Turbulence

Enhancement of minor ions such as $^3$He and heavy ions in flare-associated solar energetic particle (SEP) events remains one of the major puzzles in heliophysics. In this work, we use 3D hybrid simulations (kinetic protons and fluid electrons) to investigate particle energization in a turbulent low-beta environment similar to solar flares. It is shown that in this regime the injected large-amplitude Alfvén waves develop into compressible and anisotropic turbulence, which efficiently heats thermal ions of different species. We find that temperature increase of heavy ions is inversely proportional to the charge to mass ratio, which is consistent with observations of impulsive SEP events. Further analysis reveals that ions are energized by interacting with nearly perpendicular magnetosonic waves near proton inertial scale.

preprint2019arXiv

The Acceleration and Confinement of Energetic Electrons by a Termination Shock in a Magnetic Trap: An Explanation for Nonthermal Loop-top Sources during Solar Flares

Nonthermal loop-top sources in solar flares are the most prominent observational signature that suggests energy release and particle acceleration in the solar corona. Although several scenarios for particle acceleration have been proposed, the origin of the loop-top sources remains unclear. Here we present a model that combines a large-scale magnetohydrodynamic simulation of a two-ribbon flare with a particle acceleration and transport model for investigating electron acceleration by a fast-mode termination shock at the looptop. Our model provides spatially resolved electron distribution that evolves in response to the dynamic flare geometry. We find a concave-downward magnetic structure located below the flare termination shock, induced by the fast reconnection downflows. It acts as a magnetic trap to confine the electrons at the looptop for an extended period of time. The electrons are energized significantly as they cross the shock front, and eventually build up a power-law energy spectrum extending to hundreds of keV. We suggest that this particle acceleration and transport scenario driven by a flare termination shock is a viable interpretation for the observed nonthermal loop-top sources.

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.

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

Dissipation in PIC simulations of moderate to low \b{eta} plasma turbulence

We simulate decaying turbulence in electron-positron pair plasmas using a fully- kinetic particle-in-cell (PIC) code. We run two simulations with moderate-to-low plasma beta. The energy decay rate is found to be similar in both the cases. The perpendicular wave-number spectrum of magnetic energy shows a slope of k^-1.3 in both the cases. The particle energy distribution function shows the formation of a non-thermal feature in the case of lower plasma beta, with a slope close to E^-1. The role of thin turbulent current sheets in this process is investigated. The heating by E_{\parallel}.J_{\parallel} term dominates the E_{\perp}.J_{\perp} term. Regions of strong E_{\parallel}.J_{\parallel} are spatially well-correlated with regions of intense current sheets, which also appear correlated with regions of strong E_{\parallel} in the low beta simulation, suggesting an important role of magnetic reconnection in the dissipation of low beta plasma turbulence.

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

Electron Acceleration at a Coronal Shock Propagating Through a Large-scale Streamer-like Magnetic Field

With a test-particle simulation, we investigate the effect of large-scale coronal magnetic fields on electron acceleration at an outward-propagating coronal shock with a circular front. The coronal field is approximated by an analytical solution with a streamer-like magnetic field featured by partially open magnetic field and a current sheet at the equator atop the closed region. We show that the large-scale shock-field configuration, especially the relative curvature of the shock and the magnetic field line across which the shock is sweeping, plays an important role in the efficiency of electron acceleration. At low shock altitudes, when the shock curvature is larger than that of magnetic field lines, the electrons are mainly accelerated at the shock flanks; at higher altitudes, when the shock curvature is smaller, the electrons are mainly accelerated at the shock nose around the top of closed field lines. The above process reveals the shift of efficient electron acceleration region along the shock front during its propagation. It is also found that in general the electron acceleration at the shock flank is not so efficient as that at the top of closed field since at the top a collapsing magnetic trap can be formed. In addition, we find that the energy spectra of electrons is power-law like, first hardening then softening with the spectral index varying in a range of -3 to -6. Physical interpretations of the results and implications on the study of solar radio bursts are discussed.

preprint2016arXiv

Observation of a Metric Type N Solar Radio Burst

Type III and type-III-like radio bursts are produced by energetic electron beams guided along coronal magnetic fields. As a variant of type III bursts, Type N bursts appear as the letter "N" in the radio dynamic spectrum and reveal a magnetic mirror effect in coronal loops. Here, we report a well-observed N-shaped burst consisting of three successive branches at metric wavelength with both fundamental and harmonic components and a high brightness temperature ($>$10$^9$ K). We verify the burst as a true type N burst generated by the same electron beam from three aspects of the data. First, durations of the three branches at a given frequency increase gradually, may due to the dispersion of the beam along its path. Second, the flare site, as the only possible source of non-thermal electrons, is near the western feet of large-scale closed loops. Third, the first branch and the following two branches are localized at different legs of the loops with opposite sense of polarization. We also find that the sense of polarization of the radio burst is in contradiction to the O-mode and there exists a fairly large time delay ($\sim$3-5 s) between the fundamental and harmonic components. Possible explanations accounting for these observations are presented. Assuming the classical plasma emission mechanism, we can infer coronal parameters such as electron density and magnetic field near the radio source and make diagnostics on the magnetic mirror process.

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.

preprint2015arXiv

Nonthermally Dominated Electron Acceleration during Magnetic Reconnection in a Low-beta Plasma

By means of fully kinetic simulations, we investigate electron acceleration during magnetic reconnection in a nonrelativistic proton--electron plasma with conditions similar to solar corona and flares. We demonstrate that reconnection leads to a nonthermally dominated electron acceleration with a power-law energy distribution in the nonrelativistic low-$β$ regime but not in the high-$β$ regime, where $β$ is the ratio of the plasma thermal pressure and the magnetic pressure. The accelerated electrons contain most of the dissipated magnetic energy in the low-$β$ regime. A guiding-center current description is used to reveal the role of electron drift motions during the bulk nonthermal energization. We find that the main acceleration mechanism is a \textit{Fermi}-type acceleration accomplished by the particle curvature drift motion along the electric field induced by the reconnection outflows. Although the acceleration mechanism is similar for different plasma $β$, low-$β$ reconnection drives fast acceleration on Alfvénic timescales and develops power laws out of thermal distribution. The nonthermally dominated acceleration resulting from magnetic reconnection in low-$β$ plasma may have strong implications for the highly efficient electron acceleration in solar flares and other astrophysical systems.

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

Polarization Swings Reveal Magnetic Energy Dissipation in Blazars

The polarization signatures of the blazar emissions are known to be highly variable. In addition to small fluctuations of the polarization angle around a mean value, sometimes large (> 180^o) polarization angle swings are observed. We suggest that such p henomena can be interpreted as arising from light-travel-time effects within an underlying axisymmetric emission region. We present the first simultaneous fitting of the multi-wavelength spectrum, variability and time-dependent polarization features of a correlated optical and gamma-ray flaring event of the prominent blazar 3C279, which was accompanied by a drastic change of its polarization signatures. This unprecedented combination of spectral, variability, and polarization information in a coherent physical model allows us to place stringent constraints on the particle acceleration and magnetic-field topology in the relativistic jet of a blazar, strongly favoring a scenario in which magnetic energy dissipation is the primary driver of the flare event.

preprint2015arXiv

The Acceleration of Electrons at Collisionless Shocks Moving Through a Turbulent Magnetic Field

We perform a numerical-simulation study of the acceleration of electrons at shocks that propagate through a prespecified, kinematically defined turbulent magnetic field. The turbulence consists of broadband magnetic fluctuations that are embedded in the plasma and cover a range of wavelengths, the smallest of which is larger than the gyroadii of electrons that are initially injected into the system. We find that when the variance of the turbulent component of the upstream magnetic field is sufficiently large -- $σ^2 \sim$ 10 $B_0^2$, where $B_0$ is the strength of the background magnetic field -- electrons can be efficiently accelerated at a collisionless shock regardless of the orientation of the mean upstream magnetic field relative to the shock-normal direction. Since the local angle between the incident magnetic-field vector and the shock-normal vector can be quite large, electrons can be accelerated through shock-drift acceleration at the shock front. In the upstream region, electrons are mirrored back to the shock front leading to multiple shock encounters. Eventually the accelerated electrons are energetic enough that their gyroradii are of the same order as the wavelength of waves that are included in our description of the turbulent magnetic field. Our results are consistent with recent \textit{in situ} observations at Saturn's bow shock. The study may help understand the acceleration of electrons at shocks in space and astrophysical systems.

preprint2014arXiv

A solar type II radio burst from CME-coronal ray interaction: simultaneous radio and EUV imaging

Simultaneous radio and extreme ultraviolet (EUV)/white-light imaging data are examined for a solar type II radio burst occurring on 2010 March 18 to deduce its source location. Using a bow-shock model, we reconstruct the 3-dimensional EUV wave front (presumably the type-II emitting shock) based on the imaging data of the two STEREO spacecraft. It is then combined with the Nançay radio imaging data to infer the 3-dimensional position of the type II source. It is found that the type II source coincides with the interface between the CME EUV wave front and a nearby coronal ray structure, providing evidence that the type II emission is physically related to the CME-ray interaction. This result, consistent with those of previous studies, is based on simultaneous radio and EUV imaging data for the first time.

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

Imaging and spectroscopic observations of magnetic reconnection and chromospheric evaporation in a solar flare

Magnetic reconnection is believed to be the dominant energy release mechanism in solar flares. The standard flare model predicts both downward and upward outflow plasmas with speeds close to the coronal Alfvén speed. Yet, spectroscopic observations of such outflows, especially the downflows, are extremely rare. With observations of the newly launched Interface Region Imaging Spectrograph (IRIS), we report the detection of greatly redshifted ($\sim$125 km s$^{-1}$ along line of sight) Fe {\sc{xxi}} 1354.08Å emission line with a $\sim$100 km s$^{-1}$ nonthermal width at the reconnection site of a flare. The redshifted Fe {\sc{xxi}} feature coincides spatially with the loop-top X-Ray source observed by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). We interpret this large redshift as the signature of downward-moving reconnection outflow/hot retracting loops. Imaging observations from both IRIS and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) also reveal the eruption and reconnection processes. Fast downward-propagating blobs along these loops are also found from cool emission lines (e.g., Si {\sc{iv}}, O {\sc{iv}}, C {\sc{ii}}, Mg {\sc{ii}}) and images of AIA and IRIS. Furthermore, the entire Fe {\sc{xxi}} line is blueshifted by $\sim$260 km s$^{-1}$ at the loop footpoints, where the cool lines mentioned above all exhibit obvious redshift, a result that is consistent with the scenario of chromospheric evaporation induced by downward-propagating nonthermal electrons from the reconnection site.

preprint2014arXiv

Possible role of coronal streamer as magnetically-closed structure in shock-induced energetic electrons and metric type II radio bursts

Two solar type II radio bursts, separated by ~24 hours in time, are examined together. Both events are associated with coronal mass ejections (CMEs) erupting from the same active region (NOAA 11176) beneath a well-observed helmet streamer. We find that the type II emissions in both events ended once the CME/shock fronts passed the white-light streamer tip, which is presumably the magnetic cusp of the streamer. This leads us to conjecture that the closed magnetic arcades of the streamer may play a role in electron acceleration and type II excitation at coronal shocks. To examine such a conjecture, we conduct a test-particle simulation for electron dynamics within a large-scale partially-closed streamer magnetic configuration swept by a coronal shock. We find that the closed field lines play the role of an electron trap, via which the electrons are sent back to the shock front for multiple times, and therefore accelerated to high energies by the shock. Electrons with an initial energy of 300eV can be accelerated to tens of keV concentrating at the loop apex close to the shock front with a counter-streaming distribution at most locations. These electrons are energetic enough to excite Langmuir waves and radio bursts. Considering the fact that most solar eruptions originate from closed field regions, we suggest that the scenario may be important to the generation of more metric type IIs. This study also provides an explanation to the general ending frequencies of metric type IIs at or above 20-30 MHz and the disconnection issue between metric and interplanetary type IIs.

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.

preprint2013arXiv

Small-scale Gradients of Charged Particles in the Heliospheric Magnetic Field

Using numerical simulations of charged-particles propagating in the heliospheric magnetic field, we study small-scale gradients, or "dropouts", in the intensity of solar energetic particles seen at 1 AU. We use two turbulence models, the foot-point random motion model (Jokipii & Parker 1969; Giacalone et al. 2006) and two-component model (Matthaeus et al. 1990), to generate fluctuating magnetic fields similar to spacecraft observations at 1 AU. The turbulence models include a Kolmogorov-like magnetic field power spectrum containing a broad range of spatial scales from those that lead to large-scale field-line random walk to small scales leading to resonant pitch-angle scattering of energetic particles. We release energetic protons (20 keV - 10 MeV) from a spatially compact and instantaneous source. The trajectories of energetic charged particles in turbulent magnetic fields are numerically integrated. Spacecraft observations are mimicked by collecting particles in small windows when they pass the windows at a distance of 1 AU. We show that small-scale gradients in the intensity of energetic particles and velocity dispersions observed by spacecraft can be reproduced using the foot-point random motion model. However, no dropouts is seen in simulations using the two-component magnetic turbulence model. We also show that particle scattering in the solar wind magnetic field needs to be infrequent for intensity dropouts to form.

preprint2013arXiv

The Acceleration of Thermal Protons at Parallel Collisionless Shocks: Three-dimensional Hybrid Simulations

We present three-dimensional hybrid simulations of collisionless shocks that propagate parallel to the background magnetic field to study the acceleration of protons that forms a high-energy tail on the distribution. We focus on the initial acceleration of thermal protons and compare it with results from one-dimensional simulations. We find that for both one- and three-dimensional simulations, particles that end up in the high-energy tail of the distribution later in the simulation gained their initial energy right at the shock. This confirms previous results but is the first to demonstrate this using fully three-dimensional fields. The result is not consistent with the "thermal leakage" model. We also show that the gyrocenters of protons in the three-dimensional simulation can drift away from the magnetic field lines on which they started due to the removal of ignorable coordinates that exist in one- and two-dimensional simulations. Our study clarifies the injection problem for diffusive shock acceleration.

preprint2012arXiv

A broken solar type II radio burst induced by a coronal shock propagating across the streamer boundary

We discuss an intriguing type II radio burst that occurred on 2011 March 27. The dynamic spectrum was featured by a sudden break at about 43 MHz on the well-observed harmonic branch. Before the break, the spectrum drifted gradually with a mean rate of about -0.05 MHz/s. Following the break, the spectrum jumped to lower frequencies. The post-break emission lasted for about three minutes. It consisted of an overall slow drift which appeared to have a few fast drift sub-bands. Simultaneous observations from the Solar TErrestrial RElations Observatory (STEREO) and the Solar Dynamics Observatory (SDO) were also available and are examined for this event. We suggest that the slow-drift period before the break was generated inside a streamer by a coronal eruption driven shock, and the spectral break as well as the relatively wide spectrum after the break is a consequence of the shock crossing the streamer boundary where density drops abruptly. It is suggested that this type of radio bursts can be taken as a unique diagnostic tool for inferring the coronal density structure, as well as the radio emitting source region.

preprint2012arXiv

Effects of Turbulent Magnetic Fields on the Transport and Acceleration of Energetic Charged Particles: Numerical Simulations with Application to Heliospheric Physics

After introduction we focus on: the transport of charged particles, the acceleration of ions at shocks, and the acceleration of electrons at shocks. Chapter 2 studies the propagation of solar energetic particles(SEPs) in turbulent magnetic fields. Particle trajectories in turbulent magnetic fields are numerically integrated. The turbulence includes a Kolmogorov-like power spectrum containing a broad range of scales. Small-scale variations in particle intensities(dropouts) and velocity dispersions can be reproduced. The result gives a constraint on the error of onset analysis for inferring SEP informations. We find that dropouts are rarely produced using the two-component model(Matthaeus et al., 1990). The result questions the turbulence model. Chapter 3 studies the acceleration of ions. We use 3-D hybrid simulations to study the acceleration of low-energy particles at parallel shocks. We find that particles gain energy by reflection at the shock. The protons can move off field lines in 3-D electric and magnetic fields. We also use a stochastic integration method to study diffusive shock acceleration including large-scale magnetic variations. The results can explain the observations of anomalous cosmic rays by Voyager 1. Chapter 4 studies electron acceleration at a shock in a turbulent magnetic field by combining hybrid simulations and test-particle simulations. The acceleration is enhanced by including large-scale turbulence. Since electrons mainly follow field lines, the field-line braiding allows electrons interacting with shock many times. Ripples also contribute to acceleration by mirroring electrons. The process favors perpendicular shocks. We discuss the implication to SEPs by comparing the acceleration of electrons with that of protons. The intensity correlation of electrons and ions in SEPs implies perpendicular shocks play important roles in accelerating particles.

preprint2012arXiv

Particle Acceleration at a Flare Termination Shock: Effect of Large-scale Magnetic Turbulence

We investigate the acceleration of charged particles (both electrons and protons) at collisionless shocks predicted to exist in the vicinity of solar flares. The existence of standing termination shocks has been examined by flare models and numerical simulations e.g., Shibata,Forbes. We study electron energization by numerically integrating the equations of motion of a large number of test-particle electrons in the time-dependent two-dimensional electric and magnetic fields generated from hybrid simulations (kinetic ions and fluid electron) using parameters typical of the solar flare plasma environment. The shock is produced by injecting plasma flow toward a rigid piston. Large-scale magnetic fluctuations -- known to exist in plasmas and known to have important effects on the nonthermal electron acceleration at shocks -- are also included in our simulations. For the parameters characteristic of the flaring region, our calculations suggest that the termination shock formed in the reconnection outflow region (above post-flare loops) could accelerate electrons to a kinetic energy of a few MeV within 100 ion cyclotron periods, which is of the order of a millisecond. Given a sufficient turbulence amplitude level ($δB^2/B_0^2 \sim 0.3$), about 10% of thermal test-particle electrons are accelerated to more than 15 keV. We find that protons are also accelerated, but not to as high energy in the available time and the energy spectra are considerably steeper than that of the electrons for the parameters used in our simulations. Our results are qualitatively consistent with the observed hard X-ray emissions in solar flares.

preprint2012arXiv

The Acceleration of Electrons at Perpendicular Shocks and its Implication for Solar Energetic Particle events

We present a study of the acceleration of electrons at a perpendicular shock that propagates through a turbulent magnetic field. The energization process of electrons is investigated by utilizing a combination of hybrid (kinetic ions and fluid electron) simulations and test-particle electron simulations. In this method, the motions of the test-particle electrons are numerically integrated in the time-dependent electric and magnetic fields generated by two-dimensional hybrid simulations. We show that large-scale magnetic fluctuations effect electrons in a number of ways and lead to efficient and rapid energization at the shock front. Since the electrons mainly follow along magnetic lines of force, the large-scale braiding of field lines in space allows the fast-moving electrons to interact with the shock front and get accelerated multiple times. Ripples in the shock front occurring at various scales will also contribute to the acceleration by mirroring the electrons. Our calculation shows that this process favors electron acceleration at perpendicular shocks. The acceleration efficiency is critically dependent on the turbulence amplitude and coherence length. We also discuss the implication of this study for solar energetic particles (SEPs) by comparing the acceleration of electrons with that of protons. Their correlation indicates that perpendicular shocks play an important role in SEP events.

preprint2011arXiv

On the Amplification of Magnetic Field by a Supernova Blast Shock Wave in a Turbulent Medium

We have performed extensive two-dimensional magnetohydrodynamic simulations to study the amplification of magnetic fields when a supernova blast wave propagates into a turbulent interstellar plasma. The blast wave is driven by injecting high pressure in the simulation domain. The interstellar magnetic field can be amplified by two different processes, occurring in different regions. One is facilitated by the fluid vorticity generated by the ``rippled" shock front interacting with the background turbulence. The resulting turbulent flow keeps amplifying the magnetic field, consistent with earlier work \citep{Giacalone2007}. The other process is facilitated by the growth of the Rayleigh-Taylor instability at the contact discontinuity between the ejecta and the shocked medium. This can efficiently amplify the magnetic field and tends to produce the highest magnetic field. We investigate the dependence of the amplification on numerical parameters such as grid-cell size and on various physical parameters. We show the magnetic field has a characteristic radial profile that the downstream magnetic field gets progressively stronger away from the shock. This is because the downstream magnetic field needs a finite time to reach the efficient amplification, and will get further amplified in the Rayleigh-Taylor region. In our simulation we do not observe a systematic strong magnetic field within a small distance to the shock. This indicates that if the magnetic-field amplification in supernova remnants indeed occurs near the shock front, other processes such as three-dimensional instabilities, plasma kinetics and/or cosmic ray effect may need to be considered to explain the strong magnetic field in supernova remnants.

preprint2010arXiv

The Effect of Large Scale Magnetic Turbulence on the Acceleration of Electrons by Perpendicular Collisionless Shocks

We study the physics of electron acceleration at collisionless shocks that move through a plasma containing large-scale magnetic fluctuations. We numerically integrate the trajectories of a large number of electrons, which are treated as test particles moving in the time dependent electric and magnetic fields determined from 2-D hybrid simulations (kinetic ions, fluid electron). The large-scale magnetic fluctuations effect the electrons in a number of ways and lead to efficient and rapid energization at the shock front. Since the electrons mainly follow along magnetic lines of force, the large-scale braiding of field lines in space allows the fast-moving electrons to cross the shock front several times, leading to efficient acceleration. Ripples in the shock front occuring at various scales will also contribute to the acceleration by mirroring the electrons. Our calculation shows that this process favors electron acceleration at perpendicular shocks. The current study is also helpful in understanding the injection problem for electron acceleration by collisionless shocks. It is also shown that the spatial distribution of energetic electrons is similar to in-situ observations (e.g., Bale et al. 1999; Simnett et al. 2005). The process may be important to our understanding of energetic electrons in planetary bow shocks and interplanetary shocks, and explaining herringbone structures seen in some type II solar radio bursts.