Very recently, the Xenon1T collaboration has reported an intriguing electron recoil excess, which may imply for light dark matter. In order to interpret this anomaly, we proposed the the atmospheric dark matter (AMD) from the inelastic collision of cosmic rays with the atmosphere. Thanks to the boost effect of high energy CRs, we show that a ${\cal O}(10)$ keV ADM can be fast-moving and successfully fit the observed electron recoil spectra observed through the ADM-electron scattering process. Meanwhile, our ADM predict the scattering cross sections $\sigma_e \sim {\cal O}(10^{-37}- 10^{-38}$) cm$^{2}$, and thus can evade other direct detection constraints. The search for light meson rare decays, such as $\eta \to \pi + \slashed E_T$, would provide a complementary probe of our ADM in the future.

We perform two-dimensional hydrodynamic simulations to study the non-rotating gas flow in the region from 1-1000 parsec. This region connects the galaxy scale and the accretion disk scale. The gas is irradiated by the photons generated by the accretion disk system near the central black hole with $10^8$ solar mass. We assume that the luminosity of the central accretion system is constant and equals to $0.3$ Eddington luminosity. Gas with density higher than $10^{-24} {\rm g \ cm^{-3}}$ is found to be thermally unstable. Two phase, cold and hot, gas flow will form due to thermal instability. We calculated the ratio of cold gas accretion rate to hot gas accretion rate. This ratio is a function of gas density and generally increases with the increase of gas density. The maximum value of this ratio is $16$, when gas density is $10^{-21} {\rm g \ cm^{-3}} $. Gas with density lower than $10^{-24} {\rm g \ cm^{-3}}$ is found to be thermally stable and the gas flow is one-phase. The applications of the results are briefly discussed.

The recent measurements of the cosmological parameter $H_0$ from the direct local observations and the inferred value from the Cosmic Microwave Background show $\sim 4 \sigma$ discrepancy. We demonstrate that a keV gravitino dark matter, which has a small fraction of non-thermal component (e.g. from the late decay of NLSP bino) may reconcile this tension. Furthermore, we discuss the implied collider signatures and point out that this scenario can be tested by searching for the dilepton plus missing energy events at the LHC in near future.

The Tianlai Cylinder Pathfinder is a radio interferometer array designed to test techniques for 21 cm intensity mapping in the post-reionization Universe, with the ultimate aim of mapping the large scale structure and measuring cosmological parameters such as the dark energy equation of state. Each of its three parallel cylinder reflectors is oriented in the north-south direction, and the array has a large field of view. As the Earth rotates, the northern sky is observed by drift scanning. The array is located in Hongliuxia, a radio-quiet site in Xinjiang, and saw its first light in September 2016. In this first data analysis paper for the Tianlai cylinder array, we discuss the sub-system qualification tests, and present basic system performance obtained from preliminary analysis of the commissioning observations during 2016-2018. We show typical interferometric visibility data, from which we derive the actual beam profile in the east-west direction and the frequency band-pass response. We describe also the calibration process to determine the complex gains for the array elements, either using bright astronomical point sources, or an artificial on site calibrator source, and discuss the instrument response stability, crucial for transit interferometry. Based on this analysis, we find a system temperature of about 90 K, and we also estimate the sensitivity of the array.

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

Detection of non-gravitational interactions of massive dark matter with visible sector so far have given null results. The DM may communicate with the ordinary matter only through gravitational interactions. Besides, the majority of traditional direct detections have poor sensitivity of searching for light DM particle because of the small recoil energy. Thanks to the high energy cosmic rays (CRs), the light DM can be boosted by scattering with CRs and thus may be detected in those existing experiments. In this work, we first derive the exclusion limits of the cosmic ray boosted sub-GeV Dark Matter with gravitational mediator from the Xenon1T data.

