Axion dark matter (DM) may convert to radio-frequency electromagnetic radiation in the strong magnetic fields around neutron stars. The radio signature of such a process would be an ultra-narrow spectral peak at a frequency determined by the mass of the axion particle. We analyze data we collected from the Robert C. Byrd Green Bank Telescope in the L-band and the Effelsberg 100-m Telescope in the L-Band and S-band from a number of sources expected to produce bright signals of axion-photon conversion, including the Galactic Center of the Milky Way and the nearby isolated neutron stars RX J0720.4-3125 and RX J0806.4-4123. We find no evidence for axion DM and are able to set some of the strongest constraints to-date on the existence of axion DM in the highly-motivated mass range between ~5-11 $\mu$eV.

We discuss how reactive and dissipative non-linearities affect the intrinsic response of superconducting thin-film resonators. We explain how most, if not all, of the complex phenomena commonly seen can be described by a model in which the underlying resonance is a single-pole Lorentzian, but whose centre frequency and quality factor change as external parameters, such as readout power and frequency, are varied. What is seen during a vector-network-analyser measurement is series of samples taken from an ideal Lorentzian that is shifting and spreading as the readout frequency is changed. According to this model, it is perfectly proper to refer to, and measure, the resonant frequency and quality factor of the underlying resonance, even though the swept-frequency curves appear highly distorted and hysteretic. In those cases where the resonance curve is highly distorted, the specific shape of the trajectory in the Argand plane gives valuable insights into the second-order physical processes present. We discuss the formulation and consequences of this approach in the case of non-linear kinetic inductance, two-level-system loss, quasiparticle generation, and a generic model based on a power-law form. The generic model captures the key features of specific dissipative non-linearities, but additionally leads to insights into how general dissipative processes create characteristic forms in the Argand plane. We provide detailed formulations in each case, and indicate how they lead to the wide variety of phenomena commonly seen in experimental data. We also explain how the properties of the underlying resonance can be extracted from this data. Overall, our paper provides a self-contained compendium of behaviour that will help practitioners interpret and determine important parameters from distorted swept-frequency measurements.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

Recent studies have indicated that an anisotropic cosmic expansion may exist. In this paper, we use three datasets of type Ia supernovae (SNe Ia) to probe the isotropy of cosmic acceleration. For the Union2.1 dataset, the direction and magnitude of the dipole are $(l=309.3^{\circ} {}^{+ 15.5^{\circ}}_{-15.7^{\circ}} ,\ b = -8.9^{\circ} {}^{ + 11.2^{\circ}}_{-9.8^{\circ}} ),\ A=(1.46 \pm 0.56) \times 10^{-3}$. For the Constitution dataset, the results are $(l=67.0^{\circ}{}^{+ 66.5^{\circ}}_{-66.2^{\circ}},\ b=-0.6^{\circ}{}^{+ 25.2^{\circ}}_{-26.3^{\circ}}),\ A=(4.4 \pm 5.0) \times 10^{-4}$. For the JLA dataset, no significant dipolar deviation is found. We also explore the effects of anisotropic distributions of coordinates and redshifts on the results using Monte-Carlo simulations. We find that the anisotropic distribution of coordinates can cause dipole directions and make dipole magnitude larger. Anisotropic distribution of redshifts is found to have no significant effect on dipole fitting results.

We present several methods towards construction of precursors, which show great promise towards early predictions, of solar flare events in this paper. A data pre-processing pipeline is built to extract useful data from multiple sources, Geostationary Operational Environmental Satellites (GOES) and Solar Dynamics Observatory (SDO)/Helioseismic and Magnetic Imager (HMI), to prepare inputs for machine learning algorithms. Two classification models are presented: classification of flares from quiet times for active regions and classification of strong versus weak flare events. We adopt deep learning algorithms to capture both the spatial and temporal information from HMI magnetogram data. Effective feature extraction and feature selection with raw magnetogram data using deep learning and statistical algorithms enable us to train classification models to achieve almost as good performance as using active region parameters provided in HMI/Space-Weather HMI-Active Region Patch (SHARP) data files. Case studies show a significant increase in the prediction score around 20 hours before strong solar flare events.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

It has been suggested that radio telescopes may be sensitive to axion dark matter that resonantly converts to radio photons in the magnetospheres surrounding neutron stars (NSs). In this work, we closely examine this possibility by calculating the radiated power from and projected sensitivity to axion dark matter conversion in ensembles of NSs within astrophysical systems like galaxies and globular clusters. We use population synthesis and evolution models to describe the spatial distributions of NSs within these systems and the distributions of NS properties. Focusing on three specific targets for illustration, the Galactic Center of the Milky Way, the globular cluster M54 in the Sagittarius dwarf galaxy, and the Andromeda galaxy, we show that narrow-band radio observations with telescopes such as the Green Bank Telescope and the future Square Kilometer Array may be able to probe the quantum chromodynamics axion over roughly two orders of magnitude in mass, starting at a fraction of a $\mu$eV.

