We present the detection of VHE gamma-ray emission above 100 TeV from HAWC J2227+610 with the HAWC observatory. Combining our observations with previously published results by VERITAS, we interpret the gamma-ray emission from HAWC J2227+610 as emission from protons with a lower limit in their cutoff energy of 800 TeV. The most likely source of the protons is the associated supernova remnant G106.3+2.7, making it a good candidate for a Galactic PeVatron. However, a purely leptonic origin of the observed emission cannot be excluded at this time.

In the light of the latest measurements on the total $e^+ + e^-$ flux by CALET and DAMPE experiments, we revisit the multicomponent leptonically decaying dark matter (DM) explanations to the cosmic-ray electron/positron excesses observed previously. Especially, we use the single and double-component DM models to explore the compatibility of the AMS-02 positron fraction with the new CALET or DAMPE data. It turns out that neither single nor double-component DM models are able to fit the AMS-02 positron fraction and DAMPE total $e^+ + e^-$ flux data simultaneously. On the other hand, for the combined AMS-02 and CALET dataset, both the single and double-component DM models can provide reasonable fits. If we further take into the diffuse $\gamma$-ray constraints from Fermi-LAT, only the double-component DM models are allowed.

Due to the high energies and long distances to the sources, astrophysical observations provide a unique opportunity to test possible signatures of Lorentz invariance violation (LIV). Superluminal LIV enables the decay of photons at high energy. The High Altitude Water Cherenkov (HAWC) Observatory is among the most sensitive gamma-ray instruments currently operating above 10 TeV. HAWC finds evidence of 100 TeV photon emission from at least four astrophysical sources. These observations exclude, for the strongest of the limits set, the LIV energy scale to $2.2\times10^{31}$ eV, over 1800 times the Planck energy and an improvement of 1 to 2 orders of magnitude over previous limits.

The Large Underground Xenon (LUX) dark matter search was a 250-kg active mass dual-phase time projection chamber that operated by detecting light and ionization signals from particles incident on a xenon target. In December 2015, LUX reported a minimum 90% upper C.L. of 6e-46 cm^2 on the spin-independent WIMP-nucleon elastic scattering cross section based on a 1.4e4 kg*day exposure in its first science run. Tension between experiments and the absence of a definitive positive detection suggest it would be prudent to search for WIMPs outside the standard spin-independent/spin-dependent paradigm. Recent theoretical work has identified a complete basis of 14 independent effective field theory (EFT) operators to describe WIMP-nucleon interactions. In addition to spin-independent and spin-dependent nuclear responses, these operators can produce novel responses such as angular-momentum-dependent and spin-orbit couplings. Here we report on a search for all 14 of these EFT couplings with data from LUX's first science run. Limits are placed on each coupling as a function of WIMP mass.

In face of the small-scale structure problems of the collisionless cold dark matter (DM) paradigm, a popular remedy is to introduce a strong DM self-interaction which can be generated nonperturbatively by a MeV-scale light mediator. However, if such the mediator is unstable and decays into SM particles, the model is severely constrained by the DM direct and indirect detection experiments. In the present paper, we study a model of a self-interacting fermionic DM, endowed with a light stable scalar mediator. In this model, the DM relic abundance is dominated by the fermionic DM particle which is generated mainly via the freeze-out of its annihilations to the stable mediator. Since this channel is invisible, the DM indirect detection constraints should be greatly relaxed. Furthermore, the direct detection signals are suppressed to an unobservable level since fermionic DM scatterings with a nucleon appear at one-loop level. By further studying the bounds from the CMB and BBN on the visible channels involving the dark sector, we show that there is a large parameter space which can generate appropriate DM self-interactions at dwarf galaxy scales, while remaining compatible with other experimental constraints.

Cosmic rays, along with stellar radiation and magnetic fields, are known to make up a significant fraction of the energy density of galaxies such as the Milky Way. When cosmic rays interact in the interstellar medium, they produce gamma-ray emission which provides an important indication of how the cosmic rays propagate. Gamma rays from the Andromeda Galaxy (M31), located 785 kpc away, provide a unique opportunity to study cosmic-ray acceleration and diffusion in a galaxy with a structure and evolution very similar to the Milky Way. Using 33 months of data from the High Altitude Water Cherenkov Observatory, we search for TeV gamma rays from the galactic plane of M31. We also investigate past and present evidence of galactic activity in M31 by searching for Fermi Bubble-like structures above and below the galactic nucleus. No significant gamma-ray emission is observed, so we use the null result to compute upper limits on the energy density of cosmic rays $>10$ TeV in M31. The computed upper limits are approximately ten times higher than expected from the extrapolation of the Fermi LAT results.

