Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

The Kepler mission observed thousands of transiting exoplanet candidates around hundreds of thousands of FGK dwarf stars. He, Ford, & Ragozzine (2019) (arXiv:1907.07773) applied forward modelling to infer the distribution of intrinsic architectures of planetary systems, developed a clustered Poisson point process model for exoplanetary systems (SysSim) to reproduce the marginal distributions of the observed Kepler population, and showed that orbital periods and planet radii are clustered within a given planetary system. Here, we extend the clustered model to explore correlations between planetary systems and their host star properties. We split the sample of Kepler FGK dwarfs into two halves and model the fraction of stars with planets (between 0.5-10 $R_\oplus$ and 3--300 d), $f_{\rm swpa}$, as a linear function of the Gaia DR2 $b_p-r_p-E^*$ colour, where $E^*$ is a reddening correction. We find that the occurrence of these planetary systems rises significantly towards later type (redder or higher $b_p-r_p-E^*$) stars, with a slope of $df_{\rm swpa}/d(b_p-r_p-E^*) = 0.84_{-0.35}^{+0.37}$. The fraction of stars with planets increases from $f_{\rm swpa} = 0.34_{-0.12}^{+0.12}$ for F2V dwarfs to $f_{\rm swpa} = 0.92_{-0.15}^{+0.08}$ for mid K-dwarfs. About half ($f_{\rm swpa} = 0.57_{-0.10}^{+0.15}$) of all solar-type (G2V) dwarfs harbour a planetary system between 3 and 300 d. We considered three related models and find that the rise in $f_{\rm swpa}$ with $b_p-r_p-E^*$ is robust. While this linear $f_{\rm swpa}(b_p-r_p-E^*)$ model is simple, it can closely match the observed multiplicity distributions of both bluer and redder halves in our sample, suggesting that the architectures of planetary systems around stars of different spectral types may be similar aside from a shift in the overall fraction of planet hosting stars.

Low-surface-brightness galaxies(LSBGs) are defined as galaxies that are fainter than dark night sky and are important for studying our universe. Particularly, edge-on galaxies are useful for the study of rotational velocity and dynamical properties of galaxies. Hence here we focus on searching for edge-on LSBGs. In order to find these edge-on dim galaxies, a series of effects caused by inclination, including the surface brightness profile, internal extinction, and scale length, have been corrected. In this work, we present a catalog of 281 edge-on LSBG candidates, which are selected from the cross-match between SDSS DR7 and the 40% ALFALFA catalog. We also present the properties of these edge-on LSBG candidates including absolute magnitude, central surface brightness, B-V color, scale length, and relative thickness. Our result suggests that the correction of inclination effects is very important for obtaining a complete sample of LSBGs.

It has been argued that the mixed Higgs-$R^2$ model acts as the UV extension of the Higgs inflation, pushing up its cut-off scale in the vacuum close up to the Planck scale. In this letter, we study the inflaton oscillation stage after inflation, focusing on the effective mass of the phase direction of the Higgs field, which can cause a violent preheating process. We find that the "spikes" in the effective mass of the phase direction observed in the Higgs inflation still appear in the mixed Higgs-$R^2$ model. While the spikes appear above the cut-off scale in the Higgs-only case, they appear below the cut-off scale when the model is extended with $R^2$ term though reheating cannot be completed in the violent particle production regime since the spikes get milder.

Observations of exoplanetary systems provide clues about the intrinsic distribution of planetary systems, their architectures, and how they formed. We develop a forward modelling framework for generating populations of planetary systems and "observed" catalogues by simulating the Kepler detection pipeline (SysSim). We compare our simulated catalogues to the Kepler DR25 catalogue of planet candidates, updated to include revised stellar radii from Gaia DR2. We constrain our model based on the observed 1D marginal distributions of orbital periods, period ratios, transit depths, transit depth ratios, transit durations, transit duration ratios, and transit multiplicities. Models assuming planets with independent periods and sizes do not adequately account for the properties of the multiplanet systems. Instead, a clustered point process model for exoplanet periods and sizes provides a significantly better description of the Kepler population, particularly the observed multiplicity and period ratio distributions. We find that $0.56^{+0.18}_{-0.15}$ of FGK stars have at least one planet larger than $0.5 R_\oplus$ between 3 and 300 d. Most of these planetary systems ($\sim 98\%$) consist of one or two clusters with a median of three planets per cluster. We find that the Kepler dichotomy is evidence for a population of highly inclined planetary systems and is unlikely to be solely due to a population of intrinsically single planet systems. We provide a large ensemble of simulated physical and observed catalogues of planetary systems from our models, as well as publicly available code for generating similar catalogues given user-defined parameters.

