366 pubs · 15,117 cites · h-index 69
Cherenkov telescope Trinity One detects neutrinos from 1 PeV to 10 EeV in Utah, highlighting its sensitivity and design.
Investigates intrinsic characteristics of AlGaN-based avalanche photodiodes, suggesting design enhancements.
We developed a modular silicon photomultiplier camera to detect Earth-skimming PeV to EeV tau neutrinos with the imaging atmospheric Cherenkov technique. We built two cameras, a 256-pixel camera with S14161-6050HS SiPMs for the Trinity Demonstrator located on Frisco Peak, Utah, and a 512-pixel camera with S14521-6050AN SiPMs for the EUSO-SPB2 Cherenkov Telescope. The front-end electronics are based on the eMUSIC ASIC, and the camera signals are sampled and digitized with the 100MS/s and 12-bit AGET system. Both cameras are liquid-cooled. We detail the camera concept and the results from characterizing the SiPMs, bench testing, and calibrating the two cameras.
This work reports high-performance homojunction gallium nitride (GaN) p-i-n (PIN) avalanche photodiodes (APDs) grown on a low-defect GaN substrate and fabricated with an <tex-math notation="LaTeX">~</tex-math> 1.8 <tex-math notation="LaTeX">∘</tex-math> shallow-bevel-mesa edge termination. High-quality p-i-n GaN epitaxial layers were grown using metalorganic chemical vapor deposition (MOCVD), and the device edge termination was achieved by multistep pseudo gray-scale lithography followed by photoresist reflow and dry etching of the bevel mesa. Fabricated GaN APDs showed an ultralow dark current density <tex-math notation="LaTeX">$<$</tex-math> 10 <tex-math notation="LaTeX">⁻⁹</tex-math> A/cm <tex-math notation="LaTeX">²</tex-math> in the low-mid reverse bias voltage region prior to avalanche breakdown and achieved a current-limited photocurrent gain of <tex-math notation="LaTeX">$>$</tex-math> 10 <tex-math notation="LaTeX">⁷</tex-math> under deep-ultraviolet (DUV) illumination at <tex-math notation="LaTeX">λ</tex-math> <tex-math notation="LaTeX">$=$</tex-math> 280 nm and a reverse bias <tex-math notation="LaTeX">$-$</tex-math> 107 V. The total current density without illumination was in the order of 10 <tex-math notation="LaTeX">⁻⁵</tex-math> A/cm <tex-math notation="LaTeX">²</tex-math> at the onset of avalanche breakdown. The avalanche breakdown characteristic was confirmed by temperature-dependent J – V studies, with a derived positive temperature coefficient of 0.02 V/K for the breakdown voltage (BV). The analysis of the temperature-dependent reverse-bias characteristics indicated that the reverse leakage current components before avalanche could be attributed to trap-assisted tunneling (TAT) and phonon-assisted tunneling (PAT) processes.
This is a collection of papers presented by the JEM-EUSO Collaboration at the 38th International Cosmic Ray Conference (Nagoya, Japan, July 26-August 3, 2023)
A deep-ultraviolet Al0.6Ga0.4N p–i–n avalanche photodiode (APD) structure was grown on a (0001) AlN bulk substrate by metalorganic chemical vapor deposition. The wafer was fabricated into 20 μm diameter mesa APD devices both with and without ion-implantation with nitrogen ions on the periphery of the p-type region of the diode mesa and tested. The dark current density vs bias, photoresponse, and the optical gain of the APDs with and without ion implantation were compared. The devices fabricated with ion implantation showed improved performance, exhibiting lower dark current densities of ∼1 × 10−9 A/cm2 and a higher optical gain of ∼5.2 × 105 at a current density limit of 0.3 A/cm2. The average temperature coefficients of the reverse-bias breakdown voltage were also compared. Although the data showed negative coefficients for APDs fabricated both with and without ion implantation, the ion-implanted APDs showed an improvement relative to the devices fabricated without ion-implantation.
We present the low-temperature Geiger-mode characteristics of GaN (gallium nitride) p-i-n avalanche photodiodes (APDs). The devices have a breakdown voltage of −95V and a temperature dependence of~0.0159 ± 0.0034 V/K near 300K. The room-temperature (300 K) dark-count rate (DCR) is 23.8 MHz for a75× 75\,\,μ m²device biased at 1 % overvoltage. The DCR halves when lowering the temperature by 50°C. Based on the temperature-dependent characteristics of the DCR, we identify band-to-band tunneling as the dominant DCR generation mechanism. At 4.65 % overvoltage and 375 nm, the photon detection efficiency (PDE) is 0.82 %- limited by a low breakdown probability of 1.7 %. We discuss the measurement setup and the method to extract count rates, which is based on the Poisson distribution of the time intervals between Geiger-mode breakdowns of the APD. The setup includes a custom circuit to bias the diode and amplify its signals, a steady ultraviolet (UV) light source, and a system to control the temperature of the APD with a thermoelectric element in the range from −40 to 20°C.
Review outlines future tau neutrino detection across accelerator and astrophysical frontiers, highlighting opportunities to test standard model physics.
We report on a long-lasting, elevated gamma-ray flux state from VER J0521+211 observed by VERITAS, MAGIC, and Fermi-LAT in 2013 and 2014. The peak integral flux above 200 GeV measured with the nightly binned light curve is (8.8 ± 0.4) × 10 −7 photons m −2 s −1 , or ∼37% of the Crab Nebula flux. Multiwavelength observations from X-ray, UV, and optical instruments are also presented. A moderate correlation between the X-ray and TeV gamma-ray fluxes was observed, and the X-ray spectrum appeared harder when the flux was higher. Using the gamma-ray spectrum and four models of the extragalactic background light (EBL), a conservative 95% confidence upper limit on the redshift of the source was found to be z ≤ 0.31. Unlike the gamma-ray and X-ray bands, the optical flux did not increase significantly during the studied period compared to the archival low-state flux. The spectral variability from optical to X-ray bands suggests that the synchrotron peak of the spectral energy distribution (SED) may become broader during flaring states, which can be adequately described with a one-zone synchrotron self-Compton model varying the high-energy end of the underlying particle spectrum. The synchrotron peak frequency of the SED and the radio morphology of the jet from the MOJAVE program are consistent with the source being an intermediate-frequency-peaked BL Lac object.
Gamma-ray astronomy has become one of the main experimental ways to test the modified dispersion relations (MDRs) of photons in vacuum, obtained in some attempts to formulate a theory of quantum gravity. The MDRs in use imply time delays that depend on the energy and that increase with distance following some function of redshift. The use of transient, or variable, distant and highly energetic sources already allows us to set stringent limits on the energy scale related to this phenomenon, usually thought to be of the order of the Planck energy, but robust conclusions on the existence of MDR-related propagation effects still require the analysis of a large population of sources. In order to gather the biggest sample of sources possible for MDR searches at teraelectronvolt energies, the H.E.S.S., MAGIC, and VERITAS collaborations enacted a joint task force to combine all their relevant data to constrain the quantum gravity energy scale. In the present article, the likelihood method used to combine the data and provide a common limit is described in detail and tested through simulations of recorded data sets for a gamma-ray burst, three flaring active galactic nuclei, and two pulsars. Statistical and systematic errors are assessed and included in the likelihood as nuisance parameters. In addition, a comparison of two different formalisms for distance dependence of the time lags is performed for the first time. In a second article, to appear later, the method will be applied to all relevant data from the three experiments.
A top-illuminated deep-ultraviolet Al0.6Ga0.4N p-i-n avalanche photodiode (APD) structure was designed and grown by metalorganic chemical vapor deposition on an AlN bulk substrate and on two different quality AlN/sapphire templates, and APDs were fabricated and tested. The APD devices with a circular diameter of 20 μm have demonstrated a distinctive reverse-bias breakdown behavior. The reverse breakdown voltage of the APDs is approximately −140 V, which corresponds to a breakdown electric field of 6–6.2 MV/cm for the Al0.6Ga0.4N material as estimated by Silvaco TCAD simulation. The APDs grown on the AlN bulk substrate show the lowest leakage current density of <1 × 10−8 A/cm2 (at low reverse bias) compared to that of the devices grown on the AlN templates. From the photocurrent measurement, a maximum gain (current limited) of 1.2 × 104 is calculated. The average temperature coefficients of the breakdown voltage are negative for APD devices fabricated from both the AlN bulk substrate and the AlN templates, but these data show that the coefficient is the least negative for the APD devices grown on the low-dislocation-density AlN bulk substrate.
Flat-spectrum radio quasars (FSRQs) are the most luminous blazars at GeV energies but only rarely emit detectable fluxes of TeV gamma rays, typically during bright GeV flares. We explore the gamma-ray variability and spectral characteristics of three FSRQs that have been observed at GeV and TeV energies by Fermi-LAT and VERITAS, making use of almost 100 hr of VERITAS observations spread over 10 yr: 3C 279, PKS 1222+216, and Ton 599. We explain the GeV flux distributions of the sources in terms of a model derived from a stochastic differential equation describing fluctuations in the magnetic field in the accretion disk and estimate the timescales of magnetic flux accumulation and stochastic instabilities in their accretion disks. We identify distinct flares using a procedure based on Bayesian blocks and analyze their daily and subdaily variability and gamma-ray energy spectra. Using observations from VERITAS, as well as Fermi, Swift, and the Steward Observatory, we model the broadband spectral energy distributions of PKS 1222+216 and Ton 599 during very high energy (VHE)–detected flares in 2014 and 2017, respectively, strongly constraining the jet Doppler factors and gamma-ray emission region locations during these events. Finally, we place theoretical constraints on the potential production of PeV-scale neutrinos during these VHE flares.
HESS J0632+057 belongs to a rare subclass of binary systems that emit gamma rays above 100 GeV. It stands out for its distinctive high-energy light curve, which features a sharp “primary” peak and broader “secondary” peak. We present the results of contemporaneous observations by NuSTAR and VERITAS during the secondary peak between 2019 December and 2020 February, when the orbital phase ( ϕ ) is between 0.55 and 0.75. NuSTAR detected X-ray spectral evolution, while VERITAS detected TeV emission. We fit a leptonic wind-collision model to the multiwavelength spectra data obtained over the four NuSTAR and VERITAS observations, constraining the pulsar spin-down luminosity and the magnetization parameter at the shock. Despite long-term monitoring of the source from 2019 October to 2020 March, the MDM observatory did not detect significant variation in H α and H β line equivalent widths, an expected signature of Be-disk interaction with the pulsar. Furthermore, fitting folded Swift-XRT light-curve data with an intrabinary shock model constrained the orbital parameters, suggesting two orbital phases (at ϕ D = 0.13 and 0.37), where the pulsar crosses the Be-disk, as well as phases for the periastron ( ϕ 0 = 0.30) and inferior conjunction ( ϕ IFC = 0.75). The broadband X-ray spectra with Swift-XRT and NuSTAR allowed us to measure a higher neutral hydrogen column density at one of the predicted disk-passing phases.
The results of gamma-ray observations of the binary system HESS J0632 + 057 collected during 450 hr over 15 yr, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the Hα emission line. A significant detection of the modulation of the very high-energy gamma-ray fluxes with a period of 316.7 ± 4.4 days is reported, consistent with the period of 317.3 ± 0.7 days obtained with a refined analysis of X-ray data. The analysis of data from four orbital cycles with dense observational coverage reveals short-timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over a timescale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance, but cannot find any correlation of optical Hα parameters with fluxes at X-ray or gamma-ray energies in simultaneous observations. The key finding is that the emission of HESS J0632 + 057 in the X-ray and gamma-ray energy bands is highly variable on different timescales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system, and discussed in the context of results reported on other gamma-ray binary systems.
The Galactic Center (GC) region hosts a variety of powerful astronomical sources and rare astrophysical processes that emit a large flux of nonthermal radiation. The inner 375 pc × 600 pc region, called the Central Molecular Zone, is home to the supermassive black hole Sagittarius A*, massive cloud complexes, and particle accelerators such as supernova remnants (SNRs). We present the results of our improved analysis of the very-high-energy gamma-ray emission above 2 TeV from the GC using 125 hr of data taken with the Very Energetic Radiation Imaging Telescope Array System imaging-atmospheric Cerenkov telescope between 2010 and 2018. The central source VER J1745–290, consistent with the position of Sagittarius A*, is detected at a significance of 38 standard deviations above the background level (38 σ ), and we report its spectrum and light curve. Its differential spectrum is consistent with a power law with exponential cutoff, with a spectral index of , a flux normalization at 5.3 TeV of TeV −1 cm −2 s −1 , and cutoff energy of TeV. We also present results on the diffuse emission near the GC, obtained by combining data from multiple regions along the GC ridge, which yield a cumulative significance of 9.5 σ . The diffuse GC ridge spectrum is best fit by a power law with a hard index of 2.19 ± 0.20, showing no evidence of a cutoff up to 40 TeV. This strengthens the evidence for a potential accelerator of PeV cosmic rays being present in the GC. We also provide spectra of the other sources in our field of view with significant detections, composite SNR G0.9+0.1, and HESS J1746–285.
The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to accurately observe ultra-high-energy cosmic rays (UHECRs) and cosmic neutrinos from space with sensitivity over the full celestial sky. POEMMA will observe the air fluorescence produced by extensive air showers (EASs) from UHECRs and potentially UHE neutrinos above 20 EeV. Additionally, POEMMA has the ability to observe the Cherenkov signal from upward-moving EASs induced by Earth-interacting tau neutrinos above 20 PeV. The POEMMA spacecraft are designed to quickly re-orientate to follow up transient neutrino sources and obtain currently unparalleled neutrino flux sensitivity. Developed as a NASA Astrophysics Probe-class mission, POEMMA consists of two identical satellites flying in loose formation in 525 km altitude orbits. Each POEMMA instrument incorporates a wide field-of-view (45 ∘ ) Schmidt telescope with an optical collecting area of over 6 m 2 . The hybrid focal surface of each telescope includes a fast (1 μs) near-ultraviolet camera for EAS fluorescence observations and an ultrafast (10 ns) optical camera for Cherenkov EAS observations. In a 5-year mission, POEMMA will provide measurements that open new multi-messenger windows onto the most energetic events in the universe, enabling the study of new astrophysics and particle physics at these extreme energies.
We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of 2010 February, when an extraordinary flare reaching a level of ∼27 Crab Units above 1 TeV was measured in very high energy (VHE) γ-rays with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) observatory. This is the highest flux state for Mrk 421 ever observed in VHE γ-rays. Data are analyzed from a coordinated campaign across multiple instruments, including VHE γ-ray (VERITAS, Major Atmospheric Gamma-ray Imaging Cherenkov), high-energy γ-ray (Fermi-LAT), X-ray (Swift, Rossi X-ray Timing Experiment, MAXI), optical (including the GASP-WEBT collaboration and polarization data), and radio (Metsähovi, Owens Valley Radio Observatory, University of Michigan Radio Astronomy Observatory). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare “decline” epochs. The main flare statistics allow 2 minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of ∼25–55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor (δ ≳ 33) and the size of the emission region ( ) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10 minute binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux–flux relationship, from linear to quadratic to lack of correlation to anticorrelation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain with the classic single-zone synchrotron self-Compton model.
HESS J1943+213 is a very high energy (VHE; >100 GeV) γ-ray source in the direction of the Galactic plane. Studies exploring the classification of the source are converging toward its identification as an extreme synchrotron BL Lac object. Here we present 38 hr of VERITAS observations of HESS J1943+213 taken over 2 yr. The source is detected with a significance of ∼20 standard deviations, showing a remarkably stable flux and spectrum in VHE γ-rays. Multifrequency Very Long Baseline Array (VLBA) observations of the source confirm the extended, jet-like structure previously found in the 1.6 GHz band with the European VLBI Network and detect this component in the 4.6 and 7.3 GHz bands. The radio spectral indices of the core and the jet and the level of polarization derived from the VLBA observations are in a range typical for blazars. Data from VERITAS, Fermi-LAT, Swift-XRT, the FLWO 48″ telescope, and archival infrared and hard X-ray observations are used to construct and model the spectral energy distribution (SED) of the source with a synchrotron self-Compton model. The well-measured γ-ray peak of the SED with VERITAS and Fermi-LAT provides constraining upper limits on the source redshift. Possible contribution of secondary γ-rays from ultra-high-energy cosmic-ray-initiated electromagnetic cascades to the γ-ray emission is explored, finding that only a segment of the VHE spectrum can be accommodated with this process. A variability search is performed across X-ray and γ-ray bands. No statistically significant flux or spectral variability is detected.
On 2017 September 22, the IceCube Neutrino Observatory reported the detection of the high-energy neutrino event IC 170922A, of potential astrophysical origin. It was soon determined that the neutrino direction was consistent with the location of the gamma-ray blazar TXS 0506+056 (3FGL J0509.4+0541), which was in an elevated gamma-ray emission state as measured by the Fermi satellite. Very Energetic Radiation Imaging Telescope Array System (VERITAS) observations of the neutrino/blazar region started on 2017 September 23 in response to the neutrino alert and continued through 2018 February 6. While no significant very-high-energy (VHE; E > 100 GeV) emission was observed from the blazar by VERITAS in the two-week period immediately following the IceCube alert, TXS 0506+056 was detected by VERITAS with a significance of 5.8 standard deviations ( σ ) in the full 35 hr data set. The average photon flux of the source during this period was (8.9 ± 1.6) × 10 −12 cm −2 s −1 , or 1.6% of the Crab Nebula flux, above an energy threshold of 110 GeV, with a soft spectral index of 4.8 ± 1.3.
Combined with measurements made by very-long-baseline interferometry, the observations of fast TeV gamma-ray flares probe the structure and emission mechanism of blazar jets. However, only a handful of such flares have been detected to date, and only within the last few years have these flares been observed from lower-frequency-peaked BL Lac objects and flat-spectrum radio quasars. We report on a fast TeV gamma-ray flare from the blazar BL Lacertae observed by the Very Energetic Radiation Imaging Telescope Array System (VERITAS). with a rise time of ∼2.3 hr and a decay time of ∼36 min. The peak flux above 200 GeV is (4.2 ± 0.6) × 10 −6 photon m −2 s −1 measured with a 4-minute-binned light curve, corresponding to ∼180% of the flux that is observed from the Crab Nebula above the same energy threshold. Variability contemporaneous with the TeV gamma-ray flare was observed in GeV gamma-ray, X-ray, and optical flux, as well as in optical and radio polarization. Additionally, a possible moving emission feature with superluminal apparent velocity was identified in Very Long Baseline Array observations at 43 GHz, potentially passing the radio core of the jet around the time of the gamma-ray flare. We discuss the constraints on the size, Lorentz factor, and location of the emitting region of the flare, and the interpretations with several theoretical models that invoke relativistic plasma passing stationary shocks.