Primordial Black Holes (PBHs) may have been created by density fluctuations in the early Universe and could be as massive as $> 10^9$ solar masses or as small as the Planck mass. It has been postulated that a black hole has a temperature inversely-proportional to its mass and will thermally emit all species of fundamental particles via Hawking Radiation. PBHs with initial masses of $\sim 5 \times 10^{14}$ g (approximately one gigaton) should be expiring today with bursts of high-energy gamma radiation in the GeV--TeV energy range. The High Altitude Water Cherenkov (HAWC) Observatory is sensitive to gamma rays with energies of $\sim$300 GeV to past 100 TeV, which corresponds to the high end of the PBH burst spectrum. With its large instantaneous field-of-view of $\sim 2$ sr and a duty cycle over 95%, the HAWC Observatory is well suited to perform an all-sky search for PBH bursts. We conducted a search using 959 days of HAWC data and exclude the local PBH burst rate density above $3400~\mathrm{pc^{-3}~yr^{-1}}$ at 99% confidence, the strongest limit on the local PBH burst rate density from any existing electromagnetic measurement.

We present the first catalog of gamma-ray sources emitting above 56 and 100 TeV with data from the High Altitude Water Cherenkov (HAWC) Observatory, a wide field-of-view observatory capable of detecting gamma rays up to a few hundred TeV. Nine sources are observed above 56 TeV, all of which are likely Galactic in origin. Three sources continue emitting past 100 TeV, making this the highest-energy gamma-ray source catalog to date. We report the integral flux of each of these objects. We also report spectra for three highest-energy sources and discuss the possibility that they are PeVatrons.

We present a novel analysis technique for liquid xenon time projection chambers that allows for a lower threshold by relying on events with a prompt scintillation signal consisting of single detected photons. The energy threshold of the LUX dark matter experiment is primarily determined by the smallest scintillation response detectable, which previously required a 2-fold coincidence signal in its photomultiplier arrays, enforced in data analysis. The technique presented here exploits the double photoelectron emission effect observed in some photomultiplier models at vacuum ultraviolet wavelengths. We demonstrate this analysis using an electron recoil calibration dataset and place new constraints on the spin-independent scattering cross section of weakly interacting massive particles (WIMPs) down to 2.5 GeV/c$^2$ WIMP mass using the 2013 LUX dataset. This new technique is promising to enhance light WIMP and astrophysical neutrino searches in next-generation liquid xenon experiments.

LUX-ZEPLIN (LZ) is a next generation dark matter direct detection experiment that will operate 4850 feet underground at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. Using a two-phase xenon detector with an active mass of 7 tonnes, LZ will search primarily for low-energy interactions with Weakly Interacting Massive Particles (WIMPs), which are hypothesized to make up the dark matter in our galactic halo. In this paper, the projected WIMP sensitivity of LZ is presented based on the latest background estimates and simulations of the detector. For a 1000 live day run using a 5.6 tonne fiducial mass, LZ is projected to exclude at 90% confidence level spin-independent WIMP-nucleon cross sections above $1.6 \times 10^{-48}$ cm$^{2}$ for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP. Additionally, a $5\sigma$ discovery potential is projected reaching cross sections below the existing and projected exclusion limits of similar experiments that are currently operating. For spin-dependent WIMP-neutron(-proton) scattering, a sensitivity of $2.7 \times 10^{-43}$ cm$^{2}$ ($8.1 \times 10^{-42}$ cm$^{2}$) for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP is expected. With underground installation well underway, LZ is on track for commissioning at SURF in 2020.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

