We present a first calculation of the rate for plasmon production in semiconductors from nuclei recoiling against dark matter. The process is analogous to bremsstrahlung of transverse photon modes, but with a longitudinal plasmon mode emitted instead. For dark matter in the 10 MeV - 1 GeV mass range, we find that the plasmon bremsstrahlung rate is 4-5 orders of magnitude smaller than that for elastic scattering, but 4-5 orders of magnitude larger than the transverse bremsstrahlung rate. Because the plasmon can decay into electronic excitations and has characteristic energy given by the plasma frequency $\omega_p$, with $\omega_p \approx 16$ eV in Si crystals, plasmon production provides a distinctive signature and new method to detect nuclear recoils from sub-GeV dark matter.

Dark matter could have a small electromagnetic charge, provided the charge-to-mass ratio is much less than that of electrons or protons. This candidate, commonly known as millicharged dark matter (mDM), would form a plasma and interact with the interstellar medium and electromagnetic fields within galaxies. In general, understanding the dynamics of mDM requires consideration of collective plasma effects. It has been proposed that mDM can be accelerated in supernova remnants, forming a dark cosmic ray population that would leave distinct experimental signatures. In this work, we study a microphysical model where mDM is shocked by a supernova remnant and isotropized in the frame of the expanding fluid. We find that for $|q_\chi/m_\chi| \gtrsim 10^{-13} e / {\rm MeV}$, the isotropization length for electromagnetic plasma instabilities is much shorter than the size of the supernova remnant. This is a necessary, though not sufficient, first step for formation of a Fermi-accelerated mDM component, and determining the size of this component requires further study. We discuss additional implications of mDM interactions in supernova remnants.

The radio source 3C 264, hosted by the giant elliptical galaxy NGC 3862, was observed with VERITAS between February 2017 and May 2019. These deep observations resulted in the discovery of very-high-energy (VHE; E $>100$ GeV) $\gamma$-ray emission from this active galaxy. An analysis of $\sim$57 hours of quality-selected live time yields a detection at the position of the source, corresponding to a statistical significance of 7.8 standard deviations above background. The observed VHE flux is variable on monthly time scales, with an elevated flux seen in 2018 observations. The VHE emission during this elevated state is well-characterized by a power-law spectrum with a photon index $\Gamma = 2.20 \pm 0.27$ and flux F($>315$ GeV) = ($7.6\pm 1.2_{\mathrm stat} \pm 2.3_{\mathrm syst})\times 10^{-13}$ cm$^{-2}$ s$^{-1}$, or approximately 0.7% of the Crab Nebula flux above the same threshold. 3C 264 ($z = 0.0217$) is the most distant radio galaxy detected at VHE, and the elevated state is thought to be similar to that of the famously outbursting jet in M 87. Consequently, extensive contemporaneous multi-wavelength data were acquired in 2018 at the time of the VHE high state. An analysis of these data, including VLBA, VLA, HST, Chandra and Swift observations in addition to the VERITAS data, is presented, along with a discussion of the resulting spectral energy distribution.

Compilation of papers presented by the VERITAS Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin.

We conduct a search for periodic emission in the very high-energy gamma-ray band (VHE; E > 100 GeV) from a total of 13 pulsars in an archival VERITAS data set with a total exposure of over 450 hours. The set of pulsars includes many of the brightest young gamma-ray pulsars visible in the Northern Hemisphere. The data analysis resulted in non-detections of pulsed VHE gamma rays from each pulsar. Upper limits on a potential VHE gamma-ray flux are derived at the 95% confidence level above three energy thresholds using two methods. These are the first such searches for pulsed VHE emission from each of the pulsars, and the obtained limits constrain a possible flux component manifesting at VHEs as is seen for the Crab pulsar.

These lectures provide an introduction to models and direct detection of dark matter. We summarize the general features and motivations for candidates in the full dark matter mass range, and then restrict to the keV-TeV mass window. Candidates in this window can be produced by thermal mechanisms in the standard cosmology, and are an important target for experimental searches. We then turn to sub-GeV dark matter (light dark matter) and dark sectors, an area where many new models and experiments are currently being proposed. We discuss the cosmology of dark sectors, specific portal realizations, and some of the prospects for detection. The final parts of these lectures focus on the theory for direct detection, both reviewing the fundamentals for nuclear recoils of WIMPs and describing new directions for sub-GeV candidates. A version of these lectures was originally presented at the TASI 2018 summer school on "Theory in an Era of Data".

The angular size of a star is a critical factor in determining its basic properties. Direct measurement of stellar angular diameters is difficult: at interstellar distances stars are generally too small to resolve by any individual imaging telescope. This fundamental limitation can be overcome by studying the diffraction pattern in the shadow cast when an asteroid occults a star, but only when the photometric uncertainty is smaller than the noise added by atmospheric scintillation. Atmospheric Cherenkov telescopes used for particle astrophysics observations have not generally been exploited for optical astronomy due to the modest optical quality of the mirror surface. However, their large mirror area makes them well suited for such high-time-resolution precision photometry measurements. Here we report two occultations of stars observed by the VERITAS Cherenkov telescopes with millisecond sampling, from which we are able to provide a direct measurement of the occulted stars' angular diameter at the $\leq0.1$ milliarcsecond scale. This is a resolution never achieved before with optical measurements and represents an order of magnitude improvement over the equivalent lunar occultation method. We compare the resulting stellar radius with empirically derived estimates from temperature and brightness measurements, confirming the latter can be biased for stars with ambiguous stellar classifications.

In a recent work, we numerically studied the radiative properties of the reverberation phase of pulsar wind nebulae (PWNe), i.e., when the reverse shock created by the supernova explosion travels back towards the pulsar, compressing the wind bubble. We focused on several well-characterized PWNe and used them as examples for introducing the concept of superefficiency. The latter is a period of the PWN evolution, happening within reverberation, where the luminosity in a given band exceeds the spin-down power at the time. Here, we explore a broad range of PWN models to study their reverberation and superefficiency phases in a systematic way. Armed with these models we consider two aspects: On the one hand, we analyze via Monte Carlo simulations how many Galactic PWNe are expected to be reverberating or in a superefficiency stage at any given time, providing the first such estimations. On the other hand, we focus on searching for observational signatures of such periods. We analyze archival observations and check for the existence of possible candidates for superefficient PWNe. We also provide predictions for the future evolution of the magnetar nebula J1834.9-0846 (which we consider to be starting its reverberation period) along the next 50 years. Using our simulations as input we study how sensitive current and future X-ray satellites (like eXTP or Athena) will be to observe such evolution, concluding that they will be able to track it in detail.

