Source author record

G. Giacinti

G. Giacinti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2026arXiv

An Enigmatic PeVatron in an Area around HII Region G35.6$-$0.5

Identifying Galactic PeVatrons (PeV particle accelerators) from the ultra-high-energy (UHE, >100 TeV) $γ$-ray sources plays a crucial role in revealing the origin of Galactic cosmic rays. The UHE source 1LHAASO J1857+0203u is suggested to be associated with HESS J1858+020, which may be attributed to the possible PeVatron candidate supernova remnant (SNR) G35.6$-$0.4 or HII region G35.6$-$0.5. We perform detailed analysis on the very-high-energy and UHE $γ$-ray emissions towards this region with data from the Large High Altitude Air Shower Observatory (LHAASO). 1LHAASO J1857+0203u is detected with a significance of 11.6$σ$ above 100 TeV, indicating the presence of a PeVatron. It has an extension of $\sim 0.18^\circ$ with a power-law (PL) spectral index of $\sim$2.5 in 1-25 TeV and a point-like emission with a PL spectral index of $\sim$3.2 above 25 TeV. Using the archival CO and HI data, we identify some molecular and atomic clouds that may be associated with the TeV $γ$-ray emissions. Our modelling indicates that the TeV $γ$-ray emissions are unlikely to arise from the clouds illuminated by the protons that escaped from SNR G35.6$-$0.4. In the scenario that HII region G35.6$-$0.5 could accelerate particles to the UHE band, the observed GeV-TeV $γ$-ray emission could be well explained by a hadronic model with a PL spectral index of $\sim$2.0 and cutoff energy of $\sim$450 TeV. However, an evolved pulsar wind nebula origin cannot be ruled out.

preprint2026arXiv

An Ultrahigh-energy $γ$-ray Bubble Powered by a Super PeVatron

