Longitudinal oscillations of solar filament have been investigated via numerical simulations continuously, but mainly in one dimension (1D), where the magnetic field line is treated as a rigid flux tube. Whereas those one-dimensional simulations can roughly reproduce the observed oscillation periods, implying that gravity is the main restoring force for filament longitudinal oscillations, the decay time in one-dimensional simulations is generally longer than in observations. In this paper, we perform a two-dimensional (2D) non-adiabatic magnetohydrodynamic simulation of filament longitudinal oscillations, and compare it with the 2D adiabatic case and 1D adiabatic and non-adiabatic cases. It is found that, whereas both non-adiabatic processes (radiation and heat conduction) can significantly reduce the decay time, wave leakage is another important mechanism to dissipate the kinetic energy of the oscillating filament when the magnetic field is weak so that gravity is comparable to Lorentz force. In this case, our simulations indicate that the pendulum model might lead to an error of ~100% in determining the curvature radius of the dipped magnetic field using the longitudinal oscillation period when the gravity to Lorentz force ratio is close to unity.

Magnetic flux ropes play a key role in triggering solar flares in the solar atmosphere. In this paper, we investigate the evolution of active region NOAA 12268 within 36 hours from 2015 January 29 to 30, during which a flux rope was formed and three M-class and three C-class flares were triggered without coronal mass ejections. During the evolution of the active region, the flare emission seen in the H$\alpha$ and ultraviolet wavebands changed from a circular shape (plus an adjacent conjugated ribbon and a remote ribbon) to three relatively straight and parallel ribbons. Based on a series of reconstructed nonlinear force-free fields, we find sheared or twisted magnetic field lines and a large-scale quasi-separatrix layer (QSL) associated with 3D null points in a quadrupolar magnetic field. These features always existed and constantly evolved during the two days. The twist of the flux rope was gradually accumulated that eventually led to its instability. Around the flux rope, there were some topological structures, including a bald patch, a hyperbolic flux tube and a torus QSL. We discuss how the particular magnetic structure and its evolution produce the flare emission. In particular, the bifurcation of the flux rope can explain the transition of the flares from circular to parallel ribbons. We propose a two-stage evolution of the magnetic structure and its associated flares. In the first stage, sheared arcades under the dome-like large-scale QSL were gradually transformed into a flux rope through magnetic reconnection, which produced the circular ribbon flare. In the second stage, the flux rope bifurcated to form the three relatively straight and parallel flare ribbons.

With the rapid development of telescopes, both temporal cadence and the spatial resolution of observations are increasing. This in turn generates vast amount of data, which can be efficiently searched only with automated detections in order to derive the features of interest in the observations. A number of automated detection methods and algorithms have been developed for solar activities, based on the image processing and machine learning techniques. In this paper, after briefly reviewing some automated detection methods, we describe our efficient and versatile automated detection method for solar filaments. It is able not only to recognize filaments, determine the features such as the position, area, spine, and other relevant parameters, but also to trace the daily evolution of the filaments. It is applied to process the full disk H-alpha data observed in nearly three solar cycles, and some statistic results are presented.

Solar prominences are subject to all kinds of perturbations during their lifetime, and frequently demonstrate oscillations. The study of prominence oscillations provides an alternative way to investigate their internal magnetic and thermal structures as the oscillation characteristics depend on their interplay with the solar corona. Prominence oscillations can be classified into longitudinal and transverse types. We perform three-dimensional ideal magnetohydrodynamic simulations of prominence oscillations along a magnetic flux rope, with the aim to compare the oscillation periods with those predicted by various simplified models and to examine the restoring force. We find that the longitudinal oscillation has a period of about 49 minutes, which is in accordance with the pendulum model where the field-ligned component of gravity serves as the restoring force. In contrast, the horizontal transverse oscillation has a period of about 10 minutes and the vertical transverse oscillation has a period of about 14 minutes, and both of them can be nicely fitted with a two-dimensional slab model. We also find that the magnetic tension force dominates most of the time in transverse oscillations, except for the first minute when magnetic pressure overwhelms.

White-light flares are the flares with emissions visible in the optical continuum. They are thought to be rare and pose the most stringent requirements in energy transport and heating in the lower atmosphere. Here we present a nearly circular white-light flare on 2015 March 10 that was well observed by the Optical and Near-infrared Solar Eruption Tracer and Solar Dynamics Observatory. In this flare, there appear simultaneously both impulsive and gradual white-light kernels. The generally accepted thick-target model would be responsible for the impulsive kernels but not sufficient to interpret the gradual kernels. Some other mechanisms including soft X-ray backwarming or downward-propagating Alfven waves, acting jointly with electron beam bombardment, provide a possible interpretation. However, the origin of this kind of white-light kernels is still an open question that induces more observations and researches in the future to decipher it.

The large-amplitude longitudinal oscillations of solar filaments have been observed and explored for more than ten years. Previous studies are mainly based on the one-dimensional rigid flux tube model with a single magnetic dip. However, it is noticed that there might be two magnetic dips, and hence two threads, along one magnetic field line. Following the previous work, we intend to investigate the kinematics of the filament longitudinal oscillations when two threads are magnetically connected, which is done by solving one-dimensional radiative hydrodynamic equations with the numerical code MPI-AMRVAC. Two different types of perturbations are considered, and the difference from previous works resulting from the filament thread-thread interaction is investigated. We find that even with the inclusion of the thread-thread interaction, the oscillation period is modified weakly, by at most 20% compared to the traditional pendulum model with one thread. However, the damping timescale is significantly affected by the thread-thread interaction. Hence, we should take it into account when applying the consistent seismology to the filaments where two threads are magnetically connected.

