Theoretically, CME kinematics are related to magnetic reconnection processes in the solar corona. However, the current quantitative understanding of this relationship is based on the analysis of only a handful of events. Here we report a statistical study of 60 CME-flare events from August 2010 to December 2013. We investigate kinematic properties of CMEs and magnetic reconnection in the low corona during the early phase of the eruptions, by combining limb observations from STEREO with simultaneous on-disk views from SDO. For a subset of 42 events with reconnection rate evaluated by the magnetic fluxes swept by the flare ribbons on the solar disk observed from SDO, we find a strong correlation between the peak CME acceleration and the peak reconnection rate. Also, the maximum velocities of relatively fast CMEs (> 600 km/s) are positively correlated with the reconnection flux, but no such correlation is found for slow CMEs. A time-lagged correlation analysis suggests that the distribution of the time lag of CME acceleration relative to reconnection rate exhibits three peaks, approximately 10 minutes apart, and on average, acceleration-lead events have smaller reconnection rates. We further compare the CME total mechanical energy with the estimated energy in the current sheet. The comparison suggests that, for small-flare events, reconnection in the current sheet alone is insufficient to fuel CMEs. Results from this study suggest that flare reconnection may dominate the acceleration of fast CMEs, but for events of slow CMEs and weak reconnection, other mechanisms may be more important.

Aims. We investigate the physical and chemical conditions of molecular gas in the circumnuclear disk (CND) region of NGC 1068. Methods. We carried out a spectral line survey with the IRAM 30m telescope toward the center of NGC 1068 and mainly focused on the 2 mm band with a frequency coverage of 160.7-168.6 GHz and 176.5-184.3 GHz. Results. Fifteen lines are detected in NGC 1068, eight of which are new detections for this galaxy. We derive the rotation temperatures and column densities of fourteen molecular species. Conclusions. Based on the [HCO+ (2 - 1)]/[HOC+ (2 - 1)] ratio, we obtain a high ionization degree in the CND of NGC 1068. It is found that HC3N is concentrated in the east knot, while 13CCH, CH3CN, SO, HOC+, CS, CH3CCH, and H2CO are concentrated in the west knot. Compared to the star-forming galaxies M 82 and NGC 253, the chemistry of NGC 1068 might be less strongly affected by the UV radiation field, and its kinetic temperature might be lower.

In past decades, much progress has been achieved on the origin and evolution of coronal mass ejections (CMEs). In-situ observations of the counterparts of CMEs, especially magnetic clouds (MCs) near the Earth, have provided measurements of the structure and total flux of CME flux ropes. However, it has been difficult to measure these properties in the erupting CME flux rope, in particular in the pre-existing flux rope. In this work, we propose a model to estimate the toroidal flux of the pre-existing flux rope by subtracting the flux contributed by magnetic reconnection during the eruption from the flux measured in the MC. The flux by the reconnection is derived from geometric properties of two-ribbon flares based on a quasi-2D reconnection model. We then apply the model to four CME/flare events and find that the ratio of toroidal flux in the pre-existing flux rope to that of the associated MC lies in the range of 0.40--0.88. It indicates that the toroidal flux of the pre-existing flux rope has an important contribution to that of the CME flux rope and is usually at least as large as the flux arising from the eruption process for the selected events.

It is well accepted that a magnetic flux rope (MFR) is a critical component of many coronal mass ejections (CMEs), yet how it evolves toward eruption remains unclear. Here we investigate the continuous evolution of a pre-existing MFR, which is rooted in strong photospheric magnetic fields and electric currents. The evolution of the MFR is observed by the Solar Terrestrial Relations Observatory (STEREO) and the Solar Dynamics Observatory (SDO) from multiple viewpoints. From STEREO's perspective, the MFR starts to rise slowly above the limb five hours before it erupts as a halo CME on 2012 June 14. In SDO observations, conjugate dimmings develop on the disk, simultaneously with the gradual expansion of the MFR, suggesting that the dimmings map the MFR's feet. The evolution comprises a two-stage gradual expansion followed by another stage of rapid acceleration/eruption. Quantitative measurements indicate that magnetic twist of the MFR increases from 1.0 +/- 0.5 to 2.0 +/- 0.5 turns during the five-hour expansion, and further increases to about 4.0 turns per AU when detected as a magnetic cloud at 1 AU two day later. In addition, each stage is preceded by flare(s), implying reconnection is actively involved in the evolution and eruption of the MFR. The implications of these measurements on the CME initiation mechanisms are discussed.

This is the eighth in a series of papers reporting on a large reverberation mapping campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). We employ the recently developed dynamical modeling approach for broad-line regions (BLRs) based on the method of Pancoast et al. to analyze the reverberation mapping dataset of Mrk 142 observed in the first monitoring season. In this approach, continuum variations are reconstructed using a damped random walk process, and BLR structure is delineated using a flexible disk-like geometry, in which BLR clouds move around the central BH with Keplerian orbits or inflow/outflow motion. The approach also includes the possibilities of anisotropic emission of BLR clouds, non-linear response of the line emission to the continuum, and different long-term trends in the continuum and emission-line variations. We implement the approach in a Bayesian framework that is apt for parallel computation and use a Markov Chain Monte Carlo technique to recover the parameters and uncertainties for the modeling, including mass of the central BH. We apply three BLR models with different prescriptions of BLR clouds distributions and find that the best model for fitting the data of Mrk 142 is a two-zone BLR model, consistent with the theoretical BLR model surrounding slim accretion disks. The best model yields a BH mass of $\log (M_\bullet/M_\odot)=6.23_{-0.45}^{+0.26}$, resulting in a virial factor of $\log f=-0.36_{-0.54}^{+0.33}$ for the full width at half maximum of the H$\beta$ line measured from the mean spectrum. The virial factors for the other measures of the H$\beta$ line width are also presented.

As one of the series of papers reporting on a large reverberation mapping campaign, we apply the maximum entropy method (MEM) to 9 narrow-line Seyfert 1 galaxies with super-Eddington accretion rates observed during 2012-2013 for the velocity-delay maps of their H\beta\ and H\gamma\ emission lines. The maps of 6 objects are reliably reconstructed using MEM. The maps of H\beta\ and H\gamma\ emission lines of Mrk 335 indicate that the gas of its broad-line region (BLR) is infalling. For Mrk 142, its H\beta\ and H\gamma\ lines show signatures of outflow. The H\beta\ and H\gamma\ maps of Mrk 1044 demonstrate complex kinematics -- a virialized motion accompanied by an outflow signature, and the H\beta\ map of IRAS F12397+3333 is consistent with a disk or a spherical shell. The H\beta\ maps of Mrk 486 and MCG +06-26-012 suggest the presence of an inflow and outflow, respectively. These super-Eddington accretors show diverse geometry and kinematics. Brief discussions of their BLRs are provided for each individual object.

Spectroscopic observations of a solar limb flare recorded by SUMER on SOHO reveal, for the first time, hot fast magnetic reconnection outflows in the corona. As the reconnection site rises across the SUMER spectrometer slit, significant blue- and red-shift signatures are observed in sequence in the Fe XIX line, reflecting upflows and downflows of hot plasma jets, respectively. With the projection effect corrected, the measured outflow speed is between 900-3500 km/s, consistent with theoretical predictions of the Alfvenic outflows in magnetic reconnection region in solar impulsive events. Based on theoretic models, the magnetic field strength near the reconnection region is estimated to be 19-37 Gauss.

We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where the dense gas is traced by the HCN J=4-3 and the HCO^+ J=4-3 emission, measured in our sample of spatially-resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, L_IR/L'_dense, shows systematic variations with L_IR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between L_IR/L'_dense ratio and the warm-dust temperature gauged by the 70 \mu m/100 \mu m flux ratio. We find the luminosity ratios of IR/HCN(4-3) and IR/HCO^+(4-3), which can be taken as a proxy for the efficiency of star formation in the dense molecular gas (SFE_dense), appears to be nearly independent of the dense-gas fraction (f_dense) for our sample of galaxies. The SFE of the total molecular gas (SFE_mol) is found to increase substantially with f_dense when combining our data with that on local (ultra)luminous infrared galaxies and high-z quasars. The mean L'_HCN(4-3)/L'_HCO^+(4-3) line ratio measured for the six targeted galaxies is 0.9+/-0.6. No significant correlation is found for the L'_HCN(4-3)/L'_HCO^+(4-3) ratio with the SFR as traced by L_IR, nor with the warm-dust temperature, for the different populations of galaxies.

