This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

The expansion of hot electrons in flaring magnetic loops is crucial to understanding the dynamics of solar flares. In this paper we investigate, for the first time, the transport of hot electrons in a magnetic mirror field based on a 1-D particle-in-cell (PIC) simulation. The hot electrons with small pitch angle transport into the cold plasma, which leads to the generation of Langmuir waves in the cold plasma and ion acoustic waves in the hot plasma. The large pitch angle electrons can be confined by the magnetic mirror, resulting in the different evolution time scale between electron parallel and perpendicular temperature. This will cause the formation of electron temperature anisotropy, which can generate the whistler waves near the interface between hot electrons and cold electrons. The whistler waves can scatter the large pitch angle electrons to smaller value through the cyclotron resonance, leading to electrons escaping from the hot region. These results indicate that the whistler waves may play an important role in the transport of electrons in flaring magnetic loops. The findings from this study provide some new insights to understand the electron dynamics of solar flares.

Magnetic field are transported and tangled by turbulence, even as they lose identity due to nonideal or resistive effects. On balance field lines undergo stretch-twist-fold processes. The curvature field, a scalar that measures the tangling of the magnetic field lines, is studied in detail here, in the context of magnetohydrodynamic turbulence. A central finding is that the magnitudes of the curvature and the magnetic field are anti-correlated. High curvature co-locates with low magnetic field, which gives rise to power-law tails of the probability density function of the curvature field. The curvature drift term that converts magnetic energy into flow and thermal energy, largely depends on the curvature field behavior, a relationship that helps to explain particle acceleration due to curvature drift. This adds as well to evidence that turbulent effects most likely play an essential role in particle energization since turbulence drives stronger tangled field configurations, and therefore curvature.

In the context of space and astrophysical plasma turbulence and particle heating, several vocabularies emerge for estimating turbulent energy dissipation rate, including Kolmogorov-Yaglom third-order law and, in its various forms, $\boldsymbol{j}\cdot\boldsymbol{E}$ (work done by the electromagnetic field on particles), and $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$ (pressure-strain interaction), to name a couple. It is now understood that these energy transfer channels, to some extent, are correlated with coherent structures. In particular, we find that different energy dissipation proxies, although not point-wise correlated, are concentrated in proximity to each other, for which they decorrelate in a few $d_i$(s). However, the energy dissipation proxies dominate at different scales. For example, there is an inertial range over which the third-order law is meaningful. Contributions from scale bands stemming from scale-dependent spatial filtering show that, the energy exchange through $\boldsymbol{j}\cdot\boldsymbol{E}$ mainly results from large scales, while the energy conversion from fluid flow to internal through $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$ dominates at small scales.

In this paper, with two-dimensional (2-D) hybrid simulations, we study the generation mechanism of filamentary structures downstream of a quasi-parallel shock. The results show that in the downstream both the amplitude of magnetic field and number density exhibit obvious filamentary structures, and the magnetic field and number density are anticorrelated. Detailed analysis find that these downstream compressive waves propagate almost perpendicular to the magnetic field, and the dominant wave number is around the inverse of ion kinetic scale. Their parallel and perpendicular components roughly satisfies(where and represent the parallel and in-plane perpendicular components of magnetic field, is the wave number in the perpendicular direction, and in the ion gyroradius), and their Alfven ratio also roughly agree with the analytical relation (where and indicate the Alfven ratio and plasma beta), while the corresponding cross helicity and compressibility show good agreement with previous theoretical calculations. All these properties are consistent with those of kinetic slow waves (KSWs). Therefore, we conclude that the filamentary structures in the downstream of a quasi-parallel shock are produced due to the excitation of KSWs.

We carry out two-dimensional global particle-in-cell simulations of the interaction between the solar wind and a dipole field to study the formation of the bow shock and magnetosphere. A self-reforming bow shock ahead of a dipole field is presented by using relatively high temporal-spatial resolutions. We find that (1) the bow shock and the magnetosphere are formed and reach a quasi-stable state after several ion cyclotron periods, and (2) under the Bz southward solar wind condition the bow shock undergoes a self-reformation for low \b{eta}i and high MA. Simultaneously, a magnetic reconnection in the magnetotail is found. For high \b{eta}i and low MA, the shock becomes quasi-stationary, and the magnetotail reconnection disappears. In addition, (3) the magnetopause deflects the magnetosheath plasmas. The sheath particles injected at the quasi-perpendicular region of the bow shock can be convected to downstream of an oblique shock region. A fraction of these sheath particles can leak out from the magnetosheath at the wings of the bow shock. Hence, the downstream situation is more complicated than that for a planar shock produced in local simulations.