We report the detection of a $γ$-ray bubble spanning at least 100$\rm deg^2$ in ultra high energy (UHE) up to a few PeV in the direction of the star-forming region Cygnus X, implying the presence Super PeVatron(s) accelerating protons to at least 10 PeV. A log-parabola form with the photon index $Γ(E) = (2.71 \pm 0.02) + (0.11 \pm 0.02) \times \log_{10} (E/10 \ {\rm TeV})$ is found fitting the gamma-ray energy spectrum of the bubble well. UHE sources, `hot spots' correlated with very massive molecular clouds, and a quasi-spherical amorphous $γ$-ray emitter with a sharp central brightening are observed in the bubble. In the core of $\sim 0.5^{\circ}$, spatially associating with a region containing massive OB association (Cygnus OB2) and a microquasar (Cygnus X-3), as well as previously reported multi-TeV sources, an enhanced concentration of UHE $γ$-rays are observed with 2 photons at energies above 1 PeV. The general feature of the bubble, the morphology and the energy spectrum, are reasonably reproduced by the assumption of a particle accelerator in the core, continuously injecting protons into the ambient medium.

preprint2026arXiv

Constraining the Cosmic-ray Energy Based on Observations of Nearby Galaxy Clusters by LHAASO

Galaxy clusters act as reservoirs of high-energy cosmic rays (CRs). As CRs propagate through the intracluster medium, they generate diffuse $γ$-rays detectable by arrays such as LHAASO. These $γ$-rays result from proton-proton ($pp$) collisions of very high-energy cosmic rays (VHECRs) or inverse Compton (IC) scattering of positron-electron pairs created by $pγ$ interactions of ultra-high-energy cosmic rays (UHECRs). We analyzed diffuse $γ$-ray emission from the Coma, Perseus, and Virgo clusters using LHAASO data. Diffuse emission was modeled as a disk of radius $R_{500}$ for each cluster while accounting for point sources. No significant diffuse emission was detected, yielding 95\% confidence level (C.L.) upper limits on the $γ$-ray flux: for WCDA (1-25~TeV) and KM2A ($>25$~TeV), less than $(49.4, 13.7, 54.0)$ and $(1.34, 1.14, 0.40) \times 10^{-14}$~ph~cm$^{-2}$~s$^{-1}$ for Coma, Perseus, and Virgo, respectively. The $γ$-ray upper limits can be used to derive model-independent constraints on the integral energy of CRp above 10~TeV (corresponding to the LHAASO observational range $>1$~TeV under the $pp$ scenario) to be less than $(1.96, 0.59, 0.08) \times 10^{61}$~erg. The absence of detectable annuli/ring-like structures, indicative of cluster accretion or merging shocks, imposes further constraints on models in which the UHECRs are accelerated in the merging shocks of galaxy clusters.

preprint2026arXiv

Deep view of Composite SNR CTA1 with LHAASO in $γ$-rays up to 300 TeV

The ultra-high-energy (UHE) gamma-ray source 1LHAASO J0007+7303u is positionally associated with the composite SNR CTA1 that is located at high Galactic Latitude $b\approx 10.5^\circ$. This provides a rare opportunity to spatially resolve the component of the pulsar wind nebula (PWN) and supernova remnant (SNR) at UHE. This paper conducted a dedicated data analysis of 1LHAASO J0007+7303u using the data collected from December 2019 to July 2023. This source is well detected with significances of 21$σ$ and 17$σ$ at 8$-$100 TeV and $>$100 TeV, respectively. The corresponding extensions are determined to be 0.23$^{\circ}\pm$0.03$^{\circ}$ and 0.17$^{\circ}\pm$0.03$^{\circ}$. The emission is proposed to originate from the relativistic electrons and positrons accelerated within the PWN of PSR J0007+7303. The energy spectrum is well described by a power-law with an exponential cutoff function $dN/dE = (42.4\pm4.1)(\frac{E}{20\rm\ TeV})^{-2.31\pm0.11}\exp(-\frac{E}{110\pm25\rm\ TeV})$ $\rm\ TeV^{-1}\ cm^{-2}\ s^{-1}$in the energy range from 8 TeV to 300 TeV, implying a steady-state parent electron spectrum $dN_e/dE_e\propto (\frac{E_e}{100\rm\ TeV})^{-3.13\pm0.16}\exp[(\frac{-E_e}{373\pm70\rm\ TeV})^2]$ at energies above $\approx 50 \rm\ TeV$. The cutoff energy of the electron spectrum is roughly equal to the expected current maximum energy of particles accelerated at the PWN terminal shock. Combining the X-ray and gamma-ray emission, the current space-averaged magnetic field can be limited to $\approx 4.5\rm\ μG$. To satisfy the multi-wavelength spectrum and the $γ$-ray extensions, the transport of relativistic particles within the PWN is likely dominated by the advection process under the free-expansion phase assumption.

preprint2026arXiv

Energy calibration of LHAASO-KM2A using the cosmic ray Moon shadow

We present a precise measurement of the westward, rigidity-dependent shift of the Moon's shadow using three and a half years of cosmic-ray data collected by the Kilometer Square Array (KM2A) of the Large High Altitude Air Shower Observatory (LHAASO). These measurements enable us to calibrate the detector energy response in the range 20-260 TeV, with results showing excellent agreement with the response derived from Monte Carlo (MC) simulations of the KM2A detector. We also measure a best-fit parameter $ε= 0.015 \pm 0.08$, corresponding to a 95% confidence interval of [-14%, +17%] for the energy-scale estimation. This result establishes the exceptional accuracy of the KM2A-MC in simulating the detector's response within this energy range.

preprint2026arXiv

Energy-Dependent Shifts of Medium-Scale Anisotropies in Very-High-Energy Cosmic Rays Observed by LHAASO-KM2A

Small deviations from isotropy in the arrival directions of Galactic cosmic rays serve as a unique probe of the local magnetic environment. In this Letter, we report observations of medium-scale anisotropies (MSA) at energies above 10 TeV using the LHAASO-KM2A array. Our analysis identifies four regions of excess and four regions of deficit, each spanning angular scales of approximately ten degrees. Crucially, we detect significant energy-dependent shifts in the centroids of two excess regions: Region B and the newly identified Region $\mathrm{\widetilde{D}}$. We also characterize the energy evolution of the fractional relative intensity across both excess and deficit regions. These findings imply that the observed anisotropies are shaped by the specific realization of the local turbulent magnetic field within the cosmic ray scattering length. Such energy-dependent behaviors impose strict constraints on local turbulence models and cosmic ray propagation theories.

preprint2026arXiv

Evidence for particle acceleration approaching PeV energies in the W51 complex

The $γ$-ray emission from the W51 complex is widely acknowledged to be attributed to the interaction between the cosmic rays (CRs) accelerated by the shock of supernova remnant (SNR) W51C and the dense molecular clouds in the adjacent star-forming region, W51B. However, the maximum acceleration capability of W51C for CRs remains elusive. Based on observations conducted with the Large High Altitude Air Shower Observatory (LHAASO), we report a significant detection of $γ$ rays emanating from the W51 complex, with energies from 2 TeV to 200 TeV. The LHAASO measurements, for the first time, extend the $γ$-ray emission from the W51 complex beyond 100 TeV and reveal a significant spectrum bending at tens of TeV. By combining the ``$π^0$-decay bump" featured data from Fermi-LAT, the broadband $γ$-ray spectrum of the W51 region can be well-characterized by a simple pp-collision model. The observed spectral bending feature suggests an exponential cutoff at $\sim400$~TeV or a power-law break at $\sim200$~TeV in the CR proton spectrum, most likely providing the first evidence of SNRs serving as CR accelerators approaching the PeV regime. Additionally, two young star clusters within W51B could also be theoretically viable to produce the most energetic $γ$ rays observed by LHAASO. Our findings strongly support the presence of extreme CR accelerators within the W51 complex and provide new insights into the origin of Galactic CRs.

preprint2026arXiv

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622 + 3749 Observed by LHAASO-KM2A

We report the discovery of an extended very-high-energy (VHE) gamma-ray source around the location of the middle-aged (207.8 kyr) pulsar PSR J0622+3749 with the Large High Altitude Air Shower Observatory (LHAASO). The source is detected with a significance of $8.2σ$ for $E>25$~TeV assuming a Gaussian template. The best-fit location is (R.A., Dec.)$=(95^{\circ}\!.47\pm0^{\circ}\!.11,\,37^{\circ}\!.92 \pm0^{\circ}\!.09)$, and the extension is $0^{\circ}\!.40\pm0^{\circ}\!.07$. The energy spectrum can be described by a power-law spectrum with an index of ${-2.92 \pm 0.17_{\rm stat} \pm 0.02_{\rm sys} }$. No clear extended multi-wavelength counterpart of the LHAASO source has been found from the radio to sub-TeV bands. The LHAASO observations are consistent with the scenario that VHE electrons escaped from the pulsar, diffused in the interstellar medium, and scattered the interstellar radiation field. If interpreted as the pulsar halo scenario, the diffusion coefficient, inferred for electrons with median energies of $\sim160$~TeV, is consistent with those obtained from the extended halos around Geminga and Monogem and much smaller than that derived from cosmic ray secondaries. The LHAASO discovery of this source thus likely enriches the class of so-called pulsar halos and confirms that high-energy particles generally diffuse very slowly in the disturbed medium around pulsars.

preprint2026arXiv

LHAASO Detection of Ultra-High-Energy Gamma-Ray Emission toward the Giant Molecular Clouds

The $γ$-ray from Giant molecular clouds (GMCs) is regarded as the most ideal tool to perform in-situ measurement of cosmic ray (CR) density and spectra in our Galaxy. We report the first detection of $γ$-ray emissions in the very-high-energy (VHE) domain from the five nearby GMCs with a stacking analysis based on a 4.5-year $γ$-ray observation with the Large High Altitude Air Shower Observatory (LHAASO) experiment. The spectral energy distributions derived from the GMCs are consistent with the expected $γ$-ray flux produced via CR interacting with the ISM in the energy interval 1 - 100 $~\rm$ TeV. In addition, we investigate the presence of the CR spectral `knee' by introducing a spectral break in the $γ$-ray data. While no significant evidence for the CR knee is found, the current KM2A measurements from GMCs strongly favor a proton CR knee located above 0.9$~\rm$ PeV, which is consistent with the latest measurement of the CR spectrum by ground-based experiments.

