Showing 1–40 of 40 results for author: Xi, S

Constraints on Ultra Heavy Dark Matter Properties from Dwarf Spheroidal Galaxies with LHAASO Observations

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes… ▽ More In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes of astrophysical $γ$-ray background while large amount of dark matter. By analyzing more than 700 days observational data at LHAASO, no significant dark matter signal from 1 TeV to 1 EeV is detected. Accordingly we derive the most stringent constraints on the ultra-heavy dark matter annihilation cross-section up to EeV. The constraints on the lifetime of dark matter in decay mode are also derived. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 17 pages, 12 figures, accepted by PRL

Data quality control system and long-term performance monitor of the LHAASO-KM2A

Authors: Zhen Cao, F. Aharonian, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, A. V. Bukevich, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, H. X. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen , et al. (263 additional authors not shown)

Abstract: The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To… ▽ More The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To ensure the reliability of the LHAASO-KM2A data, a three-level quality control system has been established. It is used to monitor the status of detector units, stability of reconstructed parameters and the performance of the array based on observations of the Crab Nebula and Moon shadow. This paper will introduce the control system and its application on the LHAASO-KM2A data collected from August 2021 to July 2023. During this period, the pointing and angular resolution of the array were stable. From the observations of the Moon shadow and Crab Nebula, the results achieved using the two methods are consistent with each other. According to the observation of the Crab Nebula at energies from 25 TeV to 100 TeV, the time averaged pointing errors are estimated to be $-0.003^{\circ} \pm 0.005^{\circ}$ and $0.001^{\circ} \pm 0.006^{\circ}$ in the R.A. and Dec directions, respectively. △ Less

Submitted 13 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures

Discovery of Very-high-energy Gamma-ray Emissions from the Low Luminosity AGN NGC 4278 by LHAASO

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) i… ▽ More The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) is compatible with NGC 4278 within $\sim0.03$ degree. Variation analysis shows an indication of the variability at a few months level in the TeV band, which is consistent with low frequency observations. Based on these observations, we report the detection of TeV $γ$-ray emissions from this low-luminosity AGN NGC 4278. The observations by LHAASO-WCDA during active period has a significance level of 8.8\,$σ$ with best-fit photon spectral index $\varGamma=2.56\pm0.14$ and a flux $f_{1-10\,\rm{TeV}}=(7.0\pm1.1_{\rm{sta}}\pm0.35_{\rm{syst}})\times10^{-13}\,\rm{photons\,cm^{-2}\,s^{-1}}$, or approximately $5\%$ of the Crab Nebula. The discovery of VHE from NGC 4278 indicates that the compact, weak radio jet can efficiently accelerate particles and emit TeV photons. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 11 pages, 5 figures

LHAASO-KM2A detector simulation using Geant4

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (254 additional authors not shown)

Abstract: KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with… ▽ More KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with large altitude difference (30 m) and huge coverage (1.3 km^2). In this paper, the design of the KM2A simulation code G4KM2A based on Geant4 is introduced. The process of G4KM2A is optimized mainly in memory consumption to avoid memory overffow. Some simpliffcations are used to signiffcantly speed up the execution of G4KM2A. The running time is reduced by at least 30 times compared to full detector simulation. The particle distributions and the core/angle resolution comparison between simulation and experimental data of the full KM2A array are also presented, which show good agreement. △ Less

Submitted 7 April, 2024; originally announced April 2024.

Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (256 additional authors not shown)

Abstract: We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at… ▽ More We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at $3.67 \pm 0.05 \pm 0.15$ PeV. Below the knee, the spectral index is found to be -$2.7413 \pm 0.0004 \pm 0.0050$, while above the knee, it is -$3.128 \pm 0.005 \pm 0.027$, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -$0.1200 \pm 0.0003 \pm 0.0341$. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components. △ Less

Submitted 26 March, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

Comments: 8 pages, 3 figures

Journal ref: Physical Review Letters 132, 131002 (2024)

Does or did the supernova remnant Cassiopeia A operate as a PeVatron?

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE;… ▽ More For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE; $E_γ\geq 100$~TeV) $γ$-rays. In this context, the historical SNR Cassiopeia A (Cas A) is considered one of the most promising target for UHE observations. This paper presents the observation of Cas A and its vicinity by the LHAASO KM2A detector. The exceptional sensitivity of LHAASO KM2A in the UHE band, combined with the young age of Cas A, enabled us to derive stringent model-independent limits on the energy budget of UHE protons and nuclei accelerated by Cas A at any epoch after the explosion. The results challenge the prevailing paradigm that Cas A-type SNRs are major suppliers of PeV CRs in the Milky Way. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 11 pages, 3 figures, Accepted by the APJL

Very high energy gamma-ray emission beyond 10 TeV from GRB 221009A

Authors: Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the t… ▽ More The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the trigger. The intrinsic energy spectrum of gamma-rays can be described by a power-law after correcting for extragalactic background light (EBL) absorption. Such a hard spectrum challenges the synchrotron self-Compton (SSC) scenario of relativistic electrons for the afterglow emission above several TeV. Observations of gamma-rays up to 13 TeV from a source with a measured redshift of z=0.151 hints more transparency in intergalactic space than previously expected. Alternatively, one may invoke new physics such as Lorentz Invariance Violation (LIV) or an axion origin of very high energy (VHE) signals. △ Less

Submitted 22 November, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 49pages, 11figures

Journal ref: Science Advances, 9, eadj2778 (2023) 15 November 2023

The First LHAASO Catalog of Gamma-Ray Sources

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022.… ▽ More We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022. This catalog represents the main result from the most sensitive large coverage gamma-ray survey of the sky above 1 TeV, covering declination from $-$20$^{\circ}$ to 80$^{\circ}$. In total, the catalog contains 90 sources with an extended size smaller than $2^\circ$ and a significance of detection at $> 5σ$. Based on our source association criteria, 32 new TeV sources are proposed in this study. Among the 90 sources, 43 sources are detected with ultra-high energy ($E > 100$ TeV) emission at $> 4σ$ significance level. We provide the position, extension, and spectral characteristics of all the sources in this catalog. △ Less

Submitted 27 November, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 40 pages, 13 figures, 4 tables

Journal ref: The Astrophysical Journal Supplement Series, 271 (2024) 25

Measurement of ultra-high-energy diffuse gamma-ray emission of the Galactic plane from 10 TeV to 1 PeV with LHAASO-KM2A

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer ar… ▽ More The diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer array of the Large High Altitude Air Shower Observatory (LHAASO). Diffuse emissions from the inner ($15^{\circ}<l<125^{\circ}$, $|b|<5^{\circ}$) and outer ($125^{\circ}<l<235^{\circ}$, $|b|<5^{\circ}$) Galactic plane are detected with $29.1σ$ and $12.7σ$ significance, respectively. The outer Galactic plane diffuse emission is detected for the first time in the very- to ultra-high-energy domain ($E>10$~TeV). The energy spectrum in the inner Galaxy regions can be described by a power-law function with an index of $-2.99\pm0.04$, which is different from the curved spectrum as expected from hadronic interactions between locally measured cosmic rays and the line-of-sight integrated gas content. Furthermore, the measured flux is higher by a factor of $\sim3$ than the prediction. A similar spectrum with an index of $-2.99\pm0.07$ is found in the outer Galaxy region, and the absolute flux for $10\lesssim E\lesssim60$ TeV is again higher than the prediction for hadronic cosmic ray interactions. The latitude distributions of the diffuse emission are consistent with the gas distribution, while the longitude distributions show clear deviation from the gas distribution. The LHAASO measurements imply that either additional emission sources exist or cosmic ray intensities have spatial variations. △ Less

Submitted 19 August, 2023; v1 submitted 9 May, 2023; originally announced May 2023.

Comments: 12 pages, 8 figures, 5 tables; accepted for publication in Physical Review Letters; source mask file provided as ancillary file

Journal ref: Phys. Rev. Lett. 131, 151001 (2023)

Features of the gamma-ray pulsar halo HESS J1831$-$098

Authors: Kun Fang, Shao-Qiang Xi, Li-Zhuo Bao, Xiao-Jun Bi, En-Sheng Chen