The MAJORANA DEMONSTRATOR is sensitive to rare events near its energy threshold, including bosonic dark matter, solar axions, and lightly ionizing particles. In this analysis, a novel training set of low energy small-angle Compton scatter events is used to determine the efficiency of pulse shape analysis cuts, and we present updated bosonic dark matter and solar axion results from an 11.17 kg-y dataset using a 5 keV analysis threshold.

So far most studies on the structure of coronal mass ejections (CMEs) are conducted through white-light coronagraphs, which demonstrate about one third of CMEs exhibit the typical three-part structure in the high corona (e.g., beyond 2 Rs), i.e., the bright front, the dark cavity and the bright core. In this paper, we address the CME structure in the low corona (e.g., below 1.3 Rs) through extreme-ultraviolet (EUV) passbands and find that the three-part CMEs in the white-light images can possess a similar three-part appearance in the EUV images, i.e., a leading edge, a low-density zone, and a filament or hot channel. The analyses identify that the leading edge and the filament or hot channel in the EUV passbands evolve into the front and the core later within several solar radii in the white-light passbands, respectively. What's more, we find that the CMEs without obvious cavity in the white-light images can also exhibit the clear three-part appearance in the EUV images, which means that the low-density zone in the EUV images (observed as the cavity in white-light images) can be compressed and/or transformed gradually by the expansion of the bright core and/or the reconnection of magnetic field surrounding the core during the CME propagation outward. Our study suggests that more CMEs can possess the clear three-part structure in their early eruption stage. The nature of the low-density zone between the leading edge and the filament or hot channel is discussed.

Depending on the value of the accretion rate, black hole accretion is divided into cold and hot modes. The two modes have distinctly different physics and correspond to two feedback modes. Most previous feedback works either focus only on one mode, or the accretion physics is not always properly adopted. Here, by performing hydrodynamical numerical simulations of AGN feedback in an elliptical galaxy, we show that including both is important, and gives different results from including just one or the other. We specifically focus on the wind and radiation (but neglecting the jet) feedback in the hot mode, to explore whether this particular mode of feedback can play any important feedback role. For this aim, we have run two test models. In one model, we always adopt the cold mode no matter what value the accretion rate is; in another model, we turn off the AGN once it enters into the hot mode. We have calculated the AGN light curves, black hole growth, star formation, and AGN energy duty-cycle; and compared the results with the model in which the two modes are correctly included. Important differences are found. For example, if we were to adopt only the cold mode, the total mass of newly formed stars would become two orders of magnitude smaller, and the fraction of energy ejected within the high accretion regime (i.e., $L_{\rm BH} > 2\% L_{\rm Edd}$) would be too small to be consistent with observations. We have also investigated the respective roles of wind and radiation in the hot mode.

An unexpected strong geomagnetic storm occurred on 2018 August 26, which was caused by a slow coronal mass ejection (CME) from a gradual eruption of a large quiet-region filament. We investigate the eruption and propagation characteristics of this CME in relation to the strong geomagnetic storm with remote sensing and in situ observations. Coronal magnetic fields around the filament are extrapolated and compared with EUV observations. We determine the propagation direction and tilt angle of the CME flux rope near the Sun using a graduated cylindrical shell (GCS) model and the Sun-to-Earth kinematics of the CME with wide-angle imaging observations from STEREO A. We reconstruct the flux-rope structure using a Grad-Shafranov technique based on the in situ measurements at the Earth and compare it with those from solar observations and the GCS results. Our conclusions are as follows: (1) the eruption of the filament was unusually slow and occurred in the regions with relatively low critical heights of the coronal field decay index; (2) the axis of the CME flux rope rotated in the corona as well as in interplanetary space, which tended to be aligned with the local heliospheric current sheet; (3) the CME was bracketed between slow and fast solar winds, which enhanced the magnetic field inside the CME at 1 AU; (4) the geomagnetic storm was caused by the enhanced magnetic field and a southward orientation of the flux rope at 1 AU from the rotation of the flux rope.