The analysis and combination of data from different gamma-ray instruments involves the use of collaboration proprietary software and case-by-case methods. The effort of defining a common data format for high-level data, namely event lists and instrument response functions (IRFs), has recently started for very-high-energy gamma-ray instruments, driven by the upcoming Cherenkov Telescope Array (CTA). In this work we implemented this prototypical data format for a small set of MAGIC, VERITAS, FACT, and H.E.S.S. Crab nebula observations, and we analyzed them with the open-source gammapy software package. By combining data from $Fermi$-LAT, and from four of the currently operating imaging atmospheric Cherenkov telescopes, we produced a joint maximum likelihood fit of the Crab nebula spectrum. Aspects of the statistical errors and the evaluation of systematic uncertainty are also commented upon, along with the release format of spectral measurements. The results presented in this work are obtained using open-access on-line assets that allow for a long-term reproducibility of the results.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Dark matter (DM) could couple to particles in the Standard Model (SM) through a light vector mediator. In the limit of small coupling, this portal could be responsible for producing the observed DM abundance through a mechanism known as freeze-in. Furthermore, the requisite DM-SM couplings provide a concrete benchmark for direct and indirect searches for DM. In this paper, we present updated calculations of the relic abundance for DM produced by freeze-in through a light vector mediator. We identify an additional production channel: the decay of photons that acquire an in-medium plasma mass. These plasmon decays are a dominant channel for DM production for sub-MeV DM masses, and including this channel leads to a significant reduction in the predicted signal strength for DM searches. Accounting for production from both plasmon decays and annihilations of SM fermions, the DM acquires a highly non-thermal phase space distribution which impacts the cosmology at later times; these cosmological effects will be explored in a companion paper.

We show that polar materials are excellent targets for direct detection of sub-GeV dark matter due to the presence of gapped optical phonons as well as acoustic phonons with high sound speed. We take the example of Gallium Arsenide (GaAs), which has the properties needed for experimental realization, and where many results can be estimated analytically. We find GaAs has excellent reach to dark photon absorption, can completely cover the freeze-in benchmark for scattering via an ultralight dark photon, and is competitive with other proposals to detect sub-MeV dark matter scattering off nuclei.

We report on observations of the pulsar / Be star binary system PSR J2032+4127 / MT91 213 in the energy range between 100 GeV and 20 TeV with the VERITAS and MAGIC imaging atmospheric Cherenkov telescope arrays. The binary orbit has a period of approximately 50 years, with the most recent periastron occurring on 2017 November 13. Our observations span from 18 months prior to periastron to one month after. A new, point-like, gamma-ray source is detected, coincident with the location of PSR J2032+4127 / MT91 213. The gamma-ray light curve and spectrum are well-characterized over the periastron passage. The flux is variable over at least an order of magnitude, peaking at periastron, thus providing a firm association of the TeV source with the pulsar / Be star system. Observations prior to periastron show a cutoff in the spectrum at an energy around 0.5 TeV. This result adds a new member to the small population of known TeV binaries, and it identifies only the second source of this class in which the nature and properties of the compact object are firmly established. We compare the gamma-ray results with the light curve measured with the X-ray Telescope (XRT) on board the Neil Gehrels \textit{Swift} Observatory and with the predictions of recent theoretical models of the system. We conclude that significant revision of the models is required to explain the details of the emission we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130.

We numerically study the radiative properties of the reverberation phase of pulsar wind nebulae. Reverberation brings a significant evolution in a short period of time. We show that even the Crab nebula, associated to the more energetic pulsar of the sample we consider, has a period in its future time evolution where the X-ray luminosity will exceed the spin-down power at the time. In fact, all nebulae in our sample are expected to have a period of radio, X-ray, and GeV superefficiency, and most will also have a period of TeV superefficiency. We analyze and characterize these superefficient phases.

If a component of the dark matter has dissipative interactions, it could collapse to form a thin dark disk in our Galaxy that is coplanar with the baryonic disk. It has been suggested that dark disks could explain a variety of observed phenomena, including periodic comet impacts. Using the first data release from the $Gaia$ space observatory, we search for a dark disk via its effect on stellar kinematics in the Milky Way. Our new limits disfavor the presence of a thin dark matter disk, and we present updated measurements on the total matter density in the Solar neighborhood.

The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.

Cosmic-ray electrons and positrons (CREs) at GeV-TeV energies are a unique probe of our local Galactic neighborhood. CREs lose energy rapidly via synchrotron radiation and inverse-Compton scattering processes while propagating within the Galaxy and these losses limit their propagation distance. For electrons with TeV energies, the limit is on the order of a kiloparsec. Within that distance there are only a few known astrophysical objects capable of accelerating electrons to such high energies. It is also possible that the CREs are the products of the annihilation or decay of heavy dark matter (DM) particles. VERITAS, an array of imaging air Cherenkov telescopes in southern Arizona, USA, is primarily utilized for gamma-ray astronomy, but also simultaneously collects CREs during all observations. We describe our methods of identifying CREs in VERITAS data and present an energy spectrum, extending from 300 GeV to 5 TeV, obtained from approximately 300 hours of observations. A single power-law fit is ruled out in VERITAS data. We find that the spectrum of CREs is consistent with a broken power law, with a break energy at 710 $\pm$ 40$_{stat}$ $\pm$ 140$_{syst}$ GeV.

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

