Small-scale magnetic flux ropes (SFRs) are a type of structures in the solar wind that possess helical magnetic field lines. In a recent report (Chen & Hu 2020), we presented the radial variations of the properties of SFR from 0.29 to 8 au using in situ measurements from the Helios, ACE/Wind, Ulysses, and Voyager spacecraft. With the launch of the Parker Solar Probe (PSP), we extend our previous investigation further into the inner heliosphere. We apply a Grad-Shafranov-based algorithm to identify SFRs during the first two PSP encounters. We find that the number of SFRs detected near the Sun is much less than that at larger radial distances, where magnetohydrodynamic (MHD) turbulence may act as the local source to produce these structures. The prevalence of Alfvenic structures significantly suppresses the detection of SFRs at closer distances. We compare the SFR event list with other event identification methods, yielding a dozen well-matched events. The cross-section maps of two selected events confirm the cylindrical magnetic flux rope configuration. The power-law relation between the SFR magnetic field and heliocentric distances seems to hold down to 0.16 au.

The transonic phenomenon of black hole accretion and the existence of the photon sphere are the characteristics of strong gravitational fields near a black hole horizon. In this work, we study spherical accretion flow onto a general parametrized spherically symmetric black hole spacetimes. For this purpose, we analyze the accretion process of various perfect fluids, such as the isothermal fluid of ultra-stiff, ultra-relativistic, and sub-relativistic types and polytropic fluid, respectively. The influences of extra parameters beyond the Schwarzschild black hole in the general parameterized spherically symmetric black hole on the flow behaviors of the above-mentioned test fluids are studied in detail. In addition, by studying the accretion of ideal photon gas, we further discuss the correspondence between the sonic radius of accreting photon gas and the photon sphere for the general parameterized spherically symmetric black hole. Some possible future extensions of our analysis are also discussed.

The next generation of exoplanet-hunting spectrographs should deliver up to an order of magnitude improvement in radial velocity precision over the standard 1 m/s state of the art. This advance is critical for enabling the detection of Earth-mass planets around Sun-like stars. New calibration techniques such as laser frequency combs and stabilized etalons ensure that the instrumental stability is well characterized. However, additional sources of error include stellar noise, undetected short-period planets, and telluric contamination. To understand and ultimately mitigate error sources, the contributing terms in the error budget must be isolated to the greatest extent possible. Here, we introduce a new high cadence radial velocity program, the EXPRES 100 Earths program, which aims to identify rocky planets around bright, nearby G and K dwarfs. We also present a benchmark case: the 62-d orbit of a Saturn-mass planet orbiting the chromospherically quiet star, HD 3651. The combination of high eccentricity (0.6) and a moderately long orbital period, ensures significant dynamical clearing of any inner planets. Our Keplerian model for this planetary orbit has a residual RMS of 58 cm/s over a $\sim 6$ month time baseline. By eliminating significant contributors to the radial velocity error budget, HD 3651 serves as a standard for evaluating the long term precision of extreme precision radial velocity (EPRV) programs.

As the first successful technique used to detect exoplanets orbiting distant stars, the Radial Velocity Method aims to detect a periodic Doppler shift in a star's spectrum. We introduce a new, mathematically rigorous, approach to detect such a signal that accounts for functional relationships of neighboring wavelengths, minimizes the role of wavelength interpolation, accounts for heteroskedastic noise, and easily allows for statistical inference. Using Hermite-Gaussian functions, we show that the problem of detecting a Doppler shift in the spectrum can be reduced to linear regression in many settings. A simulation study demonstrates that the proposed method is able to accurately estimate an individual spectrum's radial velocity with precision below 0.3 m/s. Furthermore, the new method outperforms the traditional Cross-Correlation Function approach by reducing the root mean squared error up to 15 cm/s. The proposed method is also demonstrated on a new set of observations from the EXtreme PREcision Spectrometer (EXPRES) for the star 51 Pegasi, and successfully recovers estimates that agree well with previous studies of this planetary system. Data and Python3 code associated with this work can be found at https://github.com/parkerholzer/hgrv_method. The method is also implemented in the open source R package rvmethod.

Evidence for the central engine of gamma-ray bursts (GRBs) has been collected in the Neil Gehrels \emph{Swift} data. For instance, some GRBs show an internal X-ray plateau followed by very steep decay, which is difficult to be interpreted within the framework of a black hole (BH) central engine, but are consistent within a rapidly spinning magnetar engine picture. The very steep decay at the end of the plateau suggests a sudden cessation of the central engine, which is explained as the collapse of a supra-massive magnetar into a black hole when it spins down. Here we propose that some additional evidence, such as a second X-ray plateau feature would show up, if the fall-back accretion could activate the newborn BH and sufficient energy could be transferred from the newborn BH to the GRB blast wave. With a systematic data analysis for all long GRBs, we find three candidates in \emph{Swift} sample, i.e., GRBs 070802, 090111, and 120213A, whose X-ray afterglow light curves contain two plateaus, with the first one being an internal plateau. We find that in a fairly loose and reasonable parameter space, the second X-ray plateau data for all 3 GRBs could be well interpreted with our proposed model. Future observations are likely to discover more similar events, which could offer more information of the properties of the magnetar as well as the newborn BH.