One of the main purposes of the DArk Matter Particle Explorer (DAMPE) is to measure the cosmic ray nuclei up to several tens of TeV or beyond, whose origin and propagation remains a hot topic in astrophysics. The Plastic Scintillator Detector (PSD) on top of DAMPE is designed to measure the charges of cosmic ray nuclei from H to Fe and serves as a veto detector for discriminating gamma-rays from charged particles. We propose in this paper a charge reconstruction procedure to optimize the PSD performance in charge measurement. Essentials of our approach, including track finding, alignment of PSD, light attenuation correction, quenching and equalization correction are described detailedly in this paper after a brief description of the structure and operational principle of the PSD. Our results show that the PSD works very well and almost all the elements in cosmic rays from H to Fe are clearly identified in the charge spectrum.

DArk Matter Particle Explorer (DAMPE) is a space-borne apparatus for detecting the high-energy cosmic-rays like electrons, $\gamma$-rays, protons and heavy-ions. Plastic Scintillator Detector (PSD) is the top-most sub-detector of the DAMPE. The PSD is designed to measure the charge of incident high-energy particles and it also serves as a veto detector for discriminating $\gamma$-rays from charged particles. In this paper, PSD on-orbit calibration procedure is described, which includes five steps of pedestal, dynode correlation, response to minimum-ionizing particles (MIPs), light attenuation function and energy reconstruction. A method for reconstructing the charge of incident high energy cosmic-ray particles is introduced. The detection efficiency of each PSD strip is verified to be above 99.5%, the total efficiency of the PSD for charged particles is above 99.99%.

The Plastic Scintillator Detector (PSD) of the DArk Matter Particle Explorer (DAMPE) is designed to measure cosmic ray charge (Z) and to act as a veto detector for gamma-ray identification. In order to fully exploit the charge identification potential of the PSD and to enhance its capability to identify the gamma ray events, we develop a PSD detector alignment method. The path length of a given track in the volume of a PSD bar is derived taking into account the shift and rotation alignment corrections. By examining energy spectra of corner-passing events and fully contained events, position shifts and rotations of all PSD bars are obtained, and found to be on average about 1mm and 0.0015 radian respectively. To validate the alignment method, we introduce the artificial shifts and rotations of PSD bars in the detector simulation. These shift and rotation parameters can be recovered successfully by the alignment procedure. As a result of the PSD alignment procedure, the charge resolution of the PSD is improved from $4\%$ to $8\%$ depending on the nuclei.

In this paper, we study the anisotropy of cosmic acceleration the using Pantheon sample, which includes 1048 spectroscopically confirmed Type Ia supernovae (SNe Ia) covering the redshift range $0.01 < z < 2.3$. In hemisphere comparison method, we find the dipole direction is $(l = 37 \pm 40^{\circ}, b = 33 \pm 16^{\circ})$ with the maximum anisotropy level of $\delta=0.136 {}^{+0.009}_{-0.005}$. From the dipole fitting method, we find that the magnitude of anisotropy is $A = (3.7 {}^{+2.5}_{-3.7}) \times 10^{-4}$, and the direction of the dipole $(l = 329^{\circ}{}^{+ 101^{\circ}}_{-28^{\circ}}, b = 37^{\circ}{}^{+ 52^{\circ}}_{-21^{\circ}})$ in the galactic coordinate system. The result is weakly dependent on redshift from the redshift tomography analysis. The anisotropy is small and the isotropic cosmological model is an excellent approximation.

Two photometric follow-up transit (primary eclipse) observations on WASP-43 b and four observations on TrES-3 b are performed using the Xuyi Near-Earth Object Survey Telescope. After differential photometry and light curve analysis, the physical parameters of the two systems are obtained and are in good match with the literature. Combining with transit data from a lot of literature, the residuals ($O-C$) of transit observations of both systems are fitted with the linear and quadratic functions. With the linear fitting, the periods and transit timing variations (TTVs) of the planets are obtained, and no obvious periodic TTV signal is found in both systems after an analysis. The maximum mass of a perturbing planet located at the 1:2 mean motion resonance (MMR) for WASP-43 b and TrES-3 b is estimated to be 1.826 and 1.504 Earth mass, respectively. By quadratic fitting, it is confirmed that WASP-43 b may have a long-term TTV which means an orbital decay. The decay rate is shown to be $\dot{P} =(-0.005248 \pm 0.001714)$~s$\cdot$yr$^{-1}$, and compared with the previous results. Based on this, the lower limit of the stellar tidal quality parameter of WASP-43 is calculated to be $Q'_{\rm *} \ge 1.5\times10^5$, and the remaining lifetimes of the planets are presented for the different $Q'_{\rm *}$ values of the two systems, correspondingly.