We analyze and model a C5.7 two-ribbon solar flare observed by SDO, Hinode and GOES on 2011 December 26. The flare is made of many loops formed and heated successively over one and half hours, and their footpoints are brightened in the UV 1600 A before enhanced soft X-ray and EUV missions are observed in flare loops. Assuming that anchored at each brightened UV pixel is a half flaring loop, we identify more than 6,700 half flaring loops, and infer the heating rate of each loop from the UV light curve at the foot-point. In each half loop, the heating rate consists of two phases, an intense impulsive heating followed by a low-rate heating persistent for more than 20 minutes. Using these heating rates, we simulate the evolution of their coronal temperatures and densities with the model of "enthalpy-based thermal evolution of loops" (EBTEL). In the model, suppression of thermal conduction is also considered. This model successfully reproduces total soft X-ray and EUV light curves observed in fifteen pass-bands by four instruments GOES, AIA, XRT, and EVE. In this flare, a total energy of 4.9x10^30 ergs is required to heat the corona, around 40% of this energy is in the slow heating phase. About two fifth of the total energy used to heat the corona is radiated by the coronal plasmas, and the other three fifth transported to the lower atmosphere by thermal conduction.

As one of the most popular applications of lidar systems, the atmospheric visibility is defined to be inversely proportional to the atmospheric extinction coefficient at 0.55 um. Since the laser at 1.5 um shows the highest maximum permissible exposure in the wavelength ranging from 0.3 um to 10 um, the eye-safe 1.5 um lidar can be deployed in urban areas. In such a case, the measured extinction coefficient at 1.5 um should be converted to that at 0.55 um for visibility retrieval. Although several models have been established since 1962, the accurate wavelength conversion remains a challenge. An adaptive inversion algorithm for 1.5 um visibility lidar is proposed and demonstrated by using the in situ Angstrom wavelength exponent, which is derived from either a sun photometer or an aerosol spectrometer. The impact of the particle size distribution of atmospheric aerosols and the Rayleigh backscattering of atmospheric molecules are taken into account. In comparison, the Angstrom wavelength exponent derived from the sun photometer is 7.7% higher than that derived from the aerosol spectrometer. Then, using 1.5 um visibility lidar, the visibility with a temporal resolution of 5 min is detected over 48 hours in Hefei (31.83 N, 117.25 E). The average visibility error between the new method and a visibility sensor (Vaisala, PWD52) is 5.2% with the R-square value of 0.96, while the relative error between another reference visibility lidar at 532 nm and the visibility sensor is 6.7% with the R-square value of 0.91. All results agree with each other well, demonstrating the accuracy and stability of the algorithm.

Seismic maps of the Sun's far hemisphere, computed from Doppler data from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) are now being used routinely to detect strong magnetic regions on the far side of the Sun (http:/jsoc.stanford.edu/data/farside/). To test the reliability of this technique, the helioseismically inferred active region detections are compared with far-side observation of solar activity from the Solar TErrestrial RElations Observatory (STEREO), using brightness in extreme ultraviolet light (EUV) as a proxy for magnetic fields. Two approaches are used to analyze nine months of STEREO and HMI data. In the first approach, we determine whether or not new large east-limb active regions are detected seismically on the far side before they appear Earth side and study how the detectability of these regions relates to their EUV intensity. We find that, while there is a range of EUV intensities for which far-side regions may or may not be detected seismically, there appears to be an intensity level above which they are almost always detected and an intensity level below which they are never detected. In the second approach, we analyze concurrent extreme ultraviolet and helioseismic far-side observations. We find that 100% (22) of the far-side seismic regions correspond to an extreme ultraviolet plage; 95% of these either became a NOAA-designated magnetic region when reaching the east limb or were one before crossing to the far side. A low but significant correlation is found between the seismic signature strength and the EUV intensity of a farside region.

We propose a mechanism for quasi-periodic oscillations of both coronal mass ejections (CMEs) and flare loops as related to magnetic reconnection in eruptive solar flares. We perform two-dimensional numerical MHD simulations of magnetic flux rope eruption, with three different values of the global Lundquist number. In the low Lundquist number run, no oscillatory behavior is found. In the moderate Lunquist number run, on the other hand, quasi-periodic oscillations are excited both at the bottom of the flux rope and at the flare loop-top. In the high Lundquist number run, quasi-periodic oscillations are also excited; in the meanwhile, the dynamics become turbulent due to the formation of multiple plasmoids in the reconnection current sheet. In high and moderate Lundquist number runs, thin reconnection jet collide with the flux rope bottom or flare loop-top and dig them deeply. Steep oblique shocks are formed as termination shocks where reconnection jet is bent (rather than decelerated) in horizontal direction, resulting in supersonic back-flows. The structure becomes unstable, and quasi-periodic oscillation of supersonic back-flows appear at locally confined high-beta region at both the flux rope bottom and flare loop-top. We compare the observational characteristics of quasi-periodic oscillations in erupting flux ropes, post-CME current sheets, flare ribbons and light curves, with corresponding dynamical structures found in our simulation.

We present observations of distinct UV spectral properties at different locations during an atypical X-shaped flare (SOL2014-11-09T15:32) observed by the Interface Region Imaging Spectrograph (IRIS). In this flare, four chromospheric ribbons appear and converge at an X-point where a separator is anchored. Above the X-point, two sets of non-coplanar coronal loops approach laterally and reconnect at the separator. The IRIS slit was located close to the X-point, cutting across some of the flare ribbons and loops. Near the location of the separator, the Si IV 1402.77 A line exhibits significantly broadened line wings extending to 200 km/s but an unshifted line core. These spectral features suggest the presence of bidirectional flows possibly related to the separator reconnection. While at the flare ribbons, the hot Fe XXI 1354.08 A line shows blueshifts and the cool Si IV 1402.77 A, C II 1335.71 A, and Mg II 2803.52 A lines show evident redshifts up to a velocity of 80 km/s, which are consistent with the scenario of chromospheric evaporation/condensation.

We have compiled a sample of 26 metal-poor galaxies with 12 + log(O/H) < 8.1 with both infrared continuum and 1.4 GHz radio continuum data. By comparing to galaxies at higher metallicity, we have investigated the dependence on the metallicity of the IR-radio relationship at 24 um, 70 um, 100 um and 160 um bands as well as the integrated FIR luminosity. It is found that metal-poor galaxies have on average lower qIR than metal-rich ones with larger offsets at longer IR wavelengths, from -0.06 dex in q24um to -0.6 dex in q160um. The qIR of all galaxies as a whole at 160 um show positive trends with the metallicity and IR-to-FUV ratio, and negative trends with the IR color, while those at lower IR wavelengths show weaker correlations. We proposed a mechanism that invokes combined effects of low obscured-SFR/total-SFR fraction and warm dust temperature at low metallicity to interpret the above behavior of qIR, with the former reducing the IR radiation and the latter further reducing the IR emission at longer IR wavelength. Other mechanisms that are related to the radio emission including the enhanced magnetic field strength and increased thermal radio contribution are unable to reconcile the IR-wavelength-dependent differences of qIR between metal-poor and metal- rich galaxies. In contrast to qIR, the mean total-SFR/radio ratio of metal-poor galaxies is the same as the metal-rich one, indicating the 1.4 GHz radio emission is still an effective tracer of SFRs at low metallicity.

We report observations of a two-stage coronal dimming in an eruptive event of a two-ribbon flare and a fast coronal mass ejection (CME). Weak gradual dimming persists for more than half an hour before the onset of the two-ribbon flare and the fast rise of the CME. It is followed by abrupt rapid dimming. The two-stage dimming occurs in a pair of conjugate dimming regions adjacent to the two flare ribbons, and the flare onset marks the transition between the two stages of dimming. At the onset of the two-ribbon flare, transient brightenings are also observed inside the dimming regions, before rapid dimming occurs at the same places. These observations suggest that the CME structure, most probably anchored at the twin dimming regions, undergoes a slow rise before the flare onset, and its kinematic evolution has significantly changed at the onset of flare reconnection. We explore diagnostics of the CME evolution in the early phase with analysis of the gradual dimming signatures prior to the CME eruption.

We present an analysis of the apparent elongation motion of flare ribbons along the polarity inversion line (PIL) as well as the shear of flare loops in several two-ribbon flares. Flare ribbons and loops spread along the PIL at a speed ranging from a few to a hundred km/s. The shear measured from conjugate foot-points is consistent with the measurement from flare loops, and both show the decrease of shear toward a potential field as a flare evolves and ribbons and loops spread along the PIL. Flares exhibiting fast bi-directional elongation appear to have a strong shear, which may indicate a large magnetic guide field relative to the reconnection field in the coronal current sheet. We discuss how the analysis of ribbon motion could help infer properties in the corona where reconnection takes place.

Solar flares are one of the most energetic events in the solar atmosphere. It is widely accepted that flares are powered by magnetic reconnection in the corona. An eruptive flare is usually accompanied by a coronal mass ejection, both of which are probably driven by the eruption of a magnetic flux rope (MFR). Here we report an eruptive flare on 2016 March 23 observed by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory. The extreme-ultraviolet imaging observations exhibit the clear rise and eruption of an MFR. In particular, the observations reveal solid evidence for magnetic reconnection from both the corona and chromosphere during the flare. Moreover, weak reconnection is observed before the start of the flare. We find that the preflare weak reconnection is of tether-cutting type and helps the MFR to rise slowly. Induced by a further rise of the MFR, strong reconnection occurs in the rise phases of the flare, which is temporally related to the MFR eruption. We also find that the magnetic reconnection is more of 3D-type in the early phase, as manifested in a strong-to-weak shear transition in flare loops, and becomes more 2D-like in the later phase, as shown by the apparent rising motion of an arcade of flare loops.