An observation-based Galactic Cosmic Ray (GCR) spectral model for heavy nuclei is developed. Zhao and Qin (J. Geophys. Res. Space Phys.118, 1837(2013)) proposed an empirical elemental GCR spectra model for nuclear charge 5-28 over the energy range from 30 to 500 MeV/nuc, which is proved to be successful in predicting yearly averaged GCR heavy nuclei spectra.Based on the latest highly statistically precise measurements from ACE/CRIS,a further elemental GCR model with monthly averaged spectra is presented. The model can reproduce the past and predict the futureGCR intensity monthly by correlating model parameters with thecontinuous sunspot number (SSN) record. The effects of solar activity on GCR modulation are considered separately for odd and even solar cycles. Compared with other comprehensive GCR models, our modeling results are satisfyingly consistent with the GCR spectral measurements from ACE/SIS and IMP-8, and have comparable prediction accuracy as the Badhwar & O'Neill 2014 model.A detailed error analysis is also provided.Finally, the GCR carbon and iron nuclei fluxes for the subsequent two solar cycles (SC 25 and 26) are predicted and they show a potential trend in reduced flux amplitude, which is suspected to be relevant to possible weak solar cycles.

We report detections of atomic species in the atmosphere of MASCARA-2 b, using the first transit observations obtained with the newly commissioned EXPRES spectrograph. EXPRES is a highly stabilised optical echelle spectrograph, designed to detect stellar reflex motions with amplitudes down to 30 cm/s, and was recently deployed at the Lowell Discovery Telescope. By analysing the transmission spectrum of the ultra-hot Jupiter MASCARA-2 b using the cross-correlation method, we confirm previous detections of Fe I, Fe II and Na I, which likely originate in the upper regions of the inflated atmosphere. In addition, we report significant detections of Mg I and Cr II. The absorption strengths change slightly with time, possibly indicating different temperatures and chemistry in the day-side and night-side terminators. Using the effective stellar line-shape variation induced by the transiting planet, we constrain the projected spin-orbit misalignment of the system to $1.6\pm3.1$ degrees, consistent with an aligned orbit. We demonstrate that EXPRES joins a suite of instruments capable of phase-resolved spectroscopy of exoplanet atmospheres.

The nature of dark matter is one of the greatest mysteries in modern physics and astronomy. A wide variety of experiments have been carried out worldwide to search for the evidence of particle dark matter. Chinese physicists started experimental search for dark matter about ten years ago, and have produced results with high scientific impact. In this paper, we present an overview of the dark matter program in China, and discuss recent results and future directions.

The EXtreme PREcision Spectrograph (EXPRES) is an environmentally stabilized, fiber-fed, $R=137,500$, optical spectrograph. It was recently commissioned at the 4.3-m Lowell Discovery Telescope (LDT) near Flagstaff, Arizona. The spectrograph was designed with a target radial-velocity (RV) precision of 30$\mathrm{~cm~s^{-1}}$. In addition to instrumental innovations, the EXPRES pipeline, presented here, is the first for an on-sky, optical, fiber-fed spectrograph to employ many novel techniques---including an "extended flat" fiber used for wavelength-dependent quantum efficiency characterization of the CCD, a flat-relative optimal extraction algorithm, chromatic barycentric corrections, chromatic calibration offsets, and an ultra-precise laser frequency comb for wavelength calibration. We describe the reduction, calibration, and radial-velocity analysis pipeline used for EXPRES and present an example of our current sub-meter-per-second RV measurement precision, which reaches a formal, single-measurement error of 0.3$\mathrm{~m~s^{-1}}$ for an observation with a per-pixel signal-to-noise ratio of 250. These velocities yield an orbital solution on the known exoplanet host 51 Peg that matches literature values with a residual RMS of 0.895$\mathrm{~m~s^{-1}}$.

The EXtreme PREcision Spectrograph (EXPRES) is a new Doppler spectrograph designed to reach a radial velocity measurement precision sufficient to detect Earth-like exoplanets orbiting nearby, bright stars. We report on extensive laboratory testing and on-sky observations to quantitatively assess the instrumental radial velocity measurement precision of EXPRES, with a focused discussion of individual terms in the instrument error budget. We find that EXPRES can reach a single-measurement instrument calibration precision better than 10 cm/s, not including photon noise from stellar observations. We also report on the performance of the various environmental, mechanical, and optical subsystems of EXPRES, assessing any contributions to radial velocity error. For atmospheric and telescope related effects, this includes the fast tip-tilt guiding system, atmospheric dispersion compensation, and the chromatic exposure meter. For instrument calibration, this includes the laser frequency comb (LFC), flat-field light source, CCD detector, and effects in the optical fibers. Modal noise is mitigated to a negligible level via a chaotic fiber agitator, which is especially important for wavelength calibration with the LFC. Regarding detector effects, we empirically assess the impact on radial velocity precision due to pixel-position non-uniformities (PPNU) and charge transfer inefficiency (CTI). EXPRES has begun its science survey to discover exoplanets orbiting G-dwarf and K-dwarf stars, in addition to transit spectroscopy and measurements of the Rossiter-McLaughlin effect.