Within the framework that primordial black holes are formed by the direct gravitational collapse of large primordial density perturbations in the radiation dominated stage, we derive the threshold of the density contrast for the formation of rotating primordial black holes based on the simple Jeans criterion. It is found that the threshold value increases in proportion to the square of the angular momentum. We then apply the recently refined analysis on the formation threshold for non-rotating black holes to the case of rotating black holes, and contrast the derived threshold with the former. Caveats and effects ignored in our analysis are also presented, which suggests that the uncertainties of our result can be addressed only by means of numerical relativity.

Ultra-diffuse galaxies (UDGs) are objects which have very extended morphology and faint central surface brightness. Most UDGs are discovered in galaxy clusters and groups, but also some are found in low density environments. The diffuse morphology and faint surface brightness make them difficult to distinguish from the sky background. Several previous works have suggested that at least some UDGs are consistent with exponential surface brightness profiles (S\'{e}rsic n ~ 1). The surface brightness of exponential disks is enhanced in edge-on systems, so searching for edge-on systems may be an efficient way to select UDGs. In this paper, we focus on searching for edge-on HI-bearing ultra-diffuse sources (HUDS) from the 40% ALFALFA catalog, based on SDSS g- and r-band images. After correcting the observed central surface brightness to a face-on perspective, we discover 11 edge-on HUDS candidates. All these newly discovered HUDS candidates are blue and HI-bearing, similar to other HUDS in 70% ALFALFA catalog, and different from UDGs in clusters.

We present a narrow $\rm H\alpha$-band imaging survey of 357 low surface brightness galaxies (LSBGs) that are selected from the spring sky region of the 40% Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA) HI Survey. All the $\rm H\alpha$ images are obtained from the 2.16 m telescope, operated by Xinglong Observatory of the National Astronomical Observatories, Chinese Academy of Sciences. We provide the $\rm H\alpha$ fluxes and derive the global star formation rates (SFRs) of LSBGs after the Galactic extinction, internal extinction, and [NII] contamination correction. Comparing to normal star-forming galaxies, LSBGs have a similar distribution in the HI surface density ($\rm \Sigma_{HI}$), but their SFRs and star formation surface density ($\rm \Sigma_{SFR}$) are much lower. Our results show that the gas-rich LSBGs selected from the ALFALFA survey obviously deviate from the Kennicutt-Schmidt law, in the relation between the star formation surface density ($\rm \Sigma_{SFR}$) and the gas surface density ($\rm \Sigma_{gas}$). However, they follow the extended Schmidt law well when taking the stellar mass of the galaxy into consideration.

To study the disc central surface brightness ($\mu_0$) distribution in optical and near-infrared bands, we select 708 disc-dominated galaxies within a fixed distance of 57 Mpc from SDSS DR7 and UKIDSS DR10. Then we fit $\mu_0$ distribution by using single and double Gaussian profiles with an optimal bin size for the final sample of 538 galaxies in optical $griz$ bands and near-infrared $YJHK$ bands. Among the 8 bands, we find that $\mu_{0}$ distribution in optical bands can not be much better fitted with double Gaussian profiles. However, for all the near-infrared bands, the evidence of being better fitted by using double Gaussian profiles is positive. Especially for $K$ band, the evidence of a double Gaussian profile being better than a single Gaussian profile for $\mu_{0}$ distribution is very strong, the reliability of which can be approved by 1000 times test for our sample. No dust extinction correction is applied. The difference of $\mu_{0}$ distribution between optical and near-infrared bands could be caused by the effect of dust extinction in optical bands. Due to the sample selection criteria, our sample is not absolutely complete. However, the sample incompleteness does not change the double Gaussian distribution of $\mu_{0}$ in $K$ band. Furthermore, we discuss some possible reasons for the fitting results of $\mu_{0}$ distribution in $K$ band. Conclusively, the double Gaussian distribution of $\mu_{0}$ in $K$ band for our sample may depend on bulge-to-disk ratio, color and disk scalelength, rather than the inclination of sample galaxies, bin size and statistical fluctuations.