As the observational signature of the footprints of solar magnetic field lines open into the heliosphere, coronal holes provide a critical measure of the structure and evolution of these lines. Using a combination of Solar and Heliospheric Observatory / Extreme ultraviolet Imaging Telescope (SOHO/EIT), Solar Dynamics Observatory / Atmospheric Imaging Assembly (SDO/AIA) and Solar Terrestrial Relations Observatory / Extreme Ultraviolet Imager (STEREO/EUVI A/B) extreme ultraviolet (EUV) observations spanning 1996-2015 (nearly two solar cycles), coronal holes are automatically detected and characterized. Coronal hole area distributions show distinct behavior in latitude, defining the domain of polar and low-latitude coronal holes. The northern and southern polar regions show a clear asymmetry, with a lag between hemispheres in the appearance and disappearance of polar coronal holes.

A compact X-class flare on 2004-Feb-26 showed a concentrated source of hard X-rays at the tops of the flare's loops. This was analyzed in previous work (Longcope et al. 2010), and interpreted as plasma heated and compressed by slow magnetosonic shocks generated during post-reconnection retraction of the flux. That work used analytic expressions from a thin flux tube (TFT) model, which neglected many potentially important factors such as thermal conduction and chromospheric evaporation. Here we use a numerical solution of the TFT equations to produce a more comprehensive and accurate model of the same flare, including those effects previously omitted. These simulations corroborate the prior hypothesis that slow mode shocks persist well after the retraction has ended, thus producing a compact, loop-top source instead of an elongated jet, as steady reconnection models predict. Thermal conduction leads to densities higher than analytic estimates had predicted, and evaporation enhances the density still higher, but at lower temperatures. X-ray light curves and spectra are synthesized by convolving the results from a single TFT simulation with the rate at which flux is reconnected, as measured through motion of flare ribbons, for example. These agree well with light curves observed by RHESSI and GOES and spectra from RHESSI. An image created from a superposition of TFT model runs resembles one produced from RHESSI observations. This suggests that the HXR loop-top source, at least the one observed in this flare, could be the result of slow magnetosonic shocks produced in fast reconnection models like Petschek's.

An all-fiber, micro-pulse and eye-safe high spectral resolution wind lidar (HSRWL) at 1550nm is proposed and demonstrated by using a pair of upconversion single-photon detectors and a fiber Fabry-Perot scanning interferometer (FFP-SI). In order to improve the optical detection efficiency, both the transmission spectrum and the reflection spectrum of the FFP-SI are used for spectral analyses of the aerosol backscatter and the reference laser pulse. The reference signal is tapped from the outgoing laser and served as a zero velocity indicator. The Doppler shift is retrieved from a frequency response function Q, which is defined as the ratio of difference of the transmitted signal and the reflected signal to their sum. Taking advantages of high signal-to-noise ratio of the detectors and high spectral resolution of the FFP-SI, the Q spectra of the aerosol backscatter are reconstructed along the line-of-sight (LOS) of the telescope. By applying a least squares fit procedure to the measured Q spectra, the center frequencies and the bandwidths are obtained simultaneously. And then the Doppler shifts are determined relative to the center frequency of the reference signal. To eliminate the influence of temperature fluctuations on the FFP-SI, the FFP-SI is cased in a chamber with temperature stability of 0.001 during the measurement. Continuous LOS wind observations are carried out on two days at Hefei (31.843 N, 117.265 E), China. In the meantime, LOS wind measurements from the HSRWL show good agreement with the results from an ultrasonic wind sensor (Vaisala windcap WMT52). Due to the computational expensive of the convolution operation of the Q function, an empirical method is adopted to evaluate the quality of the measurements. The standard deviation of the wind speed is 0.76 m/s at the 1.8 km. The standard deviation of the retrieved bandwidth variation is 2.07 MHz at the 1.8 km.

For the first time, a versatile, eyesafe, compact and direct detection Doppler lidar is developed using upconversion single-photon detection method. An all-fiber and polarization maintaining architecture is realized to guarantee the high optical coupling efficiency and the system stability. Using integrated-optic components, the conservation of etendue of the optical receiver is achieved by manufacturing a fiber-coupled periodically poled Lithium niobate waveguide and an all-fiber Fabry-Perot interferometer (FPI). The so-called double-edge direct detection is implemented using a single-channel FPI and a single upconversion detector, incorporating time-division multiplexing method. The relative error of the system is lower than 0.1% over 9 weeks. To show the robust of the system, atmospheric wind and visibility over 48 hours are detected in the boundary layer. In the intercomparison experiments, lidar shows good agreement with the ultrasonic wind sensor (Vaisala windcap WMT52), with standard deviation of 1.04 m/s in speed and 12.3{\deg} in direction.

We report evolution of an atypical X-shaped flare ribbon which provides novel observational evidence of three-dimensional (3D) magnetic reconnection at a separator. The flare occurred on 2014 November 9. High-resolution slit-jaw 1330 A images from the Interface Region Imaging Spectrograph reveal four chromospheric flare ribbons that converge and form an X-shape. Flare brightening in the upper chromosphere spreads along the ribbons toward the center of the "X" (the X-point), and then spreads outward in a direction more perpendicular to the ribbons. These four ribbons are located in a quadrupolar magnetic field. Reconstruction of magnetic topology in the active region suggests the presence of a separator connecting to the X-point outlined by the ribbons. The inward motion of flare ribbons in the early stage therefore indicates 3D magnetic reconnection between two sets of non-coplanar loops that approach laterally, and reconnection proceeds downward along a section of vertical current sheet. Coronal loops are also observed by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory confirming the reconnection morphology illustrated by ribbon evolution.

This paper reports results of the third-year campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs) between 2014-2015. Ten new targets were selected from quasar sample of Sloan Digital Sky Survey (SDSS), which are generally more luminous than the SEAMBH candidates in last two years. H$\beta$ lags ($\tau_{_{\rm H\beta}}$) in five of the 10 quasars have been successfully measured in this monitoring season. We find that the lags are generally shorter, by large factors, than those of objects with same optical luminosity, in light of the well-known $R_{_{\rm H\beta}}-L_{5100}$ relation. The five quasars have dimensionless accretion rates of $\dot{\mathscr{M}}=10-10^3$. Combining measurements of the previous SEAMBHs, we find that the reduction of H$\beta$ lags tightly depends on accretion rates, $\tau_{_{\rm H\beta}}/\tau_{_{R-L}}\propto\dot{\mathscr{M}}^{-0.42}$, where $\tau_{_{R-L}}$ is the H$\beta$ lag from the normal $R_{_{\rm H\beta}}-L_{5100}$ relation. Fitting 63 mapped AGNs, we present a new scaling relation for the broad-line region: $R_{_{\rm H\beta}}=\alpha_1\ell_{44}^{\beta_1}\,\min\left[1,\left(\dot{\mathscr{M}}/\dot{\mathscr{M}}_c\right)^{-\gamma_1}\right]$, where $\ell_{44}=L_{5100}/10^{44}\,\rm erg~s^{-1}$ is 5100 \AA\ continuum luminosity, and coefficients of $\alpha_1=(29.6_{-2.8}^{+2.7})$ lt-d, $\beta_1=0.56_{-0.03}^{+0.03}$, $\gamma_1=0.52_{-0.16}^{+0.33}$ and $\dot{\mathscr{M}}_c=11.19_{-6.22}^{+2.29}$. This relation is applicable to AGNs over a wide range of accretion rates, from $10^{-3}$ to $10^3$. Implications of this new relation are briefly discussed.

Coronal mass ejections (CMEs) are often accompanied by coronal dimming evident in extreme ultraviolet (EUV) and soft X-ray observations. The locations of dimming are sometimes considered to map footpoints of the erupting flux rope. As the emitting material expands in the corona, the decreased plasma density leads to reduced emission observed in spectral and irradiance measurements. Therefore, signatures of dimming may reflect properties of CMEs in the early phase of its eruption. In this study, we analyze the event of flare, CME, and coronal dimming on December 26, 2011. We use the data from the Atmospheric Imaging Assembly (AIA) on Solar Dynamics Observatories (SDO) for disk observations of the dimming, and analyze images taken by EUVI, COR1, and COR2 onboard the Solar Terrestrial Relations Observatories to obtain the height and velocity of the associated CMEs observed at the limb. We also measure magnetic reconnection rate from flare observations. Dimming occurs in a few locations next to the flare ribbons, and it is observed in multiple EUV passbands. Rapid dimming starts after the onset of fast reconnection and CME acceleration, and its evolution well tracks the CME height and flare reconnection. The spatial distribution of dimming exhibits cores of deep dimming with a rapid growth, and their light curves are approximately linearly scaled with the CME height profile. From the dimming analysis, we infer the process of the CME expansion, and estimate properties of the CME.

