Particle Acceleration in the Presence of Weak Turbulence at an X-Type Neutral Point

We simulate the likely noisy situation near a reconnection region by superposing many 2D linear reconnection eigenmodes. The superposition of modes on the steady state X-type magnetic field creates multiple X- and O-type neutral points close to the original neutral point and so increases the size of the non-adiabatic region. We study test particle trajectories of initially thermal protons in these fields. Protons become trapped in this region and are accelerated by the turbulent electric field to energies up to 1 MeV in time scales relevant to solar flares. Higher energies are achieved due to the interaction of particles with increasingly turbulent electric and magnetic fields.     

Stochastic Acceleration of Energetic Particles in the Magnetosphere of Saturn

Voyager's plasma probe observations suggest that there are at least three fundamentally different plasma regimes in Saturn: the hot outer magnetosphere, the extended plasma sheet, and the inner plasma torus. At the outer regions of the inner torus some ions have been accelerated to reach energies of the order of 43 keV. We develop a model that calculates the acceleration of charged particles in the Saturn's magnetosphere. We propose that the stochastic electric field associated to the observed magnetic field fluctuations is responsible of such acceleration. A random electric field is derived from the fluctuating magnetic field – via a Monte Carlo simulation – which then is applied to the momentum equation of charged particles seeded in the magnetosphere. Taking different initial conditions, like the source of charged particles and the distribution function of their velocities, we find that particles injected with very low energies ranging from 0.129 eV to 5.659 keV can be strongly accelerated to reach much higher energies ranging from 22.220 eV to 9.711 keV as a result of 125,000 hitting events (the latter are used in the numerical code to produce the particle acceleration over a predetermined distance).     

Criticism of reconnection models of the magnetosphere

Reconnection involves singular lines called X-lines on the day and night sides of the magnetosphere, and the reconnection rate is proportional to the component of the electric field along the X-line. Although there is some indirect support for this model, nevertheless direct support is totally lacking. However, there are two distinct pieces of clearly contradictory observational evidence on the dayside. First is the failure to account for the implied energy dissipation by the magnetopause current, over 10¹¹ W, which should be easily observable as heating or enhanced flow of the plasma near the magnetopause. In marked contrast to this prediction, HEOS-2 satellite data reveal a plasma with decreased energy density and reduced flow. Second, the boundary of closed magnetic field lines is in the wrong location. In the reconnection process the plasma outflow would cut across open field lines toward higher latitudes; there should be a band of open field lines equatorward of the cleft. Observations of trapped energetic particles indicate closed field lines within the entry layer and cleft. Either one of these pieces of evidence is sufficient by itself to require drastic revision, even rejection, of the reconnection model. There is also contradictory evidence on the night side. The last closed field line capable of trapping energetic particles is poleward of auroral arcs. The implication is that the X-line is at the distant magnetopause, and not in the plasma sheet. Consequently, even if the reconnection process were operative at the nightside X-line, it would be isolated from steady state plasma sheet and auroral processes. On the other hand, substorm phenomena, in which stored magnetic energy is converted into particle kinetic energy, necessarily involve an induced electric field; that is excluded in theories of the reconnection process in which it is assumed that curl E = 0. Nevertheless, the observed easy access of energetic solar flare particles to the polar caps, and especially the preservation of interplanetary anisotropies as differences between the two polar caps, argues strongly for an open magnetosphere, with interconnection between geomagnetic and inter-planetary magnetic field lines. It is suggested that the resolution of this apparent paradox involves electric fields parallel to the magnetic field lines somewhere on the dawn and dusk sides of the magnetosphere, with an equipotential dayside magnetopause.     

Modeling of Proton Acceleration in a Magnetic Island Inside the Ripple of the Heliospheric Current Sheet