Abstract: Gamma-ray pulsar halos are ideal indicators of cosmic-ray propagation in localized regions of the Galaxy and electron injection from pulsar wind nebulae. HESS~J1831$-$098 is a candidate pulsar halo observed by both H.E.S.S. and HAWC experiments. We adopt the flux map of the H.E.S.S. Galactic plane survey and the spectrum measurements of H.E.S.S. and \textit{Fermi}-LAT to study HESS~J1831$-$098. We… ▽ More Gamma-ray pulsar halos are ideal indicators of cosmic-ray propagation in localized regions of the Galaxy and electron injection from pulsar wind nebulae. HESS~J1831$-$098 is a candidate pulsar halo observed by both H.E.S.S. and HAWC experiments. We adopt the flux map of the H.E.S.S. Galactic plane survey and the spectrum measurements of H.E.S.S. and \textit{Fermi}-LAT to study HESS~J1831$-$098. We find that HESS~J1831$-$098 meets all the criteria for a pulsar halo. The diffusion coefficient inside the halo and the conversion efficiency from the pulsar spin-down energy to the electron energy are both similar to the Geminga halo, a canonical pulsar halo. The injection spectrum can be well described by an exponentially-cutoff power law. However, the needed power-law term is very hard with $p\lesssim1$ if the diffusion coefficient is spatially and temporally independent. Considering the possible origins of the slow-diffusion environment, we adopt the two-zone diffusion model and the time-delayed slow-diffusion model. Both the models can interpret the H.E.S.S. and \textit{Fermi}-LAT results with a milder $p$. A modified injection time profile may have a similar effect. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 15 pages (one column), 4 figures, 1 table

Flux Variations of Cosmic Ray Air Showers Detected by LHAASO-KM2A During a Thunderstorm on 10 June 2021

Authors: LHAASO Collaboration, F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Zhe Cao, Zhen Cao, J. Chang, J. F. Chang, E. S. Chen, Liang Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, S. H. Chen, S. Z. Chen, T. L. Chen, X. J. Chen , et al. (248 additional authors not shown)

Abstract: The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations… ▽ More The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations of trigger rates (increases or decreases) are found to be strongly dependent on the primary zenith angle. The flux of secondary particles increases significantly, following a similar trend with that of the shower events. To better understand the observed behavior, Monte Carlo simulations are performed with CORSIKA and G4KM2A (a code based on GEANT4). We find that the experimental data (in saturated negative fields) are in good agreement with simulations, assuming the presence of a uniform upward electric field of 700 V/cm with a thickness of 1500 m in the atmosphere above the observation level. Due to the acceleration/deceleration and deflection by the atmospheric electric field, the number of secondary particles with energy above the detector threshold is modified, resulting in the changes in shower detection rate. △ Less

Submitted 6 December, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 18 pages, 11 figures

Journal ref: Chinese Phys. C 47 015001 (2023)

Peta-electron volt gamma-ray emission from the Crab Nebula

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Zhe Cao, J. Chang, J. F. Chang, B. M. Chen, E. S. Chen, J. Chen, Liang Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen , et al. (250 additional authors not shown)

Abstract: The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ pet… ▽ More The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ petaelectronvolt (PeV). The ultra-high-energy photons exhibit the presence of a PeV electron accelerator (a pevatron) with an acceleration rate exceeding 15% of the absolute theoretical limit. Assuming that unpulsed $γ$-rays are produced at the termination of the pulsar's wind, we constrain the pevatron's size, between $0.025$ and $0.1$ pc, and the magnetic field $\approx 110 μ$G. The production rate of PeV electrons, $2.5 \times 10^{36}$ erg $\rm s^{-1}$, constitutes 0.5% of the pulsar's spin-down luminosity, although we do not exclude a non-negligible contribution of PeV protons to the production of the highest energy $γ$-rays. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 43 pages, 13 figures, 2 tables; Published in Science

Journal ref: Science, 2021, Vol 373, Issue 6553, pp. 425-430

Fermi-LAT detection of extended gamma-ray emission in the vicinity of SNR G045.7$-$00.4: evidence for escaping cosmic rays interacting with the surrounding molecular clouds

Authors: Hai-Ming Zhang, Ruo-Yu Liu, Yang Su, Hui Zhu, Shao-Qiang Xi, Xiang-Yu Wang

Abstract: We present the analysis of Fermi Large Area Telecope (LAT) data of the gamma-ray emission in the vicinity of a radio supernova remnant (SNR), G045.7$-$00.4. To study the origin of the gamma-ray emission, we also make use of the CO survey data of Milky Way Imaging Scroll Painting to study the massive molecular gas complex that surrounds the SNR. The whole size of the GeV emission is significantly l… ▽ More We present the analysis of Fermi Large Area Telecope (LAT) data of the gamma-ray emission in the vicinity of a radio supernova remnant (SNR), G045.7$-$00.4. To study the origin of the gamma-ray emission, we also make use of the CO survey data of Milky Way Imaging Scroll Painting to study the massive molecular gas complex that surrounds the SNR. The whole size of the GeV emission is significantly larger than that of the radio morphology. Above 3 GeV, the GeV emission is resolved into two sources: one is spatially consistent with the position of the SNR with a size comparable to that of the radio emission, and the other is located outside of the western boundary of the SNR and spatially coincident with the densest region of the surrounding molecular cloud. We suggest that the GeV emission of the western source may arise from cosmic rays (CRs) which have escaped the SNR and illuminated the surrounding molecular cloud. We find that the gamma-ray spectra of the western source can be consistently explained by this scenario with a total energy of $\sim 10^{50}{\rm erg}$ in escaping CRs assuming the escape is isotropic. △ Less

Submitted 15 November, 2021; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: 11 pages, 7 figures and 1 table -Accepted for publication in ApJ

The physical properties of $γ$-ray quiet flat-spectrum radio quasars: why are they undetected by $Fermi$-LAT?

Authors: Xue-Jiao Deng, Rui Xue, Ze-Rui Wang, Shao-Qiang Xi, Hu-Bing Xiao, Lei-Ming Du, Zhao-Hua Xie

Abstract: During a decade of the $Fermi$-Large Area Telescope (LAT) operation, thousands of blazars have been detected in the $γ$-ray band. However, there are still numbers of blazars that have not been detected in the $γ$-ray band. In this work, we focus on investigating why some flat-spectrum radio quasars (FSRQs) are undetected by $Fermi$-LAT. By cross-matching the Candidate Gamma-ray Blazars Survey cata… ▽ More During a decade of the $Fermi$-Large Area Telescope (LAT) operation, thousands of blazars have been detected in the $γ$-ray band. However, there are still numbers of blazars that have not been detected in the $γ$-ray band. In this work, we focus on investigating why some flat-spectrum radio quasars (FSRQs) are undetected by $Fermi$-LAT. By cross-matching the Candidate Gamma-ray Blazars Survey catalog with the Fourth Catalog of Active Galactic Nuclei Detected by the $Fermi$-LAT, we select 11 $γ$-ray undetected ($γ$-ray quiet) FSRQs as our sample whose quasi-simultaneous multi-wavelength data are collected. In the framework of the conventional one-zone leptonic model, we investigate their underlying physical properties and study the possibility that they are undetected with $γ$-ray by modeling their quasi-simultaneous spectral energy distributions. In contrast to a smaller bulk Lorentz factor suggested by previous works, our results suggest that the dissipation region located relatively far away from the central super-massive black hole is more likely to be the cause of some $γ$-ray quiet FSRQs being undetected by $Fermi$-LAT. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: Accepted for publication by MNRAS. Comments welcome! (11 pages, 2 figures, 3 tables)

Self-consistent interpretations of the multi-wavelength gamma-ray spectrum of LHAASO J0621$+$3755

Authors: Kun Fang, Shao-Qiang Xi, Xiao-Jun Bi