The scattering of dark matter (DM) particles with sub-GeV masses off nuclei is difficult to detect using liquid xenon-based DM search instruments because the energy transfer during nuclear recoils is smaller than the typical detector threshold. However, the tree-level DM-nucleus scattering diagram can be accompanied by simultaneous emission of a Bremsstrahlung photon or a so-called "Migdal" electron. These provide an electron recoil component to the experimental signature at higher energies than the corresponding nuclear recoil. The presence of this signature allows liquid xenon detectors to use both the scintillation and the ionization signals in the analysis where the nuclear recoil signal would not be otherwise visible. We report constraints on spin-independent DM-nucleon scattering for DM particles with masses of 0.4-5 GeV/c$^2$ using 1.4$\times10^4$ kg$\cdot$day of search exposure from the 2013 data from the Large Underground Xenon (LUX) experiment for four different classes of mediators. This analysis extends the reach of liquid xenon-based DM search instruments to lower DM masses than has been achieved previously.

List of proceedings from the HAWC Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA.

The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is a ground-based air shower array designed to detect Cherenkov light produced in water by secondary particles from atmospheric air showers. In order to improve the sensitivity at the highest energies, especially for the shower cores falling outside the main array, 345 smaller Water Cherenkov Detectors (WCDs) were installed around the main array, the outrigger array. This extension increased the instrumented area of HAWC by a factor of four. With the increased size of the array, and the ability to detect shower particles further away from the core, understanding of the time structure of the shower front is crucial for accurate direction reconstruction and mandates proper modeling. In this contribution, we present a model of the shower front as expected to be observed by the outriggers obtained from Monte-Carlo simulations. Applying this model to shower reconstruction, the improvements on air shower parameter are studied.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

We report results from an extensive set of measurements of the \b{eta}-decay response in liquid xenon.These measurements are derived from high-statistics calibration data from injected sources of both $^{3}$H and $^{14}$C in the LUX detector. The mean light-to-charge ratio is reported for 13 electric field values ranging from 43 to 491 V/cm, and for energies ranging from 1.5 to 145 keV.

Due to the universality of gravitational interactions, it is generally expected that a stochastic gravitational wave (GW) background could form during the reheating period when the inflaton perturbatively decays with the emission of gravitons. Previously, only models in which the inflaton dominantly decays into a pair of light scalar and/or fermion particles were considered in the literature. In the present paper, we focus on the cases with a vector particle pair in the final decay product. The differential decay rates for the three-body gravitational inflaton decays are presented for two typical couplings between the inflaton and vector fields, from which we predict their respective GW frequency spectra. It turns out that, similar to the scalar and fermion cases, the obtained GW spectra is too high in frequency to be observed by the current and near-future GW detection experiments and calls for a new design of high-frequency GW detectors.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

Various dark matter models predict annual and diurnal modulations of dark matter interaction rates in Earth-based experiments as a result of the Earth's motion in the halo. Observation of such features can provide generic evidence for detection of dark matter interactions. This paper reports a search for both annual and diurnal rate modulations in the LUX dark matter experiment using over 20 calendar months of data acquired between 2013 and 2016. This search focuses on electron recoil events at low energies, where leptophilic dark matter interactions are expected to occur and where the DAMA experiment has observed a strong rate modulation for over two decades. By using the innermost volume of the LUX detector and developing robust cuts and corrections, we obtained a stable event rate of 2.3$\pm$0.2~cpd/keV$_{\text{ee}}$/tonne, which is among the lowest in all dark matter experiments. No statistically significant annual modulation was observed in energy windows up to 26~keV$_{\text{ee}}$. Between 2 and 6~keV$_{\text{ee}}$, this analysis demonstrates the most sensitive annual modulation search up to date, with 9.2$\sigma$ tension with the DAMA/LIBRA result. We also report no observation of diurnal modulations above 0.2~cpd/keV$_{\text{ee}}$/tonne amplitude between 2 and 6~keV$_{\text{ee}}$.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

It is well known that a gravitational wave (GW) experiences the damping effect when it propagates in a fluid with nonzero shear viscosity. In this paper, we propose a new method to constrain the GW damping rate and thus the fluid shear viscosity. By defining the effective distance which incorporates damping effects, we can transform the GW strain expression in a viscous Universe into the same form as that in a perfect fluid. Therefore, the constraints of the luminosity distances from the observed GW events by LIGO and Virgo can be directly applied to the effective distances in our formalism. We exploit the lognormal likelihoods for the available GW effective distances and a Gaussian likelihood for the luminosity distance inferred from the electromagnetic radiation observation of the binary neutron star merger event GW170817. Our fittings show no obvious damping effects in the current GW data, and the upper limit on the damping rate with the combined data is $6.75 \times 10^{-4}\,{\rm Mpc}^{-1}$ at 95\% confidence level. By assuming that the dark matter self-scatterings are efficient enough for the hydrodynamic description to be valid, we find that a GW event from its source at a luminosity distance $D\gtrsim 10^4\;\rm Mpc$ can be used to put a constraint on the dark matter self-interactions.