HESS J1943+213 is a very-high-energy (VHE; $>$100 GeV) $\gamma$-ray source in the direction of the Galactic Plane. Studies exploring the classification of the source are converging towards its identification as an extreme synchrotron BL Lac object. Here we present 38 hours of VERITAS observations of HESS J1943+213 taken over two years. The source is detected with $\sim$20 standard deviations significance, showing a remarkably stable flux and spectrum in VHE $\gamma$-rays. Multi-frequency very-long-baseline array (VLBA) observations of the source confirm the extended, jet-like structure previously found in the 1.6 GHz band with European VLBI Network and detect this component in the 4.6 GHz and the 7.3 GHz bands. The radio spectral indices of the core and the jet and the level of polarization derived from the VLBA observations are in a range typical for blazars. Data from VERITAS, $Fermi$-LAT, $Swift$-XRT, FLWO 48$''$ telescope, and archival infrared and hard X-ray observations are used to construct and model the spectral energy distribution (SED) of the source with a synchrotron-self-Compton model. The well-measured $\gamma$-ray peak of the SED with VERITAS and $Fermi$-LAT provides constraining upper limits on the source redshift. Possible contribution of secondary $\gamma$-rays from ultra-high-energy cosmic ray-initiated electromagnetic cascades to the $\gamma$-ray emission is explored, finding that only a segment of the VHE spectrum can be accommodated with this process. A variability search is performed across X-ray and $\gamma$-ray bands. No statistically significant flux or spectral variability is detected.

We present results from deep observations towards the Cygnus region using 300 hours of very-high-energy (VHE) $\gamma$-ray data taken with the VERITAS Cherenkov telescope array and over seven years of high-energy $\gamma$-ray data taken with the Fermi satellite at an energy above 1 GeV. As the brightest region of diffuse $\gamma$-ray emission in the northern sky, the Cygnus region provides a promising area to probe the origins of cosmic rays. We report the identification of a potential Fermi-LAT counterpart to VER J2031+415 (TeV J2032+4130), and resolve the extended VHE source VER J2019+368 into two source candidates (VER J2018+367* and VER J2020+368*) and characterize their energy spectra. The Fermi-LAT morphology of 3FGL 2021.0+4031e (the Gamma-Cygni supernova remnant) was examined and a region of enhanced emission coincident with VER J2019+407 was identified and jointly fit with the VERITAS data. By modeling 3FGL J2015.6+3709 as two sources, one located at the location of the pulsar wind nebula CTB 87 and one at the quasar QSO J2015+371, a continuous spectrum from 1 GeV to 10 TeV was extracted for VER J2016+371 (CTB 87). An additional 71 locations coincident with Fermi-LAT sources and other potential objects of interest were tested for VHE $\gamma$-ray emission, with no emission detected and upper limits on the differential flux placed at an average of 2.3% of the Crab Nebula ux. We interpret these observations in a multiwavelength context and present the most detailed $\gamma$-ray view of the region to date.

We present the results of the analysis of 8 years of Fermi-LAT data of the pulsar/pulsar wind nebula complex PSR J0205+6449/3C 58. Using a contemporaneous ephemeris, we carried out a detailed analysis of PSR J0205+6449 both during its off-peak and on-peak phase intervals. 3C 58 is significantly detected during the off-peak phase interval. We show that the spectral energy distribution at high energies is the same disregarding the phases considered, and thus that this part of the spectrum is most likely dominated by the nebula radiation. We present results of theoretical models of the nebula and the magnetospheric emission that confirm this interpretation. Possible high-energy flares from 3C 58 were searched for, but none was unambiguously identified.

On 2015 March 23, VERITAS responded to a $Swift$-BAT detection of a gamma-ray burst, with observations beginning 270 seconds after the onset of BAT emission, and only 135 seconds after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40 minute integration corresponds to about 1% of the prompt fluence. Our limit is particularly significant since the very-high-energy (VHE) observation started only $\sim$2 minutes after the prompt emission peaked, and $Fermi$-LAT observations of numerous other bursts have revealed that the high-energy emission is typically delayed relative to the prompt radiation and lasts significantly longer. Also, the proximity of GRB~150323A ($z=0.593$) limits the attenuation by the extragalactic background light to $\sim 50$ % at 100-200 GeV. We conclude that GRB 150323A had an intrinsically very weak high-energy afterglow, or that the GeV spectrum had a turnover below $\sim100$ GeV. If the GRB exploded into the stellar wind of a massive progenitor, the VHE non-detection constrains the wind density parameter to be $A\gtrsim 3\times 10^{11}$ g cm$^{-1}$, consistent with a standard Wolf-Rayet progenitor. Alternatively, the VHE emission from the blast wave would be weak in a very tenuous medium such as the ISM, which therefore cannot be ruled out as the environment of GRB 150323A.

Combined with very-long-baseline interferometry measurements, the observations of fast TeV gamma-ray flares probe the structure and emission mechanism of blazar jets. However, only a handful of such flares have been detected to date, and only within the last few years have these flares been observed from lower-frequency-peaked BL~Lac objects and flat-spectrum radio quasars. We report on a fast TeV gamma-ray flare from the blazar BL~Lacertae observed by VERITAS, with a rise time of $\sim$2.3~hr and a decay time of $\sim$36~min. The peak flux above 200 GeV is $(4.2 \pm 0.6) \times 10^{-6} \;\text{photon} \;\text{m}^{-2}\; \text{s}^{-1}$ measured with a 4-minute-binned light curve, corresponding to $\sim$180\% of the flux which is observed from the Crab Nebula above the same energy threshold. Variability contemporaneous with the TeV gamma-ray flare was observed in GeV gamma-ray, X-ray, and optical flux, as well as in optical and radio polarization. Additionally, a possible moving emission feature with superluminal apparent velocity was identified in VLBA observations at 43 GHz, potentially passing the radio core of the jet around the time of the gamma-ray flare. We discuss the constraints on the size, Lorentz factor, and location of the emitting region of the flare, and the interpretations with several theoretical models which invoke relativistic plasma passing stationary shocks.

On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC observatory were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns.

Imaging Atmospheric Cherenkov Telescopes study the highest energy (up to tens of TeV) photon emission coming from nearby and distant astrophysical sources, thus providing valuable results from searches for Lorentz Invariance Violation (LIV) effects. Highly variable, energetic and distant sources such as Pulsars and AGNs are the best targets for the Time-of-Flight LIV studies. However, the limited number of observations of AGN flares or of high-energy pulsed emission greatly restricts the potential of such studies, especially any potential LIV effects as a function of redshift. To address these issues, an inter-experiment working group has been established by the three major collaborations taking data with Imaging Atmospheric Cherenkov Telescopes (H.E.S.S., MAGIC and VERITAS) with the aim to increase sensitivity to any effects of LIV, together with an improved control of systematic uncertainties, by sharing data samples and developing joint analysis methods. This will allow an increase in the number of available sources and to perform a sensitive search for redshift dependencies. This presentation reviews the first combined maximum likelihood method analyses using simu- lations of published source observations done in the past with H.E.S.S., MAGIC and VERITAS. The results from analyses based on combined maximum likelihood methods, the strategies to deal with data from different types of sources and instruments, as well as future plans will be presented.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