preprint2026arXiv

Measurement of attenuation length of the muon content in extensive air showers from 0.3 to 30 PeV with LHAASO

The attenuation length of the muon content in extensive air showers provides important information regarding the generation and development of air showers. This information can be used not only to improve the description of such showers but also to test fundamental models of hadronic interactions. Using data from the LHAASO-KM2A experiment, the development of the muon content in high-energy air showers was studied. The attenuation length of muon content in the air showers was measured from experimental data in the energy range from 0.3 to 30 PeV using the constant intensity cut method. By comparing the attenuation length of the muon content with predictions from high-energy hadronic interaction models (QGSJET-II-04, SIBYLL 2.3d, and EPOS-LHC), it is evident that LHAASO results are significantly shorter than those predicted by the first two models (QGSJET-II-04 and SIBYLL 2.3d) but relatively close to those predicted by the third model (EPOS-LHC). Thus, the LHAASO data favor the EPOS-LHC model over the other two models. The three interaction models confirmed an increasing trend in the attenuation length as the cosmic-ray energy increases.

preprint2026arXiv

Measurement of Very-high-energy Diffuse Gamma-ray Emissions from the Galactic Plane with LHAASO-WCDA

The diffuse Galactic gamma-ray emission is a very important tool used to study the propagation and interaction of cosmic rays in the Milky Way. In this work, we report the measurements of the diffuse emission from the Galactic plane, covering Galactic longitudes from $15^{\circ}$ to $235^{\circ}$ and latitudes from $-5^{\circ}$ to $+5^{\circ}$, in an energy range of 1 TeV to 25 TeV, with the Water Cherenkov Detector Array (WCDA) of the Large High Altitude Air Shower Observatory (LHAASO). After masking the sky regions of known sources, the diffuse emission is detected with $24.6σ$ and $9.1σ$ significance in the inner Galactic plane and outer Galactic plane, respectively. The WCDA spectra in both regions can be well described by a power-law function, with spectral indices of $-2.67\pm0.05_{\rm stat}$ in the inner region and $-2.83\pm0.19_{\rm stat}$ in the outer region, respectively. Combined with the Square Kilometer Array (KM2A) measurements at higher energies, a clear softening of the spectrum is found in the inner region, with change of spectral indices by $\sim0.5$ at a break energy around $30$ TeV. The fluxes of the diffuse emission are higher by a factor of $1.5-2.7$ than the model prediction assuming local CR spectra and the gas column density, which are consistent with those measured by the KM2A. Along Galactic longitude, the spatial distribution of the diffuse emission shows deviation from that of the gas column density. The spectral shape of the diffuse emission are possibly variation in different longitude region. The WCDA measurements bridge the gap between the low-energy measurements by space detectors and the ultra-high-energy observations by LHAASO-KM2A and other experiments. These results suggest that improved modeling of the wide-band diffuse emission is required.

preprint2026arXiv

Transient Large-Scale Anisotropy in TeV Cosmic Rays due to an Interplanetary Coronal Mass Ejection

Large- or medium-scale cosmic ray anisotropy at TeV energies has not previously been confirmed to vary with time. Transient anisotropy changes have been observed below 150 GeV, especially near the passage of an interplanetary shock and coronal mass ejection containing a magnetic flux rope ejected by a solar storm, which can trigger a geomagnetic storm with practical consequences. In such events, cosmic rays provide remote sensing of the magnetic field properties. Here we report the observation of transient large-scale anisotropy in TeV cosmic ray ions using data from the Large High Altitude Air Shower Observatory (LHAASO). We analyze hourly skymaps of the transient cosmic ray intensity excess or deficit, the gradient of which indicates the direction and magnitude of transient large-scale anisotropy across the field of view. We observe enhanced anisotropy above typical hourly fluctuations with $>$5$σ$ significance during some hours of November 4, 2021, in separate data sets for four primary cosmic ray energy ranges of median energy from $E$=0.7 to 3.1 TeV. The gradient varies with energy as $E^γ$, where $γ\approx-0.5$. At a median energy $\leq$1.0 TeV, this gradient corresponds to a dipole anisotropy of at least 1\%, or possibly a weaker anisotropy of higher order. This new type of observation opens the opportunity to study interplanetary magnetic structures using air shower arrays around the world, complementing existing in situ and remote measurements of plasma properties.

