The Doppler tracking data of the Chang'e 3 lunar mission is used to constrain the stochastic background of gravitational wave in cosmology within the 1 mHz to 0.05 Hz frequency band. Our result improves on the upper bound on the energy density of the stochastic background of gravitational wave in the 0.02 Hz to 0.05 Hz band obtained by the Apollo missions, with the improvement reaching almost one order of magnitude at around 0.05 Hz. Detailed noise analysis of the Doppler tracking data is also presented, with the prospect that these noise sources will be mitigated in future Chinese deep space missions. A feasibility study is also undertaken to understand the scientific capability of the Chang'e 4 mission, due to be launched in 2018, in relation to the stochastic gravitational wave background around 0.01 Hz. The study indicates that the upper bound on the energy density may be further improved by another order of magnitude from the Chang'e 3 mission, which will fill the gap in the frequency band from 0.02 Hz to 0.1 Hz in the foreseeable future.

We present a non-supersymmetric scenario in which the R-parity symmetry $R_P = (-1)^{3(B-L)+2s}$ arises as a result of spontaneous gauge symmetry breaking, leading to a viable Dirac fermion WIMP dark matter candidate. Direct detection in nuclear recoil experiments probes dark matter masses around $2-5$ TeV for $M_{Z^{\prime}} \sim 3-4$ TeV consistent with searches at the LHC, while lepton flavor violation rates and flavor changing neutral currents in neutral meson systems lie within reach of upcoming experiments.

In the centenary year of Einstein's General Theory of Relativity, this paper reviews the current status of gravitational wave astronomy across a spectrum which stretches from attohertz to kilohertz frequencies. Sect. 1 of this paper reviews the historical development of gravitational wave astronomy from Einstein's first prediction to our current understanding the spectrum. It is shown that detection of signals in the audio frequency spectrum can be expected very soon, and that a north-south pair of next generation detectors would provide large scientific benefits. Sect. 2 reviews the theory of gravitational waves and the principles of detection using laser interferometry. The state of the art Advanced LIGO detectors are then described. These detectors have a high chance of detecting the first events in the near future. Sect. 3 reviews the KAGRA detector currently under development in Japan, which will be the first laser interferometer detector to use cryogenic test masses. Sect. 4 of this paper reviews gravitational wave detection in the nanohertz frequency band using the technique of pulsar timing. Sect. 5 reviews the status of gravitational wave detection in the attohertz frequency band, detectable in the polarisation of the cosmic microwave background, and discusses the prospects for detection of primordial waves from the big bang. The techniques described in sects. 1-5 have already placed significant limits on the strength of gravitational wave sources. Sects. 6 and 7 review ambitious plans for future space based gravitational wave detectors in the millihertz frequency band. Sect. 6 presents a roadmap for development of space based gravitational wave detectors by China while sect. 7 discusses a key enabling technology for space interferometry known as time delay interferometry.

Part of a review paper entitled "Gravitational wave astronomy: the current status.", appeared in " Science China Physics, Mechanics & Astronomy 58.12 (2015): 1-41.

The present work reports on a feasibility study commissioned by the Chinese Academy of Sciences of China to explore various possible mission options to detect gravitational waves in space alternative to that of the eLISA/LISA mission concept. Based on the relative merits assigned to science and technological viability, a few representative mission options descoped from the ALIA mission are considered. A semi-analytic Monte Carlo simulation is carried out to understand the cosmic black hole merger histories starting from intermediate mass black holes at high redshift as well as the possible scientific merits of the mission options considered in probing the light seed black holes and their coevolution with galaxies in early Universe. The study indicates that, by choosing the armlength of the interferometer to be three million kilometers and shifting the sensitivity floor to around one-hundredth Hz, together with a very moderate improvement on the position noise budget, there are certain mission options capable of exploring light seed, intermediate mass black hole binaries at high redshift that are not readily accessible to eLISA/LISA, and yet the technological requirements seem to within reach in the next few decades for China.

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85{\deg} to +20{\deg}, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

A multipath mechanism similar to that used in Australia sixty years ago by the Sea-cliff Interferometer is shown to generate correlations between the periods of oscillations observed by two distant radio telescopes pointed to the Sun. The oscillations are the result of interferences between the direct wave detected in the main antenna lobe and its reflection on ground detected in a side lobe. A model is made of such oscillations in the case of two observatories located at equal longitudes and opposite tropical latitudes, respectively in Ha Noi (Viet Nam) and Learmonth (Australia), where similar radio telescopes are operated at 1.4 GHz. Simple specular reflection from ground is found to give a good description of the observed oscillations and to explain correlations that had been previously observed and for which no satisfactory interpretation, instrumental or other, had been found.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013

After a brief reminder of the mechanism of gravitational lensing of extended sources, the particular case of the host galaxy of QSO RXJ0911, a high redshift (z~2.8) quadruply imaged quasar, is explored. The non linearity of the problem, together with the proximity of the source to a cusp of the lens inner caustic, have important consequences on the dependence of the image appearance on the size and shape of the source. Their expected main features and their interpretation in terms of source extension and shape are investigated in a spirit of simplicity and in preparation for the analysis of high sensitivity and spatial resolution images that will soon be within reach with the completion of the Atacama Large Millimeter/submillimeter Array (ALMA). In particular, the information on source size carried by relative image brightness is discussed. Extension of the results to other types of quadruply imaged quasars is briefly considered.

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

The development of extensive air showers at extreme energies is studied using a simulation model much simpler and cruder, but also more transparent and flexible, than existing sophisticated codes. Evidence for its satisfactory performance is presented. As an illustration, shower elongation rates are evaluated in the $10^{18}$ to $10^{20}$ eV region and compared with recently published data. Lateral distribution functions of both muons and electrons/photons are also briefly discussed. Reliable results are obtained in the comparison between proton-induced and iron-induced showers.

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malarg\"ue and averaged monthly models, the utility of the GDAS data is shown.

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.

We present the results of an analysis of data recorded at the Pierre Auger Observatory in which we search for groups of directionally-aligned events (or `multiplets') which exhibit a correlation between arrival direction and the inverse of the energy. These signatures are expected from sets of events coming from the same source after having been deflected by intervening coherent magnetic fields. The observation of several events from the same source would open the possibility to accurately reconstruct the position of the source and also measure the integral of the component of the magnetic field orthogonal to the trajectory of the cosmic rays. We describe the largest multiplets found and compute the probability that they appeared by chance from an isotropic distribution. We find no statistically significant evidence for the presence of multiplets arising from magnetic deflections in the present data.