We study a vector dark matter (VDM) model in which the dark sector couples to the Standard Model sector via a Higgs portal. If the portal coupling is small enough the VDM can be produced via the freeze-in mechanism. It turns out that the electroweak phase transition have a substantial impact on the prediction of the VDM relic density. We further assume that the dark Higgs boson which gives the VDM mass is so light that it can induce strong VDM self-interactions and solve the small-scale structure problems of the Universe. As illustrated by the latest LUX data, the extreme smallness of the Higgs portal coupling required by the freeze-in mechanism implies that the dark matter direct detection bounds are easily satisfied. However, the model is well constrained by the indirect detections of VDM from BBN, CMB, AMS-02, and diffuse $\gamma$/X-rays. Consequently, only when the dark Higgs boson mass is at most of ${\cal O}({\rm keV})$ does there exist a parameter region which leads to a right amount of VDM relic abundance and an appropriate VDM self-scattering while satisfying all other constraints simultaneously.

By the end of 2014, a cosmic ray muon telescope was installed at Zhongshan Station in Antarctic and has been continuously collecting data since then. It is the first surface muon telescope to be built in Antarctic. In June 2015, five CMEs were ejected towards the Earth initiating a big large Forbush decrease (FD) event. We conduct a comprehensive study of the galactic cosmic ray intensity fluctuations during the FD using the data from cosmic ray detectors of multiple stations (Zhongshan, McMurdo, South Polar and Nagoya) and he solar wind measurements from ACE and WIND. A pre-increase before the shock arrival was observed. Distinct differences exist in the timelines of the galactic cosmic ray recorded by the neutron monitors and the muon telescopes. FD onset for Zhongshan muon telescope is delayed (2.5h) with respect to SSC onset. This FD had a profile of four-step decrease. The traditional one- or two-step classification of FDs was inadequate to explain this FD.

We report a new search of weakly interacting massive particles (WIMPs) using the combined low background data sets in 2016 and 2017 from the PandaX-II experiment in China. The latest data set contains a new exposure of 77.1 live day, with the background reduced to a level of 0.8$\times10^{-3}$ evt/kg/day, improved by a factor of 2.5 in comparison to the previous run in 2016. No excess events were found above the expected background. With a total exposure of 5.4$\times10^4$ kg day, the most stringent upper limit on spin-independent WIMP-nucleon cross section was set for a WIMP with mass larger than 100 GeV/c$^2$, with the lowest exclusion at 8.6$\times10^{-47}$ cm$^2$ at 40 GeV/c$^2$.

We report an absolute calibration of the ionization yields($\textit{Q$_y$})$ and fluctuations for electronic recoil events in liquid xenon at discrete energies between 186 eV and 33.2 keV. The average electric field applied across the liquid xenon target is 180 V/cm. The data are obtained using low energy $^{127}$Xe electron capture decay events from the 95.0-day first run from LUX (WS2013) in search of Weakly Interacting Massive Particles (WIMPs). The sequence of gamma-ray and X-ray cascades associated with $^{127}$I de-excitations produces clearly identified 2-vertex events in the LUX detector. We observe the K- (binding energy, 33.2 keV), L- (5.2 keV), M- (1.1 keV), and N- (186 eV) shell cascade events and verify that the relative ratio of observed events for each shell agrees with calculations. The N-shell cascade analysis includes single extracted electron (SE) events and represents the lowest-energy electronic recoil $\textit{in situ}$ measurements that have been explored in liquid xenon.

The first searches for axions and axion-like particles with the Large Underground Xenon (LUX) experiment are presented. Under the assumption of an axio-electric interaction in xenon, the coupling constant between axions and electrons, gAe is tested, using data collected in 2013 with an exposure totalling 95 live-days $\times$ 118 kg. A double-sided, profile likelihood ratio statistic test excludes gAe larger than 3.5 $\times$ 10$^{-12}$ (90% C.L.) for solar axions. Assuming the DFSZ theoretical description, the upper limit in coupling corresponds to an upper limit on axion mass of 0.12 eV/c$^{2}$, while for the KSVZ description masses above 36.6 eV/c$^{2}$ are excluded. For galactic axion-like particles, values of gAe larger than 4.2 $\times$ 10$^{-13}$ are excluded for particle masses in the range 1-16 keV/c$^{2}$. These are the most stringent constraints to date for these interactions.

