We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

Rapid response and short time latency are very important for Time Domain Astronomy, such as the observations of Gamma-ray Bursts (GRBs) and electromagnetic (EM) counterparts of gravitational waves (GWs). Based on the near real-time Fermi/GBM data, we developed a low-latency pipeline to automatically calculate the temporal and spectral properties of GRBs. With this pipeline, some important parameters can be obtained, such as T90 and fluence, within ~20 minutes after the GRB trigger. For ~90% GRBs, T90 and fluence are consistent with the GBM catalog results within 2 sigma errors. This pipeline has been used by the Gamma-ray Bursts Polarimeter (POLAR) and the Insight Hard X-ray Modulation Telescope (Insight-HXMT) to follow up the bursts of interest. For GRB 170817A, the first EM counterpart of GW events detected by Fermi/GBM and INTEGRAL/SPI-ACS, the pipeline gave T90 and spectral information in 21 minutes after the GBM trigger, providing important information for POLAR and Insight-HXMT observations.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

We present a concept of surface decomposition extended from double Fourier series to nonnegative sinusoidal wave surfaces, on the basis of which linear ion sources apply to the ultra-precision fabrication of complex surfaces and diffractive optics. It is the first time that we have a surface descriptor for building a relationship between the fabrication process of optical surfaces and the surface characterization based on PSD analysis, which akin to Zernike polynomials used for mapping the relationship between surface errors and Seidel aberrations. Also, we demonstrate that the one-dimensional scanning of linear ion source is applicable to the removal of surface errors caused by small-tool polishing in raster scan mode as well as the fabrication of beam sampling grating of high diffractive uniformity without a post-processing procedure. The simulation results show that, in theory, optical fabrication with linear ion source is feasible and even of higher output efficiency compared with the conventional approach.

We propose the use of the kinetic Sunyaev-Zel'dovich (kSZ) effect to probe the circumgalactic medium (CGM), with the aid of a spectroscopic survey covering the same area of a SZ survey. One can design an optimal estimator of the kSZ effect of the CGM with a matched filter, and construct the cross correlation between the estimator and the peculiar velocity recovered from the galaxy survey, which can be measured by stacking a number of galaxies. We investigate two compelling profiles for the CGM, the MB profile (Maller & Bullock 2004) and the $\beta$ profile, and estimate the detectability against the synergy of a fiducial galaxy survey with number density $10^{-3}h^3\,$ Mpc$^{-3}$ and an ACT-like SZ survey. We show that the shape of the filter does not change much with redshift for the $\beta$ profile, while there are significant side lobes at $z<0.1$ for the MB profile. By stacking $\sim 10^4$ Milky Way-size halos around z $\sim 0.5$, one can get $\gtrsim$ 1 $\sigma$ signal to noise (S/N) for the both profiles. The S/N increases with decreasing redshift before it reaches a maximum ($\sim$ 7.5 at z $\simeq$ 0.15 for the MB profile, $\sim 19$ at $z\simeq 0.03$ for the $\beta$ profile). Due to the large beam size, a Planck-like CMB survey can marginally detect the kSZ signal by stacking the same number of galaxies at $z<0.1$. The search for the CGM in realistic surveys will involve dividing the galaxies into subsamples with similar redshift and mass of host halos, and scaling the results presented here to obtain the S/N.

We report on a clear solar-cycle variation of the Sun's shadow in the 10 TeV cosmic-ray flux observed by the Tibet air shower array during a full solar cycle from 1996 to 2009. In order to clarify the physical implications of the observed solar cycle variation, we develop numerical simulations of the Sun's shadow, using the Potential Field Source Surface (PFSS) model and the Current Sheet Source Surface (CSSS) model for the coronal magnetic field. We find that the intensity deficit in the simulated Sun's shadow is very sensitive to the coronal magnetic field structure, and the observed variation of the Sun's shadow is better reproduced by the CSSS model. This is the first successful attempt to evaluate the coronal magnetic field models by using the Sun's shadow observed in the TeV cosmic-ray flux.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

In this work we study the constraints on the dark matter interaction with the standard model particles, from the observations of dark matter relic density, the direct detection experiments of CDMS and XENON, and the indirect detection of the antiproton-to-proton ratio by PAMELA. A model independent way is adopted in the study by constructing the effective interaction operators between dark matter and standard model particles. The most general 4-fermion operators are investigated. We find that the constraints from different observations are complementary with each other. Especially the spin independent scattering gives very strong constraints for corresponding operators. In some cases the indirect detection of antiproton-to-proton data can actually be more sensitive than the direct detection or relic density for light dark matter (less than 70 GeV).

