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.

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.

Lijiang 2.4-meter Telescope(LJT), the largest common-purpose optical telescope in China, has been applied to the world-wide astronomers since 2008. It is located at Gaomeigu site, Lijiang Observatory(LJO), the southwest of China. The site has very good observational conditions. Since 10-year operation, several instruments have been equipped on the LJT. Astronomers can perform both photometric and spectral observations. The main scientific goals of LJT include photometric and spectral evolution of supernova, reverberation mapping of active galactic nucleus, physical properties of binary star and near-earth object(comet and asteroid), identification of exoplanet, and all kinds of transients. Until now, the masses of 41 high accretion rate black holes have been measured, and more than 168 supernova have been identified by the LJT. More than 190 papers related to the LJT have been published. In this paper, the general observation condition of the Gaomeigu site is introduced at first. Then, the LJT structure is described in detail, including the optical, mechanical, motion and control system. The specification of all the instruments, and some detailed parameters of the YFOSC is also presented. Finally, some important scientific results and future expectations are summarized.

Planetary systems can survive the stellar evolution, as evidenced by the atmospheric metal pollution and dusty disks of single white dwarfs. Recent observations show that 1 to 4 percent of single white dwarfs are accompanied by dusty disks, while the occurrence rate of metal pollution is about 25 to 50 percent . The dusty disks and metal pollution have been associated with accretion of remnant planetary systems around white dwarfs, yet the relation between these two phenomena is still unclear. Here we suggest an evolutionary scenario to link the two observational phenomena. By analyzing a sample of metal polluted white dwarfs, we find that the mass accretion rate onto the white dwarf generally follows a broken power law decay, which matches well with the theoretical prediction, if assuming dust accretion is primarily driven by Poynting-Robertson drag and the dust source is primarily delivered via dynamically falling asteroids perturbed by a Jovian planet. The presence of disks is mainly at the early stage (about 0.1 to 0.7 billion years) of the whole process of metal pollution, which is detectable until about 8 billion years, naturally explaining the fraction (about 2-16 percent) of metal-polluted white dwarfs having dusty disks. The success of this scenario also implies that the configuration of an asteroid belt with an outer gas giant might be common around stars of several solar masses.

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.

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.

Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferometer operation has been installed and the first test run was accomplished in March and April of 2016 with a rather simple configuration. The initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we summarize the construction of KAGRA, including the study of the advantages and challenges of building an underground detector and the operation of the iKAGRA interferometer together with the geophysics interferometer that has been constructed in the same tunnel.

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.

In cosmological N-body simulations, the baryon effects on the cold dark matter (CDM) halos can be used to solve the small scale problems in $\Lambda$CDM cosmology, such as cusp-core problem and missing satellites problem. It turns out that the resultant total density profiles (baryons plus CDM), for halos with mass ranges from dwarf galaxies to galaxy clusters, can match the observations of the rotation curves better than NFW profile. In our previous work, however, we found that such density profiles fail to match the most recent strong gravitational lensing observations. In this paper, we do the converse: we fit the most recent strong lensing observations with the predicted lensing probabilities based on the so-called $(\alpha,\beta,\gamma)$ double power-law profile, and use the best-fit parameters ($\alpha=3.04, \beta=1.39, \gamma=1.88$) to calculate the rotation curves. We find that, at outer parts for a typical galaxy, the rotation curve calculated with our fitted density profile is much lower than observations and those based on simulations, including the NFW profile. This again verifies and strengthen the conclusions in our previous works: in $\Lambda$CDM paradigm, it is difficult to reconcile the contradictions between the observations for rotation curves and strong gravitational lensing.

Simulations are expected to be the powerful tool to investigate the baryon effects on dark matter (DM) halos. Recent high resolution, cosmological hydrodynamic simulations (\citealt{Cintio14}, DC14) predict that the inner density profiles of DM halos depend systematically on the ratio of stellar to DM mass ($M_{\ast}/M_{\rm halo}$) which is thought to be able to provide good fits to the observed rotation curves of galaxies. The DC14 profile is fitted from the simulations which are confined to $M_{\rm halo}\le 10^{12}M_{\sun}$, in order to investigate the physical processes that may affect all halos, we extrapolate it to much larger halo mass, including that of galaxy clusters. The inner slope of DC14 profile is flat for low halo mass, it approaches 1 when the halo mass increases towards $10^{12}M_{\sun}$ and decreases rapidly after that mass. We use DC14 profile for lenses and find that it predicts too few lenses compared with the most recent strong lensing observations SQLS (\citealt{Inada12}). We also calculate the strong lensing probabilities for a simulated density profile which continues the halo mass from the mass end of DC14 ($\sim 10^{12}M_{\sun}$) to the mass that covers the galaxy clusters (\citealt{Schaller15}, Schaller15), and find that this Schaller15 model predict too many lenses compared with other models and SQLS observations. Interestingly, Schaller15 profile has no core, however, like DC14, the rotation curves of the simulated halos are in excellent agreement with observational data. Furthermore, we show that the standard two-population model SIS+NFW cannot match the most recent SQLS observations for large image separations.

