The He I 10830 \AA\ triplet is a very informative indicator of chromospheric activities as the helium is the second most abundant element in the solar atmosphere. Taking advantage of the high resolution of the 1.6 m Goode Solar Telescope (GST) at Big Bear Solar Observatory (BBSO), previous observations have shown clear evidence of the enhanced absorption, instead of typically-observed emission, for two M-class flares. In this study, we analyze the evolution of the He I 10830 10830 \AA\ emission in numerical models and compare it with observations. The models represent the RADYN simulation results obtained from the F-CHROMA database. We consider the models with the injected electron spectra parameters close to observational estimates for the 2013-August-17 flare event ($\delta=8$, $E_c = \{15,20\}$ keV, $F=\{1\times 10^{11}, 3\times{}10^{11}\}$ erg cm$^{-2}$) in detail, as well as other available models. The modeling results agree well with observations, in the sense of both the maximum intensity decrease (-17.1%, compared to the observed value of -13.7%) and the trend of temporal variation (initial absorption phase followed by the emission). All models demonstrate the increased number densities and decreased ratio of the upper and lower level populations of He I 10830 10830 \AA\ transition in the initial phase, which enhances the opacity and forms an absorption feature. Models suggest that the temperatures and free electron densities at heights of 1.3-1.5 Mm should be larger than $\sim 10^4$ K and $6\times 10^{11}$ cm$^{-3}$ thresholds for the line to start being in emission.

We propose a new machine learning approach to Stokes inversion based on a convolutional neural network (CNN) and the Milne-Eddington (ME) method. The Stokes measurements used in this study were taken by the Near InfraRed Imaging Spectropolarimeter (NIRIS) on the 1.6 m Goode Solar Telescope (GST) at the Big Bear Solar Observatory. By learning the latent patterns in the training data prepared by the physics-based ME tool, the proposed CNN method is able to infer vector magnetic fields from the Stokes profiles of GST/NIRIS. Experimental results show that our CNN method produces smoother and cleaner magnetic maps than the widely used ME method. Furthermore, the CNN method is 4~6 times faster than the ME method, and is able to produce vector magnetic fields in near real-time, which is essential to space weather forecasting. Specifically, it takes ~50 seconds for the CNN method to process an image of 720 x 720 pixels comprising Stokes profiles of GST/NIRIS. Finally, the CNN-inferred results are highly correlated to the ME-calculated results and are closer to the ME's results with the Pearson product-moment correlation coefficient (PPMCC) being closer to 1 on average than those from other machine learning algorithms such as multiple support vector regression and multilayer perceptrons (MLP). In particular, the CNN method outperforms the current best machine learning method (MLP) by 2.6% on average in PPMCC according to our experimental study. Thus, the proposed physics-assisted deep learning-based CNN tool can be considered as an alternative, efficient method for Stokes inversion for high resolution polarimetric observations obtained by GST/NIRIS.

We have studied numerically the shadows of a non-Kerr rotating compact object with quadrupole mass moment, which belongs to Manko-Novikov family. The non-integrable photon motion caused by quadrupole mass moment affects sharply the shadow of the compact object. As the deviation parameter related to quadrupole mass moment is negative, the shadow of compact object is prolate and there are two disconnected main shadows with eyebrows located symmetrically on both sides of the equatorial plane. As the deviation parameter is positive, the shadow becomes oblate and the main shadow is joined together in the equatorial plane. Moreover, in this positive cases, there is a disorder region in the left of shadow which increases with the quadrupole-deviation parameter. Interestingly, we also find that Einstein ring is broken as the deviation from Kerr metric is larger than a certain critical value. This critical value decreases with the rotation parameter of black hole. Especially, the observer on the direction of rotation axis will find some concentric bright rings in the black disc. Finally, supposing that the gravitational field of the supermassive central object of the galaxy described by this metric, we estimated the numerical values of the observables for the black hole shadow.

A recent laboratory experiment of ideal magnetohydrodynamic (MHD) instabilities reveals four distinct eruption regimes readily distinguished by the torus instability (TI) and helical kink instability (KI) parameters \citep{Myers2015}. To establish its observational counterpart, we collect 38 solar flares (stronger than GOES class M5 in general) that took place within 45$^{\circ}$ of disk center during 2011$-$2017, 26 of which are associated with a halo or partial halo coronal mass ejection (CME) (i.e., ejective events), while the others are CMEless (i.e., confined events). This is a complete sample of solar events satisfying our selection criteria detailed in the paper. For each event, we calculate decay index $n$ of the potential strapping field above the magnetic flux rope (MFR) in and around the flaring magnetic polarity inversion line (a TI parameter), and the unsigned twist number $T_w$ of the non-linear force-free (NLFF) field lines forming the same MFR (a KI parameter). We then construct a $n-T_w$ diagram to investigate how the eruptiveness depends on these parameters. We find: (1) $T_w$ appears to play little role in discriminating between confined and ejective events; (2) the events with $n\gtrsim0.8$ are all ejective and all confined events have $n\lesssim0.8$. However, $n\gtrsim0.8$ is not a necessary condition for eruption, because some events with $n\lesssim0.8$ also erupted. In addition, we investigate the MFR's geometrical parameters, apex height and distance between footpoints, as a possible factor for the eruptiveness. We briefly discuss the difference of the present result for solar eruptions with that of the laboratory result in terms of the role played by magnetic reconnection.

Polar Crown Filaments (PCFs) form above the magnetic polarity inversion line, which separates the unipolar polar fields and the nearest dispersed fields from trailing part of active regions with opposite polarity. The statistical properties of PCFs are correlated with the solar cycle. Therefore, study of PCFs plays an important role in understanding the variation of solar cycle, especially the prolonged cycle 23 and the current "abnormal" solar cycle 24. In this study, we investigate PCFs using full disk H{\alpha} data from 1973 to early 2018, recorded by Kanzelhohe Solar Observatory (KSO) and Big Bear Solar Observatory (BBSO), in digital form from 1997 to 2018 and in 35 mm film (digitized) from 1973 to 1996. PCFs are identified manually because their segmented shape and close-to-limb location were not handled well by automatical detections in several previous studies. Our results show that the PCFs start to move poleward at the beginning of each solar cycle. When the PCFs approach to the maximum latitude, the polar field strength reduces to zero followed by a reversal. The migration rates are about 0.4 to 0.7 degree per Carrington rotation, with clear N-S asymmetric pattern. In cycles 21 and 23, the PCFs in the northern hemisphere migrate faster than those in the southern hemisphere. However, in the "abnormal" cycle 24, the southern PCFs migrate faster, which is consistent with other observations of magnetic fields and radio emission. In addition, there are more days in cycle 23 and 24 without PCFs than in the previous cycles.

As one of the most violent eruptions on the Sun, flares are believed to be powered by magnetic reconnection. The fundamental physics involving the release, transfer and deposition of energy have been studied extensively. Taking advantage of the unprecedented resolution provided by the 1.6-m Goode Solar Telescope, here we show a sudden rotation of vector magnetic fields, about 12$^{\circ}$-20$^{\circ}$ counterclockwise, associated with a flare. Unlike the permanent changes reported previously, the azimuth-angle change is transient and co-spatial/temporal with H$\alpha$ emission. The measured azimuth angle becomes closer to that in potential fields suggesting untwist of flare loops. The magnetograms were obtained in the near infrared at 1.56~$\mu$m, which is minimally affected by flare emission and no intensity profile change was detected. We believe that these transient changes are real and discuss the possible explanations in which the high energy electron beams or $Alfv\acute{e}n$ waves play a crucial role.

