We report results from searches for new physics with low-energy electronic recoil data recorded with the XENON1T detector, with an exposure of 0.65 t-y and an unprecedentedly low background rate of $76 \pm 2_{stat}$ events/(t y keV) between 1-30 keV. An excess over known backgrounds is observed below 7 keV, rising towards lower energies and prominent between 2-3 keV. The solar axion model has a 3.5$\sigma$ significance, and a three-dimensional 90% confidence surface is reported for axion couplings to electrons, photons, and nucleons. This surface is inscribed in the cuboid defined by $g_{ae} < 3.7 \times 10^{-12}$, $g_{ae}g_{an}^{eff} < 4.6 \times 10^{-18}$, and $g_{ae}g_{a\gamma} < 7.6\times10^{-22}~{GeV}^{-1}$, and excludes either $g_{ae}=0$ or $g_{ae}g_{a\gamma}=g_{ae}g_{an}^{eff}=0$. The neutrino magnetic moment signal is similarly favored over background at 3.2$\sigma$ and a confidence interval of $\mu_{\nu} \in (1.4,~2.9)\times10^{-11}\mu_B$ (90% C.L.) is reported. Both results are in strong tension with stellar constraints. The excess can also be explained by $\beta$ decays of tritium at 3.2$\sigma$ significance with a corresponding tritium concentration in xenon of $(6.2 \pm 2.0) \times 10^{-25}$ mol/mol. Such a trace amount can be neither confirmed nor excluded with current knowledge of production and reduction mechanisms. The significances of the solar axion and neutrino magnetic moment hypotheses are decreased to 2.1$\sigma$ and 0.9$\sigma$, respectively, if an unconstrained tritium component is included in the fitting. With respect to bosonic dark matter, the excess favors a monoenergetic peak at ($2.3 \pm 0.2$) keV (68% C.L.) with a 3.0$\sigma$ global (4.0$\sigma$ local) significance over background. This analysis sets the most restrictive direct constraints to date on pseudoscalar and vector bosonic dark matter for most masses between 1 and 210 keV/c$^2$.

We detail the sensitivity of the liquid xenon (LXe) DARWIN observatory to solar neutrinos via elastic electron scattering. We find that DARWIN will have the potential to measure the fluxes of five solar neutrino components: $pp$, $^7$Be, $^{13}$N, $^{15}$O and $pep$. The precision of the $^{13}$N, $^{15}$O and $pep$ components is hindered by the double-beta decay of $^{136}$Xe and, thus, would benefit from a depleted target. A high-statistics observation of $pp$ neutrinos would allow us to infer the values of the weak mixing angle, $\sin^2\theta_w$, and the electron-type neutrino survival probability, $P_e$, in the electron recoil energy region from a few keV up to 200 keV for the first time, with relative precision of 5% and 4%, respectively, at an exposure of 300 ty. An observation of $pp$ and $^7$Be neutrinos would constrain the neutrino-inferred solar luminosity down to 0.2%. A combination of all flux measurements would distinguish between the high (GS98) and low metallicity (AGS09) solar models with 2.1-2.5$\sigma$ significance, independent of external measurements from other experiments or a measurement of $^8$B neutrinos through coherent elastic neutrino-nucleus scattering in DARWIN. Finally, we demonstrate that with a depleted target DARWIN may be sensitive to the neutrino capture process of $^{131}$Xe.

Gamma-ray bursts (GRBs) might be powered by a black hole (BH) hyperaccretion systems via the Blandford-Znajek (BZ) mechanism or neutrino annihilation from neutrino-dominated accretion flows (NDAFs). Magnetic coupling (MC) between the inner disc and BH can transfer angular momentum and energy from the fast-rotating BH to the disc. The neutrino luminosity and neutrino annihilation luminosity are both efficiently enhanced by the MC process. In this paper, we study the structure, luminosity, MeV neutrinos, and gravitational waves (GWs) of magnetized NDAFs (MNDAFs) under the assumption that both the BZ and MC mechanisms are present. The results indict that the BZ mechanism will compete with the neutrino annihilation luminosity to trigger jets under the different partitions of the two magnetic mechanisms. The typical neutrino luminosity and annihilation luminosity of MNDAFs are definitely higher than those of NDAFs. The typical peak energy of neutrino spectra of MNDAFs is higher than that of NDAFs, but similar to those of core-collapse supernovae. Moreover, if the MC process is dominant, then the GWs originating from the anisotropic neutrino emission will be stronger particularly for discs with high accretion rates.