We present constraints on the annihilation cross section of WIMP dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We report on the results of $\sim$230 hours of observations of five dwarf galaxies and the joint statistical analysis of four of the dwarf galaxies. We find no evidence of gamma-ray emission from any individual dwarf nor in the joint analysis. The derived upper limit on the dark matter annihilation cross section from the joint analysis is $1.35\times 10^{-23} {\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the bottom quark ($b\bar{b}$) final state, $2.85\times 10^{-24}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the tau lepton ($\tau^{+}\tau^{-}$) final state and $1.32\times 10^{-25}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the gauge boson ($\gamma\gamma$) final state.

Semiconductors are by now well-established targets for direct detection of MeV to GeV dark matter via scattering off electrons. We show that semiconductor targets can also detect significantly lighter dark matter via an absorption process. When the dark matter mass is above the band gap of the semiconductor (around an eV), absorption proceeds by excitation of an electron into the conduction band. Below the band gap, multi-phonon excitations enable absorption of dark matter in the 0.01 eV to eV mass range. Energetic dark matter particles emitted from the sun can also be probed for masses below an eV. We derive the reach for absorption of a relic kinetically mixed dark photon or pseudoscalar in germanium and silicon, and show that existing direct detection results already probe new parameter space. With only a moderate exposure, low-threshold semiconductor target experiments can exceed current astrophysical and terrestrial constraints on sub-keV bosonic dark matter.

We present a search for magnetically broadened gamma-ray emission around active galactic nuclei (AGN), using VERITAS observations of seven hard-spectrum blazars. A cascade process occurs when multi-TeV gamma rays from AGN interact with extragalactic background light (EBL) photons to produce electron-positron pairs, which then interact with cosmic microwave background (CMB) photons via inverse-Compton scattering to produce gamma rays. Due to the deflection of the electron-positron pairs, a non-zero intergalactic magnetic field (IGMF) would potentially produce detectable effects on the angular distribution of the cascade emission. In particular, an angular broadening compared to the unscattered emission could occur. Through non-detection of angularly broadened emission from 1ES 1218+304, the source with the largest predicted cascade fraction, we exclude a range of IGMF strengths around $10^{-14}$G at the 95% confidence level. The extent of the exclusion range varies with the assumptions made about the intrinsic spectrum of 1ES 1218+304 and the EBL model used in the simulation of the cascade process. All of the sources are used to set limits on the flux due to extended emission.

We present time-resolved spectral analysis of the steep decay segments of 29 bright X-ray flares of gamma-ray bursts (GRBs) observed with the Swift/X-ray telescope, and model their lightcurves and spectral index evolution behaviors with the curvature effect model. Our results show that the observed rapid flux decay and strong spectral index evolution with time can be well fit with this model, and the derived characteristic timescales ($t_c$) are in the range of $33\sim 264$ seconds. Using an empirical relation between the peak luminosity and the Lorentz factor derived from the prompt gamma-rays, we estimate the Lorentz factors of the flares ($\Gamma_{\rm X}$). We obtain $\Gamma_{\rm X}=17\sim 87$ with a median value of $52$, which is smaller than the initial Lorentz factors of prompt gamma-ray fireballs. With the derived $t_c$ and $\Gamma_{\rm X}$, we constrain the radiating regions of 13 X-ray flares, yielding $R_{\rm X}=(0.2\sim 1.1)\times 10^{16}$ cm, which are smaller than the radii of the afterglow fireballs at the peak times of the flares. A long evolution feature from prompt gamma-ray phase to the X-ray epoch is found by incorporating our results with a sample of GRBs whose initial Lorentz factors are available in literatures, i.e., $\Gamma\propto [t_{p}/(1+z)]^{-0.69\pm 0.06}$. These results may shed lights on the long term evolution of GRB central engines.

Imaging atmospheric Cherenkov telescopes (IACTs) are sensitive to rare gamma-ray photons, buried in the background of charged cosmic-ray (CR) particles, the flux of which is several orders of magnitude greater. The ability to separate gamma rays from CR particles is important, as it is directly related to the sensitivity of the instrument. This gamma-ray/CR-particle classification problem in IACT data analysis can be treated with the rapidly-advancing machine learning algorithms, which have the potential to outperform the traditional box-cut methods on image parameters. We present preliminary results of a precise classification of a small set of muon events using a convolutional neural networks model with the raw images as input features. We also show the possibility of using the convolutional neural networks model for regression problems, such as the radius and brightness measurement of muon events, which can be used to calibrate the throughput efficiency of IACTs.

We examine in depth a recent proposal to utilize superfluid helium for direct detection of sub-MeV mass dark matter. For sub-keV recoil energies, nuclear scattering events in liquid helium primarily deposit energy into long-lived phonon and roton quasiparticle excitations. If the energy thresholds of the detector can be reduced to the meV scale, then dark matter as light as ~MeV can be reached with ordinary nuclear recoils. If, on the other hand, two or more quasiparticle excitations are directly produced in the dark matter interaction, the kinematics of the scattering allows sensitivity to dark matter as light as ~keV at the same energy resolution. We present in detail the theoretical framework for describing excitations in superfluid helium, using it to calculate the rate for the leading dark matter scattering interaction, where an off-shell phonon splits into two or more higher-momentum excitations. We validate our analytic results against the measured and simulated dynamic response of superfluid helium. Finally, we apply this formalism to the case of a kinetically mixed hidden photon in the superfluid, both with and without an external electric field to catalyze the processes.

This report, based on the Dark Sectors workshop at SLAC in April 2016, summarizes the scientific importance of searches for dark sector dark matter and forces at masses beneath the weak-scale, the status of this broad international field, the important milestones motivating future exploration, and promising experimental opportunities to reach these milestones over the next 5-10 years.

Transient X-ray binaries produce major outbursts in which the X-ray flux can increase over the quiescent level by factors as large as $10^7$. The low-mass X-ray binary V 404 Cyg and the high-mass system 4U 0115+634 underwent such major outbursts in June and October 2015, respectively. We present here observations at energies above hundreds of GeV with the VERITAS observatory taken during some of the brightest X-ray activity ever observed from these systems. No gamma-ray emission has been detected by VERITAS in 2.5 hours of observations of the microquasar V 404 Cyg from 2015, June 20-21. The upper flux limits derived from these observations on the gamma-ray flux above 200 GeV of F $< 4.4\times 10^{-12}$ cm$^{-2}$ s$^{-1}$ correspond to a tiny fraction (about $10^{-6}$) of the Eddington luminosity of the system, in stark contrast to that seen in the X-ray band. No gamma rays have been detected during observations of 4U 0115+634 in the period of major X-ray activity in October 2015. The flux upper limit derived from our observations is F $< 2.1\times 10^{-12}$ cm$^{-2}$ s$^{-1}$ for gamma rays above 300 GeV, setting an upper limit on the ratio of gamma-ray to X-ray luminosity of less than 4%.

Motivated by gauge coupling unification and dark matter, we present an extension to the Standard Model where both are achieved by adding an extra new matter multiplet. Such considerations lead to a Grand Unified Theory with very heavy WIMPzilla dark matter, which has mass greater than ~10^7 GeV and must be produced before reheating ends. Naturally, we refer to this scenario as GUTzilla dark matter. Here we present a minimal GUTzilla model, adding a vector-like quark multiplet to the Standard Model. Proton decay constraints require the new multiplet to be both color and electroweak charged, which prompts us to include a new confining SU(3) gauge group that binds the multiplet into a neutral composite dark matter candidate. Current direct detection constraints are evaded due to the large dark matter mass; meanwhile, next-generation direct detection and proton decay experiments will probe much of the parameter space. The relic abundance is strongly dependent on the dynamics of the hidden confining sector, and we show that dark matter production during the epoch of reheating can give the right abundance.

Superconducting targets have recently been proposed for the direct detection of dark matter as light as a keV, via elastic scattering off conduction electrons in Cooper pairs. Detecting such light dark matter requires sensitivity to energies as small as the superconducting gap of O(meV). Here we show that these same superconducting devices can detect much lighter DM, of meV to eV mass, via dark matter absorption on a conduction electron, followed by emission of an athermal phonon. We demonstrate the power of this setup for relic kinetically mixed hidden photons, pseudoscalars, and scalars, showing the reach can exceed current astrophysical and terrestrial constraints with only a moderate exposure.

The F-type star KIC 8462852 has recently been identified as an exceptional target for SETI (search for extraterrestrial intelligence) observations. We describe an analysis methodology for optical SETI, which we have used to analyse nine hours of serendipitous archival observations of KIC 8462852 made with the VERITAS gamma-ray observatory between 2009 and 2015. No evidence of pulsed optical beacons, above a pulse intensity at the Earth of approximately 1 photon per m^2, is found. We also discuss the potential use of imaging atmospheric Cherenkov telescope arrays in searching for extremely short duration optical transients in general.

Observations of very early multi-wavelength afterglows are critical to reveal the properties of the radiating fireball and its environment as well as the central engine of gamma-ray bursts (GRBs). We report our optical observations of GRB 111228A from 95 sec to about 50 hours after the burst trigger and investigate its properties of the prompt gamma-rays and the ambient medium using our data and the data observed with {\em Swift} and {\em Fermi} missions. Our joint optical and X-ray spectral fits to the afterglow data show that the ambient medium features as low dust-to-gas ratio. Incorporating the energy injection effect, our best fit to the afterglow lightcurves with the standard afterglow model via the Markov Chain Monte Carlo (MCMC) technique shows that $\epsilon_e=(6.9\pm 0.3)\times 10^{-2}$, $\epsilon_B=(7.73\pm 0.62)\times 10^{-6}$, $E_{\rm K}=(6.32\pm 0.86)\times 10^{53}\rm erg$, $n=0.100\pm 0.014$ cm$^{-3}$. The low medium density likely implies that the afterglow jet may be in a halo or in a hot ISM. Achromatic shallow decay segment observed in the optical and X-ray bands is well explained with the long-lasting energy injection from the central engine, which would be a magnetar with a period of about 1.92 ms inferred from the data. The $E_p$ of its time-integrated prompt gamma-ray spectrum is $\sim 26$ KeV. Using the initial Lorentz factor ($\Gamma_0=476^{+225}_{-237}$) derived from our afterglow model fit, it is found that GRB 111228A satisfies the $L_{\rm iso}-E_{\rm p,z}-\Gamma_0$ relation and bridges the typical GRBs and low luminosity GRBs in this relation.

We select a sample of 34 gamma-ray bursts (GRBs) whose $\Gamma_0$ values are derived with the onset peaks observed in the afterglow lightcurves (except for GRB 060218 whose $\Gamma_0$ is estimated with its radio data), and investigate the correlations among $\Gamma_0$, the isotropic peak luminosity ($L_{\rm iso}$), and the peak energy ($E_{\rm p,z}$) of the $\nu f_\nu$ spectrum in the cosmological rest frame. An analysis of pair correlations among these observables well confirms the results reported by the previous papers. More interestingly, a tight correlation among $L_{\rm iso}$, $E_{\rm p,z}$, and $\Gamma_0$ is found from a multiple regression analysis, which takes the form of $L_{\rm iso} \propto E_{\rm p,z}^{1.34\pm 0.14} \Gamma_0^{1.32\pm 0.19}$ or $E_{\rm p,z} \propto L_{\rm iso}^{0.55\pm 0.06}\Gamma_0^{-0.50\pm 0.17}$. Nine other GRBs whose $\Gamma_0$ are derived via the pair production opacity constraint also follow such a correlation. Excluding GRB 060218, the $L_{\rm iso}-E_{\rm p,z}-\Gamma_0$ correlation is valid, and it even holds in the jet co-moving frame. However, GRB 060218 deviates the $L^{'}_{\rm iso}-E^{'}_{\rm p}$ relation of typical GRBs in the jet co-moving frame with $3\sigma$. We argue that the $L_{\rm iso} - E_{\rm p, z} - \Gamma_0$ correlation may be more physical than the $L_{\rm iso} - E_{\rm p,z}$ correlation, since physically the relationship between the observed $L_{\rm iso}$ and $E_{\rm p,z}$ not only depends on radiation physics, but also depends on the bulk motion of the jet. We explore the possible origins of this correlation and discuss its physical implications for understanding GRB jet composition and radiation mechanism.

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

A fermion dark matter candidate with a relic abundance set by annihilation through a pseudoscalar can evade constraints from direct detection experiments. We present simplified models that realize this fact by coupling a fermion dark sector to a two-Higgs doublet model. These models are generalizations of mixed bino-Higgsino dark matter in the MSSM, with more freedom in the couplings and scalar spectra. Annihilation near a pseudoscalar resonance allows a significant amount of parameter space for thermal relic dark matter compared to singlet-doublet dark matter, in which the fermions couple only to the SM Higgs doublet. In a general two-Higgs doublet model, there is also freedom for the pseudoscalar to be relatively light and it is possible to obtain thermal relic dark matter candidates even below 100 GeV. In particular, we find ample room to obtain dark matter with mass around 50 GeV and fitting the Galactic Center excess in gamma-rays. This region of parameter space can be probed by LHC searches for heavy pseudoscalars or electroweakinos, and possibly by other new collider signals.

Thermal relic dark matter particles with a mass of 31-40 GeV and that dominantly annihilate to bottom quarks have been shown to provide an excellent description of the excess gamma rays observed from the center of the Milky Way. Flavored dark matter provides a well-motivated framework in which the dark matter can dominantly couple to bottom quarks in a flavor-safe manner. We propose a phenomenologically viable model of bottom flavored dark matter that can account for the spectral shape and normalization of the gamma-ray excess while naturally suppressing the elastic scattering cross sections probed by direct detection experiments. This model will be definitively tested with increased exposure at LUX and with data from the upcoming high-energy run of the Large Hadron Collider (LHC).

Dark-matter annihilation into electroweak gauge and Higgs bosons results in $\gamma$-ray emission. We use observational upper limits on the fluxes of both line and continuum $\gamma$-rays from the Milky Way Galactic Center and from Milky Way dwarf companion galaxies to set exclusion limits on allowed dark-matter masses. (Generally, Galactic Center $\gamma$-ray line search limits from the Fermi-LAT and the H.E.S.S. experiments are most restrictive.) Our limits apply under the following assumptions: a) the dark matter species is a cold thermal relic with present mass density equal to the measured dark-matter density of the universe; b) dark-matter annihilation to standard-model particles is described in the non-relativistic limit by a single effective operator ${\cal O} \propto J_{DM}\cdot J_{SM}$, where $J_{DM}$ is a standard-model singlet current consisting of dark-matter fields (Dirac fermions or complex scalars), and $J_{SM}$ is a standard-model singlet current consisting of electroweak gauge and Higgs bosons; and c) the dark-matter mass is in the range 5 GeV to 20 TeV. We consider, in turn, the 34 possible operators with mass dimension 8 or lower with non-zero s-wave annihilation channels satisfying the above assumptions. Our limits are presented in a large number of figures, one for each of the 34 possible operators; these limits can be grouped into 13 classes determined by the field content and structure of the operators. We also identify three classes of operators (coupling to the Higgs and $SU(2)_L$ gauge bosons) that can supply a 130 GeV line with the desired strength to fit the putative line signal in Fermi data, while saturating the relic density and satisfying all other indirect constraints we consider.