The Parker Solar Probe (PSP) observed an interplanetary coronal mass ejection (ICME) event during its first orbit around the sun, among many other events. This event is analyzed by applying a wavelet analysis technique to obtain the reduced magnetic helicity, cross helicity, and residual energy, the first two of which are magnetohydrodynamics (MHD) invariants. Our results show that the ICME, as a large scale magnetic flux rope, possesses high magnetic helicity, very low cross helicity, and highly negative residual energy, thus pointing to a magnetic fluctuation dominated structure. Using the same technique, we also search for small-scale coherent magnetic flux rope structures during the period from 2018/10/22--2018/11/21, which are intrinsic to quasi-2D MHD turbulence in the solar wind. Multiple structures with duration between 8 and 300 minutes are identified from PSP in-situ spacecraft measurements. The location and scales of these structures are characterized by wavelet spectrograms of the normalized reduced magnetic helicity, normalized cross helicity and normalized residual energy. Transport theory suggests that these small-scale magnetic flux ropes may contribute to the acceleration of charged particles through magnetic reconnection processes, and the dissipation of these structures may be important for understanding the coronal heating processes.

Based on the early-year observations from Neil Gehrels \emph{Swift} Observatory, \cite{liang07} performed a systematic analysis for the shallow decay component of gamma-ray bursts (GRBs) X-ray afterglow, in order to explore its physical origin. Here we revisit the analysis with an updated sample (with Swift/XRT GRBs between February 2004 and July 2017). We find that with a larger sample, 1) the distributions of the characteristic properties of the shallow decay phase (e.g. $t_{b}$ , $S_{X}$, $\Gamma_{X,1}$, and $\alpha_{X,1}$) still accords with normal or lognormal distribution; 2) $\Gamma_{X,1}$ and $\Gamma_{\gamma}$ still show no correlation, but the tentative correlations of durations, energy fluences, and isotropic energies between the gamma-ray and X-ray phases still exist; 3) for most GRBs, there is no significant spectral evolution between the shallow decay segment and its follow-up segment, and the latter is usually consistent with the external-shock models; 4) assuming that the central engine has a power-law luminosity release history as $L(t)=L_{0}(\frac{t}{t_{0}})^{-q}$, we find that the value $q$ is mainly distributed between -0.5 and 0.5, with an average value of 0.16$\pm$ 0.12; 5) the tentative correlation between $E_{\rm{iso},X}$ and $t'_{b}$ disappears, so that the global 3-parameter correlation ($E_{\rm{iso},X}-E'_{p}-t'_{b}$) becomes less significant; 6) the anti-correlation between $L_{X}$ and $t'_{b}$ and the three-parameter correlation ($E_{\rm{iso},\gamma}-L_{X}-t_{b}$) indeed exist with a high confidence level. Overall, our results are generally consistent with \cite{liang07}, confirming their suggestion that the shallow decay segment in most bursts is consistent with an external forward shock origin, probably due to a continuous energy injection from a long-lived central engine.

We calculate the interplanetary magnetic field path lengths traveled by electrons in solar electron events detected by the WIND 3DP instrument from $1994$ to $2016$. The velocity dispersion analysis method is applied for electrons at energies of $\sim$ $27$ keV to $310$ keV. Previous velocity dispersion analyses employ the onset times, which are often affected by instrumental effects and the pre-existing background flux, leading to large uncertainties. We propose a new method here. Instead of using the peak or onset time, we apply the velocity dispersion analysis to the times that correspond to the rising phase of the fluxes that are a fraction, $\eta$, of the peak flux. We perform statistical analysis on selected events whose calculated path lengths have uncertainties smaller than $0.1$ AU. The mean and standard deviation, ($\mu$, $\sigma$), of the calculated path lengths corresponding to $\eta=$ $3/4$, $1/2$, and $1/3$ of the peak flux is ($1.17$ AU, $0.17$ AU), ($1.11$ AU, $0.14$ AU), and ($1.06$ AU, $0.15$ AU). The distribution of the calculated path lengths is also well fitted by a Gaussian distribution for the $\eta=3/4$ and $1/3$ cases. These results suggest that in these electron events the interplanetary magnetic field topology is close to the nominal Parker spiral with little field line meandering. Our results have important implications for particles' perpendicular diffusion.