Ellerman bombs (EBs) are small brightening events in the solar lower atmosphere. By original definition, the main EB's characteristic is the two emission bumps in both wings of chromospheric lines, such as H$\alpha$ and Ca II 8542 {\AA} lines. Up to now, most authors found that the temperature increase of EBs around the temperature minimum region is in the range of 600K-3000K. However, with recent IRIS observations, some authors proposed that the temperature increase of EBs could be more than 10000K. Using non-LTE semi-empirical modeling, we investigate the line profiles, continuum emission and the radiative losses for the EB models with different temperature increases, and compare them with observations. Our result indicates that if the EB maximum temperature attains more than 10000K around the temperature minimum region, then the resulted H$\alpha$ and Ca II 8542 {\AA} line profiles and the continuum emission would be much stronger than that of EB observations. Moreover, due to the high radiative losses, the high temperature EB would have a very short lifetime, which is not comparable with the observations. Thus, our study does not support the proposal that the EB temperatures are higher than 10000K.

It is important to study the fine structures of solar filaments with high-resolution observations since it can help us understand the magnetic and thermal structures of the filaments and their dynamics. In this paper, we study a newly-formed filament located inside the active region NOAA 11762, which was observed by the 1.6 m New Solar Telescope (NST) at Big Bear Solar Observatory (BBSO) from 16:40:19 UT to 17:07:58 UT on 2013 June 5. As revealed by the H$\alpha$ filtergrams, cool material is seen to be injected into the filament spine with a speed of 5--10 km s$^{-1}$. At the source of the injection, brightenings are identified in the chromosphere, which is accompanied by magnetic cancellation in the photosphere, implying the importance of magnetic reconnection in replenishing the filament with plasmas from the lower atmosphere. Counter-streamings are detected near one endpoint of the filament, with the plane-of-the-sky speed being 7--9 km s$^{-1}$ in the H$\alpha$ red-wing filtergrams and 9--25 km s$^{-1}$ in the blue-wing filtergrams. The observations are indicative of that this active region filament is supported by a sheared arcade without magnetic dips, and the counter-streamings are due to unidirectional flows with alternative directions, rather than due to the longitudinal oscillations of filament threads as in many other filaments.

Solar filaments can be formed via chromospheric evaporation followed by condensation in the corona or by the direct injection of cool plasma from the chromosphere to the corona. In this paper, with high-resolution H$\alpha$ data observed by the 1.6 m New Solar Telescope of the Big Bear Solar Observatory on 2015 August 21, we confirmed that an active-region filament is maintained by the continuous injection of cold chromospheric plasma. We find that the filament is rooted along a bright ridge in H$\alpha$, which corresponds to the intersection of a magnetic quasi-separatrix layer with the solar surface. This bright ridge consists of many small patches and the sizes of these patches are comparable to the width of the filament threads. It is found that upflows originate from the brighter patches of the ridge, whereas the downflows move toward the weaker patches of the ridge. The whole filament is composed of two opposite directional streams, implying that longitudinal oscillations are not the only cause of the counterstreamings, and unidirectional siphon flows with alternative directions are another possibility.

It has been revealed that the magnetic topology in the solar atmosphere displays hemispheric preference, i.e., negative/positive helicity in the northern/southern hemisphere, respectively. However, the strength of the hemispheric rule and its cyclic variation are controversial. In this paper, we apply a new method based on filament drainage to 571 erupting filaments from 2010 May to 2015 December in order to determine the filament chirality and its hemispheric preference. It is found that $91.6\%$ of our sample of erupting filaments follow the hemispheric rule of helicity sign. It is found that the strength of the hemispheric preference of the quiescent filaments decreases slightly from $\sim$$97\%$ in the rising phase to $\sim$$85\%$ in the declining phase of solar cycle 24, whereas the strength of the intermediate filaments keeps a high value around $96\pm4\%$ all the time. Only the active-region filaments show significant variations. Their strength of the hemispheric rule rises from $\sim$$63\%$ to $\sim$$95\%$ in the rising phase, and keeps a high value $82\pm5\%$ during the declining phase. Furthermore, during a half-year period around the solar maximum, their hemispheric preference totally vanishes. Besides, we also diagnose the magnetic configurations of the filaments based on our indirect method, and found that in our sample of erupting events, $89\%$ are inverse-polarity filaments with a flux rope magnetic configuration, whereas $11\%$ are normal-polarity filaments with a sheared arcade configuration.

The variations of the 3-D coronal magnetic fields associated with the X3.4-class flare of active region 10930 are studied in this paper. The coronal magnetic field data are reconstructed from the photospheric vector magnetograms obtained by the Hinode satellite and using the nonlinear force-free field extrapolation method developed in our previous work (He et al., 2011). The 3-D force-free factor $\alpha$, 3-D current density, and 3-D magnetic energy density are employed to analyze the coronal data. The distributions of $\alpha$ and current density reveal a prominent magnetic connectivity with strong negative $\alpha$ values and strong current density before the flare. This magnetic connectivity extends along the main polarity inversion line and is found to be totally broken after the flare. The distribution variation of magnetic energy density reveals the redistribution of magnetic energy before and after the flare. In the lower space of the modeling volume the increase of magnetic energy dominates, and in the higher space the decrease of energy dominates. The comparison with the flare onset imaging observation exhibits that the breaking site of the magnetic connectivity and site with the highest values of energy density increase coincide with the location of flare initial eruption. We conclude that a cramped positive $\alpha$ region appearing in the photosphere causes the breaking of the magnetic connectivity. A scenario for flare initial eruption is proposed in which the Lorentz force acting on the isolated electric current at the magnetic connectivity breaking site lifts the associated plasmas and causes the initial ejection.

