The basic physics and recent progresses in theoretical and particle-in-cell (PIC) simulation studies of particle acceleration in multi-island magnetic reconnection are briefly reviewed. Particle acceleration in multi-island magnetic reconnection is considered a plausible mechanism for the acceleration of energetic particles in solar flares and the solar wind. Theoretical studies have demonstrated that such a mechanism can produce the observed power-law energy distribution of energetic particles if the particle motion is sufficiently randomized in the reconnection event. However, PIC simulations seem to suggest that the first-order Fermi acceleration mechanism is unable to produce a power-law particle energy distribution function in mildly relativistic multi-island magnetic reconnections. On the other hand, while simulations of highly relativistic reconnections appear to be able to produce a power-law energy spectrum, the spectral indices obtained are generally harder than the soft power-law spectra with indices $\sim -5$ commonly observed in the solar wind and solar flare events. In addition, the plasma heating due to kinetic instabilities in 3D magnetic reconnection may "thermalize" the power-law particles, making it even more difficult for multi-island reconnections to generate a power-law spectrum. We discuss the possible reasons that may lead to these problems.

Magnetic reconnection (MR) is considered as a major source of particle energization in astrophysical plasma. In the past, analysis of MR often assumes the magnetostatic condition, i.e. $\partial_t = 0$. We show that under the Sweet-Parker-Petschek framework, steady state is an over-constraint and is not achievable. On the other hand, the quasi-steady state defined as $\partial_t \mathbf{E} = 0$ but $\partial_t \mathbf{B} \neq 0$ or equivalently $\partial_t\mathbf{j}\neq 0$ better describes the asymptotic behaviour of MR without turbulence. The upper limit of MR rate for quasi-steady MR is found to be $\sim 1/3\sqrt{3} \sim 0.19$. The limit does not apply to impulsive or turbulent MR of which $\partial_t\mathbf{B} \neq 0$ and $\partial_t\mathbf{E} \neq 0$. In impulsive MR the rate can be higher or lower than 0.19 depending on the state of the turbulence. Our results may explain the apparent discrepancy in observations of solar flare MR rates. The analysis is independent of mass ratio and thus the results are applicable to pair plasma.

Observations of the solar corona and the solar wind discover that the solar wind is unsteady and originates from the impulsive events near the surface of the Sun's atmosphere. How solar coronal activities affect the properties of the solar wind is a fundamental issue in heliophysics. We report a simulation and theoretical investigation of how nanoflare accelerated electron beams affect the kinetic-scale properties of the solar wind and generate coherent radio emission. We show that nanoflare-accelerated electron beams can trigger a nonlinear electron two stream instability, which generates kinetic Alfv\'en and whistler waves, as well as a non-Maxwellian electron velocity distribution function, consistent with observations of the solar wind. The plasma coherent emission produced in our model agrees well with the observations of Type III, J and V solar radio bursts. Open questions in the kinetic solar wind model are also discussed.

Whether turbulence induced anomalous resistivity (AR) can facilitate a fast magnetic reconnection in collisionless plasma is a subject of active debate for decades. Recent space observations suggest that the reconnection rate can be higher than the Hall-reconnection rate and turbulent dissipation is required. In this paper, using particle-in-cell simulations, we present a case study of how AR produced by Buneman instability accelerates magnetic reconnection. We first show that the AR/drag produced by Buneman instability in a thin electron current layer 1) can dissipate magnetic energy stored in the current layer through dissipation of the kinetic energy of electron beams; 2) The inhomogeneous drag caused by wave couplings spontaneously breaks the magnetic field lines and causes impulsive fast non-Hall magnetic reconnection on electron-scales with a mean rate reaching 0.6 $V_A$. We then show that a Buneman instability driven by intense electron beams around the x-point in a 3D magnetic reconnection significantly enhances the dissipation of the magnetic energy. Electron-scale magnetic reconnections driven by the inhomogeneous drag around the x-line enhances the reconnection electric field and the in-plane perpendicular magnetic field. About 40\% of the released magnetic energy is converted into the electron thermal energy by AR while 50\% is converted into the kinetic energy of the electron beams through the acceleration by the reconnection electric field. The enhanced magnetic energy dissipation is balanced by a net Poynting flux in-flow. About 10\% of the released magnetic energy is brought out by an enhanced Poynting flux out-flow. These results suggest that AR with sufficient intensity and electron-scale inhomogeneity can significantly accelerate magnetic reconnection.