We have developed a new automated small-scale magnetic flux rope (SSMFR) detection algorithm based on the Grad-Shafranov (GS) reconstruction technique. We have applied this detection algorithm to the Wind spacecraft in-situ measurements during 1996 - 2016, covering two solar cycles, and successfully detected a total number of 74,241 small-scale magnetic flux rope events with duration from 9 to 361 minutes. This large number of small-scale magnetic flux ropes has not been discovered by any other previous studies through this unique approach. We perform statistical analysis of the small-scale magnetic flux rope events based on our newly developed database, and summarize the main findings as follows. (1) The occurrence of small-scale flux ropes has strong solar cycle dependency with a rate of a few hundreds per month on average. (2) The small-scale magnetic flux ropes in the ecliptic plane tend to align along the Parker spiral. (3) In low speed ($<$ 400 km/s) solar wind, the flux ropes tend to have lower proton temperature and higher proton number density, while in high speed ($\ge$ 400 km/s) solar wind, they tend to have higher proton temperature and lower proton number density. (4) Both the duration and scale size distributions of the small-scale magnetic flux ropes obey a power law. (5) The waiting time distribution of small-scale magnetic flux ropes can be fitted by an exponential function (for shorter waiting times) and a power law function (for longer waiting times). (6) The wall-to-wall time distribution obeys double power laws with the break point at 60 minutes (corresponding to the correlation length). (7) The small-scale magnetic flux ropes tend to accumulate near the heliospheric current sheet (HCS). The entire database is available at \url{http://fluxrope.info} and in machine readable format in this article.

Extensive resonance structures are manifest in R-Matrix (RM) calculations. However, there exist a large number of highly excited electronic configurations that may contribute to background non-resonant bound-free opacity in high-temperature plasmas. Since RM calculations are very complex, and not essential for background contributions, the Relativistic Distorted Wave (RDW) method is utilized to complement ("top-up") photoionization cross sections of Fe XVII obtained using Close-Coupling Breit-Pauli R-Matrix (CC-BPRM) method. There is good agreement between RDW and BPRM for background cross sections where resonances are not present, and individual fine structure levels can be correctly matched spectroscopically, though resonances are neglected in the RDW. To ensure completeness, a high energy range up to 500 Ry above the ionization threshold for each level is considered. Interestingly, the hydrogenic Kramer's approximation also shows the same energy behavior as the RDW. Grouping separately, the BPRM configurations consist of 454 bound levels with resonances corresponding to configurations $1s^22s^22p^4nln'l'$ (n $\leq$ 3, n' $\leq$ 10); including RDW configurations there are 51,558 levels in total. The topup contribution results in $\sim$20\% increment, in addition to the 35\% enhancement from BPRM calculations over the Opacity Project value for the Rosseland Mean Opacity at the Z-temperature of 2.11 $\times 10^6$K (Pradhan and Nahar 2017).

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.

We use more than a decade of radial velocity measurements for $\alpha$ Cen A, B, and Proxima Centauri from HARPS, CHIRON, and UVES to identify the $M \sin i$ and orbital periods of planets that could have been detected if they existed. At each point in a mass-period grid, we sample a simulated, Keplerian signal with the precision and cadence of existing data and assess the probability that the signal could have been produced by noise alone. Existing data places detection thresholds in the classically defined habitable zones at about $M \sin i$ of 53 M$_{\oplus}$ for $\alpha$ Cen A, 8.4 M$_{\oplus}$ for $\alpha$ Cen B, and 0.47 M$_{\oplus}$ for Proxima Centauri. Additionally, we examine the impact of systematic errors, or "red noise" in the data. A comparison of white- and red-noise simulations highlights quasi-periodic variability in the radial velocities that may be caused by systematic errors, photospheric velocity signals, or planetary signals. For example, the red-noise simulations show a peak above white-noise simulations at the period of Proxima Centauri b. We also carry out a spectroscopic analysis of the chemical composition of the $\alpha$ Centauri stars. The stars have super-solar metallicity with ratios of C/O and Mg/Si that are similar to the Sun, suggesting that any small planets in the $\alpha$ Cen system may be compositionally similar to our terrestrial planets. Although the small projected separation of $\alpha$ Cen A and B currently hampers extreme-precision radial velocity measurements, the angular separation is now increasing. By 2019, $\alpha$ Cen A and B will be ideal targets for renewed Doppler planet surveys.

We report here the results of searching for inelastic scattering of dark matter (initial and final state dark matter particles differ by a small mass splitting) with nucleon with the first 79.6-day of PandaX-II data (Run 9). We set the upper limits for the spin independent WIMP-nucleon scattering cross section up to a mass splitting of 300 keV/c$^2$ at two benchmark dark matter masses of 1 and 10 TeV/c$^2$.

We report new searches for the solar axions and galactic axion-like dark matter particles, using the first low-background data from PandaX-II experiment at China Jinping Underground Laboratory, corresponding to a total exposure of about $2.7\times 10^4$ kg$\cdot$day. No solar axion or galactic axion-like dark matter particle candidate has been identified. The upper limit on the axion-electron coupling ($g_{Ae}$) from the solar flux is found to be about $4.35 \times 10^{-12}$ in mass range from $10^{-5}$ to 1 keV/$c^2$ with 90\% confidence level, similar to the recent LUX result. We also report a new best limit from the $^{57}$Fe de-excitation. On the other hand, the upper limit from the galactic axions is on the order of $10^{-13}$ in the mass range from 1 keV/$c^2$ to 10 keV/$c^2$ with 90\% confidence level, slightly improved compared with the LUX.

We report a new search of weakly interacting massive particles (WIMPs) using the combined low background data sets in 2016 and 2017 from the PandaX-II experiment in China. The latest data set contains a new exposure of 77.1 live day, with the background reduced to a level of 0.8$\times10^{-3}$ evt/kg/day, improved by a factor of 2.5 in comparison to the previous run in 2016. No excess events were found above the expected background. With a total exposure of 5.4$\times10^4$ kg day, the most stringent upper limit on spin-independent WIMP-nucleon cross section was set for a WIMP with mass larger than 100 GeV/c$^2$, with the lowest exclusion at 8.6$\times10^{-47}$ cm$^2$ at 40 GeV/c$^2$.