A new air-shower core-detector array (YAC : Yangbajing Air-shower Core-detector array) has been developed to measure the primary cosmic-ray composition at the "knee" energies in Tibet, China, focusing mainly on the light components. The prototype experiment (YAC-I) consisting of 16 detectors has been constructed and operated at Yangbajing (4300 m a.s.l.) in Tibet since May 2009. YAC-I is installed in the Tibet-III AS array and operates together. In this paper, we performed a Monte Carlo simulation to check the sensitivity of YAC-I+Tibet-III array to the cosmic-ray light component of cosmic rays around the knee energies, taking account of the observation conditions of actual YAC-I+Tibet-III array. The selection of light component from others was made by use of an artificial neural network (ANN). The simulation shows that the light-component spectrum estimated by our methods can well reproduce the input ones within 10\% error, and there will be about 30\% systematic errors mostly induced by the primary and interaction models used. It is found that the full-scale YAC and the Tibet-III array is powerful to study the cosmic-ray composition, in particular, to obtain the energy spectra of protons and helium nuclei around the knee energies.

Recently, AMS-02 reported their observed results of cosmic rays(CRs). In addition to the AMS-02 data, we add HESS data to estimate the spectra of CR electrons and the diffuse gamma rays above TeV. In the conventional diffusion model, a global analysis is performed on the spectral features of CR electrons and the diffuse gamma rays by GALRPOP package. The results show that the spectrum structure of the primary component of CR electrons can not be fully reproduced by a simple power law and the relevant break is around hundred GeV. At 99\% C.L., the injection indices above the break decrease from 2.54 to 2.35, but the ones below the break are only in the range 2.746 - 2.751. The spectrum of CR electrons does not need to add TeV cutoff to match the features of HESS data too. Based on the difference between the fluxes of CR electrons and the primary component of them, the predicted excess of CR positrons is consistent with the interpretations as pulsar or dark matter. In the analysis of the Galactic diffuse gamma rays with the indirect constraint of AMS-02 and HESS data, it is found that the fluxes of Galactic diffuse gamma rays are consistent with GeV data of Fermi-LAT in the high latitude regions. The results indicate that the inverse Compton scattering(IC) is the dominant component in the range of the hundred GeV to tens of TeV respectively from the high latitude regions to the low ones, and in the all regions of Galaxy the flux of diffuse gamma rays is less than CR electrons at the energy scale of 20 TeV.

Aiming at the observation of cosmic-ray chemical composition at the "knee" energy region, we have been developinga new type air-shower core detector (YAC, Yangbajing Air shower Core detector array) to be set up at Yangbajing (90.522$^\circ$ E, 30.102$^\circ$ N, 4300 m above sea level, atmospheric depth: 606 g/m$^2$) in Tibet, China. YAC works together with the Tibet air-shower array (Tibet-III) and an underground water cherenkov muon detector array (MD) as a hybrid experiment. Each YAC detector unit consists of lead plates of 3.5 cm thick and a scintillation counter which detects the burst size induced by high energy particles in the air-shower cores. The burst size can be measured from 1 MIP (Minimum Ionization Particle) to $10^{6}$ MIPs. The first phase of this experiment, named "YAC-I", consists of 16 YAC detectors each having the size 40 cm $\times$ 50 cm and distributing in a grid with an effective area of 10 m$^{2}$. YAC-I is used to check hadronic interaction models. The second phase of the experiment, called "YAC-II", consists of 124 YAC detectors with coverage about 500 m$^2$. The inner 100 detectors of 80 cm $\times $ 50 cm each are deployed in a 10 $\times$ 10 matrix from with a 1.9 m separation and the outer 24 detectors of 100 cm $\times$ 50 cm each are distributed around them to reject non-core events whose shower cores are far from the YAC-II array. YAC-II is used to study the primary cosmic-ray composition, in particular, to obtain the energy spectra of proton, helium and iron nuclei between 5$\times$$10^{13}$ eV and $10^{16}$ eV covering the "knee" and also being connected with direct observations at energies around 100 TeV. We present the design and performance of YAC-II in this paper.

