In 2018-2019 the VERITAS VHE gamma-ray observatory was augmented with highspeed optical instrumentation and continuous data recording electronics to create a sensitive Stellar Intensity Interferometry (SII) observatory, VERITAS-SII. The primary science goal of VERITAS-SII is to perform stellar diameter measurements and image analysis in the visible wavebands on a selection of bright (m< 6), hot (O/B/A) stars. The VERITAS Collaboration has agreed to the deployment and operation of VERITAS-SII during several days each month around the full moon period when VERITAS does not perform VHE gamma-ray observations. The VERITAS-SII augmentation employs custom high-speed/low-noise focal plane instrumentation using high quantum efficiency photomultiplier tubes, and a battery-powered, fiber-optic controlled High Voltage supply. To reduce engineering time, VERITAS-SII uses commercially available high-speed (250 MS/sec), continuously streaming electronics to record the time dependence of the intensity fluctuations at each VERITAS telescope. VERITAS-SII also uses fast ( < 100 psec) data acquisition clock synchronization over inter-telescope distances (greater than 100 m) using a commercially available White Rabbit based timing solution. VERITAS-SII is now in full operation at the VERITAS observatory, F.L.Whipple Observatory, Amado, AZ USA. This paper describes the design of the instrumentation hardware used for VERITAS-SII augmentation of the VERITAS observatory, the status of initial VERITAS-SII observations, and plans for future improvements to VERITAS-SII.

The discovery of the electromagnetic counterparts to the binary neutron star merger GW170817 has opened the era of GW+EM multi-messenger astronomy. Exploiting this breakthrough requires increasing samples to explore the diversity of kilonova behaviour and provide more stringent constraints on the Hubble constant, and tests of fundamental physics. LSST can play a key role in this field in the 2020s, when the gravitational wave detector network is expected to detect higher rates of merger events involving neutron stars ($\sim$10s per year) out to distances of several hundred Mpc. Here we propose comprehensive target-of-opportunity (ToOs) strategies for follow-up of gravitational-wave sources that will make LSST the premiere machine for discovery and early characterization for neutron star mergers and other gravitational-wave sources.

This project aims at exploiting the wide-field and limiting-magnitude capabilities of the LSST to fully characterise the resolved stellar populations in/around six Local Group stellar systems of different morphological type at ~30 to ~400 kpc distance from us. We selected targets that host red giant branch (RGB) stars which are within the reach of Gaia and not yet (all) saturated with the LSST. We will use RR Lyrae stars, Cepheids, SX Phoenicis, delta Scuti stars and Long Period Variables, along with the Color Magnitude Diagram of the resolved stellar populations in these 6 systems to: i) trace their different stellar generations over a spatial extension and with a depth that only the LSST can achieve; ii) measure their distances using variable stars of different type/parent stellar population and the Tip of the RGB; iii) map their 3D structures up to the periphery of their halos; iv) search for tidal streams; and v) study their Star Formation Histories over unprecedented large fractions of their bodies. Our ultimate goals are to provide a complete picture of these nearby stellar systems all the way through to their periphery, and to directly link and cross-calibrate the Gaia and LSST projects.

The Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC) is an open data challenge to classify simulated astronomical time-series data in preparation for observations from the Large Synoptic Survey Telescope (LSST), which will achieve first light in 2019 and commence its 10-year main survey in 2022. LSST will revolutionize our understanding of the changing sky, discovering and measuring millions of time-varying objects. In this challenge, we pose the question: how well can we classify objects in the sky that vary in brightness from simulated LSST time-series data, with all its challenges of non-representativity? In this note we explain the need for a data challenge to help classify such astronomical sources and describe the PLAsTiCC data set and Kaggle data challenge, noting that while the references are provided for context, they are not needed to participate in the challenge.