New constraints are presented on the spin-dependent WIMP-nucleon interaction from the PandaX-II experiment, using a data set corresponding to a total exposure of 3.3$\times10^4$ kg-days. Assuming a standard axial-vector spin-dependent WIMP interaction with $^{129}$Xe and $^{131}$Xe nuclei, the most stringent upper limits on WIMP-neutron cross sections for WIMPs with masses above 10 GeV/c$^{2}$ are set in all dark matter direct detection experiments. The minimum upper limit of $4.1\times 10^{-41}$ cm$^2$ at 90\% confidence level is obtained for a WIMP mass of 40 GeV/c$^{2}$. This represents more than a factor of two improvement on the best available limits at this and higher masses. These improved cross-section limits provide more stringent constraints on the effective WIMP-proton and WIMP-neutron couplings.

Solar filaments/prominences are one of the most common features in the corona, which may lead to energetic coronal mass ejections (CMEs) and flares when they erupt. Filaments are about one hundred times cooler and denser than the coronal material, and physical understanding of their material origin remains controversial. Two types of scenarios have been proposed: one argues that the filament plasma is brought into the corona from photosphere or chromosphere through a siphon or evaporation/injection process, while the other suggests that the material condenses from the surrounding coronal plasma due to thermal instability. The elemental abundance analysis is a reasonable clue to constrain the models, as the siphon or evaporation/injection model would predict that the filament material abundances are close to the photospheric or chromospheric ones, while the condensation model should have coronal abundances. In this letter, we analyze the elemental abundances of a magnetic cloud that contains the ejected filament material. The corresponding filament eruption occurred on 1998 April 29, accompanying an M6.8 class soft X-ray flare located at the heliographic coordinates S18E20 (NOAA 08210) and a fast halo CME with the linear velocity of 1374 km s$^{-1}$ near the Sun. We find that the abundance ratios of elements with low and high First Ionization Potential such as Fe/O, Mg/O, and Si/O are 0.150, 0.050, and 0.070, respectively, approaching their corresponding photospheric values 0.065, 0.081, and 0.066, which does not support the coronal origin of the filament plasma.

In large SEP events, ions can be accelerated at CME-driven shocks to very high energies. Spectra of heavy ions in many large SEP events show features such as roll-overs or spectral breaks. In some events when the spectra are plotted in energy/nucleon they can be shifted relative to each other to make the spectral breaks align. The amount of shift is charge-to-mass ratio (Q/A) dependent and varies from event to event. This can be understood if the spectra of heavy ions are organized by the diffusion coefficients (Cohen et al., 2005). In the work of Li et al. (2009), the Q/A dependences of the scaling is related to shock geometry when the CME-driven shock is close to the Sun. For events where multiple in-situ spacecraft observations exist, one may expect that different spacecraft are connected to different portions of the CME-driven shock that have different shock geometries, therefore yielding different Q/A dependence. In this work, we examine one SEP event which occurred on 2013 November 4. We study the Q/A dependence of the energy scaling for heavy ion spectra using Helium, oxygen, and iron ions. Observations from STEREO-A, STEREO-B and ACE are examined. We find that the scalings are different for different spacecraft. We suggest that this is because ACE, STEREO-A and STEREO- B are connected to different parts of the shock that have different shock geometries. Our analysis indicates that studying the Q/A scaling of in-situ particle spectra can serve as a powerful tool to remotely examine the shock geometry for large SEP events.

Magnetic clouds (MCs) are the interplanetary counterpart of coronal magnetic flux ropes. They can provide valuable information to reveal the flux rope characteristics at their eruption stage in the corona, which are unable to be explored in situ at present. In this paper, we make a comprehensive survey of the average iron charge state (<Q>Fe) distributions inside 96 MCs for solar cycle 23 using ACE (Advanced Composition Explorer) data. As the <Q>Fe in the solar wind are typically around 9+ to 11+, the Fe charge state is defined as high when the <Q>Fe is larger than 12+, which implies the existence of a considerable amount of Fe ions with high charge states (e.g., \geq 16+). The statistical results show that the <Q>Fe distributions of 92 (~ 96%) MCs can be classified into four groups with different characteristics. In group A (11 MCs), the <Q>Fe shows a bimodal distribution with both peaks higher than 12+. Group B (4 MCs) presents a unimodal distribution of <Q>Fe with its peak higher than 12+. In groups C (29 MCs) and D (48 MCs), the <Q>Fe remains higher and lower than 12+ throughout ACE passage through the MC, respectively. Possible explanations to these distributions are discussed.

The Large High Altitude Air Shower Observatory (LHAASO) is designed for high energy gamma ray and cosmic ray detection. A Water Cherenkov Detector Array which is sensitive to gamma ray showers above a few hundred GeV is proposed to survey gamma ray sources. The WCDA consists of 3600 PhotoMultiplier Tubes (PMT) which collect the Cherenkov light produced by the shower particles in water. Both high precision time and charge measurement are required over a large dynamic range from 1 photo electron (P.E.) to 4000 P.E. Prototype of an analogue front-end Application Specific Integrated Circuit (ASIC) fabricated in Chartered 0.35 {\mu}m CMOS technology is designed to read out PMT signal in the WCDA. This ASIC employs leading edge discrimination and RC4 shaping structure; combined with the following Time-to-Digital Converter (TDC) and Analog-to-Digital Converter (ADC), both the arrival time and charge of the PMT signal can be measured. Initial test results indicate that time resolution is better than 350 ps and charge resolution is better than 10% at 1 P.E. and better than 1% with large input signals (300 P.E. to 4000 P.E.). Besides, this ASIC has a small channel-to-channel crosstalk and low ambient temperature dependency.