Currently, the Japanese gravitational wave laser interferometer KAGRA is under construction in the Kamioka mine. As one main feature, it will employ sapphire mirrors operated at a temperature of 20K to reduce the impact from thermal noise. To reduce seismic noise, the mirrors will also be suspended from multi-stage pendulums. Thus the heat load deposited in the mirrors by absorption of the circulating laser light as well as heat load from thermal radiation will need to be extracted through the last suspension stage. This stage will consist of four thin sapphire fibers with larger heads necessary to connect the fibers to both the mirror and the upper stage. In this paper, we discuss heat conductivity measurements on different fiber candidates. While all fibers had a diameter of 1.6mm, different surface treatments and approaches to attach the heads were analyzed. Our measurements show that fibers fulfilling the basic KAGRA heat conductivity requirement of $\kappa\geq $5000W/m/K at 20K are technologically feasible.

PSR J0218$+$4232 is a millisecond pulsar (MSP) with a flux density $\sim$ 0.9 mJy at 1.4 GHz. It is very bright in the high-energy X-ray and $\gamma$-ray domains. We conducted an astrometric program using the European VLBI Network (EVN) at 1.6 GHz to measure its proper motion and parallax. A model-independent distance would also help constrain its $\gamma$-ray luminosity. We achieved a detection of signal-to-noise ratio S/N > 37 for the weak pulsar in all five epochs. Using an extragalactic radio source lying 20 arcmin away from the pulsar, we estimate the pulsar's proper motion to be $\mu_{\alpha}\cos\delta=5.35\pm0.05$ mas yr$^{-1}$ and $\mu_{\delta}=-3.74\pm 0.12$ mas yr$^{-1}$, and a parallax of $\pi=0.16\pm0.09$ mas. The very long baseline interferometry (VLBI) proper motion has significantly improved upon the estimates from long-term pulsar timing observations. The VLBI parallax provides the first model-independent distance constraints: $d=6.3^{+8.0}_{-2.3}$ kpc, with a corresponding $3\sigma$ lower-limit of $d=2.3$ kpc. This is the first pulsar trigonometric parallax measurement based solely on EVN observations. Using the derived distance, we believe that PSR J0218$+$4232 is the most energetic $\gamma$-ray MSP known to date. The luminosity based on even our 3$\sigma$ lower-limit distance is high enough to pose challenges to the conventional outer gap and slot gap models.

In interferometric cryogenic gravitational wave detectors, there are plans to cool mirrors and their suspension systems (payloads) in order to reduce thermal noise, that is, one of the fundamental noise sources. Because of the large payload masses (several hundred kg in total) and their thermal isolation, a cooling time of several months is required. Our calculation shows that a high-emissivity coating (e.g. a diamond-like carbon (DLC) coating) can reduce the cooling time effectively by enhancing radiation heat transfer. Here, we have experimentally verified the effect of the DLC coating on the reduction of the cooling time.

We demonstrate how observations of pulsars can be used to help navigate a spacecraft travelling in the solar system. We make use of archival observations of millisecond pulsars from the Parkes radio telescope in order to demonstrate the effectiveness of the method and highlight issues, such as pulsar spin irregularities, which need to be accounted for. We show that observations of four millisecond pulsars every seven days using a realistic X-ray telescope on the spacecraft throughout a journey from Earth to Mars can lead to position determinations better than approx. 20km and velocity measurements with a precision of approx. 0.1m/s.

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 EAS hybrid experiment has been designed by constructing a YAC (Yangbajing Air shower Core) detector array inside the existing Tibet-III air shower array. The first step of YAC, called "YAC-I", consists of 16 plastic scintillator units (4 rows times 4 columns) each with an area of 40 cm * 50 cm which is used to check hadronic interaction models used in AS simulations. A Monte Carlo study shows that YAC-I can record high energy electromagnetic component in the core region of air showers induced by primary particles of several tens TeV energies where the primary composition is directly measured by space experiments. It may provide a direct check of the hadronic interaction models currently used in the air shower simulations in the corresponding energy region. In present paper, the method of the observation and the sensitivity of the characteristics of the observed events to the different interaction models are discussed.

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.

Pulsed gamma-ray emission from millisecond pulsars (MSPs) has been detected by the sensitive Fermi, which sheds light on studies of the emission region and mechanism. In particular, the specific patterns of radio and gamma-ray emission from PSR J0101-6422 challenge the popular pulsar models, e.g. outer gap and two-pole caustic (TPC) models. Using the three dimension (3D) annular gap model, we have jointly simulated radio and gamma-ray light curves for three representative MSPs (PSR J0034-0534, PSR J0101-6422 and PSR J0437-4715) with distinct radio phase lags and present the best simulated results for these MSPs, particularly for PSR J0101-6422 with complex radio and gamma-ray pulse profiles and for PSR J0437-4715 with a radio interpulse. It is found that both the gamma-ray and radio emission originate from the annular gap region located in only one magnetic pole, and the radio emission region is not primarily lower than the gamma-ray region in most cases. In addition, the annular gap model with a small magnetic inclination angle instead of "orthogonal rotator" can account for MSPs' radio interpulse with a large phase separation from the main pulse. The annular gap model is a self-consistent model not only for young pulsars but also MSPs, and multi-wavelength light curves can be fundamentally explained by this model.