We examine cosmological, astrophysical and collider constraints on thermal dark matter (DM) with mass mX in the range 1 MeV to 10 GeV. Cosmic microwave background (CMB) observations, which severely constrain light symmetric DM, can be evaded if the DM relic density is sufficiently asymmetric. CMB constraints require the present anti-DM to DM ratio to be less than 2*10^{-6} (10^{-1}) for DM mass mX = 1 MeV (10 GeV) with ionizing efficiency factor f ~ 1. We determine the minimum annihilation cross section for achieving these asymmetries subject to the relic density constraint; these cross sections are larger than the usual thermal annihilation cross section. On account of collider constraints, such annihilation cross sections can only be obtained by invoking light mediators. These light mediators can give rise to significant DM self-interactions, and we derive a lower bound on the mediator mass from elliptical DM halo shape constraints. We find that halo shapes require a mediator with mass mphi > 4 * 10^{-2} MeV (40 MeV) for mX = 1 MeV (10 GeV). We map all of these constraints to the parameter space of DM-electron and DM-nucleon scattering cross sections for direct detection. For DM-electron scattering, a significant fraction of the parameter space is already ruled out by beam-dump and supernova cooling constraints.

High-precision measurements of the temperature and polarization anisotropies of the cosmic microwave background radiation have been previously employed to set robust constraints on dark matter annihilation during recombination. In this work we improve and generalize these constraints to apply to energy deposition during the recombination era with arbitrary redshift dependence. Our approach also provides more rigorous and model-independent bounds on dark matter annihilation and decay scenarios. We employ principal component analysis to identify a basis of weighting functions for the energy deposition. The coefficients of these weighting functions parameterize any energy deposition model and can be constrained directly by experiment. For generic energy deposition histories that are currently allowed by WMAP7 data, up to 3 principal component coefficients are measurable by Planck and up to 5 coefficients are measurable by an ideal cosmic variance limited experiment. For WIMP dark matter, our analysis demonstrates that the effect on the CMB is described well by a single (normalization) parameter and a "universal" redshift dependence for the energy deposition history. We give WMAP 7 constraints on both generic energy deposition histories and the universal WIMP case.