We report the results of a weakly-interacting massive particle (WIMP) dark matter search using the full 80.1\;live-day exposure of the first stage of the PandaX experiment (PandaX-I) located in the China Jin-Ping Underground Laboratory. The PandaX-I detector has been optimized for detecting low-mass WIMPs, achieving a photon detection efficiency of 9.6\%. With a fiducial liquid xenon target mass of 54.0\,kg, no significant excess event were found above the expected background. A profile likelihood analysis confirms our earlier finding that the PandaX-I data disfavor all positive low-mass WIMP signals reported in the literature under standard assumptions. A stringent bound on the low mass WIMP is set at WIMP mass below 10\,GeV/c$^2$, demonstrating that liquid xenon detectors can be competitive for low-mass WIMP searches.

We report on the first dark-matter (DM) search results from PandaX-I, a low threshold dual-phase xenon experiment operating at the China Jinping Underground Laboratory. In the 37-kg liquid xenon target with 17.4 live-days of exposure, no DM particle candidate event was found. This result sets a stringent limit for low-mass DM particles and disfavors the interpretation of previously-reported positive experimental results. The minimum upper limit, $3.7\times10^{-44}$\,cm$^2$, for the spin-independent isoscalar DM-particle-nucleon scattering cross section is obtained at a DM-particle mass of 49\,GeV/c$^2$ at 90\% confidence level.

On 2012 July 11, two solar filaments were observed in the northeast of the solar disk and their eruptions due to the interaction between them are studied by using the data from the Solar Dynamics Observatory (SDO), Solar TErrestrial RElations Observatory (STEREO) and Global Oscillation Network Group (GONG). The eastern filament (F1) first erupted toward the northeast. During the eruption of F1, some plasma from F1 fell down and was injected to the North-East part of another filament (F2), and some plasma of F1 fell down to the northern region close to F2 and caused the plasma to brighten. Meanwhile, the North-East part of F2 first started to be active and rise, but did not erupt finally. Then the South-West part of F2 erupted successfully. Therefore, the F2's eruption is a partial filament eruption. Two associated CMEs related to the eruptions were observed by STEREO/COR1. We find two possible reasons that lead to the instability and the eruption of F2. One main reason is that the magnetic loops overlying the two filaments were partially opened by the eruptive F1 and resulted in the instability of F2. The other is that the downflows from F1 might break the stability of F2.

Seeking to establish whether active region upflow material contributes to the slow solar wind, we examine in detail the plasma upflows from Active Region (AR)10978, which crossed the Sun's disc in the interval 8 to 16 December, 2007 during Carrington rotation (CR)2064. In previous work, using data from the Hinode/EUV Imaging Spectrometer, upflow velocity evolution was extensively studied as the region crossed the disc while a linear force-free magnetic extrapolation was used to confirm aspects of the velocity evolution and to establish the presence of quasi-separatrix layers at the upflow source areas. The plasma properties, temperature, density and first ionisation potential bias (FIP-bias) were measured with the spectrometer during the disc passage of the active region. Global potential field source surface (PFSS) models showed that AR 10978 was completely covered by the closed field of a helmet streamer that is part of the streamer belt. Thus it is not clear how any of the upflowing AR-associated plasma could reach the source surface at 2.5 R(Sun)and contribute to the slow solar wind. However a detailed examination of solar-wind in-situ data obtained by the Advanced Composition Explorer (ACE) spacecraft at the L1 point shows that the increase in O^7+/O^6+, C^6+/C^5+ and Fe/O - a FIP-bias proxy - are present before the heliospheric current sheet crossing. These increases, along with an accompanying reduction in proton velocity and an increase in density are characteristic of both AR and slow-wind plasma. Finally we describe a two-step reconnection process by which some of the upflowing plasma from the AR could reach the heliosphere.

During the recent solar minimum between cycles 23 and 24 (solar minimum $P_{23/24}$) the intensity of Galactic Cosmic Rays (GCRs) measured at the Earth was the highest ever recorded since space age. It is the purpose of this paper to resolve the most plausible mechanism for this unusually high intensity. A GCR transport model in three-dimensional heliosphere based on a simulation of Markov stochastic process is used to find the relation of cosmic ray modulation to various transport parameters, including solar wind (SW) speed, distance of heliospheric boundary, magnitude of interplanetary magnetic field (IMF) at the Earth, tilt angle of heliospheric current sheet (HCS), values of parallel and perpendicular diffusion coefficients. We calculate GCR proton energy spectra at the Earth for the last three solar minima $P_{21/22}$, $P_{22/23}$, and $P_{23/24}$, with the transport parameters obtained from observations. Besides weak IMF magnitude and slow SW speed, we find that a possible low magnetic turbulence, which increases the parallel diffusion and reduces the perpendicular diffusion in the polar direction, might be an additional possible mechanism for the high GCR intensity in the solar minimum $P_{23/24}$.