We investigate a global extreme-ultraviolet (EUV) wave associated with a coronal mass ejection (CME)-driven shock on 2017 September 10. The EUV wave is transmitted by north- and south-polar coronal holes (CHs), which is observed by the Solar Dynamics Observatory (SDO) and Solar Terrestrial Relations Observatory A (STEREO-A) from opposite sides of the Sun. We obtain key findings on how the EUV wave interacts with multiple coronal structures, and on its connection with the CME-driven shock: (1) the transmitted EUV wave is still connected with the shock that is incurvated to the Sun, after the shock has reached the opposite side of the eruption; (2) the south CH transmitted EUV wave is accelerated inside an on-disk, low-density region with closed magnetic fields, which implies that an EUV wave can be accelerated in both open and closed magnetic field regions; (3) part of the primary EUV wavefront turns around a bright point (BP) with a bipolar magnetic structure when it approaches a dim, low-density filament channel near the BP; (4) the primary EUV wave is diffused and apparently halted near the boundaries of remote active regions (ARs) that are far from the eruption, and no obvious AR related secondary waves are detected; (5) the EUV wave extends to an unprecedented scale of ~360{\deg} in latitudes, which is attributed to the polar CH transmission. These results provide insights into the effects of coronal density and magnetic field distributions on the evolution of an EUV wave, and into the connection between the EUV wave and the associated CME-driven shock.

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.

It is paramount from both scientific and societal perspectives to understand the generation of extreme space weather. We discuss the formation of solar superstorms based on a comparative study of the 2012 July 23 and 2017 July 23 eruptions. The first one is Carrington-class, and the second could rival the 1989 March event that caused the most intense geomagnetic storm of the space age. Observations of these events in the historically weak solar cycle 24 indicate that a solar superstorm can occur in any solar cycle and at any phase of the cycle. Recurrent patterns are identified in both cases, including the long-lived eruptive nature of the active region, a complex event composed of successive eruptions from the same active region, and in-transit interaction between the successive eruptions resulting in exceptionally strong ejecta magnetic fields at 1 AU. Each case also shows unique characteristics. Preconditioning of the upstream solar wind leading to unusually high solar wind speeds at 1 AU is observed in the first case whereas not in the latter. This may suggest that the concept of "preconditioning" appears to be necessary for making a Carrington-class storm. We find a considerable deflection by nearby coronal holes in the second case but not in the first. On the basis of these results, we propose a hypothesis for further investigation that superstorms are "perfect storms" in nature, i.e., a combination of circumstances that results in an event of unusual magnitude. Historical records of some extreme events seem to support our hypothesis.

A powerful coronal mass ejection (CME) occurred on 2017 September 10 near the end of the declining phase of the historically weak solar cycle 24. We obtain new insights concerning the geometry and kinematics of CME-driven shocks in relation to their heliospheric impacts from the optimal, multi-spacecraft observations of the eruption. The shock, which together with the CME driver can be tracked from the early stage to the outer corona, shows a large oblate structure produced by the vast expansion of the ejecta. The expansion speeds of the shock along the radial and lateral directions are much larger than the translational speed of the shock center, all of which increase during the flare rise phase, peak slightly after the flare maximum and then decrease. The near simultaneous arrival of the CME-driven shock at the Earth and Mars, which are separated by 156.6$^{\circ}$ in longitude, is consistent with the dominance of expansion over translation observed near the Sun. The shock decayed and failed to reach STEREO A around the backward direction. Comparison between ENLIL MHD simulations and the multi-point in situ measurements indicates that the shock expansion near the Sun is crucial for determining the arrival or non-arrival and space weather impact at certain heliospheric locations. The large shock geometry and kinematics have to be taken into account and properly treated for accurate predictions of the arrival time and space weather impact of CMEs.