The discovery of stationary "EIT waves" about 16 years ago posed a big challenge to the then favorite fast-mode wave model for coronal "EIT waves". It encouraged the proposing of various non-wave models, and played an important role in approaching the recent converging viewpoint, {\it i.e.} there are two types of EUV waves. However, it was recently discovered that a stationary wave front can also be generated when a fast-mode wave passes through a magnetic quasi-separatrix layer (QSL). In this paper, we perform a magnetohydrodynamic (MHD) numerical simulation of the interaction between a fast-mode wave and a magnetic QSL, and a stationary wave front is reproduced. The analysis of the numerical results indicates that near the plasma beta $\sim 1$ layer in front of the magnetic QSL, part of the fast-mode wave is converted to a slow-mode MHD wave, which is then trapped inside the magnetic loops, forming a stationary wave front. Our research implies that we have to be cautious in identifying the nature of a wave since there may be mode conversion during the propagation of the waves driven by solar eruptions.

Magnetic reconnection is a rapid energy release process that is believed to be responsible for flares on the Sun and stars. Nevertheless, such flare-related reconnection is mostly detected to occur in the corona, while there have been few studies concerning the reconnection in the chromosphere or photosphere. Here we present both spectroscopic and imaging observations of magnetic reconnection in the chromosphere leading to a microflare. During the flare peak time, chromospheric line profiles show significant blueshifted/redshifted components on the two sides of the flaring site, corresponding to upflows and downflows with velocities of $\pm$(70--80) km s$^{-1}$, comparable with the local Alfv\'{e}n speed as expected by the reconnection in the chromosphere. The three-dimensional nonlinear force-free field configuration further discloses twisted field lines (a flux rope) at a low altitude, cospatial with the dark threads in He I 10830 \r{A} images. The instability of the flux rope may initiate the flare-related reconnection. These observations provide clear evidence of magnetic reconnection in the chromosphere and show the similar mechanisms of a microflare to those of major flares.

We improve our filament automated detection method which was proposed in our previous works. It is then applied to process the full disk H$\alpha$ data mainly obtained by Big Bear Solar Observatory (BBSO) from 1988 to 2013, spanning nearly 3 solar cycles. The butterfly diagrams of the filaments, showing the information of the filament area, spine length, tilt angle, and the barb number, are obtained. The variations of these features with the calendar year and the latitude band are analyzed. The drift velocities of the filaments in different latitude bands are calculated and studied. We also investigate the north-south (N-S) asymmetries of the filament numbers in total and in each subclass classified according to the filament area, spine length, and tilt angle. The latitudinal distribution of the filament number is found to be bimodal. About 80% of all the filaments have tilt angles within [0{\deg}, 60{\deg}]. For the filaments within latitudes lower (higher) than 50{\deg} the northeast (northwest) direction is dominant in the northern hemisphere and the southeast (southwest) direction is dominant in the southern hemisphere. The latitudinal migrations of the filaments experience three stages with declining drift velocities in each of solar cycles 22 and 23, and it seems that the drift velocity is faster in shorter solar cycles. Most filaments in latitudes lower (higher) than 50{\deg} migrate toward the equator (polar region). The N-S asymmetry indices indicate that the southern hemisphere is the dominant hemisphere in solar cycle 22 and the northern hemisphere is the dominant one in solar cycle 23.

Two-ribbon brightenings are one of the most remarkable characteristics of an eruptive solar flare and are often used for predicting the occurrence of coronal mass ejections (CMEs). Nevertheless, it was called in question recently whether all two-ribbon flares are eruptive. In this paper, we investigate a two ribbon-like white-light (WL) flare that is associated with a failed magnetic flux rope (MFR) eruption on 2015 January 13, which has no accompanying CME in the WL coronagraph. Observations by \textit{Optical and Near-infrared Solar Eruption Tracer} and \textit{Solar Dynamics Observatory} reveal that, with the increase of the flare emission and the acceleration of the unsuccessfully erupting MFR, two isolated kernels appear at the WL 3600 {\AA} passband and quickly develop into two elongated ribbon-like structures. The evolution of the WL continuum enhancement is completely coincident in time with the variation of \textit{Fermi} hard X-ray 26--50 keV flux. Increase of continuum emission is also clearly visible at the whole FUV and NUV passbands observed by \textit{Interface Region Imaging Spectrograph}. Moreover, in one WL kernel, the \ion{Si}{4}, \ion{C}{2}, and \ion{Mg}{2} h/k lines display significant enhancement and non-thermal broadening. However, their Doppler velocity pattern is location-dependent. At the strongly bright pixels, these lines exhibit a blueshift; while at moderately bright ones, the lines are generally redshifted. These results show that the failed MFR eruption is also able to produce a two-ribbon flare and high-energy electrons that heat the lower atmosphere, causing the enhancement of the WL and FUV/NUV continuum emissions and chromospheric evaporation.