Abstract: LHAASO J0621$+$3755 is a TeV gamma-ray halo newly identified by LHAASO-KM2A. It is likely to be generated by electrons trapped in a slow-diffusion zone around PSR J0622$+$3749 through inverse Compton scattering. When the gamma-ray spectrum of LHAASO-KM2A is fitted, the GeV fluxes derived by the commonly used one-zone normal diffusion model for electron propagation are significantly higher than the… ▽ More LHAASO J0621$+$3755 is a TeV gamma-ray halo newly identified by LHAASO-KM2A. It is likely to be generated by electrons trapped in a slow-diffusion zone around PSR J0622$+$3749 through inverse Compton scattering. When the gamma-ray spectrum of LHAASO-KM2A is fitted, the GeV fluxes derived by the commonly used one-zone normal diffusion model for electron propagation are significantly higher than the upper limits (ULs) of Fermi-LAT. In this work, we respectively adopt the one-zone superdiffusion and two-zone normal diffusion models to solve this conflict. For the superdiffusion scenario, we find that a model with superdiffusion index $α\lesssim1.2$ can meet the constraints of Fermi-LAT observation. For the two-zone diffusion scenario, the size of the slow-diffusion zone is required to be smaller than $\sim50$ pc, which is consistent with theoretical expectations. Future precise measurements of the Geminga halo may further distinguish between these two scenarios for the electron propagation in pulsar halos. △ Less

Submitted 5 July, 2021; originally announced July 2021.

Comments: 13 pages (one column), 4 figures, 1 table

Calibration of the Air Shower Energy Scale of the Water and Air Cherenkov Techniques in the LHAASO experiment

Authors: F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Z. Cao Z. Cao, J. Chang, J. F. Chang, X. C. Chang, B. M. Chen, J. Chen, L. Chen, L. Chen, L. Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (233 additional authors not shown)

Abstract: The Wide Field-of-View Cherenkov Telescope Array (WFCTA) and the Water Cherenkov Detector Arrays (WCDA) of LHAASO are designed to work in combination for measuring the energy spectra of various cosmic ray species over a very wide energy range from a few TeV to 10 PeV. The energy calibration of WCDA can be achieved with a proven technique of measuring the westward shift of the Moon shadow of galact… ▽ More The Wide Field-of-View Cherenkov Telescope Array (WFCTA) and the Water Cherenkov Detector Arrays (WCDA) of LHAASO are designed to work in combination for measuring the energy spectra of various cosmic ray species over a very wide energy range from a few TeV to 10 PeV. The energy calibration of WCDA can be achieved with a proven technique of measuring the westward shift of the Moon shadow of galactic cosmic rays due to the geomagnetic field. This deflection angle $Δ$ is inversely proportional to the energy of the cosmic rays. The precise measurements of the shifts by WCDA allows us to calibrate its energy scale for energies as high as 35 TeV. The energy scale measured by WCDA can be used to cross calibrate the energy reconstructed by WFCTA, which spans the whole energy range up to 10 PeV. In this work, we will demonstrate the feasibility of the method using the data collected from April 2019 to January 2020 by the WFCTA array and WCDA-1 detector, the first of the three water Cherenkov ponds, already commissioned at LHAASO site. △ Less

Submitted 13 April, 2021; v1 submitted 11 April, 2021; originally announced April 2021.

Journal ref: Phys. Rev. D 104, 062007 (2021)

Construction and On-site Performance of the LHAASO WFCTA Camera

Authors: F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Z. Cao, Z. Cao, J. Chang, J. F. Chang, X. C. Chang, B. M. Chen, J. Chen, L. Chen, L. Chen, L. Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (234 additional authors not shown)

Abstract: The focal plane camera is the core component of the Wide Field-of-view Cherenkov/fluorescence Telescope Array (WFCTA) of the Large High-Altitude Air Shower Observatory (LHAASO). Because of the capability of working under moonlight without aging, silicon photomultipliers (SiPM) have been proven to be not only an alternative but also an improvement to conventional photomultiplier tubes (PMT) in this… ▽ More The focal plane camera is the core component of the Wide Field-of-view Cherenkov/fluorescence Telescope Array (WFCTA) of the Large High-Altitude Air Shower Observatory (LHAASO). Because of the capability of working under moonlight without aging, silicon photomultipliers (SiPM) have been proven to be not only an alternative but also an improvement to conventional photomultiplier tubes (PMT) in this application. Eighteen SiPM-based cameras with square light funnels have been built for WFCTA. The telescopes have collected more than 100 million cosmic ray events and preliminary results indicate that these cameras are capable of working under moonlight. The characteristics of the light funnels and SiPMs pose challenges (e.g. dynamic range, dark count rate, assembly techniques). In this paper, we present the design features, manufacturing techniques and performances of these cameras. Finally, the test facilities, the test methods and results of SiPMs in the cameras are reported here. △ Less

Submitted 4 July, 2021; v1 submitted 29 December, 2020; originally announced December 2020.

Comments: 45 pages, 21 figures, article

Journal ref: Eur. Phys. J. C 81, 657 (2021)

Neutrino Emission from an Off-Axis Jet Driven by the Tidal Disruption Event AT2019dsg

Authors: Ruo-Yu Liu, Shao-Qiang Xi, Xiang-Yu Wang

Abstract: Recently, a high-energy muon neutrino event was detected in association with a tidal disruption event (TDE) AT2019dsg at the time about 150 days after the peak of the optical/UV luminosity. We propose that such a association could be interpreted as arising from hadronic interactions between relativistic protons accelerated in the jet launched from the TDE and the intense radiation field of TDE ins… ▽ More Recently, a high-energy muon neutrino event was detected in association with a tidal disruption event (TDE) AT2019dsg at the time about 150 days after the peak of the optical/UV luminosity. We propose that such a association could be interpreted as arising from hadronic interactions between relativistic protons accelerated in the jet launched from the TDE and the intense radiation field of TDE inside the optical/UV photosphere, if we are observing the jet at a moderate angle (i.e., approximately 10-30 degree) with respect to the jet axis. Such an off-axis viewing angle leads to a high gas column density in the line of sight which provides a high opacity for the photoionization and the Bethe-Heitler process, {and allows the existence of an intrinsic long-term X-ray radiation of comparatively high emissivity}. As a result, the cascade emission accompanying the neutrino production, which would otherwise overshoot the flux limits in X-ray and/or GeV band, is significantly obscured or absorbed. Since the jets of TDEs are supposed to be randomly oriented in the sky, the source density rate of TDE with an off-axis jet is significantly higher than that of TDE with an on-axis jet. Therefore, an off-axis jet is naturally expected in a nearby TDE being discovered, supporting the proposed scenario. △ Less

Submitted 7 November, 2020; originally announced November 2020.

Comments: 8 pages, 4 figures, published in PRD, 102, 083028

The observation of the Crab Nebula with LHAASO-KM2A for the performance study

Authors: F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Z. Cao, Z. Cao, J. Chang, J. F. Chang, X. C. Chang, B. M. Chen, J. Chen, L. Chen, L. Chen, L. Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (234 additional authors not shown)

Abstract: As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the… ▽ More As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the first observation on the Crab Nebula, a standard candle in very high energy gamma-ray astronomy. We detect gamma-ray signals from the Crab Nebula in both energy ranges of 10$-$100 TeV and $>$100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance including angular resolution, pointing accuracy and cosmic ray background rejection power. The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE =(1.13$\pm$0.05$_{stat}$$\pm$0.08$_{sys}$)$\times$10$^{-14}$$\cdot$(E/20TeV)$^{-3.09\pm0.06_{stat}\pm0.02_{sys}}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$. It is consistent with previous measurements by other experiments. This opens a new window of gamma-ray astronomy above 0.1 PeV through which ultrahigh-energy gamma-ray new phenomena, such as cosmic PeVatrons, might be discovered. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: 13 pages, 15 figures,submitted to CPC

Journal ref: Chinese Phys. C 45 025002 (2021)

Limiting Superluminal Neutrino Velocity and Lorentz Invariance Violation by Neutrino Emission from the Blazar TXS 0506+056

Authors: Kai Wang, Shao-Qiang Xi, Lijing Shao, Ruo-Yu Liu, Zhuo Li, Zhong-Kai Zhang