Studying the stellar mass, age, luminosity, star-formation rate, and impact parameter of quasar absorber host galaxies can aid in the understanding of galaxy formation and evolution as well as in testing their models. We derive the Spectral Energy Distribution (SED) and impact parameter limits of low redshift ($z_{abs} = 0.37 - 0.55$) Mg II absorbers and of higher redshift ($z_{abs} = 1.0 - 2.5$) 2175 \AA\ dust absorbers (2DAs). We use an imaging stacking technique to statistically boost the signal-to-noise ratio (SNR) to increase detection of the absorber host galaxies. The point spread function of the background quasar is modeled with Principal Component Analysis (PCA). This method efficiently reduces the uncertainty of traditional PSF modeling. Our SED for Mg II absorbers indicates that low redshift Mg II absorber host galaxies are likely star-forming galaxies transitioning into red quiescent galaxies, with a low star formation rate of 2.2 $M_\odot$ $yr^{-1}$. From the stacked images and simulations, we show that the average impact parameter of 2DAs is > 5 times smaller than that of Mg II absorbers, at < 7 kpc instead of Mg II absorbers' 48 kpc, indicating that 2DAs are likely associated with disk components of high redshift galaxies. This means that 2DAs are likely good probes to study precursors to the Milky Way.

The solar eruption on 2012 January 27 resulted in a wide-spread solar energetic particle (SEP) event observed by STEREO A and the near-Earth spacecraft (separated by 108\degree). The event was accompanied by an X-class flare, extreme-ultraviolet (EUV) wave and fast coronal mass ejection (CME). We investigate the particle release by comparing the release times of particles at the spacecraft and the times when magnetic connectivity between the source and the spacecraft was established. The EUV wave propagating to the magnetic footpoint of the spacecraft in the lower corona and the shock expanding to the open field line connecting the spacecraft in the upper corona are thought to be responsible for the particle release. We track the evolution of the EUV wave and model the propagation of the shock using EUV and white-light observations. No obvious evidence indicates that the EUV wave reached the magnetic footpoint of either STEREO A or L1-observers. Our shock modeling shows that the release time of the particles observed at L1 was consistent with the time when the shock first established contact with the magnetic field line connecting L1-observers. The release of the particles observed by STEREO A was delayed relative to the time when the shock was initially connected to STEREO A via the magnetic field line. We suggest that the particle acceleration efficiency of the portion of the shock connected to the spacecraft determines the release of energetic particles at the spacecraft.

The \textsc{Majorana Demonstrator} is an ultra low-background experiment searching for neutrinoless double-beta decay in $^{76}$Ge. The heavily shielded array of germanium detectors, placed nearly a mile underground at the Sanford Underground Research Facility in Lead, South Dakota, also allows searches for new exotic physics. Free, relativistic, lightly-ionizing particles with electrical charges less than $e$ are forbidden by the standard model but predicted by some of its extensions. If such particles exist, they might be detected in the \textsc{Majorana Demonstrator} by searching for multiple- detector events with individual-detector energy depositions down to 1 keV. This search is background free and no candidate events have been found in 285 days of data taking. New direct-detection limits are set for the flux of lightly ionizing particles for charges as low as $e$/1000.

Recently, DAMPE has released its first results on the high-energy cosmic-ray electrons and positrons (CREs) from about $25$ GeV to $4.6$ TeV, which directly detect a break at $\sim 1$ TeV. This result gives us an excellent opportunity to study the source of the CREs excess. In this work, we used the data fo proton and helium flux (from AMS-02 and CREAM), $\bar{\mathrm{p}}/\mathrm{p}$ ratio (from AMS-02), positron flux (from AMS-02) and CREs flux (from DAMPE without the peak signal point at $\sim 1.4$ TeV) to do global fitting simultaneously, which can account the influence from the propagation model, the nuclei and electron primary source injection and the secondary lepton production precisely. For extra source to interpret the excess in lepton spectrum, we consider two separate scenarios (pulsar and dark matter annihilation via leptonic channels) to construct the bump ($\gtrsim 100$ GeV) and the break at $\sim 1$ TeV. The result shows: (i) in pulsar scenario, the spectral index of the injection should be $\nu_{\mathrm{psr}} \sim 0.65$ and the cut-off should be $R_{c} \sim 650$ GV; (ii) in dark matter scenario, the dark matter particle's mass is $m_{\chi} \sim 1208$ GeV and the cross section is $\langle \sigma v \rangle \sim 1.48 \times 10^{-23} \mathrm{cm}^{3} \mathrm{s}^{-1}$. Moreover, in the dark matter scenario, the $\tau \bar{\tau}$ annihilation channel is highly suppressed, and a DM model is built to satisfy the fitting results.