Classification of transient and variable light curves is an essential step in using astronomical observations to develop an understanding of their underlying physical processes. However, upcoming deep photometric surveys, including the Large Synoptic Survey Telescope (LSST), will produce a deluge of low signal-to-noise data for which traditional labeling procedures are inappropriate. Probabilistic classification is more appropriate for the data but are incompatible with the traditional metrics used on deterministic classifications. Furthermore, large survey collaborations intend to use these classification probabilities for diverse science objectives, indicating a need for a metric that balances a variety of goals. We describe the process used to develop an optimal performance metric for an open classification challenge that seeks probabilistic classifications and must serve many scientific interests. The Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC) is an open competition aiming to identify promising techniques for obtaining classification probabilities of transient and variable objects by engaging a broader community both within and outside astronomy. Using mock classification probability submissions emulating archetypes of those anticipated of PLAsTiCC, we compare the sensitivity of metrics of classification probabilities under various weighting schemes, finding that they yield qualitatively consistent results. We choose as a metric for PLAsTiCC a weighted modification of the cross-entropy because it can be meaningfully interpreted. Finally, we propose extensions of our methodology to ever more complex challenge goals and suggest some guiding principles for approaching the choice of a metric of probabilistic classifications.

In the analysis of VERITAS and other IACT data, it is expected to get a distribution of statistical significances in sky map bins with mean of zero and width of unity in the absence of a $\gamma$-ray signal. However, it is not uncommon to see significance distributions of width greater than unity, indicating that the background is poorly estimated and the significances in the region of interest are incorrect. This work explores the origins of these wider significance distributions and develop solutions to this issue and test these solutions on samples of VERITAS data.

In the cosmological paradigm, cold dark matter (DM) dominates the mass content of the Universe and is present at every scale. Candidates for DM include many extensions of the Standard Model with weakly interacting massive particles (WIMPs) in the mass range from ~10 GeV to greater than 10 TeV. The self-annihilation or decay of WIMPs in astrophysical regions of high DM density can produce secondary particles including very-high-energy (VHE; E>100 GeV) gamma rays. VERITAS, an array of atmospheric Cherenkov telescopes, sensitive to VHE gamma rays in the 85 GeV-30 TeV energy range, has been utilized for indirect DM searches. Astrophysical candidates for indirect DM detection by VERITAS are dwarf spheroidal galaxies (dSphs) of the Local Group and the Galactic Center, among others. This is a report on the current status of VERITAS observations of dark matter targets and the results of a combined search of four dSphs.

We report results of a wideband search for periodic gravitational waves from isolated neutron stars within the Orion spur towards both the inner and outer regions of our Galaxy. As gravitational waves interact very weakly with matter, the search is unimpeded by dust and concentrations of stars. One search disk (A) is $6.87^\circ$ in diameter and centered on $20^\textrm{h}10^\textrm{m}54.71^\textrm{s}+33^\circ33'25.29"$, and the other (B) is $7.45^\circ$ in diameter and centered on $8^\textrm{h}35^\textrm{m}20.61^\textrm{s}-46^\circ49'25.151"$. We explored the frequency range of 50-1500 Hz and frequency derivative from $0$ to $-5\times 10^{-9}$ Hz/s. A multi-stage, loosely coherent search program allowed probing more deeply than before in these two regions, while increasing coherence length with every stage. Rigorous followup parameters have winnowed initial coincidence set to only 70 candidates, to be examined manually. None of those 70 candidates proved to be consistent with an isolated gravitational wave emitter, and 95% confidence level upper limits were placed on continuous-wave strain amplitudes. Near $169$ Hz we achieve our lowest 95% CL upper limit on worst-case linearly polarized strain amplitude $h_0$ of $6.3\times 10^{-25}$, while at the high end of our frequency range we achieve a worst-case upper limit of $3.4\times 10^{-24}$ for all polarizations and sky locations.

The very-high-energy (VHE; E > 100 GeV) gamma-ray experiment, VERITAS, detected exceptionally bright flares from the high-mass X-ray binary LS I +61$^\circ$ 303 during the period October-December 2014. LS I +61$^\circ$ 303 is a known VHE gamma-ray source, the flux from which varies strongly with the orbital period of ~26.5 days. The maximum VHE flux is found around apastron (orbital phase ~0.6) at a level typically corresponding to 10-15% of the Crab Nebula flux (>300 GeV). During these most recent observations, relatively short (day scale), bright TeV flares were observed from the source around apastron in two orbital cycles (October and November). Both cases exhibited peak fluxes above 25% of the Crab Nebula flux (>300 GeV), making these the brightest VHE flares ever detected from this source. In the last orbital cycle observed (December), the source had returned to its historical level of activity. The results of these VERITAS observations from 2014 will be presented.