New laboratory measurements using an Electron Beam Ion Trap (EBIT) and an x-ray microcalorimeter are presented for the n=3 to n=2 Fe XVII emission lines in the 15 {\AA} to 17 {\AA} range, along with new theoretical predictions for a variety of electron energy distributions. This work improves upon our earlier work on these lines by providing measurements at more electron impact energies (seven values from 846 to 1185 eV), performing an in situ determination of the x-ray window transmission, taking steps to minimize the ion impurity concentrations, correcting the electron energies for space charge shifts, and estimating the residual electron energy uncertainties. The results for the 3C/3D and 3s/3C line ratios are generally in agreement with the closest theory to within 10%, and in agreement with previous measurements from an independent group to within 20%. Better consistency between the two experimental groups is obtained at the lowest electron energies by using theory to interpolate, taking into account the significantly different electron energy distributions. Evidence for resonance collision effects in the spectra is discussed. Renormalized values for the absolute cross sections of the 3C and 3D lines are obtained by combining previously published results, and shown to be in agreement with the predictions of converged R-matrix theory. This work establishes consistency between results from independent laboratories and improves the reliability of these lines for astrophysical diagnostics. Factors that should be taken into account for accurate diagnostics are discussed, including electron energy distribution, polarization, absorption/scattering, and line blends.