Crossings of the heliospheric current sheet (HCS) at the Earth’s orbit are often associated with observations of anisotropic beams of energetic protons accelerated to energies from hundreds of keV to several MeV and above. A connection between this phenomenon and the occurrence of small-scale magnetic islands (SMIs) near reconnecting current sheets has recently been found. This study shows how pre-accelerated protons can be energized additionally due to oscillations of multiple SMIs inside the ripple of the reconnecting HCS. A model of the electromagnetic field of an oscillating 3D SMI with a characteristic size of ~0.001 AU is developed. A SMI is supposed to be bombarded by protons accelerated by magnetic reconnection at the HCS to energies from ~1keV to tens of keV. Numerical simulations have demonstrated that the resulting longitudinal inductive electric fields can additionally reaccelerate protons injected into a SMI. It is shown that there is a local “acceleration” region within the island in which particles gain energy most effectively. As a result, their average escape energies range from hundreds of keV to 2 MeV and above. There is almost no particle acceleration outside the region. It is shown that energies gained by protons significantly depend on the initial phase and the place of their entry into a SMI but weakly depend on the initial energy. Therefore, low-energy particles can be accelerated more efficiently than high-energy particles, and all particles can reach the total energy limit upon their escape from a SMI. It is also found that the escape velocity possesses a strong directional anisotropy. The results are consistent with observations in the solar wind plasma.

     

Field-Guided Proton Acceleration at Reconnecting x-Points in Flares

An explicitly energy-conserving full orbit code CUEBIT, developed originally to describe energetic particle effects in laboratory fusion experiments, has been applied to the problem of proton acceleration in solar flares. The model fields are obtained from solutions of the linearised MHD equations for reconnecting modes at an X-type neutral point, with the additional ingredient of a longitudinal magnetic field component. To accelerate protons to the highest observed energies on flare timescales, it is necessary to invoke anomalous resistivity in the MHD solution. It is shown that the addition of a longitudinal field component greatly increases the efficiency of ion acceleration, essentially because it greatly reduces the magnitude of drift motions away from the vicinity of the X-point, where the accelerating component of the electric field is largest. Using plasma parameters consistent with flare observations, we obtain proton distributions extending up to -ray-emitting energies (>1 MeV). In some cases the energy distributions exhibit a bump-on-tail in the MeV range. In general, the shape of the distribution is sensitive to the model parameters.     

On the electric field screening by electron–positron pairs in a pulsar magnetosphere

We present a steady one-dimensional model for a pulsar polar cap accelerator, where the field-aligned electric field and flow are solved self-consistently with a given current density. It is assumed that no particles return to the star. It is known that the space-charge-limited flow is accelerated to energies high enough to create electron–positron pairs if the assumed current density is high enough. We find that when pairs are created in such a space-charge-limited flow, the accelerating electric field is screened out within a short distance after pair creation, if the pair particle flux is larger than a critical value. We also find that a space charge density wave is excited in the screening region.
We find that a pair flux larger than the critical value M _c=10³–10⁵ must be reached in a layer with thickness equal to the braking distance for the decelerating component. Therefore, the required multiplicity – the number of pairs created by one primary particle – is too large to be realized in the actual pulsar magnetosphere. We suggest that in order to obtain a localized potential drop along the polar cap magnetic flux, one needs to take into account additional effects such as wave–particle interaction or quasi-periodic pair creation.     

Test particle acceleration in torsional spine magnetic reconnection

Three-dimensional (3D) magnetic reconnection is taking place commonly in astrophysical and space plasmas, especially in solar flares which are rich sources of highly energetic particles. One of the proposed mechanisms for steady-state 3D magnetic reconnection is “torsional spine reconnection”. By using the magnetic and electric fields for “torsional spine reconnection”, we numerically investigate the features of test particle acceleration with input parameters for the solar corona. We show that efficient acceleration of a relativistic proton is possible near the null point where it can gain up to 100 MeV of kinetic energy within a few milliseconds. However, varying the injection position results in different scenarios for proton acceleration. A proton is most efficiently accelerated when it is injected at the point where the magnetic field lines change their curvature in the fan plane. Moreover, a proton injected far away from the null point cannot be accelerated and, even in some cases, it is trapped in the magnetic field. In addition, adopting either spatially uniform or non-uniform localized plasma resistivity does not much influence the features of trajectory.     

Energies of Electrons Accelerated in Turbulent Reconnecting Current Sheets in Solar Flares

