-
List of contributions from the Cherenkov Telescope Array Consortium presented
at the 35th International Cosmic Ray Conference, July 12-20 2017, Busan, Korea.
-
List of contributions from the Cherenkov Telescope Array (CTA) Consortium
presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy
(Gamma 2016), July 11-15, 2016, in Heidelberg, Germany.
-
The binary millisecond radio pulsar PSR J1023+0038 exhibits many
characteristics similar to the gamma-ray binary system PSR B1259--63/LS 2883,
making it an ideal candidate for the study of high-energy non-thermal emission.
It has been the subject of multi-wavelength campaigns following the
disappearance of the pulsed radio emission in 2013 June, which revealed the
appearance of an accretion disk around the neutron star. We present the results
of very high-energy gamma-ray observations carried out by VERITAS before and
after this change of state. Searches for steady and pulsed emission of both
data sets yield no significant gamma-ray signal above 100 GeV, and upper limits
are given for both a steady and pulsed gamma-ray flux. These upper limits are
used to constrain the magnetic field strength in the shock region of the PSR
J1023+0038 system. Assuming that very high-energy gamma rays are produced via
an inverse-Compton mechanism in the shock region, we constrain the shock
magnetic field to be greater than $\sim$2 G before the disappearance of the
radio pulsar and greater than $\sim$10 G afterwards.
-
The BL Lacertae object 1ES 1440+122 was observed in the energy range from 85
GeV to 30 TeV by the VERITAS array of imaging atmospheric Cherenkov telescopes.
The observations, taken between 2008 May and 2010 June and totalling 53 hours,
resulted in the discovery of $\gamma$-ray emission from the blazar, which has a
redshift $z$=0.163. 1ES 1440+122 is detected at a statistical significance of
5.5 standard deviations above the background with an integral flux of
(2.8$\pm0.7_{\mathrm{stat}}\pm0.8_{\mathrm{sys}}$) $\times$ 10$^{-12}$
cm$^{-2}$ s$^{-1}$ (1.2\% of the Crab Nebula's flux) above 200 GeV. The
measured spectrum is described well by a power law from 0.2 TeV to 1.3 TeV with
a photon index of 3.1 $\pm$ 0.4$_{\mathrm{stat}}$ $\pm$ 0.2$_{\mathrm{sys}}$.
Quasi-simultaneous multi-wavelength data from the Fermi Large Area Telescope
(0.3--300 GeV) and the Swift X-ray Telescope (0.2--10 keV) are additionally
used to model the properties of the emission region. A synchrotron self-Compton
model produces a good representation of the multi-wavelength data. Adding an
external-Compton or a hadronic component also adequately describes the data.
-
The very high energy (VHE; E $>$ 100 GeV) blazar Markarian 501 was observed
between April 17 and May 5 (MJD 54938--54956), 2009, as part of an extensive
multi-wavelength campaign from radio to VHE. Strong VHE $\gamma$-ray activity
was detected on May 1st with Whipple and VERITAS, when the flux (E $>$ 400 GeV)
increased to 10 times the pre-flare baseline flux ($3.9{\times 10^{-11}}~{\rm
ph~cm^{-2}~s^{-1}}$), reaching five times the flux of the Crab Nebula. This
coincided with a decrease in the optical polarization and a rotation of the
polarization angle by 15$^{\circ}$. This VHE flare showed a fast flux variation
with an increase of a factor $\sim$4 in 25 minutes, and a falling time of
$\sim$50 minutes. We present the observations of the quiescent state previous
to the flare and of the high state after the flare, focusing on the flux and
spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined
with optical and radio data.
-
We present results from multiwavelength observations of the BL Lacertae
object 1ES 1741+196, including results in the very-high-energy $\gamma$-ray
regime using the Very Energetic Radiation Imaging Telescope Array System
(VERITAS). The VERITAS time-averaged spectrum, measured above 180 GeV, is
well-modelled by a power law with a spectral index of
$2.7\pm0.7_{\mathrm{stat}}\pm0.2_{\mathrm{syst}}$. The integral flux above 180
GeV is $(3.9\pm0.8_{\mathrm{stat}}\pm1.0_{\mathrm{syst}})\times 10^{-8}$
m$^{-2}$ s$^{-1}$, corresponding to 1.6% of the Crab Nebula flux on average.
The multiwavelength spectral energy distribution of the source suggests that
1ES 1741+196 is an extreme-high-frequency-peaked BL Lacertae object. The
observations analysed in this paper extend over a period of six years, during
which time no strong flares were observed in any band. This analysis is
therefore one of the few characterizations of a blazar in a non-flaring state.
-
Between the beginning of its full-scale scientific operations in 2007 and
2012, the VERITAS Cherenkov telescope array observed more than 130 blazars; of
these, 26 were detected as very-high-energy (VHE; E>100 GeV) {\gamma}-ray
sources. In this work, we present the analysis results of a sample of 114
undetected objects. The observations constitute a total live-time of ~570
hours. The sample includes several unidentified Fermi-Large Area Telescope
(LAT) sources (located at high Galactic latitude) as well as all the sources
from the second Fermi-LAT catalog which are contained within the field of view
of the VERITAS observations. We have also performed optical spectroscopy
measurements in order to estimate the redshift of some of these blazars that do
not have a spectroscopic distance estimate. We present new optical spectra from
the Kast instrument on the Shane telescope at the Lick observatory for 18
blazars included in this work, which allowed for the successful measurement or
constraint on the redshift of four of them. For each of the blazars included in
our sample we provide the flux upper limit in the VERITAS energy band. We also
study the properties of the significance distributions and we present the
result of a stacked analysis of the data-set, which shows a 4 {\sigma} excess.
-
The Galactic Center Ridge has been observed extensively in the past by both
GeV and TeV gamma-ray instruments revealing a wealth of structure, including a
diffuse component as well as the point sources G0.9+0.1 (a composite supernova
remnant) and Sgr A* (believed to be associated with the supermassive black hole
located at the center of our Galaxy). Previous very high energy (VHE) gamma-ray
observations with the H.E.S.S. experiment have also detected an extended TeV
gamma-ray component along the Galactic plane in the >300 GeV gamma-ray regime.
Here we report on observations of the Galactic Center Ridge from 2010-2014 by
the VERITAS telescope array in the >2 TeV energy range. From these observations
we 1.) provide improved measurements of the differential energy spectrum for
Sgr A* in the >2 TeV gamma-ray regime, 2.) provide a detection in the >2 TeV
gamma-ray emission from the composite SNR G0.9+0.1 and an improved
determination of its multi-TeV gamma-ray energy spectrum, 3.) report on the
detection of VER J1746-289, a localized enhancement of >2 TeV gamma-ray
emission along the Galactic plane.
-
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.
-
The TeV binary system LS I +61$^\circ$ 303 is known for its regular,
non-thermal emission pattern which traces the orbital period of the compact
object in its 26.5 day orbit around its B0 Ve star companion. The system
typically presents elevated TeV emission around apastron passage with flux
levels between 5% and 15% of the steady flux from the Crab Nebula (> 300 GeV).
In this article, VERITAS observations of LS I +61$^\circ$ 303 taken in late
2014 are presented, during which bright TeV flares around apastron at flux
levels peaking above 30% of the Crab Nebula flux were detected. This is the
brightest such activity from this source ever seen in the TeV regime. The
strong outbursts have rise and fall times of less than a day. The short
timescale of the flares, in conjunction with the observation of 10 TeV photons
from LS I +61$^\circ$ 303 during the flares, provides constraints on the
properties of the accelerator in the source.
-
Outbursts from gamma-ray quasars provide insights on the relativistic jets of
active galactic nuclei and constraints on the diffuse radiation fields that
fill the Universe. The detection of significant emission above 100 GeV from a
distant quasar would show that some of the radiated gamma rays escape
pair-production interactions with low-energy photons, be it the extragalactic
background light (EBL), or the radiation near the supermassive black hole lying
at the jet's base. VERITAS detected gamma-ray emission up to 200 GeV from PKS
1441+25 (z=0.939) during April 2015, a period of high activity across all
wavelengths. This observation of PKS 1441+25 suggests that the emission region
is located thousands of Schwarzschild radii away from the black hole. The
gamma-ray detection also sets a stringent upper limit on the near-ultraviolet
to near-infrared EBL intensity, suggesting that galaxy surveys have resolved
most, if not all, of the sources of the EBL at these wavelengths.
-
We present coordinated multiwavelength observations of the bright, nearby BL
Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift,
NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments,
providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR
yielded previously unattainable sensitivity in the 3-79 keV range, revealing
that the spectrum softens when the source is dimmer until the X-ray spectral
shape saturates into a steep power law with a photon index of approximately 3,
with no evidence for an exponential cutoff or additional hard components up to
about 80 keV. For the first time, we observed both the synchrotron and the
inverse-Compton peaks of the spectral energy distribution (SED) simultaneously
shifted to frequencies below the typical quiescent state by an order of
magnitude. The fractional variability as a function of photon energy shows a
double-bump structure which relates to the two bumps of the broadband SED. In
each bump, the variability increases with energy which, in the framework of the
synchrotron self-Compton model, implies that the electrons with higher energies
are more variable. The measured multi-band variability, the significant
X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in
both bands, the lack of correlation between optical/UV and X-ray flux, the low
degree of polarization and its significant (random) variations, the short
estimated electron cooling time, and the significantly longer variability
timescale observed in the NuSTAR light curves point toward in-situ electron
acceleration, and suggest that there are multiple compact regions contributing
to the broadband emission of Mrk 421 during low-activity states.
-
Compilation of papers presented by the VERITAS Collaboration at the 34th
International Cosmic Ray Conference (ICRC), held July 30 through August 6, 2015
in The Hague, The Netherlands.
-
The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is a
ground-based array located at the Fred Lawrence Whipple Observatory in southern
Arizona and is one of the world's most sensitive gamma-ray instruments at
energies of 85 GeV to $>$30 TeV. VERITAS has a wide scientific reach that
includes the study of extragalactic and Galactic objects as well as the search
for astrophysical signatures of dark matter and the measurement of cosmic rays.
In this paper, we will summarize the current status of the VERITAS observatory
and present some of the scientific highlights from the last two years, focusing
in particular on those results shown at the 2015 ICRC in The Hague,
Netherlands.
-
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.
-
List of contributions from the CTA Consortium presented at the 34th
International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The
Netherlands.
-
The VERITAS array of Cherenkov telescopes has carried out a deep
observational program on the nearby dwarf spheroidal galaxy Segue 1. We report
on the results of nearly 48 hours of good quality selected data, taken between
January 2010 and May 2011. No significant $\gamma$-ray emission is detected at
the nominal position of Segue 1, and upper limits on the integrated flux are
derived. According to recent studies, Segue 1 is the most dark matter-dominated
dwarf spheroidal galaxy currently known. We derive stringent bounds on various
annihilating and decaying dark matter particle models. The upper limits on the
velocity-weighted annihilation cross-section are $\mathrm{<\sigma v >^{95% CL}
\lesssim 10^{-23} cm^{3} s^{-1}}$, improving our limits from previous
observations of dwarf spheroidal galaxies by at least a factor of two for dark
matter particle masses $\mathrm{m_{\chi}\gtrsim 300 GeV}$. The lower limits on
the decay lifetime are at the level of $\mathrm{\tau^{95% CL} \gtrsim 10^{24}
s}$. Finally, we address the interpretation of the cosmic ray lepton anomalies
measured by ATIC and PAMELA in terms of dark matter annihilation, and show that
the VERITAS observations of Segue 1 disfavor such a scenario.
-
During moonlit nights, observations with ground-based Cherenkov telescopes at
very high energies (VHE, $E>100$ GeV) are constrained since the photomultiplier
tubes (PMTs) in the telescope camera are extremely sensitive to the background
moonlight. Observations with the VERITAS telescopes in the standard
configuration are performed only with a moon illumination less than 35$\%$ of
full moon. Since 2012, the VERITAS collaboration has implemented a new
observing mode under bright moonlight, by either reducing the voltage applied
to the PMTs (reduced-high-voltage configuration, RHV), or by utilizing
UV-transparent filters. While these operating modes result in lower sensitivity
and increased energy thresholds, the extension of the available observing time
is useful for monitoring variable sources such as blazars and sources requiring
spectral measurements at the highest energies. In this paper we report the
detection of $\gamma$-ray flaring activity from the BL Lac object 1ES 1727+502
during RHV observations. This detection represents the first evidence of VHE
variability from this blazar. The integral flux is
$(1.1\pm0.2)\times10^{-11}\mathrm{cm^{-2}s^{-1}}$ above 250 GeV, which is about
five times higher than the low-flux state. The detection triggered additional
\veritas\ observations during standard dark-time. Multiwavelength observations
with the FLWO 48" telescope, and the Swift and Fermi satellites are presented
and used to produce the first spectral energy distribution (SED) of this object
during $\gamma$-ray flaring activity. The SED is then fitted with a standard
synchrotron-self-Compton model, placing constraints on the properties of the
emitting region and of the acceleration mechanism at the origin of the
relativistic particle population in the jet.
-
We summarize broadband observations of the TeV-emitting blazar 1ES 1959+650,
including optical R-band observations by the robotic telescopes Super-LOTIS and
iTelescope, UV observations by Swift UVOT, X-ray observations by the Swift
X-ray Telescope (XRT), high-energy gamma-ray observations with the Fermi Large
Area Telescope (LAT) and very-high-energy (VHE) gamma-ray observations by
VERITAS above 315 GeV, all taken between 17 April 2012 and 1 June 2012 (MJD
56034 and 56079). The contemporaneous variability of the broadband spectral
energy distribution is explored in the context of a simple synchrotron self
Compton (SSC) model. In the SSC emission scenario, we find that the parameters
required to represent the high state are significantly different than those in
the low state. Motivated by possible evidence of gas in the vicinity of the
blazar, we also investigate a reflected-emission model to describe the observed
variability pattern. This model assumes that the non-thermal emission from the
jet is reflected by a nearby cloud of gas, allowing the reflected emission to
re-enter the blob and produce an elevated gamma-ray state with no simultaneous
elevated synchrotron flux. The model applied here, although not required to
explain the observed variability pattern, represents one possible scenario
which can describe the observations. As applied to an elevated VHE state of 66%
of the Crab Nebula flux, observed on a single night during the observation
period, the reflected-emission scenario does not support a purely leptonic
non-thermal emission mechanism. The reflected emission model does, however,
predict a reflected photon field with sufficient energy to enable elevated
gamma-ray emission via pion production with protons of energies between 10 and
100 TeV.
-
We present results from VERITAS observations of the BL Lac object PG 1553+113
spanning the years 2010, 2011, and 2012. The time-averaged spectrum, measured
between 160 and 560\,GeV, is well described by a power law with a spectral
index of $4.33 \pm 0.09$. The time-averaged integral flux above $200\,$GeV
measured for this period was $(1.69 \pm 0.06) \times 10^{-11} \, \mathrm{ph} \,
\mathrm{cm}^{-2} \, \mathrm{s}^{-1}$, corresponding to 6.9\% of the Crab Nebula
flux. We also present the combined $\gamma$-ray spectrum from the Fermi Large
Area Telescope and VERITAS covering an energy range from 100~MeV to 560~GeV.
The data are well fit by a power law with an exponential cutoff at $\rm {101.9
\pm 3.2 \, \mathrm{GeV}} $. The origin of the cutoff could be intrinsic to
PG~1553+113 or be due to the $\gamma$-ray opacity of our universe through pair
production off the extragalactic background light (EBL). Given lower limits to
the redshift of $\rm z \negthinspace > \negthinspace 0.395$ based on optical/UV
observations of PG~1553+113, the cutoff would be dominated by EBL absorption.
Conversely, the small statistical uncertainties of the VERITAS energy spectrum
have allowed us to provide a robust upper limit on the redshift of PG 1553+113
of $z \negthinspace \leq \negthinspace 0.62$. A strongly-elevated mean flux of
$(2.50 \pm 0.14) \times 10^{-11} \, \mathrm{ph} \, \mathrm{cm}^{-2} \,
\mathrm{s}^{-1}$ (10.3\% of the Crab Nebula flux) was observed during 2012,
with the daily flux reaching as high as $(4.44 \pm 0.71) \times 10^{-11} \,
\mathrm{ph} \, \mathrm{cm}^{-2} \, \mathrm{s}^{-1}$ (18.3\% of the Crab Nebula
flux) on MJD 56048. The light curve measured during the 2012 observing season
is marginally inconsistent with a steady flux, giving a $\chi^2$ probability
for a steady flux of 0.03\%.
-
We report on multifrequency observations performed during 2012 December-2013
August of the first narrow-line Seyfert 1 galaxy detected in gamma rays, PMN
J0948+0022 ($z$ = 0.5846). A gamma-ray flare was observed by the Large Area
Telescope on board Fermi during 2012 December-2013 January, reaching a daily
peak flux in the 0.1-100 GeV energy range of (155 $\pm$ 31) $\times$10$^{-8}$
ph cm$^{-2}$ s$^{-1}$ on 2013 January 1, corresponding to an apparent isotropic
luminosity of about 1.5$\times$10$^{48}$ erg s$^{-1}$. The gamma-ray flaring
period triggered Swift and VERITAS observations in addition to radio and
optical monitoring by OVRO, MOJAVE, and CRTS. A strong flare was observed in
optical, UV, and X-rays on 2012 December 30, quasi-simultaneously to the
gamma-ray flare, reaching a record flux for this source from optical to gamma
rays. VERITAS observations at very high energy (E > 100 GeV) during 2013
January 6-17 resulted in an upper limit of F (> 0.2 TeV) <
4.0$\times$10$^{-12}$ ph cm$^{-2}$ s$^{-1}$. We compared the spectral energy
distribution (SED) of the flaring state in 2013 January with that of an
intermediate state observed in 2011. The two SEDs, modelled as synchrotron
emission and an external Compton scattering of seed photons from a dust torus,
can be modelled by changing both the electron distribution parameters and the
magnetic field.
-
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.
-
Prompt emission from the very fluent and nearby (z=0.34) gamma-ray burst GRB
130427A was detected by several orbiting telescopes and by ground-based,
wide-field-of-view optical transient monitors. Apart from the intensity and
proximity of this GRB, it is exceptional due to the extremely long-lived
high-energy (100 MeV to 100 GeV) gamma-ray emission, which was detected by the
Large Area Telescope on the Fermi Gamma-ray Space Telescope for ~70 ks after
the initial burst. The persistent, hard-spectrum, high-energy emission suggests
that the highest-energy gamma rays may have been produced via synchrotron
self-Compton processes though there is also evidence that the high-energy
emission may instead be an extension of the synchrotron spectrum. VERITAS, a
ground-based imaging atmospheric Cherenkov telescope array, began follow-up
observations of GRB 130427A ~71 ks (~20 hr) after the onset of the burst. The
GRB was not detected with VERITAS; however, the high elevation of the
observations, coupled with the low redshift of the GRB, make VERITAS a very
sensitive probe of the emission from GRB 130427A for E > 100 GeV. The
non-detection and consequent upper limit derived place constraints on the
synchrotron self-Compton model of high-energy gamma-ray emission from this
burst.
-
We report on deep observations of the extended TeV gamma-ray source MGRO
J1908+06 made with the VERITAS very high energy (VHE) gamma-ray observatory.
Previously, the TeV emission has been attributed to the pulsar wind nebula
(PWN) of the Fermi-LAT pulsar PSR J1907+0602. We detect MGRO J1908+06 at a
significance level of 14 standard deviations (14 sigma) and measure a photon
index of 2.20 +/- 0.10_stat +/- 0.20_sys. The TeV emission is extended,
covering the region near PSR J1907+0602 and also extending towards SNR
G40.5--0.5. When fitted with a 2-dimensional Gaussian, the intrinsic extension
has a standard deviation of sigma_src = 0.44 +/- 0.02 degrees. In contrast to
other TeV PWNe of similar age in which the TeV spectrum softens with distance
from the pulsar, the TeV spectrum measured near the pulsar location is
consistent with that measured at a position near the rim of G40.5--0.5, 0.33
degrees away.
-
We present deep VERITAS observations of the blazar PKS 1424+240, along with
contemporaneous Fermi Large Area Telescope, Swift X-ray Telescope and Swift UV
Optical Telescope data between 2009 February 19 and 2013 June 8. This blazar
resides at a redshift of $z\ge0.6035$, displaying a significantly attenuated
gamma-ray flux above 100 GeV due to photon absorption via pair-production with
the extragalactic background light. We present more than 100 hours of VERITAS
observations from three years, a multiwavelength light curve and the
contemporaneous spectral energy distributions. The source shows a higher flux
of (2.1$\pm0.3$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV in 2009 and
2011 as compared to the flux measured in 2013, corresponding to
(1.02$\pm0.08$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV. The measured
differential very high energy (VHE; $E\ge100$ GeV) spectral indices are
$\Gamma=$3.8$\pm$0.3, 4.3$\pm$0.6 and 4.5$\pm$0.2 in 2009, 2011 and 2013,
respectively. No significant spectral change across the observation epochs is
detected. We find no evidence for variability at gamma-ray opacities of greater
than $\tau=2$, where it is postulated that any variability would be small and
occur on longer than year timescales if hadronic cosmic-ray interactions with
extragalactic photon fields provide a secondary VHE photon flux. The data
cannot rule out such variability due to low statistics.