The major scientific goals of DArk Matter Particle Explorer (DAMPE) are to study cosmic-ray electrons (including positrons) and gamma rays from 5 GeV to 10 TeV and nuclei from Z = 1 to 26 up to 100 TeV. The deposited energy measured by the Bismuth Germanate Oxide (BGO) calorimeter of DAMPE is affected by fluorescence attenuation in BGO crystals that are 600 mm long. In this work, an in-orbit attenuation calibration method is reported, and energy correction of the sensitive detector unit of the BGO calorimeter is also presented.

Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{\beta\beta}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $\sigma/\mu$ is as low as (0.79$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{\beta\beta}$. We achieve this by using a signal correction method to rectify the saturation effects of the signal readout system. The very good energy resolution from keV to MeV energies demonstrated in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events.

Direct dark matter detection experiments based on a liquid xenon target are leading the search for dark matter particles with masses above $\sim$ 5 GeV/c$^2$, but have limited sensitivity to lighter masses because of the small momentum transfer in dark matter-nucleus elastic scattering. However, there is an irreducible contribution from inelastic processes accompanying the elastic scattering, which leads to the excitation and ionization of the recoiling atom (the Migdal effect) or the emission of a Bremsstrahlung photon. In this letter, we report on a probe of low-mass dark matter with masses down to about 85 MeV/c$^2$ by looking for electronic recoils induced by the Migdal effect and Bremsstrahlung, using data from the XENON1T experiment. Besides the approach of detecting both scintillation and ionization signals, we exploit an approach that uses ionization signals only, which allows for a lower detection threshold. This analysis significantly enhances the sensitivity of XENON1T to light dark matter previously beyond its reach.

We report constraints on light dark matter (DM) models using ionization signals in the XENON1T experiment. We mitigate backgrounds with strong event selections, rather than requiring a scintillation signal, leaving an effective exposure of $(22 \pm 3)$ tonne-days. Above $\sim\!0.4\,\mathrm{keV}_\mathrm{ee}$, we observe $<1 \, \text{event}/(\text{tonne} \times \text{day} \times \text{keV}_\text{ee})$, which is more than one thousand times lower than in similar searches with other detectors. Despite observing a higher rate at lower energies, no DM or CEvNS detection may be claimed because we cannot model all of our backgrounds. We thus exclude new regions in the parameter spaces for DM-nucleus scattering for DM masses $m_\chi$ within $3-6\,\mathrm{GeV}/\mathrm{c}^2$, DM-electron scattering for $m_\chi > 30\,\mathrm{MeV}/\mathrm{c}^2$, and absorption of dark photons and axion-like particles for $m_\chi$ within $0.186 - 1 \, \mathrm{keV}/\mathrm{c}^2$.

As progenitors of gamma-ray bursts (GRBs), core collapse of massive stars and coalescence of compact object binaries are believed to be powerful sources of gravitational waves (GWs). In the collapsar scenario, a rotating stellar-mass black hole (BH) surrounded by a hyperaccretion disk might be running in the center of a massive collapsar, which is one of the plausible central engines of long GRBs. Such a BH hyperaccretion disk would be in a state of a neutrino-dominated accretion flow (NDAF) at the initial stage of the accretion process; meanwhile, the jets attempt to break out from the envelope and circumstellar medium to power GRBs. In addition to collapsars, the BH hyperaccretion systems are important sources of neutrinos and GWs. In this paper, we investigated the GW emission generated by the anisotropic neutrino emission from NDAFs in the collapsar scenarios. As the results indicate, the typical frequency of GWs is $\sim$ 1-100 Hz, and the masses and metallicities of the progenitor stars have slight effects on the GW strains. The GWs from NDAFs might be detected by operational or planned detectors at the distance of 10 kpc. Moreover, comparisons of the detectable GWs from collapsars, NDAFs, and GRB jets (internal shocks) are displayed. By combining the electromagnetic counterparts, neutrinos, and GWs, one may constrain the characteristics of collapsars and central BH accretion systems.

Precise measurements of the energy spectra of cosmic rays (CRs) show various kinds of features deviating from single power-laws, which give very interesting and important implications on their origin and propagation. Previous measurements from a few balloon and space experiments indicate the existence of spectral softenings around 10 TV for protons (and probably also for Helium nuclei). Very recently, the DArk Matter Particle Explorer (DAMPE) measurement about the proton spectrum clearly reveals such a softening with a high significance. Here we study the implications of these new measurements, as well as the groundbased indirect measurements, on the origin of CRs. We find that a single component of CRs fails to fit the spectral softening and the air shower experiment data simultaneously. In the framework of multiple components, we discuss two possible scenarios, the multiple source population scenario and the background plus nearby source scenario. Both scenarios give reasonable fits to the wide-band data from TeV to 100 PeV energies. Considering the anisotropy observations, the nearby source model is favored.

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.

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.

As the plausible central engine of gamma-ray bursts (GRBs), a black hole (BH) hyperaccretion disk should be in a state of neutrino-dominated accretion flow (NDAF) if the accretion rate is larger than the ignition rate of an NDAF. A rotating stellar-mass BH surrounded by a hyperaccretion disk might be born in the center of a massive core collapsar. In the initial hundreds of seconds of the accretion process, the mass supply rate of the massive progenitor is generally higher than the ignition accretion rate, but the jets are generally choked in the envelope. Thus, neutrinos can be emitted from the center of a core collapsar. In this paper, we study the effects of the masses and metallicities of progenitor stars on the time-integrated spectra of electron neutrinos from NDAFs. The peak energies of the calculated spectra are approximately 10-20 MeV. The mass of a collapsar has little influence on the neutrino spectrum, and a low metallicity is beneficial to the production of low-energy ($\lesssim$ 1 MeV) neutrinos. We also investigate the differences in the electron neutrino spectra between NDAFs and proto-neutron stars. Combining with the electromagnetic counterparts and multi-messenger astronomy, one may verify the possible remnants of the core collapse of massive stars with future neutrino detectors.

We report the first experimental results on spin-dependent elastic weakly interacting massive particle (WIMP) nucleon scattering from the XENON1T dark matter search experiment. The analysis uses the full ton year exposure of XENON1T to constrain the spin-dependent proton-only and neutron-only cases. No significant signal excess is observed, and a profile likelihood ratio analysis is used to set exclusion limits on the WIMP-nucleon interactions. This includes the most stringent constraint to date on the WIMP-neutron cross section, with a minimum of $6.3\times10^{-42}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level. The results are compared with those from collider searches and used to exclude new parameter space in an isoscalar theory with an axial-vector mediator.

Understanding the sources of lunar water is crucial for studying the history of lunar evolution, and also the solar wind interaction with the Moon and other airless bodies. Recent observations revealed lunar hydration is very likely a surficial dynamic process driven by solar wind. Solar wind is shielded over a period of 3-5 days as the Moon passes through the Earth's magnetosphere, during which a significant loss of hydration is expected from previous works.Here we study lunar hydration inside the magnetosphere using orbital spectral data, which unexpectedly found that the polar surficial OH/H2O abundance remains at the same level when in the solar wind and in the magnetosphere. We suggest that particles from the magnetosphere (Earth wind, naturally different from solar wind) contribute to lunar hydration. From lunar orbital plasma observations, we find the existence of optimal energy ranges, other than 1 keV as previously thought, for surface hydration formation. These optimal energy ranges deduced from space observations may provide strong implications for laboratory experiments simulating lunar hydration processes.

We present first results on the scalar WIMP-pion coupling from 1 t$\times$yr of exposure with the XENON1T experiment. This interaction is generated when the WIMP couples to a virtual pion exchanged between the nucleons in a nucleus. In contrast to most non-relativistic operators, these pion-exchange currents can be coherently enhanced by the total number of nucleons, and therefore may dominate in scenarios where spin-independent WIMP-nucleon interactions are suppressed. Moreover, for natural values of the couplings, they dominate over the spin-dependent channel due to their coherence in the nucleus. Using the signal model of this new WIMP-pion channel, no significant excess is found, leading to an upper limit cross section of $6.4\times10^{-46}$ cm$^2$ (90 % confidence level) at 30 GeV/c$^2$ WIMP mass.

The DArk Matter Particle Explorer (DAMPE) developed in China was designed to search for evidence of dark matter particles by observing primary cosmic rays and gamma rays in the energy range from 5 GeV to 10 TeV. Since its launch in December 2015, a large quantity of data has been recorded. With the data set acquired during more than a year of operation in space, a precise time-dependent calibration for the energy measured by the BGO ECAL has been developed. In this report, the instrumentation and development of the BGO Electromagnetic Calorimeter (BGO ECAL) are briefly described. The calibration on orbit, including that of the pedestal, attenuation length, minimum ionizing particle peak, and dynode ratio, is discussed, and additional details about the calibration methods and performance in space are presented.

We report on a search for Weakly Interacting Massive Particles (WIMPs) using 278.8 days of data collected with the XENON1T experiment at LNGS. XENON1T utilizes a liquid xenon time projection chamber with a fiducial mass of $(1.30 \pm 0.01)$ t, resulting in a 1.0 t$\times$yr exposure. The energy region of interest, [1.4, 10.6] $\mathrm{keV_{ee}}$ ([4.9, 40.9] $\mathrm{keV_{nr}}$), exhibits an ultra-low electron recoil background rate of $(82\substack{+5 \\ -3}\textrm{ (sys)}\pm3\textrm{ (stat)})$ events/$(\mathrm{t}\times\mathrm{yr}\times\mathrm{keV_{ee}})$. No significant excess over background is found and a profile likelihood analysis parameterized in spatial and energy dimensions excludes new parameter space for the WIMP-nucleon spin-independent elastic scatter cross-section for WIMP masses above 6 GeV/c${}^2$, with a minimum of $4.1\times10^{-47}$ cm$^2$ at 30 GeV/c${}^2$ and 90% confidence level.

In this paper, we describe the XENON100 data analyses used to assess the target-intrinsic background sources radon ($^{222}$Rn), thoron ($^{220}$Rn) and krypton ($^{85}$Kr). We detail the event selections of high-energy alpha particles and decay-specific delayed coincidences. We derive distributions of the individual radionuclides inside the detector and quantify their abundances during the main three science runs of the experiment over a period of $\sim$ 4 years, from January 2010 to January 2014. We compare our results to external measurements of radon emanation and krypton concentrations where we find good agreement. We report an observed reduction in concentrations of radon daughters that we attribute to the plating-out of charged ions on the negatively biased cathode.

We report on the response of liquid xenon to low energy electronic recoils below 15 keV from beta decays of tritium at drift fields of 92 V/cm, 154 V/cm and 366 V/cm using the XENON100 detector. A data-to-simulation fitting method based on Markov Chain Monte Carlo is used to extract the photon yields and recombination fluctuations from the experimental data. The photon yields measured at the two lower fields are in agreement with those from literature; additional measurements at a higher field of 366 V/cm are presented. The electronic and nuclear recoil discrimination as well as its dependence on the drift field and photon detection efficiency are investigated at these low energies. The results provide new measurements in the energy region of interest for dark matter searches using liquid xenon.

The DArk Matter Particle Explorer (DAMPE), a high energy cosmic ray and $\gamma$-ray detector in space, has recently reported the new measurement of the total electron plus positron flux between 25 GeV and 4.6 TeV. A spectral softening at $\sim0.9$ TeV and a tentative peak at $\sim1.4$ TeV have been reported. We study the physical implications of the DAMPE data in this work. The presence of the spectral break significantly tightens the constraints on the model parameters to explain the electron/positron excesses. The spectral softening can either be explained by the maximum acceleration limits of electrons by astrophysical sources, or a breakdown of the common assumption of continuous distribution of electron sources at TeV energies in space and time. The tentive peak at $\sim1.4$ TeV implies local sources of electrons/positrons with quasi-monochromatic injection spectrum. We find that the cold, ultra-relativistic $e^+e^-$ winds from pulsars may give rise to such a structure. The pulsar is requird to be middle-aged, relatively slowly-rotated, mildly magnetized, and isolated in a density cavity. The annihilation of DM particles ($m_{\chi}\sim1.5$ TeV) into $e^+e^-$ pairs in a nearby clump or an over-density region may also explain the data. In the DM scenario, the inferred clump mass (or density enhancement) is about $10^7-10^8$ M$_\odot$ (or $17-35$ times of the canonical local density) assuming a thermal production cross section, which is relatively extreme compared with the expectation from numerical simulations. A moderate enhancement of the annihilation cross section via, e.g., the Sommerfeld mechanism or non-thermal production, is thus needed.

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.

We report the first dark matter search results from XENON1T, a $\sim$2000-kg-target-mass dual-phase (liquid-gas) xenon time projection chamber in operation at the Laboratori Nazionali del Gran Sasso in Italy and the first ton-scale detector of this kind. The blinded search used 34.2 live days of data acquired between November 2016 and January 2017. Inside the (1042$\pm$12) kg fiducial mass and in the [5, 40] $\mathrm{keV}_{\mathrm{nr}}$ energy range of interest for WIMP dark matter searches, the electronic recoil background was $(1.93 \pm 0.25) \times 10^{-4}$ events/(kg $\times$ day $\times \mathrm{keV}_{\mathrm{ee}}$), the lowest ever achieved in a dark matter detector. A profile likelihood analysis shows that the data is consistent with the background-only hypothesis. We derive the most stringent exclusion limits on the spin-independent WIMP-nucleon interaction cross section for WIMP masses above 10 GeV/c${}^2$, with a minimum of 7.7 $\times 10^{-47}$ cm${}^2$ for 35-GeV/c${}^2$ WIMPs at 90% confidence level.

We present the first search for dark matter-induced delayed coincidence signals in a dual-phase xenon time projection chamber, using the 224.6 live days of the XENON100 science run II. This very distinct signature is predicted in the framework of magnetic inelastic dark matter which has been proposed to reconcile the modulation signal reported by the DAMA/LIBRA collaboration with the null results from other direct detection experiments. No candidate event has been found in the region of interest and upper limits on the WIMP's magnetic dipole moment are derived. The scenarios proposed to explain the DAMA/LIBRA modulation signal by magnetic inelastic dark matter interactions of WIMPs with masses of 58.0 GeV/c$^2$ and 122.7 GeV/c$^2$ are excluded at 3.3 $\sigma$ and 9.3 $\sigma$, respectively.

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.

In this paper, we presented a study of temperature effects on BGO calorimeters using proton MIP's collected in the first year operation of DAMPE. By directly comparing MIP calibration constants used by DAMPE data production pipe line, we found an experimental relation between temperature and signal amplitudes of each BGO bar: a general deviation of -1.162%/$^{\circ}$C,and -0.47%/$^{\circ}$C to -1.60%/$^{\circ}$C statistically for each detector element. During 2016, DAMPE's temperature changed by about 7 degrees due to solar elevation angle and the corresponding energy scale bias is about 8%. By frequent MIP calibration operation, this kind of bias is eliminated to an acceptable value.

We present results of searches for vector and pseudo-scalar bosonic super-WIMPs, which are dark matter candidates with masses at the keV-scale, with the XENON100 experiment. XENON100 is a dual-phase xenon time projection chamber operated at the Laboratori Nazionali del Gran Sasso. A profile likelihood analysis of data with an exposure of 224.6 live days $\times$ 34\,kg showed no evidence for a signal above the expected background. We thus obtain new and stringent upper limits in the $(8-125)$\,keV/c$^2$ mass range, excluding couplings to electrons with coupling constants of $g_{ae} > 3\times10^{-13}$ for pseudo-scalar and $\alpha'/\alpha > 2\times10^{-28}$ for vector super-WIMPs, respectively. These limits are derived under the assumption that super-WIMPs constitute all of the dark matter in our galaxy.

The XENON1T experiment at the Laboratori Nazionali del Gran Sasso (LNGS) is the first WIMP dark matter detector operating with a liquid xenon target mass above the ton-scale. Out of its 3.2t liquid xenon inventory, 2.0t constitute the active target of the dual-phase time projection chamber. The scintillation and ionization signals from particle interactions are detected with low-background photomultipliers. This article describes the XENON1T instrument and its subsystems as well as strategies to achieve an unprecedented low background level. First results on the detector response and the performance of the subsystems are also presented.

A software system has been developed for the DArk Matter Particle Explorer (DAMPE) mission, a satellite-based experiment. The DAMPE software is mainly written in C++ and steered using Python script. This article presents an overview of the DAMPE offline software, including the major architecture design and specific implementation for simulation, calibration and reconstruction. The whole system has been successfully applied to DAMPE data analysis, based on which some results from simulation and beam test experiments are obtained and presented.

We report on WIMP search results in the XENON100 detector using a non-relativistic effective field theory approach. The data from science run II (34 kg $\times$ 224.6 live days) was re-analyzed, with an increased recoil energy interval compared to previous analyses, ranging from $(6.6 - 240)~\mathrm{keV_\mathrm{nr}}$. The data is found to be compatible with the background-only hypothesis. We present 90% confidence level exclusion limits on the coupling constants of WIMP-nucleon effective operators using a binned profile likelihood method. We also consider the case of inelastic WIMP scattering, where incident WIMPs may up-scatter to a higher mass state, and set exclusion limits on this model as well.

We present the first constraints on the spin-dependent, inelastic scattering cross section of Weakly Interacting Massive Particles (WIMPs) on nucleons from XENON100 data with an exposure of 7.64$\times$10$^3$\,kg\,day. XENON100 is a dual-phase xenon time projection chamber with 62\,kg of active mass, operated at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy and designed to search for nuclear recoils from WIMP-nucleus interactions. Here we explore inelastic scattering, where a transition to a low-lying excited nuclear state of $^{129}$Xe is induced. The experimental signature is a nuclear recoil observed together with the prompt de-excitation photon. We see no evidence for such inelastic WIMP-$^{129}$Xe interactions. A profile likelihood analysis allows us to set a 90\% C.L. upper limit on the inelastic, spin-dependent WIMP-nucleon cross section of $3.3 \times 10^{-38}$\,cm$^{2}$ at 100\,GeV/c$^2$. This is the most constraining result to date, and sets the pathway for an analysis of this interaction channel in upcoming, larger dual-phase xenon detectors.

The Hamamatsu R11410-21 photomultiplier tube is the photodetector of choice for the XENON1T dual-phase time projection chamber. The device has been optimized for a very low intrinsic radioactivity, a high quantum efficiency and a high sensitivity to single photon detection. A total of 248 tubes are currently operated in XENON1T, selected out of 321 tested units. In this article the procedures implemented to evaluate the large number of tubes prior to their installation in XENON1T are described. The parameter distributions for all tested tubes are shown, with an emphasis on those selected for XENON1T, of which the impact on the detector performance is discussed. All photomultipliers have been tested in a nitrogen atmosphere at cryogenic temperatures, with a subset of the tubes being tested in gaseous and liquid xenon, simulating their operating conditions in the dark matter detector. The performance and evaluation of the tubes in the different environments is reported and the criteria for rejection of PMTs are outlined and quantified.

We report on WIMP search results of the XENON100 experiment, combining three runs summing up to 477 live days from January 2010 to January 2014. Data from the first two runs were already published. A blind analysis was applied to the last run recorded between April 2013 and January 2014 prior to combining the results. The ultra-low electromagnetic background of the experiment, ~$5 \times 10^{-3}$ events/(keV$_{\mathrm{ee}}\times$kg$\times$day) before electronic recoil rejection, together with the increased exposure of 48 kg $\times$ yr improves the sensitivity. A profile likelihood analysis using an energy range of (6.6 - 43.3) keV$_{\mathrm{nr}}$ sets a limit on the elastic, spin-independent WIMP-nucleon scattering cross section for WIMP masses above 8 GeV/$c^2$, with a minimum of 1.1 $\times 10^{-45}$ cm$^2$ at 50 GeV/$c^2$ and 90% confidence level. We also report updated constraints on the elastic, spin-dependent WIMP-nucleon cross sections obtained with the same data. We set upper limits on the WIMP-neutron (proton) cross section with a minimum of 2.0 $\times 10^{-40}$ cm$^2$ (52$\times 10^{-40}$ cm$^2$) at a WIMP mass of 50 GeV/$c^2$, at 90% confidence level.

We report on a search for electronic recoil event rate modulation signatures in the XENON100 data accumulated over a period of 4 years, from January 2010 to January 2014. A profile likelihood method, which incorporates the stability of the XENON100 detector and the known electronic recoil background model, is used to quantify the significance of periodicity in the time distribution of events. There is a weak modulation signature at a period of $431^{+16}_{-14}$ days in the low energy region of $(2.0-5.8)$ keV in the single scatter event sample, with a global significance of $1.9\,\sigma$, however no other more significant modulation is observed. The expected annual modulation of a dark matter signal is not compatible with this result. Single scatter events in the low energy region are thus used to exclude the DAMA/LIBRA annual modulation as being due to dark matter electron interactions via axial vector coupling at $5.7\,\sigma$.

We perform a low-mass dark matter search using an exposure of 30\,kg$\times$yr with the XENON100 detector. By dropping the requirement of a scintillation signal and using only the ionization signal to determine the interaction energy, we lowered the energy threshold for detection to 0.7\,keV for nuclear recoils. No dark matter detection can be claimed because a complete background model cannot be constructed without a primary scintillation signal. Instead, we compute an upper limit on the WIMP-nucleon scattering cross section under the assumption that every event passing our selection criteria could be a signal event. Using an energy interval from 0.7\,keV to 9.1\,keV, we derive a limit on the spin-independent WIMP-nucleon cross section that excludes WIMPs with a mass of 6\,GeV/$c^2$ above $1.4 \times 10^{-41}$\,cm$^2$ at 90\% confidence level.

DARk matter WImp search with liquid xenoN (DARWIN) will be an experiment for the direct detection of dark matter using a multi-ton liquid xenon time projection chamber at its core. Its primary goal will be to explore the experimentally accessible parameter space for Weakly Interacting Massive Particles (WIMPs) in a wide mass-range, until neutrino interactions with the target become an irreducible background. The prompt scintillation light and the charge signals induced by particle interactions in the xenon will be observed by VUV sensitive, ultra-low background photosensors. Besides its excellent sensitivity to WIMPs above a mass of 5 GeV/c2, such a detector with its large mass, low-energy threshold and ultra-low background level will also be sensitive to other rare interactions. It will search for solar axions, galactic axion-like particles and the neutrinoless double-beta decay of 136-Xe, as well as measure the low-energy solar neutrino flux with <1% precision, observe coherent neutrino-nucleus interactions, and detect galactic supernovae. We present the concept of the DARWIN detector and discuss its physics reach, the main sources of backgrounds and the ongoing detector design and R&D efforts.

We present monitoring observations of the active T Tauri star RW Aur, from 2010 October to 2015 January, using optical high-resolution (R>10000) spectroscopy with CFHT-ESPaDOnS. Optical photometry in the literature shows bright, stable fluxes over most of this period, with lower fluxes (by 2-3 mag.) in 2010 and 2014. In the bright period our spectra show clear photospheric absorption, complicated variation in the Ca II 8542 A emission}profile shapes, and a large variation in redshifted absorption in the O I 7772 and 8446 A and He I 5876 A lines, suggesting unstable mass accretion during this period. In contrast, these line profiles are relatively uniform during the faint periods, suggesting stable mass accretion. During the faint periods the photospheric absorption lines are absent or marginal, and the averaged Li I profile shows redshifted absorption due to an inflow. We discuss (1) occultation by circumstellar material or a companion and (2) changes in the activity of mass accretion to explain the above results, together with near-infrared and X-ray observations from 2011-2015. Neither scenario can simply explain all the observed trends, and more theoretical work is needed to further investigate their feasibilities.