Continuous plasma coherent emission is maintained by repetitive Langmuir collapse driven by the nonlinear evolution of a strong electron two-stream instability. The Langmuir waves are modulated by solitary waves in the linear stage, and by electrostatic whistler waves in the nonlinear stage. Modulational instability leads to Langmuir collapse and electron heating that fills in cavitons. The high pressure is released via excitation of a short wavelength ion acoustic mode that is damped by electrons and that re-excites small-scale Langmuir waves---this process closes a feedback loop that maintains the continuous coherent emission.

It is well known that electron beams accelerated in solar flares can drive two-stream instability and produce radio bursts in the solar corona as well as in the interplanetary medium. Recent observations show that the solar wind likely originates from nanoflare-like events near the surface of the Sun where locally heated plasma escapes along open field lines into space. Recent numerical simulations and theoretical studies show that electron two-stream instability (ETSI) driven by nanoflare-accelerated electron beams can produce the observed nanoflare-type radio bursts, the non-Maxwellian electron velocity distribution function of the solar wind, and the kinetic scale turbulence in solar wind. This brief review focus on the basic theoretical framework and recent progress in the nonlinear evolution of ETSI, including the formation of electron holes, Langmuir wave generation in warm plasma, and the nonlinear modulation instability and Langmuir collapse. Potential applications in heliophysics and astrophysics are discussed.

We summarize our recent studies on the origin of solar wind kinetic scale turbulence and electron halo in the electron velocity distribution function. Increasing observations of nanoflares and microscopic type III radio bursts strongly suggest that nanoflares and accelerated electron beams are common in the corona. Based on particle-in-cell simulations, we show that both the core-halo feature and kinetic scale turbulence observed in the solar wind can be produced by the nonlinear evolution of electron two-stream instability driven by nanoflare accelerated electron beams. The energy exchange between waves and particles reaches equilibrium in the inner corona and the key features of the turbulence and velocity distribution are preserved as the solar wind escapes into interplanetary space along open magnetic field lines. Observational tests of the model and future theoretical work are discussed.

The formation of the observed core-halo feature in the solar wind electron velocity distribution function is a long-time puzzle. In this letter based on the current knowledge of nanoflares we show that the nanoflare-accelerated electron beams are likely to trigger a strong electron two-stream instability that generates kinetic Alfv\'en wave and whistler wave turbulence, as we demonstrated in a previous paper. We further show that the core-halo feature produced during the origin of kinetic turbulence is likely to originate in the inner corona and can be preserved as the solar wind escapes to space along open field lines. We formulate a set of equations to describe the heating processes observed in the simulation and show that the core-halo temperature ratio of the solar wind is insensitive to the initial conditions in the corona and is related to the core-halo density ratio of the solar wind and to the quasi-saturation property of the two-stream instability at the time when the exponential decay ends. This relation can be extended to the more general core-halo-strahl feature in the solar wind. The temperature ratio between the core and hot components is nearly independent of the heliospheric distance to the sun. We show that the core-halo relative drift previously reported is a relic of the fully saturated two stream instability. Our theoretical results are consistent with the observations while new tests for this model are provided.

The observed sub-proton scale turbulence spectrum in the solar wind raises the question of how that turbulence originates. Observations of keV energetic electrons during solar quite-time suggest them as possible source of free energy to drive the turbulence. Using particle-in-cell simulations, we explore how free energy in energetic electrons, released by an electron two-stream instability drives Weibel-like electromagnetic waves that excite wave-wave interactions. Consequently, both kinetic Alfv\'enic and whistler waves are excited that evolve through inverse and forward magnetic energy cascades.