We attempt to propose a method for automatically detecting the solar filament chirality and barb bearing. We first introduce the unweighted undirected graph concept and adopt the Dijkstra shortest-path algorithm to recognize the filament spine. Then, we use the polarity inversion line (PIL) shift method for measuring the polarities on both sides of the filament, and employ the connected components labeling method to identify the barbs and calculate the angle between each barb and the spine to determine the bearing of the barbs, i.e., left or right. We test the automatic detection method with H-alpha filtergrams from the Big Bear Solar Observatory (BBSO) H-alpha archive and magnetograms observed with the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). Four filaments are automatically detected and illustrated to show the results. The barbs in different parts of a filament may have opposite bearings. The filaments in the southern hemisphere (northern hemisphere) mainly have left-bearing (right-bearing) barbs and positive (negative) magnetic helicity, respectively. The tested results demonstrate that our method is efficient and effective in detecting the bearing of filament barbs. It is demonstrated that the conventionally believed one-to-one correspondence between filament chirality and barb bearing is not valid. The correct detection of the filament axis chirality should be done by combining both imaging morphology and magnetic field observations.

Magnetic reconnection, a change of magnetic field connectivity, is a fundamental physical process in which magnetic energy is released explosively. It is responsible for various eruptive phenomena in the universe. However, this process is difficult to observe directly. Here, the magnetic topology associated with a solar reconnection event is studied in three dimensions (3D) using the combined perspectives of two spacecraft. The sequence of extreme ultraviolet (EUV) images clearly shows that two groups of oppositely directed and non-coplanar magnetic loops gradually approach each other, forming a separator or quasi-separator and then reconnecting. The plasma near the reconnection site is subsequently heated from $\sim$1 to $\ge$5 MK. Shortly afterwards, warm flare loops ($\sim$3 MK) appear underneath the hot plasma. Other observational signatures of reconnection, including plasma inflows and downflows, are unambiguously revealed and quantitatively measured. These observations provide direct evidence of magnetic reconnection in a 3D configuration and reveal its origin.

Ellerman bombs (EBs) are tiny brightenings often observed near sunspots. The most impressive characteristic of the EB spectra is the two emission bumps in both wings of the H$\alpha$ and \ion{Ca}{II} 8542 {\AA} lines. High-resolution spectral data of three small EBs were obtained on 2013 June 6 with the largest solar telescope, the 1.6 meter New Solar Telescope (NST), at the Big Bear Solar Observatory. The characteristics of these EBs are analyzed. The sizes of the EBs are in the range of 0.3\arcsec\--0.8\arcsec\ and their durations are only 3--5 minutes. Our semi-empirical atmospheric models indicate that the heating occurs around the temperature minimum region with a temperature increase of 2700--3000 K, which is surprisingly higher than previously thought. The radiative and kinetic energies are estimated to be as high as 5$\times$10$^{25}$--3.0$\times$10$^{26}$ ergs despite the small size of these EBs. Observations of the magnetic field show that the EBs appeared just in a parasitic region with mixed polarities and accompanied by mass motions. Nonlinear force-free field extrapolation reveals that the three EBs are connected with a series of magnetic field lines associated with bald patches, which strongly implies that these EBs should be produced by magnetic reconnection in the solar lower atmosphere. According to the lightcurves and the estimated magnetic reconnection rate, we propose that there is a three phase process in EBs: pre-heating, flaring and cooling phases.

Helical magnetic flux rope (MFR) is a fundamental structure of corona mass ejections (CMEs) and has been discovered recently to exist as a sigmoidal channel structure prior to its eruption in the extreme ultraviolet (EUV) high temperature passbands of the Atmospheric Imaging Assembly (AIA). However, when and where the MFR is built up are still elusive. In this paper, we investigate two MFRs (MFR1 and MFR2) in detail, whose eruptions produced two energetic solar flares and CMEs on 2014 April 18 and 2014 September 10, respectively. The AIA EUV images reveal that for a long time prior to their eruption, both MFR1 and MFR2 are under formation, which is probably through magnetic reconnection between two groups of sheared arcades driven by the shearing and converging flows in the photosphere near the polarity inversion line. At the footpoints of the MFR1, the \textit{Interface Region Imaging Spectrograph} Si IV, C II, and Mg II lines exhibit weak to moderate redshifts and a non-thermal broadening in the pre-flare phase. However, a relatively large blueshift and an extremely strong non-thermal broadening are found at the formation site of the MFR2. These spectral features consolidate the proposition that the reconnection plays an important role in the formation of MFRs. For the MFR1, the reconnection outflow may propagate along its legs, penetrating into the transition region and the chromosphere at the footpoints. For the MFR2, the reconnection probably takes place in the lower atmosphere and results in the strong blueshift and non-thermal broadening for the Mg II, C II, and Si IV lines.

