skip to main content


This content will become publicly available on September 1, 2025

Title: Observation of the Galactic Center PeVatron beyond 100 TeV with HAWC
Abstract

We report an observation of ultrahigh-energy (UHE) gamma rays from the Galactic center (GC) region, using 7 yr of data collected by the High-Altitude Water Cherenkov (HAWC) Observatory. The HAWC data are best described as a point-like source (HAWC J1746-2856) with a power-law spectrum (dN/dE=ϕE/26TeVγ), whereγ= −2.88 ± 0.15stat− 0.1sysandϕ= 1.5 × 10−15(TeV cm2s)−1±0.3stat0.13sys+0.08sysextending from 6 to 114 TeV. We find no evidence of a spectral cutoff up to 100 TeV using HAWC data. Two known point-like gamma-ray sources are spatially coincident with the HAWC gamma-ray excess: Sgr A* (HESS J1745-290) and the Arc (HESS J1746-285). We subtract the known flux contribution of these point sources from the measured flux of HAWC J1746-2856 to exclude their contamination and show that the excess observed by HAWC remains significant (>5σ), with the spectrum extending to >100 TeV. Our result supports that these detected UHE gamma rays can originate via hadronic interaction of PeV cosmic-ray protons with the dense ambient gas and confirms the presence of a proton PeVatron at the GC.

 
more » « less
Award ID(s):
1914549 1806408 1912708 2310104 2209103
PAR ID:
10552328
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Publisher / Repository:
American Astronomical Society
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
973
Issue:
1
ISSN:
2041-8205
Page Range / eLocation ID:
L34
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta (βp= 0.25) collisionless ion–electron shocks with mass ratiomi/me= 200, fast Mach numberMms=1–4, and upstream magnetic field angleθBn= 55°–85° from the shock normalnˆ. It is known that shock electron heating is described by an ambipolar,B-parallel electric potential jump, Δϕ, that scales roughly linearly with the electron temperature jump. Our simulations haveΔϕ/(0.5miush2)0.1–0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measureϕ, including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate ofϕin our low-βpshocks. We further focus on twoθBn= 65° shocks: aMs=4(MA=1.8) case with a long, 30diprecursor of whistler waves alongnˆ, and aMs=7(MA=3.2) case with a shorter, 5diprecursor of whistlers oblique to bothnˆandB;diis the ion skin depth. Within the precursors,ϕhas a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of theMs=4,θBn= 65° case,ϕshows a weak dependence on the electron plasma-to-cyclotron frequency ratioωpece, andϕdecreases by a factor of 2 asmi/meis raised to the true proton–electron value of 1836.

     
    more » « less
  2. Abstract

    We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over 3 yr (from 2018 October to 2021 September). We infer the black hole (BH) mass (MBH) with host galaxy scaling relations, showing that the sampleMBHranges from 105.1Mto 108.2M. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frameg-band luminosity function can be well described by a broken power law ofϕ(Lg)Lg/Lbk0.3+Lg/Lbk2.61, withLbk= 1043.1erg s−1. In the BH mass regime of 105.3≲ (MBH/M) ≲ 107.3, the TDE mass function followsϕ(MBH)MBH0.25, which favors a flat local BH mass function (dnBH/dlogMBHconstant). We confirm the significant rate suppression at the high-mass end (MBH≳ 107.5M), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass ofMgal∼ 1010M, the average optical TDE rate is ≈3.2 × 10−5galaxy−1yr−1. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6] × 10−5galaxy−1yr−1in galaxies with red, green, and blue colors.

     
    more » « less
  3. A<sc>bstract</sc>

    A search for the fully reconstructed$$ {B}_s^0 $$Bs0→ μ+μγdecay is performed at the LHCb experiment using proton-proton collisions at$$ \sqrt{s} $$s= 13 TeV corresponding to an integrated luminosity of 5.4 fb1. No significant signal is found and upper limits on the branching fraction in intervals of the dimuon mass are set$$ {\displaystyle \begin{array}{cc}\mathcal{B}\left({B}_s^0\to {\mu}^{+}{\mu}^{-}\gamma \right)<4.2\times {10}^{-8},& m\left({\mu}^{+}{\mu}^{-}\right)\in \left[2{m}_{\mu },1.70\right]\textrm{GeV}/{c}^2,\\ {}\mathcal{B}\left({B}_s^0\to {\mu}^{+}{\mu}^{-}\gamma \right)<7.7\times {10}^{-8},&\ m\left({\mu}^{+}{\mu}^{-}\right)\in \left[\textrm{1.70,2.88}\right]\textrm{GeV}/{c}^2,\\ {}\mathcal{B}\left({B}_s^0\to {\mu}^{+}{\mu}^{-}\gamma \right)<4.2\times {10}^{-8},& m\left({\mu}^{+}{\mu}^{-}\right)\in \left[3.92,{m}_{B_s^0}\right]\textrm{GeV}/{c}^2,\end{array}} $$BBs0μ+μγ<4.2×108,mμ+μ2mμ1.70GeV/c2,BBs0μ+μγ<7.7×108,mμ+μ1.70, 2.88GeV/c2,BBs0μ+μγ<4.2×108,mμ+μ3.92mBs0GeV/c2,

    at 95% confidence level. Additionally, upper limits are set on the branching fraction in the [2mμ,1.70] GeV/c2dimuon mass region excluding the contribution from the intermediateϕ(1020) meson, and in the region combining all dimuon-mass intervals.

     
    more » « less
  4. Abstract

    A star completely destroyed in a tidal disruption event (TDE) ignites a luminous flare that is powered by the fallback of tidally stripped debris to a supermassive black hole (SMBH) of massM. We analyze two estimates for the peak fallback rate in a TDE, one being the “frozen-in” model, which predicts a strong dependence of the time to peak fallback rate,tpeak, on both stellar mass and age, with 15 days ≲tpeak≲ 10 yr for main sequence stars with masses 0.2 ≤M/M≤ 5 andM= 106M. The second estimate, which postulates that the star is completely destroyed when tides dominate the maximum stellar self-gravity, predicts thattpeakis very weakly dependent on stellar type, withtpeak=23.2±4.0daysM/106M1/2for 0.2 ≤M/M≤ 5, whiletpeak=29.8±3.6daysM/106M1/2for a Kroupa initial mass function truncated at 1.5M. This second estimate also agrees closely with hydrodynamical simulations, while the frozen-in model is discrepant by orders of magnitude. We conclude that (1) the time to peak luminosity in complete TDEs is almost exclusively determined by SMBH mass, and (2) massive-star TDEs power the largest accretion luminosities. Consequently, (a) decades-long extra-galactic outbursts cannot be powered by complete TDEs, including massive-star disruptions, and (b) the most highly super-Eddington TDEs are powered by the complete disruption of massive stars, which—if responsible for producing jetted TDEs—would explain the rarity of jetted TDEs and their preference for young and star-forming host galaxies.

     
    more » « less
  5. Abstract

    We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizationsσ≲ 10−4and pair-loading factorsZ±≲ 10 are studied, whereZ±is the number of loaded electrons and positrons per ion. We find the following: (1) The shock becomes mediated by the ion Larmor gyration in the mean field whenσexceeds a critical valueσLthat decreases withZ±. AtσσLthe shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease withZ±, leading to lowerσL. (2) The energy fraction carried by the post-shock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per post-shock electron scales asE¯eZ±+11. (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low asσ≈ 5 × 10−6. The ions then become essentially thermal with mean energyE¯i, while electrons form a nonthermal tail, extending fromEZ±+11E¯itoE¯i. Whenσ= 0, particle acceleration is enhanced by the formation of intense magnetic cavities that populate the precursor during the late stages of shock evolution. Here, the maximum energy of the nonthermal ions and electrons keeps growing over the duration of the simulation. Alongside the simulations, we develop theoretical estimates consistent with the numerical results. Our findings have important implications for models of early gamma-ray burst afterglows.

     
    more » « less