preprint2025arXiv

Observation of the $γ$-ray Emission from W43 with LHAASO

In this paper, we report the detection of the very-high-energy (VHE, $ 100{\rm\ GeV} < E < 100{\rm\ TeV} $) and ultra-high-energy (UHE, $E > 100\rm\ TeV$) $γ$-ray emissions from the direction of the young star-forming region W43, observed by the Large High Altitude Air Shower Observation (LHAASO). The extended $γ$-ray source was detected with a significance of ${\sim}16\,σ$ by KM2A and ${\sim}17\,σ$ by WCDA, respectively. The angular extension of this $γ$-ray source is about 0.5 degrees, corresponding to a physical size of about 50 pc. We discuss the origin of the $γ$-ray emission and possible cosmic ray acceleration in the W43 region using multi-wavelength data. Our findings suggest that W43 is likely another young star cluster capable of accelerating cosmic rays (CRs) to at least several hundred TeV.

preprint2025arXiv

Study of Ultra-High-Energy Gamma-Ray Source 1LHAASO J0056+6346u and Its Possible Origins

We report a dedicated study of the newly discovered extended UHE $γ$-ray source 1LHAASO J0056+6346u. Analyzing 979 days of LHAASO-WCDA data and 1389 days of LHAASO-KM2A data, we observed a significant excess of $γ$-ray events with both WCDA and KM2A. Assuming a point power-law source with a fixed spectral index, the significance maps reveal excesses of ${\sim}12.65\,σ$, ${\sim}22.18\,σ$, and ${\sim}10.24\,σ$ in the energy ranges of 1--25 TeV, 25--100 TeV, and $> 100$ TeV, respectively. We use a 3D likelihood algorithm to derive the morphological and spectral parameters, and the source is detected with significances of $12.65\,σ$ by WCDA and $25.27\,σ$ by KM2A. The best-fit positions derived from WCDA and KM2A data are (R.A. = $13.96^\circ\pm0.09^\circ$, Decl. = $63.92^\circ\pm0.05^\circ$) and (R.A. = $14.00^\circ\pm0.05^\circ$, Decl. = $63.79^\circ\pm0.02^\circ$), respectively. The angular size ($r_{39}$) of 1LHAASO J0056+6346u is $0.34^\circ\pm0.04^\circ$ at 1--25 TeV and $0.24^\circ\pm0.02^\circ$ at $> 25$ TeV. The differential flux of this UHE $γ$-ray source can be described by an exponential cutoff power-law function: $(2.67\pm0.25) \times 10^{-15} (E/20\,\text{TeV})^{-1.97\pm0.10} e^{-E/(55.1\pm7.2)\,\text{TeV}} \,\text{TeV}^{-1}\,\text{cm}^{-2}\,\text{s}^{-1}$. To explore potential sources of $γ$-ray emission, we investigated the gas distribution around 1LHAASO J0056+6346u. 1LHAASO J0056+6346u is likely to be a TeV PWN powered by an unknown pulsar, which would naturally explain both its spatial and spectral properties. Another explanation is that this UHE $γ$-ray source might be associated with gas content illuminated by a nearby CR accelerator, possibly the SNR candidate G124.0+1.4.

preprint2022arXiv

Galactic Science with the Southern Wide-field Gamma-ray Observatory

The Southern Wide-field Gamma-ray Observatory (SWGO) is a proposed ground-based gamma-ray detector that will be located in the Southern Hemisphere and is currently in its design phase. In this contribution, we will outline the prospects for Galactic science with this Observatory. Particular focus will be given to the detectability of extended sources, such as gamma-ray halos around pulsars; optimisation of the angular resolution to mitigate source confusion between known TeV sources; and studies of the energy resolution and sensitivity required to study the spectral features of PeVatrons at the highest energies. Such a facility will ideally complement contemporaneous observatories in studies of high energy astrophysical processes in our Galaxy.

preprint2022arXiv

HAWC Study of the Ultra-High-Energy Spectrum of MGRO J1908+06

We report TeV gamma-ray observations of the ultra-high-energy source MGRO J1908+06 using data from the High Altitude Water Cherenkov (HAWC) Observatory. This source is one of the highest-energy known gamma-ray sources, with emission extending past 200 TeV. Modeling suggests that the bulk of the TeV gamma-ray emission is leptonic in nature, driven by the energetic radio-faint pulsar PSR J1907+0602. Depending on what assumptions are included in the model, a hadronic component may also be allowed. Using the results of the modeling, we discuss implications for detection prospects by multi-messenger campaigns.

preprint2022arXiv

Signatures of anisotropic diffusion around PeVatrons in 100 TeV gamma-ray data

The Tibet AS$γ$ collaboration has reported a diffuse $γ$-ray emission signal from the Galactic Plane. We consider the possibility that the diffuse emission from the outer Galactic Plane at the highest energies is produced by cosmic rays spreading from a single supernova-type source either in the Local or Perseus arm of the Milky Way. We show that anisotropic diffusion of multi-PeV cosmic rays along the Galactic magnetic field can produce an extended source spanning ten(s) of degrees on the sky, with a flux-per-unit-solid-angle consistent with Tibet AS$γ$ measurements. Observations of this new type of very extended sources, and measurements of their morphologies, can be used to characterize the anisotropic diffusion of PeV cosmic rays in the Galactic magnetic field, and to constrain the locations and properties of past PeVatrons.