We present a study of the waiting time distributions (WTDs) of solar energetic particle (SEP) events observed with the spacecraft $WIND$ and $GOES$. Both the WTDs of solar electron events (SEEs) and solar proton events (SPEs) display a power-law tail $\sim \Delta t^{-\gamma}$. The SEEs display a broken power-law WTD. The power-law index is $\gamma_{1} =$ 0.99 for the short waiting times ($<$70 hours) and $\gamma_{2} =$ 1.92 for large waiting times ($>$100 hours). The break of the WTD of SEEs is probably due to the modulation of the corotating interaction regions (CIRs). The power-law index $\gamma \sim$ 1.82 is derived for the WTD of SPEs that is consistent with the WTD of type II radio bursts, indicating a close relationship between the shock wave and the production of energetic protons. The WTDs of SEP events can be modeled with a non-stationary Poisson process which was proposed to understand the waiting time statistics of solar flares (Wheatland 2000; Aschwanden $\&$ McTiernan 2010). We generalize the method and find that, if the SEP event rate $\lambda = 1/\Delta t$ varies as the time distribution of event rate $f(\lambda) = A \lambda^{-\alpha}exp(-\beta \lambda)$, the time-dependent Poisson distribution can produce a power-law tail WTD $\sim \Delta t^{\alpha - 3}$, where $0 \leq \alpha < 2$.

We present a study of 7 large solar proton events (SPEs) of current solar cycle 24 (from 2009 January up to date). They were recorded by GOES spacecraft with highest proton fluxes over 200 pfu for energies $>$10 MeV. In situ particle measurements show that: (1) The profiles of the proton fluxes are highly dependent of the locations of their solar sources, namely flares or coronal mass ejections (CMEs); (2) The solar particle release (SPR) times fall in the decay phase of the flare emission, and are in accordance with the times when the CMEs travel to an average height of 7.9 solar radii; (3) The time differences between the SPR and the flare peak are also dependent of the locations of the solar active regions (ARs). The results tend to support the concept of proton acceleration by the CME-driven shock, even though there exists a possibility of particle acceleration at flare site with subsequent perpendicular diffusion of accelerated particles in the interplanetary magnetic field (IMF). We derive the integral time-of-maximum (TOM) spectra of solar protons in two forms: a single power-law distribution and a power law broken with an exponential tail. It is found that the unique Ground Level Enhancement (GLE) event on 2012 May 17 displays a hardest spectrum and a largest broken energy that may explain why the this event could extend to relativistic energy.

We study the H$\alpha$ and Ca II 8542 \r{A} line spectra of four typical Ellerman bombs (EBs) in active region NOAA 11765 on 2013 June 6, observed with the Fast Imaging Solar Spectrograph installed at the 1.6 meter New Solar Telescope at Big Bear Solar Observatory. Considering that EBs may occur in a restricted region in the lower atmosphere, and that their spectral lines show particular features, we propose a two-cloud model to fit the observed line profiles. The lower cloud can account for the wing emission, and the upper cloud is mainly responsible for the absorption at line center. After choosing carefully the free parameters, we get satisfactory fitting results. As expected, the lower cloud shows an increase of the source function, corresponding to a temperature increase of 400--1000 K in EBs relative to the quiet Sun. This is consistent with previous results deduced from semi-empirical models and confirms that a local heating occurs in the lower atmosphere during the appearance of EBs. We also find that the optical depths can increase to some extent in both the lower and upper clouds, which may result from either a direct heating in the lower cloud, or illumination by an enhanced radiation on the upper cloud. The velocities derived from this method, however, are different from those obtained using the traditional bisector method, implying that one should be cautious when interpreting this parameter. The two-cloud model can thus be used as an efficient method to deduce the basic physical parameters of EBs.

The dynamics of a filament channel are observed with imaging and spectroscopic telescopes before and during the filament eruption on 2011 January 29. The extreme ultraviolet (EUV) spectral observations reveal that there is no EUV counterparts of the H$\alpha$ counter-streamings in the filament channel, implying that the ubiquitous H$\alpha$ counter-streamings found by previous research are mainly due to longitudinal oscillations of filament threads, which are not in phase between each other. However, there exist larger-scale patchy counter-streamings in EUV along the filament channel from one polarity to the other, implying that there is another component of uni-directional flow (in the range of $\pm$10 km s$^{-1}$) inside each filament thread in addition to the implied longitudinal oscillation. Our results suggest that the flow direction of the larger-scale patchy counter-streaming plasma in the EUV is related to the intensity of the plage or active network, with the upflows being located at brighter areas of the plage and downflows at the weaker areas. Besides, we propose a new method to determine the chirality of an erupting filament based on the skewness of the conjugate filament drainage sites. It suggests that the right-skewed drainage corresponds to sinistral chirality, whereas the left-skewed drainage corresponds to dextral chirality.

High-resolution H$\alpha$ observations indicate that filaments consist of an assembly of thin threads. In quiescent filaments, the threads are generally short, whereas in active region filaments, the threads are generally long. In order to explain these observational features, we performed one-dimensional radiative hydrodynamic simulations of filament formation along a dipped magnetic flux tube in the framework of the chromospheric evaporation-coronal condensation model. The geometry of a dipped magnetic flux tube is characterized by three parameters, i.e., the depth ($D$), the half-width ($w$), and the altitude ($h$) of the magnetic dip. The parameter survey in the numerical simulations shows that allowing the filament thread to grow in 5 days, the maximum length ($L_{th}$) of the filament thread increases linearly with $w$, and decreases linearly with $D$ and $h$. The dependence is fitted into a linear function $L_{th}=0.84w-0.88D-2.78h+17.31$ Mm. Such a relation can qualitatively explain why quiescent filaments have shorter threads and active region filaments have longer threads.