In this paper, a new version of the interacting model of new agegraphic dark energy (INADE) is proposed and analyzed in detail. The interaction between dark energy and dark matter is reconsidered. The interaction term $Q=bH_0\rho_{\rm de}^\alpha\rho_{\rm dm}^{1-\alpha}$ is adopted, which abandons the Hubble expansion rate $H$ and involves both $\rho_{\rm de}$ and $\rho_{\rm dm}$. Moreover, the new initial condition for the agegraphic dark energy is used, which solves the problem of accommodating baryon matter and radiation in the model. The solution of the model can be given using an iterative algorithm. A concrete example for the calculation of the model is given. Furthermore, the model is constrained by using the combined Planck data (Planck+BAO+SNIa+$H_0$) and the combined WMAP-9 data (WMAP+BAO+SNIa+$H_0$). Three typical cases are considered: (A) $Q=bH_0\rho_{\rm de}$, (B) $Q=bH_0\sqrt{\rho_{\rm de}\rho_{\rm dm}}$, and (C) $Q=bH_0\rho_{\rm dm}$, which correspond to $\alpha=1$, 1/2, and 0, respectively. The departures of the models from the $\Lambda$CDM model are measured by the $\Delta$BIC and $\Delta$AIC values. It is shown that the INADE model is better than the NADE model in the fit, and the INADE(A) model is the best in fitting data among the three cases.

We report observations of the acceleration and trapping of energetic ions and electrons between a pair of corotating interaction regions (CIRs). The event occurred in Carrington Rotation 2060. Observed at spacecraft STEREO-B, the two CIRs were separated by less than 5 days. In contrast to other CIR events, the fluxes of energetic ions and electrons in this event reached their maxima between the trailing-edge of the first CIR and the leading edge of the second CIR. The radial magnetic field (Br) reversed its sense and the anisotropy of the flux also changed from sunward to anti-sunward between the two CIRs. Furthermore, there was an extended period of counter-streaming suprathermal electrons between the two CIRs. Similar observations for this event were also obtained for ACE and STEREO-A. We conjecture that these observations were due to a "U-shape" large scale magnetic field topology connecting the reverse shock of the first CIR and the forward shock of the second CIR. Such a disconnected U-shaped magnetic field topology may have formed due to magnetic reconnection in the upper corona.

We report observations of the acceleration and trapping of energetic ions and electrons between a pair of corotating interaction regions (CIRs). The event occurred in Carrington Rotation 2060. Observed at spacecraft STEREO-B, the two CIRs were separated by less than 5 days. In contrast to other CIR events, the fluxes of energetic ions and electrons in this event reached their maxima between the trailing-edge of the first CIR and the leading edge of the second CIR. The radial magnetic field (Br) reversed its sense and the anisotropy of the flux also changed from sunward to anti-sunward between the two CIRs. Furthermore, there was an extended period of counter-streaming suprathermal electrons between the two CIRs. Similar observations for this event were also obtained for ACE and STEREO-A. We conjecture that these observations were due to a "U-shape" large scale magnetic field topology connecting the reverse shock of the first CIR and the forward shock of the second CIR. Such a disconnected U-shaped magnetic field topology may have formed due to magnetic reconnection in the upper corona.

Cosmic Rays (CRs) include Galactic Cosmic Rays (GCRs) and Anomalous Cosmic Rays (ACRs). The CR flux data of protons and heavy-ions observed with spacecraft are often seriously contaminated by Solar Energetic Particle (SEP) events. In this work, we separate SEPs from CRs of ACE/SIS spacecraft observations with an automatic despiking algorithm, so we are able to study the different variations of SEPs and CRs over a solar cycle. In particular, we study the elemental ratio, first ionization potential dependence and elemental dependence, and information entropy of SEPs and CRs. So that we can gain new insights into energetic particles' different compositions, origins, and transport processes, etc.

The LHAASO (Large High Altitude Air Shower Observatory) experiment is proposed for very high energy gamma ray source survey, in which the WCDA (Water Cherenkov Detector Array) is the one of the major components. In the WCDA, a total of 3600 PMTs are placed under water in four ponds, each with a size of 150 m x 150 m. Precise time and charge measurement is required for the PMT signals, over a large signal amplitude range from single P.E. (photo electron) to 4000 P.E. To fulfill the high requirement of signal measurement in so many front end nodes scattered in a large area, special techniques are developed, such as multiple gain readout, hybrid transmission of clocks, commands, and data, precise clock phase alignment, and new trigger electronics. We present the readout electronics architecture for the WCDA and several prototype modules, which are now under test in the laboratory.

To explore the difference between the most two recent solar minima, we analyze the in-situ ACE and ULYSSES observations and examine the distributions of the three types of solar wind (streamer-stalk-associated wind, wind from outside the streamer stalk that can be associated, in part, with coronal holes, and interplanetary coronal mass ejections). We use the taxonomy provided by Zhao et al. (2009) to identify the three types of solar wind. We then map the in-situ observations to the 2.5 solar radii surface. With the aid of the potential-field-source surface model (PFSS), we calculate the normal distance from the solar wind "foot point" to the local helisopheric current sheet on that surface. We find that the source region of the streamer stalk wind is narrower compared to the previous minimum. The area outside the streamer stalk is accordingly larger, but the magnetic field strength is observed to be lower, with the result that the total amount of the magnetic open flux from the outside of streamer stalk region is conserved in the two successive solar minima. The implications of the conservation of open magnetic flux for models of the behavior of the solar magnetic field are discussed.