Low-scale supersymmetry breaking in string motivated theories implies the presence of O(100) TeV scale moduli, which generically lead to a significant modification of the history of the universe prior to Big Bang Nucleosynthesis. Such an approach implies a non-thermal origin for dark matter resulting from scalar decay, where the lightest supersymmetric particle can account for the observed dark matter relic density. We study the further effect of the decay on the baryon asymmetry of the universe, and find that this can satisfactorily address the problem of the over-production of the baryon asymmetry by the Affleck-Dine mechanism in the MSSM. Remarkably, there is a natural connection between the baryon and dark matter abundances today, which leads to a solution of the `Cosmic Coincidence Problem'.

We propose that the fermionic superpartner of a weak-scale Goldstone boson can be a natural WIMP candidate. The p-wave annihilation of this `Goldstone fermion' into pairs of Goldstone bosons automatically generates the correct relic abundance, whereas the XENON100 direct detection bounds are evaded due to suppressed couplings to the Standard Model. Further, it is able to avoid indirect detection constraints because the relevant s-wave annihilations are small. The interactions of the Goldstone supermultiplet can induce non-standard Higgs decays and novel collider phenomenology.

The thermal Sunyaev-Zel'dovich (tSZ) effect directly measures the thermal pressure of free electrons integrated along the line of sight and thus contains valuable information on the thermal history of the universe. However, the redshift information is entangled in the projection along the line of sight. This projection effect severely degrades the power of the tSZ effect to reconstruct the thermal history. We investigate the tSZ tomography technique to recover this otherwise lost redshift information by cross correlating the tSZ effect with galaxies of known redshifts, or alternatively with matter distribution reconstructed from weak lensing tomography. We investigate in detail the 3D distribution of the gas thermal pressure and its relation with the matter distribution, through our adiabatic hydrodynamic simulation and the one with additional gastrophysics including radiative cooling, star formation and supernova feedback. (1) We find a strong correlation between the gas pressure and matter distribution, with a typical cross correlation coefficient r ~ 0.7 at k . 3h/Mpc and z < 2. This tight correlation will enable robust cross correlation measurement between SZ surveys such as Planck, ACT and SPT and lensing surveys such as DES and LSST, at ~20-100{\sigma} level. (2) We propose a tomography technique to convert the measured cross correlation into the contribution from gas in each redshift bin to the tSZ power spectrum. Uncertainties in gastrophysics may affect the reconstruction at ~ 2% level, due to the ~ 1% impact of gastrophysics on r, found in our simulations. However, we find that the same gastrophysics affects the tSZ power spectrum at ~ 40% level, so it is robust to infer the gastrophysics from the reconstructed redshift resolved contribution.

The kinetic Sunyaev Zel'dovich effect (kSZ) effect is a potentially powerful probe to the missing baryons. However, the kSZ signal is overwhelmed by various contaminations and the cosmological application is hampered by loss of redshift information due to the projection effect. We propose a kSZ tomography method to alleviate these problems, with the aid of galaxy spectroscopic redshift surveys. We propose to estimate the large scale peculiar velocity through the 3D galaxy distribution, weigh it by the 3D galaxy density and adopt the product projected along the line of sight with a proper weighting as an estimator of the true kSZ temperature fluctuation $\Theta$. We thus propose to measure the kSZ signal through the $\Hat{\Theta}$-$\Theta$ cross correlation. This approach has a number of advantages (see details in the abstract of the paper). We test the proposed kSZ tomography against non-adiabatic and adiabatic hydrodynamical simulations. We confirm that $\hat{\Theta}$ is indeed tightly correlated with $\Theta$ at $k\la 1h/$Mpc, although nonlinearities in the density and velocity fields and nonlinear redshift distortion do weaken the tightness of the $\hat{\Theta}$-$\Theta$ correlation. We further quantify the reconstruction noise in $\Hat{\Theta}$ from galaxy distribution shot noise. Based on these results, we quantify the applicability of the proposed kSZ tomography for future surveys. We find that, in combination with the BigBOSS-N spectroscopic redshift survey, the PLANCK CMB experiment will be able to detect the kSZ with an overall significance of $\sim 50\sigma$ and further measure its redshift distribution at many redshift bins over $0<z<2$.

We address the cosmological moduli/gravitino problems and the issue of too little thermal but excessive non-thermal dark matter from the decays of moduli. The main examples we study are the G2-MSSM models arising from M theory compactifications, which allow for a precise calculation of moduli decay rates and widths. We find that the late decaying moduli satisfy both BBN constraints and avoid the gravitino problem. The non-thermal production of wino LSPs, which is a prediction of G2-MSSM models, gives a relic density of about the right order of magnitude.