Abstract: The detection of high-energy neutrino coincident with the blazar TXS 0506+056 provides a unique opportunity to test Lorentz invariance violation (LIV) in the neutrino sector. Thanks to the precisely measured redshift, i.e., $z=0.3365$, the comoving distance of the neutrino source is determined. In this work, we obtain and discuss the constraints on the superluminal neutrino velocity $δ_ν$ and the… ▽ More The detection of high-energy neutrino coincident with the blazar TXS 0506+056 provides a unique opportunity to test Lorentz invariance violation (LIV) in the neutrino sector. Thanks to the precisely measured redshift, i.e., $z=0.3365$, the comoving distance of the neutrino source is determined. In this work, we obtain and discuss the constraints on the superluminal neutrino velocity $δ_ν$ and the LIV by considering the energy loss of superluminal neutrino during propagation. Given superluminal electron velocity ($δ_e \ge 0$), a very stringent constraint on superluminal neutrino velocity can be reached, i.e., $δ_ν\lesssim 1.3\times 10^{-18}$, corresponding to the quantum gravity (QG) scale $M_{\rm QG,1} \gtrsim 5.7 \times 10^{3} M_{\rm Pl}$ and $M_{\rm QG,2} \gtrsim 9.3 \times 10^{-6} M_{\rm Pl}$ for linear (quadratic) LIV, which are $\sim 12$ orders of magnitude tighter for linear LIV and $\sim 9$ orders tighter for quadratic LIV compared to the time-of-flight constraint from MeV neutrinos of SN 1987A. While given the subluminal electron velocity, a weaker constraint on the superluminal neutrino velocity is obtained, i.e., $δ_ν\lesssim 8 \times 10^{-17}$, which is consistent with the conclusions of previous works. We also study the neutrino detection probability due to the distortion of neutrino spectral shape during propagation, which gives slightly weaker constraints than above by a factor of $\sim2$. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Comments: 4 pages, 1 table, accepted by Phys. Rev. D

Journal ref: Phys. Rev. D 102 (2020) 063027

The physical properties of $Fermi$-4LAC flat spectrum radio quasars

Authors: Can Tan, Rui Xue, Lei-Ming Du, Shao-Qiang Xi, Ze-Rui Wang, Zhao-Hua Xie

Abstract: In this work, we collect quasi-simultaneous infrared, optical, X-ray and $γ$-ray data of 60 $Fermi$-4LAC flat spectrum radio quasars (FSRQs). In the framework of the conventional one-zone leptonic model, we investigate the physical properties of $Fermi$-4LAC FSRQs' jets by modeling their quasi-simultaneous spectral energy distributions (SEDs). Our main results are summarized as follows. (1) There… ▽ More In this work, we collect quasi-simultaneous infrared, optical, X-ray and $γ$-ray data of 60 $Fermi$-4LAC flat spectrum radio quasars (FSRQs). In the framework of the conventional one-zone leptonic model, we investigate the physical properties of $Fermi$-4LAC FSRQs' jets by modeling their quasi-simultaneous spectral energy distributions (SEDs). Our main results are summarized as follows. (1) There is a linear correlation between synchrotron peak frequency and curvature of the electron energy distribution. As suggested by previous works, the slope of the best linear fitting equation of this correlation is consistent with statistic acceleration which needs a fluctuation of fractional acceleration gain. (2) The gamma-ray dissipation regions are located at the range from 0.1 to 10 pc away from the super-massive black hole, and located outside the broad-line region (BLR) and within the dusty torus (DT). (3) A size relation $P_{\rm e}$ (the kinetic power carried in relativistic electrons) $\sim$ $P_{\rm B}$ (Poynting flux) $\leq$ $P_{\rm r}$ (the radiative power ) $<$ $P_{\rm p}$ (the kinetic power in cold protons) is found in our modeling. Among them, $P_{\rm e}\sim P_{\rm B}$ suggests that SEDs of almost all FSRQs with parameters are close to equipartition between the magnetic field and the relativistic electrons. The $P_{\rm e} < P_{\rm r}$ suggest that the most energy of the relativistic electrons are dissipated by EC radiation for FSRQs. (4) There is an anti-correlation between the peak energy of SEDs ($γ_{\rm peak}$) and the jet power ($P_{\rm jet}$), which is consistent with the blazar sequence. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: Accepted for publication in ApJS

GeV $γ$-ray emission from M33 and Arp~299

Authors: Shao-Qiang Xi, Hai-Ming Zhang, Ruo-Yu Liu, Xiang-Yu Wang

Abstract: Star-forming galaxies are huge reservoirs of cosmic rays (CRs) and these CRs convert a significant fraction of their energy into $γ$-rays by colliding with the interstellar medium (ISM). Several nearby star-forming galaxies have been detected in GeV-TeV $γ$-rays. It is also found that the $γ$-ray luminosities in 0.1-100 GeV correlate well with indicators of star formation rates of the galaxies, su… ▽ More Star-forming galaxies are huge reservoirs of cosmic rays (CRs) and these CRs convert a significant fraction of their energy into $γ$-rays by colliding with the interstellar medium (ISM). Several nearby star-forming galaxies have been detected in GeV-TeV $γ$-rays. It is also found that the $γ$-ray luminosities in 0.1-100 GeV correlate well with indicators of star formation rates of the galaxies, such as the total infrared (IR) luminosity. In this paper, we report a systematic search for possible $γ$-ray emission from galaxies in the IRAS Revised Bright Galaxies Sample, using 11.4 years of $γ$-ray data taken by the Fermi Large Area Telescope (LAT). Two new galaxies, M33 and Arp 299, are detected significantly. The two galaxies are consistent with the empirical correlation between the $γ$-ray luminosity and total infrared luminosity, suggesting that their $γ$-ray emissions should mainly originate from CRs interacting with ISM. Nevertheless, there is a tentative evidence that the flux of the $γ$-ray emission from Arp~299 is variable. If the variability is true, part of the emission from Arp 299 should originate from the obscured AGN in this interacting galaxy system. In addition, we find that the $γ$-ray excess from M33 is located at the northeast region of the galaxy, where a supergiant H II region, NGC604, resides. This indicates that some bright star-forming regions in spiral galaxies could play a dominant role in the galaxy in producing $γ$-ray emission. △ Less

Submitted 1 July, 2020; v1 submitted 17 March, 2020; originally announced March 2020.

Comments: 9 pages, 9 figures, 2 tables; ApJ accepted

A serendipitous discovery of GeV gamma-ray emission from supernova 2004dj in a survey of nearby star-forming galaxies with Fermi-LAT

Authors: Shao-Qiang Xi, Ruo-Yu Liu, Xiang-Yu Wang, Rui-Zhi Yang, Qiang Yuan, Bing Zhang

Abstract: The interaction between a supernova ejecta and the circum-stellar medium drives a strong shock wave which accelerates particles (i.e., electrons and protons). The radio and X-ray emission observed after the supernova explosion constitutes the evidence of the electron acceleration. The accelerated protons are expected to produce GeV-TeV gamma-ray emission via $pp$ collisions, but the flux is usuall… ▽ More The interaction between a supernova ejecta and the circum-stellar medium drives a strong shock wave which accelerates particles (i.e., electrons and protons). The radio and X-ray emission observed after the supernova explosion constitutes the evidence of the electron acceleration. The accelerated protons are expected to produce GeV-TeV gamma-ray emission via $pp$ collisions, but the flux is usually low since only a small fraction of the supernova kinetic energy is converted into the shock energy at the very early time. The low gamma-ray flux of the nearest supernova explosion, SN 1987A, agrees with this picture. Here we report a serendipitous discovery of a fading GeV gamma-ray source in spatial coincidence with the second nearest supernova--SN 2004dj from our gamma-ray survey of nearby star-forming galaxies with Fermi-LAT. The total gamma-ray energy released by SN 2004dj is about $6\times10^{47}{\rm erg}$. We interpret this gamma-ray emission arising from the supernova ejecta interacting with a surrounding high-density shell, which decelerates the ejecta and converts ~1% of the ejecta's kinetic energy to relativistic protons. In addition, our gamma-ray survey of nearby star-forming galaxies discovers GeV emissions from two star-forming galaxies, i.e., Arp 299 and M33, for the first time. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: 8 pages, 5 figures, 1 table, comments are welcome

Constraints on the intergalactic magnetic field from $γ$-ray observations of GRB 190114C

Authors: Ze-Rui Wang, Shao-Qiang Xi, Ruo-Yu Liu, Rui Xue, Xiang-Yu Wang