In 1989, the Whipple 10m telescope achieved the first indisputable detection of a TeV gamma-ray source, the Crab Nebula. Until its decommissioning in 2011, the Whipple telescope took regular measurements of the nebula. With the recent discovery of GeV gamma-ray flaring activity in the Crab Nebula, it is an opportune time to return to the Whipple telescope data set and search its extensive archive for evidence of TeV flares. A data set on the Crab Nebula spanning ten years, 2000 - 2010, is compiled and searched for day-scale flaring activity using a Bayesian-block binning algorithm. No evidence for significant flaring activity is found. Monte Carlo simulations show that low levels of flux increase on short timescales are difficult to detect. Assuming a flare duration of seven days, 99% confidence level upper limits are calculated for the possible frequency of five-fold, two-fold and 1.5-fold flares in the data set. An upper limit of 0.02 flares per year is found for the five-fold flare, and a limit of 0.27 flares per year is placed on the two-fold flare. The detection of the 1.5-fold flare is consistent with the false-positive rate of the method, and so cannot be excluded.

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

We present optical, X-ray, high energy ($\lessapprox 30$ GeV) and very high energy ($\gtrapprox 100$ GeV; VHE) observations of the high-frequency peaked blazar Mrk 421 taken between 2008 May 24 and June 23. A high energy $\gamma$-ray signal was detected by AGILE with \sqrt{TS}=4.5 on June 9--15, with $F(E>100 \mathrm{MeV})= 42^{+14}_{-12}\times 10^{-8}$ photons cm$^{-2}$ s$^{-1}$. This flaring state is brighter than the average flux observed by EGRET by a factor of $\sim$3, but still consistent with the highest EGRET flux. In hard X-rays (20-60 keV) SuperAGILE resolved a 5-day flare (June 9-15) peaking at $\sim$ 55 mCrab. SuperAGILE, RXTE/ASM and Swift/BAT data show a correlated flaring structure between soft and hard X-rays. Hints of the same flaring behavior are also detected in the simultaneous optical data provided by the GASP-WEBT. A Swift/XRT observation near the flaring maximum revealed the highest 2-10 keV flux ever observed from this source, of 2.6 $\times 10^{-9}$ erg cm$^{-2}$ s$^{-1}$ (i.e. > 100 mCrab). A peak synchrotron energy of $\sim$3 keV was derived, higher than typical values of $\sim$0.5-1 keV. VHE observations with MAGIC and VERITAS on June 6-8 show the flux peaking in a bright state, well correlated with the X-rays. This extraordinary set of simultaneous data, covering a twelve-decade spectral range, allowed for a deep analysis of the spectral energy distribution as well as of correlated light curves. The $\gamma$-ray flare can be interpreted within the framework of the synchrotron self-Compton model in terms of a rapid acceleration of leptons in the jet.

In May 2002, the BL Lacertae object (BL Lac), 1ES1959+650, was detected in the Very High Energy (VHE) regime by the VERITAS collaboration using the Whipple $\gamma$-ray observatory. During the following two months, the VERITAS collaboration observed episodes of flaring; at times, the nightly average VHE signal was three times as bright as the Crab Nebula, the standard candle in VHE astronomy. 1ES1959+650 was also monitored in the optical and X-ray bands. The X-ray light curve showed periodic behavior and correlation with the VHE and the optical bands.

We present the first strong detection of very high energy gamma-rays from the close (z=0.048) X-ray selected BL Lacertae object 1ES1959+650. Observations were made with the Whipple 10m telescope on Mt. Hopkins, Arizona, using the atmospheric Cherenkov imaging technique. The flux between May and July 2002 was highly variable, with a mean of 0.64 +/- 0.03 times the steady flux from the Crab Nebula and reaching a maximum of five Crab, with variability on timescales as short as seven hours.