preprint2022arXiv

The first days of type II-P core collapse supernovae in the gamma-ray range

Type II-P supernovæ~(SNe), the most common core-collapse SNe type, result from the explosions of red supergiant stars. Their detection in the radio domain testifies of the presence of relativistic electrons, and shows that they are potentially efficient energetic particle accelerators. If hadrons can also be accelerated, these energetic particles are expected to interact with the surrounding medium to produce a gamma-ray signal even in the multi--TeV range. The intensity of this signal depends on various factors, but an essential one is the density of the circumstellar medium. Such a signal should however be limited by electron-positron pair production arising from the interaction of the gamma-ray photons with optical photons emitted by the supernova photosphere, which can potentially degrade the gamma-ray signal by over ten orders of magnitude in the first days/weeks following the explosion. We calculate the gamma-gamma opacity from a detailed modelling of the time evolution of the forward shock and supernova photosphere, taking a full account of the non-isotropy of the photon interactions. We discuss the time-dependent gamma-ray TeV emission from type II-P SNe as a function of the stellar progenitor radius and mass-loss rate, as well as the explosion energy and mass of the ejected material. We evaluate the detectability of the SNe with the next generation of Cherenkov telescopes. We find that, while most extragalactic events may be undetectable, type II-P SNe exploding in our Galaxy or in the Magellanic Clouds should be detected by gamma-ray observatories such as the upcoming Cherenkov Telescope Array.

preprint2021arXiv

3HWC: The Third HAWC Catalog of Very-High-Energy Gamma-ray Sources

We present a new catalog of TeV gamma-ray sources using 1523 days of data from the High Altitude Water Cherenkov (HAWC) observatory. The catalog represents the most sensitive survey of the Northern gamma-ray sky at energies above several TeV, with three times the exposure compared to the previous HAWC catalog, 2HWC. We report 65 sources detected at $\geq$ 5 sigma significance, along with the positions and spectral fits for each source. The catalog contains eight sources that have no counterpart in the 2HWC catalog, but are within $1^\circ$ of previously detected TeV emitters, and twenty sources that are more than $1^\circ$ away from any previously detected TeV source. Of these twenty new sources, fourteen have a potential counterpart in the fourth \textit{Fermi} Large Area Telescope catalog of gamma-ray sources. We also explore potential associations of 3HWC sources with pulsars in the ATNF pulsar catalog and supernova remnants in the Galactic supernova remnant catalog.

preprint2021arXiv

Cosmic-ray current-driven instabilities -- revisiting environmental conditions

The growth of magneto-hydrodynamic fluctuations relevant to cosmic ray confinement in and near their sources, and the effects of local plasma conditions is revisited. We consider cases where cosmic rays penetrate a medium which may contain a fraction of neutral particles, and explore the possible effects of high-order cosmic-ray anisotropies. An algorithm for calculating the dispersion relation for arbitrary distributions, and anisotropies is presented, and a general solution for power-law cosmic-ray distributions is provided. Implications for the resulting instabilities near to strong Galactic cosmic-ray sources are discussed. We argue that cosmic-ray streaming in weakly ionised plasmas eliminates the need for the existence of an evanescent band in the dispersion relation, a conclusion which may be confirmed by gamma-ray observations. The necessity for additional multi-scale numerical simulations is highlighted, as understanding the non-linear behaviour is crucial.

preprint2021arXiv

Spectrum and Morphology of the Very-High-Energy Source HAWC J2019+368

The MGRO J2019+37 region is one of the brightest sources in the sky at TeV energies. It was detected in the 2 year HAWC catalog as 2HWC J2019+367 and here we present a detailed study of this region using data from HAWC. This analysis resolves the region into two sources: HAWC J2019+368 and HAWC J2016+371. We associate HAWC J2016+371 with the evolved supernova remnant CTB 87, although its low significance in this analysis prevents a detailed study at this time. An investigation of the morphology (including possible energy dependent morphology) and spectrum for HAWC J2019+368 is the focus of this work. We associate HAWC J2019+368 with PSR J2021+3651 and its X-ray pulsar wind nebula, the Dragonfly nebula. Modeling the spectrum measured by HAWC and Suzaku reveals a $\sim$7 kyr pulsar and nebula system producing the observed emission at X-ray and $γ$-ray energies.

preprint2021arXiv

Ultra-high energy Inverse Compton emission from Galactic electron accelerators

It is generally held that >100 TeV emission from astrophysical objects unambiguously demonstrates the presence of PeV protons or nuclei, due to the unavoidable Klein-Nishina suppression of inverse Compton emission from electrons. However, in the presence of inverse Compton dominated cooling, hard high-energy electron spectra are possible. We show that the environmental requirements for such spectra can naturally be met in spiral arms, and in particular in regions of enhanced star formation activity, the natural locations for the most promising electron accelerators: powerful young pulsars. Our scenario suggests a population of hard ultra-high energy sources is likely to be revealed in future searches, and may also provide a natural explanation for the 100 TeV sources recently reported by HAWC.

preprint2020arXiv

HAWC J2227+610 and its association with G106.3+2.7, a new potential Galactic PeVatron