Flare emissions in X-ray and EUV wavelengths have previously been modeled as the plasma response to impulsive heating from magnetic reconnection. Some flares exhibit gradually evolving X-ray and EUV light curves, which are believed to result from superposition of an extended sequence of impulsive heating events occurring in different adjacent loops or even unresolved threads within each loop. In this paper, we apply this approach to a long duration two-ribbon flare SOL2011-09-13T22 observed by the Atmosphere Imaging Assembly (AIA). We find that to reconcile with observed signatures of flare emission in multiple EUV wavelengths, each thread should be heated in two phases, an intense impulsive heating followed by a gradual, low-rate heating tail that is attenuated over 20-30 minutes. Each AIA resolved single loop may be composed of several such threads. The two-phase heating scenario is supported by modeling with both a zero-dimensional and a 1D hydrodynamic code. We discuss viable physical mechanisms for the two-phase heating in a post-reconnection thread.

As a natural consequence of cosmological hierarchical structure formation, sub-parsec supermassive black hole binaries (SMBHBs) should be common in galaxies but thus far have eluded spectroscopic identification. Based on four decades of optical spectroscopic monitoring, we report that the nucleus of NGC 5548, a nearby Seyfert galaxy long suspected to have experienced a major merger about one billion years ago, exhibits long-term variability with a period of 14 years in the optical continuum and broad Hbeta emission line. Remarkably, the double-peaked profile of Hbeta shows systematic velocity changes with a similar period. These pieces of observations plausibly indicate that a SMBHB resides in the center of NGC 5548. The complex, secular variations in the line profiles can be explained by orbital motion of a binary with equal mass and a semi-major axis of 22 light-days (corresponding to 18 milli-parsec). At a distance of 75 Mpc, NGC 5548 is one of the nearest sub-parsec SMBHB candidates that offers an ideal laboratory for gravitational wave detection.

In the sixth of the series of papers reporting on a large reverberation mapping (RM) campaign of active galactic nuclei (AGNs) with high accretion rates, we present velocity-resolved time lags of H$\beta$ emission lines for nine objects observed in the campaign during 2012$-$2013. In order to correct the line-broadening caused by seeing and instruments before the analysis of velocity-resolved RM, we adopt Richardson-Lucy deconvolution to reconstruct their H$\beta$ profiles. The validity and effectiveness of the deconvolution are checked out by Monte Carlo simulation. Five among the nine objects show clear dependence of time delay on velocity. Mrk 335 and Mrk 486 show signatures of gas inflow whereas the clouds in the broad-line regions (BLRs) of Mrk 142 and MCG +06-26-012 tend to be radial outflowing. Mrk 1044 is consistent with the case of virialized motions. The lags of the rest four are not velocity-resolvable. The velocity-resolved RM of super-Eddington accreting massive black holes (SEAMBHs) shows that they have diversity of the kinematics in their BLRs. Comparing with the AGNs with sub-Eddington accretion rates, we do not find significant differences in the BLR kinematics of SEAMBHs.

Chromospheric evaporation refers to dynamic mass motions in flare loops as a result of rapid energy deposition in the chromosphere. These have been observed as blueshifts in X-ray and extreme-ultraviolet (EUV) spectral lines corresponding to upward motions at a few tens to a few hundreds of km/s. Past spectroscopic observations have also revealed a dominant stationary component, in addition to the blueshifted component, in emission lines formed at high temperatures (~10 MK). This is contradictory to evaporation models predicting predominant blueshifts in hot lines. The recently launched Interface Region Imaging Spectrograph (IRIS) provides high resolution imaging and spectroscopic observations that focus on the chromosphere and transition region in the UV passband. Using the new IRIS observations, combined with coordinated observations from the EUV Imaging Spectrometer, we study the chromospheric evaporation process from the upper chromosphere to corona during an X1.0 flare on 2014 March 29. We find evident evaporation signatures, characterized by Doppler shifts and line broadening, at two flare ribbons separating from each other, suggesting that chromospheric evaporation takes place in successively formed flaring loops throughout the flare. More importantly, we detect dominant blueshifts in the high temperature Fe XXI line (~10 MK), in agreement with theoretical predictions. We also find that, in this flare, gentle evaporation occurs at some locations in the rise phase of the flare, while explosive evaporation is detected at some other locations near the peak of the flare. There is a conversion from gentle to explosive evaporation as the flare evolves.

The bulk of the radiative output of a solar flare is emitted from the chromosphere, which produces enhancements in the optical and UV continuum, and in many lines, both optically thick and thin. We have, until very recently, lacked observations of two of the strongest of these lines: the Mg II h & k resonance lines. We present a detailed study of the response of these lines to a solar flare. The spatial and temporal behaviour of the integrated intensities, k/h line ratios, line of sight velocities, line widths and line asymmetries were investigated during an M class flare (SOL2014-02-13T01:40). Very intense, spatially localised energy input at the outer edge of the ribbon is observed, resulting in redshifts equivalent to velocities of ~15-26km/s, line broadenings, and a blue asymmetry in the most intense sources. The characteristic central reversal feature that is ubiquitous in quiet Sun observations is absent in flaring profiles, indicating that the source function increases with height during the flare. Despite the absence of the central reversal feature, the k/h line ratio indicates that the lines remain optically thick during the flare. Subordinate lines in the Mg II passband are observed to be in emission in flaring sources, brightening and cooling with similar timescales to the resonance lines. This work represents a first analysis of potential diagnostic information of the flaring atmosphere using these lines, and provides observations to which synthetic spectra from advanced radiative transfer codes can be compared.

We present imaging and spectroscopic observations from the Interface Region Imaging Spectrograph (IRIS) of the evolution of the flare ribbon in the SOL2014-04-18T13:03 M-class flare event, at high spatial resolution and time cadence. These observations reveal small-scale substructure within the ribbon, which manifests as coherent quasi-periodic oscillations in both position and Doppler velocities. We consider various alternative explanations for these oscillations, including modulation of chromospheric evaporation flows. Among these we find the best support for some form of wave localized to the coronal current sheet, such as a tearing mode or Kelvin-Helmholtz instability.

We have completed two years of photometric and spectroscopic monitoring of a large number of active galactic nuclei (AGNs) with very high accretion rates. In this paper, we report on the result of the second phase of the campaign, during 2013--2014, and the measurements of five new H$\beta$ time lags out of eight monitored AGNs. All five objects were identified as super-Eddington accreting massive black holes (SEAMBHs). The highest measured accretion rates for the objects in this campaign are $\dot{\mathscr{M}}\gtrsim 200$, where $\dot{\mathscr{M}}= \dot{M}_{\bullet}/L_{\rm Edd}c^{-2}$, $\dot{M}_{\bullet}$ is the mass accretion rates, $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. We find that the H$\beta$ time lags in SEAMBHs are significantly shorter than those measured in sub-Eddington AGNs, and the deviations increase with increasing accretion rates. Thus, the relationship between broad-line region size ($R_{_{\rm H\beta}}$) and optical luminosity at 5100\AA, $R_{_{\rm H\beta}}-L_{5100}$, requires accretion rate as an additional parameter. We propose that much of the effect may be due to the strong anisotropy of the emitted slim-disk radiation. Scaling $R_{_{\rm H\beta}}$ by the gravitational radius of the black hole, we define a new radius-mass parameter ($Y$) and show that it saturates at a critical accretion rate of $\dot{\mathscr{M}}_c=6\sim 30$, indicating a transition from thin to slim accretion disk and a saturated luminosity of the slim disks. The parameter $Y$ is a very useful probe for understanding the various types of accretion onto massive black holes. We briefly comment on implications to the general population of super-Eddington AGNs in the universe and applications to cosmology.

This is the third in a series of papers reporting on a large reverberation-mapping campaign aimed to study the properties of active galactic nuclei (AGNs) with high accretion rates. We present new results on the variability of the optical Fe II emission lines in 10 AGNs observed by the Yunnan Observatory 2.4m telescope during 2012--2013. We detect statistically significant time lags, relative to the AGN continuum, in nine of the sources. This accurate measurement is achieved by using a sophisticated spectral fitting scheme that allows for apparent flux variations of the host galaxy, and several narrow lines, due to the changing observing conditions. Six of the newly detected lags are indistinguishable from the Hbeta lags measured in the same sources. Two are significantly longer and one is slightly shorter. Combining with Fe II lags reported in previous studies, we find a Fe II radius--luminosity relationship similar to the one for Hbeta, although our sample by itself shows no clear correlation. The results support the idea that Fe II emission lines originate in photoionized gas which, for the majority of the newly reported objects, is indistinguishable from the Hbeta-emitting gas. We also present a tentative correlation between the lag and intensity of Fe II and Hbeta and comment on its possible origin.