The well-established notion of a "common population" of the accelerated electrons simultaneously producing the hard X-ray (HXR) and the microwave (MW) emission during the flare impulsive phase has been challenged by some studies reporting the discrepancies between the HXR-inferred and the MW-inferred electron energy spectra. The traditional methods of their spectral inversion have some problems that can be mainly attributed to the unrealistic and the oversimplified treatment of the flare emission. To properly address this problem, we use a Non-linear Force Free Field (NLFFF) model extrapolated from an observed photospheric magnetogram as input to the three-dimensional, multi-wavelength modeling platform GX Simulator, and create a unified electron population model that can simultaneously reproduce the observed HXR and MW observations. We model the end of the impulsive phase of the 2015-06-22 M6.5 flare, and constrain the modeled electron spatial and energy parameters using observations made by the highest-resolving instruments currently available in two wavelengths, the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) for HXR and the Expanded Owens Valley Solar Array (EOVSA) for MW. Our results suggest that the HXR-emitting electron population model fits the standard flare model with a broken power-law spectrum (E_break ~ 200 keV) that simultaneously produces the HXR footpoint emission and the MW high frequency emission. The model also includes an "HXR invisible" population of nonthermal electrons that are trapped in a large volume of magnetic field above the HXR-emitting loops, which is observable by its gyrosynchrotron (GS) radiation emitting mainly in MW low frequency range.

We analyze the high-resolution observations from the 1.6m telescope at Big Bear Solar Observatory that cover an active region filament. Counter-streaming motions are clearly observed in the filament. The northern end of the counter-streaming motions extends to a light bridge, forming a spectacular circulation pattern around a sunspot, with clockwise motion in the blue wing and counterclockwise motion in the red wing as observed in Halpha off-bands. The apparent speed of the flow is around 10 to 60 km/s in the filament, decreasing to 5 to 20 km/s in the light bridge. The most intriguing results are the magnetic structure and the counter-streaming motions in the light bridge. Similar to those in the filament, magnetic fields show a dominant transverse component in the light bridge. However, the filament is located between opposite magnetic polarities, while the light bridge is between strong fields of the same polarity. We analyze the power of oscillations with the image sequences of constructed Dopplergrams, and find that the filament's counter-streaming motion is due to physical mass motion along fibrils, while the light bridge's counter-streaming motion is due to oscillation in the direction along the line-of-sight. The oscillation power peaks around 4 minutes. However, the section of the light bridge next to the filament also contains a component of the extension of the filament in combination with the oscillation, indicating that some strands of the filament are extended to and rooted in that part of the light bridge.

In the solar corona, magnetic helicity slowly and continuously accumulates in response to plasma flows tangential to the photosphere and magnetic flux emergence through it. Analyzing this transfer of magnetic helicity is key for identifying its role in the dynamics of active regions (ARs). The connectivity-based helicity flux density method was recently developed for studying the 2D and 3D transfer of magnetic helicity in ARs. The method takes into account the 3D nature of magnetic helicity by explicitly using knowledge of the magnetic field connectivity, which allows it to faithfully track the photospheric flux of magnetic helicity. Because the magnetic field is not measured in the solar corona, modeled 3D solutions obtained from force-free magnetic field extrapolations must be used to derive the magnetic connectivity. Different extrapolation methods can lead to markedly different 3D magnetic field connectivities, thus questioning the reliability of the connectivity-based approach in observational applications. We address these concerns by applying this method to the isolated and internally complex AR 11158 with different magnetic field extrapolation models. We show that the connectivity-based calculations are robust to different extrapolation methods, in particular with regards to identifying regions of opposite magnetic helicity flux. We conclude that the connectivity-based approach can be reliably used in observational analyses and is a promising tool for studying the transfer of magnetic helicity in ARs and relate it to their flaring activity.

It is now a common practice to use local magnetic parameters such as magnetic decay index for explaining solar eruptions from active regions, but there can be an alternative view that the global properties of the source region should be counted as a more important factor. We discuss this issue based on {\it Solar Dynamics Observatory} (SDO) observations of the three successive eruptions within 1.5 hours from the NOAA active region 11444 and the magnetic parameters calculated using the nonlinear force-free field (NLFFF) model. Two violent eruptions occurred in the regions with relatively high magnetic twist number (0.5--1.5) and high decay index (0.9--1.1) at the nominal height of the filament (12$''$) and otherwise a mild eruption occurred, which supports the local parameter paradigm. Our main point is that the time sequence of the eruptions did not go with these parameters. It is argued that an additional factor, in the form of stabilizing force, should operate to determine the onset of the first eruption and temporal behaviors of subsequent eruptions. As supporting evidence, we report that the heating and fast plasma flow continuing for a timescale of an hour was the direct cause for the first eruption, and that the unidirectional propagation of the disturbance determined the timing of subsequent eruptions. Both of these factors are associated with the overall magnetic structure rather than local magnetic properties of the active region.

How eruption can recur from a confined magnetic structure is discussed based on the {\it Solar Dynamics Observatory} (SDO) observations of the NOAA active region 11444, which produced three eruptions within 1.5 hours on March 27, 2012. The active region had the positive polarity magnetic fields in the center surrounded by the negative polarity fields around. Since such a distribution of magnetic polarity tends to form a dome-like magnetic fan structure confined over the active region, the multiple eruptions was puzzling. Our investigation reveals that this event exhibits several properties distinct from other eruptions associated with magnetic fan structures: (i) a long filament encircling the active region was present before the eruptions; (ii) expansion of the open-closed boundary of the field lines after each eruption suggestive of the growing fan-dome structure, and (iii) the ribbons inside the closed magnetic polarity inversion line evolving in response to the expanding open-closed boundary. It thus appears that in spite of multiple eruptions the fan-dome structure remained undamaged, and the closing back field lines after each eruption rather reinforced the fan-dome structure. We argue that the multiple eruptions could occur in this active region in spite of its confined magnetic structure because the filament encircling the active region was adequate for slipping through the magnetic separatrix to minimize the damage to its overlying fan-dome structure. The result of this study provides a new insight into the productivity of eruptions from a confined magnetic structure.

We report the intriguing large-scale dynamic phenomena associated with the M6.5 flare~(SOL2015-06-22T18:23) in NOAA active region 12371, observed by RHESSI, Fermi, and the Atmospheric Image Assembly (AIA) and Magnetic Imager (HMI) on the Solar Dynamic Observatory (SDO). The most interesting feature of this event is a third ribbon (R3) arising in the decay phase, propagating along a dimming channel (seen in EUV passbands) towards a neighboring sunspot. The propagation of R3 occurs in the presence of hard X-ray footpoint emission, and is broadly visible at temperatures from 0.6~MK to over 10~MK through the Differential Emission Measure (DEM) analysis. The coronal loops then undergo an apparent slipping motion following the same path of R3, after a $\sim$ 80~min delay. To understand the underlying physics, we investigate the magnetic configuration and the thermal structure of the flaring region. Our results are in favor of a slipping-type reconnection followed by the thermodynamic evolution of coronal loops. In comparison with those previously reported slipping reconnection events, this one proceeds across a particularly long distance ($\sim$60 Mm) over a long period of time ($\sim$50 min), and shows two clearly distinguished phases: the propagation of the footpoint brightening driven by nonthermal particle injection and the apparent slippage of loops governed by plasma heating and subsequent cooling.