We present experimental constraints on the spin-dependent WIMP-nucleon elastic cross sections from the total 129.5 kg-year exposure acquired by the Large Underground Xenon experiment (LUX), operating at the Sanford Underground Research Facility in Lead, South Dakota (USA). A profile likelihood ratio analysis allows 90% CL upper limits to be set on the WIMP-neutron (WIMP-proton) cross section of $\sigma_n$ = 1.6$\times 10^{-41}$ cm$^{2}$ ($\sigma_p$ = 5$\times 10^{-40}$ cm$^{2}$) at 35 GeV$c^{-2}$, almost a sixfold improvement over the previous LUX spin-dependent results. The spin-dependent WIMP-neutron limit is the most sensitive constraint to date.

In this Technical Design Report (TDR) we describe the LZ detector to be built at the Sanford Underground Research Facility (SURF). The LZ dark matter experiment is designed to achieve sensitivity to a WIMP-nucleon spin-independent cross section of three times ten to the negative forty-eighth square centimeters.

We report constraints on spin-independent weakly interacting massive particle (WIMP)-nucleon scattering using a 3.35e4 kg-day exposure of the Large Underground Xenon (LUX) experiment. A dual-phase xenon time projection chamber with 250 kg of active mass is operated at the Sanford Underground Research Facility under Lead, South Dakota (USA). With roughly fourfold improvement in sensitivity for high WIMP masses relative to our previous results, this search yields no evidence of WIMP nuclear recoils. At a WIMP mass of 50 GeV/c^2, WIMP-nucleon spin-independent cross sections above 2.2e-46 cm^2 are excluded at the 90% confidence level. When combined with the previously reported LUX exposure, this exclusion strengthens to 1.1e-46 cm^2 at 50 GeV/c^2.

The Large Underground Xenon (LUX) experiment is a dual-phase liquid xenon time projection chamber (TPC) operating at the Sanford Underground Research Facility in Lead, South Dakota. A calibration of nuclear recoils in liquid xenon was performed $\textit{in situ}$ in the LUX detector using a collimated beam of mono-energetic 2.45 MeV neutrons produced by a deuterium-deuterium (D-D) fusion source. The nuclear recoil energy from the first neutron scatter in the TPC was reconstructed using the measured scattering angle defined by double-scatter neutron events within the active xenon volume. We measured the absolute charge ($Q_{y}$) and light ($L_{y}$) yields at an average electric field of 180 V/cm for nuclear recoil energies spanning 0.7 to 74 keV and 1.1 to 74 keV, respectively. This calibration of the nuclear recoil signal yields will permit the further refinement of liquid xenon nuclear recoil signal models and, importantly for dark matter searches, clearly demonstrates measured ionization and scintillation signals in this medium at recoil energies down to $\mathcal{O}$(1 keV).

We propose a new technique for the calibration of nuclear recoils in large noble element dual-phase time projection chambers used to search for WIMP dark matter in the local galactic halo. This technique provides an $\textit{in situ}$ measurement of the low-energy nuclear recoil response of the target media using the measured scattering angle between multiple neutron interactions within the detector volume. The low-energy reach and reduced systematics of this calibration have particular significance for the low-mass WIMP sensitivity of several leading dark matter experiments. Multiple strategies for improving this calibration technique are discussed, including the creation of a new type of quasi-monoenergetic 272 keV neutron source. We report results from a time-of-flight based measurement of the neutron energy spectrum produced by an Adelphi Technology, Inc. DD108 neutron generator, confirming its suitability for the proposed nuclear recoil calibration.

We present constraints on weakly interacting massive particles (WIMP)-nucleus scattering from the 2013 data of the Large Underground Xenon dark matter experiment, including $1.4\times10^{4}\;\mathrm{kg\; day}$ of search exposure. This new analysis incorporates several advances: single-photon calibration at the scintillation wavelength, improved event-reconstruction algorithms, a revised background model including events originating on the detector walls in an enlarged fiducial volume, and new calibrations from decays of an injected tritium $\beta$ source and from kinematically constrained nuclear recoils down to 1.1 keV. Sensitivity, especially to low-mass WIMPs, is enhanced compared to our previous results which modeled the signal only above a 3 keV minimum energy. Under standard dark matter halo assumptions and in the mass range above 4 $\mathrm{GeV}\,c^{-2}$, these new results give the most stringent direct limits on the spin-independent WIMP-nucleon cross section. The 90% C.L. upper limit has a minimum of 0.6 zb at 33 $\mathrm{GeV}\,c^{-2}$ WIMP mass.