Yohkoh observations strongly suggest that electron acceleration in solar flares occurs in magnetic reconnection regions in the corona above the soft X-ray flare loops. Unfortunately, models for particle acceleration in reconnecting current sheets predict electron energy gains in terms of the reconnection electric field and the thickness of the sheet, both of which are extremely difficult to measure. It can be shown, however, that application of Ohm's law in a turbulent current sheet, combined with energy and Maxwell's equations, leads to a formula for the electron energy gain in terms of the flare power output, the magnetic field strength, the plasma density and temperature in the sheet, and its area. Typical flare parameters correspond to electron energies between a few tens of keV and a few MeV. The calculation supports the viewpoint that electrons that generate the continuum gamma-ray and hard X-ray emissions in impulsive solar flares are accelerated in a large-scale turbulent current sheet above the soft X-ray flare loops.     

Gamma-Ray Emission from Pulsar Outer Magnetospheres

We investigate a stationary pair production cascade in the outer magnetosphere of an isolated, spinning neutron star. The charge depletion due to global flows of charged particles, causes a large electric field along the magnetic field lines. Migratory electrons and/or positrons are accelerated by this field to radiate gamma-rays via curvature and inverse-Compton processes. Some of such gamma-rays collide with the X-rays to materialize as pairs in the gap. The replenished charges partially screen the electric field, which is self-consistently solved together with the energy distribution of particles and gamma-rays at each point along the field lines. By solving the set of Maxwell and Boltzmann equations, we demonstrate that an external injection of charged particles at nearly Goldreich-Julian rate does not quench the gap but shifts its position and that the particle energy distribution cannot be described by a power-law. The injected particles are accelerated in the gap and escape from it with large Lorentz factors. We show that such escaping particles migrating outside of the gap contribute significantly to the gamma-ray luminosity for young pulsars and that the soft gamma-ray spectrum between 100 MeV and 3 GeV observed for the Vela pulsar can be explained by this component. We also discuss that the luminosity of the gamma-rays emitted by the escaping particles is naturally proportional to the square root of the spin-down luminosity.     

Particle Acceleration by a Time-Varying Electric Field in Merging Magnetic Fields

Traditional models for particle acceleration by magnetic reconnection in solar flares assumed a constant electric field in a steady reconnecting magnetic field. Although this assumption may be justified during the gradual phase of flares, the situation is different during the impulsive phase. Observed rapid variations in flare emissions imply that reconnection is non-steady and a time-varying electric field is present in a reconnecting current sheet. This paper describes exploratory calculations of charged particle orbits in an oscillating electric field present either at a neutral plane or a neutral line of two-dimensional magnetic field. A simple analytical model makes it possible to explain the effects of particle trapping and resonant acceleration previously noted by Petkaki and MacKinnon in a numerical simulation. As an application, electron acceleration to X-ray generating energies in impulsive solar flares is discussed within the context of the model.     

Particle Interactions with Single or Multiple 3D Solar Reconnecting Current Sheets

The acceleration of charged particles (electrons and protons) in flaring solar active regions is analyzed using numerical experiments. The acceleration is modeled as a stochastic process taking place by the interaction of the particles with local magnetic reconnection sites via multiple steps. Two types of local reconnecting topologies are studied: the Harris-type and the X-point. A formula for the maximum kinetic energy gain in a Harris-type current sheet, found in a previous work of ours, fits well the numerical data for a single step of the process. A generalization is then given approximating the kinetic energy gain through an X-point. In the case of the multiple step process, in both topologies the kinetic energy distribution of the particles is found to acquire a practically invariant form after a small number of steps. This tendency is interpreted theoretically. Other characteristics of the acceleration process are given, such as the mean acceleration time and the pitch angle distributions of the particles.     

High-energy gamma-ray emission from solar flares: Constraining the accelerated proton spectrum

Using a multi-component model to describe the γ-ray emission, we investigate the flares of December 16, 1988 and March 6, 1989 which exhibited unambiguous evidence of neutral pion decay. The observations are then combined with theoretical calculations of pion production to constrain the accelerated proton spectra. The detection of π⁰ emisson alone can indicate much about the energy distribution and spectral variation of the protons accelerated to pion producing energies. Here both the intensity and detailed spectral shape of the Doppler-broadened π⁰ decay feature are used to determine the spectral form of the accelerated proton energy distribution. The Doppler width of this γ-ray emission provides a unique diagnostic of the spectral shape at high energies, independent of any normalisation. To our knowledge, this is the first time that this diagnostic has been used to constrain the proton spectra. The form of the energetic proton distribution is found to be severely limited by the observed intensity and Doppler width of the π⁰ decay emission, demonstrating effectively the diagnostic capabilities of the π⁰ decay γ-rays. The spectral index derived from the γ-ray intensity is found to be much harder than that derived from the Doppler width. To reconcile this apparent discrepancy we investigate the effects of introducing a high-energy cut-off in the accelerated proton distribution. With cut-off energies of around 0.5–0.8 GeV and relatively hard spectra, the observed intensities and broadening can be reproduced with a single energetic proton distribution above the pion production threshold.     