We present the detection of VHE gamma-ray emission above 100 TeV from HAWC J2227+610 with the HAWC observatory. Combining our observations with previously published results by VERITAS, we interpret the gamma-ray emission from HAWC J2227+610 as emission from protons with a lower limit in their cutoff energy of 800 TeV. The most likely source of the protons is the associated supernova remnant G106.3+2.7, making it a good candidate for a Galactic PeVatron. However, a purely leptonic origin of the observed emission cannot be excluded at this time.

preprint2019arXiv

On the TeV Halo Fraction in gamma-ray bright Pulsar Wind Nebulae

The discovery of extended TeV emission around the Geminga and PSR B0656+14 pulsars, with properties consistent with free particle propagation in the interstellar medium (ISM), has sparked considerable discussion on the possible presence of such halos in other systems. Here we make an assessment of the current TeV source population associated with energetic pulsars, in terms of size and estimated energy density. Based on two alternative estimators we conclude that a large majority of the known TeV sources have emission originating in the zone energetically and dynamically dominated by the pulsar (i.e. the pulsar wind nebula), rather than from a halo of particles diffusing in to the ISM. Furthermore, whilst the number of established halos will surely increase in the future, we find that it is unlikely that such halos contribute significantly to the total TeV $γ$-ray luminosity from electrons accelerated in PWN.

preprint2015arXiv

Collisionless Shocks and TeV Neutrinos before Supernova Shock Breakout from an Optically Thick Wind

During a supernova explosion, a radiation-dominated shock (RDS) travels through its progenitor. A collisionless shock (CS) is usually assumed to replace it during shock breakout (SB). We demonstrate here that for some realistic progenitors enshrouded in optically thick winds, such as possibly SN 2008D, a CS forms deep inside the wind, soon after the RDS leaves the core, and therefore significantly before SB. The RDS does not survive the transition from the core to the thick wind when the wind close to the core is not sufficiently dense to compensate for the $r^{-2}$ dilution of photons due to shock curvature. This typically happens when the shock velocity is $\lesssim 0.1 {\rm c} \, (\frac{u_{\rm w}}{10\,{\rm km/s}}) (\frac{\dot{M}}{5 \cdot 10^{-4} \, {\rm M}_\odot {\rm /yr}})^{-1} (\frac{r_\ast}{10^{13}\,{\rm cm}})$, where $u_{\rm w}$, $\dot{M}$ and $r_\ast$ are respectively the wind velocity, mass-loss rate and radius of the progenitor star. The radiative CS results in a hard spectrum of the photon flash at breakout, which would produce an X-ray flash. Cosmic ray acceleration would start before SB, for such progenitors. A fraction of secondary TeV neutrinos can reach the observer up to more than ten hours before the first photons from breakout, providing information on the invisible layers of the progenitor.

preprint2015arXiv

Escape model for Galactic cosmic rays and an early extragalactic transition

We show that the cosmic ray (CR) knee can be entirely explained by energy-dependent CR leakage from the Milky Way, with an excellent fit to all existing data. We test this hypothesis calculating the trajectories of individual CRs in the Galactic magnetic field. We find that the CR escape time $τ_{\rm esc}(E)$ exhibits a knee-like structure around $E/Z={\rm few}\times 10^{15}$ eV for small coherence lengths and strengths of the turbulent magnetic field. The resulting intensities for different groups of nuclei are consistent with the ones determined by KASCADE and KASCADE-Grande, using simple power-laws as injection spectra. The transition from Galactic to extragalactic CRs is terminated at $\approx 2\times 10^{18}$ eV, while extragalactic CRs contribute sizeable to the subdominant proton flux already for $\gtrsim 2\times 10^{16}$ eV. The natural source of extragalactic CRs in the intermediate energy region up to the ankle are in this model normal and starburst galaxies. The escape model provides a good fit to $\ln(A)$ data; it predicts that the phase of the CR dipole varies strongly in the energy range between $1\times 10^{17}$ and $3\times 10^{18}$ eV, while our estimate for the dipole magnitude is consistent with observations.

preprint2015arXiv

Unified model for cosmic rays above $10^{17}$ eV and the diffuse gamma-ray and neutrino backgrounds

We investigate how the extragalactic proton component derived within the "escape model" can be explained by astrophysical sources. We consider as possible cosmic ray (CR) sources normal/starburst galaxies and radio-loud active galactic nuclei (AGN). We find that the contribution to the total extragalactic proton flux from normal and starburst galaxies is only subdominant and does not fit the spectral shape deduced in the escape model. In the case of radio-loud AGN, we show that the complete extragalactic proton spectrum can be explained by a single source population, BL Lac/FR I, for any of the potential acceleration sites in these sources. We calculate the diffuse neutrino and $γ$-ray fluxes produced by these CR protons interacting with gas inside their sources. For a spectral slope of CRs close to $α=2.1-2.2$ as suggested by shock acceleration, we find that these UHECR sources contribute the dominant fraction of both the isotropic $γ$-ray background and of the extragalactic part of the astrophysical neutrino signal observed by IceCube.

preprint2014arXiv

Cosmic ray propagation in the interstellar magnetic fields

The propagation of TeV-PeV cosmic rays (CR) in our Galaxy can be described as a diffusive process. We discuss here two effects, with important observational consequences, that cannot be predicted by the diffusion approximation in its usual form. First, we present an explanation for the CR anisotropies observed at small angular scales on the sky. We show that the local magnetic field configuration within a CR mean free path from Earth naturally results in CR flux anisotropies at small and medium scales. Second, we point out that TeV-PeV CRs should be expected to diffuse strongly anisotropically in the interstellar medium on scales smaller than the maximum scale of spatial fluctuations of the field, $\sim 100$ pc. This notably questions the usual assumptions on CR diffusion around sources.