We present measurements of the electron-recoil (ER) response of the LUX dark matter detector based upon 170,000 highly pure and spatially-uniform tritium decays. We reconstruct the tritium energy spectrum using the combined energy model and find good agreement with expectations. We report the average charge and light yields of ER events in liquid xenon at 180 V/cm and 105 V/cm and compare the results to the NEST model. We also measure the mean charge recombination fraction and its fluctuations, and we investigate the location and width of the LUX ER band. These results provide input to a re-analysis of the LUX Run3 WIMP search.

Monte Carlo simulations are performed to study the correlation between the ground cosmic ray intensity and near-earth thunderstorms electric field at YBJ (4300 m a.s.l., Tibet, China). The variations of the secondary cosmic ray intensity are found to be highly dependent on the strength and polarity of the electric field. In negative fields and in positive fields greater than 600 V/cm, the total number of ground comic ray positrons and electrons increases with increasing electric field strength. And these values increase more obviously when involving a shower with lower primary energy or a higher zenith angle. While in positive fields ranging from 0 to 600 V/cm, the total number of ground comic ray positrons and electrons declines and the amplitude is up to 3.1% for vertical showers. A decrease of intensity occurs for inclined showers in positive fields less than 500 V/cm, which is accompanied by smaller amplitudes. In this paper, the intensity changes are discussed, especially concerning the decreases in positive electric fields. Our simulation results are in good agreement with ground-based experimental results obtained from ARGO-YBJ and the Carpet air shower array. These results could be helpful in understanding the acceleration mechanisms of secondary charged particles caused by an atmospheric electric field.

ARGO-YBJ, located at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l., Tibet, China), is a full coverage air shower array, with an energy threshold of 300 GeV for gamma-ray astronomy. Most of the recorded events are single front showers, satisfying the trigger requirement of at least 20 particles detected in a given time window. However, in 13% of the events, two randomly arriving showers may be recorded in the same time window, and the second one, in generally smaller, does not need to satisfy the trigger condition. These events are called double front shower events. By using these small showers, well under the trigger threshold, the detector primary energy threshold can be lowered to a few tens of GeV. In this paper, the angular resolution that can be achieved with these events is evaluated by a full Monte Carlo simulation. The ARGO-YBJ sensitivity in detecting gamma-ray bursts (GRBs) by using double front shower events is also studied for various cutoff energies, time durations, and zenith angles of GRBs in the field view of ARGO.

The design and performance of the LUX-ZEPLIN (LZ) detector is described as of March 2015 in this Conceptual Design Report. LZ is a second-generation dark-matter detector with the potential for unprecedented sensitivity to weakly interacting massive particles (WIMPs) of masses from a few GeV/c2 to hundreds of TeV/c2. With total liquid xenon mass of about 10 tonnes, LZ will be the most sensitive experiment for WIMPs in this mass region by the end of the decade. This report describes in detail the design of the LZ technical systems. Expected backgrounds are quantified and the performance of the experiment is presented. The LZ detector will be located at the Sanford Underground Research Facility in South Dakota. The organization of the LZ Project and a summary of the expected cost and current schedule are given.

We revisit the multi-component leptonically decaying dark matter (DM) scenario to explain the possible electron/positron excesses with the recently updated AMS-02 data. We find that both the single- and two-component DM models can fit the positron fraction and $e^+/e^-$ respective fluxes, in which the two-component ones provide better fits. However, for the single-component models, the recent AMS-02 data on the positron fraction limit the DM cutoff to be smaller than 1 TeV, which conflicts with the high-energy behavior of the AMS-02 total $e^++e^-$ flux spectrum, while the two-component DM models do not possess such a problem. We also discuss the constraints from the Fermi-LAT measurement of the diffuse $\gamma$-ray spectrum. We show that the two-component DM models are consistent with the current DM lifetime bounds. In contrast, the best-fit DM lifetimes in the single-component models are actually excluded.