Studying the dynamics of filaments at pre-eruption phase can shed light on the precursor of eruptive events. Such studies in high-resolution (in the order of 0.1") are highly desirable yet very rare so far. In this work, we present a detailed observation of a pre-eruption evolution of a filament obtained by the 1.6 m New Solar Telescope (NST) at Big Bear Solar Observatory (BBSO). One end of the filament is anchored at the sunspot in NOAA active region (AR) 11515, which is well observed by NST H$\alpha$ off-bands four hours before till one hour after the filament eruption. A M1.6 flare is associated with the eruption. We observed persistent downflowing materials along the H$\alpha$ multi-threaded component of the loop towards the AR end during the pre-eruption phase. We traced the trajectories of plasma blobs along the H$\alpha$ threads and obtained the plane-of-sky velocity of 45 km s$^{-1}$ on average. We further estimated the real velocities of the downflows and the altitude of the filament by matching the observed H$\alpha$ threads with magnetic field lines extrapolated from a nonlinear force-free field (NLFFF) model. Observation of chromospheric brightenings (BZs) at the footpoints of the falling plasma blobs is also presented in the paper. The lower limit of the kinetic energy per second of the downflows through the BZs is found to be $\sim$ 10$^{21}$ erg. Larger FOV observations from BBSO full disk H$\alpha$ images show that the AR end of the filament started ascending four hours before the flare. We attribute the observed downflows at the AR end of the filament to the draining effect of the filament rising prior to its eruption. During the slow-rise phase, the downflows continuously drained away $\sim$ 10$^{15}$g mass from the filament over a few hours, which is believed to be essential for the instability at last, and could be an important precursor of eruptive events.

Sunspots are concentrations of magnetic field visible on the solar surface (photosphere). It was considered implausible that solar flares, as resulted from magnetic reconnection in the tenuous corona, would cause a direct perturbation of the dense photosphere involving bulk motion. Here we report the sudden flare-induced rotation of a sunspot using the unprecedented spatiotemporal resolution of the 1.6 m New Solar Telescope, supplemented by magnetic data from the Solar Dynamics Observatory. It is clearly observed that the rotation is non-uniform over the sunspot: as the flare ribbon sweeps across, its different portions accelerate (up to 50 deg per hr) at different times corresponding to peaks of flare hard X-ray emission. The rotation may be driven by the surface Lorentz-force change due to the back reaction of coronal magnetic restructuring and is accompanied by a downward Poynting flux. These results have direct consequences for our understanding of energy and momentum transportation in the flare-related phenomena.

Solar flares are generally believed to be powered by free magnetic energy stored in the corona, but the build up of coronal energy alone may be insufficient for the imminent flare occurrence. The flare onset mechanism is a critical but less understood problem, insights into which could be gained from small-scale energy releases known as precursors, which are observed as small pre-flare brightenings in various wavelengths, and also from certain small-scale magnetic configurations such as the opposite polarity fluxes, where magnetic orientation of small bipoles is opposite to that of the ambient main polarities. However, high-resolution observations of flare precursors together with the associated photospheric magnetic field dynamics are lacking. Here we study precursors of a flare using unprecedented spatiotemporal resolution of the 1.6 m New Solar Telescope, complemented by novel microwave data. Two episodes of precursor brightenings are initiated at a small-scale magnetic channel (a form of opposite polarity fluxes) with multiple polarity inversions and enhanced magnetic fluxes and currents, lying near the footpoints of sheared magnetic loops. The low-atmospheric origin of these precursor emissions is corroborated by microwave spectra. We propose that the emerging magnetic channel field interacts with the sheared arcades to cause precursor brightenings at the main flare core region. These high-resolution results provide evidence of low-atmospheric small-scale energy release and possible relationship to the onset of the main flare.

In this study we synthesize the results of four previous studies on the global energetics of solar flares and associated coronal mass ejections (CMEs), which include magnetic, thermal, nonthermal, and CME energies in 399 solar M and X-class flare events observed during the first 3.5 years of the Solar Dynamics Observatory (SDO) mission. Our findings are: (1) The sum of the mean nonthermal energy of flare-accelerated particles ($E_{\mathrm{nt}}$), the energy of direct heating ($E_{\mathrm{dir}}$), and the energy in coronal mass ejections ($E_{\mathrm{CME}}$), which are the primary energy dissipation processes in a flare, is found to have a ratio of $(E_{\mathrm{nt}}+E_{\mathrm{dir}}+ E_{\mathrm{CME}})/E_{\mathrm{mag}} = 0.87 \pm 0.18$, compared with the dissipated magnetic free energy $E_{\mathrm{mag}}$, which confirms energy closure within the measurement uncertainties and corroborates the magnetic origin of flares and CMEs; (2) The energy partition of the dissipated magnetic free energy is: $0.51\pm0.17$ in nonthermal energy of $\ge 6$ keV electrons, $0.17\pm0.17$ in nonthermal $\ge 1$ MeV ions, $0.07\pm0.14$ in CMEs, and $0.07\pm0.17$ in direct heating; (3) The thermal energy is almost always less than the nonthermal energy, which is consistent with the thick-target model; (4) The bolometric luminosity in white-light flares is comparable with the thermal energy in soft X-rays (SXR); (5) Solar Energetic Particle (SEP) events carry a fraction $\approx 0.03$ of the CME energy, which is consistent with CME-driven shock acceleration; and (6) The warm-target model predicts a lower limit of the low-energy cutoff at $e_c \approx 6$ keV, based on the mean differential emission measure (DEM) peak temperature of $T_e=8.6$ MK during flares. This work represents the first statistical study that establishes energy closure in solar flare/CME events.

Solar flares signify the sudden release of magnetic energy and are sources of so called space weather. The fine structures (below 500 km) of flares are rarely observed and are accessible to only a few instruments world-wide. Here we present observation of a solar flare using exceptionally high resolution images from the 1.6~m New Solar Telescope (NST) equipped with high order adaptive optics at Big Bear Solar Observatory (BBSO). The observation reveals the process of the flare in unprecedented detail, including the flare ribbon propagating across the sunspots, coronal rain (made of condensing plasma) streaming down along the post-flare loops, and the chromosphere's response to the impact of coronal rain, showing fine-scale brightenings at the footpoints of the falling plasma. Taking advantage of the resolving power of the NST, we measure the cross-sectional widths of flare ribbons, post-flare loops and footpoint brighenings, which generally lie in the range of 80-200 km, well below the resolution of most current instruments used for flare studies. Confining the scale of such fine structure provides an essential piece of information in modeling the energy transport mechanism of flares, which is an important issue in solar and plasma physics.

Solar flares are sudden flashes of brightness on the Sun and are often associated with coronal mass ejections and solar energetic particles which have adverse effects in the near Earth environment. By definition, flares are usually referred to bright features resulting from excess emission. Using the newly commissioned 1.6~m New Solar Telescope at Big Bear Solar Observatory, here we show a striking "negative" flare with a narrow, but unambiguous "dark" moving front observed in He I 10830 \AA, which is as narrow as 340 km and is associated with distinct spectral characteristics in H-alpha and Mg II lines. Theoretically, such negative contrast in He I 10830 \AA\ can be produced under special circumstances, by nonthermal-electron collisions, or photoionization followed by recombination. Our discovery, made possible due to unprecedented spatial resolution, confirms the presence of the required plasma conditions and provides unique information in understanding the energy release and radiative transfer in astronomical objects.

We present a set of tools for detecting small-scale solar magnetic cancellations and the disk counterpart of type II spicules (the so-called Rapid Blueshifted Excursions (RBEs)), using line-of-sight photospheric magnetograms and chromospheric spectroscopic observations, respectively. For tracking magnetic cancellation, we improve the Southwest Automatic Magnetic Identification Suite (SWAMIS) so that it is able to detect certain obscure cancellations that can be easily missed. For detecting RBEs, we use a normalized reference profile to reduce false-positive detections caused by the non-uniform background and seeing condition. Similar to the magnetic feature tracking in SWAMIS, we apply a dual-threshold method to enhance the accuracy of RBE detection. These tools are employed to analyze our coordinated observations using the Interferometric BIdimensional Spectrometer at Dunn Solar Telescope (DST) of the National Solar Observatory (NSO) and Hinode. We present the statistical properties of magnetic cancellations and RBEs, and explore their correlation using this data set.

Chromospheric rapid blueshifted excursions (RBEs) are suggested to be the disk counterparts of type II spicules at the limb and believed to contribute to the coronal heating process. Previous identification of RBEs was mainly based on feature detection using Dopplergrams. In this paper, we study RBEs on 2011 October 21 in a very quiet region at the disk center, which were observed with the high-cadence imaging spectroscopy of the Ca II 8542 A line from the Interferometric Bidimensional Spectrometer (IBIS). By using an automatic spectral analysis algorithm, a total of 98 RBEs are identified during a 11 minute period. Most of these RBEs have either a round or elongated shape, with an average area of 1.2 arcsec^2. The detailed temporal evolution of spectra from IBIS makes possible a quantitative determination of the velocity (~16 km/s) and acceleration (~400 m/s^2) of Ca II 8542 RBEs, and reveal an additional deceleration (~-160 m/s^2) phase that usually follows the initial acceleration. In addition, we also investigate the association of RBEs with the concomitant photospheric magnetic field evolution, using coordinated high-resolution and high-sensitivity magnetograms made by Hinode. Clear examples are found where RBEs appear to be associated with the preceding magnetic flux emergence and/or the subsequent flux cancellation. However, a further analysis with the aid of the Southwest Automatic Magnetic Identification Suite does not yield a significant statistical association between these RBEs and magnetic field evolution. We discuss the implications of our results in the context of understanding the driving mechanism of RBEs.

We present the first part of a project on the global energetics of solar flares and coronal mass ejections (CMEs) that includes about 400 M- and X-class flares observed with AIA and HMI onboard SDO. We calculate the potential energy, free energy, and the flare-dissipated magnetic energy. We calculate these magnetic parameters using two different NLFFF codes: The COR-NLFFF code uses the line-of-sight magnetic field component $B_z$ from HMI to define the potential field, and the 2D coordinates of automatically detected coronal loops in 6 coronal wavelengths from AIA to measure the helical twist of coronal loops caused by vertical currents, while the PHOT-NLFFF code extrapolates the photospheric 3D vector fields. We find agreement between the two codes in the measurement of free energies and dissipated energies within a factor of $ \approx 3$. The size distributions of magnetic parameters exhibit powerlaw slopes that are approximately consistent with the fractal-diffusive self-organized criticality model. The magnetic parameters exhibit scaling laws for the nonpotential energy, $E_{np} \propto E_p^{1.02}$, for the free energy, $E_{free} \propto E_p^{1.7}$ and $E_{free} \propto B_{\varphi}^{1.0} L^{1.5}$, for the dissipated energy, $E_{diss} \propto E_p^{1.6}$ and $E_{diss} \propto E_{free}^{0.9}$, and the energy dissipation volume, $V \propto E_{diss}^{1.2}$. The potential energies vary in the range of $E_p = 1 \times 10^{31} - 4 \times 10^{33}$ erg, while the free energy has a ratio of $E_{free}/E_p \approx 1%-25%$. The Poynting flux amounts to $F_{flare} \approx 5 \times 10^{8} - 10^{10}$ erg cm$^{-2}$ s$^{-1}$ during flares, which averages to $F_{AR} \approx 6 \times 10^6$ erg cm$^{-2}$ s$^{-1}$ during the entire observation period and is comparable with the coronal heating rate requirement in active regions.

Photospheric vector magnetograms from Helioseismic and Magnetic Imager on board the Solar Dynamic Observatory are utilized as the boundary conditions to extrapolate both non-linear force-free and potential magnetic fields in solar corona. Based on the extrapolations, we are able to determine the free magnetic energy (FME) stored in active regions (ARs). Over 3000 vector magnetograms in 61 ARs were analyzed. We compare FME with ARs' flare index (FI) and find that there is a weak correlation ($<60\%$) between FME and FI. FME shows slightly improved flare predictability relative to total unsigned magnetic flux of ARs in the following two aspects: (1) the flare productivity predicted by FME is higher than that predicted by magnetic flux and (2) the correlation between FI and FME is higher than that between FI and magnetic flux. However, this improvement is not significant enough to make a substantial difference in time-accumulated FI, rather than individual flare, predictions.

Large, complex, active regions may produce multiple flares within a certain period of one or two days. These flares could occur in the same location with similar morphologies, commonly referred to as homologous flares. In 2011 September, active region NOAA 11283 produced a pair of homologous flares on the 6th and 7th, respectively. Both of them were white-light (WL) flares, as captured by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory in visible continuum at 617.3 nm which is believed to originate from the deep solar atmosphere.We investigate the WL emission of these X-class flares with HMIs seeing-free imaging spectroscopy. The durations of impulsive peaks in the continuum are about 4 minutes. We compare the WL with hard X-ray (HXR) observations for the September 6 flare and find a good correlation between the continuum and HXR both spatially and temporally. In absence of RHESSI data during the second flare on September 7, the derivative of the GOES soft X-ray is used and also found to be well correlated temporally with the continuum. We measure the contrast enhancements, characteristic sizes, and HXR fluxes of the twin flares, which are similar for both flares, indicating analogous triggering and heating processes. However, the September 7 flare was associated with conspicuous sunquake signals whereas no seismic wave was detected during the flare on September 6. Therefore, this comparison suggests that the particle bombardment may not play a dominant role in producing the sunquake events studied in this paper.

We present the observation of a major solar eruption that is associated with fast sunspot rotation. The event includes a sigmoidal filament eruption, a coronal mass ejection, and a GOES X2.1 flare from NOAA active region 11283. The filament and some overlying arcades were partially rooted in a sunspot. The sunspot rotated at $\sim$10$^\circ$ per hour rate during a period of 6 hours prior to the eruption. In this period, the filament was found to rise gradually along with the sunspot rotation. Based on the HMI observation, for an area along the polarity inversion line underneath the filament, we found gradual pre-eruption decreases of both the mean strength of the photospheric horizontal field ($B_h$) and the mean inclination angle between the vector magnetic field and the local radial (or vertical) direction. These observations are consistent with the pre-eruption gradual rising of the filament-associated magnetic structure. In addition, according to the Non-Linear Force-Free-Field reconstruction of the coronal magnetic field, a pre-eruption magnetic flux rope structure is found to be in alignment with the filament, and a considerable amount of magnetic energy was transported to the corona during the period of sunspot rotation. Our study provides evidences that in this event sunspot rotation plays an important role in twisting, energizing, and destabilizing the coronal filament-flux rope system, and led to the eruption. We also propose that the pre-event evolution of $B_h$ may be used to discern the driving mechanism of eruptions.

This Letter reports two rarely observed three-ribbon flares (M1.9 and C9.2) on 2012 July 6 in NOAA AR 11515, which we found with Halpha observations of 0.1" resolution from the New Solar Telescope and CaII H images from Hinode. The flaring site is characterized with an intriguing "fish-bone-like" morphology evidenced by both Halpha images and a nonlinear force-free field (NLFFF) extrapolation, where two semi-parallel rows of low-lying, sheared loops connect an elongated, parasitic negative field with the sandwiching positive fields. The NLFFF model also shows that the two rows of loops are asymmetric in height and have opposite twists, and are enveloped by large-scale field lines including open fields. The two flares occurred in succession in half an hour and are located at the two ends of the flaring region. The three ribbons of each flare run parallel to the PIL, with the outer two lying in the positive field and the central one in the negative field. Both flares show surge-like flows in Halpha apparently toward the remote region, while the C9.2 flare is also accompanied by EUV jets possibly along the open field lines. Interestingly, the 12-25 keV hard X-ray sources of the C9.2 flare first line up with the central ribbon then shift to concentrate on the top of the higher branch of loops. These results are discussed in favor of reconnection along the coronal null-line producing the three flare ribbons and the associated ejections.

Rapid, irreversible changes of magnetic topology and sunspot structure associated with flares have been systematically observed in recent years. The most striking features include the increase of horizontal field at the polarity inversion line (PIL) and the co-spatial penumbral darkening. A likely explanation of the above phenomenon is the back reaction to the coronal restructuring after eruptions: a coronal mass ejection carries the upward momentum while the downward momentum compresses the field lines near the PIL. Previous studies could only use low resolution (above 1") magnetograms and white-light images. Therefore, the changes are mostly observed for X-class flares. Taking advantage of the 0.1" spatial resolution and 15s temporal cadence of the New Solar Telescope at Big Bear Solar Observatory, we report in detail the rapid formation of sunspot penumbra at the PIL associated with the C7.4 flare on 2012 July 2. It is unambiguously shown that the solar granulation pattern evolves to alternating dark and bright fibril structure, the typical pattern of penumbra. Interestingly, the appearance of such a penumbra creates a new delta sunspot. The penumbral formation is also accompanied by the enhancement of horizontal field observed using vector magnetograms from the Helioseismic and Magnetic Imager. We explain our observations as due to the eruption of a flux rope following magnetic cancellation at the PIL. Subsequently the re-closed arcade fields are pushed down towards the surface to form the new penumbra. NLFFF extrapolation clearly shows both the flux rope close to the surface and the overlying fields.

We present an unprecedented high-resolution \ha\ imaging spectroscopic observation of a C4.1 flare taken with IBIS on 2011 October 22. The flare consists of a main circular ribbon that occurred in a parasitic magnetic configuration and a remote ribbon that was observed by the IBIS. Such a circular-ribbon flare with a remote brightening is predicted in 3D fan-spine reconnection but so far has been rarely observed. During the flare impulsive phase, we define "core" and "halo" structures in the observed ribbon. Examining the \ha\ emission spectra averaged in the flare core and halo areas, we find that only those from the flare cores show typical nonthermal electron beam heating characteristics. These characteristics include: broad and centrally reversed emission spectra, excess emission in the red wing with regard to the blue wing (i.e., red asymmetry), and redshifted bisectors of the emission spectra. We also observe rather quick timescales for the heating (30 s) and cooling (14--33 s) in the flare core locations. Therefore, we suggest that the flare cores revealed by IBIS track the sites of electron beam precipitation with exceptional spatial and temporal resolution. The flare cores show two-stage motion (a parallel motion along the ribbon followed by an expansion motion perpendicular to the ribbon) during the two impulsive phases of the flare. Some cores jump quickly (30 \kms) between discrete magnetic elements implying reconnection involving different flux tubes. We observe a very high temporal correlation ($\gtrsim0.9$) between the integrated \ha\ and HXR emission during the flare impulsive phase. A short time delay (4.6 s) is also found in the \ha\ emission spikes relative to HXR bursts. The ionization timescale of the cool chromosphere and the extra time taken for the electrons to travel to the remote ribbon site may contribute to this delay.

Numerical simulations suggest that kink and torus instabilities are two potential contributors to the initiation and prorogation of eruptive events. A magnetic parameter named decay index (i.e., the coronal magnetic gradient of the overlying fields above the eruptive flux ropes) could play an important role in controlling kinematics of eruptions. Previous studies have identified a threshold range of the decay index that distinguishes between eruptive and confined configurations. Here we advance the study by investigating if there is a clear correlation between the decay index and CME speed. 38 CMEs associated with filament eruptions and/or two-ribbon flares are selected using the Halpha data from the Global Halpha Network. The filaments and flare ribbons observed in Halpha associated with the CMEs help to locate the magnetic polarity inversion line, along which the decay index is calculated based on the potential field extrapolation using MDI magnetograms as boundary conditions. The speeds of CMEs are obtained from the LASCO C2 CME catalog available online. We find that the mean decay index increases with CME speed for those CMEs with a speed below 1000 km/s, and stays flat around 2.2 for the CMEs with higher speeds. In addition, we present a case study of a partial filament eruption, in which the decay indexes show different values above the erupted/non-erupted part.

We study the strong gravitational lensing on the equatorial plane of a quasi-Kerr compact object with arbitrary quadrupole moments which can be used to model the super-massive central object of the galaxy. We find that, when the quadrupolar correction parameter $\xi$ takes the positive (negative) value, the photon-sphere radius $r_{ps}$, the minimum impact parameter $u_{ps}$, the coefficient $\bar{b}$, the relative magnitudes $r_m$ and the angular position of the relativistic images $\theta_{\infty}$ are larger (smaller) than the results obtained in the Kerr black hole, but the coefficient $\bar{a}$, the deflection angle $\alpha(\theta)$ and the angular separation $s$ are smaller (larger) than that in the Kerr black hole. These features may offer a way to probe special properties for some rotating compact objects by the astronomical instruments in the future.

The properties of the ergosphere and energy extraction by Penrose process in a rotating non-Kerr black hole are investigated. It is shown that the ergosphere is sensitive to the deformation parameter $\epsilon$ and the shape of the ergosphere becomes thick with increase of the parameter $\epsilon$. It is of interest to note that, comparing with the Kerr black hole, the deformation parameter $\epsilon$ can enhance the maximum efficiency of the energy extraction process greatly. Especially, for the case of $a>M$, the non-Kerr metric describes a superspinning compact object and the maximum efficiency can exceed 60%, while it is only 20.7% for the extremal Kerr black hole.

We study the strong gravitational lensing in the background of a rotating non-Kerr compact object with a deformed parameter $\epsilon$ and an unbound rotation parameter $a$. We find that the photon sphere radius and the deflection angle depend sharply on the parameters $\epsilon$ and $a$. For the case in which the black hole is more prolate than a Kerr black hole, the photon sphere exists only in the regime $\epsilon\leq\epsilon_{max}$ for prograde photon. The upper limit $\epsilon_{max}$ is a function of the rotation parameter $a$. As $\epsilon>\epsilon_{max}$, the deflection angle of the light ray closing very to the naked singularity is a positive finite value, which is different from those in both the usual Kerr black hole spacetime and in the rotating naked singularity described by Janis-Newman-Winicour metric. For the oblate black hole and the retrograde photon, there does not exist such a threshold value. Modelling the supermassive central object of the Galaxy as a rotating non-Kerr compact object, we estimated the numerical values of the coefficients and observables for gravitational lensing in the strong field limit.

We study the accretion process in the thin disk around a rotating non-Kerr black hole with a deformed parameter and an unbound rotation parameter. Our results show that the presence of the deformed parameter $\epsilon$ modifies the standard properties of the disk. For the case in which the black hole is more oblate than a Kerr black hole, the larger deviation leads to the smaller energy flux, the lower radiation temperature and the fainter spectra luminosity in the disk. For the black hole with positive deformed parameter, we find that the effect of the deformed parameter on the disk becomes more complicated. It depends not only on the rotation direction of the black hole and the orbit particles, but also on the sign of the difference between the deformed parameter $\epsilon$ and a certain critical value $\epsilon_{c}$. These significant features in the mass accretion process may provide a possibility to test the no-hair theorem in the strong field regime in future astronomical observations.

The rapid, irreversible change of the photospheric magnetic field has been recognized as an important element of the solar flare process. This Letter reports such a rapid change of magnetic fields during the 2011 February 13 M6.6 flare in NOAA AR 11158 that we found from the vector magnetograms of the Helioseismic and Magnetic Imager with 12-min cadence. High-resolution magnetograms of Hinode that are available at ~-5.5, -1.5, 1.5, and 4 hrs relative to the flare maximum are used to reconstruct three-dimensional coronal magnetic field under the nonlinear force-free field (NLFFF) assumption. UV and hard X-ray images are also used to illuminate the magnetic field evolution and energy release. The rapid change is mainly detected by HMI in a compact region lying in the center of the magnetic sigmoid, where the mean horizontal field strength exhibited a significant increase by 28%. The region lies between the initial strong UV and hard X-ray sources in the chromosphere, which are cospatial with the central feet of the sigmoid according to the NLFFF model. The NLFFF model further shows that strong coronal currents are concentrated immediately above the region, and that more intriguingly, the coronal current system underwent an apparent downward collapse after the sigmoid eruption. These results are discussed in favor of both the tether-cutting reconnection producing the flare and the ensuing implosion of the coronal field resulting from the energy release.

We have studied the quasinormal modes of a massive scalar field coupling to Einstein's tensor in the spacelike stretched $AdS_3$ black hole spacetime. We find that both the right-moving and left-moving quasinormal frequencies depend not only on the warped parameter $v$ of black hole, but also on the coupling between the scalar field and Einstein's tensor. Moreover, we also discuss the warped $AdS/CFT$ correspondence from the quasinormal modes and probe the effects of the coupling on the left and right conformal weights $h_L$ and $h_R$ of the operators dual to the scalar field in the boundary.

We study the accretion process in the thin disk around a squashed Kaluza-Klein black hole and probe the effects of the extra dimensional scale $\rho_0$ on the physical properties of the disk. Our results show that with the increase of the parameter $\rho_0$, the energy flux, the conversion efficiency, the radiation temperature, the spectra luminosity and the spectra cut-off frequency of the thin accretion disk decrease, but the inner border of the disk increases. This implies that the extra dimension scale imprints in the mass accretion process in the disk.

We study the influence of the spacetime noncommutative parameter on the strong field gravitational lensing in the noncommutative Reissner-Nordstr\"{o}m black-hole spacetime. Supposing that the gravitational field of the supermassive central object of the Galaxy is described by this metric, we estimate the numerical values of the coefficients and observables for strong gravitational lensing. Our results show that with the increase of the parameter $\sqrt{\vartheta}$, the observables $\theta_{\infty}$ and $r_m$ decrease, while $s$ increases. Our results also show that i) if $\sqrt{\vartheta}$ is strong, the observables are close to those of the noncommutative Schwarzschild black hole lensing; ii) if $\sqrt{\vartheta}$ is weak, the observables are close to those of the commutative Reissner-Nordstr\"{o}m black hole lensing; iii) the detectable scope of $\vartheta$ in a noncommutative Reissner-Nordstr\"{o}m black hole lensing is $0.12\leq\sqrt{\vartheta}\leq0.26$, which is wider than that in a noncommutative Schwarzschild black hole lensing, $0.18\leq\sqrt{\vartheta}\leq0.26$. This may offer a way to probe the spacetime noncommutative constant $\vartheta$ by the astronomical instruments in the future.

We study the collision of two particles with the different rest masses moving in the equatorial plane of a Kerr-Taub-NUT spacetime and get the center-of-mass (CM) energy for the particles. We find that the CM energy depends not only on the rotation parameter, $a$, but also on the NUT charge of the Kerr-Taub-NUT spacetime, $n$. Especially, for the extremal Kerr-Taub-NUT spacetime, an unlimited CM energy can be approached if the parameter $a$ is in the range $[1,\sqrt{2}]$, which is different from that of the Kerr and Kerr-Newman black holes.

We use rotation stereoscopy to estimate the height of a steady-state solar feature relative to the photosphere, based on its apparent motion in the image plane recorded over several days of observation. The stereoscopy algorithm is adapted to work with either one- or two-dimensional data (i.e. from images or from observations that record the projected position of the source along an arbitrary axis). The accuracy of the algorithm is tested on simulated data, and then the algorithm is used to estimate the coronal radio source heights associated with the active region NOAA 10956, based on multifrequency imaging data over 7 days from the Siberian Solar Radio Telescope near 5.7 GHz, the Nobeyama Radio Heliograph at 17 GHz, as well as one-dimensional scans at multiple frequencies spanning the 5.98--15.95 GHz frequency range from the RATAN-600 instrument. The gyroresonance emission mechanism, which is sensitive to the coronal magnetic field strength, is applied to convert the estimated radio source heights at various frequencies, h(f), to information about magnetic field vs. height B(h), and the results are compared to a magnetic field extrapolation derived from photospheric magnetic field observations obtained by Hinode and MDI. We found that the gyroresonant emission comes from the heights exceeding location of the third gyrolayer irrespectively on the magnetic extrapolation method; implications of this finding for the coronal magnetography and coronal plasma physics are discussed.

We investigate the strong gravitational lensing in a squashed Kaluza-Klein black hole immersed in the G\"{o}del universe with global rotation. Our result show that the strong gravitational lensing in the squashed Kaluza-Klein G\"{o}del black hole spacetime has some distinct behaviors from that in the Kerr case. In the squashed Kaluza-Klein G\"{o}del black hole spacetime, the photon sphere radius, the minimum impact parameter, the coefficient $\bar{a}$, $\bar{b}$ and the deflection angle $\alpha(\theta)$ in the $\phi$ direction are independent of whether the photon goes with or against the global rotation of the G\"{o}del Universe. While in the Kerr black hole, the values of these quantities for the prograde photons are different from those for the retrograde photons. Moreover, the coefficient of $\bar{b}$ increases with $j$ in the squashed Kaluza-Klein G\"{o}del black hole, but decreases with $a$ in the Kerr case. We also probe the influence of the squashed effect on the strong gravitational lensing in this black hole and find that in the extremely squashed case $\rho_0=0$, the coefficient $\bar{a}$ is a constant 1 and is independent of the global rotation of the G\"{o}del Universe. Furthermore, we assume that the gravitational field of the supermassive central object of the Galaxy can be described by this metric and estimate the numerical values of the coefficients and the main observables in the strong gravitational lensing.

We present the first interferometric observation of a zebra-pattern radio burst with simultaneous high spectral (~ 1 MHz) and high time (20 ms) resolution. The Frequency-Agile Solar Radiotelescope (FASR) Subsystem Testbed (FST) and the Owens Valley Solar Array (OVSA) were used in parallel to observe the X1.5 flare on 14 December 2006. By using OVSA to calibrate the FST the source position of the zebra pattern can be located on the solar disk. With the help of multi-wavelength observations and a nonlinear force-free field (NLFFF) extrapolation, the zebra source is explored in relation to the magnetic field configuration. New constraints are placed on the source size and position as a function of frequency and time. We conclude that the zebra burst is consistent with a double-plasma resonance (DPR) model in which the radio emission occurs in resonance layers where the upper hybrid frequency is harmonically related to the electron cyclotron frequency in a coronal magnetic loop.

We study collision of two general geodesic particles around the Kerr-Newman (KN) black hole and get the center-of-mass (CM) energy of the non-marginally and marginally bound critical particles in the direct collision and LSO collision scenarios. We find the constraint conditions that arbitrarily high CM energy can be obtained for the near-horizon collision of two general geodesic particles in the extremal KN black hole, and note that the charge decreases the value of the latitude in which arbitrarily high CM energy can occurs. We also interpret why the high-velocity collision belt centers to the equator with the increase of the charge.

Noncommutative geometry may be a starting point to a quantum gravity. We study the influence of the spacetime noncommutative parameter on the strong field gravitational lensing in the noncommutative Schwarzschild black-hole spacetime and obtain the angular position and magnification of the relativistic images. Supposing that the gravitational field of the supermassive central object of the galaxy described by this metric, we estimate the numerical values of the coefficients and observables for strong gravitational lensing. Comparing to the Reissner-Norstr\"{om} black hole, we find that the influences of the spacetime noncommutative parameter is similar to those of the charge, just these influences are much smaller. This may offer a way to distinguish a noncommutative black hole from a Reissner-Norstr\"{om} black hole, and may probe the spacetime noncommutative constant $\vartheta$ \cite{foot} by the astronomical instruments in the future.

We study the second-order phase transition (SOPT) for the spherically symmetric Kehagias-Sfetsos (KS) black hole in the deformed H\v{o}rava-Lifshitz gravity by applying the methods of equilibrium and non-equilibrium fluctuations. We find that, although the KS black hole has only one mass parameter as the usual Schwarzschild ones, the SOPT will take place if the mass of the KS black hole changes across the critical point $\frac{\sqrt{5+\sqrt{33}}(\sqrt{33}-1)} {16\sqrt{\omega}} $. The result show us that there is difference between the H\v{o}rava-Lifshitz gravity and the Einstein's gravity theory.

Filament eruptions and hard X-ray (HXR) source motions are commonly observed in solar flares, which provides critical information on the coronal magnetic reconnection. This Letter reports an event on 2005 January 15, in which we found an asymmetric filament eruption and a subsequent coronal mass ejection together with complicated motions of HXR sources during the GOES-class X2.6 flare. The HXR sources initially converge to the magnetic polarity inversion line (PIL), and then move in directions either parallel or perpendicular to the PIL depending on the local field configuration. We distinguish the evolution of the HXR source motion in four phases and associate each of them with distinct regions of coronal magnetic fields as reconstructed using a non-linear force-free field extrapolation. It is found that the magnetic reconnection proceeds along the PIL toward the regions where the overlying field decreases with height more rapidly. It is also found that not only the perpendicular but the parallel motion of the HXR sources correlates well with the HXR lightcurve. These results are discussed in favor of the torus instability as an important factor in the eruptive process.

To study the three-dimensional (3D) magnetic field topology and its long-term evolution associated with the X3.4 flare of 2006 December 13, we investigate the coronal relative magnetic helicity in the flaring active region (AR) NOAA 10930 during the time period of December 8-14. The coronal helicity is calculated based on the 3D nonlinear force-free magnetic fields reconstructed by the weighted optimization method of Wiegelmann, and is compared with the amount of helicity injected through the photospheric surface of the AR. The helicity injection is determined from the magnetic helicity flux density proposed by Pariat et al. using Solar and Heliospheric Observatory/Michelson Doppler Imager magnetograms. The major findings of this study are the following. (1) The time profile of the coronal helicity shows a good correlation with that of the helicity accumulation by injection through the surface. (2) The coronal helicity of the AR is estimated to be -4.3$\times$10$^{43}$ Mx$^{2}$ just before the X3.4 flare. (3) This flare is preceded not only by a large increase of negative helicity, -3.2$\times$10$^{43}$ Mx$^{2}$, in the corona over ~1.5 days but also by noticeable injections of positive helicity though the photospheric surface around the flaring magnetic polarity inversion line during the time period of the channel structure development. We conjecture that the occurrence of the X3.4 flare is involved with the positive helicity injection into an existing system of negative helicity.

A magnetic channel - a series of polarity reversals separating elongated flux threads with opposite polarities - may be a manifestation of a highly non-potential magnetic configuration in active regions. To understand its formation we have carried out a detailed analysis of the magnetic channel in AR 10930 using data taken by the Solar Optical Telescope/Hinode. As a result, we found upflows (-0.5 to -1.0 km/s) and downflows (+1.5 to +2.0 km/s) inside and at both tips of the thread respectively, and a pair of strong vertical currents of opposite polarity along the channel. Moreover, our analysis of the nonlinear force-free fields constructed from the photospheric magnetic field indicates that the current density in the lower corona may have gradually increased as a result of the continuous emergence of the highly sheared flux along the channel. With these results, we suggest that the magnetic channel originates from the emergence of the twisted flux tube that has formed below the surface before the emergence.

In this letter, we study in detail the evolution of the post-flare loops on 2005 January 15 that occurred between two consecutive solar eruption events, both of which generated a fast halo CME and a major flare. The post-flare loop system, formed after the first CME/flare eruption, evolved rapidly, as manifested by the unusual accelerating rise motion of the loops. Through nonlinear force-free field (NLFFF) models, we obtain the magnetic structure over the active region. It clearly shows that the flux rope below the loops also kept rising accompanied with increasing twist and length. Finally, the post-flare magnetic configuration evolved to a state that resulted in the second CME/flare eruption. This is an event in which the post-flare loops can re-flare in a short period of $\sim$16 hr following the first CME/flare eruption. The observed re-flaring at the same location is likely driven by the rapid evolution of the flux rope caused by the magnetic flux emergence and the rotation of the sunspot. This observation provides valuable information on CME/flare models and their prediction.

We have investigated the variation of magnetic helicity over a span of several days around the times of 11 X-class flares which occurred in seven active regions (NOAA 9672, 10030, 10314, 10486, 10564, 10696, and 10720) using the magnetograms taken by the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO). As a major result we found that each of these major flares was preceded by a significant helicity accumulation over a long period (0.5 to a few days). Another finding is that the helicity accumulates at a nearly constant rate and then becomes nearly constant before the flares. This led us to distinguish the helicity variation into two phases: a phase of monotonically increasing helicity and the following phase of relatively constant helicity. As expected, the amount of helicity accumulated shows a modest correlation with time-integrated soft X-ray flux during flares. However, the average helicity change rate in the first phase shows even stronger correlation with the time-integrated soft X-ray flux. We discuss the physical implications of this result and the possibility that this characteristic helicity variation pattern can be used as an early warning sign for solar eruptions.

The relation between the quasinormal modes (QNMs) and the second order thermodynamic phase transition (SOTPT) for the Reissner-Nordstr\"om (RN) black hole is studied. It is shown that the quasinormal frequencies of the RN black hole start to get a spiral-like shape in the complex $\omega$ plane and both the real and imaginary parts become the oscillatory functions of the charge if the real part of the quasinormal frequencies arrives at its maximum at the second order phase transition point of Davies for given overtone number and angular quantum number. That is to say, we can find out the SOTPT point from the QNMs of the RN black hole. The fact shows that the quasinormal frequencies carry the thermodynamical information of the RN black hole.

The quasinormal modes (QNMs) and the late-time behavior of arbitrary spin fields are studied in the background of a Schwarzschild black hole with a global monopole (SBHGM). It has been shown that the real part of the QNMs for a SBHGM decreases as the symmetry breaking scale parameter $H$ increases but imaginary part increases instead. For large overtone number $n$, these QNMs become evenly spaced and the spacing for the imaginary part equals to $-i(1-H)^{3/2}/(4M)$ which is dependent of $H$ but independent of the quantum number $l$. It is surprisingly found that the late-time behavior is dominated by an inverse power-law tail $t^{-2[1+\sqrt{(s+1/2)^{2}+ (l-s)(l+s+1)/(1-H)}]}$ for each $l$, and as $H\to0$ it reduces to the Schwarzschild case $t^{-(2l+3)}$ which is independent of the spin number $s$.

The rare phenomenon of ribbon-like hard X-ray (HXR) sources up to 100 keV found in the 2005 May 13 M8.0 flare observed with the \textit{Reuven Ramaty High Energy Solar Spectroscopic Imager} provides detailed information on the spatial distribution of flare HXR emission. In this Letter, we further investigate the characteristics of HXR emission in this event using imaging spectroscopy, from which we obtain spatially resolved HXR spectral maps during the flare impulsive phase. As a result we found, along a flare ribbon, an anticorrelation relationship between the local HXR flux and the local HXR spectral index. We suggest that this can be regarded as a spatial analog of the well-known temporal soft-hard-soft spectral evolution pattern of the integrated HXR flux. We also found an anticorrelation between HXR spectral index and local electric field along the ribbon, which suggests the electron acceleration by the electric field during flares.

The Dirac quasinormal modes (QNMs) of the Kerr-Newman black hole are investigated using continued fraction approach. It is shown that the quasinormal frequencies in the complex $\omega$ plane move counterclockwise as the charge or angular momentum per unit mass of the black hole increases. They get a spiral-like shape, moving out of their Schwarzschild or Reissner-Nordstr\"om values and "looping in" towards some limiting frequencies as the charge and angular momentum per unit mass tend to their extremal values. The number of the spirals increases as the overtone number increases but decreases as the angular quantum number increases. It is also found that both the real and imaginary parts are oscillatory functions of the angular momentum per unit mass, and the oscillation becomes faster as the overtone number increases but slower as the angular quantum number increases.

The late-time evolution of the charged massive Dirac fields in the background of a Reissner-Norstr\"om (RN) black hole is studied. It is found that the intermediate late-time behavior is dominated by an inverse power-law decaying tail without any oscillation in which the dumping exponent depends not only on the multiple number of the wave mode but also on the field parameters. It is also found that, at very late times, the oscillatory tail has the decay rate of $t^{-5/6}$ and the oscillation of the tail has the period $2\pi/\mu$ which is modulated by two types of long-term phase shifts.

The quasinormal modes (QNMs) associated with the decay of Dirac field perturbation around a Schwarzschild black hole is investigated by using continued fraction and Hill-determinant approaches. It is shown that the fundamental quasinormal frequencies become evenly spaced for large angular quantum number and the spacing is given by $\omega_{\lambda+1}- \omega_{\lambda}=0.38490-0.00000i$. The angular quantum number has the surprising effect of increasing real part of the quasinormal frequencies, but it almost does not affect imaginary part, especially for low overtones. In addition, the quasinormal frequencies also become evenly spaced for large overtone number and the spacing for imaginary part is $Im(\omega_{n+1})-Im(\omega_n)\approx -i/4M$ which is same as that of the scalar, electromagnetic, and gravitational perturbations.

The quasinormal modes (QNMs) of Dirac field perturbations of a Reissner-Nordstr\"om black hole in an asymptotically Anti-de Sitter spacetime are investigated. We find that both the real and imaginary parts of the fundamental quasinormal frequencies for large black holes are the linear functions of the Hawking temperature, and the slope of the lines for the real parts decreases while that for the magnitude of the imaginary parts increases as the black hole charge increases. According to the Anti-de Sitter/Conformal Field Theory (AdS/CFT) correspondence, the fact shows that different charge presents different time scale in three-dimensional CFT. Another interesting result is that the quasinormal frequencies become evenly spaced for high overtone number, and in the spacing expressions the real part decreases while the magnitude of the imaginary part increases as the charge increases. We also study the relation between quasinormal frequencies and angular quantum number and find that the real part increases while the magnitude of the imaginary part decreases as the angular quantum number increases.

The quasinormal modes (QNMs) associated with the decay of Dirac field perturbation around a Schwarzschild-anti-de Sitter (SAdS) black hole is investigated by using Horowitz-Hubeny approach. We find that both the real and the imaginary parts of the fundamental quasinormal frequencies for large black holes are linear functions of the Hawking temperature, and the real part of the fundamental quasinormal frequencies for intermediate and small black holes approximates a temperature curve and the corresponding imaginary part is almost a linear function of the black hole radius. The angular quantum number has the surprising effect of increasing the damping time scale and decreasing the oscillation time scale for intermediate and small black holes. We also find that the rates between quasinormal frequencies and black hole radii for large, intermediate and small black holes are linear functions of the overtone number, i. e., all the modes are evenly spaced.

We investigate the quasinormal mode frequencies for the massless Dirac field in static four dimensional $AdS$ space-time. The separation of the Dirac equation is achieved for the first time in $AdS$ space. Besides the relevance that this calculation can have in the framework of the $AdS/CFT$ correspondence between M-theory on $AdS_4\times S^7$ and SU(N) super Yang-Mills theory on $M_3$, it also serves to fill in a gap in the literature, which has only been concerned with particles of integral spin $0,1,2$.

The late-time tail behavior of massive Dirac fields is investigated in the Schwarzschild black-hole geometry and the result is compared with that of the massive scalar fields. It is shown that in the intermediate times there are three kinds of differences between the massive Dirac and scalar fields, (I) the asymptotic behavior of massive Dirac fields is dominated by a decaying tail without any oscillation, but the massive scalar field by a oscillatory inverse power-law decaying tail, (II) the dumping exponent for the massive Dirac field depends not only on the multiple number of the wave mode but also on the mass of the Dirac field, while that for the massive scalar field depends on the multiple number only, and (III) the decay of the massive Dirac field is slower than that of the massive scalar field.

The late-time tail behavior of massive scalar fields is studied analytically in a stationary axisymmetric EMDA black hole geometry. It is shown that the asymptotic behavior of massive perturbations is dominated by the oscillatory inverse power-law decaying tail $ t^{-(l+3/2)}\sin(\mu t)$ at the intermediate late times, and by the asymptotic tail $ t^{-5/6}\sin(\mu t)$ at asymptotically late times. Our result seems to suggest that the intermediate tails $ t^{-(l+3/2)}\sin(\mu t)$ and the asymptotically tails $t^{-5/6} \sin(\mu t)$ may be quite general features for evolution of massive scalar fields in any four dimensional asymptotically flat rotating black hole backgrounds.

The quasinormal modes of the Reissner-Nordstr\"om de Sitter black hole for the massless Dirac fields are studied using the P\"oshl-Teller potential approximation. We find that the magnitude of the imaginary part of the quasinormal frequencies decreases as the cosmological constant or the orbital angular momentum increases, but it increases as the charge or the overtone number increases. An interesting feature is that the imaginary part is almost linearly related to the real part as the cosmological constant changes for fixed charge, and the linearity becomes better as the orbital angular momentum increases. We also prove exactly that the Dirac quasinormal frequencies are the same for opposite chirality.

By using the 't Hooft's "brick wall" model and the Pauli-Villars regularization scheme we calculate the statistical-mechanical entropies arising from the quantum scalar field in different coordinate settings, such as the Painlev\'{e} and Lemaitre coordinates. At first glance, it seems that the entropies would be different from that in the standard Schwarzschild coordinate since the metrics in both the Painlev\'{e} and Lemaitre coordinates do not possess the singularity at the event horizon as that in the Schwarzschild-like coordinate. However, after an exact calculation we find that, up to the subleading correction, the statistical-mechanical entropies in these coordinates are equivalent to that in the Schwarzschild-like coordinate. The result is not only valid for black holes and de Sitter spaces, but also for the case that the quantum field exerts back reaction on the gravitational field provided that the back reaction does not alter the symmetry of the spacetime.