We present results for Chandra observations of comets, 17P/Holmes (17P) and 8P/Tuttle (8P). 17P was observed for 30 ksec right after its major outburst, on 31 Oct 2007 (10:07 UT) and comet 8P/Tuttle was observed in 2008 January for 47 ksec. During the two Chandra observations, 17P was producing at least 100 times more water than 8P but was 2.2 times further away from the Sun. Also, 17P is the first comet observed at high latitude (+19.1 degrees) during solar minimum, while 8P was observed at a lower solar latitude (3.4 degrees). The X-ray spectrum of 17P is unusually soft with little significant emission at energies above 500 eV. Depending on our choice of background, we derive a 300 to 1000 eV flux of 0.5 to 4.5 x 10^-13 ergs/cm2/sec, with over 90% of the emission in the 300 to 400 eV range. This corresponds to an X-ray luminosity between 0.4 to 3.3 x 10^15 ergs/sec. 17P's lack of X-rays in the 400 to 1000 eV range, in a simple picture, may be attributed to the polar solar wind, which is depleted in highly charged ions. 8P/Tuttle was much brighter, with an average count rate of 0.20 counts/s in the 300 to 1000 eV range. We derive an average X-ray flux in this range of 9.4 x 10^-13 ergs/cm2/sec and an X-ray luminosity for the comet of 1.7 x 10^14 ergs/sec. The light curve showed a dramatic decrease in flux of over 60% between observations on January 1st and 4th. When comparing outer regions of the coma to inner regions, its spectra showed a decrease in ratios of CVI/CV, OVIII/OVII, as predicted by recent solar wind charge exchange emission models. There are remarkable differences between the X-ray emission from these two comets, further demonstrating the qualities of cometary X-ray observations, and solar wind charge exchange emission in more general as a means of remote diagnostics of the interaction of astrophysical plasmas.

We investigate the physical meaning of the prominent blue shifts of Ne VIII, which is observed to be associated with quiet-Sun network junctions (boundary intersections), through data analyses combining force-free-field extrapolations with EUV spectroscopic observations. For a middle-latitude region, we reconstruct the magnetic funnel structure in a sub-region showing faint emission in EIT-Fe 195. This funnel appears to consist of several smaller funnels that originate from network lanes, expand with height and finally merge into a single wide open-field region. However, the large blue shifts of Ne VIII are generally not associated with open fields, but seem to be associated with the legs of closed magnetic loops. Moreover, in most cases significant upflows are found in both of the funnel-shaped loop legs. These quasi-steady upflows are regarded as signatures of mass supply to the coronal loops rather than the solar wind. Our observational result also reveals that in many cases the upflows in the upper transition region (TR) and the downflows in the middle TR are not fully cospatial. Based on these new observational results, we suggest different TR structures in coronal holes and in the quiet Sun.

The merging cluster of galaxies A2255 is covered by the Sloan Digital Sky Survey (SDSS) survey. In this paper we perform a morphological classification on the basis of the SDSS imaging and spectral data, and investigate the morphological dependence of the star formation rates (SFRs) for these member galaxies. As we expect, a tight correlation between the normalized SFR by stellar mass (SFR/M$_*$) and the H$\alpha$ equivalent width is found for the late-type galaxies in A2255. The correlation of SFR/M$_*$ with the continuum break strength at 4000 \AA is also confirmed. The SFR/M$_*$ - M$_*$ correlation is found for both the early- and late-type galaxies, indicating that the star formation activity tends to be suppressed when the assembled stellar mass M$_*$) increases, and this correlation is tighter and steeper for the late-type cluster galaxies. Compared with the mass range of field spiral galaxies, only two massive late-type galaxies with M$_*>10^{11}$ M$_{\odot}$ are survived in A2255, suggesting that the gas disks of massive spiral galaxies could have been tidally stripped during cluster formation. Additionally, the SFR variation with the projected radial distance are found to be heavily dependent upon galaxy morphology: the early-type galaxies have a very weak inner decrease in SFR/M$_*$, while the inner late-type galaxies tend to have higher SFR/M$_*$ values than the outer late-types. This may suggest that the galaxy-scale turbulence stimulated by the merging of subclusters might have played different roles on early- and late-type galaxies, which leads to a suppression of the star formation activity for E/S0 galaxies and a SFR enhancement for spiral and irregular galaxies.

Rapid, efficient deactivation of the photoexcited PAH cations accounts for their remarkable photostability and have important implications for astrochemistry, as these cations are the leading candidates for the species responsible for the diffuse interstellar bands (DIB) observed throughout the Galaxy.Ultrafast relaxation dynamics for photoexcited PAH cations isolated in boric acid glass have been studied using femtosecond and picosecond transient grating spectroscopy. With the exception of perylene+, the recovery kinetics for the ground doublet (D0) states of these radical cations are biexponential, containing a fast (< 200 fs) and a slow (3-20 ps) components. No temperature dependence or isotope effect was observed for the fast component, whereas the slow component exhibits both the H/D isotope effect (1.1-1.3) and strong temperature dependence (15 to 300 K). We suggest that the fast component is due to internal Dn to D0 conversion and the slow component is due to vibrational energy transfer (VET) from a hot D0 state to the glass matrix.