Warped accretion disks have attracted intensive attention because of their critical role on shaping the spin of supermassive massive black holes (SMBHs) through the Bardeen-Petterson effect, a general relativistic effect that leads to final alignments or anti-alignments between black holes and warped accretion disks. We study such alignment processes by explicitly taking into account the finite sizes of accretion disks and the episodic lifetimes of AGNs that delineate the duration of gas fueling onto accretion disks. We employ an approximate global model to simulate the evolution of accretion disks, allowing to determine the gravitomagnetic torque that drives the alignments in a quite simple way. We then track down the evolutionary paths for mass and spin of black holes both in a single activity episode and over a series of episodes. Given with randomly and isotropically oriented gas fueling over episodes, we calculate the spin evolution with different episodic lifetimes and find that it is quite sensitive to the lifetimes. We therefore propose that spin distribution of SMBHs can place constraints on the episodic lifetimes of AGNs and vice versa. Applications of our results on the observed spin distributions of SMBHs and the observed episodic lifetimes of AGNs are discussed, although both the measurements at present are yet ambiguous to draw a firm conclusion. Our prescription can be easily incorporated into semi-analytic models for black hole growth and spin evolution.

Coronal holes are regions on the Sun's surface that map the foot-prints of open magnetic field lines. We have developed an automated routine to detect and track boundaries of long-lived coronal holes using full-disk EUV images obtained by SoHO:EIT, SDO:AIA, and STEREO:EUVI. We measure coronal hole areas and magnetic flux in these holes, and compare the measurements with calculations by the PFSS model. It is shown that, from 1996 through 2010, the total area of coronal holes measured with EIT images varies between 5$\%$ and 17$\%$ of the total solar surface area, and the total unsigned open flux varies between $2-5 \times 10^{22}$ Mx. The solar cycle dependence of these measurements are similar to the PFSS results, but the model yields larger hole areas and greater open flux than observed by EIT. The AIA/EUVI measurements from 2010-2013 show coronal hole area coverage of 5-10$\%$ of the total surface area, with significant contribution from low latitudes, which is under-represented by EIT. AIA/EUVI have measured much enhanced open magnetic flux in the range of $2-4 \times 10^{22}$ Mx, which is about twice the flux measured by EIT, and matches with the PFSS calculated open flux, with discrepancies in the location and strength of coronal holes. A detailed comparison between the three measurements (by EIT, AIA-EUVI, and PFSS) indicates that coronal holes in low latitudes contribute significantly to the total open magnetic flux. These low-latitude coronal holes are not well measured with either the He I 10830 line in previous studies, or EIT EUV images; neither are they well captured by the static PFSS model. The enhanced observations from AIA/EUVI allow a more accurate measure of these low latitude coronal holes, and their contribution to open magnetic flux.

We report on the detailed and systematic study of field-line twist and length distributions within magnetic flux ropes embedded in Interplanetary Coronal Mass Ejections (ICMEs). The Grad-Shafranov reconstruction method is utilized together with a constant-twist nonlinear force-free (Gold-Hoyle) flux rope model to reveal the close relation between the field-line twist and length in cylindrical flux ropes, based on in-situ Wind spacecraft measurements. We show that the field-line twist distributions within interplanetary flux ropes are inconsistent with the Lundquist model. In particular we utilize the unique measurements of magnetic field-line lengths within selected ICME events as provided by Kahler et al. (2011) based on energetic electron burst observations at 1 AU and the associated type III radio emissions detected by the Wind spacecraft. These direct measurements are compared with our model calculations to help assess the flux-rope interpretation of the embedded magnetic structures. By using the different flux-rope models, we show that the in-situ direct measurements of field-line lengths are consistent with a flux-rope structure with spiral field lines of constant and low twist, largely different from that of the Lundquist model, especially for relatively large-scale flux ropes.

Supermassive black holes in active galactic nuclei (AGNs) undergo a wide range of accretion rates, which lead to diversity of appearance. We consider the effects of anisotropic radiation from accretion disks on the broad-line region (BLR), from the Shakura-Sunyaev regime to slim disks with super-Eddington accretion rates. The geometrically thick funnel of the inner region of slim disks produces strong self-shadowing effects that lead to very strong anisotropy of the radiation field. We demonstrate that the degree of anisotropy of the radiation fields grows with increasing accretion rate. As a result of this anisotropy, BLR clouds receive different spectral energy distributions depending on their location relative to the disk, resulting in diverse observational appearance of the BLR. We show that the self-shadowing of the inner parts of the disk naturally produces two dynamically distinct regions of the BLR, depending on accretion rate. These two regions manifest themselves as kinematically distinct components of the broad H$\beta$ line profile with different line widths and fluxes, which jointly account for the Lorentzian profile generally observed in narrow-line Seyfert 1 galaxies. In the time domain, these two components are expected reverberate with different time lags with respect to the varying ionizing continuum, depending on the accretion rate and the viewing angle of the observer. The diverse appearance of the BLR due to the anisotropic ionizing energy source can be tested by reverberation mapping of H$\beta$ and other broad emission lines (e.g., \feii), providing a new tool to diagnose the structure and dynamics of the BLR. Other observational consequences of our model are also explored.

This is the second in a series of papers reporting on a large reverberation mapping (RM) campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). The goal is to identify super-Eddington accreting massive black holes (SEAMBHs) and to use their unique properties to construct a new method for measuring cosmological distances. Based on theoretical models, the saturated bolometric luminosity of such sources is proportional to the BH mass which can be used to obtain their distance. Here we report on five new RM measurements and show that in four of the cases we can measure the BH mass and three of these sources are SEAMBHs. Together with the three sources from our earlier work, we now have six new sources of this type. We use a novel method based on a minimal radiation efficiency to identify nine additional SEAMBHs from earlier RM-based mass measurements. We use a Bayesian analysis to determine the parameters of the new distance expression, and the method uncertainties, from the observed properties of the objects in the sample. The ratio of the newly measured distances to the standard cosmological ones has a mean scatter of 0.14 dex, indicating that SEAMBHs can be use as cosmological distance probes. With their high luminosity, long period of activity and large numbers at high redshifts, SEAMBHs have a potential to extend the cosmic distance ladder beyond the range now explored by type Ia supernovae.

Magnetic reconnection is essential to release the flux rope during its ejection. The question remains: how does the magnetic reconnection change the flux rope structure? Following the original study of \citet{Qiu2007}, we compare properties of ICME/MC flux ropes measured at 1 AU and properties of associated solar progenitors including flares, filaments, and CMEs. In particular, the magnetic field-line twist distribution within interplanetary magnetic flux ropes is systematically derived and examined. Our analysis shows that for most of these events, the amount of twisted flux per AU in MCs is comparable with the total reconnection flux on the Sun, and the sign of the MC helicity is consistent with the sign of helicity of the solar source region judged from the geometry of post-flare loops. Remarkably, we find that about one half of the 18 magnetic flux ropes, most of them being associated with erupting filaments, have a nearly uniform and relatively low twist distribution from the axis to the edge, and the majority of the other flux ropes exhibit very high twist near the axis, of up to $\gtrsim 5$ turns per AU, which decreases toward the edge. The flux ropes are therefore not linear force free. We also conduct detailed case studies showing the contrast of two events with distinct twist distribution in MCs as well as different flare and dimming characteristics in solar source regions, and discuss how reconnection geometry reflected in flare morphology may be related to the structure of the flux rope formed on the Sun.

On June 17, 2012, we observed the evolution of a C-class flare associated with the eruption of a filament near a large sunspot in the active region NOAA 11504. We obtained high spatial resolution filtergrams using the 1.6 m New Solar Telescope at the Big Bear Solar Observatory in broad-band TiO at 706 nm (bandpass:10 \AA) and He I 10830 \AA\ narrow-band (bandpass: 0.5 \AA, centered 0.25 \AA\ to the blue). We analyze the spatio-temporal behavior of the He I 10830 \AA\ data, which were obtained over a 90" X 90" field of view with a cadence of 10 sec. We also analyze simultaneous data from the Atmospheric Imaging Assembly and Extreme Ultraviolet Variability Experiment instruments on board the Solar Dynamics Observatory spacecraft, and data from Reuven Ramaty High Energy Solar Spectroscopic Imager and GOES spacecrafts. Non-thermal effects are ignored in this analysis. Several quantitative aspects of the data, as well as models derived using the "0D" Enthalpy-Based Thermal Evolution of Loops model (EBTEL: Klimchuk et al. 2008) code, indicate that the triplet states of the 10830 \AA\ multiplet are populated by photoionization of chromospheric plasma followed by radiative recombination. Surprisingly, the He II 304 \AA\ line is reasonably well matched by standard emission measure calculations, along with the C IV emission which dominates the AIA 1600 \AA\ channel during flares. This work lends support to some of our previous work combining X-ray, EUV and UV data of flares to build models of energy transport from corona to chromosphere.