Abstract: Very high energy photons from cosmological gamma-ray bursts (GRBs) are expected to interact with extragalactic background light (EBL) and produce electron-positron pairs when they propagate through intergalactic medium (IGM). These relativistic pairs will then up-scatter cosmic microwave background (CMB) photons and emit secondary GeV emission. Meanwhile, the motion of these pairs are deflected by… ▽ More Very high energy photons from cosmological gamma-ray bursts (GRBs) are expected to interact with extragalactic background light (EBL) and produce electron-positron pairs when they propagate through intergalactic medium (IGM). These relativistic pairs will then up-scatter cosmic microwave background (CMB) photons and emit secondary GeV emission. Meanwhile, the motion of these pairs are deflected by intergalactic magnetic field (IGMF), so the secondary GeV photons arrive later than the primary emission. It has been suggested that the properties of the secondary GeV emission can be used to constrain IGMF. Recently, TeV gamma-ray emission has been detected, for the first time, from a GRB (GRB 190114C) by the MAGIC telescope and its steep ${\rm γ-ray}$ spectrum shows a clear evidence of absorption by EBL. We then constrain the IGMF with the GeV flux limit obtained from the $Fermi$-LAT observations. We find a limit of $>10^{-19.5}$ G for the coherence length of $λ\leq 1$ Mpc. Although this limit is weaker than that obtained by using blazars, it represents the first limit from ${\rm γ-ray}$ observations of GRBs, which provides an independent constraint on IGMF. We also find that, for transient ${\rm γ-ray}$ sources, one can choose a favorable time window to search for the echo emission at a particular energy. △ Less

Submitted 19 March, 2020; v1 submitted 5 January, 2020; originally announced January 2020.

Comments: 6 pages, 5 figures, Accepted by PRD

Journal ref: Phys. Rev. D 101, 083004 (2020)

Discovery of a spatially extended GeV source in the vicinity of the TeV halo candidate 2HWC J1912+099: a TeV halo or supernova remnant ?

Authors: Hai-Ming Zhang, Shao-Qiang Xi, Ruo-Yu Liu, Yu-Liang Xin, Siming Liu, Xiang-Yu Wang

Abstract: Observations by HAWC and Milagro have detected spatially extended TeV sources surrounding middle-aged ($t\sim100-400 \,{\rm kyr}$) pulsars like Geminga and PSR B0656+14, which have been named "TeV Halos", representing very extended TeV pulsar wind nebulae (PWNe) powered by relatively old pulsars. A few more HAWC-detected sources have been suggested to be TeV halo candidates. In this paper, we sear… ▽ More Observations by HAWC and Milagro have detected spatially extended TeV sources surrounding middle-aged ($t\sim100-400 \,{\rm kyr}$) pulsars like Geminga and PSR B0656+14, which have been named "TeV Halos", representing very extended TeV pulsar wind nebulae (PWNe) powered by relatively old pulsars. A few more HAWC-detected sources have been suggested to be TeV halo candidates. In this paper, we search for possible GeV counterparts of three TeV halo candidates with Fermi Large Area Telescopes. We detect a new spatially extended GeV source in the vicinity of the TeV halo candidate 2HWC J1912+099, which is also detected by HESS (HESS J1912+101). We find that the size of the GeV source is significantly larger than that of the TeV emission measured by HESS, and a spatial template characteristic of a PWN can fit the GeV data. We suggest that the GeV source is an extended PWN powered by the central middle-aged pulsar PSR J1913+1011. This discovery favors the TeV halo scenario for the TeV source 2HWC J1912+099 (HESS J1912+101), although the possible shell-like morphology measured by HESS challenges this interpretation. Alternatively, the TeV emission could be dominated by a supernova remnant (SNR) via the hadronic process. Future multi-wavelength observations of the source and more precise measurements of the spatial profile of the TeV emission will be useful to distinguish between the two scenarios. △ Less

Submitted 22 November, 2019; v1 submitted 28 September, 2019; originally announced September 2019.

Comments: Accepted by ApJ, 7 pages, 6 figures

Synchrotron self-Compton emission from external shocks as the origin of the sub-TeV emission in GRB 180720B and GRB 190114C

Authors: Xiang-Yu Wang, Ruo-Yu Liu, Hai-Ming Zhang, Shao-Qiang Xi, Bing Zhang

Abstract: Recently, very high-energy photons above 100 GeV were reported to be detected from GRB~190114C and GRB~180720B at, respectively, 100-1000 s and 10 hours after the burst. We model the available broad-band data of both GRBs with the synchrotron plus synchrotron self-Compton (SSC) emission of the afterglow shocks. We find that the sub-TeV emission of GRB~180720B can be interpreted as the SSC emission… ▽ More Recently, very high-energy photons above 100 GeV were reported to be detected from GRB~190114C and GRB~180720B at, respectively, 100-1000 s and 10 hours after the burst. We model the available broad-band data of both GRBs with the synchrotron plus synchrotron self-Compton (SSC) emission of the afterglow shocks. We find that the sub-TeV emission of GRB~180720B can be interpreted as the SSC emission from afterglow shocks expanding in a constant density circum-burst medium. The SSC emission of GRB~190114C dominates over the synchrotron component from GeV energies at $\sim100$ s, which can explain the possible hard spectrum of the GeV emission at this time. The extrapolated flux of this SSC component to sub-TeV energies can explain the high-significance detection of GRB~190114C by the MAGIC telescope. The parameter values (such as the circum-burst density and shock microphysical parameters) in the modeling are not unusual for both GRBs, implying that the detection of sub-TeV photons from these two bursts should be attributed to their large burst energies and low redshifts. △ Less

Submitted 6 September, 2019; v1 submitted 27 May, 2019; originally announced May 2019.

Comments: 7 pages, 3 figures, accepted by ApJ

Tentative evidence of spatially extended GeV emission from SS433/W50

Authors: Xiao-Na Sun, Rui-Zhi Yang, Bing Liu, Shao-Qiang Xi, Xiang-Yu Wang

Abstract: We analyze 10 years of Fermi-LAT data towards the SS433/W50 region. With the latest source catalog and diffuse background models, the gamma-ray excess from SS433/W50 is detected with a significance of 6σ in the photon energy range of 500 MeV - 10 GeV. Our analysis indicates that an extended flat disk morphology is preferred over a point-source description, suggesting that the GeV emission region i… ▽ More We analyze 10 years of Fermi-LAT data towards the SS433/W50 region. With the latest source catalog and diffuse background models, the gamma-ray excess from SS433/W50 is detected with a significance of 6σ in the photon energy range of 500 MeV - 10 GeV. Our analysis indicates that an extended flat disk morphology is preferred over a point-source description, suggesting that the GeV emission region is much larger than that of the TeV emission detected by HAWC. The size of the GeV emission is instead consistent with the extent of the radio nebula W50, a supernova remnant being distorted by the jets, so we suggest that the GeV emission may originate from this supernova remnant. The spectral result of the GeV emission is also consistent with an supernova remnant origin. We also derive the GeV flux upper limits on the TeV emission region, which put moderate constrains on the leptonic models to explain the multiwavelength data. △ Less

Submitted 15 May, 2019; v1 submitted 10 April, 2019; originally announced April 2019.

Comments: 7 pages, 4 figures, accepted for publication in A&A

Journal ref: A&A 626, A113 (2019)

Comparative study of gamma-ray emission from molecular clouds and star-forming galaxies

Authors: Fang-Kun Peng, Shao-Qiang Xi, Xiang-Yu Wang, Qi-Jun Zhi, Di Li