Influence of the Variations of Current Sheet Parameters on the Acceleration of Electrons

By the test particle method, we have investigated the kinematic characteristics of the electrons in the reconnecting current sheet with a guiding magnetic field B_z after they are accelerated by the supper-Dreicer electric field E_z. Firstly, the influence of the guiding magnetic field B_z on the particle acceleration is discussed under the assumption that B_z is constant in magnitude but different in orientation with respect to the electric field. In this case, the variation of the B_z direction directly leads to the variation of electron trajectories and makes electrons leave the current sheet along different paths. If B_z is parallel to E_z, the pitch angles of the accelerated electrons are close to 180°. If B_z is anti-parallel to E_z, the pitch angles of the accelerated electrons are close to 0°. The orientation of the guiding magnetic field just makes the electric field accelerate selectively the electrons in different regions, but does not change the energy distribution of electrons, and the finally derived energy spectrum is the common power-law spectrum E^?γ. In typical coronal conditions, γ is about 2.9. The further study indicates that the magnitude of γ depends on the strengths of the guiding magnetic field and reconnecting electric field, as well as the scale of the current sheet. Then, the kinematic characteristics of the accelerated electrons in the current sheet with multiple X-points and O-points are also studied. The result indicates that the existences of the X-points and O-points have the particles constrained in the accelerating region to obtain the maximum acceleration, and the final energy spectrum has the characteristics of multi-power law spectra.     

The termination shock in a striped pulsar wind

The origin of radio emission from plerions is considered. Recent observations suggest that radio-emitting electrons are presently accelerated rather than having been injected at early stages of the plerion evolution. The observed flat spectra without a low-frequency cut-off imply an acceleration mechanism that raises the average particle energy by orders of magnitude but leaves most of the particles at an energy of less than approximately a few hundred MeV. It is suggested that annihilation of the alternating magnetic field at the pulsar wind termination shock provides the necessary mechanism. Toroidal stripes of opposite magnetic polarity are formed in the wind emanating from an obliquely rotating pulsar magnetosphere (the striped wind). At the termination shock, the flow compresses and the magnetic field annihilates by driven reconnection. Jump conditions are obtained for the shock in a striped wind. It is shown that the post-shock magnetohydrodynamic parameters of the flow are the same as if the energy of the alternating field had already been converted into plasma energy upstream of the shock. Therefore, the available estimates of the ratio of the Poynting flux to the matter energy flux, σ, should be attributed not to the total upstream Poynting flux but only to that associated with the average magnetic field. A simple model for the particle acceleration in the shocked striped wind is presented.     

Particle distributions and X-ray spectra in single or multiple solar current sheets

The acceleration of charged particles in a site of magnetic reconnection is analysed by detailed numerical simulations. Single or multiple encounters of the particles with Harris-type reconnecting current sheets (RCSs) are modelled as an overall stochastic process taking place within an active region. RCS physical parameters are selected in a parameter space relevant to solar flares. Initially, the charged particles form a thermal (Maxwellian) distribution corresponding to coronal temperature  ≃2 × 10⁶ K  . Our main goal is to investigate how the acceleration process changes the shape of the particles' kinetic energy distribution. The evolution of the kinetic energy distribution, calculated numerically after one encounter of the particles with a single RCS, is found to be in good agreement with our previously published analytical formulae. In the case of consecutive encounters, we find that the kinetic distribution tends to converge to a practically invariant form after a relatively small number of encounters. We construct a discrete stochastic process that reproduces the numerical distributions and we provide a theoretical interpretation of the asymptotic convergence of the energy distribution. We finally compute the theoretical X-ray spectra that would be emitted by the simulated particles in a thick target model of radiation.     

Particle Acceleration Due to Coronal Non-null Magnetic Reconnection