We show that axion dark matter (DM) may be detectable through narrow radio lines emitted from neutron stars. The neutron star magnetosphere hosts a strong magnetic field and a plasma frequency that increases towards the neutron star surface. As the axions pass through the magnetosphere, they can resonantly convert into radio photons in a narrow region around the radius at which the plasma frequency equals the axion mass. The bandwidth of the signal is set by the small DM velocity dispersion far away from the neutron star. We solve the axion-photon mixing equations, including a full treatment of the magnetized plasma and associated anisotropic dielectric tensor, to obtain the conversion probability. We discuss possible neutron-star targets and how they may probe the QCD axion parameter space in the mass range of ~0.2-40 $\mu$eV.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

he DArk Matter Particle Explorer (DAMPE) is a general purposed satellite-borne high energy $\gamma-$ray and cosmic ray detector, and among the scientific objectives of DAMPE are the searches for the origin of cosmic rays and an understanding of Dark Matter particles. As one of the four detectors in DAMPE, the Plastic Scintillator Detector (PSD) plays an important role in the particle charge measurement and the photons/electrons separation. The PSD has 82 modules, each consists of a long organic plastic scintillator bar and two PMTs at both ends for readout, in two layers and covers an overall active area larger than 82 cm $\times$ 82 cm. It can identify the charge states for relativistic ions from H to Fe, and the detector efficiency for Z=1 particles can reach 0.9999. The PSD has been successfully launched with DAMPE on Dec. 17, 2015. In this paper, the design, the assembly, the qualification tests of the PSD and some of the performance measured on the ground have been described in detail.

In this work, we extensively investigate the formation of near 4:2:1 mean motion resonances (MMRs) configuration by performing two sets of N-body simulations. We model the eccentricity damping, gas drag, type I and type II planetary migration of planetesimals, planetary embryos and giant planets in the first sets. For the simulations of giant planets with type II migration, the massive terrestrial planets, with a mass up to several Earth masses, are likely produced in the systems. We further show that by shepherding and/or scattering mechanisms through Jovian planet's type II migration, the terrestrial planets and giant planets in the systems can be evolved into a near 4:2:1 MMRs. Moreover, the models are applicable to the formation of Kepler-238 and 302 systems. In the second set, we study the 4:2:1 MMRs formation in the terrestrial planetary systems, where the planets undergo type I migration and eccentricity damping. By considering type I migration, $\sim$ 17.1\% of the simulations indicate that terrestrial planets are evolved into 4:2:1 MMRs. However, this probability should depend on the initial conditions of planets. Hence, we conclude that both type I and type II migration can play a crucial role in close-in terrestrial planet formation.

We systematically study the renormalizable three-term polynomial inflation in the supersymmetric and non-supersymmetric models. The supersymmetric inflaton potentials can be realized in supergravity theory, and only have two independent parameters. We show that the general renormalizable supergravity model is equivalent to one kind of our supersymmetric models. We find that the spectral index and tensor-to-scalar ratio can be consistent with the Planck and BICEP2 results, but the running of spectral index is always out of the $2\sigma$ range. If we do not consider the BICEP2 experiment, these inflationary models can be highly consistent with the Planck observations and saturate its upper bound on the tensor-to-scalar ratio ($r \le 0.11$). Thus, our models can be tested at the future Planck and QUBIC experiments.

Whether natural factors could interpret the rise of the Earth's surface temperature is still controversial. Though numerous recent researches have reported apparent correlations between solar activity and the Earth's climate, solar activity has encountered a big problem when describing the rapid global warming after 1970s. Our investigation shows the good positive correlations between the Earth's surface Ultraviolet irradiance (280-400 nm) and the Earth's surface temperature both in temporal and spatial variations by analyzing the global surface Ultraviolet irradiance (280-400 nm) and global surface temperature data from 1980-1999. The rise of CO$_2$ cannot interpret the good positive correlations, and we could even get an opposite result to the good correlations when employing the rise of CO$_2$ to describe the relation between them. Based on the good positive correlations, we suggest a new effect, named "Highly Excited Water Vapor" (HEWV) effect, which can interpret how the Sun influences the Earth's surface temperature reasonably, including the rapid warming after 1970s.