Using observations of pulsars from the Parkes Pulsar Timing Array (PPTA) project we develop the first pulsar-based timescale that has a precision comparable to the uncertainties in international atomic timescales. Our ensemble of pulsars provides an Ensemble Pulsar Scale (EPS) analogous to the free atomic timescale Echelle Atomique Libre (EAL). The EPS can be used to detect fluctuations in atomic timescales and therefore can lead to a new realisation of Terrestrial Time, TT(PPTA11). We successfully follow features known to affect the frequency of the International Atomic Timescale (TAI) and we find marginally significant differences between TT(PPTA11) and TT(BIPM11). We discuss the various phenomena that lead to a correlated signal in the pulsar timing residuals and therefore limit the stability of the pulsar timescale.

A correlation between the peak luminosity and the peak energy has been found by Yonetoku et al. as $L_{p}\propto E_{p,i}^{2.0}$ for 11 pre-Swift long gamma-ray bursts. In this study, for a greatly expanded sample of 148 long gamma-ray bursts in the Swift era, we find that the correlation still exists, but most likely with a slightly different power-law index, i.e., $L_{p}\propto E_{p,i} ^{1.7}$. In addition, we have collected 17 short gamma-ray bursts with necessary data. It is found that the correlation of $L_{p}\propto E_{p,i} ^{1.7}$ also exists for this sample of short events. It is argued that the radiation mechanism of both long and short gamma-ray bursts should be similar, i.e., of quasi-thermal origin caused by the photosphere and the dissipation occurring very near the central engine. Some key parameters of the process are constrained. Our results suggest that the radiation process of both long and short bursts may be dominated by thermal emission, rather than the single synchrotron radiation. This might put strong physical constraints on the theoretical models.

An algorithm is proposed for constructing a group (ensemble) pulsar time based on the application of optimal Wiener filters. This algorithm makes it possible to separate the contributions of variations of the atomic time scale and of the pulsar rotation to barycentric residual deviations of the pulse arrival times. The method is applied to observations of the pulsars PSR B1855+09 and PSR B1937+21, and is used to obtain corrections to UTC relative to the group pulsar time PT$_{\rm ens}$. Direct comparison of the terrestrial time TT(BIPM06) and the group pulsar time PT$_{\rm ens}$ shows that they disagree by no more than $0.4\pm 0.17\; \mu$s. Based on the fractional instability of the time difference TT(BIPM06) -- PT$_{\rm ens}$, a new limit for the energy density of the gravitational-wave background is established at the level $\Omega_g {h}^2\sim 10^{-9}$.

We show how pulsar observations may be used to construct a time standard that is independent of terrestrial time standards. The pulsar time scale provides a method to determine the stability of terrestrial time standards over years to decades. Here, we summarise the method, provide initial results and discuss the possibilities and limitations of our pulsar time scale.

Standard cosmology has many successes on large scales, but faces some fundamental difficulties on small, galactic scales. One such difficulty is the cusp/core problem. High resolution observations of the rotation curves for dark matter dominated low surface brightness (LSB) galaxies imply that galactic dark matter halos have a density profile with a flat central core, whereas N-body structure formation simulations predict a divergent (cuspy) density profile at the center. It has been proposed that this problem can be resolved by stellar feedback driving turbulent gas motion that erases the initial cusp. However, strong gravitational lensing prefers a cuspy density profile for galactic halos. In this paper, we use the most recent high resolution observations of the rotation curves of LSB galaxies to fit the core size as a function of halo mass, and compare the resultant lensing probability to the observational results for the well defined combined sample of the Cosmic Lens All-Sky Survey (CLASS) and Jodrell Bank/Very Large Array Astrometric Survey (JVAS). The lensing probabilities based on such density profiles are too low to match the observed lensing in CLASS/JVAS. High baryon densities in the galaxies that dominate the lensing statistics can reconcile this discrepancy, but only if they steepen the mass profile rather than making it more shallow. This places contradictory demands upon the effects of baryons on the central mass profiles of galaxies.