The XENON1T experiment is currently in the commissioning phase at the Laboratori Nazionali del Gran Sasso, Italy. In this article we study the experiment's expected sensitivity to the spin-independent WIMP-nucleon interaction cross section, based on Monte Carlo predictions of the electronic and nuclear recoil backgrounds. The total electronic recoil background in $1$ tonne fiducial volume and ($1$, $12$) keV electronic recoil equivalent energy region, before applying any selection to discriminate between electronic and nuclear recoils, is $(1.80 \pm 0.15) \cdot 10^{-4}$ ($\rm{kg} \cdot day \cdot keV)^{-1}$, mainly due to the decay of $^{222}\rm{Rn}$ daughters inside the xenon target. The nuclear recoil background in the corresponding nuclear recoil equivalent energy region ($4$, $50$) keV, is composed of $(0.6 \pm 0.1)$ ($\rm{t} \cdot y)^{-1}$ from radiogenic neutrons, $(1.8 \pm 0.3) \cdot 10^{-2}$ ($\rm{t} \cdot y)^{-1}$ from coherent scattering of neutrinos, and less than $0.01$ ($\rm{t} \cdot y)^{-1}$ from muon-induced neutrons. The sensitivity of XENON1T is calculated with the Profile Likelihood Ratio method, after converting the deposited energy of electronic and nuclear recoils into the scintillation and ionization signals seen in the detector. We take into account the systematic uncertainties on the photon and electron emission model, and on the estimation of the backgrounds, treated as nuisance parameters. The main contribution comes from the relative scintillation efficiency $\mathcal{L}_\mathrm{eff}$, which affects both the signal from WIMPs and the nuclear recoil backgrounds. After a $2$ y measurement in $1$ t fiducial volume, the sensitivity reaches a minimum cross section of $1.6 \cdot 10^{-47}$ cm$^2$ at m$_\chi$=$50$ GeV/$c^2$.

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.