A new multi-wavelength solar telescope, Optical and Near-infrared Solar Eruption Tracer (ONSET) of Nanjing University, was constructed, being fabricated by Nanjing Institute of Astronomical Optics & Technology and run in cooperation with Yunnan Astronomical Observatory. ONSET is able to observe the Sun in three wavelength windows: He {\small I} 10830 \AA, H$\alpha$, and white-light at 3600 \AA and 4250 \AA, which are selected in order to obtain the dynamics in the corona, chromosphere, and the photosphere simultaneously. Full-disk or partial-disk solar images with a field of 10 arcmin at three wavelengths can be obtained nearly simultaneously. It is designed to trace solar eruptions with high spatial and temporal resolutions. This telescope was installed at a new solar observing site near Fuxian Lake in Yunnan Province, southwest China. The site is located at E102N24, with an altitude of 1722 m. The seeing is stable and of high quality. We give a brief description of the scientific objectives and the basic structure of the ONSET. Some preliminary results are also presented.

An X3.4 solar flare and a fast halo coronal mass ejection (CME) occurred on 2006 December 13, accompanied by a high flux of energetic particles recorded both in near-Earth space and at ground level. Our purpose is to provide evidence of flare acceleration in a major solar energetic particle (SEP) event. We first present observations from ACE/EPAM, GOES, and the Apatity neutron monitor. It is found that the initial particle release time coincides with the flare emission and that the spectrum becomes softer and the anisotropy becomes weaker during particle injection, indicating that the acceleration source changes from a confined coronal site to a widespread interplanetary CME-driven shock. We then describe a comprehensive study of the associated flare active region. By use of imaging data from HINODE/SOT and SOHO/MDI magnetogram, we infer the flare magnetic reconnection rate in the form of the magnetic flux change rate. This correlates in time with the microwave emission, indicating a physical link between the flare magnetic reconnection and the acceleration of nonthermal particles. Combining radio spectrograph data from Huairou/NOAC, Culgoora/IPS, Learmonth/RSTN, and WAVES/WIND leads to a continuous and longlasting radio burst extending from a few GHz down to several kHz. Based on the photospheric vector magnetogram from Huairou/NOAC and the nonlinear force free field (NFFF) reconstruction method, we derive the 3D magnetic field configuration shortly after the eruption. Furthermore, we also compute coronal field lines extending to a few solar radii using a potential-field source-surface (PFSS) model. Both the so-called type III-l burst and the magnetic field configuration suggest that open-field lines extend from the flare active region into interplanetary space, allowing the accelerated and charged particles escape.

We developed a method to automatically detect and trace solar filaments in H\alpha\ full-disk images. The program is able not only to recognize filaments and determine their properties, such as the position, the area, the spine, and other relevant parameters, but also to trace the daily evolution of the filaments. The program consists of three steps: First, preprocessing is applied to correct the original images; Second, the Canny edge-detection method is used to detect filaments; Third, filament properties are recognized through the morphological operators. To test the algorithm, we applied it to the observations from the Mauna Loa Solar Observatory (MLSO), and the program is demonstrated to be robust and efficient. H\alpha\ images obtained by MLSO from 1998 to 2009 are analyzed, and a butterfly diagram of filaments is obtained. It shows that the latitudinal migration of solar filaments has three trends in the Solar Cycle 23: The drift velocity was fast from 1998 to the solar maximum; After the solar maximum, it became relatively slow. After 2006, the migration became divergent, signifying the solar minimum. About 60% filaments with the latitudes larger than $50^{\circ}$ migrate towards the polar regions with relatively high velocities, and the latitudinal migrating speeds in the northern and the southern hemispheres do not differ significantly in the Solar Cycle 23.

We attempt to understand the white-light flare (WLF) that was observed on 2012 March 9 with a newly constructed multi-wavelength solar telescope called the Optical and Near-infrared Solar Eruption Tracer (ONSET). We analyzed WLF observations in radio, H-alpha, white-light, ultraviolet, and X-ray bands. We also studied the magnetic configuration of the flare via the nonlinear force-free field (NLFFF) extrapolation and the vector magnetic field observed by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). Continuum emission enhancement clearly appeared at the 3600 angstrom and 4250 angstrom bands, with peak contrasts of 25% and 12%, respectively. The continuum emission enhancement closely coincided with the impulsive increase in the hard X-ray emission and a microwave type III burst at 03:40 UT. We find that the WLF appeared at one end of either the sheared or twisted field lines or both. There was also a long-lasting phase in the H-alpha and soft X-ray bands after the white-light emission peak. In particular, a second, yet stronger, peak appeared at 03:56 UT in the microwave band. This event shows clear evidence that the white-light emission was caused by energetic particles bombarding the lower solar atmosphere. A two-step magnetic reconnection scenario is proposed to explain the entire process of flare evolution, i.e., the first-step magnetic reconnection between the field lines that are highly sheared or twisted or both, and the second-step one in the current sheet, which is stretched by the erupting flux rope. The WLF is supposed to be triggered in the first-step magnetic reconnection at a relatively low altitude.

In this paper we investigate the perturbation theory of the asymptotically safe inflation and we find that all modes of gravitational waves perturbation become ghosts in order to achieve a large enough number of e-folds. Formally we can calculate the power spectrum of gravitational waves perturbation, but we find that it is negative. It indicates that there is serious trouble with the asymptotically safe inflation.