In this paper, two-dimensional (2-D) hybrid simulations are performed to investigate ion dynamics at a rippled quasi-parallel shock. The results show that the ripples around the shock front are inherent structures of a quasi-parallel shock, and the reformation of the shock is not synchronous along the surface of the shock front. By following the trajectories of the upstream ions, we find that these ions behave differently when they interact with the shock front at different positions along the shock surface. The upstream particles are easier to transmit through the upper part of a ripple, and the bulk velocity in the corresponding downstream is larger, where a high-speed jet is formed. In the lower part of the ripple, the upstream particles tend to be reflected by the shock. For the reflected ions by the shock, they may suffer multiple stage acceleration when moving along the shock surface, or trapped between the upstream waves and the shock front. At last, these ions may escape to the further upstream or enter the downstream, therefore, the superthermal ions can be found in both the upstream and downstream.

We present the particle-in-cell (PIC) simulation results of the interaction of a high-energy lepton plasma flow with background electron-proton plasma and focus on the acceleration processes of the protons. It is found that the acceleration follows a two-stage processes. In the first stage, protons are accelerated transversely (perpendicular to the lepton flow) by the turbulent magnetic field "islands" generated via the strong Weibel-type instabilities. The accelerated protons shows a perfect inverse-power energy spectrum. As the interaction continues, a shockwave structure forms and the protons in front of the shockwave are reflected at twice of the shock speed, resulting in a quasi-monoenergetic peak located near 200MeV under the simulation parameters.

The transition between the supersonic solar wind and the subsonic heliosheath, the termination shock (TS), was observed by Voyager 2 (V2) on 2007 August 31-September 1 at a distance of 84 AU from the Sun. The data reveal multiple crossings of a complex, quasi-perpendicular supercritical shock. These experimental data are the starting point for a more sophisticated analysis that includes computer modeling of a shock in the presence of pickup ions (PUIs). here, we present two-dimensional (2-D) particle-in-cell (PIC) simulations of the TS including PUIs self-consistently. We also report the ion velocity distribution across the TS using the Faraday cup data from V2. A relatively complete plasma and magnetic field data set from V2 gives us the opportunity to do a full comparison between the experimental data and PIC simulation results. Our results show that: (1) The nonstationarity of the shock front is mainly caused by the ripples along the shock front and these ripples from even if the percentage of PUIs is high. (2) PUIs play a key role in the energy dissipation of the TS, and most of the incident ion dynamic energy is transferred to the thermal energy of PUIs instead of solar wind ions (SWIs). (3) The simulated composite heliosheath ion velocity distribution function is a superposition of a cold core formed by transmitted SWIs, the shoulders contributed by the hot reflected SWIs and directly transmitted PUIs, and the wings of the distribution dominated by the very hot reflected PUIs. (4) The V2 Faraday cups observed the cool core of the distribution, so they saw only a tip of the iceberg. For the evolution of the cool core distribution function across the TS, the computed results agree reasonably well with the V2experimental results.