Abstract: Star-forming regions on different scales, such as giant molecular clouds in our Galaxy and star-forming galaxies, emit GeV gamma-rays. These are thought to originate from hadronic interactions of cosmic-ray (CR) nuclei with the interstellar medium. It has recently been shown that the gamma-ray luminosity ($L_γ$) of star-forming galaxies is well correlated with their star formation rates (SFR). We… ▽ More Star-forming regions on different scales, such as giant molecular clouds in our Galaxy and star-forming galaxies, emit GeV gamma-rays. These are thought to originate from hadronic interactions of cosmic-ray (CR) nuclei with the interstellar medium. It has recently been shown that the gamma-ray luminosity ($L_γ$) of star-forming galaxies is well correlated with their star formation rates (SFR). We investigated \textsl{Fermi} data of eight Galactic molecular clouds in the Gould belt and found that molecular clouds do not follow the $L_γ-{\rm SFR}$ correlation of star-forming galaxies. We also compared the scaling relations of gamma-ray luminosity, SFR, and the gas mass for molecular clouds and star-forming galaxies. Using a multiple-variable regression analysis, we found different dependences of gamma-ray emission on SFR or mass for molecular clouds and star-forming galaxies. This suggests that different mechanisms may govern the production of gamma-rays in these two types of sources. Specifically, the strong dependence on mass supports that gamma-ray emission of molecular clouds primarily comes from {\em \textup{passive}} interaction by diffuse Galactic CRs, whereas the strong dependence on SFR supports that gamma-ray emission of star-forming galaxies originates from CRs that are accelerated by local {\em \textup{active}} sources. △ Less

Submitted 17 November, 2018; originally announced November 2018.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 621, A70 (2019)

GeV observations of the extended pulsar wind nebulae constrain the pulsar interpretations of the cosmic-ray positron excess

Authors: Shao-Qiang Xi, Ruo-Yu Liu, Zhi-Qiu Huang, Kun Fang, Xiang-Yu Wang

Abstract: It has long been suggested that nearby pulsars within $\sim 1 \,{\rm kpc}$ are the leading candidate of the 10-500 GeV cosmic-ray positron excess measured by PAMELA and other experiments. The recent measurement of surface brightness profile of TeV nebulae surrounding Geminga and PSR~B0656+14 by the High-Altitude Water Cherenkov Observatory (HAWC) suggests inefficient diffusion of particles from th… ▽ More It has long been suggested that nearby pulsars within $\sim 1 \,{\rm kpc}$ are the leading candidate of the 10-500 GeV cosmic-ray positron excess measured by PAMELA and other experiments. The recent measurement of surface brightness profile of TeV nebulae surrounding Geminga and PSR~B0656+14 by the High-Altitude Water Cherenkov Observatory (HAWC) suggests inefficient diffusion of particles from the sources, giving rise to a debate on the pulsar interpretation of the cosmic-ray positron excess. Here we argue that GeV observations provide more direct constraints on the positron density in the TeV nebulae in the energy range of 10-500 GeV and hence on the origin of the observed positron excess. Motivated by this, we search for GeV emission from the TeV nebulae with the \textsl{Fermi} Large Area Telescope (LAT). No spatially-extended GeV emission is detected from these two TeV nebulae in the framework of two-zone diffusion spatial templates, suggesting a relatively low density of GeV electrons/positrons in the TeV nebulae. A joint modelling of the data from HAWC and \textsl{Fermi}-LAT disfavors Geminga and PSR~B0656+14 as the dominant source of the positron excess at $\sim 50-500$ GeV for the usual Kolmogorov-type diffusion, while for an energy-independent diffusion, a dominant part of the positron excess contributed by them cannot be ruled out by the current data. △ Less

Submitted 9 May, 2019; v1 submitted 25 October, 2018; originally announced October 2018.

Comments: 26 pages, 11 figures, significantly expanded, to appear in ApJ

Is the gamma-ray emission around the Coma galaxy cluster a dark matter origin?

Authors: Yun-Feng Liang, Zi-Qing Xia, Shao-Qiang Xi, Shang Li, Zhao-Qiang Shen, Yi-Zhong Fan

Abstract: Recently, gamma-ray emission in the direction of Coma, with a TS value of $\sim 40$, has been reported. In this work we will discuss the possibility of such a residual emission coming from dark matter annihilation. Our results show that the gamma-ray emission within the Coma region is also spatially correlated to the mass distribution derived from weak gravitational lensing measurements very well.… ▽ More Recently, gamma-ray emission in the direction of Coma, with a TS value of $\sim 40$, has been reported. In this work we will discuss the possibility of such a residual emission coming from dark matter annihilation. Our results show that the gamma-ray emission within the Coma region is also spatially correlated to the mass distribution derived from weak gravitational lensing measurements very well. However the dark matter models are not supported by the spectral analysis results and constraints by observations of other targets. Thus we derive the upper limits of the dark matter annihilation cross section according to the observation of the Coma region. △ Less

Submitted 5 January, 2018; originally announced January 2018.

Comments: 8 pages, 4 figures, comments are welcome

Detection of gamma-ray emission from the Coma cluster with Fermi Large Area Telescope and tentative evidence for an extended spatial structure

Authors: Shao-Qiang Xi, Xiang-Yu Wang, Yun-Feng Liang, Fang-Kun Peng, Rui-Zhi Yang, Ruo-Yu Liu

Abstract: Many galaxy clusters have giant halos of non-thermal radio emission, indicating the presence of relativistic electrons in the clusters. Relativistic protons may also be accelerated by merger and/or accretion shocks in galaxy clusters. These cosmic-ray (CR) electrons and/or protons are expected to produce gamma-rays through inverse-Compton scatterings or inelastic $pp$ collisions respectively. Desp… ▽ More Many galaxy clusters have giant halos of non-thermal radio emission, indicating the presence of relativistic electrons in the clusters. Relativistic protons may also be accelerated by merger and/or accretion shocks in galaxy clusters. These cosmic-ray (CR) electrons and/or protons are expected to produce gamma-rays through inverse-Compton scatterings or inelastic $pp$ collisions respectively. Despite of intense efforts in searching for high-energy gamma-ray emission from galaxy clusters, conclusive evidence is still missing so far. Here we report the discovery of $\ge 200$ MeV gamma-ray emission from the Coma cluster direction with an unbinned likelihood analysis of the 9 years of {\it Fermi}-LAT Pass 8 data. The gamma-ray emission shows a spatial morphology roughly coincident with the giant radio halo, with an apparent excess at the southwest of the cluster. Using the test statistic analysis, we further find tentative evidence that the gamma-ray emission at the Coma center is spatially extended. The extended component has an integral energy flux of $\sim 2\times 10^{-12}{\rm \ erg\ cm^{-2}\ s^{-1}}$ in the energy range of 0.2 - 300 GeV and the spectrum is soft with a photon index of $\simeq-2.7$. Interpreting the gamma-ray emission as arising from CR proton interaction, we find that the volume-averaged value of the CR to thermal pressure ratio in the Coma cluster is about $\sim 2\%$. Our results show that galaxy clusters are likely a new type of GeV gamma-ray sources, and they are probably also giant reservoirs of CR protons. △ Less

Submitted 16 August, 2018; v1 submitted 25 September, 2017; originally announced September 2017.

Comments: 10 pages, 10 figures, Accepted by Physical Review D, more spatial models for the gamma-ray emission are used, systematic checks on the results are added

Journal ref: Phys. Rev. D 98, 063006 (2018)

A New Method to Test the Einstein's Weak Equivalence Principle

Authors: H. Yu, S. Q. Xi, F. Y. Wang

Abstract: The Einstein's weak equivalence principle (WEP) is one of the foundational assumptions of general relativity and some other gravity theories. In the theory of parametrized post-Newtonian (PPN), the difference between the PPN parameters $γ$ of different particles or the same type of particle with different energies, $Δγ$, represents the violation of WEP. Current constraints on $Δγ$ are derived from… ▽ More The Einstein's weak equivalence principle (WEP) is one of the foundational assumptions of general relativity and some other gravity theories. In the theory of parametrized post-Newtonian (PPN), the difference between the PPN parameters $γ$ of different particles or the same type of particle with different energies, $Δγ$, represents the violation of WEP. Current constraints on $Δγ$ are derived from the observed time delay between correlated particles of astronomical sources. However, the observed time delay is contaminated by other effects, such as the time delays due to different particle emission times, the potential Lorentz invariance violation and none-zero photon rest mass. Therefore, current constraints are only upper limits. Here we propose a new method to test WEP based on the fact that the gravitational time delay is direction-dependent while others are not. This is the first method which can naturally correct other time delay effects. Using the time delay measurements of BASTE gamma-ray burst (GRB) sample and the gravitational potential of local super galaxy cluster Laniakea, we find that the constraint on $Δγ$ of different energy photons can be as low as $10^{-14}$. In the future, if more gravitational wave events and fast radio bursts with much more precise time delay measurements are observed, this method can give a reliable and tight constraint on WEP. △ Less