There exist cislunar and trans-lunar libration points near the Moon, which are referred as the LL1 and LL2 points respectively and can generate the different types of low-energy trajectories transferring from Earth to Moon. The time-dependent analytic model including the gravitational forces from the Sun, Earth and Moon is employed to investigate the energy-minimal and practical transfer trajectories. However, different from the circular restricted three-body problem, the equivalent gravitational equilibria are defined according to the geometry of instantaneous Hill's boundary due to the gravitational perturbation from the Sun. The relationship between the altitudes of periapsis and eccentricities is achieved from the Poincar\'e mapping for all the lunar captured trajectories, which presents the statistical feature of the fuel cost and captured orbital elements rather than generating a specified Moon-captured segment. The minimum energy required by the captured trajectory on a lunar circular orbit is deduced in the spatial bi-circular model. It is presented that the asymptotical behaviors of invariant manifolds approaching to/from the libration points or halo orbits are destroyed by the solar perturbation. In fact, the energy-minimal cislunar transfer trajectory is acquired by transiting LL1 point, while the energy-minimal trans-lunar transfer trajectory is obtained by transiting LL2 point. Finally, the transfer opportunities for the practical trajectories escaped from the Earth and captured by the Moon are yielded by transiting halo orbits near LL1 and LL2 points, which can be used to generate the whole trajectories.