We highlight the progress, current status, and open challenges of QCD-driven physics, in theory and in experiment. We discuss how the strong interaction is intimately connected to a broad sweep of physical problems, in settings ranging from astrophysics and cosmology to strongly-coupled, complex systems in particle and condensed-matter physics, as well as to searches for physics beyond the Standard Model. We also discuss how success in describing the strong interaction impacts other fields, and, in turn, how such subjects can impact studies of the strong interaction. In the course of the work we offer a perspective on the many research streams which flow into and out of QCD, as well as a vision for future developments.

We investigate heating and evolution of flare loops in a C4.7 two-ribbon flare on 2011 February 13. From SDO/AIA imaging observations, we can identify two sets of loops. Hinode/EIS spectroscopic observations reveal blueshifts at the feet of both sets of loops. The evolution and dynamics of the two sets are quite different. The first set of loops exhibits blueshifts for about 25 minutes followed by redshifts, while the second set shows stronger blueshifts, which are maintained for about one hour. The UV 1600 observation by AIA also shows that the feet of the second set of loops brighten twice. These suggest that continuous heating may be present in the second set of loops. We use spatially resolved UV light curves to infer heating rates in the few tens of individual loops comprising the two loop systems. With these heating rates, we then compute plasma evolution in these loops with the "enthalpy-based thermal evolution of loops" (EBTEL) model. The results show that, for the first set of loops, the synthetic EUV light curves from the model compare favorably with the observed light curves in six AIA channels and eight EIS spectral lines, and the computed mean enthalpy flow velocities also agree with the Doppler shift measurements by EIS. For the second set of loops modeled with twice-heating, there are some discrepancies between modeled and observed EUV light curves in low-temperature bands, and the model does not fully produce the prolonged blueshift signatures as observed. We discuss possible causes for the discrepancies.

We report first results from a large project to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). Such objects may be different from other AGNs in being powered by slim accretion disks and showing saturated accretion luminosities, but both are not yet fully understood. The results are part of a large reverberation mapping (RM) campaign using the 2.4-m Shangri-La telescope at the Yunnan Observatory in China. The goals are to investigate the gas distribution near the BH and the properties of the central accretion disks, to measure BH mass and Eddington ratios, and to test the feasibility of using such objects as a new type of cosmological candles. The paper presents results for three objects, Mrk 335, Mrk 142 and IRAS F12397+3333 with H$\beta$ time lags relative to the 5100\AA\ continuum of $10.6^{+1.7}_{-2.9}$, $6.4^{+0.8}_{-2.2}$ and $11.4^{+2.9}_{-1.9}$ days, respectively. The corresponding BH masses are $(8.3_{-3.2}^{+2.6})\times 10^6M_{\odot}$, $(3.4_{-1.2}^{+0.5})\times 10^6M_{\odot}$ and $(7.5_{-4.1}^{+4.3})\times 10^6M_{\odot}$, and the lower limits on the Eddington ratios 0.6, 2.3, and 4.6 for the minimal radiative efficiency of 0.038. Mrk 142 and IRAS F12397+333 (extinction corrected) clearly deviate from the currently known relation between H$\beta$ lag and continuum luminosity. The three Eddington ratios are beyond the values expected in thin accretion disks and two of them are the largest measured so far among objects with RM-based BH masses. We briefly discuss implications for slim disks, BH growth and cosmology.

A solar flare is composed of impulsive energy release events by magnetic reconnection, which forms and heats flare loops. Recent studies have revealed a two-phase evolution pattern of UV 1600\AA\ emission at the feet of these loops: a rapid pulse lasting for a few seconds to a few minutes, followed by a gradual decay on timescales of a few tens of minutes. Multiple band EUV observations by AIA further reveal very similar signatures. These two phases represent different but related signatures of an impulsive energy release in the corona. The rapid pulse is an immediate response of the lower atmosphere to an intense thermal conduction flux resulting from the sudden heating of the corona to high temperatures (we rule out energetic particles due to a lack of significant hard X-ray emission). The gradual phase is associated with the cooling of hot plasma that has been evaporated into the corona. The observed footpoint emission is again powered by thermal conduction (and enthalpy), but now during a period when approximate steady state conditions are established in the loop. UV and EUV light curves of individual pixels may therefore be separated into contributions from two distinct physical mechanisms to shed light on the nature of energy transport in a flare. We demonstrate this technique using coordinated, spatially resolved observations of UV and EUV emission from the footpoints of a C3.2 thermal flare.

We analyze and model an M8.0 flare on 2005 May 13 observed by TRACE and RHESSI to determine the energy release rate from magnetic reconnection that forms and heats numerous flare loops. The flare exhibits two ribbons in UV 1600 {\AA} emission. Analysis shows that the UV light curve at each flaring pixel rises impulsively within a few minutes, and decays slowly with a timescale >10 min. Since the lower atmosphere (transition region and chromosphere) responds to energy deposit nearly instantaneously, the rapid UV brightening is thought to reflect the energy release process in the newly formed flare loop rooted at the footpoint. We utilize spatially resolved (down to 1 arcsec) UV light curves and thick-target hard X-ray emission to construct heating functions of a few thousand flare loops anchored at UV foot points, and compute plasma evolution in these loops using the EBTEL model. The modeled coronal temperatures and densities of these flare loops are used to calculate synthetic soft X-ray spectra and light curves, which compare favorably with those observed by RHESSI and GOES/XRS. The time-dependent transition region DEM for each loop during its decay phase is also computed with a simplified model and used to calculate the optically-thin C IV line emission, which dominates the UV 1600 {\AA} bandpass during the flare. The computed C IV line emission decays at the same rate as observed. This study presents a method to constrain heating of reconnection-formed flare loops using all available observables independently, and provides insight into the physics of energy release and plasma heating during the flare. With this method, the lower limit of the total energy used to heat the flare loops in this event is estimated to be 1.22e31 ergs, of which only 1.9e30 ergs is carried by beam-driven upflows during the impulsive phase, suggesting that the coronal plasmas are predominantly heated in situ.

We investigate the alignment processes of spinning black holes and their surrounding warped accretion disks in a frame of two different types of feeding at the outer boundaries. We consider (1) fixed flows in which gas is continually fed with a preferred angular momentum, and (2) free flows in which there is no gas supply and the disks diffuse freely at their outer edges. As expected, we find that for the cases of fixed flows the black hole disk systems always end up aligning on timescales of several 1e6 yr, irrespective of the initial inclinations. If the initial inclination angles are larger than pi/2, the black hole accretion transits from retrograde to prograde fashion, and the accreted mass onto the black holes during these two phases is comparable. On the other hand, for the cases of free flows, both alignments and anti-alignments can occur, depending on the initial inclinations and the ratios of the angular momentum of the disks to that of the black holes. In such cases, the disks will be consumed within timescales of 1e6 yr by black holes accreting at the Eddington limit. We propose that there is a close connection between the black hole spin and the lifetime for which the feeding persists, which determines the observable episodic lifetimes of active galactic nuclei. We conclude that careful inclusion of the disk feeding at the outer boundaries is crucial for modeling the evolution of the black hole spin.

In this paper, we analyze hard X-ray spikes observed by RHESSI to understand their temporal, spectral, and spatial properties. A recently developed demodulation code was applied to hard X-ray light curves in several energy bands observed by RHESSI. Hard X-ray spikes were selected from the demodulated flare light curves. We measured the spike duration, the energy-dependent time delay, and count spectral index of these spikes. We also located the hard X-ray source emitting these spikes from RHESSI mapping that was coordinated with imaging observations in visible and UV wavelengths. We identify quickly varying structures of <1 s during the rise of hard X-rays in five flares. These hard X-ray spikes can be observed at photon energies over 100 keV. They exhibit sharp rise and decay with a duration (FWHM) of less than 1 s. Energy-dependent time lags are present in some spikes. It is seen that the spikes exhibit harder spectra than underlying components, typically by 0.5 in the spectral index when they are fitted to power-law distributions. RHESSI clean maps at 25-100 keV with an integration of 2 s centered on the peak of the spikes suggest that hard X-ray spikes are primarily emitted by double foot-point sources in magnetic fields of opposite polarities. With the RHESSI mapping resolution of ~ 4 arsec, the hard X-ray spike maps do not exhibit detectable difference in the spatial structure from sources emitting underlying components. Coordinated high-resolution imaging UV and infrared observations confirm that hard X-ray spikes are produced in magnetic structures embedded in the same magnetic environment of the underlying components. The coordinated high-cadence TRACE UV observations of one event possibly reveal new structures on spatial scales <1-2 arsec at the time of the spike superposed on the underlying component. They are probably sources of hard X-ray spikes.