The evolution of the galaxy stellar mass function is especially useful to test the current model of galaxy formation. Observational data have revealed a few inconsistencies with predictions from the $\Lambda {\rm CDM}$ model. For example, most massive galaxies have already been observed at very high redshifts, and they have experienced only mild evolution since then. In conflict with this, semi-analytical models of galaxy formation predict an insufficient number of massive galaxies at high redshift and a rapid evolution between redshift 1 and 0 . In addition, there is a strong correlation between star formation rate and stellar mass for star-forming galaxies, which can be roughly reproduced with the model, but with a normalization that is too low at high redshift. Furthermore, the stellar mass density obtained from the integral of the cosmic star formation history is higher than the measured one by a factor of 2. In this paper, we study these issues using a semi-analytical model that includes: 1) cold gas accretion in massive halos at high redshift; 2) tidal stripping of stellar mass from satellite galaxies; and 3) an evolving stellar initial mass function (bottom-light) with a higher gas recycle fraction. Our results show that the combined effects from 1) and 2) can predict sufficiently massive galaxies at high redshifts and reproduce their mild evolution at low redshift, While the combined effects of 1) and 3) can reproduce the correlation between star formation rate and stellar mass for star-forming galaxies across wide range of redshifts. A bottom-light/top-heavy stellar IMF could partly resolve the conflict between the stellar mass density and cosmic star formation history.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

Cosmological numerical simulations of galaxy formation have led to the cuspy density profile of a pure cold dark matter halo toward the center, which is in sharp contradiction with the observations of the rotation curves of cold dark matter-dominated dwarf and low surface brightness disk galaxies, with the latter tending to favor mass profiles with a flat central core. Many efforts have been devoted to resolve this cusp-core problem in recent years, among them, baryon-cold dark matter interactions are considered to be the main physical mechanisms erasing the cold dark matter (CDM) cusp into a flat core in the centers of all CDM halos. Clearly, baryon-cold dark matter interactions are not customized only for CDM-dominated disk galaxies, but for all types, including giant ellipticals. We first fit the most recent high resolution observations of rotation curves with the Burkert profile, then use the constrained core size-halo mass relation to calculate the lensing frequency, and compare the predicted results with strong lensing observations. Unfortunately, it turns out that the core size constrained from rotation curves of disk galaxies cannot be extrapolated to giant ellipticals. We conclude that, in the standard cosmological paradigm, baryon-cold dark matter interactions are not universal mechanisms for galaxy formation, and therefore, they cannot be true solutions to the cusp-core problem.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

Aims. We try to identify ultra-cool dwarfs from the seventh Data Release of the Sloan Digital Sky Survey (SDSS DR7) with SDSS i-z and r-z colors. We also obtain proper motion data from SDSS, 2MASS, and UKIDSS and improve spectral typing from SDSS and 2MASS photometric colors. Methods. We selected ultra-cool dwarf candidates from the SDSS DR7 with new photometric selection criteria, which are based on a parameterization study of known L and T dwarfs. The objects are then cross-identified with the Two Micron All Sky Survey and the Fourth Data Release of the UKIRT Infrared Deep Sky Survey (UKIDSS DR4). We derive proper motion constraints by combining SDSS, 2MASS, and UKIDSS positional information. In this way we are able to assess, to some extent, the credence of our sample using a multi epoch approach, which complements spectroscopic confirmation. Some of the proper motions are affected by short baselines, but, as a general tool, this method offers great potential to confirm faint L dwarfs as UKIDSS coverage increases. In addition we derive updated color-spectral type relations for L and T dwarfs with SDSS and 2MASS magnitudes. Results. We present 59 new nearby M and L dwarfs selected from the imaging catalog of the SDSS DR7, including proper motions and spectral types calculated from the updated color-spectral type relations. and obtain proper motions from SDSS, 2MASS, and UKIDSS for all of our objects.

$f(R)$ gravity, capable of driving the late-time acceleration of the universe, is emerging as a promising alternative to dark energy. Various $f(R)$ gravity models have been intensively tested against probes of the expansion history, including type Ia supernovae (SNIa), the cosmic microwave background (CMB) and baryon acoustic oscillations (BAO). In this paper we propose to use the statistical lens sample from Sloan Digital Sky Survey Quasar Lens Search Data Release 3 (SQLS DR3) to constrain $f(R)$ gravity models. This sample can probe the expansion history up to $z\sim2.2$, higher than what probed by current SNIa and BAO data. We adopt a typical parameterization of the form $f(R)=R-\alpha H^2_0(-\frac{R}{H^2_0})^\beta$ with $\alpha$ and $\beta$ constants. For $\beta=0$ ($\Lambda$CDM), we obtain the best-fit value of the parameter $\alpha=-4.193$, for which the 95% confidence interval that is [-4.633, -3.754]. This best-fit value of $\alpha$ corresponds to the matter density parameter $\Omega_{m0}=0.301$, consistent with constraints from other probes. Allowing $\beta$ to be free, the best-fit parameters are $(\alpha, \beta)=(-3.777, 0.06195)$. Consequently, we give $\Omega_{m0}=0.285$ and the deceleration parameter $q_0=-0.544$. At the 95% confidence level, $\alpha$ and $\beta$ are constrained to [-4.67, -2.89] and [-0.078, 0.202] respectively. Clearly, given the currently limited sample size, we can only constrain $\beta$ within the accuracy of $\Delta\beta\sim 0.1$ and thus can not distinguish between $\Lambda$CDM and $f(R)$ gravity with high significance, and actually, the former lies in the 68% confidence contour. We expect that the extension of the SQLS DR3 lens sample to the SDSS DR5 and SDSS-II will make constraints on the model more stringent.