Submitted 4 May, 2018; v1 submitted 8 August, 2017; originally announced August 2017.

Comments: 27 pages, 4 figures, 4 tables, accepted for publication in ApJ

Search for GeV counterparts to fast radio bursts with Fermi

Authors: Shao-Qiang Xi, Pak-Hin Thomas Tam, Fang-Kun Peng, Xiang-Yu Wang

Abstract: The non-repeating fast radio bursts (FRBs) could arise from catastrophic stellar explosions or magnetar giant flares, so relativistic blast waves might be produced in these events. Motivated by this, we here search for GeV counterparts to all non-repeating FRBs with Fermi Large Area Telescope (LAT), including FRB 131104 that is claimed to be possibly associated with a gamma-ray transient candidate… ▽ More The non-repeating fast radio bursts (FRBs) could arise from catastrophic stellar explosions or magnetar giant flares, so relativistic blast waves might be produced in these events. Motivated by this, we here search for GeV counterparts to all non-repeating FRBs with Fermi Large Area Telescope (LAT), including FRB 131104 that is claimed to be possibly associated with a gamma-ray transient candidate detected by Swift Burst Alert Telescope (BAT). FRB 131104 enters the field of view (FoV) of LAT ~5000s after the burst time, so we are only able to search for the GeV afterglow emission during this period, but no significant GeV emission is detected. We also perform a search for GeV emissions from other FRBs, but no significant GeV emissions are detected either. Upper limit fluences in the range of (4.7-29.2)x10^(-7) erg/cm2 are obtained, and then the upper limits of the isotropic blast wave kinetic energy of about (1-200)x10^(53) erg are inferred under certain assumptions. Although the current limits on the isotropic blast wave energy are not sufficiently stringent to rule out the connection between FRBs and GRB-like transients, future more sensitive observations with Fermi or Imaging Atmospheric Cherenkov Telescopes might be able to constrain the connection. △ Less

Submitted 10 May, 2017; originally announced May 2017.

Comments: 12 pages, 2 figures, 2 tables, ApJ Letters, in press

The History of GRB Outflows: Ejection Lorentz Factor and Radiation Location of X-Ray Flares

Authors: Hui-Jun Mu, Da-Bin Lin, Shao-Qiang Xi, Ting-Ting Lin, Yuan-Zhu Wang, Yun-Feng Liang, Lian-Zhong Lv, Jin Zhang, En-Wei Liang

Abstract: We present time-resolved spectral analysis of the steep decay segments of 29 bright X-ray flares of gamma-ray bursts (GRBs) observed with the Swift/X-ray telescope, and model their lightcurves and spectral index evolution behaviors with the curvature effect model. Our results show that the observed rapid flux decay and strong spectral index evolution with time can be well fit with this model, and… ▽ More We present time-resolved spectral analysis of the steep decay segments of 29 bright X-ray flares of gamma-ray bursts (GRBs) observed with the Swift/X-ray telescope, and model their lightcurves and spectral index evolution behaviors with the curvature effect model. Our results show that the observed rapid flux decay and strong spectral index evolution with time can be well fit with this model, and the derived characteristic timescales ($t_c$) are in the range of $33\sim 264$ seconds. Using an empirical relation between the peak luminosity and the Lorentz factor derived from the prompt gamma-rays, we estimate the Lorentz factors of the flares ($Γ_{\rm X}$). We obtain $Γ_{\rm X}=17\sim 87$ with a median value of $52$, which is smaller than the initial Lorentz factors of prompt gamma-ray fireballs. With the derived $t_c$ and $Γ_{\rm X}$, we constrain the radiating regions of 13 X-ray flares, yielding $R_{\rm X}=(0.2\sim 1.1)\times 10^{16}$ cm, which are smaller than the radii of the afterglow fireballs at the peak times of the flares. A long evolution feature from prompt gamma-ray phase to the X-ray epoch is found by incorporating our results with a sample of GRBs whose initial Lorentz factors are available in literatures, i.e., $Γ\propto [t_{p}/(1+z)]^{-0.69\pm 0.06}$. These results may shed lights on the long term evolution of GRB central engines. △ Less

Submitted 22 December, 2016; v1 submitted 17 August, 2016; originally announced August 2016.

Comments: 28 pages, 33 figures, accepted for publication in ApJ

Comprehensive study of the X-ray flares from gamma-ray bursts observed by Swift

Authors: Shuang-Xi Yi, Shao-Qiang Xi, Hai Yu, F. Y. Wang, Hui-Jun Mu, Lian-Zhong Lv, En-Wei Liang

Abstract: X-ray flares are generally supposed to be produced by the later central engine activities, and may share the similar physical origin with prompt emission of gamma-ray bursts (GRBs). In this paper, we have analyzed all significant X-ray flares from the GRBs observed by {\em Swift} from April 2005 to March 2015. The catalog contains 468 bright X-ray flares, including 200 flares with redshifts. We ob… ▽ More X-ray flares are generally supposed to be produced by the later central engine activities, and may share the similar physical origin with prompt emission of gamma-ray bursts (GRBs). In this paper, we have analyzed all significant X-ray flares from the GRBs observed by {\em Swift} from April 2005 to March 2015. The catalog contains 468 bright X-ray flares, including 200 flares with redshifts. We obtain the fitting results of X-ray flares, such as start time, peak time, duration, peak flux, fluence, peak luminosity, and mean luminosity. The peak luminosity decreases with peak time, following a power-law behavior $L_p \propto T_{peak,z}^{-1.27}$. The flare duration increases with peak time. The 0.3-10 keV isotropic energy of X-ray flares distribution is a lognormal peaked at $10^{51.2}$ erg. We also study the frequency distributions of flare parameters, including energies, durations, peak fluxes, rise times, decay times and waiting times. Power-law distributions of energies, durations, peak fluxes, and waiting times are found in GRB X-ray flares and solar flares. These distributions could be well explained by a fractal-diffusive, self-organized criticality model. Some theoretical models basing on magnetic reconnection have been proposed to explain X-ray flares. Our result shows that the relativistic jets of GRBs may be Poynting-flux dominated. △ Less

Submitted 24 March, 2016; originally announced March 2016.

Comments: 36 pages, 10 figures, 2 tables, accepted for publication in ApJS

Universal Behavior of X-ray Flares from Black Hole Systems

Authors: F. Y. Wang, Z. G. Dai, S. X. Yi, S. Q. Xi

Abstract: X-ray flares have been discovered in black hole systems, such as gamma-ray bursts, the tidal disruption event Swift J1644+57, the supermassive black hole Sagittarius A$^*$ at the center of our Galaxy, and some active galactic nuclei. Their occurrences are always companied by relativistic jets. However, it is still unknown whether there is a physical analogy among such X-ray flares produced in blac… ▽ More X-ray flares have been discovered in black hole systems, such as gamma-ray bursts, the tidal disruption event Swift J1644+57, the supermassive black hole Sagittarius A$^*$ at the center of our Galaxy, and some active galactic nuclei. Their occurrences are always companied by relativistic jets. However, it is still unknown whether there is a physical analogy among such X-ray flares produced in black hole systems spanning nine orders of magnitude in mass. Here we report the observed data of X-ray flares, and show that they have three statistical properties similar to solar flares, including power-law distributions of energies, durations, and waiting times, which both can be explained by a fractal-diffusive self-organized criticality model. These statistical similarities, together with the fact that solar flares are triggered by a magnetic reconnection process, suggest that all of the X-ray flares are consistent with magnetic reconnection events, implying that their concomitant relativistic jets may be magnetically dominated. △ Less

Submitted 15 November, 2014; originally announced November 2014.

Comments: 23 pages, 4 tables, 6 figures. Accepted for publication by ApJS

Photosphere emission in the X-Ray Flares of Swift Gamma-Ray Bursts and Implications for the Fireball Properties