We investigate the coronal and interplanetary evolution of a coronal mass ejection (CME) launched on 2010 September 4 from a source region linking two active regions (ARs) 11101 and 11103, using extreme ultraviolet imaging, magnetogram, white-light and in situ observations from SDO, STEREO, SOHO, VEX and Wind. A potential-field source-surface model is employed to examine the configuration of the coronal magnetic field surrounding the source region. The graduated cylindrical shell model and a triangulation method are applied to determine the kinematics of the CME in the corona and interplanetary space. From the remote sensing and in situ observations we obtain some key results: (1) the CME was deflected in both the eastward and southward directions in the low corona by the magnetic pressure from the two ARs and possibly interacted with another ejection, which caused that the CME arrived at VEX that was longitudinally distant from the source region; (2) although VEX was closer to the Sun, the observed and derived CME arrival times at VEX are not earlier than those at Wind, which suggests the importance of determining both the frontal shape and propagation direction of the CME in interplanetary space; (3) the ICME was compressed in the radial direction while the longitudinal transverse size was extended.

We examine the structure, propagation and expansion of the shock associated with the 2012 July 23 extreme coronal mass ejection (CME). Characteristics of the shock determined from multi-point imaging observations are compared to in situ measurements at different locations and a complex radio type II burst, which according to our definition has multiple branches that may not all be fundamental-harmonic related. The white-light shock signature can be modeled reasonably well by a spherical structure and was expanding backward even on the opposite side of the Sun. The expansion of the shock, which was roughly self-similar after the first $\sim$1.5 hours from launch, largely dominated over the translation of the shock center for the time period of interest. Our study also suggests a bow-shock morphology around the nose at later times due to the outward motion in combination with the expansion of the ejecta. The shock decayed and failed to reach Mercury in the backward direction and STEREO B and Venus in the lateral directions, as indicated by the imaging and in situ observations. The shock in the nose direction, however, may persist to the far outer heliosphere, with predicted impact on Dawn around 06 UT on July 25 and on Jupiter around 23:30 UT on July 27 by an MHD model. The type II burst shows properties generally consistent with the spatial/temporal variations of the shock deduced from imaging and in situ observations. In particular, the low-frequency bands agree well with the in situ measurements of a very low density ahead of the shock at STEREO A.

Nuclear fusion reactions are the most important processes in nature to power stars and produce new elements, and lie at the center of the understanding of nucleosynthesis in the universe. It is critically important to study the reactions in full plasma environments that are close to true astrophysical conditions. By using laser-driven counter-streaming collisionless plasmas, we studied the fusion D$+$D$\rightarrow n +^3$He in a Gamow-like window around 27 keV. The results show that astrophysical nuclear reaction yield can be modulated significantly by the self-generated electromagnetic fields and the collective motion of the plasma. This plasma-version mini-collider may provide a novel tool for studies of astrophysics-interested nuclear reactions in plasma with tunable energies in earth-based laboratories.

We study the solar energetic particle (SEP) event associated with the 2012 July 23 extreme solar storm, for which STEREO and the spacecraft at L1 provide multi-point remote sensing and in situ observations. The extreme solar storm, with a superfast shock and extremely enhanced ejecta magnetic fields observed near 1 AU at STEREO A, was caused by the combination of successive coronal mass ejections (CMEs). Meanwhile, energetic particles were observed by STEREO and near-Earth spacecraft such as ACE and SOHO, suggestive of a wide longitudinal spread of the particles at 1 AU. Combining the SEP observations with in situ plasma and magnetic field measurements we investigate the longitudinal distribution of the SEP event in connection with the associated shock and CMEs. Our results underscore the complex magnetic configuration of the inner heliosphere formed by solar eruptions. The examinations of particle intensities, proton anisotropy distributions, element abundance ratios, magnetic connectivity and spectra also give important clues for the particle acceleration, transport and distribution.