Five newly observed transit light curves of the TrES-3 planetary system are presented. Together with other light curve data from literature, 23 transit light curves in total, which cover an overall timescale of 911 epochs, have been analyzed through a standard procedure. From these observational data, the system's orbital parameters are determined and possible transit timing variations are investigated. Given that a null transit-timing-variation produces a fit with reduced chi^2=1.52, our results agree with previous work, that transit timing variations might not exist in these data. However, a 1-frequency oscillating transit-timing-variation model, giving a fit with a reduced chi^2=0.93, does possess a statistically higher probability. It is, thus, concluded that future observations and dynamical simulations for this planetary system will be very important.

Through the international collaborators, we recently established a network of existing and working meter-class telescopes to look for planetary transit events. As a first step, we focus on the TrES3 system, and conclude that there could be some level of transit timing variations.

Properties of proton resonances in $^{18}$Ne have been investigated efficiently by utilizing a technique of proton resonant elastic scattering with a $^{17}$F radioactive ion (RI) beam and a thick proton target. A 4.22~MeV/nucleon $^{17}$F RI beam was produced via a projectile-fragmentation reaction, and subsequently separated by a Radioactive Ion Beam Line in Lanzhou ({\tt RIBLL}). Energy spectra of the recoiled protons were measured by two sets of $\Delta$E-E silicon telescope at center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175${^\circ}$$\pm$5${^\circ}$, $\theta_{c.m.}$$\approx$152${^\circ}$$\pm$8${^\circ}$, respectively. Several proton resonances in $^{18}$Ne were observed, and their resonant parameters have been determined by an $R$-matrix analysis of the differential cross sections in combination with the previous results. The resonant parameters are related to the reaction-rate calculation of the stellar $^{14}$O($\alpha$,$p$)$^{17}$F reaction, which was thought to be the breakout reaction from the hot CNO cycles into the $rp$-process in x-ray bursters. Here, $J^\pi$=(3$^-$, 2$^-$) are tentatively assigned to the 6.15-MeV state which was thought the key 1$^-$ state previously. In addition, a doublet structure at 7.05 MeV are tentatively identified, and its contribution to the resonant reaction rate of $^{14}$O($\alpha$,$p$)$^{17}$F could be enhanced by at least factors of about 4$\sim$6 in comparison with the previous estimation involving only a singlet. The present calculated resonant rates are much larger than those previous values, and it may imply that this breakout reaction could play a crucial role under x-ray bursters conditions.

Understanding the nature of Dark Matter and Dark Energy is one of the most pressing issues in cosmology and fundamental physics. The purpose of the DUNE (Dark UNiverse Explorer) mission is to study these two cosmological components with high precision, using a space-based weak lensing survey as its primary science driver. Weak lensing provides a measure of the distribution of dark matter in the universe and of the impact of dark energy on the growth of structures. DUNE will also include a complementary supernovae survey to measure the expansion history of the universe, thus giving independent additional constraints on dark energy. The baseline concept consists of a 1.2m telescope with a 0.5 square degree optical CCD camera. It is designed to be fast with reduced risks and costs, and to take advantage of the synergy between ground-based and space observations. Stringent requirements for weak lensing systematics were shown to be achievable with the baseline concept. This will allow DUNE to place strong constraints on cosmological parameters, including the equation of state parameter of the dark energy and its evolution from redshift 0 to 1. DUNE is the subject of an ongoing study led by the French Space Agency (CNES), and is being proposed for ESA's Cosmic Vision programme.

Numerical simulation plays an important role in the study of structure formation of the universe. However, the mass resolution in current simulations is still poor. Due to technical difficulties, it is necessary to use both greatly reduced number density of particles and greatly raised unit particle mass. Consequently, the particle masses used in cosmological simulations are about $10^{70}$ times larger than the $Gev$ candidates in particle physics. This is a huge physical bias that cannot be neglected when interpreting the results of the simulations. Here we discuss how such a bias affects Cold Dark Matter (CDM) cosmological simulations. We find that the small scale properties of the CDM particle system are changed in two aspects. 1) An upper limit is imposed on the spatial resolution of the simulation results. 2) Most importantly, an unexpected short mean free path is produced, and the corresponding two body scattering cross section is close to the value expected in the Self-Interaction Dark Matter(SIDM) model. Since the mean free path of real CDM particle systems is much longer than that in the simulations, our results imply 1) that there is probably no 'cusp problem' in real CDM halos, and 2) that a much longer time is needed to form new virialized halos in real CDM particle systems than in the simulations. This last result can help us understand the 'substructure problem'. Our discussion can also explain why the massive halos in the simulations may have smaller concentration coefficients.