We analyze a two-field inflationary model consisting of the Ricci scalar squared ($R^2$) term and the standard Higgs field non-minimally coupled to gravity in addition to the Einstein $R$ term. Detailed analysis of the power spectrum of this model with mass hierarchy is presented, and we find that one can describe this model as an effective single-field model in the slow-roll regime with a modified sound speed. The scalar spectral index predicted by this model coincides with those given by the $R^2$ inflation and the Higgs inflation implying that there is a close relation between this model and the $R^2$ inflation already in the original (Jordan) frame. For a typical value of the self-coupling of the standard Higgs field at the high energy scale of inflation, the role of the Higgs field in parameter space involved is to modify the scalaron mass, so that the original mass parameter in the $R^2$ inflation can deviate from its standard value when non-minimal coupling between the Ricci scalar and the Higgs field is large enough.

We present the observed H$\alpha$ flux and derived star formation rates (SFRs) for a fall sample of low$-$surface$-$brightness galaxies (LSBGs). The sample is selected from the fall sky region of the 40$\%$ ALFALFA {\ion{H}{1}} survey $-$ SDSS DR7 photometric data, and all the $H\alpha$ images were obtained using the 2.16 m telescope, operated by the National Astronomy Observatories, Chinese Academy of Sciences. A total of 111 LSBGs were observed and $H\alpha$ flux was measured in 92 of them. Though almost all the LSBGs in our sample are {\ion{H}{1}}$-$rich, their SFRs derived from the extinction and filter$-$transmission$-$corrected $H\alpha$ flux, are less than 1$M_{\sun}$$yr^{-1}$. LSBGs and star forming galaxies have similar {\ion{H}{1}} surface densities, but LSBGs have much lower SFRs and SFR surface densities than star$-$forming galaxies. Our results show that LSBGs deviate from the Kennicutt-Schmidt law significantly, which indicate that they have low star formation efficiency. The SFRs of LSBGs are close to average SFRs in Hubble time and support the previous arguments that most of the LSBGs are stable systems and they tend to seldom contain strong interactions or major mergers during their star formation histories.

A significant fraction of main sequence stars are part of a triple system. We study the long-term stability and dynamical outcomes of triple stellar systems using a large number of long-term direct N-body integrations with relativistic precession. We find that the previously proposed stability criteria by Eggleton & Kiseleva 1995 and Mardling & Aarseth 2001 predict the stability against ejections reasonably well for a wide range of parameters. Assuming that the triple stellar systems follow orbital and mass distributions from FGK binary stars in the field, we find that in ~1% and ~0.5% of the triple systems lead to a direct head-on collision (impact velocity ~ escape velocity) between main sequence (MS) stars and between a MS star and a stellar-mass compact object, respectively. We conclude that triple interactions are the dominant channel for direct collisions involving a MS star in the field with a rate of one event every ~100 years in the Milky Way. We estimate that the fraction of triple systems that forms short-period binaries is up to ~23% with only up to ~13% being the result of three-body interactions with tidal dissipation, which is consistent with previous work using a secular code.

We present a spectroscopic redshift catalog from the LAMOST Complete Spectroscopic Survey of Pointing Area (LaCoSSPAr) in the Southern Galactic Cap (SGC), which is designed to observe all sources (Galactic and extra-galactic) by using repeating observations with a limiting magnitude of $r=18.1~mag$ in two $20~deg^2$ fields. The project is mainly focusing on the completeness of LAMOST ExtraGAlactic Surveys (LEGAS) in the SGC, the deficiencies of source selection methods and the basic performance parameters of LAMOST telescope. In both fields, more than 95% of galaxies have been observed. A post-processing has been applied to LAMOST 1D spectrum to remove the majority of remaining sky background residuals. More than 10,000 spectra have been visually inspected to measure the redshift by using combinations of different emission/absorption features with uncertainty of $\sigma_{z}/(1+z)<0.001$. In total, there are 1528 redshifts (623 absorption and 905 emission line galaxies) in Field A and 1570 redshifts (569 absorption and 1001 emission line galaxies) in Field B have been measured. The results show that it is possible to derive redshift from low SNR galaxies with our post-processing and visual inspection. Our analysis also indicates that up to 1/4 of the input targets for a typical extra-galactic spectroscopic survey might be unreliable. The multi-wavelength data analysis shows that the majority of mid-infrared-detected absorption (91.3%) and emission line galaxies (93.3%) can be well separated by an empirical criterion of $W2-W3=2.4$. Meanwhile, a fainter sequence paralleled to the main population of galaxies has been witnessed both in $M_r$/$W2-W3$ and $M_*$/$W2-W3$ diagrams, which could be the population of luminous dwarf galaxies but contaminated by the edge-on/highly inclined galaxies ($\sim30\%$).