In this paper we investigate the space and velocity distributions of old neutron stars (aged 109 to 1010 yr) in our Galaxy. Galactic old Neutron Stars (NSs) population fills a torus-like area extending to a few tens kiloparsecs above the galactic plane. The initial velocity distribution of NSs is not well known, in this work we adopt a three component initial distribution, as given by the contribution of kick velocities, circular velocities and Maxwellian velocities. For the spatial initial distribution we use a Gamma function. We then use Monte Carlo simulations to follow the evolution of the NSs under the influence of the Paczy{\P}nski Galactic gravitational potential. Our calculations show that NS orbits have a very large Galactic radial expansion and that their radial distribution peak is quite close to their progenitors' one. We also study the NS vertical distribution and find that it can well be described by a double exponential low. Finally, we investigate the correlation of the vertical and radial distribution and study the radial dependence of scale-heights.

The propagation properties of coronal mass ejections (CMEs) are crucial to predict its geomagnetic effect. A newly developed three dimensional (3D) mask fitting reconstruction method using coronagraph images from three viewpoints has been described and applied to the CME ejected on August 7, 2010. The CME's 3D localisation, real shape and morphological evolution are presented. Due to its interaction with the ambient solar wind, the morphology of this CME changed significantly in the early phase of evolution. Two hours after its initiation, it was expanding almost self-similarly. CME's 3D localisation is quite helpful to link remote sensing observations to in situ measurements. The investigated CME was propagating to Venus with its flank just touching STEREO B. Its corresponding ICME in the interplanetary space shows a possible signature of a magnetic cloud with a preceding shock in VEX observations, while from STEREO B only a shock is observed. We have calculated three principle axes for the reconstructed 3D CME cloud. The orientation of the major axis is in general consistent with the orientation of a filament (polarity inversion line) observed by SDO/AIA and SDO/HMI. The flux rope axis derived by the MVA analysis from VEX indicates a radial-directed axis orientation. It might be that locally only the leg of the flux rope passed through VEX. The height and speed profiles from the Sun to Venus are obtained. We find that the CME speed possibly had been adjusted to the speed of the ambient solar wind flow after leaving COR2 field of view and before arriving Venus. A southward deflection of the CME from the source region is found from the trajectory of the CME geometric center. We attribute it to the influence of the coronal hole where the fast solar wind emanated from.