We study a two-component dark matter (DM) model in which the two Majorana fermionic DM components with nearly degenerate masses are stabilized by an $Z_2$ symmetry and interact with the right-handed muon and tau only via real Yukawa couplings, together with an additional $Z_2$-odd singly-charged scalar. In this setup, the decay from the heavy DM to the lighter one via the transition electric dipole yields the 3.55 keV X-ray signal observed recently. The Yukawa couplings in the dark sector are assumed to be hierarchical, so that the observed DM relic abundance can be achieved with the leading s-wave amplitudes without a fine-tuning. We also consider the constraints from flavor physics, DM direct detections and collider searches, respectively.

The identification of impact craters on planetary surfaces provides important information about their geological history. Most studies have relied on individual analysts who map and identify craters and interpret crater statistics. However, little work has been done to determine how the counts vary as a function of technique, terrain, or between researchers. Furthermore, several novel internet-based projects ask volunteers with little to no training to identify craters, and it was unclear how their results compare against the typical professional researcher. To better understand the variation among experts and to compare with volunteers, eight professional researchers have identified impact features in two separate regions of the moon. Small craters (diameters ranging from 10 m to 500 m) were measured on a lunar mare region and larger craters (100s m to a few km in diameter) were measured on both lunar highlands and maria. Volunteer data were collected for the small craters on the mare. Our comparison shows that the level of agreement among experts depends on crater diameter, number of craters per diameter bin, and terrain type, with differences of up to $\sim\pm45%$. We also found artifacts near the minimum crater diameter that was studied. These results indicate that caution must be used in most cases when interpreting small variations in crater size-frequency distributions and for craters $\le10$ pixels across. Because of the natural variability found, projects that emphasize many people identifying craters on the same area and using a consensus result are likely to yield the most consistent and robust information.

The Large Underground Xenon (LUX) dark matter experiment aims to detect rare low-energy interactions from Weakly Interacting Massive Particles (WIMPs). The radiogenic backgrounds in the LUX detector have been measured and compared with Monte Carlo simulation. Measurements of LUX high-energy data have provided direct constraints on all background sources contributing to the background model. The expected background rate from the background model for the 85.3 day WIMP search run is $(2.6\pm0.2_{\textrm{stat}}\pm0.4_{\textrm{sys}})\times10^{-3}$~events~keV$_{ee}^{-1}$~kg$^{-1}$~day$^{-1}$ in a 118~kg fiducial volume. The observed background rate is $(3.6\pm0.4_{\textrm{stat}})\times10^{-3}$~events~keV$_{ee}^{-1}$~kg$^{-1}$~day$^{-1}$, consistent with model projections. The expectation for the radiogenic background in a subsequent one-year run is presented.

The multi-component decaying dark matter (DM) scenario is investigated to explain the possible excesses in the positron fraction by PAMELA and recently confirmed by AMS-02, and in the total $e^+ +e^-$ flux observed by Fermi-LAT. By performing the $\chi^2$ fits, we find that two DM components are already enough to give a reasonable fit of both AMS-02 and Fermi-LAT data. The best-fitted results show that the heavier DM component with its mass 1.5 TeV dominantly decays through the $\mu$-channel, while the lighter one of 100 GeV mainly through the $\tau$-channel. As a byproduct, the fine structure around 100 GeV observed by AMS-02 and Fermi-LAT can be naturally explained by the dropping due to the lighter DM component. With the obtained model parameters by the fitting, we calculate the diffuse $\gamma$-ray emission spectrum in this two-component DM scenario, and find that it is consistent with the data measured by Fermi-LAT. We also construct a microscopic particle DM model to naturally realize the two-component DM scenario, and point out an interesting neutrino signal which is possible to be measured in the near future by IceCube.

LUX, the world's largest dual-phase xenon time-projection chamber, with a fiducial target mass of 118 kg and 10,091 kg-days of exposure thus far, is currently the most sensitive direct dark matter search experiment. The initial null-result limit on the spin-independent WIMP-nucleon scattering cross-section was released in October 2013, with a primary scintillation threshold of 2 phe, roughly 3 keVnr for LUX. The detector has been deployed at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, and is the first experiment to achieve a limit on the WIMP cross-section lower than $10^{-45}$ cm$^{2}$. Here we present a more in-depth discussion of the novel energy scale employed to better understand the nuclear recoil light and charge yields, and of the calibration sources, including the new internal tritium source. We found the LUX data to be in conflict with low-mass WIMP signal interpretations of other results.