Context. Hard X-ray (HXR) spikes refer to fine time structures on timescales of seconds to milliseconds in high-energy HXR emission profiles during solar flare eruptions. Aims. We present a preliminary statistical investigation of temporal and spectral properties of HXR spikes. Methods. Using a three-sigma spike selection rule, we detected 184 spikes in 94 out of 322 flares with significant counts at given photon energies, which were detected from demodulated HXR light curves obtained by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). About one fifth of these spikes are also detected at photon energies higher than 100 keV. Results. The statistical properties of the spikes are as follows. (1) HXR spikes are produced in both impulsive flares and long-duration flares with nearly the same occurrence rates. Ninety percent of the spikes occur during the rise phase of the flares, and about 70% occur around the peak times of the flares. (2) The time durations of the spikes vary from 0.2 to 2 s, with the mean being 1.0 s, which is not dependent on photon energies. The spikes exhibit symmetric time profiles with no significant difference between rise and decay times. (3) Among the most energetic spikes, nearly all of them have harder count spectra than their underlying slow-varying components. There is also a weak indication that spikes exhibiting time lags in high-energy emissions tend to have harder spectra than spikes with time lags in low-energy emissions.

We analyze and model an M1.0 flare observed by SDO/AIA and Hinode/EIS to investigate how flare loops are heated and evolve subsequently. The flare is composed of two distinctive loop systems observed in EUV images. The UV 1600 \AA emission at the feet of these loops exhibits a rapid rise, followed by enhanced emission in different EUV channels observed by AIA and EIS. Such behavior is indicative of impulsive energy deposit and the subsequent response in overlying coronal loops that evolve through different temperatures. Using the method we recently developed, we infer empirical heating functions from the rapid rise of the UV light curves for the two loop systems, respectively, treating them as two big loops of cross-sectional area 5\arcsec by 5\arcsec, and compute the plasma evolution in the loops using the EBTEL model (Klimchuk et al. 2008). We compute the synthetic EUV light curves, which, with the limitation of the model, reasonably agree with observed light curves obtained in multiple AIA channels and EIS lines: they show the same evolution trend and their magnitudes are comparable by within a factor of two. Furthermore, we also compare the computed mean enthalpy flow velocity with the Doppler shift measurements by EIS during the decay phase of the two loops. Our results suggest that the two different loops with different heating functions as inferred from their footpoint UV emission, combined with their different lengths as measured from imaging observations, give rise to different coronal plasma evolution patterns captured both in the model and observations.

We analyze high cadence high resolution observations of a C3.2 flare obtained by AIA/SDO on August 1, 2010. The flare is a long duration event with soft X-ray and EUV radiation lasting for over four hours. Analysis suggests that magnetic reconnection and formation of new loops continue for more than two hours. Furthermore, the UV 1600\AA\ observations show that each of the individual pixels at the feet of flare loops is brightened instantaneously with a timescale of a few minutes, and decays over a much longer timescale of more than 30 minutes. We use these spatially resolved UV light curves during the rise phase to construct empirical heating functions for individual flare loops, and model heating of coronal plasmas in these loops. The total coronal radiation of these flare loops are compared with soft X-ray and EUV radiation fluxes measured by GOES and AIA. This study presents a method to observationally infer heating functions in numerous flare loops that are formed and heated sequentially by reconnection throughout the flare, and provides a very useful constraint to coronal heating models.

In this paper, we present comprehensive analysis of a two-ribbon flare observed in UV 1600{\AA} by Transition Region and Coronal Explorer and in HXRs by Reuven Ramaty High Energy Solar Spectroscopic Imager. HXR (25-100 keV) imaging observations show two kernels of size (FWHM) 15?? moving along the two UV ribbons. We find the following results. (1) UV brightening is substantially enhanced wherever and whenever the compact HXR kernel is passing, and during the HXR transit across a certain region, the UV count light curve in that region is temporally correlated with the HXR total flux light curve. After the passage of the HXR kernel, the UV light curve exhibits smooth monotonical decay. (2)We measure the apparent motion speed of the HXR sources and UV ribbon fronts, and decompose the motion into parallel and perpendicular motions with respect to the magnetic polarity inversion line (PIL). It is found that HXR kernels and UV fronts exhibit similar apparent motion patterns and speeds. The parallel motion dominates during the rise of the HXR emission, and the perpendicular motion starts and dominates at the HXR peak, the apparent motion speed being 10-40 km s-1. (3) We also find that UV emission is characterized by a rapid rise correlated with HXRs, followed by a long decay on timescales of 15-30 minutes. The above analysis provides evidence that UV brightening is primarily caused by beam heating, which also produces thick-target HXR emission. The thermal origin of UV emission cannot be excluded, but would produce weaker heating by one order of magnitude. The extended UV ribbons in this event are most likely a result of sequential reconnection along the PIL, which produces individual flux tubes (post-flare loops), subsequent non-thermal energy release and heating in these flux tubes, and then the very long cooling time of the transition region at the feet of these flux tubes.

We present a quantitative model of the magnetic energy stored and then released through magnetic reconnection for a flare on 26 Feb 2004. This flare, well observed by RHESSI and TRACE, shows evidence of non-thermal electrons only for a brief, early phase. Throughout the main period of energy release there is a super-hot (T>30 MK) plasma emitting thermal bremsstrahlung atop the flare loops. Our model describes the heating and compression of such a source by localized, transient magnetic reconnection. It is a three-dimensional generalization of the Petschek model whereby Alfven-speed retraction following reconnection drives supersonic inflows parallel to the field lines, which form shocks heating, compressing, and confining a loop-top plasma plug. The confining inflows provide longer life than a freely-expanding or conductively-cooling plasma of similar size and temperature. Superposition of successive transient episodes of localized reconnection across a current sheet produces an apparently persistent, localized source of high-temperature emission. The temperature of the source decreases smoothly on a time scale consistent with observations, far longer than the cooling time of a single plug. Built from a disordered collection of small plugs, the source need not have the coherent jet-like structure predicted by steady-state reconnection models. This new model predicts temperatures and emission measure consistent with the observations of 26 Feb 2004. Furthermore, the total energy released by the flare is found to be roughly consistent with that predicted by the model. Only a small fraction of the energy released appears in the super-hot source at any one time, but roughly a quarter of the flare energy is thermalized by the reconnection shocks over the course of the flare. All energy is presumed to ultimately appear in the lower-temperature T<20 MK, post-flare loops.

n order to better understand the solar genesis of interplanetary magnetic clouds (MCs) we model the magnetic and topological properties of four large eruptive solar flares and relate them to observations. We use the three-dimensional Minimum Current Corona model \cite{Longcope1996d} and observations of pre-flare photospheric magnetic field and flare ribbons to derive values of reconnected magnetic flux, flare energy, flux rope helicity and orientation of the flux rope poloidal field. We compare model predictions of those quantities to flare and MC observations and within the estimated uncertainties of the methods used find the following. The predicted model reconnection fluxes are equal to or lower than the reconnection fluxes inferred from the observed ribbon motions. Both observed and model reconnection fluxes match the MC poloidal fluxes. The predicted flux rope helicities match the MC helicities. The predicted free energies lie between the observed energies and the estimated total flare luminosities. The direction of the leading edge of the MC's poloidal field is aligned with the poloidal field of the flux rope in the AR rather than the global dipole field. These findings compel us to believe that magnetic clouds associated with these four solar flares are formed by low-corona magnetic reconnection during the eruption, rather than eruption of pre-existing structures in the corona or formation in the upper corona with participation of the global magnetic field. We also note that since all four flares occurred in active regions without significant pre-flare flux emergence and cancellation, the energy and helicity we find are stored by shearing and rotating motions, which are sufficient to account for the observed radiative flare energy and MC helicity.

The 21 cm emission and absorption from gaseous halos around the first generation of star depend on the Wouthuysen-Field (W-F) coupling, which relates the spin temperature with the kinetic temperature of hydrogen gas via the resonant scattering between Lyman alpha photons and neutral hydrogen. Although the center object generally is a strong source of these photons, the transfer of these photons in the 21 cm region is inefficient, as the optical depth of the photons is large. Consequently, these photons from the source may not be able to transfer to the entire 21 cm region timely to provide the W-F coupling. This problem is important because the lifetime of first stars generally is short. The problem is investigated with numerical solution of the integro-differential equation, which describes the kinetics of these resonant photons in both physical and frequency spaces. We show that the photon transfer process in the physical space is actually coupled to that in the frequency space. Firstly diffusion in the frequency space provides a shortcut for the diffusion in the physical space. It makes the mean time for the escape of the resonant photon in optical depth \tau media roughly proportional to the optical depth \tau, not \tau^2. Secondly the resonant scattering is effective in bouncing photons with a frequency which is not equal to initial frequency back to the initial frequency. This process can restore initial frequency photons and establish the local Boltzmann distribution of the photon spectrum around the initial frequency. Therefore, the mechanism of 'escape via shortcut' plus 'bounce back' enables W-F coupling to be properly realized in the 21 cm region around first stars. This mechanism also works for photons injected into the 21 cm region by redshift.