A new code, named MAP, is written in Fortran language for magnetohydrodynamics (MHD) calculation with the adaptive mesh refinement (AMR) and Message Passing Interface (MPI) parallelization. There are several optional numerical schemes for computing the MHD part, namely, modified Mac Cormack Scheme (MMC), Lax-Friedrichs scheme (LF) and weighted essentially non-oscillatory (WENO) scheme. All of them are second order, two-step, component-wise schemes for hyperbolic conservative equations. The total variation diminishing (TVD) limiters and approximate Riemann solvers are also equipped. A high resolution can be achieved by the hierarchical block-structured AMR mesh. We use the extended generalized Lagrange multiplier (EGLM) MHD equations to reduce the non-divergence free error produced by the scheme in the magnetic induction equation. The numerical algorithms for the non-ideal terms, e.g., the resistivity and the thermal conduction, are also equipped in the MAP code. The details of the AMR and MPI algorithms are described in the paper.

Microflares are small activities in solar low atmosphere, some are in the low corona, and others in the chromosphere. Observations show that some of the microflares are triggered by magnetic reconnection between emerging flux and a pre-existing background magnetic field. We perform 2.5D compressible resistive MHD simulations of magnetic reconnection with gravity considered. The background magnetic field is a canopy-type configuration which is rooted at the boundary of the solar supergranule. By changing the bottom boundary conditions in the simulation, new magnetic flux emerges up at the center of the supergranule and reconnects with the canopy-type magnetic field. We successfully simulate the coronal and chromospheric microflares, whose current sheets are located at the corona and the chromosphere, respectively. The microflare of coronal origin has a bigger size and a higher temperature enhancement than that of chromospheric origin. In the microflares of coronal origin, we also found a hot jet ($\sim$$1.8 \times 10^6$ K), which is probably related to the observational EUV/SXR jets, and a cold jet ($\sim$$10^4$ K), which is similar to the observational H$\alpha$/Ca surges, whereas there is only an H$\alpha$/Ca bright point in the microflares of chromospheric origin. The study of parameter dependence shows that the size and strength of the emerging magnetic flux are the key parameters which determine the height of the reconnection location, and further determine the different observational features of the microflares.

Coronal "EIT waves" appear as EUV bright fronts propagating across a significant part of the solar disk. The intriguing phenomenon provoked continuing debates on their nature and their relation with coronal mass ejections (CMEs). In this paper, we first summarize all the observational features of "EIT waves", which should be accounted for by any successful model. The theoretical models constructed during the past 10 years are then reviewed. Finally, the implication of the "EIT wave" research to the understanding of CMEs is discussed. The necessity is pointed out to revisit the nature of CME frontal loop.

We performed three dimensional resistive magnetohydrodynamic simulations to study the magnetic reconnection using an initially shearing magnetic field configuration (force free field with a current sheet in the middle of the computational box). It is shown that there are two types of reconnection jets: the ordinary reconnection jets and fan-shaped jets, which are formed along the guide magnetic field. The fan-shaped jets are much different from the ordinary reconnection jets which are ejected by magnetic tension force. There are two driving forces for accelerating the fan-shaped jets. The one is the Lorentz force which dominates the motion of fluid elements at first and then the gas pressure gradient force accelerates the fluid elements in the later stage. The dependence on magnetic reconnection angle and resistivity value has also been studied. The formation and evolution of these jets provide a new understanding of dynamic magnetohydrodynamic jets.

Magnetic reconnection is a fundamental process in space and astrophysical plasmas in which oppositely directed magnetic fields changes its connectivity and eventually converts its energy into kinetic and thermal energy of the plasma. Recently, ubiquitous jets (for example, chromospheric anemone jets, penumbral microjets, umbral light bridge jets) have been observed by Solar Optical Telescope on board the satellite Hinode. These tiny and frequently occurring jets are considered to be a possible evidence of small-scale ubiquitous reconnection in the solar atmosphere. However, the details of three dimensional magnetic configuration are still not very clear. Here we propose a new model based on three dimensional simulations of magnetic reconnection using a typical current sheet magnetic configuration with a strong guide field. The most interesting feature is that the jets produced by the reconnection eventually move along the guide field lines. This model provides a fresh understanding of newly discovered ubiquitous jets and moreover a new observational basis for the theory of astrophysical magnetic reconnection.

With gravity, ionization, and radiation being considered, we perform 2.5D compressible resistive MHD simulations of chromospheric magnetic reconnection using the CIP-MOCCT scheme. The temperature distribution of the quiet-Sun atmospheric model VALC and the helium abundance (10%) are adopted. Our 2.5D MHD simulation reproduces qualitatively the temperature enhancement observed in chromospheric microflares. The temperature enhancement $\Delta T$ is demonstrated to be sensitive to the background magnetic field, whereas the total evolution time $\Delta t$ is sensitive to the magnitude of the anomalous resistivity. Moveover, we found a scaling law, which is described as $\Delta T/\Delta t \sim {n_H}^{-1.5} B^{2.1} {\eta_0}^{0.88}$. Our results also indicate that the velocity of the upward jet is much greater than that of the downward jet and the X-point may move up or down.