We make statistics about the number distribution of normal radio pulsars according to the pulsar spin period. And find the physical meaning of the traditional statistical method about this problem is very limited. According to the statistical method proposed by us, the distribution rule is exactly the exponentially decay with the spin period increasing, especially important, which can be interpreted as normal radio pulsar number decays as the radiative elements decay. It discloses that though belonging to two different extreme scales of the universe and having different inner active mechanisms, the normal radio pulsars and radiative elements accord with the same varying rule. And like the half decay period of the radiative elements we can also introduce the half decay spin period of the normal radio pulsars.

Based the Alfven wave oscillation model (AWOM) and relativistic precession model (RPM) for twin kHz QPOs, we estimate the emission positions of most detected kHz QPOs to be at r=18+-3 km (R/15km) except Cir X-1 at r = 30\+-5 km (R/15km). For the proposed Keplerian frequency as an upper limit to kHz QPO, the spin effects in Kerr Spacetime are discussed, which have about a 5% (2%) modification for that of the Schwarzchild case for the spin frequency of 1000 (400) Hz.The application to the four typical QPO sources, Cir X-1, Sco X-1, SAX J1808.4-3658 and XTE 1807-294, is mentioned.

With the possible detection of the fastest spinning nuclear-powered pulsar XTE J1739-285 of frequency 1122 Hz (0.8913 ms), it arouses us to constrain the mass and radius of its central compact object and to imply the stellar matter compositions: neutrons or quarks. Spun-up by the accreting materials to such a high rotating speed, the compact star should have either a small radius or short innermost stable circular orbit. By the empirical relation between the upper kHz quasi-periodic oscillation frequency and star spin frequency, a strong constraint on mass and radius is obtained as 1.51 solar masses and 10.9 km, which excludes most equations of states (EOSs) of normal neutrons and strongly hints the star promisingly to be a strange quark star. Furthermore, the star magnetic field is estimated to be about $4\times10^{7} (G) < B < 10^{9} (G) $, which reconciles with those of millisecond radio pulsars, revealing the clues of the evolution linkage of two types of astrophysical objects.