We have obtained CO(J=2-1) spectra of nine face-on low surface brightness galaxies(LSBGs) using the JCMT 15-meter telescope and observed Ha images using the 2.16-meter telescope of NAOC. As no CO has been detected, only upper limits on the H2 masses are given. The upper limits of total molecular hydrogen masses are about (1.2-82.4) x 10^7 Msun. Their star formation rates are mainly lower than 0.4 Msun/yr and star formation efficiencies are lower than 1.364 x 10^10 /yr. Our results show that the absence of molecular gas content is the direct reason for the low star formation rate. The low star formation efficiency probably resulted from the low efficiency of HI gas transforming to H2 gas.

A characteristic signature of String Gas Cosmology is primordial power spectra for scalar and tensor modes which are almost scale-invariant but with a red tilt for scalar modes but a blue tilt for tensor modes. This feature, however, can also be realized in the so-called G-inflation model, in which Horndeski operators are introduced which leads to a blue tensor tilt by softly breaking the Null Energy Condition. In this article we search for potential observational differences between these two cosmologies by performing detailed perturbation analyses based on the Effective Field Theory approach. Our results show that, although both two models produce blue tilted tensor perturbations, they behave differently in three aspects. Firstly, String Gas Cosmology predicts a specific consistency relation between the index of the scalar modes ns and that of tensor ones nt, which is hard to be reproduced by G-inflation. Secondly, String Gas Cosmology typically predicts non-Gaussianities which are highly suppressed on observable scales, while G-inflation gives rise to observationally large non-Gaussianities because the kinetic terms in the action become important during inflation. However, after finely tuning the model parameters of G-inflation it is possible to obtain a blue tensor spectrum and negligible non-Gaussianities with a degeneracy between the two models. This degeneracy can be broken by a third observable, namely the scale dependence of the nonlinearity parameter, which vanishes for G-inflation but has a blue tilt in the case of String Gas Cosmology. Therefore, we conclude that String Gas Cosmology is in principle observationally distinguishable from the single field inflationary cosmology, even allowing for modifications such as G-inflation.

Planet formation theories predict a large but still undetected population of short-period terrestrial planets orbiting brown dwarfs. Should specimens of this population be discovered transiting relatively bright and nearby brown dwarfs, the Jupiter-size and the low luminosity of their hosts would make them exquisite targets for detailed atmospheric characterisation with JWST and future ground-based facilities. The eventual discovery and detailed study of a significant sample of transiting terrestrial planets orbiting nearby brown dwarfs could prove to be useful not only for comparative exoplanetology but also for astrobiology, by bringing us key information on the physical requirements and timescale for the emergence of life. In this context, we present a search for transit-signals in archival time-series photometry acquired by the Spitzer Space Telescope for a sample of 44 nearby brown dwarfs. While these 44 targets were not particularly selected for their brightness, the high precision of their Spitzer light curves allows us to reach sensitivities below Earth-sized planets for 75% of the sample and down to Europa-sized planets on the brighter targets. We could not identify any unambiguous planetary signal. Instead, we could compute the first limits on the presence of planets on close-in orbits. We find that within a 1.28 day orbit, the occurrence rate of planets with a radius between 0.75 and 3.25 R$_\oplus$ is {\eta} < 67 $\pm$ 1%. For planets with radii between 0.75 and 1.25 R$_\oplus$, we place a 95% confident upper limit of {\eta} < 87 $\pm$ 3%. If we assume an occurrence rate of {\eta} = 27% for these planets with radii between 0.75 and 1.25 R$_\oplus$, as the discoveries of the Kepler-42b and TRAPPIST-1b systems would suggest, we estimate that 175 brown dwarfs need to be monitored in order to guarantee (95%) at least one detection.