The magnetic inelastic dark matter (MiDM) model, in which dark matter inelastically scatters off nuclei through a magnetic dipole interaction, has previously been shown to reconcile the DAMA/LIBRA annual modulation signal with null results from other experiments. In this work, we explore the unique directional detection signature of MiDM. After the dark matter scatters into its excited state, it decays with a lifetime of order 1 microsecond and emits a photon with energy ~100 keV. Both the nuclear recoil and the corresponding emitted photon can be detected by studying delayed coincidence events. The recoil track and velocity of the excited state can be reconstructed from the nuclear interaction vertex and the photon decay vertex. The angular distribution of the WIMP recoil tracks is sharply peaked and modulates daily. It is therefore possible to observe the directional modulation of WIMP-nucleon scattering without a large-volume gaseous directional detection experiment. Furthermore, current experiments such as XENON100 can immediately measure this directional modulation and constrain the MiDM parameter space with an exposure of a few thousand kg day.

The DAMA/LIBRA collaboration has detected an annual modulation of the recoil rate in NaI crystals with the phase expected for WIMP scattering events. This signal is dramatically inconsistent with upper limits from other experiments for elastically scattering weak-scale WIMPs. However, the results are compatible for the case of inelastic dark matter (iDM). The iDM theory, as implemented by Tucker-Smith and Weiner, constrains the WIMP to a tight contour in sigma_n-delta space, where delta is the mass difference between the ground state and excited WIMPs. An urgent priority in direct detection is to test this scenario. The crucial test of the iDM explanation of DAMA -- an "experimentum crucis" -- is an experiment with directional sensitivity, which can measure the daily modulation in direction. Because the contrast can be 100%, it is a sharper test than the much smaller annual modulation in the rate. We estimate the significance of such an experiment as a function of the WIMP mass, cross section, background rate, and other parameters. The proposed experiment severely constrains the DAMA/iDM scenario even with modest exposure (~1000 kg day) on gaseous xenon.

Recent cosmic ray, gamma ray, and microwave signals observed by Fermi, PAMELA, and WMAP indicate an unexpected primary source of e+e- at 10-1000 GeV. We fit these data to "standard backgrounds" plus a new source, assumed to be a separable function of position and energy. For the spatial part, we consider three cases: annihilating dark matter, decaying dark matter, and pulsars. In each case, we use GALPROP to inject energy in log-spaced energy bins and compute the expected cosmic-ray and photon signals for each bin. We then fit a linear combination of energy bins, plus backgrounds, to the data. We use a non-parametric fit, with no prior constraints on the spectrum except smoothness and non-negativity. In addition, we consider arbitrary modifications to the energy spectrum of the "ordinary" primary source function, fixing its spatial part, finding this alone to be inadequate to explain the PAMELA or WMAP signals. We explore variations in the fits due to choice of magnetic field, primary electron injection index, spatial profiles, propagation parameters, and fit regularization method. Dark matter annihilation fits well, where our fit finds a mass of ~1 TeV and a boost factor times energy fraction of ~70. While it is possible for dark matter decay and pulsars to fit the data, unconventionally high magnetic fields and radiation densities are required near the Galactic Center to counter the relative shallowness of the assumed spatial profiles. We also fit to linear combinations of these three scenarios, though the fit is much less constrained.

Context: Sh2-104 is a Galactic H ii region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observational confirmation of the collect-and-collapse model of triggered star-formation. Aims: We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In addition, we investigate the relationship between the dust emissivity index {\beta} and the dust temperature, T_dust. Methods: Using Herschel PACS and SPIRE images at 100, 160, 250, 350 and 500 {\mu}m we determine T_dust and {\beta} throughout Sh2-104, fitting the spectral energy distributions (SEDs) obtained from aperture photometry. With the SPIRE Fourier transform spectrometer (FTS) we obtained spectra at different positions in the Sh2-104 region. We detect J-ladders of CO and 13CO, with which we derive the gas temperature and column density. We also detect proxies of ionizing flux as the [NII] 3P1-3P0 and [CI] 3P2-3P1 transitions. Results: We find an average value of {\beta} ~ 1.5 throughout Sh2-104, as well as a T dust difference between the photodissociation region (PDR, ~ 25 K) and the interior (~ 40 K) of the bubble. We recover the anti-correlation between {\beta} and dust temperature reported numerous times in the literature. The relative isotopologue abundances of CO appear to be enhanced above the standard ISM values, but the obtained value is very preliminary and is still affected by large uncertainties.

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M.

Project GRAND is an extensive air shower array of proportional wire chambers. It has 64 stations in a 100m x 100m area; each station has eight planes of proportional wire chambers with a 50 mm steel absorber plate above the bottom two planes. This arrangement of planes, each 1.25 square meters of area, allow an angular measurement for each track to 0.25 degrees in each of two projections. The steel absorber plate allows a measurement of the identity of each muon track to 96% accuracy. Two data-taking triggers allow data to be simultaneously taken for a) extensive air showers (multiple coincidence station hits) at about 1 Hz and b) single muons (single tracks of identified muons) at 2000 Hz. Eight on-line computers pre-analyze the single track data and store the results on magnetic tape in compacted form with a minimum of computer dead-time. One additional computer reads data from the shower triggers and records this raw data on a separate magnetic tape with no pre-analysis.

An increased data sample of identified secondary muons is collected at detection level. Project GRAND identifies secondary cosmic ray muons from electrons utilizing a thin steel absorber and tracking (PWC) chambers. The resulting angular resolution for the primary is about +/- 5 degrees. Since there is a 1.5% probability that a 100 GeV gamma primary will produce a single muon track at ground level, the high statistics allows a search for angular enhancements due to gammas. More than 100 billion muons are identified and analyzed. A table of stellar sources is examined.