The late-time evolution behavior of the autonomous system in the SO(1,1) dark energy model with power-law potential is studied. Big Rip may be a critical point of the autonomous system. This means that such a Big Rip may be considered as the middle state between the expanding and contracting phases of phantom universe. This result is also valid for some special interactions between matter and dark energy.

We report the discovery of a massive SCUBA core with evidence of inflow and outflow motions. This core is detected by SCUBA at both 450 and 850 micro meter. Barely resolved by the telescope beam at 450 micro meter, it has a size of 10 arcsec, corresponding to 0.28 pc at a distance of 5.7 kpc. The dust temperature is estimated to be less than 29 K, the total mass is 820 solar mass and the average density is 1.1e-6 per cube-cm in a region with a radius of 5 arcsec. Follow-up spectral line observations, including HCN (3-2), HCO+ (3-2), H13CO+ (3-2) and C17O (2-1) reveal a typical blue profile which indicates that this core is collapsing. The CO (3-2) line profile is as broad as 38 km/s, indicating outflow motions in this region. This core is approximately 1.5 pc away from the known HII region G25.4NW, but there are no obvious radio, IRAS, MSX or Spitzer sources associated with it. We suggest that this core is at a very early stage of massive star or cluster formation.

In order to depict the transition from deceleration to acceleration expansion of the universe we use a power-law expansion scale factor, $a\sim t^{n_0+bt^m}$, with $n_0$, $b$ and $m$ three parameters determined by $H_0$, $q_0$ and $z_T$. For the spatially flat, isotropic and homogeneous universe, such a scale factor leads to the results that the dark energy density is slowly changing currently, and predicts the equation of state $w_X$ changes from $w_X>-1$ to $w_X<-1$.

