sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Orr, Matthew

Normalized to: Orr, M.

46 article(s) in total. 1513 co-authors, from 1 to 28 common article(s). Median position in authors list is 51,0.

[1]  oai:arXiv.org:2005.12916  [pdf] - 2102972
Pressure balance in the multiphase ISM of cosmologically simulated disk galaxies
Comments: 20 pages, 14 figures, submitted to MNRAS
Submitted: 2020-05-26
Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disk galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyze how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the midplane. Bulk flows (e.g., inflows and fountains) are important at a few disk scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total midplane pressure is well-predicted by the weight of the disk gas, and we show that it also scales linearly with the star formation rate surface density (Sigma_SFR). These results support the notion that the Kennicutt-Schmidt relation arises because Sigma_SFR and the gas surface density (Sigma_g) are connected via the ISM midplane pressure.
[2]  oai:arXiv.org:1912.04290  [pdf] - 2043561
Variations in the slope of the resolved star-forming main sequence: a tool for constraining the mass of star-forming regions
Comments: Accepted for publication in MNRAS Letters
Submitted: 2019-12-09, last modified: 2020-02-04
The correlation between galaxies' integrated stellar masses and star formation rates (the `star formation main sequence'; SFMS) is a well-established scaling relation. Recently, surveys have found a relationship between the star formation rate and stellar mass surface densities on kpc and sub-kpc scales (the `resolved SFMS'; rSFMS). In this work, we demonstrate that the rSFMS emerges naturally in FIRE-2 zoom-in simulations of Milky Way-mass galaxies. We make SFR and stellar mass maps of the simulated galaxies at a variety of spatial resolutions and star formation averaging time-scales and fit the rSFMS using multiple methods from the literature. While the absolute value of the SFMS slope depends on the fitting method, the slope is steeper for longer star formation time-scales and lower spatial resolutions regardless of the fitting method employed. We present a toy model that quantitatively captures the dependence of the simulated galaxies' rSFMS slope on spatial resolution and use it to illustrate how this dependence can be used to constrain the characteristic mass of star-forming clumps.
[3]  oai:arXiv.org:2001.08197  [pdf] - 2035880
Reproducing the CO-to-H$_2$ conversion factor in cosmological simulations of Milky Way-mass galaxies
Comments: Submitted to MNRAS. 14 pages, 10 figures
Submitted: 2020-01-22
We present models of CO(1-0) emission from Milky Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which in our models sets the attenuation of the incident UV radiation field. Commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H$_2$ abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of $\sim 3$ pc in cold gas ($T < 300$ K) for both CO and H$_{2}$, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H$_{2}$ conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model which controls the amount of radiation that penetrates giant molecular clouds.
[4]  oai:arXiv.org:1911.05251  [pdf] - 1998133
Live Fast, Die Young: GMC lifetimes in the FIRE cosmological simulations of Milky Way-mass galaxies
Comments: submitted to MNRAS
Submitted: 2019-11-12
We present the first measurement of the lifetimes of Giant Molecular Clouds (GMCs) in cosmological simulations at $z = 0$, using the Latte suite of FIRE-2 simulations of Milky Way-mass galaxies. We track GMCs with total gas mass $\gtrsim 10^5$ M$_\odot$ at high spatial ($\sim1$ pc), mass ($7100$ M$_{\odot}$), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the Milky Way and nearby galaxies. We find GMC lifetimes of $5-7$ Myr, or 1-2 freefall times, on average, with less than 1\% of clouds living longer than 20 Myr. We find increasing GMC lifetimes with galactocentric radius, implying that environment affects the evolutionary cycle of GMCs. However, our GMC lifetimes show no systematic dependence on GMC mass or amount of star formation. These results are broadly consistent with inferences from the literature and provide an initial investigation into ultimately understanding the physical processes that govern GMC lifetimes in a cosmological setting.
[5]  oai:arXiv.org:1911.00020  [pdf] - 2115122
Swirls of FIRE: Spatially Resolved Gas Velocity Dispersions and Star Formation Rates in FIRE-2 Disk Environments
Comments: 18 pages, 12 figures, submitted to MNRAS
Submitted: 2019-10-31
We study the spatially resolved (sub-kpc) gas velocity dispersion ($\sigma$)--star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way mass disk galaxies at late times. In agreement with observations, we find a relatively flat relationship, with $\sigma \approx 15-30$ km/s in neutral gas across 3 dex in SFRs. We show that higher dense gas fractions (ratios of dense gas to neutral gas) and SFRs are correlated at constant $\sigma$. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher $\sigma$ at constant SFR. The limits of the $\sigma$-$\Sigma_{\rm SFR}$ relation correspond to the onset of strong outflows. We see evidence of "on-off" cycles of star formation in the simulations, corresponding to feedback injection timescales of 10-100 Myr, where SFRs oscillate about equilibrium SFR predictions. Finally, SFRs and velocity dispersions in the simulations agree well with feedback-regulated and marginally stable gas disk (Toomre's $Q =1$) model predictions, and the data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e., ${\rm 1\%}$ per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of $\sigma$, it appears to be strongly subdominant to stellar feedback in the simulated galaxy disks.
[6]  oai:arXiv.org:1810.09460  [pdf] - 1876406
A Simple Non-equilibrium Feedback Model for Galaxy-Scale Star Formation: Delayed Feedback and SFR Scatter
Comments: 15 pages, 6 figures, accepted by MNRAS (04/25/2019)
Submitted: 2018-10-22, last modified: 2019-04-30
We explore a class of simple non-equilibrium star formation models within the framework of a feedback-regulated model of the ISM, applicable to kiloparsec-scale resolved star formation relations (e.g. Kennicutt-Schmidt). Combining a Toomre-Q-dependent local star formation efficiency per free-fall time with a model for delayed feedback, we are able to match the normalization and scatter of resolved star formation scaling relations. In particular, this simple model suggests that large ($\sim$dex) variations in star formation rates (SFRs) on kiloparsec scales may be due to the fact that supernova feedback is not instantaneous following star formation. The scatter in SFRs at constant gas surface density in a galaxy then depends on the properties of feedback and when we observe its star-forming regions at various points throughout their collapse/star formation "cycles". This has the following important observational consequences: (1) the scatter and normalization of the Kennicutt-Schmidt relation are relatively insensitive to the local (small-scale) star formation efficiency, (2) but gas depletion times and velocity dispersions are; (3) the scatter in and normalization of the Kennicutt-Schmidt relation is a sensitive probe of the feedback timescale and strength; (4) even in a model where $\tilde Q_{\rm gas}$ deterministically dictates star formation locally, time evolution, variation in local conditions (e.g., gas fractions and dynamical times), and variations between galaxies can destroy much of the observable correlation between SFR and $\tilde Q_{\rm gas}$ in resolved galaxy surveys. Additionally, this model exhibits large scatter in SFRs at low gas surface densities, in agreement with observations of flat outer HI disk velocity dispersion profiles.
[7]  oai:arXiv.org:1812.03997  [pdf] - 2025359
Cosmic Rays or Turbulence can Suppress Cooling Flows (Where Thermal Heating or Momentum Injection Fail)
Comments: 22 pages, 16 figures
Submitted: 2018-12-10
The quenching `maintenance' and `cooling flow' problems are important from the Milky Way through massive cluster elliptical galaxies. Previous work has shown that some source of energy beyond that from stars and pure magnetohydrodynamic processes is required, perhaps from AGN, but even the qualitative form of this energetic input remains uncertain. Different scenarios include thermal `heating,' direct wind or momentum injection, cosmic ray heating or pressure support, or turbulent `stirring' of the intra-cluster medium (ICM). We investigate these in $10^{12}-10^{14}\,{\rm M}_{\odot}$ halos using high-resolution non-cosmological simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, including simplified toy energy-injection models, where we arbitrarily vary the strength, injection scale, and physical form of the energy. We explore which scenarios can quench without violating observational constraints on energetics or ICM gas. We show that turbulent stirring in the central $\sim100\,$kpc, or cosmic-ray injection, can both maintain a stable low-SFR halo for $>$Gyr timescales with modest energy input, by providing a non-thermal pressure which stably lowers the core density and cooling rates. In both cases, associated thermal-heating processes are negligible. Turbulent stirring preserves cool-core features while mixing condensed core gas into the hotter halo and is by far the most energy efficient model. Pure thermal heating or nuclear isotropic momentum injection require vastly larger energy, are less efficient in lower-mass halos, easily over-heat cores, and require fine-tuning to avoid driving unphysical temperature gradients or gas expulsion from the halo center.
[8]  oai:arXiv.org:1702.06148  [pdf] - 1782667
FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
Comments: 64 pages, 40 figures. Simulation animations and visualizations available at http://www.tapir.caltech.edu/~phopkins/Site/animations and http://fire.northwestern.edu . Paper includes complete FIRE algorithms and public ICs (http://www.tapir.caltech.edu/~phopkins/publicICs). Updated to match published version
Submitted: 2017-02-20, last modified: 2018-11-11
The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.
[9]  oai:arXiv.org:1809.09120  [pdf] - 1897847
The failure of stellar feedback, magnetic fields, conduction, and morphological quenching in maintaining red galaxies
Comments: 16 pages, 12 figures
Submitted: 2018-09-24
The quenching "maintenance'" and related "cooling flow" problems are important in galaxies from Milky Way mass through clusters. We investigate this in halos with masses $\sim 10^{12}-10^{14}\,{\rm M}_{\odot}$, using non-cosmological high-resolution hydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model. We specifically focus on physics present without AGN, and show that various proposed "non-AGN" solution mechanisms in the literature, including Type Ia supernovae, shocked AGB winds, other forms of stellar feedback (e.g. cosmic rays), magnetic fields, Spitzer-Braginskii conduction, or "morphological quenching" do not halt or substantially reduce cooling flows nor maintain "quenched" galaxies in this mass range. We show that stellar feedback (including cosmic rays from SNe) alters the balance of cold/warm gas and the rate at which the cooled gas within the galaxy turns into stars, but not the net baryonic inflow. If anything, outflowing metals and dense gas promote additional cooling. Conduction is important only in the most massive halos, as expected, but even at $\sim 10^{14}\,{\rm M}_{\odot}$ reduces inflow only by a factor $\sim 2$ (owing to saturation effects and anisotropic suppression). Changing the morphology of the galaxies only slightly alters their Toomre-$Q$ parameter, and has no effect on cooling (as expected), so has essentially no effect on cooling flows or maintaining quenching. This all supports the idea that additional physics, e.g., AGN feedback, must be important in massive galaxies.
[10]  oai:arXiv.org:1701.01788  [pdf] - 1691943
What FIREs Up Star Formation: the Emergence of the Kennicutt-Schmidt Law from Feedback
Comments: 23 pages, 16 figures, accepted to MNRAS, fourth revision
Submitted: 2017-01-06, last modified: 2018-05-30
We present an analysis of the global and spatially-resolved Kennicutt-Schmidt (KS) star formation relation in the FIRE (Feedback In Realistic Environments) suite of cosmological simulations, including halos with $z = 0$ masses ranging from $10^{10}$ -- $10^{13}$ M$_{\odot}$. We show that the KS relation emerges and is robustly maintained due to the effects of feedback on local scales regulating star-forming gas, independent of the particular small-scale star formation prescriptions employed. We demonstrate that the time-averaged KS relation is relatively independent of redshift and spatial averaging scale, and that the star formation rate surface density is weakly dependent on metallicity and inversely dependent on orbital dynamical time. At constant star formation rate surface density, the `Cold \& Dense' gas surface density (gas with $T < 300$~K and $n > 10$~cm$^{-3}$, used as a proxy for the molecular gas surface density) of the simulated galaxies is $\sim$0.5~dex less than observed at $\sim$kpc scales. This discrepancy may arise from underestimates of the local column density at the particle-scale for the purposes of shielding in the simulations. Finally, we show that on scales larger than individual giant molecular clouds, the primary condition that determines whether star formation occurs is whether a patch of the galactic disk is thermally Toomre-unstable (not whether it is self-shielding): once a patch can no longer be thermally stabilized against fragmentation, it collapses, becomes self-shielding, cools, and forms stars, regardless of epoch or environment.
[11]  oai:arXiv.org:1712.02795  [pdf] - 1728835
Discrete Effects in Stellar Feedback: Individual Supernovae, Hypernovae, and IMF Sampling in Dwarf Galaxies
Comments: 9 pages, 4 figures
Submitted: 2017-12-07
Using high-resolution simulations from the FIRE-2 (Feedback In Realistic Environments) project, we study the effects of discreteness in stellar feedback processes on the evolution of galaxies and the properties of the interstellar medium (ISM). We specifically consider the discretization of supernovae (SNe), including hypernovae (HNe), and sampling the initial mass function (IMF). We study these processes in cosmological simulations of dwarf galaxies with $z=0$ stellar masses $M_{\ast}\sim 10^{4}-3\times10^{6}\,M_\odot$ (halo masses $\sim 10^{9}-10^{10}\,M_\odot$). We show that the discrete nature of individual SNe (as opposed to a model in which their energy/momentum deposition is continuous over time, similar to stellar winds) is crucial in generating a reasonable ISM structure and galactic winds and in regulating dwarf stellar masses. However, once SNe are discretized, accounting for the effects of IMF sampling on continuous mechanisms such as radiative feedback and stellar mass-loss (as opposed to adopting IMF-averaged rates) has weak effects on galaxy-scale properties. We also consider the effects of rare HNe events with energies $\sim 10^{53}\,{\rm erg}$. The effects of HNe are similar to the effects of clustered explosions of SNe -- which are already captured in our default simulation setup -- and do not quench star formation (provided that the HNe do not dominate the total SNe energy budget), which suggests that HNe yield products should be observable in ultra-faint dwarfs today.
[12]  oai:arXiv.org:1711.00178  [pdf] - 1608560
SDSS-IV MaNGA: Constraints on the Conditions for Star Formation in Galaxy Discs
Comments: 12 pages, 6 figures, accepted for publication in MNRAS
Submitted: 2017-10-31, last modified: 2017-11-19
Regions of disc galaxies with widespread star formation tend to be both gravitationally unstable and self-shielded against ionizing radiation, whereas extended outer discs with little or no star formation tend to be stable and unshielded on average. We explore what drives the transition between these two regimes, specifically whether discs first meet the conditions for self-shielding (parameterized by dust optical depth, $\tau$) or gravitational instability (parameterized by a modified version of Toomre's instability parameters, $Q_{\rm thermal}$, which quantifies the stability of a gas disc that is thermally supported at $T=10^4$ K). We first introduce a new metric formed by the product of these quantities, $Q_{\rm thermal}\tau$, which indicates whether the conditions for disk instability or self-shielding are easier to meet in a given region of a galaxy, and we discuss how $Q_{\rm thermal}\tau$ can be constrained even in the absence of direct gas information. We then analyse a sample of 13 galaxies with resolved gas measurements and find that on average galaxies will reach the threshold for disk instabilities ($Q_{\rm thermal}<1$) before reaching the threshold for self-shielding ($\tau>1$). Using integral field spectroscopic observations of a sample of 236 galaxies from the MaNGA survey, we find that the value of $Q_{\rm thermal}\tau$ in star-forming discs is consistent with similar behavior. These results support a scenario where disc fragmentation and collapse occurs before self-shielding, suggesting that gravitational instabilities are the primary condition for widespread star formation in galaxy discs. Our results support similar conclusions based on recent galaxy simulations.
[13]  oai:arXiv.org:1709.10099  [pdf] - 1589040
Stacked star formation rate profiles of bursty galaxies exhibit 'coherent' star formation
Comments: 7 pages, 4 figures, accepted for publication in ApJL
Submitted: 2017-09-28
In a recent work based on 3200 stacked H$\alpha$ maps of galaxies at $z \sim 1$, Nelson et al.~find evidence for `coherent star formation': the stacked SFR profiles of galaxies above (below) the 'star formation main sequence' (MS) are above (below) that of galaxies on the MS at all radii. One might interpret this result as inconsistent with highly bursty star formation and evidence that galaxies evolve smoothly along the MS rather than crossing it many times. We analyze six simulated galaxies at $z\sim1$ from the Feedback in Realistic Environments (FIRE) project in a manner analogous to the observations to test whether the above interpretations are correct. The trends in stacked SFR profiles are qualitatively consistent with those observed. However, SFR profiles of individual galaxies are much more complex than the stacked profiles: the former can be flat or even peak at large radii because of the highly clustered nature of star formation in the simulations. Moreover, the SFR profiles of individual galaxies above (below) the MS are not systematically above (below) those of MS galaxies at all radii. We conclude that the time-averaged coherent star formation evident stacks of observed galaxies is consistent with highly bursty, clumpy star formation of individual galaxies and is not evidence that galaxies evolve smoothly along the MS.
[14]  oai:arXiv.org:1608.01569  [pdf] - 1498394
Very High Energy outburst of Markarian 501 in May 2009
Aliu, E.; Archambault, S.; Archer, A.; Arlen, T.; Aune, T.; Barnacka, A.; Behera, B.; Beilicke, M.; Benbow, W.; Berger, K.; Bird, R.; Bouvier, A.; Böttcher, M.; Buchovecky, M.; Buckley, J. H.; Bugaev, V.; Cardenzana, J. V; Cerruti, M.; Cesarini, A.; Chen, X.; Ciupik, L.; Collins-Hughes, E.; Connolly, M. P.; Cui, W.; Dumm, J.; Eisch, J. D.; Falcone, A.; Federici, S.; Feng, Q.; Finley, J. P.; Fleischhack, H.; Fortin, P.; Fortson, L.; Furniss, A.; Galante, N.; Gall, D.; Gillanders, G. H.; Griffin, S.; Griffiths, S. T.; Grube, J.; Gyuk, G.; Hütten, M.; Håkansson, N.; Holder, J.; Hughes, G.; Humensky, T. B.; Johnson, C. A.; Kaaret, P.; Kar, P.; Kelley-Hoskins, N.; Kertzman, M.; Khassen, Y.; Kieda, D.; Krause, M.; Krawczynski, H.; Krennrich, F.; Lang, M. J.; Madhavan, A. S; Maier, G.; McArthur, S.; McCann, A.; Meagher, K.; Millis, J.; Moriarty, P.; Mukherjee, R.; Nieto, D.; de Bhróithe, A. O'Faoláin; Ong, R. A.; Orr, M.; Otte, A. N.; Pandel, D.; Park, N.; Pelassa, V.; Perkins, J. S.; Pichel, A.; Pohl, M.; Popkow, A.; Quinn, J.; Ragan, K.; Reyes, L. C.; Reynolds, P. T.; Roache, E.; Rousselle, J.; Rovero, A. C.; Saxon, D. B.; Sembroski, G. H.; Shahinyan, K.; Sheidaei, F.; Skole, C.; Smith, A. W.; Staszak, D.; Telezhinsky, I.; Theiling, M.; Todd, N. W.; Tucci, J. V.; Tyler, J.; Varlotta, A.; Vassiliev, V. V.; Vincent, S.; Wakely, S. P.; Weiner, O. M.; Weinstein, A.; Welsing, R.; Wilhelm, A.; Williams, D. A.; Zitzer, B.; Baring, M. G.; Gónzalez, J. Becerra; Cillis, A. N.; Horan, D.; Paneque, D.
Comments: 13 pages, 10 figures. Accepted for publication in Astronomy and Astrophysics
Submitted: 2016-08-04
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.
[15]  oai:arXiv.org:1202.2144  [pdf] - 1934751
VERITAS Deep Observations of the Dwarf Spheroidal Galaxy Segue 1
Comments: Added recently published erratum at the end of the original article
Submitted: 2012-02-09, last modified: 2015-07-07
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.
[16]  oai:arXiv.org:1411.1439  [pdf] - 1223034
VERITAS Observations of the BL Lac Object PG 1553+113
Comments: in press in ApJ, 2014
Submitted: 2014-11-05
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\%.
[17]  oai:arXiv.org:1409.1210  [pdf] - 1216673
Photon-Dominated Region Modeling of the [C I],[C II], and CO Line Emission from a Boundary in the Taurus Molecular Cloud
Comments: Accepted for publication in the Astrophysical Journal (Sept. 1, 2014), 11 pages, 7 figures
Submitted: 2014-09-03
We present [Ci] and [Cii] observations of a linear edge region in the Taurus molecular cloud, and model this region as a cylindrically symmetric PDR exposed to a low-intensity UV radiation field. The sharp, long profile of the linear edge makes it an ideal case to test PDR models and determine cloud parameters. We compare observations of the [C i], 3P1 -> 3P0 (492 GHz), [C i] 3P2 -> 3P1 (809 GHz), and [Cii] 2P3/2 -> 2P1/2 (1900 GHz) transitions, as well as the lowest rotational transitions of 12CO and 13CO, with line intensities produced by the RATRAN radiative transfer code from the results of the Meudon PDR code. We constrain the density structure of the cloud by fitting a cylindrical density function to visual extinction data. We study the effects of variation of the FUV field, 12C/13C isotopic abundance ratio, sulfur depletion, cosmic ray ionization rate, and inclination of the filament relative to the sky-plane on the chemical network of the PDR model and resulting line emission. We also consider the role of suprathermal chemistry and density inhomogeneities. We find good agreement between the model and observations, and that the integrated line intensities can be explained by a PDR model with an external FUV field of 0.05 G0, a low ratio of 12C to 13C ~ 43, a highly depleted sulfur abundance (by a factor of at least 50), a cosmic ray ionization rate (3 - 6) x 10-17 s^-1, and without significant effects from inclination, clumping or suprathermal chemistry.
[18]  oai:arXiv.org:1310.8150  [pdf] - 1180359
Observation of Markarian 421 in TeV gamma rays over a 14-year time span
Comments: 29 pages, 15 figures. Accepted for publication in Astroparticle Physics
Submitted: 2013-10-30
The variability of the blazar Markarian 421 in TeV gamma rays over a 14-year time period has been explored with the Whipple 10 m telescope. It is shown that the dynamic range of its flux variations is large and similar to that in X-rays. A correlation between the X-ray and TeV energy bands is observed during some bright flares and when the complete data sets are binned on long timescales. The main database consists of 878.4 hours of observation with the Whipple telescope, spread over 783 nights. The peak energy response of the telescope was 400 GeV with 20% uncertainty. This is the largest database of any TeV-emitting active galactic nucleus (AGN) and hence was used to explore the variability profile of Markarian 421. The time-averaged flux from Markarian 421 over this period was 0.446$\pm$0.008 Crab flux units. The flux exceeded 10 Crab flux units on three separate occasions. For the 2000-2001 season the average flux reached 1.86 Crab units, while in the 1996-1997 season the average flux was only 0.23 Crab units.
[19]  oai:arXiv.org:1310.5662  [pdf] - 737916
Snowmass Cosmic Frontiers 6 (CF6) Working Group Summary --The Bright Side of the Cosmic Frontier: Cosmic Probes of Fundamental Physics
Comments: Snowmass 2013 CF6 Working Group Report; author list corrected
Submitted: 2013-10-21, last modified: 2013-10-25
Report of the CF6 Working Group at Snowmass 2013. Topics addressed include ultra-high energy cosmic rays, neutrinos, gamma rays, baryogenesis, and experiments probing the fundamental nature of spacetime.
[20]  oai:arXiv.org:1310.6498  [pdf] - 1180199
Long term observations of B2 1215+30 with VERITAS
Comments: 34 pages, 5 figures, accepted for publication in ApJ
Submitted: 2013-10-24
We report on VERITAS observations of the BL Lac object B2 1215+30 between 2008 and 2012. During this period, the source was detected at very high energies (VHE; E > 100 GeV) by VERITAS with a significance of $8.9\sigma$ and showed clear variability on time scales larger than months. In 2011, the source was found to be in a relatively bright state and a power-law fit to the differential photon spectrum yields a spectral index of $3.6 \pm 0.4_{\mathrm{stat}} \pm 0.3_{\mathrm{syst}}$ with an integral flux above 200 GeV of $(8.0 \pm 0.9_{\mathrm{stat}} \pm 3.2_{\mathrm{syst}}) \times 10^{-12}\, \mathrm{cm}^{-2} \mathrm{s}^{-1}$. No short term variability could be detected during the bright state in 2011. Multi-wavelength data were obtained contemporaneous with the VERITAS observations in 2011 and cover optical (Super-LOTIS, MDM, Swift-UVOT), X-ray (Swift-XRT), and gamma-ray (Fermi-LAT) frequencies. These were used to construct the spectral energy distribution (SED) of B2 1215+30. A one-zone leptonic model is used to model the blazar emission and the results are compared to those of MAGIC from early 2011 and other VERITAS-detected blazars. The SED can be well reproduced with model parameters typical for VHE-detected BL Lacs.
[21]  oai:arXiv.org:1307.2232  [pdf] - 699247
CTA contributions to the 33rd International Cosmic Ray Conference (ICRC2013)
Consortium, The CTA; :; Abril, O.; Acharya, B. S.; Actis, M.; Agnetta, G.; Aguilar, J. A.; Aharonian, F.; Ajello, M.; Akhperjanian, A.; Alcubierre, M.; Aleksic, J.; Alfaro, R.; Aliu, E.; Allafort, A. J.; Allan, D.; Allekotte, I.; Aloisio, R.; Amato, E.; Ambrosi, G.; Ambrosio, M.; Anderson, J.; Angüner, E. O.; Antonelli, L. A.; Antonuccio, V.; Antonucci, M.; Antoranz, P.; Aravantinos, A.; Argan, A.; Arlen, T.; Aramo, C.; Armstrong, T.; Arnaldi, H.; Arrabito, L.; Asano, K.; Ashton, T.; Asorey, H. G.; Aune, T.; Awane, Y.; Baba, H.; Babic, A.; Baby, N.; Bähr, J.; Bais, A.; Baixeras, C.; Bajtlik, S.; Balbo, M.; Balis, D.; Balkowski, C.; Ballet, J.; Bamba, A.; Bandiera, R.; Barber, A.; Barbier, C.; Barceló, M.; Barnacka, A.; Barnstedt, J.; de Almeida, U. Barres; Barrio, J. A.; Basili, A.; Basso, S.; Bastieri, D.; Bauer, C.; Baushev, A.; Becciani, U.; Becerra, J.; Becerra, J.; Becherini, Y.; Bechtol, K. C.; Tjus, J. Becker; Beckmann, V.; Bednarek, W.; Behera, B.; Belluso, M.; Benbow, W.; Berdugo, J.; Berge, D.; Berger, K.; Bernard, F.; Bernardino, T.; Bernlöhr, K.; Bertucci, B.; Bhat, N.; Bhattacharyya, S.; Biasuzzi, B.; Bigongiari, C.; Biland, A.; Billotta, S.; Bird, T.; Birsin, E.; Bissaldi, E.; Biteau, J.; Bitossi, M.; Blake, S.; Bigas, O. Blanch; Blasi, P.; Bobkov, A.; Boccone, V.; Böttcher, M.; Bogacz, L.; Bogart, J.; Bogdan, M.; Boisson, C.; Gargallo, J. Boix; Bolmont, J.; Bonanno, G.; Bonardi, A.; Bonev, T.; Bonifacio, P.; Bonnoli, G.; Bordas, P.; Borgland, A.; Borkowski, J.; Bose, R.; Botner, O.; Bottani, A.; Bouchet, L.; Bourgeat, M.; Boutonnet, C.; Bouvier, A.; Brau-Nogué, S.; Braun, I.; Bretz, T.; Briggs, M.; Brigida, M.; Bringmann, T.; Britto, R.; Brook, P.; Brun, P.; Brunetti, L.; Bruno, P.; Bucciantini, N.; Buanes, T.; Buckley, J.; Bühler, R.; Bugaev, V.; Bulgarelli, A.; Bulik, T.; Busetto, G.; Buson, S.; Byrum, K.; Cailles, M.; Cameron, R.; Camprecios, J.; Canestrari, R.; Cantu, S.; Capalbi, M.; Caraveo, P.; Carmona, E.; Carosi, A.; Carosi, R.; Carr, J.; Carter, J.; Carton, P. -H.; Caruso, R.; Casanova, S.; Cascone, E.; Casiraghi, M.; Castellina, A.; Catalano, O.; Cavazzani, S.; Cazaux, S.; Cerchiara, P.; Cerruti, M.; Chabanne, E.; Chadwick, P.; Champion, C.; Chaves, R.; Cheimets, P.; Chen, A.; Chiang, J.; Chiappetti, L.; Chikawa, M.; Chitnis, V. R.; Chollet, F.; Christof, A.; Chudoba, J.; Cieślar, M.; Cillis, A.; Cilmo, M.; Codino, A.; Cohen-Tanugi, J.; Colafrancesco, S.; Colin, P.; Colome, J.; Colonges, S.; Compin, M.; Conconi, P.; Conforti, V.; Connaughton, V.; Conrad, J.; Contreras, J. L.; Coppi, P.; Coridian, J.; Corona, P.; Corti, D.; Cortina, J.; Cossio, L.; Costa, A.; Costantini, H.; Cotter, G.; Courty, B.; Couturier, S.; Covino, S.; Crimi, G.; Criswell, S. J.; Croston, J.; Cusumano, G.; Dafonseca, M.; Dale, O.; Daniel, M.; Darling, J.; Davids, I.; Dazzi, F.; de Angelis, A.; De Caprio, V.; De Frondat, F.; Pino, E. M. de Gouveia Dal; de la Calle, I.; De La Vega, G. A.; Lopez, R. de los Reyes; de Lotto, B.; De Luca, A.; de Naurois, M.; de Oliveira, Y.; Wilhelmi, E. de Oña; de Palma, F.; de Souza, V.; Decerprit, G.; Decock, G.; Deil, C.; Delagnes, E.; Deleglise, G.; Delgado, C.; della Volpe, D.; Demange, P.; Depaola, G.; Dettlaff, A.; Di Girolamo, T.; Di Giulio, C.; Di Paola, A.; Di Pierro, F.; di Sciascio, G.; Díaz, C.; Dick, J.; Dickherber, R.; Dickinson, H.; Diez-Blanco, V.; Digel, S.; Dimitrov, D.; Disset, G.; Djannati-Ataï, A.; Doert, M.; Dohmke, M.; Domainko, W.; Prester, D. Dominis; Donat, A.; Dorner, D.; Doro, M.; Dournaux, J. -L.; Drake, G.; Dravins, D.; Drury, L.; Dubois, F.; Dubois, R.; Dubus, G.; Dufour, C.; Dumas, D.; Dumm, J.; Durand, D.; Dwarkadas, V.; Dyks, J.; Dyrda, M.; Ebr, J.; Edy, E.; Egberts, K.; Eger, P.; Einecke, S.; Eleftheriadis, C.; Elles, S.; Emmanoulopoulos, D.; Engelhaupt, D.; Enomoto, R.; Ernenwein, J. -P.; Errando, M.; Etchegoyen, A.; Evans, P. A.; Falcone, A.; Faltenbacher, A.; Fantinel, D.; Farakos, K.; Farnier, C.; Farrell, E.; Fasola, G.; Favill, B. W.; Fede, E.; Federici, S.; Fegan, S.; Feinstein, F.; Ferenc, D.; Ferrando, P.; Fesquet, M.; Fetfatzis, P.; Fiasson, A.; Fillin-Martino, E.; Fink, D.; Finley, C.; Finley, J. P.; Fiorini, M.; Curcoll, R. Firpo; Flandrini, E.; Fleischhack, H.; Flores, H.; Florin, D.; Focke, W.; Föhr, C.; Fokitis, E.; Font, L.; Fontaine, G.; Fornasa, M.; Förster, A.; Fortson, L.; Fouque, N.; Franckowiak, A.; Franco, F. J.; Frankowski, A.; Fransson, C.; Fraser, G. W.; Frei, R.; Fresnillo, L.; Fruck, C.; Fugazza, D.; Fujita, Y.; Fukazawa, Y.; Fukui, Y.; Funk, S.; Gäbele, W.; Gabici, S.; Gabriele, R.; Gadola, A.; Galante, N.; Gall, D.; Gallant, Y.; Gámez-García, J.; Garczarczyk, M.; García, B.; López, R. Garcia; Gardiol, D.; Gargano, F.; Garrido, D.; Garrido, L.; Gascon, D.; Gaug, M.; Gaweda, J.; Gebremedhin, L.; Geffroy, N.; Gerard, L.; Ghedina, A.; Ghigo, M.; Ghislain, P.; Giannakaki, E.; Gianotti, F.; Giarrusso, S.; Giavitto, G.; Giebels, B.; Giglietto, N.; Gika, V.; Giomi, M.; Giommi, P.; Giordano, F.; Girard, N.; Giro, E.; Giuliani, A.; Glanzman, T.; Glicenstein, J. -F.; Godinovic, N.; Golev, V.; Berisso, M. Gomez; Gómez-Ortega, J.; Gonzalez, M. M.; González, A.; González, F.; Muñoz, A. González; Gothe, K. S.; Grabarczyk, T.; Gougerot, M.; Graciani, R.; Grandi, P.; Grañena, F.; Granot, J.; Grasseau, G.; Gredig, R.; Green, A.; Greenshaw, T.; Grégoire, T.; Grillo, A.; Grimm, O.; Grondin, M. -H.; Grube, J.; Grudzinska, M.; Gruev, V.; Grünewald, S.; Grygorczuk, J.; Guarino, V.; Gunji, S.; Gyuk, G.; Hadasch, D.; Hagedorn, A.; Hagiwara, R.; Hahn, J.; Hakansson, N.; Hallgren, A.; Heras, N. Hamer; Hara, S.; Hardcastle, M. J.; Harezlak, D.; Harris, J.; Hassan, T.; Hatanaka, K.; Haubold, T.; Haupt, A.; Hayakawa, T.; Hayashida, M.; Heller, R.; Henault, F.; Henri, G.; Hermann, G.; Hermel, R.; Herrero, A.; Hervet, O.; Hidaka, N.; Hinton, J. A.; Hirotani, K.; Hoffmann, D.; Hofmann, W.; Hofverberg, P.; Holder, J.; Hörandel, J. R.; Horns, D.; Horville, D.; Houles, J.; Hrabovsky, M.; Hrupec, D.; Huan, H.; Huber, B.; Huet, J. -M.; Hughes, G.; Humensky, T. B.; Huovelin, J.; Huppert, J. -F.; Ibarra, A.; Ikawa, D.; Illa, J. M.; Impiombato, D.; Incorvaia, S.; Inoue, S.; Inoue, Y.; Iocco, F.; Ioka, K.; Israel, G. L.; Jablonski, C.; Jacholkowska, A.; Jacquemier, J.; Jamrozy, M.; Janiak, M.; Jean, P.; Jeanney, C.; Jimenez, J. J.; Jogler, T.; Johnson, C.; Johnson, T.; Journet, L.; Juffroy, C.; Jung, I.; Kaaret, P.; Kabuki, S.; Kagaya, M.; Kakuwa, J.; Kalkuhl, C.; Kankanyan, R.; Karastergiou, A.; Kärcher, K.; Karczewski, M.; Karkar, S.; Kasperek, J.; Kastana, D.; Katagiri, H.; Kataoka, J.; Katarzyński, K.; Katz, U.; Kawanaka, N.; Kazanas, D.; Kelley-Hoskins, N.; Kellner-Leidel, B.; Kelly, H.; Kendziorra, E.; Khélifi, B.; Kieda, D. B.; Kifune, T.; Kihm, T.; Kishimoto, T.; Kitamoto, K.; Kluźniak, W.; Knapic, C.; Knapp, J.; Knödlseder, J.; Köck, F.; Kocot, J.; Kodani, K.; Köhne, J. -H.; Kohri, K.; Kokkotas, K.; Kolitzus, D.; Komin, N.; Kominis, I.; Konno, Y.; Köppel, H.; Korohoda, P.; Kosack, K.; Koss, G.; Kossakowski, R.; Koul, R.; Kowal, G.; Koyama, S.; Kozioł, J.; Krähenbühl, T.; Krause, J.; Krawzcynski, H.; Krennrich, F.; Krepps, A.; Kretzschmann, A.; Krobot, R.; Krueger, P.; Kubo, H.; Kudryavtsev, V. A.; Kushida, J.; Kuznetsov, A.; La Barbera, A.; La Palombara, N.; La Parola, V.; La Rosa, G.; Lacombe, K.; Lamanna, G.; Lande, J.; Languignon, D.; Lapington, J. S.; Laporte, P.; Laurent, B.; Lavalley, C.; Flour, T. Le; Padellec, A. Le; Lee, S. -H.; Lee, W. H.; Lefèvre, J. -P.; Leich, H.; de Oliveira, M. A. Leigui; Lelas, D.; Lenain, J. -P.; Leoni, R.; Leopold, D. J.; Lerch, T.; Lessio, L.; Leto, G.; Lieunard, B.; Lieunard, S.; Lindemann, R.; Lindfors, E.; Liolios, A.; Lipniacka, A.; Lockart, H.; Lohse, T.; Lombardi, S.; Longo, F.; Lopatin, A.; Lopez, M.; López-Coto, R.; López-Oramas, A.; Lorca, A.; Lorenz, E.; Louis, F.; Lubinski, P.; Lucarelli, F.; Lüdecke, H.; Ludwin, J.; Luque-Escamilla, P. L.; Lustermann, W.; Luz, O.; Lyard, E.; Maccarone, M. C.; Maccarone, T. J.; Madejski, G. M.; Madhavan, A.; Mahabir, M.; Maier, G.; Majumdar, P.; Malaguti, G.; Malaspina, G.; Maltezos, S.; Manalaysay, A.; Mancilla, A.; Mandat, D.; Maneva, G.; Mangano, A.; Manigot, P.; Mannheim, K.; Manthos, I.; Maragos, N.; Marcowith, A.; Mariotti, M.; Marisaldi, M.; Markoff, S.; Marszałek, A.; Martens, C.; Martí, J.; Martin, J. -M.; Martin, P.; Martínez, G.; Martínez, F.; Martínez, M.; Massaro, F.; Masserot, A.; Mastichiadis, A.; Mathieu, A.; Matsumoto, H.; Mattana, F.; Mattiazzo, S.; Maurer, A.; Maurin, G.; Maxfield, S.; Maya, J.; Mazin, D.; Comb, L. Mc; McCann, A.; McCubbin, N.; McHardy, I.; McKay, R.; Meagher, K.; Medina, C.; Melioli, C.; Melkumyan, D.; Melo, D.; Mereghetti, S.; Mertsch, P.; Meucci, M.; Meyer, M.; Michałowski, J.; Micolon, P.; Mihailidis, A.; Mineo, T.; Minuti, M.; Mirabal, N.; Mirabel, F.; Miranda, J. M.; Mirzoyan, R.; Mistò, A.; Mizuno, T.; Moal, B.; Moderski, R.; Mognet, I.; Molinari, E.; Molinaro, M.; Montaruli, T.; Monte, C.; Monteiro, I.; Moore, P.; Olaizola, A. Moralejo; Mordalska, M.; Morello, C.; Mori, K.; Morlino, G.; Morselli, A.; Mottez, F.; Moudden, Y.; Moulin, E.; Mrusek, I.; Mukherjee, R.; Munar-Adrover, P.; Muraishi, H.; Murase, K.; Murphy, A. StJ.; Nagataki, S.; Naito, T.; Nakajima, D.; Nakamori, T.; Nakayama, K.; Naumann, C.; Naumann, D.; Naumann-Godo, M.; Nayman, P.; Nedbal, D.; Neise, D.; Nellen, L.; Neronov, A.; Neustroev, V.; Neyroud, N.; Nicastro, L.; Nicolau-Kukliński, J.; Niedźwiecki, A.; Niemiec, J.; Nieto, D.; Nikolaidis, A.; Nishijima, K.; Nishikawa, K. -I.; Noda, K.; Nolan, S.; Northrop, R.; Nosek, D.; Nowak, N.; Nozato, A.; Oakes, L.; O'Brien, P. T.; Ohira, Y.; Ohishi, M.; Ohm, S.; Ohoka, H.; Okuda, T.; Okumura, A.; Olive, J. -F.; Ong, R. A.; Orito, R.; Orr, M.; Osborne, J. P.; Ostrowski, M.; Otero, L. A.; Otte, N.; Ovcharov, E.; Oya, I.; Ozieblo, A.; Padilla, L.; Pagano, I.; Paiano, S.; Paillot, D.; Paizis, A.; Palanque, S.; Palatka, M.; Pallota, J.; Palatiello, M.; Panagiotidis, K.; Panazol, J. -L.; Paneque, D.; Panter, M.; Panzera, M. R.; Paoletti, R.; Papayannis, A.; Papyan, G.; Paredes, J. M.; Pareschi, G.; Parraud, J. -M.; Parsons, D.; Pauletta, G.; Arribas, M. Paz; Pech, M.; Pedaletti, G.; Pelassa, V.; Pelat, D.; Perez, M. d. C.; Persic, M.; Petrucci, P. -O.; Peyaud, B.; Pichel, A.; Pieloth, D.; Pierre, E.; Pita, S.; Pivato, G.; Pizzolato, F.; Platino, M.; Platos, Ł.; Platzer, R.; Podkladkin, S.; Pogosyan, L.; Pohl, M.; Pojmanski, G.; Ponz, J. D.; Potter, W.; Poutanen, J.; Prandini, E.; Prast, J.; Preece, R.; Profeti, F.; Prokoph, H.; Prouza, M.; Proyetti, M.; Puerto-Giménez, I.; Pühlhofer, G.; Puljak, I.; Punch, M.; Pyzioł, R.; Quel, E. J.; Quesada, J.; Quinn, J.; Quirrenbach, A.; Racero, E.; Rainò, S.; Rajda, P. J.; Rameez, M.; Ramon, P.; Rando, R.; Rannot, R. C.; Rataj, M.; Raue, M.; Ravignani, D.; Reardon, P.; Reimann, O.; Reimer, A.; Reimer, O.; Reitberger, K.; Renaud, M.; Renner, S.; Reville, B.; Rhode, W.; Ribó, M.; Ribordy, M.; Richards, G.; Richer, M. G.; Rico, J.; Ridky, J.; Rieger, F.; Ringegni, P.; Ripken, J.; Ristori, P. R.; Rivière, A.; Rivoire, S.; Rob, L.; Rodeghiero, G.; Roeser, U.; Rohlfs, R.; Rojas, G.; Romano, P.; Romaszkan, W.; Romero, G. E.; Rosen, S. R.; Lees, S. Rosier; Ross, D.; Rouaix, G.; Rousselle, J.; Rousselle, S.; Rovero, A. C.; Roy, F.; Royer, S.; Rudak, B.; Rulten, C.; Rupiński, M.; Russo, F.; Ryde, F.; Saavedra, O.; Sacco, B.; Saemann, E. O.; Saggion, A.; Sahakian, V.; Saito, K.; Saito, T.; Saito, Y.; Sakaki, N.; Sakonaka, R.; Salini, A.; Sanchez, F.; Sanchez-Conde, M.; Sandoval, A.; Sandaker, H.; Sant'Ambrogio, E.; Santangelo, A.; Santos, E. M.; Sanuy, A.; Sapozhnikov, L.; Sarkar, S.; Sartore, N.; Sasaki, H.; Satalecka, K.; Sawada, M.; Scalzotto, V.; Scapin, V.; Scarcioffolo, M.; Schafer, J.; Schanz, T.; Schlenstedt, S.; Schlickeiser, R.; Schmidt, T.; Schmoll, J.; Schovanek, P.; Schroedter, M.; Schubert, A.; Schultz, C.; Schultze, J.; Schulz, A.; Schure, K.; Schussler, F.; Schwab, T.; Schwanke, U.; Schwarz, J.; Schwarzburg, S.; Schweizer, T.; Schwemmer, S.; Schwendicke, U.; Schwerdt, C.; Segreto, A.; Seiradakis, J. -H.; Sembroski, G. H.; Servillat, M.; Seweryn, K.; Sharma, M.; Shayduk, M.; Shellard, R. C.; Shi, J.; Shibata, T.; Shibuya, A.; Shore, S.; Shum, E.; Sideras-Haddad, E.; Sidoli, L.; Sidz, M.; Sieiro, J.; Sikora, M.; Silk, J.; Sillanpää, A.; Singh, B. B.; Sironi, G.; Sitarek, J.; Skole, C.; Smareglia, R.; Smith, A.; Smith, D.; Smith, J.; Smith, N.; Sobczyńska, D.; Sol, H.; Sottile, G.; Sowiński, M.; Spanier, F.; Spiga, D.; Spyrou, S.; Stamatescu, V.; Stamerra, A.; Starling, R. L. C.; Stawarz, Ł.; Steenkamp, R.; Stegmann, C.; Steiner, S.; Stella, C.; Stergioulas, N.; Sternberger, R.; Sterzel, M.; Stinzing, F.; Stodulski, M.; Stolarczyk, Th.; Straumann, U.; Strazzeri, E.; Stringhetti, L.; Suarez, A.; Suchenek, M.; Sugawara, R.; Sulanke, K. -H.; Sun, S.; Supanitsky, A. D.; Suric, T.; Sutcliffe, P.; Sykes, J. M.; Szanecki, M.; Szepieniec, T.; Szostek, A.; Tagliaferri, G.; Tajima, H.; Takahashi, H.; Takahashi, K.; Takalo, L.; Takami, H.; Talbot, G.; Tammi, J.; Tanaka, M.; Tanaka, S.; Tasan, J.; Tavani, M.; Tavernet, J. -P.; Tejedor, L. A.; Telezhinsky, I.; Temnikov, P.; Tenzer, C.; Terada, Y.; Terrier, R.; Teshima, M.; Testa, V.; Tezier, D.; Thayer, J.; Thuermann, D.; Tibaldo, L.; Tibaldo, L.; Tibolla, O.; Tiengo, A.; Timpanaro, M. C.; Tluczykont, M.; Peixoto, C. J. Todero; Tokanai, F.; Tokarz, M.; Toma, K.; Tonachini, A.; Torii, K.; Tornikoski, M.; Torres, D. F.; Torres, M.; Toscano, S.; Toso, G.; Tosti, G.; Totani, T.; Toussenel, F.; Tovmassian, G.; Travnicek, P.; Treves, A.; Trifoglio, M.; Troyano, I.; Tsinganos, K.; Ueno, H.; Umana, G.; Umehara, K.; Upadhya, S. S.; Usher, T.; Uslenghi, M.; Vagnetti, F.; Valdes-Galicia, J. F.; Vallania, P.; Vallejo, G.; van Driel, W.; van Eldik, C.; Vandenbrouke, J.; Vanderwalt, J.; Vankov, H.; Vasileiadis, G.; Vassiliev, V.; Veberic, D.; Vegas, I.; Vercellone, S.; Vergani, S.; Verzi, V.; Vettolani, G. P.; Veyssière, C.; Vialle, J. P.; Viana, A.; Videla, M.; Vigorito, C.; Vincent, P.; Vincent, S.; Vink, J.; Vlahakis, N.; Vlahos, L.; Vogler, P.; Voisin, V.; Vollhardt, A.; von Gunten, H. -P.; Vorobiov, S.; Vuerli, C.; Waegebaert, V.; Wagner, R.; Wagner, R. G.; Wagner, S.; Wakely, S. P.; Walter, R.; Walther, T.; Warda, K.; Warwick, R. S.; Wawer, P.; Wawrzaszek, R.; Webb, N.; Wegner, P.; Weinstein, A.; Weitzel, Q.; Welsing, R.; Werner, M.; Wetteskind, H.; White, R. J.; Wierzcholska, A.; Wiesand, S.; Wilhelm, A.; Wilkinson, M. I.; Williams, D. A.; Willingale, R.; Winde, M.; Winiarski, K.; Wischnewski, R.; Wiśniewski, Ł.; Wojcik, P.; Wood, M.; Wörnlein, A.; Xiong, Q.; Yadav, K. K.; Yamamoto, H.; Yamamoto, T.; Yamazaki, R.; Yanagita, S.; Yebras, J. M.; Yelos, D.; Yoshida, A.; Yoshida, T.; Yoshikoshi, T.; Yu, P.; Zabalza, V.; Zacharias, M.; Zajczyk, A.; Zampieri, L.; Zanin, R.; Zdziarski, A.; Zech, A.; Zhao, A.; Zhou, X.; Zietara, K.; Ziolkowski, J.; Ziółkowski, P.; Zitelli, V.; Zurbach, C.; Zychowski, P.
Comments: Index of CTA conference proceedings at the ICRC2013, Rio de Janeiro (Brazil). v1: placeholder with no arXiv links yet, to be replaced once individual contributions have been all submitted. v2: final with arXiv links to all CTA contributions and full author list
Submitted: 2013-07-08, last modified: 2013-07-29
Compilation of CTA contributions to the proceedings of the 33rd International Cosmic Ray Conference (ICRC2013), which took place in 2-9 July, 2013, in Rio de Janeiro, Brazil
[22]  oai:arXiv.org:1307.6772  [pdf] - 712876
Multiwavelength observations and modeling of 1ES 1959+650 in a low flux state
Comments: 9 pages, 7 figures. This is an author-created, un-copyedited version of an article accepted for publication in The Astrophysical Journal. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it
Submitted: 2013-07-25
We report on the VERITAS observations of the high-frequency peaked BL Lac object 1ES 1959+650 in the period 2007-2011. This source is detected at TeV energies by VERITAS at 16.4 standard deviation (sigma) significance in 7.6 hours of observation in a low flux state. A multiwavelength spectral energy distribution (SED) is constructed from contemporaneous data from VERITAS, Fermi-LAT, RXTE PCA, and Swift UVOT. Swift XRT data is not included in the SED due to a lack of simultaneous observations with VERITAS. In contrast to the orphan gamma-ray flare exhibited by this source in 2002, the X-ray flux of the source is found to vary by an order of magnitude, while other energy regimes exhibit less variable emission. A quasi-equilibrium synchrotron self-Compton model with an additional external radiation field is used to describe three SEDs corresponding to the lowest, highest, and average X-ray states. The variation in the X-ray spectrum is modeled by changing the electron injection spectral index, with minor adjustments of the kinetic luminosity in electrons. This scenario produces small-scale flux variability of order >~2 in the HE (E>1 MeV) and VHE (E>100 GeV) gamma-ray regimes, which is corroborated by the Fermi-LAT, VERITAS, and Whipple 10m telescope light curves.
[23]  oai:arXiv.org:1305.6508  [pdf] - 1171624
Discovery of TeV Gamma-ray Emission Toward Supernova Remnant SNR G78.2+2.1
Comments: 18 pages, 5 figures, accepted for publication in ApJ
Submitted: 2013-05-28
We report the discovery of an unidentified, extended source of very-high-energy (VHE) gamma-ray emission, VER J2019+407, within the radio shell of the supernova remnant SNR G78.2+2.1, using 21.4 hours of data taken by the VERITAS gamma-ray observatory in 2009. These data confirm the preliminary indications of gamma-ray emission previously seen in a two-year (2007-2009) blind survey of the Cygnus region by VERITAS. VER J2019+407, which is detected at a post-trials significance of 7.5 standard deviations in the 2009 data, is localized to the northwestern rim of the remnant in a region of enhanced radio and X-ray emission. It has an intrinsic extent of 0.23^{\circ} \pm 0.03^{\circ} (stat)+0.04^{\circ}_{-0.02}^{\circ}(sys) and its spectrum is well-characterized by a differential power law (dN/dE = N_0 \times (E/TeV)^{-\Gamma}) with a photon index of {\Gamma} = 2.37 \pm 0.14 (stat) \pm 0.20 (sys) and a flux normalization of N0 = 1.5 \pm 0.2 (stat) \pm 0.4(sys) \times 10^-12 ph TeV^{-1} cm^{-2} s^{-1}. This yields an integral flux of 5.2 \pm 0.8 (stat) \pm 1.4 (sys) \times 10^-12 ph cm^{-2} s^{-1} above 320 GeV, corresponding to 3.7% of the Crab Nebula flux. We consider the relationship of the TeV gamma-ray emission with the GeV gamma-ray emission seen from SNR G78.2+2.1 as well as that seen from a nearby cocoon of freshly accelerated cosmic rays. Multiple scenarios are considered as possible origins for the TeV gamma-ray emission, including hadronic particle acceleration at the supernova remnant shock.
[24]  oai:arXiv.org:1304.8057  [pdf] - 657808
The Extragalactic Background Light (EBL): A Probe of Fundamental Physics and a Record of Structure Formation in the Universe
Comments: Input for Snowmass 2013
Submitted: 2013-04-30
A brief white paper discussing the connections between the Extragalactic Background Light (EBL) and fundamental physics through gamma-ray observations is provided.
[25]  oai:arXiv.org:1303.2142  [pdf] - 1934854
VERITAS Observations of High-Frequency-Peaked BL Lac Objects
Comments: 2012 Fermi Symposium proceedings - eConf C121028
Submitted: 2013-03-08, last modified: 2013-03-15
Here we present highlights from VERITAS observations of high-frequency-peaked BL Lac objects (HBLs). We discuss the key science motivations for observing these sources -- including performing multiwavelength campaigns critical to understanding the emission mechanisms at work in HBLs, constraining the intensity and spectra shape of the extragalactic background light, and placing limits on the strength of the intergalactic magnetic field.
[26]  oai:arXiv.org:1212.4739  [pdf] - 1158614
Discovery of TeV Gamma-ray Emission from CTA 1 by VERITAS
Comments: 10 pages, 9 figures, accepted for publication in ApJ
Submitted: 2012-12-19
We report the discovery of TeV gamma-ray emission coincident with the shell-type radio supernova remnant (SNR) CTA 1 using the VERITAS gamma-ray observatory. The source, VER J0006+729, was detected as a 6.5 standard deviation excess over background and shows an extended morphology, approximated by a two-dimensional Gaussian of semi-major (semi-minor) axis 0.30 degree (0.24 degree) and a centroid 5' from the Fermi gamma-ray pulsar PSR J0007+7303 and its X-ray pulsar wind nebula (PWN). The photon spectrum is well described by a power-law dN/dE = N_0 (E/3 TeV)^(-\Gamma), with a differential spectral index of \Gamma = 2.2 +- 0.2_stat +- 0.3_sys, and normalization N_0 = (9.1 +- 1.3_stat +- 1.7_sys) x 10^(-14) cm^(-2) s^(-1) TeV^(-1). The integral flux, F_\gamma = 4.0 x 10^(-12) erg cm^(-2) s^(-1) above 1 TeV, corresponds to 0.2% of the pulsar spin-down power at 1.4 kpc. The energetics, co-location with the SNR, and the relatively small extent of the TeV emission strongly argue for the PWN origin of the TeV photons. We consider the origin of the TeV emission in CTA 1.
[27]  oai:arXiv.org:1211.3073  [pdf] - 1157743
Rapid TeV Gamma-Ray Flaring of BL Lacertae
Comments: Accepted for publication in The Astrophysical Journal
Submitted: 2012-11-13
We report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6-minute exposure, when the integral flux above 200 GeV reached $(3.4\pm0.6) \times 10^{-6} \;\text{photons}\;\text{m}^{-2}\text{s}^{-1}$, roughly 125% of the Crab Nebula flux measured by VERITAS. The light curve indicates that the observations missed the rising phase of the flare but covered a significant portion of the decaying phase. The exponential decay time was determined to be $13\pm4$ minutes, making it one of the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spectrum of BL Lacertae during the flare was soft, with a photon index of $3.6\pm 0.4$, which is in agreement with the measurement made previously by MAGIC in a lower flaring state. Contemporaneous radio observations of the source with the Very Long Baseline Array (VLBA) revealed the emergence of a new, superluminal component from the core around the time of the TeV gamma-ray flare, accompanied by changes in the optical polarization angle. Changes in flux also appear to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of the flare, although they are difficult to quantify precisely due to sparse coverage. A strong flare was seen at radio wavelengths roughly four months later, which might be related to the gamma-ray flaring activities. We discuss the implications of these multiwavelength results.
[28]  oai:arXiv.org:1210.7224  [pdf] - 1152495
VERITAS Observations of Six Bright, Hard-Spectrum Fermi-LAT Blazars
Comments:
Submitted: 2012-10-26
We report on VERITAS very-high-energy (VHE; E>100 GeV) observations of six blazars selected from the Fermi Large Area Telescope First Source Catalog (1FGL). The gamma-ray emission from 1FGL sources was extrapolated up to the VHE band, taking gamma-ray absorption by the extragalactic background light into account. This allowed the selection of six bright, hard-spectrum blazars that were good candidate TeV emitters. Spectroscopic redshift measurements were attempted with the Keck Telescope for the targets without Sloan Digital Sky Survey (SDSS) spectroscopic data. No VHE emission is detected during the observations of the six sources described here. Corresponding TeV upper limits are presented, along with contemporaneous Fermi observations and non-concurrent Swift UVOT and XRT data. The blazar broadband spectral energy distributions (SEDs) are assembled and modeled with a single-zone synchrotron self-Compton model. The SED built for each of the six blazars show a synchrotron peak bordering between the intermediate- and high-spectrum-peak classifications, with four of the six resulting in particle-dominated emission regions.
[29]  oai:arXiv.org:1210.4786  [pdf] - 1152222
Search for a correlation between very-high-energy gamma rays and giant radio pulses in the Crab pulsar
Comments: 12 pages, 12 figures, To appear in The Astrophysical Journal
Submitted: 2012-10-17
We present the results of a joint observational campaign between the Green Bank radio telescope and the VERITAS gamma-ray telescope, which searched for a correlation between the emission of very-high-energy (VHE) gamma rays ($E_{\gamma} >$ 150 GeV) and Giant Radio Pulses (GRPs) from the Crab pulsar at 8.9 GHz. A total of 15366 GRPs were recorded during 11.6 hours of simultaneous observations, which were made across four nights in December 2008 and in November and December 2009. We searched for an enhancement of the pulsed gamma-ray emission within time windows placed around the arrival time of the GRP events. In total, 8 different time windows with durations ranging from 0.033 ms to 72 s were positioned at three different locations relative to the GRP to search for enhanced gamma-ray emission which lagged, led, or was concurrent with, the GRP event. Further, we performed separate searches on main pulse GRPs and interpulse GRPs and on the most energetic GRPs in our data sample. No significant enhancement of pulsed VHE emission was found in any of the preformed searches. We set upper limits of 5-10 times the average VHE flux of the Crab pulsar on the flux simultaneous with interpulse GRPs on single-rotation-period time scales. On $\sim$8-second time scales around interpulse GRPs, we set an upper limit of 2-3 times the average VHE flux. Within the framework of recent models for pulsed VHE emission from the Crab pulsar, the expected VHE-GRP emission correlations are below the derived limits.
[30]  oai:arXiv.org:1208.0676  [pdf] - 1937794
Constraints on Cosmic Rays, Magnetic Fields, and Dark Matter from Gamma-Ray Observations of the Coma Cluster of Galaxies with VERITAS and Fermi
Comments: ApJ in press
Submitted: 2012-08-03
Observations of radio halos and relics in galaxy clusters indicate efficient electron acceleration. Protons should likewise be accelerated, suggesting that clusters may also be sources of very high-energy (VHE; E>100 GeV) gamma-ray emission. We report here on VHE gamma-ray observations of the Coma galaxy cluster with the VERITAS array of imaging Cherenkov telescopes, with complementing Fermi-LAT observations at GeV energies. No significant gamma-ray emission from the Coma cluster was detected. Integral flux upper limits at the 99% confidence level were measured to be on the order of (2-5)*10^-8\ ph. m^-2 s^-1 (VERITAS, >220 GeV} and ~2*10^-6 ph. m^-2 s^-1 (Fermi, 1-3 GeV), respectively. We use the gamma-ray upper limits to constrain CRs and magnetic fields in Coma. Using an analytical approach, the CR-to-thermal pressure ratio is constrained to be < 16% from VERITAS data and < 1.7% from Fermi data (averaged within the virial radius). These upper limits are starting to constrain the CR physics in self-consistent cosmological cluster simulations and cap the maximum CR acceleration efficiency at structure formation shocks to be <50%. Assuming that the radio-emitting electrons of the Coma halo result from hadronic CR interactions, the observations imply a lower limit on the central magnetic field in Coma of (2 - 5.5) muG, depending on the radial magnetic-field profile and on the gamma-ray spectral index. Since these values are below those inferred by Faraday rotation measurements in Coma (for most of the parameter space), this {renders} the hadronic model a very plausible explanation of the Coma radio halo. Finally, since galaxy clusters are dark-matter (DM) dominated, the VERITAS upper limits have been used to place constraints on the thermally-averaged product of the total self-annihilation cross section and the relative velocity of the DM particles, <\sigma v>. (abr.)
[31]  oai:arXiv.org:1206.4080  [pdf] - 1124213
Multiwavelength observations of the AGN 1ES 0414+009 with VERITAS, Fermi-LAT, Swift-XRT, and MDM
Comments: Accepted for publication in The Astrophysical Journal
Submitted: 2012-06-18, last modified: 2012-06-20
We present observations of the BL Lac object 1ES 0414+009 in the >200 GeV gamma-ray band by the VERITAS array of Cherenkov telescopes. 1ES 0414+009 was observed by VERITAS between January 2008 and February 2011, resulting in 56.2 hours of good quality pointed observations. These observations resulted in a detection of 822 events from the source corresponding to a statistical significance of 6.4 standard deviations (6.4 sigma) above the background. The source flux, showing no evidence for variability, is measured as 5.2 +/- 1.1_stat +/- 2.6_sys * 10^-12 photons cm^-2 s^-1 above 200 GeV, equivalent to approximately 2% of the Crab Nebula flux above this energy. The differential photon spectrum from 230 GeV to 850 GeV is well fit by a power law with an photon index of Gamma 3.4 +/- 0.5_stat +/- 0.3_sys and a flux normalization of 1.6 +/- 0.3_stat +/- 0.8_sys * 10^-11 photons cm^-2 s^-1 at 300 GeV. We also present multiwavelength results taken in the optical (MDM), X-ray (Swift-XRT), and GeV (Fermi-LAT) bands and use these results to construct a broadband spectral energy distribution (SED). Modeling of this SED indicates that homogenous one-zone leptonic scenarios are not adequate to describe emission from the system, with a lepto-hadronic model providing a better fit to the data.
[32]  oai:arXiv.org:1205.5287  [pdf] - 1123624
VERITAS Observations of the Nova in V407 Cygni
Comments: 17 pages, 4 figures, accepted for publication in ApJ
Submitted: 2012-05-23
We report on very high energy (E > 100 GeV) gamma-ray observations of V407 Cygni, a symbiotic binary that underwent a nova outburst producing 0.1-10 GeV gamma rays during 2010 March 10-26. Observations were made with the Very Energetic Radiation Imaging Telescope Array System during 2010 March 19-26 at relatively large zenith angles, due to the position of V407 Cyg. An improved reconstruction technique for large zenith angle observations is presented and used to analyze the data. We do not detect V407 Cygni and place a differential upper limit on the flux at 1.6 TeV of 2.3 \times 10^(-12) erg cm^(-2) s^(-1) (at the 95% confidence level). When considered jointly with data from Fermi-LAT, this result places limits on the acceleration of very high energy particles in the nova.
[33]  oai:arXiv.org:1204.0865  [pdf] - 1117757
Discovery of High-energy and Very High Energy Gamma-ray Emission from the Blazar RBS 0413
Comments: 14 pages, 3 figures, accepted to ApJ, corresponding author: G. D. Senturk
Submitted: 2012-04-04
We report on the discovery of high-energy (HE; E > 0.1 GeV) and very high-energy (VHE; E > 100 GeV) gamma-ray emission from the high-frequency-peaked BL Lac object RBS 0413. VERITAS, a ground-based gamma-ray observatory, detected VHE gamma rays from RBS 0413 with a statistical significance of 5.5 standard deviations (sigma) and a gamma-ray flux of (1.5 \pm 0.6stat \pm 0.7syst) \times 10^(-8) photons m^(-2) s^(-1) (\sim 1% of the Crab Nebula flux) above 250 GeV. The observed spectrum can be described by a power law with a photon index of 3.18 \pm 0.68stat \pm 0.30syst. Contemporaneous observations with the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope detected HE gamma rays from RBS 0413 with a statistical significance of more than 9 sigma, a power-law photon index of 1.57 \pm 0.12stat +0.11sys -0.12sys and a gamma-ray flux between 300 MeV and 300 GeV of (1.64 \pm 0.43stat +0.31sys -0.22sys) \times 10^(-5) photons m^(-2) s^(-1). We present the results from Fermi-LAT and VERITAS, including a spectral energy distribution modeling of the gamma-ray, quasi-simultaneous X-ray (Swift-XRT), ultraviolet (Swift-UVOT) and R-band optical (MDM) data. We find that, if conditions close to equipartition are required, both the combined synchrotron self-Compton/external-Compton and the lepto-hadronic models are preferred over a pure synchrotron self-Compton model.
[34]  oai:arXiv.org:1111.5341  [pdf] - 554253
The 2010 very high energy gamma-ray flare & 10 years of multi-wavelength observations of M 87
Collaboration, The H. E. S. S.; :; Abramowski, A.; Acero, F.; Aharonian, F.; Akhperjanian, A. G.; Anton, G.; Balzer, A.; Barnacka, A.; de Almeida, U. Barres; Becherini, Y.; Becker, J.; Behera, B.; Bernlöhr, K.; Birsin, E.; Biteau, J.; Bochow, A.; Boisson, C.; Bolmont, J.; Bordas, P.; Brucker, J.; Brun, F.; Brun, P.; Bulik, T.; Büsching, I.; Carrigan, S.; Casanova, S.; Cerruti, M.; Chadwick, P. M.; Charbonnier, A.; Chaves, R. C. G.; Cheesebrough, A.; Clapson, A. C.; Coignet, G.; Cologna, G.; Conrad, J.; Dalton, M.; Daniel, M. K.; Davids, I. D.; Degrange, B.; Deil, C.; Dickinson, H. J.; Djannati-Ataï, A.; Domainko, W.; Drury, L. O'C.; Dubus, G.; Dutson, K.; Dyks, J.; Dyrda, M.; Egberts, K.; Eger, P.; Espigat, P.; Fallon, L.; Farnier, C.; Fegan, S.; Feinstein, F.; Fernandes, M. V.; Fiasson, A.; Fontaine, G.; Förster, A.; Füßling, M.; Gallant, Y. A.; Gast, H.; Gérard, L.; Gerbig, D.; Giebels, B.; Glicenstein, J. F.; Glück, B.; Goret, P.; Göring, D.; Häffner, S.; Hague, J. D.; Hampf, D.; Hauser, M.; Heinz, S.; Heinzelmann, G.; Henri, G.; Hermann, G.; Hinton, J. A.; Hoffmann, A.; Hofmann, W.; Hofverberg, P.; Holler, M.; Horns, D.; Jacholkowska, A.; de Jager, O. C.; Jahn, C.; Jamrozy, M.; Jung, I.; Kastendieck, M. A.; Katarzyński, K.; Katz, U.; Kaufmann, S.; Keogh, D.; Khangulyan, D.; Khélifi, B.; Klochkov, D.; Kluźniak, W.; Kneiske, T.; Komin, Nu.; Kosack, K.; Kossakowski, R.; Laffon, H.; Lamanna, G.; Lennarz, D.; Lohse, T.; Lopatin, A.; Lu, C. -C.; Marandon, V.; Marcowith, A.; Masbou, J.; Maurin, D.; Maxted, N.; Mayer, M.; McComb, T. J. L.; Medina, M. C.; Méhault, J.; Moderski, R.; Moulin, E.; Naumann, C. L.; Naumann-Godo, M.; de Naurois, M.; Nedbal, D.; Nekrassov, D.; Nguyen, N.; Nicholas, B.; Niemiec, J.; Nolan, S. J.; Ohm, S.; Wilhelmi, E. de Oña; Opitz, B.; Ostrowski, M.; Oya, I.; Panter, M.; Arribas, M. Paz; Pedaletti, G.; Pelletier, G.; Petrucci, P. -O.; Pita, S.; Pühlhofer, G.; Punch, M.; Quirrenbach, A.; Raue, M.; Rayner, S. M.; Reimer, A.; Reimer, O.; Renaud, M.; Reyes, R. de los; Rieger, F.; Ripken, J.; Rob, L.; Rosier-Lees, S.; Rowell, G.; Rudak, B.; Rulten, C. B.; Ruppel, J.; Sahakian, V.; Sanchez, D. A.; Santangelo, A.; Schlickeiser, R.; Schöck, F. M.; Schulz, A.; Schwanke, U.; Schwarzburg, S.; Schwemmer, S.; Sheidaei, F.; Skilton, J. L.; Sol, H.; Spengler, G.; Stawarz, Ł.; Steenkamp, R.; Stegmann, C.; Stinzing, F.; Stycz, K.; Sushch, I.; Szostek, A.; Tavernet, J. -P.; Terrier, R.; Tluczykont, M.; Valerius, K.; van Eldik, C.; Vasileiadis, G.; Venter, C.; Vialle, J. P.; Viana, A.; Vincent, P.; Völk, H. J.; Volpe, F.; Vorobiov, S.; Vorster, M.; Wagner, S. J.; Ward, M.; White, R.; Wierzcholska, A.; Zacharias, M.; Zajczyk, A.; Zdziarski, A. A.; Zech, A.; Zechlin, H. -S.; Collaboration, The MAGIC; :; Aleksić, J.; Antonelli, L. A.; Antoranz, P.; Backes, M.; Barrio, J. A.; Bastieri, D.; González, J. Becerra; Bednarek, W.; Berdyugin, A.; Berger, K.; Bernardini, E.; Biland, A.; Blanch, O.; Bock, R. K.; Boller, A.; Bonnoli, G.; Tridon, D. Borla; Braun, I.; Bretz, T.; Cañellas, A.; Carmona, E.; Carosi, A.; Colin, P.; Colombo, E.; Contreras, J. L.; Cortina, J.; Cossio, L.; Covino, S.; Dazzi, F.; De Angelis, A.; del Pozo, E. De Cea; De Lotto, B.; Mendez, C. Delgado; Ortega, A. Diago; Doert, M.; Domínguez, A.; Prester, D. Dominis; Dorner, D.; Doro, M.; Elsaesser, D.; Ferenc, D.; Fonseca, M. V.; Font, L.; Fruck, C.; López, R. J. García; Garczarczyk, M.; Garrido, D.; Giavitto, G.; Godinović, N.; Hadasch, D.; Häfner, D.; Herrero, A.; Hildebrand, D.; Höhne-Mönch, D.; Hose, J.; Hrupec, D.; Huber, B.; Jogler, T.; Klepser, S.; Krähenbühl, T.; Krause, J.; La Barbera, A.; Lelas, D.; Leonardo, E.; Lindfors, E.; Lombardi, S.; López, M.; Lorenz, E.; Makariev, M.; Maneva, G.; Mankuzhiyil, N.; Mannheim, K.; Maraschi, L.; Mariotti, M.; Martínez, M.; Mazin, D.; Meucci, M.; Miranda, J. M.; Mirzoyan, R.; Miyamoto, H.; Moldón, J.; Moralejo, A.; Munar, P.; Nieto, D.; Nilsson, K.; Orito, R.; Oya, I.; Paneque, D.; Paoletti, R.; Pardo, S.; Paredes, J. M.; Partini, S.; Pasanen, M.; Pauss, F.; Perez-Torres, M. A.; Persic, M.; Peruzzo, L.; Pilia, M.; Pochon, J.; Prada, F.; Moroni, P. G. Prada; Prandini, E.; Puljak, I.; Reichardt, I.; Reinthal, R.; Rhode, W.; Ribó, M.; Rico, J.; Rügamer, S.; Saggion, A.; Saito, K.; Saito, T. Y.; Salvati, M.; Satalecka, K.; Scalzotto, V.; Scapin, V.; Schultz, C.; Schweizer, T.; Shayduk, M.; Shore, S. N.; Sillanpää, A.; Sitarek, J.; Sobczynska, D.; Spanier, F.; Spiro, S.; Stamerra, A.; Steinke, B.; Storz, J.; Strah, N.; Surić, T.; Takalo, L.; Takami, H.; Tavecchio, F.; Temnikov, P.; Terzić, T.; Tescaro, D.; Teshima, M.; Thom, M.; Tibolla, O.; Torres, D. F.; Treves, A.; Vankov, H.; Vogler, P.; Wagner, R. M.; Weitzel, Q.; Zabalza, V.; Zandanel, F.; Zanin, R.; Collaboration, The VERITAS; :; Arlen, T.; Aune, T.; Beilicke, M.; Benbow, W.; Bouvier, A.; Bradbury, S. M.; Buckley, J. H.; Bugaev, V.; Byrum, K.; Cannon, A.; Cesarini, A.; Ciupik, L.; Connolly, M. P.; Cui, W.; Dickherber, R.; Duke, C.; Errando, M.; Falcone, A.; Finley, J. P.; Finnegan, G.; Fortson, L.; Furniss, A.; Galante, N.; Gall, D.; Godambe, S.; Griffin, S.; Grube, J.; Gyuk, G.; Hanna, D.; Holder, J.; Huan, H.; Hui, C. M.; Kaaret, P.; Karlsson, N.; Kertzman, M.; Khassen, Y.; Kieda, D.; Krawczynski, H.; Krennrich, F.; Lang, M. J.; LeBohec, S.; Maier, G.; McArthur, S.; McCann, A.; Moriarty, P.; Mukherjee, R.; Nuñez, P. D.; Ong, R. A.; Orr, M.; Otte, A. N.; Park, N.; Perkins, J. S.; Pichel, A.; Pohl, M.; Prokoph, H.; Ragan, K.; Reyes, L. C.; Reynolds, P. T.; Roache, E.; Rose, H. J.; Ruppel, J.; Schroedter, M.; Sembroski, G. H.; Şentürk, G. D.; Telezhinsky, I.; Tešić, G.; Theiling, M.; Thibadeau, S.; Varlotta, A.; Vassiliev, V. V.; Vivier, M.; Wakely, S. P.; Weekes, T. C.; Williams, D. A.; Zitzer, B.; de Almeida, U. Barres; Cara, M.; Casadio, C.; Cheung, C. C.; McConville, W.; Davies, F.; Doi, A.; Giovannini, G.; Giroletti, M.; Hada, K.; Hardee, P.; Harris, D. E.; Junor, W.; Kino, M.; Lee, N. P.; Ly, C.; Madrid, J.; Massaro, F.; Mundell, C. G.; Nagai, H.; Perlman, E. S.; Steele, I. A.; Walker, R. C.; Wood, D. L.
Comments: 19 pages, 5 figures; Corresponding authors: M. Raue, L. Stawarz, D. Mazin, P. Colin, C. M. Hui, M. Beilicke; Fig. 1 lightcurve data available online: http://www.desy.de/~mraue/m87/
Submitted: 2011-11-22, last modified: 2012-02-20
Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.
[35]  oai:arXiv.org:1112.4518  [pdf] - 1092506
VERITAS Observations of day-scale flaring of M87 in 2010 April
Comments: Accepted for publication in ApJ; 17 pages, 5 figures; Corresponding author: C. M. Hui
Submitted: 2011-12-19
VERITAS has been monitoring the very-high-energy (VHE; >100GeV) gamma-ray activity of the radio galaxy M87 since 2007. During 2008, flaring activity on a timescale of a few days was observed with a peak flux of (0.70 +- 0.16) X 10^{-11} cm^{-2} s^{-1} at energies above 350GeV. In 2010 April, VERITAS detected a flare from M87 with peak flux of (2.71 +- 0.68) X 10^{-11} cm^{-2} s^{-1} for E>350GeV. The source was observed for six consecutive nights during the flare, resulting in a total of 21 hr of good quality data. The most rapid flux variation occurred on the trailing edge of the flare with an exponential flux decay time of 0.90^{+0.22}_{-0.15} days. The shortest detected exponential rise time is three times as long, at 2.87^{+1.65}_{-0.99} days. The quality of the data sample is such that spectral analysis can be performed for three periods: rising flux, peak flux, and falling flux. The spectra obtained are consistent with power-law forms. The spectral index at the peak of the flare is equal to 2.19 +- 0.07. There is some indication that the spectrum is softer in the falling phase of the flare than the peak phase, with a confidence level corresponding to 3.6 standard deviations. We discuss the implications of these results for the acceleration and cooling rates of VHE electrons in M87 and the constraints they provide on the physical size of the emitting region.
[36]  oai:arXiv.org:1109.0050  [pdf] - 1083747
VERITAS Observations of Gamma-Ray Bursts Detected by Swift
Comments: 11 pages, 6 figures
Submitted: 2011-08-31, last modified: 2011-11-25
We present the results of sixteen Swift-triggered GRB follow-up observations taken with the VERITAS telescope array from January, 2007 to June, 2009. The median energy threshold and response time of these observations was 260 GeV and 320 s, respectively. Observations had an average duration of 90 minutes. Each burst is analyzed independently in two modes: over the whole duration of the observations and again over a shorter time scale determined by the maximum VERITAS sensitivity to a burst with a t^-1.5 time profile. This temporal model is characteristic of GRB afterglows with high-energy, long-lived emission that have been detected by the Large Area Telescope (LAT) on-board the Fermi satellite. No significant VHE gamma-ray emission was detected and upper limits above the VERITAS threshold energy are calculated. The VERITAS upper limits are corrected for gamma-ray extinction by the extragalactic background light (EBL) and interpreted in the context of the keV emission detected by Swift. For some bursts the VHE emission must have less power than the keV emission, placing constraints on inverse Compton models of VHE emission.
[37]  oai:arXiv.org:1111.2390  [pdf] - 436777
VERITAS Collaboration Contributions to the 32nd International Cosmic Ray Conference
Comments: 1 html page. ICRC 2011, Beijing, China. A full list of authors and collaborating institutions can be found at http://veritas.sao.arizona.edu/conferences/authors?icrc2011
Submitted: 2011-11-09
Compilation of papers contributed by the VERITAS Collaboration to the 32nd International Cosmic Ray Conference, held 11-18 August 2011 in Beijing, China.
[38]  oai:arXiv.org:1110.5949  [pdf] - 1085180
Multiwavelength Observations of the Previously Unidentified Blazar RX J0648.7+1516
Comments:
Submitted: 2011-10-26, last modified: 2011-11-08
We report on the VERITAS discovery of very-high-energy (VHE) gamma- ray emission above 200 GeV from the high-frequency-peaked BL Lac object RXJ0648.7+1516 (GBJ0648+1516), associated with 1FGLJ0648.8+1516. The photon spectrum above 200 GeV is fit by a power law dN/dE = F0(E/E0)-{\Gamma} with a photon index {\Gamma} of 4.4 {\pm} 0.8stat {\pm}0.3syst and a flux normalization F0 of (2.3 {\pm}0.5stat {\pm} 1.2sys) {\times} 10-11 TeV-1cm-2s-1 with E0 = 300 GeV. No VHE vari- ability is detected during VERITAS observations of RXJ0648.7+1516 between 2010 March 4 and April 15. Following the VHE discovery, the optical identifica- tion and spectroscopic redshift were obtained using the Shane 3-m Telescope at the Lick Observatory, showing the unidentified object to be a BL Lac type with a redshift of z = 0.179. Broadband multiwavelength observations contemporaneous with the VERITAS exposure period can be used to sub-classify the blazar as a high-frequency-peaked BL Lac (HBL) object, including data from the MDM ob- servatory, Swift-UVOT and XRT, and continuous monitoring at photon energies above 1 GeV from the Fermi Large Area Telescope (LAT). We find that in the absence of undetected, high-energy rapid variability, the one-zone synchrotron self-Compton model (SSC) overproduces the high-energy gamma-ray emission measured by the Fermi-LAT over 2.3 years. The SED can be parameterized sat- isfactorily with an external-Compton or lepto-hadronic model, which have two and six additional free parameters, respectively, compared to the one-zone SSC model.
[39]  oai:arXiv.org:1111.1225  [pdf] - 1579386
VERITAS: Status and Highlights
Comments: Presented at the 32nd ICRC, Beijing, 2011
Submitted: 2011-11-04
The VERITAS telescope array has been operating smoothly since 2007, and has detected gamma-ray emission above 100 GeV from 40 astrophysical sources. These include blazars, pulsar wind nebulae, supernova remnants, gamma-ray binary systems, a starburst galaxy, a radio galaxy, the Crab pulsar, and gamma-ray sources whose origin remains unidentified. In 2009, the array was reconfigured, greatly improving the sensitivity. We summarize the current status of the observatory, describe some of the scientific highlights since 2009, and outline plans for the future.
[40]  oai:arXiv.org:1110.4643  [pdf] - 1934669
VERITAS Observations of the BL Lac Object PG 1553+113 Between May 2010 and May 2011
Comments: 4 pages, to appear in proceedings of the 32nd ICRC, Beijing, China
Submitted: 2011-10-20
VERITAS, an array of imaging atmospheric-Cherenkov telescopes, studies blazars in the energy range between ~100 GeV and ~30 TeV. With its excellent sensitivity at these energies, and ever-deepening source exposures, VERITAS is in a position to probe distant blazars for clear absorption signatures in their very-high-energy gamma-ray spectra due to interactions with the extragalactic background light (EBL). We discuss results from recent VERITAS observations of PG 1553+113 (z > 0.4) which have resulted in the most significant very-high-energy detection ever obtained for this source. The most recent VERITAS spectral measurements are used to place an upper limit on the source redshift of z < 0.5 at the 95% confidence level. Also discussed are the prospects of using these observations, along with those of other hard- spectrum blazars, to place constraints on the EBL.
[41]  oai:arXiv.org:1110.1036  [pdf] - 1084659
Automated analysis of Fermi-LAT data to trigger ground-based gamma-ray observations
Comments: 4 pages, 4 figures, to appear in proceedings of the 32nd ICRC, Beijing, China
Submitted: 2011-10-05
The Fermi Gamma-ray Space Telescope regularly surveys the entire sky in the energy range between 0.3 and 100 GeV with an homogeneous coverage. This makes Fermi a very useful guide for ground-based Cherenkov-telescope arrays like VERITAS, that are sensitive at energies above 100 GeV. The VERITAS collaboration uses information from Fermi-LAT to select potential targets and identify flaring objects in the GeV band. We describe three different analysis pipelines that automatically process and analyze Fermi-LAT data on a daily basis to look for flaring objects: daily light curves of selected targets, an all-sky analysis to identify new flaring sources, and an algorithm that looks for clusters of high-energy photons in time and arrival direction.
[42]  oai:arXiv.org:1108.3797  [pdf] - 428636
Detection of Pulsed Gamma Rays Above 100 GeV from the Crab Pulsar
Comments: 17 pages, corresponding authors: Nepomuk Otte nepomuk.otte@gmail.com; Andrew McCann mccann@hep.physics.mcgill.ca; Martin Schroedter schroedter@veritas.sao.arizona.edu
Submitted: 2011-08-18
We report the detection of pulsed gamma rays from the Crab pulsar at energies above 100 Gigaelectronvolts (GeV) with the VERITAS array of atmospheric Cherenkov telescopes. The detection cannot be explained on the basis of current pulsar models. The photon spectrum of pulsed emission between 100 Megaelectronvolts (MeV) and 400 GeV is described by a broken power law that is statistically preferred over a power law with an exponential cutoff. It is unlikely that the observation can be explained by invoking curvature radiation as the origin of the observed gamma rays above 100 GeV. Our findings require that these gamma rays be produced more than 10 stellar radii from the neutron star.
[43]  oai:arXiv.org:1107.1738  [pdf] - 1077861
VERITAS Observations of the Unusual Extragalactic Transient Swift J164449.3+573451
Comments:
Submitted: 2011-07-08
We report on very-high-energy ($>$100 GeV) gamma-ray observations of Swift J164449.3+573451, an unusual transient object first detected by the {\it Swift} Observatory and later detected by multiple radio, optical and X-ray observatories. A total exposure of 28 hours was obtained on Swift J164449.3+573451 with VERITAS during 2011 March 28 -- April 15. We do not detect the source and place a differential upper limit on the emission at 500 GeV during these observations of $1.4 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ (99% confidence level). We also present time-resolved upper limits and use a flux limit averaged over the X-ray flaring period to constrain various emission scenarios that can accommodate both the radio-through-X-ray emission detected from the source and the lack of detection by VERITAS.
[44]  oai:arXiv.org:1105.0449  [pdf] - 1076364
VERITAS Observations of the TeV Binary LS I +61 303 During 2008-2010
Comments: Accepted for Publication in The Astrophysical Journal, 21 pages, 3 tables, 6 figures
Submitted: 2011-05-02
We present the results of observations of the TeV binary LS I +61 303 with the VERITAS telescope array between 2008 and 2010, at energies above 300 GeV. In the past, both ground-based gamma-ray telescopes VERITAS and MAGIC have reported detections of TeV emission near the apastron phases of the binary orbit. The observations presented here show no strong evidence for TeV emission during these orbital phases; however, during observations taken in late 2010, significant emission was detected from the source close to the phase of superior conjunction (much closer to periastron passage) at a 5.6 standard deviation (5.6 sigma) post-trials significance. In total, between October 2008 and December 2010 a total exposure of 64.5 hours was accumulated with VERITAS on LS I +61 303, resulting in an excess at the 3.3 sigma significance level for constant emission over the entire integrated dataset. The flux upper limits derived for emission during the previously reliably active TeV phases (i.e. close to apastron) are less than 5% of the Crab Nebula flux in the same energy range. This result stands in apparent contrast to previous observations by both MAGIC and VERITAS which detected the source during these phases at >10% of the Crab Nebula flux. During the two year span of observations, a large amount of X-ray data were also accrued on LS I +61 303 by the Swift X-ray Telescope (XRT) and the Rossi X-ray Timing Explorer Timing (RXTE) Proportional Counter Array (PCA). We find no evidence for a correlation between emission in the X-ray and TeV regimes during 20 directly overlapping observations. We also comment on data obtained contemporaneously by the Fermi Large Area Telescope (LAT).
[45]  oai:arXiv.org:1102.3871  [pdf] - 966390
Discovery of TeV Gamma Ray Emission from Tycho's Supernova Remnant
Comments: 6 pages, 3 figures. Accepted in ApJL
Submitted: 2011-02-18, last modified: 2011-04-18
We report the discovery of TeV gamma-ray emission from the Type Ia supernova remnant (SNR) G120.1+1.4, known as Tycho's supernova remnant. Observations performed in the period 2008-2010 with the VERITAS ground-based gamma-ray observatory reveal weak emission coming from the direction of the remnant, compatible with a point source located at $00^{\rm h} \ 25^{\rm m} \ 27.0^{\rm s},\ +64^{\circ} \ 10^{\prime} \ 50^{\prime\prime}$ (J2000). The TeV photon spectrum measured by VERITAS can be described with a power-law $dN/dE = C(E/3.42\;\textrm{TeV})^{-\Gamma}$ with $\Gamma = 1.95 \pm 0.51_{stat} \pm 0.30_{sys}$ and $C = (1.55 \pm 0.43_{stat} \pm 0.47_{sys}) \times 10^{-14}$ cm$^{-2}$s$^{-1}$TeV$^{-1}$. The integral flux above 1 TeV corresponds to $\sim 0.9%$ percent of the steady Crab Nebula emission above the same energy, making it one of the weakest sources yet detected in TeV gamma rays. We present both leptonic and hadronic models which can describe the data. The lowest magnetic field allowed in these models is $\sim 80 \mu$G, which may be interpreted as evidence for magnetic field amplification.
[46]  oai:arXiv.org:1101.3498  [pdf] - 956131
Strong New Constraints on the Extragalactic Background Light in the Near- to Mid-IR
Comments: 37 pages, 13 figures, updated subject headings, accepted by ApJ
Submitted: 2011-01-18, last modified: 2011-03-29
Direct measurements of the extragalactic background light (EBL) in the near-IR to mid-IR waveband are extremely difficult due to an overwhelming foreground from the zodiacal light that outshines the faint cosmological diffuse radiation field by more than an order of magnitude. Indirect constraints on the EBL are provided by gamma-ray observations of AGN. Using the combination of the Fermi Gamma-Ray Space Telescope together with the current generation of ground-based air Cherenkov telescopes (H.E.S.S., MAGIC, and VERITAS) provides unprecedented sensitivity and spectral coverage for constraining the EBL in the near- to mid-IR. In this paper we present new limits on the EBL based on the analysis of the broad-band spectra of a select set of gamma-ray blazars covering 200 MeV to several TeV. The EBL intensity at 15 microns is constrained to be 1.36 +/- 0.58 nW m^-2 sr^-1. We find that the fast evolution and baseline EBL models of Stecker et al. (2006), as well as the model of Kneiske et al. (2004), predict significantly higher EBL intensities in the mid-IR (15 microns) than is allowed by the constraints derived here. In addition, the model of Franceschini et al. (2008) and the fiducial model of Dominguez et al. (2011) predict near- to mid-IR ratios smaller than that predicted by our analysis. Namely, their intensities in the near-IR are too low while their intensities in the mid-IR are marginally too high. All of the aforementioned models are inconsistent with our analysis at the >3 sigma level.