The nonlinear development of the strong Buneman instability and the associated fast electron heating in thin current layers with $\Omega_e/\omega_{pe} <1$ are explored. Phase mixing of the electrons in wave potential troughs and a rapid increase in temperature are observed during the saturation of the instability. We show that the motion of trapped electrons can be described using a Hamiltonian formalism in the adiabatic approximation. The process of separatrix crossing as electrons are trapped and de-trapped is irreversible and guarantees that the resulting electron energy gain is a true heating process.

Buneman instability is often driven in magnetic reconnection. Understanding how velocity shear in the beams driving the Buneman instability affects the growth and saturation of waves is relevant to turbulence, heating, and diffusion in magnetic reconnection. Using a Mathieu-equation analysis for weak cosine velocity shear together with Vlasov simulations, the effects of shear on the kinetic Buneman instability are studied in a plasma consisting of strongly magnetized electrons and cold unmagnetized ions. In the linearly unstable phase, shear enhances the coupling between oblique waves and the sheared electron beam, resulting in a wider range of unstable eigenmodes with common lower growth rates. The wave couplings generate new features of the electric fields in space, which can persist into the nonlinear phase when electron holes form. Lower hybrid instabilities simultaneously occur at $k_{\shortparallel}/k_{\perp} \sim \sqrt{m_e/m_i}$ with a much lower growth rate, and are not affected by the velocity shear.

Using particle-in-cell simulations and kinetic theory, we explore the current-driven turbulence and associated electron heating in the dissipation region during 3D magnetic reconnection with a guide field. At late time the turbulence is dominated by the Buneman and lower hybrid instabilities. Both produce electron holes that co-exist but have very different propagation speeds. The associated scattering of electrons by the holes enhances electron heating in the dissipation region.

The nonlinear development of streaming instabilities in the current layers formed during magnetic reconnection with a guide field is explored. Theory and 3-D particle-in-cell simulations reveal two distinct phases. First, the parallel Buneman instability grows and traps low velocity electrons. The remaining electrons then drive two forms of turbulence: the parallel electron-electron two-stream instability and the nearly-perpendicular lower hybrid instability. The high velocity electrons resonate with the turbulence and transfer momentum to the ions and low velocity electrons.

If Gamma-ray Bursts (GRBs) are generally associated with supernovae like 1998bw, a relatively wide intrinsic luminosity function is implied, which indicates the existence of a large amount of undetected dim bursts, and a much higher event-rate than is often assumed. If it is assumed that the intrinsic luminosity function of GRBs is a power-law: ${\phi (L)} \propto {L^{-\beta}}$ ($\beta > 0$, $L_{min} \leq L \leq L_{max}$), data from the BATSE 4B catalog can be used to constrain slope index $\beta$ and the dynamic range width Log ${L_{max}\over L_{min}}$. Using a K-S test comparison with the observational Log $N$ - Log $P$, we find constraints on the GRB fireball model, GRB jets, and the possible GRB contribution to cosmic gamma-ray background. We find the acceptable dynamic range for $10^2<L_{max}/L_{min} < 10^7$.Our results show that jet model is more likely to be related more highly energetic explosion than fireball model.Our studies also show that the luminosity function provided by a purely special relativistic effect on a jet is outside of the K-S test acceptable range. Were intrinsic beaming to confine the jet to $\theta_{max} < 1/\gamma$ however, the effects of relativistic beaming would not dominate.

Recent optical afterglow observations of gamma-ray bursts indicate a setting and distance scale that many relate to star-formation regions. In this paper, we use $<V/V_{max}>$ and a set of artificial trigger thresholds to probe several potential GRB source density evolutionary scenarios. In particular, we compare a uniform subset of BATSE 4B data to cosmological scenarios where GRBs evolve as the comoving density, the star formation rate, the QSO rate, and the SN Type Ic rate. Standard candle bursts with power-law spectra and a universe without vacuum energy were assumed. Our results significantly favor a comoving density model, implying that GRB source density evolution is weaker than expected in these evolutionary scenarios. GRB density might still follow star-formation rates given proper concurrent GRB luminosity evolution, significant beaming, significant error in standard candle assumptions, or were a significant modification of star formation rate estimates to occur.