The Large Underground Xenon (LUX) experiment, a dual-phase xenon time-projection chamber operating at the Sanford Underground Research Facility (Lead, South Dakota), was cooled and filled in February 2013. We report results of the first WIMP search dataset, taken during the period April to August 2013, presenting the analysis of 85.3 live-days of data with a fiducial volume of 118 kg. A profile-likelihood analysis technique shows our data to be consistent with the background-only hypothesis, allowing 90% confidence limits to be set on spin-independent WIMP-nucleon elastic scattering with a minimum upper limit on the cross section of $7.6 \times 10^{-46}$ cm$^{2}$ at a WIMP mass of 33 GeV/c$^2$. We find that the LUX data are in strong disagreement with low-mass WIMP signal interpretations of the results from several recent direct detection experiments.

The very high energy (VHE; E > 100 GeV) blazar Markarian 501 has a well-studied history of extreme spectral variability and is an excellent laboratory for studying the physical processes within the jets of active galactic nuclei. However, there are few detailed multiwavelength studies of Markarian 501 during its quiescent state, due to its low luminosity. A short-term multiwavelength study of Markarian 501 was coordinated in March 2009, focusing around a multi-day observation with the Suzaku X-ray satellite and including {\gamma}-ray data from VERITAS, MAGIC, and the Fermi Gamma-ray Space Telescope with the goal of providing a well-sampled multiwavelength baseline measurement of Markarian 501 in the quiescent state. The results of these quiescent-state observations are compared to the historically extreme outburst of April 16, 1997, with the goal of examining variability of the spectral energy distribution between the two states. The derived broadband spectral energy distribution shows the characteristic double-peaked profile. We find that the X-ray peak shifts by over two orders of magnitude in photon energy between the two flux states while the VHE peak varies little. The limited shift in the VHE peak can be explained by the transition to the Klein-Nishina regime. Synchrotron self-Compton models are matched to the data and the implied Klein-Nishina effects are explored.

We report on observations of very high-energy gamma rays from the shell-type supernova remnant Cassiopeia A with the VERITAS stereoscopic array of four imaging atmospheric Cherenkov telescopes in Arizona. The total exposure time for these observations is 22 hours, accumulated between September and November of 2007. The gamma-ray source associated with the SNR Cassiopeia A was detected above 200 GeV with a statistical significance of 8.3 s.d. The estimated integral flux for this gamma-ray source is about 3% of the Crab-Nebula flux. The photon spectrum is compatible with a power law dN/dE ~ E^(-Gamma) with an index Gamma = 2.61 +/- 0.24(stat) +/- 0.2(sys). The data are consistent with a point-like source. We provide a detailed description of the analysis results, and discuss physical mechanisms that may be responsible for the observed gamma-ray emission.

VERITAS is the high-sensitivity instrument of latest generation. It is often used for the short AGN monitoring exposures evenly distributed over entire observational season of a source of interest. Each of these exposures is long enough to detect the source at the flux level of about 1 Crab. During the 2009 observing season a number of exposures of Mkn 501 with VERITAS revealed variable TeV gamma-ray emission at the flux level eventually exceeding 2 Crab. The spectral and flux variability measurements in TeV gamma rays for the 2009 data sample of Mkn 501 are summarized in this paper.

There is increasing evidence that conventional cold dark matter (CDM) models lead to conflicts between observations and numerical simulations of dark matter halos on sub-galactic scales. Spergel and Steinhardt showed that if the CDM is strongly self-interacting, then the conflicts disappear. However, the assumption of strong self-interaction would rule out the favored candidates for CDM, namely weakly interacting massive particles (WIMPs), such as the neutralino. In this paper we propose a mechanism of non-thermal production of WIMPs and study its implications on the power spectrum. We find that the non-vanishing velocity of the WIMPs suppresses the power spectrum on small scales compared to what it obtained in the conventional CDM model. Our results show that, in this context, WIMPs as candidates for dark matter can work well both on large scales and on sub-galactic scales.

We re-examine approximations in the analytical calculation of the primordial spectrum of cosmological perturbation produced during inflation. Taking two inflation models (chaotic inflation and natural inflation) as examples, we numerically verify the accuracy of these approximations.