We recalculate the strong lensing probability as a function of the image separation in TeVeS (tensor-vector-scalar) cosmology, which is a relativistic version of MOND (MOdified Newtonian Dynamics). The lens is modeled by the Hernquist profile. We assume an open cosmology with $\Omega_b=0.04$ and $\Omega_\Lambda=0.5$ and three different kinds of interpolating functions. Two different galaxy stellar mass functions (GSMF) are adopted: PHJ (Panter-Heavens-Jimenez, 2004) determined from SDSS data release one and Fontana (Fontana et al., 2006) from GOODS-MUSIC catalog. We compare our results with both the predicted probabilities for lenses by Singular Isothermal Sphere (SIS) galaxy halos in LCDM (lambda cold dark matter) with Schechter-fit velocity function, and the observational results of the well defined combined sample of Cosmic Lens All-Sky Survey (CLASS) and Jodrell Bank/Very Large Array Astrometric Survey (JVAS). It turns out that the interpolating function $\mu(x)=x/(1+x)$ combined with Fontana GSMF matches the results from CLASS/JVAS quite well.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

Mayall II = G1 is one of the most luminous globular clusters (GCs) known in M31. New deep, high-resolution observations with the Advanced Camera for Surveys on the {\sl Hubble Space Telescope} are used to provide accurate photometric data to the smallest radii yet. In particular, we present the precise variation of ellipticity and position angle, and of surface brightness for the core of the object. Based on these accurate photometric data, we redetermine the structural parameters of G1 by fitting a single-mass isotropic King model. We derive a core radius, $r_c=0.21\pm0.01\arcsec (=0.78\pm0.04 \rm{pc})$, a tidal radius, $r_t=21.8\pm1.1\arcsec (=80.7\pm3.9 \rm{pc})$, and a concentration index $c=\log (r_t/r_c)=2.01\pm0.02$. The central surface brightness is 13.510 mag arcsec$^{-2}$. We also calculate the half-light radius, at $r_h=1.73\pm0.07\arcsec(=6.5\pm0.3 \rm{pc})$. The results show that, within 10 core radii, a King model fits the surface brightness distribution well. We find that this object falls in the same region of the $M_V$ vs. $\log R_h$ diagram as $\omega$ Centauri, M54 and NGC 2419 in the Milky Way. All three of these objects have been claimed to be the stripped cores of now defunct dwarf galaxies. We discuss in detail whether GCs, stripped cores of dwarf spheroidals and normal dwarf galaxies form a continuous distribution in the $M_V$ versus $\log R_h$ plane, or if GCs and dwarf spheroidals constitute distinct classes of objects; we present arguments in favour of this latter view.

We calculate the strong lensing probability as a function of the image-separation $\Delta\theta$ in TeVeS (tensor-vector-scalar) cosmology, which is a relativistic version of MOND (MOdified Newtonian Dynamics). The lens, often an elliptical galaxy, is modeled by the Hernquist profile. We assume a flat cosmology with $\Omega_b=1-\Omega_\Lambda=0.04$ and the simplest interpolating function $\mu(x)={\rm min}(1,x)$. For comparison, we recalculated the probabilities for lenses by Singular Isothermal Sphere (SIS) galaxy halos in LCDM with Schechter-fit velocity function. The amplification bias is calculated based on the magnification of the second bright image rather than the total of the two brighter images. Our calculations show that the Hernquist model predicts insufficient but acceptable probabilities in flat TeVeS cosmology compared with the results of the well defined combined sample of Cosmic Lens All-Sky Survey (CLASS) and Jodrell Bank/Very Large Array Astrometric Survey (JVAS); at the same time, it predicts higher probabilities than SIS model in LCDM at small image separations.

We use the statistics of strong gravitational lenses to investigate whether mass profiles with a flat density core are supported. The probability for lensing by halos modeled by a nonsingular truncated isothermal sphere (NTIS) with image separations greater than a certain value (ranging from zero to ten arcseconds) is calculated. NTIS is an analytical model for the postcollapse equilibrium structure of virialized objects derived by Shapiro, Iliev & Raga. This profile has a soft core and matches quite well with the mass profiles of dark matter-dominated dwarf galaxies deduced from their observed rotation curves. It also agrees well with the NFW (Navarro-Frenk-White) profile at all radii outside of a few NTIS core radii. Unfortunately, comparing the results with those for singular lensing halos (NFW and SIS+NFW) and strong lensing observations, the probabilities for lensing by NTIS halos are far too low. As this result is valid for any other nonsingular density profiles (with a large core radius), we conclude that nonsingular density profiles (with a large core radius) for CDM halos are ruled out by statistics of strong gravitational lenses.