Long-term stability of deformable mirrors (DM) is a critical performance requirement for instruments requiring open-loop corrections. The effects of temperature changes in the DM performance are equally critical for such instruments. This paper investigates the long-term stability of three different Iris AO PTT111 DMs that were calibrated at different times ranging from 13 months to nearly 29 months prior to subsequent testing. Performance testing showed that only a small increase in positioning errors occurred from the initial calibration date to the test dates. The increases in errors ranged from as little as 1.38 nm rms after 18 months to 5.68 nm rms after 29 months. The paper also studies the effects of small temperature changes, up to 6.2{\deg}C around room temperature. For three different arrays, the errors ranged from 0.62-1.42 nm rms/{\deg}C. Removing the effects of packaging shows that errors are $\le$0.50 nm rms/{\deg}C. Finally, measured data showed that individual segments deformed $\le$0.11 nm rms/{\deg}C when heated.

We report the results of a weakly-interacting massive particle (WIMP) dark matter search using the full 80.1\;live-day exposure of the first stage of the PandaX experiment (PandaX-I) located in the China Jin-Ping Underground Laboratory. The PandaX-I detector has been optimized for detecting low-mass WIMPs, achieving a photon detection efficiency of 9.6\%. With a fiducial liquid xenon target mass of 54.0\,kg, no significant excess event were found above the expected background. A profile likelihood analysis confirms our earlier finding that the PandaX-I data disfavor all positive low-mass WIMP signals reported in the literature under standard assumptions. A stringent bound on the low mass WIMP is set at WIMP mass below 10\,GeV/c$^2$, demonstrating that liquid xenon detectors can be competitive for low-mass WIMP searches.

We apply one of lazy learning methods named k-nearest neighbor algorithm (kNN) to estimate the photometric redshifts of quasars, based on various datasets from the Sloan Digital Sky Survey (SDSS), UKIRT Infrared Deep Sky Survey (UKIDSS) and Wide-field Infrared Survey Explorer (WISE) (the SDSS sample, the SDSS-UKIDSS sample, the SDSS-WISE sample and the SDSS-UKIDSS-WISE sample). The influence of the k value and different input patterns on the performance of kNN is discussed. kNN arrives at the best performance when k is different with a special input pattern for a special dataset. The best result belongs to the SDSS-UKIDSS-WISE sample. The experimental results show that generally the more information from more bands, the better performance of photometric redshift estimation with kNN. The results also demonstrate that kNN using multiband data can effectively solve the catastrophic failure of photometric redshift estimation, which is met by many machine learning methods. By comparing the performance of various methods for photometric redshift estimation of quasars, kNN based on KD-Tree shows its superiority with the best accuracy for our case.

We show that, even for extremely stiff systems, explicit integration may compete in both accuracy and speed with implicit methods if algebraic methods are used to stabilize the numerical integration. The required stabilizing algebra depends on whether the system is well-removed from equilibrium or near equilibrium. This paper introduces a quantitative distinction between these two regimes and addresses the former case in depth, presenting explicit asymptotic methods appropriate when the system is extremely stiff but only weakly equilibrated. A second paper examines quasi-steady-state methods as an alternative to asymptotic methods in systems well away from equilibrium and a third paper extends these methods to equilibrium conditions in extremely stiff systems using partial equilibrium methods. All three papers present systematic evidence for timesteps competitive with implicit methods. Because an explicit method can execute a timestep faster than an implicit method, algebraically-stabilized explicit algorithms might permit integration of larger networks than have been feasible before in various disciplines.

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.

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.

Using the data taken from Tibet II High Density (HD) Array (1997 February-1999 September) and Tibet-III array (1999 November-2005 November), our previous northern sky survey for TeV $\gamma-$ray point sources has now been updated by a factor of 2.8 improved statistics. From $0.0^{\circ}$ to $60.0^{\circ}$ in declination (Dec) range, no new TeV $\gamma-$ray point sources with sufficiently high significance were identified while the well-known Crab Nebula and Mrk421 remain to be the brightest TeV $\gamma-$ray sources within the field of view of the Tibet air shower array. Based on the currently available data and at the 90% confidence level (C.L.), the flux upper limits for different power law index assumption are re-derived, which are approximately improved by 1.7 times as compared with our previous reported limits.

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 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.

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.