Various topological features, for example magnetic null points and separators, have been inferred as likely sites of magnetic reconnection and particle acceleration in the solar atmosphere. In fact, magnetic reconnection is not constrained to solely take place at or near such topological features and may also take place in the absence of such features. Studies of particle acceleration using non-topological reconnection experiments embedded in the solar atmosphere are uncommon. We aim to investigate and characterise particle behaviour in a model of magnetic reconnection which causes an arcade of solar coronal magnetic field to twist and form an erupting flux rope, crucially in the absence of any common topological features where reconnection is often thought to occur. We use a numerical scheme that evolves the gyro-averaged orbit equations of single electrons and protons in time and space, and simulate the gyromotion of particles in a fully analytical global field model. We observe and discuss how the magnetic and electric fields of the model and the initial conditions of each orbit may lead to acceleration of protons and electrons up to 2 MeV in energy (depending on model parameters). We describe the morphology of time-dependent acceleration and impact sites for each particle species and compare our findings to those recovered by topologically based studies of three-dimensional (3D) reconnection and particle acceleration. We also broadly compare aspects of our findings to general observational features typically seen during two-ribbon flare events.     

Acceleration of energetic particles of the outer regions of planetary magnetospheres: Inferences from laboratory and space experiments

High energy particles, with energies above those attainable by adiabatic or steady-state electric field acceleration, have been observed in and around the outer regions of planetary magnetospheres. Acceleration by large amplitude sporadic cross-tail electric fields over an order of magnitude greater than steady-state convection fields is proposed as a source of these particles. It is suggested that such explosive electric fields will occur intermittently in the vicinity of the tail neutral line in the expansive phases of substorms. We use laboratory Double Inverse Pinch Device (DIPD) and satellite evidence to estimate this electric potential for substorms at Earth; values of 500 kV to 2 MV are calculated, in agreement with particle observations. It is further suggested that these particles, which have been accelerated in the night side magnetosphere, drift to the dayside on closed field lines, and under certain interplanetary conditions can escape to regions upstream of the bow shock.     

The gamma ray background from large scale structure formation

Hierarchical clustering of dark matter halos is thought to describe well the large scale structure of the universe. The baryonic component of the halos is shock heated to the virial temperature while a small fraction of the energy flux through the shocks may be energized through the first order Fermi process to relativistic energy per particle. It has been proposed that the electrons accelerated in this way may upscatter the photons of the universal microwave background to gamma ray energies and indeed generate a diffuse background of gamma rays that compares well to the observations. In this paper, we calculate the spectra of the particles accelerated at the merger shocks and re-evaluate the contribution of structure formation to the extragalactic diffuse gamma ray background (EDGRB), concluding that this contribution adds up to at most 10% of the observed EDGRB.     

Breakout and Tether-Cutting Eruption Models Are Both Catastrophic (Sometimes)

We present a simplified analytic model of a quadrupolar magnetic field and flux rope to model coronal mass ejections. The model magnetic field is two-dimensional, force-free and has current only on the axis of the flux rope and within two current sheets. It is a generalization of previous models containing a single current sheet anchored to a bipolar flux distribution. Our new model can undergo quasi-static evolution either due to changes at the boundary or due to magnetic reconnection at either current sheet. We find that all three kinds of evolution can lead to a catastrophe, known as loss of equilibrium. Some equilibria can be driven to catastrophic instability either through reconnection at the lower current sheet, known as tether cutting, or through reconnection at the upper current sheet, known as breakout. Other equilibria can be destabilized through only one and not the other. Still others undergo no instability, but they evolve increasingly rapidly in response to slow steady driving (ideal or reconnective). One key feature of every case is a response to reconnection different from that found in simpler systems. In our two-current-sheet model a reconnection electric field in one current sheet causes the current in that sheet to increase rather than decrease. This suggests the possibility for the microscopic reconnection mechanism to run away.     

Formation of SCR energy spectra during stochastic acceleration with allowance for Coulomb losses

The stochastic acceleration of heavy ions by Alfvén turbulence is considered with allowance for Coulomb losses. The pattern of energy dependence of these losses gives rise to characteristic features in the energy spectra of the accelerated particles at energies of the order of several MeV nucleon^?1. The manifestation of these features in the spectra is sensitive to the temperature and density of the medium, which can serve as a basis for plasma diagnostics in the flare region. Some impulsive solar energetic particle events during which features in the spectra of ³He and ⁴He were observed are considered as an example.