-
Low-frequency radio continuum observations of edge-on galaxies are ideal to
study cosmic-ray electrons (CREs) in halos via radio synchrotron emission and
to measure magnetic field strengths. We obtained new observations of the
edge-on spiral galaxy NGC 891 at 129-163 MHz with the LOw Frequency ARray
(LOFAR) and at 13-18 GHz with the Arcminute Microkelvin Imager (AMI) and
combine them with recent high-resolution Very Large Array (VLA) observations at
1-2 GHz, enabling us to study the radio continuum emission over two orders of
magnitude in frequency. The spectrum of the integrated nonthermal flux density
can be fitted by a power law with a spectral steepening towards higher
frequencies or by a curved polynomial. Spectral flattening at low frequencies
due to free-free absorption is detected in star-forming regions of the disk.
The mean magnetic field strength in the halo is 7 +- 2 $\mu$G. The scale
heights of the nonthermal halo emission at 146 MHz are larger than those at 1.5
GHz everywhere, with a mean ratio of 1.7 +- 0.3, indicating that spectral
ageing of CREs is important and that diffusive propagation dominates. The halo
scale heights at 146 MHz decrease with increasing magnetic field strengths
which is a signature of dominating synchrotron losses of CREs. On the other
hand, the spectral index between 146 MHz and 1.5 GHz linearly steepens from the
disk to the halo, indicating that advection rather than diffusion is the
dominating CRE transport process. This issue calls for refined modelling of CRE
propagation.
-
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.
-
The Cherenkov Telescope Array, CTA, will be the major global observatory for
very high energy gamma-ray astronomy over the next decade and beyond. The
scientific potential of CTA is extremely broad: from understanding the role of
relativistic cosmic particles to the search for dark matter. CTA is an explorer
of the extreme universe, probing environments from the immediate neighbourhood
of black holes to cosmic voids on the largest scales. Covering a huge range in
photon energy from 20 GeV to 300 TeV, CTA will improve on all aspects of
performance with respect to current instruments.
The observatory will operate arrays on sites in both hemispheres to provide
full sky coverage and will hence maximize the potential for the rarest
phenomena such as very nearby supernovae, gamma-ray bursts or gravitational
wave transients. With 99 telescopes on the southern site and 19 telescopes on
the northern site, flexible operation will be possible, with sub-arrays
available for specific tasks. CTA will have important synergies with many of
the new generation of major astronomical and astroparticle observatories.
Multi-wavelength and multi-messenger approaches combining CTA data with those
from other instruments will lead to a deeper understanding of the broad-band
non-thermal properties of target sources.
The CTA Observatory will be operated as an open, proposal-driven observatory,
with all data available on a public archive after a pre-defined proprietary
period. Scientists from institutions worldwide have combined together to form
the CTA Consortium. This Consortium has prepared a proposal for a Core
Programme of highly motivated observations. The programme, encompassing
approximately 40% of the available observing time over the first ten years of
CTA operation, is made up of individual Key Science Projects (KSPs), which are
presented in this document.
-
We present radio continuum maps of 12 nearby ($D\leq 27~\rm Mpc$), edge-on
($i\geq 76^{\circ}$), late-type spiral galaxies mostly at $1.4$ and 5 GHz,
observed with the Australia Telescope Compact Array, Very Large Array,
Westerbork Synthesis Radio Telescope, Effelsberg 100-m and Parkes 64-m
telescopes. All galaxies show clear evidence of radio haloes, including the
first detection in the Magellanic-type galaxy NGC 55. In 11 galaxies, we find a
thin and a thick disc that can be better fitted by exponential rather than
Gaussian functions. We fit our SPINNAKER (SPectral INdex Numerical Analysis of
K(c)osmic-ray Electron Radio-emission) 1D cosmic-ray transport models to the
vertical model profiles of the non-thermal intensity and to the non-thermal
radio spectral index in the halo. We simultaneously fit for the advection speed
(or diffusion coefficient) and magnetic field scale height. In the thick disc,
the magnetic field scale heights range from 2 to 8 kpc with an average across
the sample of $3.0\pm 1.7~\rm kpc$; they show no correlation with either
star-formation rate (SFR), SFR surface density ($\Sigma_{\rm SFR}$) or rotation
speed ($V_{\rm rot}$). The advection speeds range from 100 to $700~\rm
km\,s^{-1}$ and display correlations of $V\propto \rm SFR^{0.36\pm 0.06}$ and
$V\propto \Sigma_{\rm SFR}^{0.39\pm 0.09}$; they agree remarkably well with the
escape velocities ($0.5\leq V/V_{\rm esc}\leq 2$), which can be explained by
cosmic-ray driven winds. Radio haloes show the presence of disc winds in
galaxies with $\Sigma_{\rm SFR} > 10^{-3}~\rm M_{\odot}\,yr^{-1}\,kpc^{-2}$
that extend over several kpc and are driven by processes related to the
distributed star formation in the disc.
-
The composition and morphology of the interstellar medium in starburst
galaxies has been well investigated, but the magnetic field properties are
still uncertain. The nearby starburst galaxy M82 provides a unique opportunity
to investigate the mechanisms leading to the amplification and reduction of
turbulent and regular magnetic fields. Possible scenarios of the contribution
of the magnetic field to the star-formation rate are evaluated. Archival data
from the VLA and WSRT were combined and re-reduced to cover the wavelength
regime between 3cm and 22cm. All observations revealed polarised emission in
the inner part of the galaxy, while extended polarised emission up to a
distance of 2kpc from the disk was only detected at 18cm and 22cm. The
observations hint at a magnetised bar in the inner part of the galaxy. We
calculate the mass inflow rate due to magnetic stress of the bar to 7.1 solar
masses per year, which can be a significant contribution to the star-formation
rate of M82 of approximately 13 solar masses per year. The halo shows polarised
emission, which might be the remnant of a regular disk field. Indications for a
helical field in the inner part of the outflow cone are provided. The coherence
length of the magnetic field in the centre is similar to the size of giant
molecular clouds. Using polarisation spectra more evidence for a close coupling
of the ionised gas and the magnetic field as well as a two-phase magnetic field
topology were found. Electron densities in the halo are similar to the ones
found in the Milky Way. The magnetic field morphology is similar to the one in
other nearby starburst galaxies with possible large-scale magnetic loops in the
halo and a helical magnetic field inside the outflow cones. The special
combination of a magnetic bar and a circumnuclear ring are able to
significantly raise the star-formation rate in this galaxy by magnetic braking.
-
List of contributions from the Cherenkov Telescope Array Consortium presented
at the 35th International Cosmic Ray Conference, July 12-20 2017, Busan, Korea.
-
Compilation of papers presented by the VERITAS Collaboration at the 35th
International Cosmic Ray Conference (ICRC), held July 12 through July 20, 2017
in Busan, South Korea.
-
Very-high-energy (VHE; $>$ 100 GeV) gamma-ray emission from the blazar RGB
J2243+203 was discovered with the VERITAS Cherenkov telescope array, during the
period between 21 and 24 December 2014. The VERITAS energy spectrum from this
source can be fit by a power law with a photon index of $4.6 \pm 0.5$, and a
flux normalization at 0.15 TeV of $(6.3 \pm 1.1) \times 10^{-10} ~
\textrm{cm}^{-2} \textrm{s}^{-1} \textrm{TeV}^{-1}$. The integrated
\textit{Fermi}-LAT flux from 1 GeV to 100 GeV during the VERITAS detection is
$(4.1 \pm 0.8) \times 10^{\textrm{-8}}
~\textrm{cm}^{\textrm{-2}}\textrm{s}^{\textrm{-1}}$, which is an order of
magnitude larger than the four-year-averaged flux in the same energy range
reported in the 3FGL catalog, ($4.0 \pm 0.1 \times 10^{\textrm{-9}} ~
\textrm{cm}^{\textrm{-2}}\textrm{s}^{\textrm{-1}}$). The detection with VERITAS
triggered observations in the X-ray band with the \textit{Swift}-XRT. However,
due to scheduling constraints \textit{Swift}-XRT observations were performed 67
hours after the VERITAS detection, not simultaneous with the VERITAS
observations. The observed X-ray energy spectrum between 2 keV and 10 keV can
be fitted with a power-law with a spectral index of $2.7 \pm 0.2$, and the
integrated photon flux in the same energy band is $(3.6 \pm 0.6) \times
10^{-13} ~\textrm{cm}^{-2} \textrm{s}^{-1}$. EBL model-dependent upper limits
of the blazar redshift have been derived. Depending on the EBL model used, the
upper limit varies in the range from z $<~0.9$ to z $<~1.1$.
-
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 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.
-
We present very-high-energy $\gamma$-ray observations of the BL Lac object
1ES 2344+514 taken by the Very Energetic Radiation Imaging Telescope Array
System (VERITAS) between 2007 and 2015. 1ES 2344+514 is detected with a
statistical significance above background of $20.8\sigma$ in $47.2$ hours
(livetime) of observations, making this the most comprehensive very-high-energy
study of 1ES 2344+514 to date. Using these observations the temporal properties
of 1ES 2344+514 are studied on short and long times scales. We fit a constant
flux model to nightly- and seasonally-binned light curves and apply a
fractional variability test, to determine the stability of the source on
different timescales. We reject the constant-flux model for the 2007-2008 and
2014-2015 nightly-binned light curves and for the long-term seasonally-binned
light curve at the $> 3\sigma$ level. The spectra of the time-averaged emission
before and after correction for attenuation by the extragalactic background
light are obtained. The observed time-averaged spectrum above 200 GeV is
satisfactorily fitted (${\chi^2/NDF = 7.89/6}$) by a power-law function with
index $\Gamma = 2.46 \pm 0.06_{stat} \pm 0.20_{sys} $ and extends to at least 8
TeV. The extragalactic-background-light-deabsorbed spectrum is adequately fit
(${\chi^2/NDF = 6.73/6}$) by a power-law function with index $\Gamma = 2.15 \pm
0.06_{stat} \pm 0.20_{sys} $ while an F-test indicates that the power-law with
exponential cutoff function provides a marginally-better fit ($\chi^2/NDF $ =
$2.56 / 5 $) at the 2.1$\sigma$ level. The source location is found to be
consistent with the published radio location and its spatial extent is
consistent with a point source.
-
The discovery of the Higgs boson by the LHC experiments ATLAS and CMS has
marked a milestone for particle physics. Yet, there are still many open
questions that cannot be answered within the Standard Model (SM). For example,
the generation of the observed matter-antimatter asymmetry in the universe
through baryogenesis can only be explained qualitatively in the SM. A simple
extension of the SM compatible with the current theoretical and experimental
constraints is given by the 2-Higgs-Doublet Model (2HDM) where a second Higgs
doublet is added to the Higgs sector. We investigate the possibility of a
strong first order electroweak phase transition in the CP-conserving 2HDM type
I and type II where either of the CP-even Higgs bosons is identified with the
SM-like Higgs boson. The renormalisation that we apply on the loop-corrected
Higgs potential allows us to efficiently scan the 2HDM parameter space and
simultaneously take into account all relevant theoretical and up-to-date
experimental constraints. The 2HDM parameter regions found to be compatible
with the applied constraints and a strong electroweak phase transition are
analysed systematically. Our results show that there is a strong interplay
between the requirement of a strong phase transition and collider phenomenology
with testable implications for searches at the LHC.
-
Imaging atmospheric Cherenkov telescopes (IACTs) are equipped with sensitive
photomultiplier tube (PMT) cameras. Exposure to high levels of background
illumination degrades the efficiency of and potentially destroys these
photo-detectors over time, so IACTs cannot be operated in the same
configuration in the presence of bright moonlight as under dark skies. Since
September 2012, observations have been carried out with the VERITAS IACTs under
bright moonlight (defined as about three times the night-sky-background (NSB)
of a dark extragalactic field, typically occurring when Moon illumination >
35%) in two observing modes, firstly by reducing the voltage applied to the
PMTs and, secondly, with the addition of ultra-violet (UV) bandpass filters to
the cameras. This has allowed observations at up to about 30 times previous NSB
levels (around 80% Moon illumination), resulting in 30% more observing time
between the two modes over the course of a year. These additional observations
have already allowed for the detection of a flare from the 1ES 1727+502 and for
an observing program targeting a measurement of the cosmic-ray positron
fraction. We provide details of these new observing modes and their performance
relative to the standard VERITAS observations.
-
High-energy gamma-ray emission from supernova remnants (SNRs) has provided a
unique perspective for studies of Galactic cosmic-ray acceleration. Tycho's SNR
is a particularly good target because it is a young, type Ia SNR that is
well-studied over a wide range of energies and located in a relatively clean
environment. Since the detection of gamma-ray emission from Tycho's SNR by
VERITAS and Fermi-LAT, there have been several theoretical models proposed to
explain its broadband emission and high-energy morphology. We report on an
update to the gamma-ray measurements of Tycho's SNR with 147 hours of VERITAS
and 84 months of Fermi-LAT observations, which represents about a factor of two
increase in exposure over previously published data. About half of the VERITAS
data benefited from a camera upgrade, which has made it possible to extend the
TeV measurements toward lower energies. The TeV spectral index measured by
VERITAS is consistent with previous results, but the expanded energy range
softens a straight power-law fit. At energies higher than 400 GeV, the
power-law index is $2.92 \pm 0.42_{\mathrm{stat}} \pm 0.20_{\mathrm{sys}}$. It
is also softer than the spectral index in the GeV energy range, $2.14 \pm
0.09_{\mathrm{stat}} \pm 0.02_{\mathrm{sys}}$, measured by this study using
Fermi--LAT data. The centroid position of the gamma-ray emission is coincident
with the center of the remnant, as well as with the centroid measurement of
Fermi--LAT above 1 GeV. The results are consistent with an SNR shell origin of
the emission, as many models assume. The updated spectrum points to a lower
maximum particle energy than has been suggested previously.
-
B2 1215+30 is a BL Lac-type blazar that was first detected at TeV energies by
the MAGIC atmospheric Cherenkov telescopes, and subsequently confirmed by the
VERITAS observatory with data collected between 2009 and 2012. In 2014 February
08, VERITAS detected a large-amplitude flare from B2 1215+30 during routine
monitoring observations of the blazar 1ES 1218+304, located in the same field
of view. The TeV flux reached 2.4 times the Crab Nebula flux with a variability
timescale of < 3.6 h. Multiwavelength observations with Fermi-LAT, Swift, and
the Tuorla observatory revealed a correlated high GeV flux state and no
significant optical counterpart to the flare, with a spectral energy
distribution where the gamma-ray luminosity exceeds the synchrotron luminosity.
When interpreted in the framework of a one-zone leptonic model, the observed
emission implies a high degree of beaming, with Doppler factor > 10, and an
electron population with spectral index < 2.3.
-
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 study the spectral energy distribution (SED) of the radio continuum
emission from the KINGFISH sample of nearby galaxies to understand the
energetics and origin of this emission. Effelsberg multi-wavelength
observations at 1.4GHz, 4.8GHz, 8.5GHz, and 10.5GHz combined with archive data
allow us, for the first time, to determine the mid-radio continuum (1-10 GHz,
MRC) bolometric luminosities and further present calibration relations vs. the
monochromatic radio luminosities. The 1-10 GHz radio SED is fitted using a
Bayesian Markov Chain Monte Carlo (MCMC) technique leading to measurements for
the nonthermal spectral index and the thermal fraction f_th with mean values of
alpha_nt=0.97+-0.16 (0.79+-0.15 for the total spectral index) and f_th= 10% +-
9% at 1.4 GHz. The MRC luminosity changes over ~3 orders of magnitude in the
sample. The thermal emission is responsible for ~23% of the MRC on average. We
also compare the extinction-corrected diagnostics of star formation rate with
the thermal and nonthermal radio tracers and derive the first star formation
calibration relations using the MRC radio luminosity. The nonthermal spectral
index flattens with increasing star formation rate surface density, indicating
the effect of the star formation feedback on the cosmic ray electron population
in galaxies. Comparing the radio and IR SEDs, we find that the FIR-to-MRC ratio
could decrease with star formation rate, due to the amplification of the
magnetic fields in star forming regions. This particularly implies a decrease
in the ratio at high redshifts, where mostly luminous/star forming galaxies are
detected.
-
Blazars are variable emitters across all wavelengths over a wide range of
timescales, from months down to minutes. It is therefore essential to observe
blazars simultaneously at different wavelengths, especially in the X-ray and
gamma-ray bands, where the broadband spectral energy distributions usually
peak.
In this work, we report on three "target-of-opportunity" (ToO) observations
of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring
event at TeV energies in 2014. These observations feature long, continuous, and
simultaneous exposures with XMM-Newton (covering X-ray and optical/ultraviolet
bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous
observations from other gamma-ray facilities (MAGIC and Fermi-LAT) and a number
of radio and optical facilities. Although neither rapid flares nor significant
X-ray/TeV correlation are detected, these observations reveal subtle changes in
the X-ray spectrum of the source over the course of a few days. We search the
simultaneous X-ray and TeV data for spectral hysteresis patterns and time
delays, which could provide insight into the emission mechanisms and the source
properties (e.g. the radius of the emitting region, the strength of the
magnetic field, and related timescales). The observed broadband spectra are
consistent with a one-zone synchrotron self-Compton model. We find that the
power spectral density distribution at $\gtrsim 4\times 10^{-4}$ Hz from the
X-ray data can be described by a power-law model with an index value between
1.2 and 1.8, and do not find evidence for a steepening of the power spectral
index (often associated with a characteristic length scale) compared to the
previously reported values at lower frequencies.
-
We describe and report the status of a neutrino-triggered program in IceCube
that generates real-time alerts for gamma-ray follow-up observations by
atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable
of monitoring the whole sky continuously, high-energy gamma-ray telescopes have
restricted fields of view and in general are unlikely to be observing a
potential neutrino-flaring source at the time such neutrinos are recorded. The
use of neutrino-triggered alerts thus aims at increasing the availability of
simultaneous multi-messenger data during potential neutrino flaring activity,
which can increase the discovery potential and constrain the phenomenological
interpretation of the high-energy emission of selected source classes (e.g.
blazars). The requirements of a fast and stable online analysis of potential
neutrino signals and its operation are presented, along with first results of
the program operating between 14 March 2012 and 31 December 2015.
-
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.
-
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%.
-
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 describe a method for spectral cleaning and timing calibration of short
voltage time series data from individual radio interferometer receivers. It
makes use of the phase differences in Fast Fourier Transform (FFT) spectra
across antenna pairs. For strong, localized terrestrial sources these are
stable over time, while being approximately uniform-random for a sum over many
sources or for noise. Using only milliseconds-long datasets, the method finds
the strongest interfering transmitters, a first-order solution for relative
timing calibrations, and faulty data channels. No knowledge of gain response or
quiescent noise levels of the receivers is required. With relatively small data
volumes, this approach is suitable for use in an online system monitoring setup
for interferometric arrays.
We have applied the method to our cosmic-ray data collection, a collection of
measurements of short pulses from extensive air showers, recorded by the LOFAR
radio telescope. Per air shower, we have collected 2 ms of raw time series data
for each receiver. The spectral cleaning has a calculated optimal sensitivity
corresponding to a power signal-to-noise ratio of 0.08 (or -11 dB) in a
spectral window of 25 kHz, for 2 ms of data in 48 antennas. This is well
sufficient for our application. Timing calibration across individual antenna
pairs has been performed at 0.4 ns precision; for calibration of signal clocks
across stations of 48 antennas the precision is 0.1 ns. Monitoring differences
in timing calibration per antenna pair over the course of the period 2011 to
2015 shows a precision of 0.08 ns, which is useful for monitoring and
correcting drifts in signal path synchronizations.
A cross-check method for timing calibration is presented, using a pulse
transmitter carried by a drone flying over the array. Timing precision is
similar, 0.3 ns.
-
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 Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black holes to cosmic voids on the largest scales. Covering a huge range in photon energy from 20 GeV to 300 TeV, CTA will improve on all aspects of performance with respect to current instruments. The observatory will operate arrays on sites in both hemispheres to provide full sky coverage and will hence maximize the potential for the rarest phenomena such as very nearby supernovae, gamma-ray bursts or gravitational wave transients. With 99 telescopes on the southern site and 19 telescopes on the northern site, flexible operation will be possible, with sub-arrays available for specific tasks. CTA will have important synergies with many of the new generation of major astronomical and astroparticle observatories. Multi-wavelength and multi-messenger approaches combining CTA data with those from other instruments will lead to a deeper understanding of the broad-band non-thermal properties of target sources. The CTA Observatory will be operated as an open, proposal-driven observatory, with all data available on a public archive after a pre-defined proprietary period. Scientists from institutions worldwide have combined together to form the CTA Consortium. This Consortium has prepared a proposal for a Core Programme of highly motivated observations. The programme, encompassing approximately 40% of the available observing time over the first ten years of CTA operation, is made up of individual Key Science Projects (KSPs), which are presented in this document.
We present radio continuum maps of 12 nearby ($D\leq 27~\rm Mpc$), edge-on ($i\geq 76^{\circ}$), late-type spiral galaxies mostly at $1.4$ and 5 GHz, observed with the Australia Telescope Compact Array, Very Large Array, Westerbork Synthesis Radio Telescope, Effelsberg 100-m and Parkes 64-m telescopes. All galaxies show clear evidence of radio haloes, including the first detection in the Magellanic-type galaxy NGC 55. In 11 galaxies, we find a thin and a thick disc that can be better fitted by exponential rather than Gaussian functions. We fit our SPINNAKER (SPectral INdex Numerical Analysis of K(c)osmic-ray Electron Radio-emission) 1D cosmic-ray transport models to the vertical model profiles of the non-thermal intensity and to the non-thermal radio spectral index in the halo. We simultaneously fit for the advection speed (or diffusion coefficient) and magnetic field scale height. In the thick disc, the magnetic field scale heights range from 2 to 8 kpc with an average across the sample of $3.0\pm 1.7~\rm kpc$; they show no correlation with either star-formation rate (SFR), SFR surface density ($\Sigma_{\rm SFR}$) or rotation speed ($V_{\rm rot}$). The advection speeds range from 100 to $700~\rm km\,s^{-1}$ and display correlations of $V\propto \rm SFR^{0.36\pm 0.06}$ and $V\propto \Sigma_{\rm SFR}^{0.39\pm 0.09}$; they agree remarkably well with the escape velocities ($0.5\leq V/V_{\rm esc}\leq 2$), which can be explained by cosmic-ray driven winds. Radio haloes show the presence of disc winds in galaxies with $\Sigma_{\rm SFR} > 10^{-3}~\rm M_{\odot}\,yr^{-1}\,kpc^{-2}$ that extend over several kpc and are driven by processes related to the distributed star formation in the disc.
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 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.
We present very-high-energy $\gamma$-ray observations of the BL Lac object 1ES 2344+514 taken by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) between 2007 and 2015. 1ES 2344+514 is detected with a statistical significance above background of $20.8\sigma$ in $47.2$ hours (livetime) of observations, making this the most comprehensive very-high-energy study of 1ES 2344+514 to date. Using these observations the temporal properties of 1ES 2344+514 are studied on short and long times scales. We fit a constant flux model to nightly- and seasonally-binned light curves and apply a fractional variability test, to determine the stability of the source on different timescales. We reject the constant-flux model for the 2007-2008 and 2014-2015 nightly-binned light curves and for the long-term seasonally-binned light curve at the $> 3\sigma$ level. The spectra of the time-averaged emission before and after correction for attenuation by the extragalactic background light are obtained. The observed time-averaged spectrum above 200 GeV is satisfactorily fitted (${\chi^2/NDF = 7.89/6}$) by a power-law function with index $\Gamma = 2.46 \pm 0.06_{stat} \pm 0.20_{sys} $ and extends to at least 8 TeV. The extragalactic-background-light-deabsorbed spectrum is adequately fit (${\chi^2/NDF = 6.73/6}$) by a power-law function with index $\Gamma = 2.15 \pm 0.06_{stat} \pm 0.20_{sys} $ while an F-test indicates that the power-law with exponential cutoff function provides a marginally-better fit ($\chi^2/NDF $ = $2.56 / 5 $) at the 2.1$\sigma$ level. The source location is found to be consistent with the published radio location and its spatial extent is consistent with a point source.
The discovery of the Higgs boson by the LHC experiments ATLAS and CMS has marked a milestone for particle physics. Yet, there are still many open questions that cannot be answered within the Standard Model (SM). For example, the generation of the observed matter-antimatter asymmetry in the universe through baryogenesis can only be explained qualitatively in the SM. A simple extension of the SM compatible with the current theoretical and experimental constraints is given by the 2-Higgs-Doublet Model (2HDM) where a second Higgs doublet is added to the Higgs sector. We investigate the possibility of a strong first order electroweak phase transition in the CP-conserving 2HDM type I and type II where either of the CP-even Higgs bosons is identified with the SM-like Higgs boson. The renormalisation that we apply on the loop-corrected Higgs potential allows us to efficiently scan the 2HDM parameter space and simultaneously take into account all relevant theoretical and up-to-date experimental constraints. The 2HDM parameter regions found to be compatible with the applied constraints and a strong electroweak phase transition are analysed systematically. Our results show that there is a strong interplay between the requirement of a strong phase transition and collider phenomenology with testable implications for searches at the LHC.
Imaging atmospheric Cherenkov telescopes (IACTs) are equipped with sensitive photomultiplier tube (PMT) cameras. Exposure to high levels of background illumination degrades the efficiency of and potentially destroys these photo-detectors over time, so IACTs cannot be operated in the same configuration in the presence of bright moonlight as under dark skies. Since September 2012, observations have been carried out with the VERITAS IACTs under bright moonlight (defined as about three times the night-sky-background (NSB) of a dark extragalactic field, typically occurring when Moon illumination > 35%) in two observing modes, firstly by reducing the voltage applied to the PMTs and, secondly, with the addition of ultra-violet (UV) bandpass filters to the cameras. This has allowed observations at up to about 30 times previous NSB levels (around 80% Moon illumination), resulting in 30% more observing time between the two modes over the course of a year. These additional observations have already allowed for the detection of a flare from the 1ES 1727+502 and for an observing program targeting a measurement of the cosmic-ray positron fraction. We provide details of these new observing modes and their performance relative to the standard VERITAS observations.
High-energy gamma-ray emission from supernova remnants (SNRs) has provided a unique perspective for studies of Galactic cosmic-ray acceleration. Tycho's SNR is a particularly good target because it is a young, type Ia SNR that is well-studied over a wide range of energies and located in a relatively clean environment. Since the detection of gamma-ray emission from Tycho's SNR by VERITAS and Fermi-LAT, there have been several theoretical models proposed to explain its broadband emission and high-energy morphology. We report on an update to the gamma-ray measurements of Tycho's SNR with 147 hours of VERITAS and 84 months of Fermi-LAT observations, which represents about a factor of two increase in exposure over previously published data. About half of the VERITAS data benefited from a camera upgrade, which has made it possible to extend the TeV measurements toward lower energies. The TeV spectral index measured by VERITAS is consistent with previous results, but the expanded energy range softens a straight power-law fit. At energies higher than 400 GeV, the power-law index is $2.92 \pm 0.42_{\mathrm{stat}} \pm 0.20_{\mathrm{sys}}$. It is also softer than the spectral index in the GeV energy range, $2.14 \pm 0.09_{\mathrm{stat}} \pm 0.02_{\mathrm{sys}}$, measured by this study using Fermi--LAT data. The centroid position of the gamma-ray emission is coincident with the center of the remnant, as well as with the centroid measurement of Fermi--LAT above 1 GeV. The results are consistent with an SNR shell origin of the emission, as many models assume. The updated spectrum points to a lower maximum particle energy than has been suggested previously.
B2 1215+30 is a BL Lac-type blazar that was first detected at TeV energies by the MAGIC atmospheric Cherenkov telescopes, and subsequently confirmed by the VERITAS observatory with data collected between 2009 and 2012. In 2014 February 08, VERITAS detected a large-amplitude flare from B2 1215+30 during routine monitoring observations of the blazar 1ES 1218+304, located in the same field of view. The TeV flux reached 2.4 times the Crab Nebula flux with a variability timescale of < 3.6 h. Multiwavelength observations with Fermi-LAT, Swift, and the Tuorla observatory revealed a correlated high GeV flux state and no significant optical counterpart to the flare, with a spectral energy distribution where the gamma-ray luminosity exceeds the synchrotron luminosity. When interpreted in the framework of a one-zone leptonic model, the observed emission implies a high degree of beaming, with Doppler factor > 10, and an electron population with spectral index < 2.3.
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 study the spectral energy distribution (SED) of the radio continuum emission from the KINGFISH sample of nearby galaxies to understand the energetics and origin of this emission. Effelsberg multi-wavelength observations at 1.4GHz, 4.8GHz, 8.5GHz, and 10.5GHz combined with archive data allow us, for the first time, to determine the mid-radio continuum (1-10 GHz, MRC) bolometric luminosities and further present calibration relations vs. the monochromatic radio luminosities. The 1-10 GHz radio SED is fitted using a Bayesian Markov Chain Monte Carlo (MCMC) technique leading to measurements for the nonthermal spectral index and the thermal fraction f_th with mean values of alpha_nt=0.97+-0.16 (0.79+-0.15 for the total spectral index) and f_th= 10% +- 9% at 1.4 GHz. The MRC luminosity changes over ~3 orders of magnitude in the sample. The thermal emission is responsible for ~23% of the MRC on average. We also compare the extinction-corrected diagnostics of star formation rate with the thermal and nonthermal radio tracers and derive the first star formation calibration relations using the MRC radio luminosity. The nonthermal spectral index flattens with increasing star formation rate surface density, indicating the effect of the star formation feedback on the cosmic ray electron population in galaxies. Comparing the radio and IR SEDs, we find that the FIR-to-MRC ratio could decrease with star formation rate, due to the amplification of the magnetic fields in star forming regions. This particularly implies a decrease in the ratio at high redshifts, where mostly luminous/star forming galaxies are detected.
Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring event at TeV energies in 2014. These observations feature long, continuous, and simultaneous exposures with XMM-Newton (covering X-ray and optical/ultraviolet bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-LAT) and a number of radio and optical facilities. Although neither rapid flares nor significant X-ray/TeV correlation are detected, these observations reveal subtle changes in the X-ray spectrum of the source over the course of a few days. We search the simultaneous X-ray and TeV data for spectral hysteresis patterns and time delays, which could provide insight into the emission mechanisms and the source properties (e.g. the radius of the emitting region, the strength of the magnetic field, and related timescales). The observed broadband spectra are consistent with a one-zone synchrotron self-Compton model. We find that the power spectral density distribution at $\gtrsim 4\times 10^{-4}$ Hz from the X-ray data can be described by a power-law model with an index value between 1.2 and 1.8, and do not find evidence for a steepening of the power spectral index (often associated with a characteristic length scale) compared to the previously reported values at lower frequencies.
We describe and report the status of a neutrino-triggered program in IceCube that generates real-time alerts for gamma-ray follow-up observations by atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable of monitoring the whole sky continuously, high-energy gamma-ray telescopes have restricted fields of view and in general are unlikely to be observing a potential neutrino-flaring source at the time such neutrinos are recorded. The use of neutrino-triggered alerts thus aims at increasing the availability of simultaneous multi-messenger data during potential neutrino flaring activity, which can increase the discovery potential and constrain the phenomenological interpretation of the high-energy emission of selected source classes (e.g. blazars). The requirements of a fast and stable online analysis of potential neutrino signals and its operation are presented, along with first results of the program operating between 14 March 2012 and 31 December 2015.
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.
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%.
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 describe a method for spectral cleaning and timing calibration of short voltage time series data from individual radio interferometer receivers. It makes use of the phase differences in Fast Fourier Transform (FFT) spectra across antenna pairs. For strong, localized terrestrial sources these are stable over time, while being approximately uniform-random for a sum over many sources or for noise. Using only milliseconds-long datasets, the method finds the strongest interfering transmitters, a first-order solution for relative timing calibrations, and faulty data channels. No knowledge of gain response or quiescent noise levels of the receivers is required. With relatively small data volumes, this approach is suitable for use in an online system monitoring setup for interferometric arrays. We have applied the method to our cosmic-ray data collection, a collection of measurements of short pulses from extensive air showers, recorded by the LOFAR radio telescope. Per air shower, we have collected 2 ms of raw time series data for each receiver. The spectral cleaning has a calculated optimal sensitivity corresponding to a power signal-to-noise ratio of 0.08 (or -11 dB) in a spectral window of 25 kHz, for 2 ms of data in 48 antennas. This is well sufficient for our application. Timing calibration across individual antenna pairs has been performed at 0.4 ns precision; for calibration of signal clocks across stations of 48 antennas the precision is 0.1 ns. Monitoring differences in timing calibration per antenna pair over the course of the period 2011 to 2015 shows a precision of 0.08 ns, which is useful for monitoring and correcting drifts in signal path synchronizations. A cross-check method for timing calibration is presented, using a pulse transmitter carried by a drone flying over the array. Timing precision is similar, 0.3 ns.
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.