Authors: Fang-Kun Peng, En-Wei Liang, Xiang-Yu Wang, Shu-Jin Hou, Shao-Qiang Xi, Rui-Jing Lu, Jin Zhang, Bing Zhang

Abstract: X-ray flares of gamma-ray bursts (GRBs) are usually observed in the soft X-ray range and the spectral coverage is limited. In this paper, we present an analysis of 32 GRB X-ray flares that are simultaneously observed by both BAT and XRT on board the Swift mission, so a joint spectral analysis with a wider spectral coverage is possible. Our results show that the joint spectra of 19 flares are fitte… ▽ More X-ray flares of gamma-ray bursts (GRBs) are usually observed in the soft X-ray range and the spectral coverage is limited. In this paper, we present an analysis of 32 GRB X-ray flares that are simultaneously observed by both BAT and XRT on board the Swift mission, so a joint spectral analysis with a wider spectral coverage is possible. Our results show that the joint spectra of 19 flares are fitted with the absorbed single power-law or the Band function models. More interestingly, the joint spectra of the other 13 X-ray flares are fitted with the absorbed single power-law model plus a black body (BB) component. Phenomenally, the observed spectra of these 13 flares are analogous to several GRBs with a thermal component, but only with a much lower temperature of $kT=1\sim 3$ keV. Assuming that the thermal emission is the photosphere emission of the GRB fireball, we derive the fireball properties of the 13 flares that have redshift measurements, such as the bulk Lorentz factor $Γ_{\rm ph}$ of the outflow. The derived $Γ_{\rm ph}$ range from $50$ to $150$ and a relation of $Γ_{\rm ph}$ to the thermal emission luminosity is found. It is consistent with the $Γ_0-L_{\rm iso}$ relation that are derived for the prompt gamma-ray emission. We discuss the physical implications of these results within the content of jet composition and radiation mechanism of GRBs and X-ray flares. △ Less

Submitted 18 September, 2014; originally announced September 2014.

Comments: 31 pages, 8 figures, and 4 tables. Accepted for publication in ApJ

Internal Energy Dissipation of Gamma-Ray Bursts Observed with Swift: Precursors, Prompt Gamma-rays, Extended emission and Late X-ray Flares

Authors: You-Dong Hu, En-Wei Liang, Shao-Qiang Xi, Fang-Kun Peng, Rui-Jing Lu, Lian-Zhong LV, Bing Zhang

Abstract: We jointly analyze the gamma-ray burst (GRB) data observed with BAT and XRT on board the Swift mission to present a global view on the internal energy dissipation processes in GRBs, including precursors, prompt gamma-ray emission, extended soft gamma-ray emission, and late X-ray flares. The Bayesian block method is utilized to analyze the BAT lightcurves to identify various emission episodes. Our… ▽ More We jointly analyze the gamma-ray burst (GRB) data observed with BAT and XRT on board the Swift mission to present a global view on the internal energy dissipation processes in GRBs, including precursors, prompt gamma-ray emission, extended soft gamma-ray emission, and late X-ray flares. The Bayesian block method is utilized to analyze the BAT lightcurves to identify various emission episodes. Our results suggest that these emission components likely share a same physical origin, which is repeated activation of the GRB central engine. What we observe in the gamma-ray band may be the tip-of-iceberg of more extended underlying activities. The precursor emission, which is detected in about 10% of {\em Swift} GRBs, is preferably detected in those GRBs that have a massive star core-collapse origin. The soft extended emission (EE) tail, on the other hand, is preferably detected in those GRBs that have a compact star merger origin. Bright X-ray emission is detected during the BAT quiescent phases prior to subsequent gamma-ray peaks, implying that X-ray emission may be detectable prior the BAT trigger time. Future GRB alert instruments with soft X-ray capability would be essential to reveal the early stage of GRB central engine activities, sheding light into jet composition and jet launching mechanism in GRBs. △ Less

Submitted 22 May, 2014; originally announced May 2014.

Comments: 15 pages, including 15 figures and 2 tables. Accepted for publication in ApJ

GRB 121027A: long-lasting, energetic X-ray flares and clues to radiation mechanism and progenitor star

Authors: Fang-kun Peng, You-Dong Hu, Shao-Qiang Xi, Xiang-Gao Wang, Rui-Jing Lu, En-Wei Liang, Bing Zhang

Abstract: GRB 121027A is un-usual with its extremely long-lasting, energetic X-ray flares. The total energy release in X-ray flares is about one order of magnitude higher than prompt gamma-rays, making it special from most long GRBs. We show that while the prompt gamma-ray emission satisfies the empirical E_{iso}-E_{p} relation of typical long GRBs, the X-ray flares, whose spectra can be fit with a cutoff-p… ▽ More GRB 121027A is un-usual with its extremely long-lasting, energetic X-ray flares. The total energy release in X-ray flares is about one order of magnitude higher than prompt gamma-rays, making it special from most long GRBs. We show that while the prompt gamma-ray emission satisfies the empirical E_{iso}-E_{p} relation of typical long GRBs, the X-ray flares, whose spectra can be fit with a cutoff-power-law model with well-constrained E_p, significantly deviate from such a relation. Nonetheless, a time-resolved spectral analysis of X-ray flares suggest that the X-ray emission is consistent with the L_{iso}-E_{p} relation of long GRBs. We constrain the minimum Lorentz factor of the X-ray flares to be ~14, which is consistent with the Gamma-L_{iso} relation. Our results imply that prompt gamma-ray emission and late X-ray flares share the similar radiation mechanism, but originate from the outflows with different Lorentz factors. We search for similar GRBs from the Swift GRB archives, and find that the z=6.29 GRB 050904 is a carbon copy of GRB 121027A. The long-lasting, energetic X-ray flares in these GRBs demand significant accretion at late times, which point towards a large-radius progenitor star. △ Less

Submitted 20 February, 2013; originally announced February 2013.

Comments: 5 pages, 4 figures, submitted to ApJ Letters

A Comprehensive Study of Gamma-Ray Burst Optical Emission: I. Flares and Early Shallow Decay Component

Authors: Liang Li, En-Wei Liang, Qing-Wen Tang, Jie-Min Chen, Shao-Qiang Xi, Hou-Jun LV, He Gao, Bing Zhang, Jin Zhang, Shuang-Xi Yi, Rui-Jing Lu, Lian-Zhong LV, Jian-Yan Wei

Abstract: Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are compiled from the literature. By empirical fitting we identify eight possible emission components and summarize the results in a "synthetic" lightcurve. Both optical flare and early shallow-decay components are likely related to long-term central engine activities. We focus on their statistical properties in this paper. Twenty-fou… ▽ More Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are compiled from the literature. By empirical fitting we identify eight possible emission components and summarize the results in a "synthetic" lightcurve. Both optical flare and early shallow-decay components are likely related to long-term central engine activities. We focus on their statistical properties in this paper. Twenty-four optical flares are obtained from 19 GRBs. The isotropic R-band energy is smaller than 1% of $E_{γ, \rm iso}$. The relation between isotropic luminosities of the flares and gamma-rays follows $L^{\rm F}_{\rm R, iso}\propto L_{γ, \rm iso}^{1.11\pm 0.27}$. Later flares tend to be wider and dimmer, i.e., $w^{\rm F}\sim t^{\rm F}_{\rm p}/2$ and $L^{\rm F}_{\rm R, iso}\propto [t^{\rm F}_{\rm p}/(1+z)]^{-1.15\pm0.15}$. The detection probability of the optical flares is much smaller than that of X-ray flares. An optical shallow decay segment is observed in 39 GRBs. The relation between the break time and break luminosity is a power-law, with an index of $-0.78\pm 0.08$, similar to that derived from X-ray flares. The X-ray and optical breaks are usually chromatic, but a tentative correlation is found. We suggest that similar to the prompt optical emission that tracks $γ$-rays, the optical flares are also related to the erratic behavior of the central engine. The shallow decay component is likely related to a long-lasting spinning-down central engine or piling up of flare materials onto the blastwave. Mixing of different emission components may be the reason of the diverse chromatic afterglow behaviors. △ Less

Submitted 10 August, 2012; v1 submitted 11 March, 2012; originally announced March 2012.

Comments: 43 pages, 13 figures, 3 tables, accepted for publication in ApJ