The high energy gamma ray sources provided by experiment EGRET on the Compton Gamma Ray Observatory are examined, testing Project GRAND's ability to detect primary gamma rays by means of studying secondary muons. There is ~1.5% chance that a 100 GeV primary cosmic gamma ray will produce a muon at detection level (Fasso & Poirier, 1999), a probability which increases with increasing gamma energy. Project GRAND has 80 square meters of muon detector which identifies secondary muons >0.1 GeV and measures their angles to 0.26 deg (projected angle in the XZ and YZ planes). Data taken during the last two years are analyzed. A table of EGRET's gamma ray sources is examined (The Third EGRET Catalog of High Energy Gamma Ray Sources (EGRET webpage, 1999). EGRET's flux (>0.1 GeV), angular information, and spectral index were extrapolated to GRAND's energy region (>30 GeV). Then the geometrical acceptance was calculated for each of these events. Thus, GRAND's sensitivity to each of EGRET's sources was predicted. A product of extrapolated flux and GRAND's detection sensitivity yields an overall parameter indicating GRAND's relative sensitivity to each source, assuming an energy dependence given by EGRET's spectral index. The 14 sources with the best extrapolated detection efficiency were selected to be examined with the data.

Project GRAND presents the results of a search for coincident high-energy gamma ray events in the direction and at the time of nine Gamma Ray Bursts (GRBs) detected by BATSE. A gamma ray has a non-negligible hadron production cross section; for each gamma ray of energy of 100 GeV, there are 0.015 muons which reach detection level (Fasso & Poirier, 1999). These muons are identified and their angles are measured in stations of eight planes of proportional wire chambers (PWCs). A 50 mm steel plate above the bottom pair of planes is used to distinguish muons from electrons. The mean angular resolution is 0.26 degrees over a +/- 61 degree range in the XZ and YZ planes. The BATSE GRB catalogue is examined for bursts which are near zenith for Project GRAND. The geometrical acceptance is calculated for each of these events. The product is then taken of the GRB flux and GRAND's geometrical acceptance. The nine sources with the best combination of detection efficiency and BATSE's intensity are selected to be examined in the data. The most significant detection of these nine sources is at a statistical significance of +3.7 sigma; this is also the GRB with the highest product of GRB flux and geometrical acceptance.

Project GRAND presents results on the atomic composition of primary cosmic rays. This is accomplished by determining the average height of primary particles that cause extensive air showers detected by Project GRAND. Particles with a larger cross sectional area, such as iron nuclei, are likely to start an extensive air shower higher in the atmosphere whereas protons, with a smaller cross section, would pass through more air before interacting and thus start showers at lower heights. Such heights can be determined by extrapolating identified muon tracks backward (upward) to determine their height of origin (Gress et al., 1997). Since muons are from the top, hadronic part of the shower, they are a good estimator for the beginning of the shower. The data for this study were taken during the previous year with 20 million shower events.

Project GRAND preents a map in right ascension and declination for single cosmic secondary tracks which have been identified as muons. These muons are measured in stations of proportional wire chambers (PWCs) which have a mean angular resolution of 0.26 deg. A 50 mm steel plate interspersed with the PWCs is used to distinguish muons from electrons. Two years of data are accumulated in a map of right ascension and declination which covers all of the northern hemisphere and a small portion of the southern hemisphere. The total number of muon tracks in this data sample are 40 billion. Significant large angular structure is found in this map; the deviations from average are small compared to the average, but because of the high statistics, the deviations from average are statistically significant to many standard deviations. Systematic effects in the data are minimized by averaging over this two year period. These results are similar to those obtained in previous data.

Project GRAND measures the asymmetries of secondary cosmic muons by calculating the differences between the muons from the east and the muons from the west for complete solar days. A similar north-south difference is also presented. These muons are identified in single track data taken with the 64 stations of Project GRAND; the data cover a two year span with a total muon sample of 30 billion muons.

Project GRAND is an extensive air shower array utilizing position sensitive detectors of proportional wire chambers. The 64 detectors deployed in a field 100 m x 100 m are located at 86.2 deg W and 41.7 deg N, at 220 m above sea level. The project was completed about two years ago and has been taking data, simultaneously on two triggers: 1) multiple-hit coincidence triggers which collect data on extensive air showers, and 2) 200 Hz triggers which collect the single tracks stored in each station during the last 5 msec.

We present nine color CCD intermediate-band spectrophotometry of a two square degree field centered on the old open cluster M67, from 3890$\rm \AA$ to nearly 1$\mu$. These observations are taken as a part of the BATC (Beijing-Arizona-Taipei-Connecticut) Color Survey of the Sky, for both scientific and calibration reasons. With these data we show that the BATC survey can reach its goal of obtaining spectrophotometry to a zero point accuracy of 0.01 mag, and down to V = 21 with 0.3 mag random error. We fit the color-magnitude diagrams (CMDs) with Worthey's theoretical models. The net result is the excellent fit of the 4.0 Gyr, [Fe/H] = $-0.10$ model to our data, including a good fit to the main sequence (MS) turn-off. Our data are consistent with a toy model with 50\% of the stars in M67 being binaries and a random distribution of binary mass-ratios, although other models with different mass-ratio distributions cannot be ruled out. The spatial distribution and mass function (MF) of stars in M67 show marked effects of dynamical evolution and evaporation of stars from the cluster. Blue stragglers and binary stars are the most condensed within the cluster, with degree of condensation depending on mass.We find M67 to have an elongated shape, oriented at an angle of $15^{\circ}$ relative to the galactic plane. Within its tidal radius, the observed MF of M67 between 1.2 $\rm M_\odot$ and $\rm 0.8 M_\odot$ has a Salpeter slope $\rm \eta = -1.93 \pm 0.66$. For stars of mass below 0.8 $\rm M_\odot$, $\rm \eta \sim 0$. It is plausible that the leveling-off of the MF at lower masses is a result of evaporation of lower mass stars in this mass range at a rate of one every $\sim 10^7$ years. If so, it is plausible that the IMF of M67 has the canonical field value of $\rm \eta = -2.0$.