AIM: To investigate several partially-erupting prominences to study their relationship with other CME-associated phenomena and to compare these observations with observables predicted by a model of partially-expelled flux ropes (Gibson & Fan, 2006a, b). METHODS: We have studied 6 selected events with partially-erupting prominences using multi wavelength observations recorded by the Extreme-ultraviolet Imaging Telescope (EIT), Transition Region and Coronal Explorer (TRACE), Mauna Loa Solar Observatory (MLSO), Big Bear Solar Observatory (BBSO) and soft X-ray telescope (SXT). The observational features associated with partially-erupting prominences were then compared with the predicted observables from the model. RESULTS: The partially-expelled-flux-rope (PEFR) model of Gibson & Fan (2006a, b) can explain the partial eruption of these prominences, and in addition predicts a variety of other CME-related observables that provide evidence for internal reconnection during eruption. We find that all of the partially-erupting prominences studied in this paper exhibit indirect evidence for internal reconnection. Moreover, all cases showed evidence of at least one observable unique to the PEFR model, e.g., dimmings external to the source region, and/or a soft X-ray cusp overlying a reformed sigmoid. CONCLUSIONS: The PEFR model provides a plausible mechanism to explain the observed evolution of partially-erupting-prominence-associated CMEs in our study.

We develop a numerical solver for the integral-differential equations, which describes the radiative transfer of photon distribution in the frequency space with resonant scattering of Lyalpha photons by hydrogen gas in the early universe. The time-dependent solutions of this equation is crucial to the estimation of the effect of the Wouthuysen-Field (WF) coupling in relation to the 21 cm emission and absorption at the epoch of reionization. The resonant scattering leads to the photon distribution in the frequency space to be piecewise smooth containing sharp changes. The weighted essentially nonoscillatory (WENO) scheme is suitable to handle this problem, as this algorithm has been found to be highly stable and robust for solving Boltzmann equation. We test this numerical solver against analytic solutions of the evolution of the photon distribution in rest background, analytic solution in expanding background without resonant scattering and formation of local Boltzmann distribution around the resonant frequency with the temperature same as that of atom for recoil. We find that evolution of photon distribution undergoes three phases; profile is similar to the initial one, a flat plateau (without recoil) or local Boltzmann distribution (with recoil) forms around the resonant frequency, and finally the distribution around the resonant frequency is saturated when the photons from the source is balanced by the redshift of the expansion. This result indicates that the onset of the W-F coupling should not be determined by the third phase, but by the time scale of the second phase. We found that the time scale of the W-F coupling is equal to about a few hundreds of the mean free flight time of photons with resonant frequency, and is independent of the Sobolev parameter if this parameter is much less than 1.

We study the Wouthuysen-Field coupling at early universe with numerical solutions of the integrodifferential equation describing the kinetics of photons undergoing resonant scattering. The numerical solver is developed based on the weighted essentially non-oscillatory (WENO) scheme for the Boltzmann-like integrodifferential equation. We focus on the time evolution of the Wouthuysen-Field (W-F) coupling in relation to the 21 cm emission and absorption at the epoch of reionization. We show that a local Boltzmann distribution will be formed if photons with frequency \sim \nu_0 have undergone a ten thousand or more times of scattering, which corresponds to the order of 10^3 yrs for neutral hydrogen density of the concordance \Lambda CDM model. The time evolution of the shape and width of the local Boltzmann distribution actually doesn't dependent on the details of atomic recoil, photon sources, or initial conditions very much. However, the intensity of photon flux at the local Boltzmann distribution is substantially time-dependent. The time scale of approaching the saturated intensity can be as long as 10^5-10^6 yrs for typical parameters of the \Lambda CDM model. The intensity of the local Boltzmann distribution at time less than 10^5 yrs is significantly lower than that of the saturation state. Therefore, it may not be always reasonable to assume that the deviation of the spin temperature of 21 cm energy states from cosmic background temperature is mainly due to the W-F coupling if first stars or their emission/absorption regions evolved with a time scale equal to or less than Myrs.

We investigate the volume growth of ionized regions around UV photon sources with the WENO algorithm, which is an effective solver of photon kinetics in the phase space described by the radiative transfer equation. We show that the volume growth rate, either of isolated ionized regions or of clustered regions in merging, generally consists of three phases: fast or relativistic growth phase at the early stage, slow growth phase at the later stage, and a transition phase between the fast and slow phases. We also show that the volume growth of ionized regions around clustered sources with intensity $\dot{E}_i$ ($i=1, 2, ...$) would have the same behavior as a single source with intensity $\dot{E}=\sum_i\dot{E}_i$, if all the distances between nearest neighbor sources $i$ and $j$ are smaller than $c(t^i_c+t^j_c)$, $t^i_c$ being the time scale $t_c$ of source $i$. Therefore, a tightly clustered UV photon sources would lead to a slow growth of ionized volume. This effect would be important for studying the redshift-dependence of 21cm signals from the reionization epoch.

We investigate the 21 cm signals from the UV ionizing sources in the reionization epoch. The formation and evolution of 21 cm emission and absorption regions depend essentially on the kinetics of photons in the physical and frequency spaces. To solve the radiative transfer equation, we use the WENO algorithm, which is effective to capture the sharp ionization profile and the cut-off at the front of light $(r=ct)$ and to handle the small fraction of neutral hydrogen and helium in the ionized sphere. We show that a spherical shell of 21 cm emission and absorption will develop around a point source once the speed of the ionization front (I-front) is significantly lower than the speed of light. The 21 cm shell extends from the I-front to the front of light; its inner part is the emission region and its outer part is the absorption region. The 21 cm emission region depends strongly on the intensity, frequency-spectrum and life-time of the UV ionizing source. For a source of short life-time, no 21 cm emission region can be formed if the source dies out before the I-front speed is significantly lower than the speed of light. Yet, a 21 cm absorption region can form and develop even after the emission of the source ceases.

We develop a numerical solver for radiative transfer problems based on the weighted essentially nonoscillatory (WENO) scheme modified with anti-diffusive flux corrections, in order to solve the temperature and ionization profiles around a point source of photons in the reionization epoch. Algorithms for such simulation must be able to handle the following two features: 1. the sharp profiles of ionization and temperature at the ionizing front (I-front) and the heating front (T-front), and 2. the fraction of neutral hydrogen within the ionized sphere is extremely small due to the stiffness of the rate equations of atom processes. The WENO scheme can properly handle these two features, as it has been shown to have high order of accuracy and good convergence in capturing discontinuities and complicated structures in fluid as well as to be significantly superior over piecewise smooth solutions containing discontinuities. With this algorithm, we show the time-dependence of the preheated shell around a UV photon source. In the first stage the I-front and T-front are coincident, and propagate with almost the speed of light. In later stage, when the frequency spectrum of UV photons is hardened, the speeds of propagation of the ionizing and heating fronts are both significantly less than the speed of light, and the heating front is always beyond the ionizing front. In the spherical shell between the I- and T-fronts, the IGM is heated, while atoms keep almost neutral. The time scale of the preheated shell evolution is dependent on the intensity of the photon source. We also find that the details of the pre-heated shell and the distribution of neutral hydrogen remained in the ionized sphere are actually sensitive to the parameters used. The WENO algorithm can provide stable and robust solutions to study these details.

We apply a new method to determine the magnetic field in coronal loops using observations of coronal loop oscillations. We analyze seven Doppler shift oscillation events detected by SUMER in the hot flare line Fe XIX to obtain oscillation periods of these events. The geometry, temperature, and electron density of the oscillating loops are measured from coordinated multi-channel soft X-ray imaging observations from SXT. All the oscillations are consistent with standing slow waves in their fundamental mode. The parameters are used to calculate the magnetic field of coronal loops based on MHD wave theory. For the seven events, the plasma $\beta$ is in the range 0.15-0.91 with a mean of 0.33$\pm$0.26, and the estimated magnetic field varies between 21-61 G with a mean of 34$\pm$14 G. With background emission subtracted, the estimated magnetic field is reduced by 9%-35%. The maximum backgroud subtraction gives a mean of 22$\pm$13 G in the range 12-51 G. We discuss measurement uncertainties and the prospect of determining coronal loop magnetic fields from future observations of coronal loops and Doppler shift oscillations.

We show that the algorithm based on the weighted essentially nonoscillatory (WENO) scheme with anti-diffusive flux corrections can be used as a solver of the radiative transfer equations. This algorithm is highly stable and robust for solving problems with both discontinuities and smooth solution structures. We test this code with the ionized sphere around point sources. It shows that the WENO scheme can reveal the discontinuity of the radiative or ionizing fronts as well as the evolution of photon frequency spectrum with high accuracy on coarse meshes and for a very wide parameter space. This method would be useful to study the details of the ionized patch given by individual source in the epoch of reionization. We demonstrate this method by calculating the evolution of the ionized sphere around point sources in physical and frequency spaces. It shows that the profile of the fraction of neutral hydrogen and the ionized radius are sensitively dependent on the intensity of the source.