This paper describes the liquid scintillator veto system for the COSINE-100 dark matter experiment and its performance. The COSINE-100 detector consists of 8 NaI(Tl) crystals immersed in 2,200~L of linear alkylbenzene based liquid scintillator. The liquid scintillator tags between 65 and 75\% of the internal $^{40}$K background in the 2--6~keV energy region. We also describe the background model for the liquid scintillator, which is primarily used to assess its energy calibration and threshold.

The surface $^{210}$Pb is one of the main background sources for dark-matter-search experiments using NaI(Tl) crystals, its spectral features associated with the beta-decay events for energies less than 60 keV depends on the depth distribution of $^{210}$Pb in the surface of an NaI(Tl) crystal. Therefore, we must understand the profile of surface $^{210}$Pb to precisely model the background measurement in the low-energy region for the low-background experiment using NaI(Tl) crystals. We estimate the depth profile of the surface $^{210}$Pb contamination by modeling the measured spectrum of the alpha emission from the decay of $^{210}$Po at the decay sequence of the surface $^{210}$Pb contamination that is obtained using an $^{222}$Rn-contaminated crystal. Moreover, we describe the energy spectra of the surface contaminations by performing Geant4 Monte Carlo simulations that are used for fitting an exponential curve, which is a function of the surface depth, with the measured data, by using the log likelihood method. The low- and high-energy events from the beta decay of surface $^{210}$Pb are also modeled to improve the depth profile for shallow depths. We simulate the energy spectra from beta decays of $^{210}$Pb that are exponentially distributed in the surface by following two exponential functions that consider mean-depth parameters in the data fitting; we observed that the energy spectra are in good agreement with the measured data.

There is a long standing debate about whether or not the annual modulation signal reported by the DAMA/LIBRA collaboration is induced by Weakly Interacting Massive Particles~(WIMP) in the galaxy's dark matter halo scattering from nuclides in their NaI(Tl) crystal target/detector. This is because regions of WIMP-mass vs. WIMP-nucleon cross-section parameter space that can accommodate the DAMA/LIBRA-phase1 modulation signal in the context of the standard WIMP dark matter galactic halo and isospin-conserving~(canonical), spin-independent~(SI) WIMP-nucleon interactions have been excluded by many of other dark matter search experiments including COSINE-100, which uses the same NaI(Tl) target/detector material. Moreover, the recently released DAMA/LIBRA-phase2 results are inconsistent with an interpretation as WIMP-nuclide scattering via the canonical SI interaction and prefer, instead, isospin-violating or spin-dependent interactions. Dark matter interpretations of the DAMA/LIBRA signal are sensitive to the NaI(Tl) scintillation efficiency for nuclear recoils, which is characterized by so-called quenching factors~(QF), and the QF values used in previous studies differ significantly from recently reported measurements, which may have led to incorrect interpretations of the DAMA/LIBRA signal. In this article, the compatibility of the DAMA/LIBRA and COSINE-100 results, in light of the new QF measurements is examined for different possible types of WIMP-nucleon interactions. The resulting allowed parameter space regions associated with the DAMA/LIBRA signal are explicitly compared with 90\% confidence level upper limits from the initial 59.5~day COSINE-100 exposure. With the newly measured QF values, the allowed 3$\sigma$ regions from the DAMA/LIBRA data are still generally excluded by the COSINE-100 data.

COSINE-100 is a direct detection dark matter search experiment that uses a 106 kg array of eight NaI(Tl) crystals that are kept underground at the Yangyang Underground Laboratory to avoid cosmogenic activation of radioisotopes by cosmic rays. Even though the cosmogenic activity is declining with time, there are still significant background rates from the remnant nuclides. In this paper, we report measurements of cosmogenic isotope contaminations with less than one year half-lives that are based on extrapolations of the time dependent activities of their characteristic energy peaks to activity rates at the time the crystals were deployed underground. For longer-lived $^{109}$Cd ($T_{1/2}=1.6$ y) and $^{22}$Na ($T_{1/2}=2.6$ y), we investigate time correlations of characteristic $\gamma$/X-ray peaks. The inferred sea-level production rates are compared with caluclations based on the ACTIVIA and MENDL-2 model calculations and experimental data. For $^{3}$H, which has a long, 12.3 year half-life, we evaluated the activity levels from the exposure times and determined a cosmogenic activation rate that is consistent with other measurements.

We present new constraints on the dark matter-induced annual modulation signal using 1.7 years, of COSINE-100 data with a total exposure of 97.7 kg$\cdot$years. The COSINE-100 experiment, consisting of 106 kg of NaI(Tl) target material, is designed to carry out a model-independent test of DAMA/LIBRA's claim of WIMP discovery by searching for the same annual modulation signal using the same NaI(Tl) target. The crystal data show a 2.7 cpd/kg/keV background rate on average in the 2--6 keV energy region of interest. Using a $\chi$-squared minimization method we observe best fit values for modulation amplitude and phase of 0.0092$\pm$0.0067 cpd/kg/keV and 127.2$\pm$45 d, respectively.

We present results from a search for solar axions with the COSINE-100 experiment. We find no evidence of solar axion events from a data-set of 6,303.9 kg$\cdot$days exposure and set a 90\,\% confidence level upper limit on the axion-electron coupling, $g_{ae}$, at 1.70~$\times$~$10^{-11}$ for an axion mass less than 1\,keV/c$^2$. This limit excludes QCD axions heavier than 0.59\,eV/c$^2$ in the DFSZ model and 168.1\,eV/c$^2$ in the KSVZ model.

A search for inelastic boosted dark matter (iBDM) using the COSINE-100 detector with 59.5 days of data is presented. This relativistic dark matter is theorized to interact with the target material through inelastic scattering with electrons, creating a heavier state that subsequently produces standard model particles, such as an electron-positron pair. In this study, we search for this electron-positron pair in coincidence with the initially scattered electron as a signature for an iBDM interaction. No excess over the predicted background event rate is observed. Therefore, we present limits on iBDM interactions under various hypotheses, one of which allows us to explore an area of the experimental search for iBDM using a terrestrial detector.

Using two volume-limited Main galaxy samples of the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we explore influences of galaxy interactions on AGN activity. It is found that in the faint volume-limited sample, paired galaxies have a slightly higher AGN fraction than isolated galaxies, whereas in the luminous volume-limited sample, an opposite trend can be observed. The significance is <1 sigma. Thus, we do not observe strong evidence that interactions or mergers likely trigger the AGN activity.

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.