Abstract We have searched for radio pulsations toward 49 Fermi Large Area Telescope (LAT) 1FGL Catalog γ -ray sources using the Green Bank Telescope at 350 MHz. We detected 18 millisecond pulsars ...(MSPs) in blind searches of the data; 10 of these were discoveries unique to our survey. 16 are binaries, with eight having short orbital periods P B < 1 day. No radio pulsations from young pulsars were detected, although three targets are coincident with apparently radio-quiet γ -ray pulsars discovered in LAT data. Here, we give an overview of the survey and present radio and γ -ray timing results for the 10 MSPs discovered. These include the only isolated MSP discovered in our survey and six short- P B binary MSPs. Of these, three have very-low-mass companions ( M c ≪ 0.1 M ⊙ ) and hence belong to the class of black widow pulsars. Two have more massive, nondegenerate companions with extensive radio eclipses and orbitally modulated X-ray emission consistent with the redback class. Significant γ -ray pulsations have been detected from nine of the discoveries. This survey and similar efforts suggest that the majority of Galactic γ -ray sources at high Galactic latitudes are either MSPs or relatively nearby nonrecycled pulsars, with the latter having on average a much smaller radio/ γ -ray beaming ratio as compared to MSPs. It also confirms that past surveys suffered from an observational bias against finding short- P B MSP systems.
Aims. We investigate the extension of the very high-energy spectral tail of the Crab Pulsar at energies above 400 GeV. Methods. We analyzed ~320 h of good-quality Crab data obtained with the MAGIC ...telescope from February 2007 to April 2014. Results. We report the most energetic pulsed emission ever detected from the Crab Pulsar reaching up to 1.5 TeV. The pulse profile shows two narrow peaks synchronized with those measured in the GeV energy range. The spectra of the two peaks follow two different power-law functions from 70 GeV up to 1.5 TeV and connect smoothly with the spectra measured above 10 GeV by the Large Area Telescope (LAT) on board the Fermi satellite. When making a joint fit of the LAT and MAGIC data above 10 GeV the photon indices of the spectra differ by 0.5 ± 0.1. Conclusions. Using data from the MAGIC telescopes we measured the most energetic pulsed photons from a pulsar to date. Such TeV pulsed photons require a parent population of electrons with a Lorentz factor of at least 5 × 106. These results strongly suggest IC scattering off low-energy photons as the emission mechanism and a gamma-ray production region in the vicinity of the light cylinder.
Context. QSO B0218+357 is a gravitationally lensed blazar located at a redshift of 0.944. The gravitational lensing splits the emitted radiation into two components that are spatially ...indistinguishable by gamma-ray instruments, but separated by a 10–12 day delay. In July 2014, QSO B0218+357 experienced a violent flare observed by the Fermi-LAT and followed by the MAGIC telescopes. Aims. The spectral energy distribution of QSO B0218+357 can give information on the energetics of z ~ 1 very high energy gamma-ray sources. Moreover the gamma-ray emission can also be used as a probe of the extragalactic background light at z ~ 1. Methods. MAGIC performed observations of QSO B0218+357 during the expected arrival time of the delayed component of the emission. The MAGIC and Fermi-LAT observations were accompanied by quasi-simultaneous optical data from the KVA telescope and X-ray observations by Swift-XRT. We construct a multiwavelength spectral energy distribution of QSO B0218+357 and use it to model the source. The GeV and sub-TeV data obtained by Fermi-LAT and MAGIC are used to set constraints on the extragalactic background light. Results. Very high energy gamma-ray emission was detected from the direction of QSO B0218+357 by the MAGIC telescopes during the expected time of arrival of the trailing component of the flare, making it the farthest very high energy gamma-ray source detected to date. The observed emission spans the energy range from 65 to 175 GeV. The combined MAGIC and Fermi-LAT spectral energy distribution of QSO B0218+357 is consistent with current extragalactic background light models. The broadband emission can be modeled in the framework of a two-zone external Compton scenario, where the GeV emission comes from an emission region in the jet, located outside the broad line region.
Abstract
Scattered sunlight from the interplanetary dust (IPD) cloud in our solar system presents a serious foreground challenge for spectrophotometric measurements of the extragalactic background ...light (EBL). In this work, we report on inferred measurements of the absolute intensity of the zodiacal light (ZL) using the novel technique of Fraunhofer line spectroscopy on the deepest 8542 Å line of the near-infrared Ca
ii
absorption triplet. The measurements are performed with the narrow band spectrometer (NBS) on board the Cosmic Infrared Background Experiment sounding rocket instrument. We use the NBS data to test the accuracy of two ZL models widely cited in the literature, the Kelsall and Wright models, which have been used in foreground removal analyses that produce high and low EBL results respectively. We find a mean reduced
χ
2
= 3.5 for the Kelsall model and
χ
2
= 2.0 for the Wright model. The best description of our data is provided by a simple modification to the Kelsall model, which includes a free ZL offset parameter. This adjusted model describes the data with a reduced
χ
2
= 1.5 and yields an inferred offset amplitude of 46 ± 19 nW m
−2
sr
−1
extrapolated to 12500 Å. These measurements elude to the potential existence of a dust cloud component in the inner solar system whose intensity does not strongly modulate with the Earth’s motion around the Sun.
Abstract
Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their ...conventional parameter space (100 GeV ≲
M
χ
≲ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM;
M
χ
≳ 100 TeV) has been suggested as an underexplored alternative to the WIMP paradigm, we search for an indirect dark matter annihilation signal in a higher mass range (up to 30 PeV) with the VERITAS
γ
-ray observatory. With 216 hr of observations of four dwarf spheroidal galaxies, we perform an unbinned likelihood analysis. We find no evidence of a
γ
-ray signal from UHDM annihilation above the background fluctuation for any individual dwarf galaxy nor for a joint-fit analysis, and consequently constrain the velocity-weighted annihilation cross section of UHDM for dark matter particle masses between 1 TeV and 30 PeV. We additionally set constraints on the allowed radius of a composite UHDM particle.
Context. During February–March 2014, the MAGIC telescopes observed the high-frequency peaked BL Lac 1ES 1011+496 (z = 0.212) in flaring state at very-high energy (VHE, E> 100 GeV). The flux reached a ...level of more than ten times higher than any previously recorded flaring state of the source. Aims. To describe the characteristics of the flare presenting the light curve and the spectral parameters of the night-wise spectra and the average spectrum of the whole period. From these data we aim to detect the imprint of the extragalactic background light (EBL) in the VHE spectrum of the source, to constrain its intensity in the optical band. Methods. We analyzed the gamma-ray data from the MAGIC telescopes using the standard MAGIC software for the production of the light curve and the spectra. To constrain the EBL, we implement the method developed by the H.E.S.S. collaboration, in which the intrinsic energy spectrum of the source is modeled with a simple function (≤4 parameters), and the EBL-induced optical depth is calculated using a template EBL model. The likelihood of the observed spectrum is then maximized, including a normalization factor for the EBL opacity among the free parameters. Results. The collected data allowed us to describe the night-wise flux changes and also to produce differential energy spectra for all nights in the observed period. The estimated intrinsic spectra of all the nights could be fitted by power-law functions. Evaluating the changes in the fit parameters, we conclude that the spectral shape for most of the nights were compatible, regardless of the flux level, which enabled us to produce an average spectrum from which the EBL imprint could be constrained. The likelihood ratio test shows that the model with an EBL density 1.07 (–0.20, +0.24)stat+sys, relative to the one in the tested EBL template, is preferred at the 4.6σ level to the no-EBL hypothesis, with the assumption that the intrinsic source spectrum can be modeled as a log-parabola. This would translate into a constraint of the EBL density in the wavelength range 0.24 μm, 4.25 μm, with a peak value at 1.4 μm of λFλ = 12.27-2.29+2.75 nW m-2 sr-1, including systematics.
Abstract We use the Very Energetic Radiation Imaging telescope Array System (VERITAS) imaging air Cherenkov telescope array to obtain the first measured angular diameter of β UMa at visual ...wavelengths using stellar intensity interferometry (SII) and independently constrain the limb-darkened angular diameter. The age of the Ursa Major moving group has been assessed from the ages of its members, including nuclear member Merak ( β UMa), an A1-type subgiant, by comparing effective temperature and luminosity constraints to model stellar evolution tracks. Previous interferometric limb-darkened angular-diameter measurements of β UMa in the near-infrared (Center for High Angular Resolution Astronomy (CHARA) Array, 1.149 ± 0.014 mas) and mid-infrared (Keck Nuller, 1.08 ± 0.07 mas), together with the measured parallax and bolometric flux, have constrained the effective temperature. This paper presents current VERITAS-SII observation and analysis procedures to derive squared visibilities from correlation functions. We fit the resulting squared visibilities to find a limb-darkened angular diameter of 1.07 ± 0.04 (stat) ± 0.05 (sys) mas, using synthetic visibilities from a stellar atmosphere model that provides a good match to the spectrum of β UMa in the optical wave band. The VERITAS-SII limb-darkened angular diameter yields an effective temperature of 9700 ± 200 ± 200 K, consistent with ultraviolet spectrophotometry, and an age of 390 ± 29 ± 32 Myr, using MESA Isochrones and Stellar Tracks. This age is consistent with 408 ± 6 Myr from the CHARA Array angular diameter.
Abstract
G106.3+2.7, commonly considered to be a composite supernova remnant (SNR), is characterized by a boomerang-shaped pulsar wind nebula (PWN) and two distinct (“head” and “tail”) regions in the ...radio band. A discovery of very-high-energy gamma-ray emission (
E
γ
> 100 GeV) followed by the recent detection of ultrahigh-energy gamma-ray emission (
E
γ
> 100 TeV) from the tail region suggests that G106.3+2.7 is a PeVatron candidate. We present a comprehensive multiwavelength study of the Boomerang PWN (100″ around PSR J2229+6114) using archival radio and Chandra data obtained two decades ago, a new NuSTAR X-ray observation from 2020, and upper limits on gamma-ray fluxes obtained by Fermi-LAT and VERITAS observatories. The NuSTAR observation allowed us to detect a 51.67 ms spin period from the pulsar PSR J2229+6114 and the PWN emission characterized by a power-law model with Γ = 1.52 ± 0.06 up to 20 keV. Contrary to the previous radio study by Kothes et al., we prefer a much lower PWN
B
-field (
B
∼ 3
μ
G) and larger distance (
d
∼ 8 kpc) based on (1) the nonvarying X-ray flux over the last two decades, (2) the energy-dependent X-ray size of the PWN resulting from synchrotron burn-off, and (3) the multiwavelength spectral energy distribution (SED) data. Our SED model suggests that the PWN is currently re-expanding after being compressed by the SNR reverse shock ∼1000 yr ago. In this case, the head region should be formed by GeV–TeV electrons injected earlier by the pulsar propagating into the low-density environment.
Abstract
Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ∼ 10–100 times greater than ordinary core-collapse supernovae. One popular model to explain the ...enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma-rays, generated in the wind nebula of the magnetar, could escape through the expanding supernova ejecta at late times (months or more after optical peak). This paper presents a search for gamma-ray emission in the broad energy band from 100 MeV to 30 TeV from two Type I SLSNe, SN2015bn, and SN2017egm, using observations from Fermi-LAT and VERITAS. Although no gamma-ray emission was detected from either source, the derived upper limits approach the putative magnetar’s spin-down luminosity. Prospects are explored for detecting very-high-energy (VHE; 100 GeV–100 TeV) emission from SLSNe-I with existing and planned facilities such as VERITAS and CTA.
Abstract
We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS ...observations of the source were triggered by the detection of a hard-spectrum GeV flare on 2015 May 15 with the Fermi-Large Area Telescope (LAT). A combined 5 hr VERITAS exposure on May 16 and 18 resulted in a strong 13
σ
detection with a differential photon spectral index, Γ = 3.8 ± 0.4, and a flux level at 9% of the Crab Nebula above 120 GeV. This also triggered target-of-opportunity observations with Swift, optical photometry, polarimetry, and radio measurements, also presented in this work, in addition to the VERITAS and Fermi-LAT data. A temporal analysis of the gamma-ray flux during this period finds evidence of a shortest variability timescale of
τ
obs
= 6.2 ± 0.9 hr, indicating emission from compact regions within the jet, and the combined gamma-ray spectrum shows no strong evidence of a spectral cutoff. An investigation into correlations between the multiwavelength observations found evidence of optical and gamma-ray correlations, suggesting a single-zone model of emission. Finally, the multiwavelength spectral energy distribution is well described by a simple one-zone leptonic synchrotron self-Compton radiation model.