Using the SNe Ia data we determine the three parameters in the power-law expanding universe model with time-dependent power \cite{W}. Inputting $H_0$ and $t_0$, then we find the $\dot{a}-t$ evolution curve with $m=5.0$ and $q_0=-0.90$ can fit very well to that from SNe observation data. The model predicts the transition redshift $z\simeq0.38$. The dark energy deduced from this model have phantom property but the universe doesn't encounter the Big Rip singularity. Assuming that this model with the three parameters is valid for the future universe, then we predict that the total energy density of the universe is decreasing and will soon reach its minimum.

A statistical study of the properties of molecular outflows is performed based on an up-to-date sample. 391 outflows were identified in published articles or preprints before February 28, 2003. The parameters of position, morphology, mass, energy, outflow dynamics and central source luminosity are presented for each outflow source. Outflow lobe polarity is known for all the sources, and 84% are found to be bipolar. The sources are divided into low mass and high mass groups according to either the available bolometric luminosity of the central source or the outflow mass.Energetic mass ejection may be a common aspect of the formation of high mass as well as low mass stars. Outflow masses are correlated strongly with bolometric luminosity of the center sources, which was obtained for the first time. There are also correlations between the central source luminosity and the parameters of mechanical luminosity and the thrust or force necessary to drive the outflow. The results show that flow mass, momentum and energy depend on the nature of the central source. Despite their similarity, there are differences between the high mass and low mass outflows on collimation, flow masses, mass-dynamical time relation and associated objects such as HH objects and water masers.Sources with characteristics of collapse or infall comprise 12% of the entire outflow sample. The spatial distribution of the outflow sources in the Galaxy is presented and the local occurrence rate is compared with the stellar birth rate.

In this paper, we show the possibility that an extended complex scalar field can be considered as an extended complex quintessence. For this model, we derive the basic equations with a parameter $\eta $ which govern the evolution of the Universe. Our model may contain a complex quintessence model for $\eta =-1$ and a real quitessence model for $\eta =0$.