Halo coronal mass ejections (CMEs) were found to be significantly faster than normal CMEs, which was a long-standing puzzle. In order to solve the puzzle, we first investigate the observed properties of 31 limb CMEs that display clearly loop-shaped frontal loops. The observational results show a strong tendency that slower CMEs are weaker in the white-light intensity. Then, we perform a Monte Carlo simulation of 20 000 artificial limb CMEs that have average velocity of $\sim$523 km s$^{-1}$. The Thomson scattering of these events is calculated when they are assumed to be observed as limb and halo events, respectively. It is found that the white-light intensity of many slow CMEs becomes remarkably reduced as they turn from being viewed as a limb event to as a halo event. When the intensity is below the background solar wind fluctuation, it is assumed that they would be missed by coronagraphs. The average velocity of "detectable" halo CMEs is $\sim$922 km s$^{-1}$, very close to the observed value. It also indicates that wider events are more likely to be recorded. The results soundly suggest that the higher average velocity of halo CMEs is due to that a majority of slow events and a part of narrow fast events carrying less material are so faint that they are blended with the solar wind fluctuations, and therefore are not observed.

The direct current (DC) electric field near the reconnection region has been proposed as an effective mechanism to accelerate protons and electrons in solar flares. A power-law energy spectrum was generally claimed in the simulations of electron acceleration by the reconnection electric field. However in most of the literature, the electric and magnetic fields were chosen independently. In this paper, we perform test-particle simulations of electron acceleration in a reconnecting magnetic field, where both the electric and magnetic fields are adopted from numerical simulations of the MHD equations. It is found that the accelerated electrons present a truncated power-law energy spectrum with an exponential tail at high energies, which is analogous to the case of diffusive shock acceleration. The influences of reconnection parameters on the spectral feature are also investigated, such as the longitudinal and transverse components of the magnetic field and the size of the current sheet. It is suggested that the DC electric field alone might not be able to reproduce the observed single or double power-law distributions.

In large gradual solar energetic particle (SEP) events, especially the ground-level enhancement (GLE) events, where and how energetic particles are accelerated is still a problem. By using imaging data from TRACE, Yohkoh/HXT, SOHO/MDI and SOHO/EIT, along with the data from the GOES, Apatity NM, and SOHO/LASCO CME catalog, the evolution of the X5.7 two-ribbon flare and the associated SEP event on 14 July 2000 are studied. It is found that the magnetic reconnection in this event consists of two parts, and the induced electric field Erec is temporally correlated with the evolution of hard X-ray and gamma-ray emission. In particular, the first hard X-ray and gamma-ray emission peak occurred at 10:22 UT, corresponding to the magnetic reconnection in the western part of the flare ribbons and the maximum Erec of 9.5 V/cm; the second emission peak at 10:27 UT, corresponding to the eastern part and the maximum Erec of 13.0 V/cm. We also analyze the SEP injection profiles as functions of time and CME-height, and find two-component injection which may result from different acceleration mechanisms. A reasonable conclusion is that reconnection electric field makes a crucial contribution to the acceleration of relativistic particles and to the impulsive component of the large gradual SEP event, while CME-driven shocks play a dominant role in the gradual component.

An X17.2 solar flare occurred on 2003 October 28, accompanied by multi-wavelength emissions and a high flux of relativistic particles observed at 1AU. We present the analytic results of the TRACE, SOHO, RHESSI, ACE, GOES, hard X-ray (INTEGRAL satellite), radio (Onderejov radio telescope), and neutron monitor data. It is found that the inferred magnetic reconnection electric field correlates well with the hard X-ray, gamma-ray, and neutron emission at the Sun. Thus the flare's magnetic reconnection probably makes a crucial contribution to the prompt relativistic particles, which could be detected at 1 AU. Since the neutrons were emitted a few minutes before the injection of protons and electrons, we propose a magnetic-field evolution configuration to explain this delay. We do not exclude the effect of CME-driven shock, which probably plays an important role in the delayed gradual phase of solar energetic particles.

There is observational evidence showing that stellar and solar flares occur in a similar circumstance, although the former are usually much more energetic. It is expected that the bombardment by high energy electrons is one of the chief heating processes of the flaring atmosphere. In this paper we study how a precipitating electron beam can influence the line profiles of Ly_alpha, H_alpha, Ca II K and Ca II 8542. We use a model atmosphere of a dMe star and make non-LTE computations taking into account the non-thermal collisional rates due to the electron beam. The results show that the four lines can be enhanced to different extents. The relative enhancement increases with increasing formation height of the lines. Varying the energy flux of the electron beam has different effects on the four lines. The wings of Ly_alpha and H_alpha become increasingly broad with the beam flux; change of the Ca II K and Ca II 8542 lines, however, is most significant in the line centre. Varying the electron energy (i.e., the low-energy cut-off for a power law beam) has a great influence on the Ly_alpha line, but little on the H_alpha and Ca II lines. An electron beam of higher energy precipitates deeper, thus producing less enhancement of the Ly_alpha line. The Ly_alpha/H_alpha flux ratio is thus sensitive to the electron energy.

We explore the possibility that the continuum emission in stellar flares is powered by nonthermal electrons accelerated during the flares. We compute the continuum spectra from an atmospheric model for a dMe star, AD Leo, at its quiescent state, when considering the nonthermal excitation and ionisation effects by precipitating electron beams. The results show that if the electron beam has an energy flux large enough, the U band brightening and, in particular, the U-B colour are roughly comparable with observed values for a typical large flare. Moreover, for electron beams with a moderate energy flux, a decrease of the emission at the Paschen continuum appears. This can explain at least partly the continuum dimming observed in some stellar flares. Adopting an atmospheric model for the flaring state can further raise the continuum flux but it yields a spectral colour incomparable with observations. This implies that the nonthermal effects may play the chief role in powering the continuum emission in some stellar flares.