The largest geomagnetic storms of solar cycle 24 so far occurred on 2015 March 17 and June 22 with $D_{\rm st}$ minima of $-223$ and $-195$ nT, respectively. Both of the geomagnetic storms show a multi-step development. We examine the plasma and magnetic field characteristics of the driving coronal mass ejections (CMEs) in connection with the development of the geomagnetic storms. A particular effort is to reconstruct the in situ structure using a Grad-Shafranov technique and compare the reconstruction results with solar observations, which gives a larger spatial perspective of the source conditions than one-dimensional in situ measurements. Key results are obtained concerning how the plasma and magnetic field characteristics of CMEs control the geomagnetic storm intensity and variability: (1) a sheath-ejecta-ejecta mechanism and a sheath-sheath-ejecta scenario are proposed for the multi-step development of the 2015 March 17 and June 22 geomagnetic storms, respectively; (2) two contrasting cases of how the CME flux-rope characteristics generate intense geomagnetic storms are found, which indicates that a southward flux-rope orientation is not a necessity for a strong geomagnetic storm; and (3) the unexpected 2015 March 17 intense geomagnetic storm resulted from the interaction between two successive CMEs plus the compression by a high-speed stream from behind, which is essentially the "perfect storm" scenario proposed by \citet[][i.e., a combination of circumstances results in an event of unusual magnitude]{liu14a}, so the "perfect storm" scenario may not be as rare as the phrase implies.

The short burst GRB 130912A was detected by Swift, Fermi satellites and several ground-based optical telescopes. Its X-ray light curve decayed with time normally. The optical emission, however, displayed a long term plateau, which is the longest one in current short GRB observations. In this work we examine the physical origin of the X-ray and optical emission of this peculiar event. We find that the canonical forward shock afterglow emission model can account for the X-ray and optical data self-consistently and the energy injection model that has been widely adopted to interpret the shallowly-decaying afterglow emission is not needed. We also find that the burst was born in a very-low density interstellar medium, consistent with the compact object merger model. Significant fractions of the energy of the forward shock have been given to accelerate the non-thermal electrons and amplify the magnetic fields (i.e., $\epsilon_{\rm e}\sim 0.37$ and $\epsilon_{\rm B}\sim 0.16$, respectively), which are much larger than those inferred in most short burst afterglow modeling and can explain why the long-lasting optical afterglow plateau is rare in short GRBs.

We collected 19.4 days of data from four 750 g TeO2 bolometers, and in three of them we were able to set the energy threshold around 3 keV using a new analysis technique. We found a background rate ranging from 25 cpd/keV/kg at 3 keV to 2 cpd/keV/kg at 25 keV, and a peak at 4.7 keV. The origin of this peak is presently unknown, but its presence is confirmed by a reanalysis of 62.7 kg.days of data from the finished CUORICINO experiment. Finally, we report the expected sensitivities of the CUORE0 (52 bolometers) and CUORE (988 bolometers) experiments to a WIMP annual modulation signal.

High-precise calculation of near-infrared absorption and transient spectra in bulk semiconductors irradiated by an intense terahertz (THz) laser shows that, the ionization rate of the ground-state exciton, probed through the dynamic Fano resonance, may decrease with the increase of the THz laser intensity. This counterintuitive effect indicates the excitons are stabilized against the field-ionization. The much lower Rydberg energy and ``atomic unit'' of laser intensity for excitons, and the possibility of creating excitons from the ``vacuum state'' allow observing the excitonic stabilization in experiments, in contrast to the case of atomic stabilization.