preprint2014arXiv

Explaining the Spectra of Cosmic Ray Groups above the Knee by Escape from the Galaxy

We investigate the possibility that the cosmic ray (CR) knee is entirely explained by the energy-dependent CR leakage from the Milky Way. We test this hypothesis calculating the trajectories of individual CRs with energies between $E/Z=10^{14}$ eV and $10^{17}$ eV propagating them in the regular and turbulent Galactic magnetic field. We find a knee-like structure of the CR escape time $τ_{\rm esc}(E)$ around $E/Z={\rm few}\times 10^{15}$ eV for a coherence length $l_{\rm c} \simeq 2$ pc of the turbulent field, while the decrease of $τ_{\rm esc}(E)$ slows down around $E/Z\simeq 10^{16}$ eV in models with a weak turbulent magnetic field. Assuming that the injection spectra of CR nuclei are power-laws, the resulting CR intensities in such a turbulence are consistent with the energy spectra of CR nuclei determined by KASCADE and KASCADE-Grande. We calculate the resulting CR dipole anisotropy as well as the source rate in this model.

preprint2013arXiv

Anisotropic Cosmic Ray Diffusion and its Implications for Gamma-Ray Astronomy

Analyses of TeV-PeV cosmic ray (CR) diffusion around their sources usually assume either isotropic diffusion or anisotropic diffusion due to the regular Galactic magnetic field. We show that none of them are adequate on distances smaller than the maximal scale Lmax ~ 100 pc of fluctuations in the turbulent interstellar magnetic field. As a result, we predict anisotropic gamma-ray emissions around CR proton and electron sources, even for uniform densities of target gas. The centers of extended emission regions may have non-negligible offsets from their sources, leading to risks of misidentification. Gamma-rays from CR filaments have steeper energy spectra than those from surrounding regions. We point out that gamma-ray telescopes can be used in the future as a new way to probe and deduce the parameters of the interstellar magnetic field.

preprint2012arXiv

Cosmic Ray Anisotropy as Signature for the Transition from Galactic to Extragalactic Cosmic Rays

We constrain the energy at which the transition from Galactic to extragalactic cosmic rays occurs by computing the anisotropy at Earth of cosmic rays emitted by Galactic sources. Since the diffusion approximation starts to loose its validity for E/Z >~ 10^(16-17) eV, we propagate individual cosmic rays (CRs) using Galactic magnetic field models and taking into account both their regular and turbulent components. The turbulent field is generated on a nested grid which allows spatial resolution down to fractions of a parsec. Assuming sufficiently frequent Galactic CR sources, the dipole amplitude computed for a mostly light or intermediate primary composition exceeds the dipole bounds measured by the Auger collaboration around E ~ 10^18 eV. Therefore, a transition at the ankle or above would require a heavy composition or a rather extreme Galactic magnetic field with strength >~ 10 muG. Moreover, the fast rising proton contribution suggested by KASCADE-Grande data between 10^17 eV and 10^18 eV should be of extragalactic origin. In case heavy nuclei dominate the flux at E >~ 10^18 eV, the transition energy can be close to the ankle, if Galactic CRs are produced by sufficiently frequent transients as e.g. magnetars.

preprint2012arXiv

Filamentary Diffusion of Cosmic Rays on Small Scales

We investigate the diffusion of cosmic rays (CR) close to their sources. Propagating individual CRs in purely isotropic turbulent magnetic fields with maximal scale of spatial variations Lmax, we find that CRs diffuse anisotropically at distances r <~ Lmax from their sources. As a result, the CR densities around the sources are strongly irregular and show filamentary structures. We determine the transition time t* to standard diffusion as t* ~ 10^4 yr (Lmax/150 pc)^b (E/PeV)^(-g) (Brms/4 muG)^g, with b ~ 2 and g = 0.25-0.5 for a turbulent field with Kolmogorov power spectrum. We calculate the photon emission due to CR interactions with gas and the resulting irregular source images.

preprint2012arXiv

Local Magnetic Turbulence and TeV-PeV Cosmic Ray Anisotropies

In the energy range from ~ 10^12 eV to ~ 10^15 eV, the Galactic cosmic ray flux has anisotropies both on large scales, with an amplitude of the order of 0.1%, and on scales between ~ 10 and ~ 30 degrees, with amplitudes smaller by a factor of a few. With a diffusion coefficient inferred from Galactic cosmic ray chemical abundances, the diffusion approximation predicts a dipolar anisotropy of comparable size, but does not explain the smaller scale anisotropies. We demonstrate here that energy dependent smaller scale anisotropies naturally arise from the local concrete realization of the turbulent magnetic field within the cosmic ray scattering length. We show how such anisotropies could be calculated if the magnetic field structure within a few tens of parsecs from Earth were known.

preprint2011arXiv

Method to Look for Imprints of Ultrahigh Energy Nuclei Sources