In the standard flat cosmological constant ($\Lambda$) cold dark matter (CDM) cosmology, a model of two populations of lens halos for strong gravitational lensing can reproduce the results of the Jodrell-Bank VLA Astrometric Survey (JVAS) and the Cosmic Lens All-Sky Survey (CLASS) radio survey. In such a model, lensing probabilities are sensitive to three parameters: the concentration parameter $c_1$, the cooling mass scale $M_\mathrm{c}$ and the value of the CDM power spectrum normalization parameter $\sigma_8$. The value ranges of these parameters are constrained by various observations. However, we found that predicted lensing probabilities are also quite sensitive to the flux density (brightness) ratio $q_{\mathrm{r}}$ of the multiple lensing images, which has been, in fact, a very important selection criterion of a sample in any lensing survey experiments. We re-examine the above mentioned model by considering the flux ratio and galactic central Super Massive Black Holes (SMBHs), in flat, low-density cosmological models with different cosmic equations of state $\omega$, and find that the predicted lensing probabilities without considering $q_{\mathrm{r}}$ are over-estimated. A low value of $q_\mathrm{r}$ can be compensated by raising the cooling mass scale $M_\mathrm{c}$ in fitting the predicted lensing probabilities to JVAS/CLASS observations. In order to determine the cosmic equation of state $\omega$, the uncertainty in $M_\mathrm{c}$ must be resolved. The effects of SMBHs cannot be detected by strong gravitational lensing method when $q_{\mathrm{r}}\leq 10$.

We investigate the linear amplitude of mass fluctuations in the universe, $\sigma_8$, and the present mass density parameter of the Universe, $\Omega_\mathrm{m}$, from the statistical strong gravitational lensing. We use the two populations of lens halos model with fixed cooling mass scale $M_\mathrm{c}=3\times 10^{13}h^{-1}M_{\sun}$ to match the observed lensing probabilities, and leave $\sigma_8$ or $\Omega_\mathrm{m}$ as a free parameter to be constrained by data. Another varying parameter is the equation of state of dark energy $\omega$, and its typical values of -1, -2/3, -1/2 and -1/3 are investigated. We find that $\sigma_8$ is degenerate with $\Omega_\mathrm{m}$ in a way similar to that suggested by present day cluster abundance as well as cosmic shear lensing measurements: $\sigma_8\Omega_\mathrm{m}^{0.6}\approx 0.33$ (Bahcall & Bode\cite{bahcall03a} and references therein). However, both $\sigma_8\leq 0.7$ and $\Omega_\mathrm{m}\leq 0.2$ can be safely ruled out, the best value is when $\sigma_8=1.0$, $\Omega_\mathrm{m}=0.3$ and $\omega=-1$. This result is different from that obtained by Bahcall & Bode (\cite{bahcall03a}), who gives $\sigma_8 =0.98\pm 0.1$ and $\Omega_m =0.17\pm 0.05$. For $\sigma_8=1.0$, higher value of $\Omega_\mathrm{m}=0.35$ requires $\omega=-2/3$ and $\Omega_\mathrm{m}=0.40$ requires $\omega=-1/2$.

We study the angular momentum profiles both for dark matter and for gas within virialized halos, using a statistical sample of halos drawn from cosmological hydrodynamics simulations. Three simulations have been analyzed, one is the ``non-radiative'' simulation, and the other two have radiative cooling. We find that the gas component on average has a larger spin and contains a smaller fraction of mass with negative angular momentum than its dark matter counterpart in the non-radiative model. As to the cooling models, the gas component shares approximately the same spin parameter as its dark matter counterpart, but the hot gas has a higher spin and is more aligned in angular momentum than dark matter, while the opposite holds for the cold gas. After the mass of negative angular momentum is excluded, the angular momentum profile of the hot gas component approximately follows the universal function originally proposed by Bullock et al. for dark matter, though the shape parameter $\mu$ is much larger for hot gas and is comfortably in the range required by observations of disk galaxies. Since disk formation is related to the distribution of hot gas that will cool, our study may explain the fact that the disk component of observed galaxies contains a smaller fraction of low angular momentum material than dark matter in halos.