We present chromospheric activity index $S\rm_{HK}$ measurements for over 13,000 F, G and K disk stars with high signal-to-noise ratio ($>$ 60) spectra in the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) spectroscopic sample. A parameter $\delta$S is defined as the difference between $S\rm_{HK}$ and a `zero' emission line fitted by several of the most inactive stars. The $S\rm_{HK}$ indices of subgiant stars tend to be much lower than dwarfs, which provide a way to distinguish dwarfs and giants with relatively low resolution spectra. Cooler stars are generally more active and display a larger scatter than hotter stars. Stars associated with the thick disk are in general less active than those of the thin disk. The fraction of K dwarfs that are active drops with vertical distance from the Galactic plane. Metallicity affects $S\rm_{HK}$ measurements differently among F, G and K dwarfs in this sample. Using the open clusters NGC 2420, M67 and NGC6791 as calibrations, ages of most field stars in this SDSS sample range from 3-8 Gyr.

We present new quasars discovered in the vicinity of the Andromeda and Triangulum galaxies with the LAMOST during the 2010 and 2011 observational seasons. Quasar candidates are selected based on the available SDSS, KPNO 4 m telescope, XSTPS optical, and WISE near infrared photometric data. We present 509 new quasars discovered in a stripe of ~135 sq. deg from M31 to M33 along the Giant Stellar Stream in the 2011 pilot survey datasets, and also 17 new quasars discovered in an area of ~100 sq. deg that covers the central region and the southeastern halo of M31 in the 2010 commissioning datasets. These 526 new quasars have i magnitudes ranging from 15.5 to 20.0, redshifts from 0.1 to 3.2. They represent a significant increase of the number of identified quasars in the vicinity of M31 and M33. There are now 26, 62 and 139 known quasars in this region of the sky with i magnitudes brighter than 17.0, 17.5 and 18.0 respectively, of which 5, 20 and 75 are newly-discovered. These bright quasars provide an invaluable collection with which to probe the kinematics and chemistry of the ISM/IGM in the Local Group of galaxies. A total of 93 quasars are now known with locations within 2.5 deg of M31, of which 73 are newly discovered. Tens of quasars are now known to be located behind the Giant Stellar Stream, and hundreds behind the extended halo and its associated substructures of M31. The much enlarged sample of known quasars in the vicinity of M31 and M33 can potentially be utilized to construct a perfect astrometric reference frame to measure the minute PMs of M31 and M33, along with the PMs of substructures associated with the Local Group of galaxies. Those PMs are some of the most fundamental properties of the Local Group.

Reconnection of the self-generated magnetic fields in laser-plasma interaction was first investigated experimentally by Nilson {\it et al.} [Phys. Rev. Lett. 97, 255001 (2006)] by shining two laser pulses a distance apart on a solid target layer. An elongated current sheet (CS) was observed in the plasma between the two laser spots. In order to more closely model magnetotail reconnection, here two side-by-side thin target layers, instead of a single one, are used. It is found that at one end of the elongated CS a fan-like electron outflow region including three well-collimated electron jets appears. The ($>1$ MeV) tail of the jet energy distribution exhibits a power-law scaling. The enhanced electron acceleration is attributed to the intense inductive electric field in the narrow electron dominated reconnection region, as well as additional acceleration as they are trapped inside the rapidly moving plasmoid formed in and ejected from the CS. The ejection also induces a secondary CS.

We aim at investigating the formation of jet-like features in the lower solar atmosphere, e.g. chromosphere and transition region, as a result of magnetic reconnection. Magnetic reconnection as occurring at chromospheric and transition regions densities and triggered by magnetic flux emergence is studied using a 2.5D MHD code. The initial atmosphere is static and isothermal, with a temperature of 20,000 K. The initial magnetic field is uniform and vertical. Two physical environments with different magnetic field strength (25 G and 50 G) are presented. In each case, two sub-cases are discussed, where the environments have different initial mass density. In the case where we have a weaker magnetic field (25 G) and higher plasma density ($N_e=2\times 10^{11}$ cm$^{-3}$), valid for the typical quiet Sun chromosphere, a plasma jet would be observed with a temperature of 2--3 $\times 10^4$ K and a velocity as high as 40 km/s. The opposite case of a medium with a lower electron density ($N_e=2\times 10^{10}$ cm$^{-3}$), i.e. more typical for the transition region, and a stronger magnetic field of 50 G, up-flows with line-of-sight velocities as high as 90 km/s and temperatures of 6 $\times$ 10$^5$ K, i.e. upper transition region -- low coronal temperatures, are produced. Only in the latter case, the low corona Fe IX 171 \AA\ shows a response in the jet which is comparable to the O V increase. The results show that magnetic reconnection can be an efficient mechanism to drive plasma outflows in the chromosphere and transition region. The model can reproduce characteristics, such as temperature and velocity for a range of jet features like a fibril, a spicule, an hot X-ray jet or a transition region jet by changing either the magnetic field strength or the electron density, i.e. where in the atmosphere the reconnection occurs.

The luminescent properties of CsI(Na) crystals are studied in this report. By using a TDS3054C oscilloscope with a sampling frequency of 5 GS/s, we find out that nuclear recoil signals are dominated by very fast light pulse with a decay time of ~20 ns, while {\gamma}-ray signals have a decay time of ~600 ns. The wavelength of nuclear recoil and {\gamma}-ray signals are also different. The study of n/{\gamma} separation shows that the rejection factor can reach an order of 10-7 with signal efficiency more than 80% at an equivalent electron recoil energy of 20 keV or more. Such a property makes CsI(Na) an ideal candidate for dark matter searches.

We report the discovery of eight new quasars in one extragalactic field (five degree centered at RA=$08^h58^m08.2^s$, Dec=$01^o32'29.7''$) with the LAMOST commissioning observations on December 18, 2009. These quasars, with $i$ magnitudes from 16.44 to 19.34 and redshifts from 0.898 to 2.773, were not identified in the SDSS spectroscopic survey, though six of them with redshifts less than 2.5 were selected as quasar targets in SDSS. Except one source without near-IR $Y$-band data, seven of these eight new quasars meet a newly proposed quasar selection criterion involving both near-IR and optical colors. Two of them were found in the 'redshift desert' for quasars ($z$ from 2.2 to 3), indicating that the new criterion is efficient for recovering the missing quasars with similar optical colors as stars. Although LAMOST met some problems during the commissioning observations, we were still able to identify other 38 known SDSS quasars in this field, with $i$ magnitudes from 16.24 to 19.10 and redshifts from 0.297 to 4.512. Our identifications imply that a substantial fraction of quasars may be missing in the previous quasar surveys. The implication of our results to the future LAMOST quasar survey is discussed.

The redshift range from 2.2 to 3, is known as the 'redshift desert' of quasars because quasars with redshift in this range have similar optical colors as normal stars and are thus difficult to be found in optical sky surveys. A quasar candidate, SDSS J085543.40-001517.7, which was selected by a recently proposed criterion involving near-IR $Y-K$ and optical $g-z$ colors, was identified spectroscopically as a new quasar with redshift of 2.427 by the LAMOST commissioning observation in December 2009 and confirmed by the observation made with the NAOC/Xinglong 2.16m telescope in March 2010. This quasar was not targeted in the SDSS spectroscopic survey because it locates in the stellar locus of the optical color-color diagrams, while it is clearly separated from stars in the $Y-K$ vs. $g-z$ diagram. Comparing with other SDSS quasars we found this new quasar with $i$ magnitude of 16.44 is apparently the brightest one in the redshift range from 2.3 to 2.7. From the spectral properties we derived its central black hole mass as $(1.4\sim3.9) \times 10^{10} M_\odot$ and the bolometric luminosity as $3.7\times 10^{48}$ \ergs, which indicates that this new quasar is intrinsically very bright and belongs to the most luminous quasars in the universe. Our identification supports that quasars in the redshift desert can be found by the quasar selection criterion involving the near-IR colors. More missing quasars are expected to be recovered by the future LAMOST spectroscopic surveys, which is important to the study of the cosmological evolution of quasars at redshift higher than 2.2.

Magnetic reconnection in the low atmosphere, e.g. chromosphere, is investigated in various physical environments. Its implications for the origination of explosive events (small--scale jets) are discussed. A 2.5-dimensional resistive magnetohydrodynamic (MHD) model in Cartesian coordinates is used. It is found that the temperature and velocity of the outflow jets as a result of magnetic reconnection are strongly dependent on the physical environments, e.g. the magnitude of the magnetic field strength and the plasma density. If the magnetic field strength is weak and the density is high, the temperature of the jets is very low (~10,000 K) as well as its velocity (~40 km/s). However, if environments with stronger magnetic field strength (20 G) and smaller density (electron density Ne=2x10^{10} cm^{-3}) are considered, the outflow jets reach higher temperatures of up to 600,000 K and a line-of-sight velocity of up to 130 km/s which is comparable with the observational values of jet-like events.

Two scenarios of possible ion heating due to finite amplitude parallel propagating Alfv\'en waves in the solar atmosphere are investigated using a 1D test particle approach. 1. An finite amplitude Alfv\'en wave is instantly introduced into a plasma (or equivalently, new ions are instantly created). 2. New ions are constantly created. In both scenarios, ions will be picked up by the Alfv\'en wave. In case 1, the wave scatters ions in the transverse direction leading to a randomization (or heating) process. This process is complete when a phase shift of $\pm \pi$ in the ion gyro-speed is produced between particles with characteristic parallel thermal speed and particles with zero parallel speed. This corresponds to $t={\pi \over k v_{th}}$ ($k$ is the wave number and $v_{th}$ is the ion thermal speed). A ring velocity distribution can be produced for large wave amplitude. The process yields a mass-proportional heating in the transverse direction, a temperature anisotropy and a bulk flow along the background magnetic field. In case 2, continuous ion creation represents a continuing phase shift in the ion gyro-speed leading to heating. New particles are picked up by the Alfv\'en wave within one ion gyroperiod. It is speculated that the mechanism may operate in the chromosphere and active regions where transient events may generate finite amplitude Alfv\'en waves.