Extended TeV Halos May Commonly Exist around Middle-Aged Pulsars
Not Available
Engineering topics
Publications and source records attributed to Fan, K. L. (ORCID:0000000282464751).
Not Available
The origin of high-energy galactic cosmic rays is yet to be understood, but some galactic cosmic-ray accelerators can accelerate cosmic rays up to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding material or radiation, resulting in the production of gamma-rays and neutrinos. To optimize for the detection of such associated production of gamma-rays and neutrinos for a given source morphology and spectrum, a multimessenger analysis that combines gamma-rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework with IceCube Maximum Likelihood Analysis software and HAWC Accelerated Likelihood to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 yr of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that, for five of the sources, the gamma-ray emission observed by HAWC cannot be produced purely from hadronic interactions. We report the limit for the fraction of gamma-rays produced by hadronic interactions for these five sources.
Active galactic nuclei (AGN) are prime candidate sources of the high-energy, astrophysical neutrinos detected by IceCube. This is demonstrated by the real-time multimessenger detection of the blazar TXS 0506+056 and the recent evidence of neutrino emission from NGC 1068 from a separate time-averaged study. However, the production mechanism of the astrophysical neutrinos in AGN is not well established, which can be resolved via correlation studies with photon observations. For neutrinos produced due to photohadronic interactions in AGN, in addition to a correlation of neutrinos with high-energy photons, there would also be a correlation of neutrinos with photons emitted at radio wavelengths. In this work, we perform an in-depth stacking study of the correlation between 15 GHz radio observations of AGN reported in the MOJAVE XV catalog, and 10 yr of neutrino data from IceCube. We also use a time-dependent approach, which improves the statistical power of the stacking analysis. No significant correlation was found for both analyses, and upper limits are reported. When compared to the IceCube diffuse flux, at 100 TeV and for a spectral index of 2.5, the upper limits derived are ~3% and ~9% for the time-averaged and time-dependent cases, respectively.
Abstract Neutrino flares in the sky are searched for in data collected by IceCube between 2011 and 2021 May. This data set contains cascade-like events originating from charged-current electron neutrino and tau neutrino interactions and all-flavor neutral-current interactions. IceCube’s previous all-sky searches for neutrino flares used data sets consisting of track-like events originating from charged-current muon neutrino interactions. The cascade data set is statistically independent of the track data sets, and while inferior in angular resolution, the low-background nature makes it competitive and complementary to previous searches. No statistically significant flare of neutrino emission was observed in an all-sky scan. Upper limits are calculated on neutrino flares of varying duration from 1 hr to 100 days. Furthermore, constraints on the contribution of these flares to the diffuse astrophysical neutrino flux are presented, showing that multiple unresolved transient sources may contribute to the diffuse astrophysical neutrino flux.
IceCube alert events are neutrinos with a moderate-to-high probability of having astrophysical origin. In this study, we analyze 11 yr of IceCube data and investigate 122 alert events and a selection of high-energy tracks detected between 2009 and the end of 2021. This high-energy event selection (alert events + high-energy tracks) has an average probability of ≥0.5 of being of astrophysical origin. We search for additional continuous and transient neutrino emission within the high-energy events' error regions. We find no evidence for significant continuous neutrino emission from any of the alert event directions. The only locally significant neutrino emission is the transient emission associated with the blazar TXS 0506+056, with a local significance of 3σ, which confirms previous IceCube studies. When correcting for 122 test positions, the global p-value is 0.156 and compatible with the background hypothesis. We constrain the total continuous flux emitted from all 122 test positions at 100 TeV to be below 1.2 × 10 –15 (TeV cm 2 s) –1 at 90% confidence assuming an E –2 spectrum. This corresponds to 4.5% of IceCube's astrophysical diffuse flux. Overall, we find no indication that alert events in general are linked to lower-energetic continuous or transient neutrino emission.
The LIGO/Virgo collaboration published the catalogs GWTC-1, GWTC-2.1, and GWTC-3 containing candidate gravitational-wave (GW) events detected during its runs O1, O2, and O3. These GW events can be possible sites of neutrino emission. In this paper, we present a search for neutrino counterparts of 90 GW candidates using IceCube DeepCore, the low-energy infill array of the IceCube Neutrino Observatory. The search is conducted using an unbinned maximum likelihood method, within a time window of 1000 s, and uses the spatial and timing information from the GW events. The neutrinos used for the search have energies ranging from a few GeV to several tens of TeV. We do not find any significant emission of neutrinos, and place upper limits on the flux and the isotropic-equivalent energy emitted in low-energy neutrinos. We also conduct a binomial test to search for source populations potentially contributing to neutrino emission. We report a nondetection of a significant neutrino-source population with this test.
We present a catalog of likely astrophysical neutrino track-like events from the IceCube Neutrino Observatory. IceCube began reporting likely astrophysical neutrinos in 2016, and this system was updated in 2019. The catalog presented here includes events that were reported in real time since 2019, as well as events identified in archival data samples starting from 2011. We report 275 neutrino events from two selection channels as the first entries in the catalog, the IceCube Event Catalog of Alert Tracks, which will see ongoing extensions with additional alerts. The Gold and Bronze alert channels respectively provide neutrino candidates with a 50% and 30% probability of being astrophysical, on average assuming an astrophysical neutrino power-law energy spectral index of 2.19. For each neutrino alert, we provide the reconstructed energy, direction, false-alarm rate, probability of being astrophysical in origin, and likelihood contours describing the spatial uncertainty in the alert's reconstructed location. We also investigate a directional correlation of these neutrino events with gamma-ray and X-ray catalogs, including 4FGL, 3HWC, TeVCat, and Swift-BAT.
The IceCube Neutrino Observatory sends realtime neutrino alerts with a high probability of being astrophysical in origin. We present a new method to correlate these events and possible candidate sources using 2089 blazars from the Fermi-LAT 4LAC-DR2 catalog and with 3413 active galactic nuclei (AGNs) from the Radio Fundamental Catalog. No statistically significant neutrino emission was found in any of the catalog searches. The result suggests that a small fraction, <1%, of the studied AGNs emit neutrinos that pass the alert criteria, and is compatible with prior evidence for neutrino emission presented by IceCube and other authors from sources such as TXS 0506 + 056 and PKS 1502 + 106. We also present cross-checks to other analyses that claim a significant correlation using similar data samples.
The understanding of novae, the thermonuclear eruptions on the surfaces of white dwarf stars in binaries, has recently undergone a major paradigm shift. Though the bolometric luminosity of novae was long thought to arise directly from photons supplied by the thermonuclear runaway, recent gigaelectronvolt (GeV) gamma-ray observations have supported the notion that a significant portion of the luminosity could come from radiative shocks. More recently, observations of novae have lent evidence that these shocks are acceleration sites for hadrons for at least some types of novae. In this scenario, a flux of neutrinos may accompany the observed gamma rays. As the gamma rays from most novae have only been observed up to a few GeV, novae have previously not been considered as targets for neutrino telescopes, which are most sensitive at and above teraelectronvolt (TeV) energies. Here, we present the first search for neutrinos from novae with energies between a few GeV and 10 TeV using IceCube-DeepCore, a densely instrumented region of the IceCube Neutrino Observatory with a reduced energy threshold. We search both for a correlation between gamma-ray and neutrino emission as well as between optical and neutrino emission from novae. We find no evidence for neutrino emission from the novae considered in this analysis and set upper limits for all gamma-ray detected novae.
Beginning in 2016, the IceCube Neutrino Observatory has sent out alerts in real time containing the information of high-energy ($E$ ≳ 100 TeV) neutrino candidate events with moderate to high (≳30%) probability of astrophysical origin. In this work, we use a recent catalog of such alert events, which, in addition to events announced in real time, includes events that were identified retroactively and covers the time period of 2011–2020. We also search for additional, lower-energy neutrinos from the arrival directions of these IceCube alerts. We show how performing such an analysis can constrain the contribution of rare populations of cosmic neutrino sources to the diffuse astrophysical neutrino flux. After searching for neutrino emission coincident with these alert events on various timescales, we find no significant evidence of either minute-scale or day-scale transient neutrino emission or of steady neutrino emission in the direction of these alert events. This study also shows how numerous a population of neutrino sources has to be to account for the complete astrophysical neutrino flux. Assuming that sources have the same luminosity, an $E$ -2.5 neutrino spectrum, and number densities that follow star formation rates, the population of sources has to be more numerous than 7 × 10 -9 Mpc -3 . This number changes to 3 × 10 -7 Mpc -3 if number densities instead have no cosmic evolution.
Core-collapse supernovae are a promising potential high-energy neutrino source class. We test for correlation between seven years of IceCube neutrino data and a catalog containing more than 1000 core-collapse supernovae of types IIn and IIP and a sample of stripped-envelope supernovae. We search both for neutrino emission from individual supernovae as well as for combined emission from the whole supernova sample, through a stacking analysis. No significant spatial or temporal correlation of neutrinos with the cataloged supernovae was found. All scenarios were tested against the background expectation and together yield an overall p-value of 93%; therefore, they show consistency with the background only. The derived upper limits on the total energy emitted in neutrinos are 1.7 × 10 48 erg for stripped-envelope supernovae, 2.8 × 10 48 erg for type IIP, and 1.3 × 10 49 erg for type IIn SNe, the latter disfavoring models with optimistic assumptions for neutrino production in interacting supernovae. We conclude that stripped-envelope supernovae and supernovae of type IIn do not contribute more than 14.6% and 33.9%, respectively, to the diffuse neutrino flux in the energy range of about [ 10 3 –10 5 ] GeV, assuming that the neutrino energy spectrum follows a power-law with an index of –2.5. Under the same assumption, we can only constrain the contribution of type IIP SNe to no more than 59.9%. Thus, core-collapse supernovae of types IIn and stripped-envelope supernovae can both be ruled out as the dominant source of the diffuse neutrino flux under the given assumptions.
We present the results of dark matter (DM) searches in a sample of 31 dwarf irregular (dIrr) galaxies within the field of view of the HAWC Observatory. dIrr galaxies are DM-dominated objects in which astrophysical gamma-ray emission is estimated to be negligible with respect to the secondary gamma-ray flux expected by annihilation or decay of weakly interacting massive particles (WIMPs). While we do not see any statistically significant DM signal in dIrr galaxies, we present the exclusion limits (95% C.L.) for annihilation cross section and decay lifetime for WIMP candidates with masses between 1 and 100 TeV. Exclusion limits from dIrr galaxies are relevant and complementary to benchmark dwarf Spheroidal (dSph) galaxies. In fact, dIrr galaxies are targets kinematically different from benchmark dSph, preserving the footprints of different evolution histories. We compare the limits from dIrr galaxies to those from ultrafaint and classical dSph galaxies previously observed with HAWC. We find that the constraints are comparable to the limits from classical dSph galaxies and ~2 orders of magnitude weaker than the ultrafaint dSph limits.
Galactic PeV cosmic-ray accelerators (PeVatrons) are Galactic sources theorized to accelerate cosmic rays up to PeV in energy. The accelerated cosmic rays are expected to interact hadronically with nearby ambient gas or the interstellar medium, resulting in $γ$-rays and neutrinos. Recently, the Large High Altitude Air Shower Observatory (LHAASO) identified 12 $γ$-ray sources with emissions above 100 TeV, making them candidates for PeVatrons. While at these high energies the Klein–Nishina effect exponentially suppresses leptonic emission from Galactic sources, evidence for neutrino emission would unequivocally confirm hadronic acceleration. Here, we present the results of a search for neutrinos from these $γ$-ray sources and stacking searches testing for excess neutrino emission from all 12 sources as well as their subcatalogs of supernova remnants and pulsar wind nebulae with 11 yr of track events from the IceCube Neutrino Observatory. No significant emissions were found. Based on the resulting limits, we place constraints on the fraction of $γ$-ray flux originating from the hadronic processes in the Crab Nebula and LHAASO J2226+6057.
Extended very-high-energy (VHE; 0.1–100 TeV) γ-ray emission has been observed around several middle-aged pulsars and referred to as "TeV halos." Their formation mechanism remains under debate. It is also unknown whether they are ubiquitous or related to a certain subgroup of pulsars. With 2321 days of observation, the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory detected VHE γ-ray emission at the location of the radio-quiet pulsar PSR J0359+5414 with >6σ significance. By performing likelihood tests with different spectral and spatial models and comparing the TeV spectrum with multiwavelength observations of nearby sources, we show that this excess is consistent with a TeV halo associated with PSR J0359+5414, though future observation of HAWC and multiwavelength follow-ups are needed to confirm this nature. This new halo candidate is located in a noncrowded region in the outer galaxy. It shares similar properties to the other halos but its pulsar is younger and radio-quiet. Our observation implies that TeV halos could commonly exist around pulsars and their formation does not depend on the configuration of the pulsar magnetosphere.
This paper reports on the γ -ray properties of the 2018 Galactic nova V392 Per, spanning photon energies ~0.1 GeV–100 TeV by combining observations from the Fermi Gamma-ray Space Telescope and the HAWC Observatory. As one of the most rapidly evolving γ -ray signals yet observed for a nova, GeV γ -rays with a power-law spectrum with an index Γ = 2.0 ± 0.1 were detected over 8 days following V392 Per’s optical maximum. HAWC observations constrain the TeV γ -ray signal during this time and also before and after. We observe no statistically significant evidence of TeV γ -ray emission from V392 Per, but present flux limits. Tests disfavor the extension of the Fermi Large Area Telescope spectrum to energies above 5 TeV by 2 standard deviations (95%) or more. We fit V392 Per’s GeV γ -rays with hadronic acceleration models, incorporating optical observations, and compare the calculations with HAWC limits.
Gamma-ray bursts (GRBs) are considered as promising sources of ultra-high-energy cosmic rays (UHECRs) due to their large power output. Observing a neutrino flux from GRBs would offer evidence that GRBs are hadronic accelerators of UHECRs. Previous IceCube analyses, which primarily focused on neutrinos arriving in temporal coincidence with the prompt gamma-rays, found no significant neutrino excess. The four analyses presented in this paper extend the region of interest to 14 days before and after the prompt phase, including generic extended time windows and targeted precursor searches. GRBs were selected between 2011 May and 2018 October to align with the data set of candidate muon-neutrino events observed by IceCube. No evidence of correlation between neutrino events and GRBs was found in these analyses. Limits are set to constrain the contribution of the cosmic GRB population to the diffuse astrophysical neutrino flux observed by IceCube. Prompt neutrino emission from GRBs is limited to ≲1% of the observed diffuse neutrino flux, and emission on timescales up to 10 4 s is constrained to 24% of the total diffuse flux.
The majority of astrophysical neutrinos have undetermined origins. The IceCube Neutrino Observatory has observed astrophysical neutrinos but has not yet identified their sources. Blazars are promising source candidates, but previous searches for neutrino emission from populations of blazars detected in ≳GeV gamma rays have not observed any significant neutrino excess. Recent findings in multimessenger astronomy indicate that high-energy photons, coproduced with high-energy neutrinos, are likely to be absorbed and reemitted at lower energies. Thus, lower-energy photons may be better indicators of TeV–PeV neutrino production. This paper presents the first time-integrated stacking search for astrophysical neutrino emission from MeV-detected blazars in the first Fermi Large Area Telescope low energy (1FLE) catalog using ten years of IceCube muon–neutrino data. The results of this analysis are found to be consistent with a background-only hypothesis. Assuming an E -2 neutrino spectrum and proportionality between the blazars MeV gamma-ray fluxes and TeV–PeV neutrino flux, the upper limit on the 1FLE blazar energy-scaled neutrino flux is determined to be 1.64 x 10 -12 TeV cm -2 s -1 at 90% confidence level. This upper limit is approximately 1% of IceCube’s diffuse muon–neutrino flux measurement.
The extragalactic background light (EBL) contains all the radiation emitted by nuclear and accretion processes in stars and compact objects since the epoch of recombination. Measuring the EBL density directly is challenging, especially in the near-to-far-infrared wave band, mainly due to the zodiacal light foreground. Instead, gamma-ray astronomy offers the possibility to indirectly set limits on the EBL by studying the effects of gamma-ray absorption in the very high energy (VHE: >100 GeV) spectra of distant blazars. The High Altitude Water Cherenkov Gamma Ray Observatory (HAWC) is one of the few instruments sensitive to gamma rays with energies above 10 TeV. This offers the opportunity to probe the EBL in the near/mid-IR region: λ = 1–100 μm. In this study, we fit physically motivated emission models to Fermi-LAT gigaelectronvolt data to extrapolate the intrinsic teraelectronvolt spectra of blazars. We then simulate a large number of absorbed spectra for different randomly generated EBL model shapes and calculate Bayesian credible bands in the EBL intensity space by comparing and testing the agreement between the absorbed spectra and HAWC extragalactic observations of two blazars. The resulting bands are in agreement with current EBL lower and upper limits, showing a downward trend toward higher wavelength values λ > 10 μm also observed in previous measurements.