Lensing probabilities of quasars with image separations greater than \Delta\theta and flux density ratios less than q_{\mathrm{r}} are calculated by foreground dark matter halos in a flat, cosmological constant dominated (\LambdaCDM) universe . The mass density of the lenses is taken to be the Navarro-Frenk-White (NFW) profile on all mass scales, plus a central point mass for low-mass halos with M<M_c=5\times 10^{13}h^{-1}M_{\odot}. We introduce a quantity M_{\mathrm{eff}}, which is a point mass ranging from 1 to 1000 times the mass M_{\bullet} of a supermassive black hole (SMBH) inhabiting the center of each galaxy, to describe the contributions of galactic central SMBHs and galactic bulges to lensing probabilities. The lensing cross section and thus the lensing probability are quite sensitive to the flux density ratio q_{\mathrm{r}} of multiple images in our calculations. It is shown that, to reproduce the lensing survey results of JVAS/CLASS for q_{\mathrm{r}}<10, about 20% of the bulge mass is needed as a point mass for each galaxy. Since there is still considerable uncertainty regarding the value of the spectrum normalization parameter \sigma_8, we investigate the effect of varying this parameter within its entire observational range (from 0.7 to 1.1), and find that low \sigma_8 values (\leq 0.7) are ruled out, and the best fit value is \sigma_8\simeq 1.0.

We investigate how central black holes (BHs) inhabited in galactic dark halos could affect strong gravitational lensing. The distribution of integral lensing probability with image separations are calculated for quasars of redshift 1.5 by foreground dark matter halos. The mass density of dark halos is taken to be the Navarro-Frenk-White (NFW) profile such that, when the mass of a halo is less than $10^{14} M_{\sun}$, its central black holes or a bulge is included as a point mass. The relationship between the masses $M_{\bullet}$ of supermassive black holes and the total gravitational mass $M_{\mathrm{DM}}$ of their host galaxy is adopted from the most recent literature. Only a flat $\Lambda$CDM model is considered here. It is shown that, while a single black hole for each galaxy contributes considerable but not sufficient lensing probabilities at small image separations compared with those without black holes, the bulges (which are about 100--1000 times larger in mass than a typical black hole) would definitely contribute enough probabilities at small image separations, although it gives too high probabilities at large separation angles compared with lensing observations.

We study the angular momentum profile of dark matter halos for a statistical sample drawn from a set of high-resolution cosmological simulations of $256^3$ particles. Two typical Cold Dark Matter (CDM) models have been analyzed, and the halos are selected to have at least $3\times 10^4$ particles in order to reliably measure the angular momentum profile. In contrast with the recent claims of Bullock et al., we find that the degree of misalignment of angular momentum within a halo is very high. About 50 percent of halos have more than 10 percent of halo mass in the mass of negative angular momentum $j$. After the mass of negative $j$ is excluded, the cumulative mass function $M(<j)$ follows approximately the universal function proposed by Bullock et al., though we still find a significant fraction of halos ($\sim 50%$) which exhibit systematic deviations from the universal function. Our results, however, are broadly in good agreement with a recent work of van den Bosch et al.. We also study the angular momentum profile of halos in a Warm Dark Matter (WDM) model and a Self-Interacting Dark Matter (SIDM) model. We find that the angular momentum profile of halos in the WDM is statistically indistinguishable from that in the CDM model, but the angular momentum of halos in the SIDM is reduced by the self-interaction of dark matter.

We explore the possibility of attributing the recent discovery of the hard X-ray variable source CXOM82 J095550.2+694047 in M82 to the gravitational magnification by an intervening stellar obeject along the line of sight as microlens. The duration of the event ($>84$ days) allows us to set robust constraints on the mass and location of the mircolensing object when combined with the dynamical properties of the Galactic halo, the M82 and typical globular clusters. Except for the extremely low probablity, the microlensing magnification by MACHO in either the Galactic halo or the M82 halo is able to explain the X-ray variability of CXOM82 J095550.2+694047. It is hoped that the lensing hypothesis can be soon tested by the light curve measurement.

We present a quantitative analysis of the time-delay effect on the cosmic background radiation (CBR) by static gravitational potential of galaxy clusters. This is primarily motivated by growing observational evidence that clusters have essentially experienced no-evolution since redshift $z\approx1$, indicating that the contribution of a time-dependent potential to CBR anisotropy discussed in literature could be rather small for the dynamically-relaxed clusters. Using the softened isothermal sphere model and the universal density profile for the mass distribution of rich clusters, we calculate the CBR anisotropy by the time-delay effect and compare it with those generated by the thermal and kinematic S-Z effects as well as by the transverse motion of clusters. While it is unlikely that the time-delay effect is detectable in the current S-Z measurement because of its small amplitude of $10^{-6}$-$10^{-7}$ and its achromaticity, it nevertheless leads to an uncertainty of $\sim10%$ in the measurement of the kinematic S-Z effect of clusters. Future cosmological application of the peculiar velocity of clusters to be measured through the S-Z effect should therefore take this uncertainty into account.