We propose a new method to search for heavy nuclei sources, on top of background, in the Ultra-High Energy Cosmic Ray data. We apply this method to the 69 events recently published by the Pierre Auger Collaboration and find a tail of events for which it reconstructs the source at a few degrees from the Virgo galaxy cluster. The reconstructed source is located at ~ 8.5 degrees from M87. The probability to have such a cluster of events in some random background and reconstruct the source position in any direction of the sky is about 7 x 10^(-3). The probability to reconstruct the source at less than 10 degrees from M87 in a data set already containing such a cluster of events is about 4 x 10^(-3). This may be a hint at the Virgo cluster as a bright ultra-high energy nuclei source. We investigate the ability of current and future experiments to validate or rule out this possibility, and discuss several alternative solutions which could explain the existing anisotropy in the Auger data.

preprint2011arXiv

Propagation of Ultrahigh Energy Nuclei in the Magnetic Field of our Galaxy

In this work, we present detailed simulations for propagation of ultra-high energy (UHE) heavy nuclei, with E > 60 EeV, within recent Galactic Magnetic Field (GMF) models. We investigate the impacts of the regular and turbulent components of the GMF. We show that with UHE heavy nuclei, there is no one-to-one correspondence between the arrival directions of cosmic rays (CR) measured at Earth and the direction of their extragalactic sources. Sources can have several distorted images on the sky. We compute images of galaxy clusters and of the supergalactic plane in recent GMF models and show the challenges, and possibilities, of "UHECR astronomy" with heavy nuclei. Finally, we present a quantitative study of the impact of the GMF on the (de-)magnification of source fluxes, due to magnetic lensing effects. We find that for 60 EeV iron nuclei, sources located in up to about one fifth of the sky would have their fluxes so strongly demagnified that they would not be detectable at Earth, even by the next generation of UHECR experiments.

preprint2011arXiv

Search for Individual Ultrahigh Energy Cosmic Ray Sources in the Future Data

We propose a new way to detect individual bright Ultra-High Energy Cosmic Ray (UHECR) sources above background if the Galactic Magnetic Field (GMF) gives the main contribution to UHECR deflections. This method can be directly applied to maps given by experiments. It consists in starting from at least two high energy events above 6x10^19 eV, and looking at lower energy tails. We test the efficiency of the method and investigate its dependence on different parameters. In case of detection, the source position and the local GMF deflection power are reconstructed. Both reconstructions are strongly affected by the turbulent GMF. With the parameters adopted in this study, for 68 % of reconstructed sources, the angular position is less than one degree from the real one. For typical turbulent field strengths of 4 μG at the Earth position and 1.5 kpc extension in the halo, one can reconstruct the deflection power with 25 % precision in 68 % of cases.

preprint2011arXiv

Search for Nuclei Sources in the Ultra-High Energy Cosmic Ray Data

We present a new method to search for heavy nuclei sources, on top of background, in the Ultra-High Energy Cosmic Ray data. We apply it to the 69 events with energies E>55 EeV published by the Pierre Auger Collaboration. We find a set of events for which the method reconstructs the source near the Virgo galaxy cluster. The probability to have a comparable set of events in some background is ~ 0.7 %. The reconstructed source is located at ~ 8.5 degrees from the active galaxy M87. The probability to reconstruct the source at less than 10 degrees from M87 for data already containing a comparable set of events is ~ 0.4 %. This may be a hint at the Virgo galaxy cluster as an ultra-high energy heavy nuclei source. We discuss the capability of current and near future experiments to test this possibility. Such a scenario gives a self-consistent description of the Auger anisotropy and composition data at the highest energies.

preprint2011arXiv

Ultrahigh Energy Nuclei in the Turbulent Galactic Magnetic Field

In this work we study how the turbulent component of the Galactic magnetic field (GMF) affects the propagation of ultrahigh energy heavy nuclei. We investigate first how the images of individual sources and of the supergalactic plane depend on the properties of the turbulent GMF. Then we present a quantitative study of the impact of the turbulent field on (de-) magnification of source fluxes, due to magnetic lensing effects. We also show that it is impossible to explain the Pierre Auger data assuming that all ultrahigh energy nuclei are coming from Cen A, even in the most favorable case of a strong, extended turbulent field in the Galactic halo.

preprint2010arXiv

Search for single sources of ultra high energy cosmic rays on the sky

In this paper, we suggest a new way to identify single bright sources of Ultra High Energy Cosmic Rays (UHECR) on the sky, on top of background. We look for doublets of events at the highest energies, E > 6 x 10^19 eV, and identify low energy tails, which are deflected by the Galactic Magnetic Field (GMF). For the sources which are detected, we can recover their angular positions on the sky within one degree from the real ones in 68% of cases. The reconstruction of the deflection power of the regular GMF is strongly affected by the value of the turbulent GMF. For typical values of 4 microG near the Earth, one can reconstruct the deflection power with 25% precision in 68% of cases.

preprint2010arXiv

Ultrahigh Energy Nuclei in the Galactic Magnetic Field

Observations are consistent with a significant fraction of heavy nuclei in the cosmic ray flux above a few times 10^19 eV. Such nuclei can be deflected considerably in the Galactic magnetic field, with important implications for the search of their sources. We perform detailed simulations of heavy nuclei propagation within recent Galactic magnetic field models. While such models are not yet sufficiently constrained to predict deflection maps in detail, we find general features of the distribution of (de-) magnified flux from sources. Since in most theoretical models sources of heavy nuclei are located in the local large scale structure of galaxies, we show examples of images of several nearby galaxy clusters and of the supergalactic plane. Such general features may be useful to develop efficient methods for source reconstruction from observed ultrahigh energy cosmic ray arrival directions.