-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Disentangling the radio emission of the supernova remnant W51C
Authors:
M. F. Zhang,
W. W. Tian,
D. A. Leahy,
H. Zhu,
X. H. Cui,
S. S. Shan
Abstract:
We simulate the evolution of supernova remnant (SNR) W51C. The simulation shows the existence of a new northeast edge. We present magnetic field structure of the W51 complex (SNR W51C and two HII regions W51A/B) by employing the 11 cm survey data of Effelsberg. This new edge is identified and overlaps with W51A along the line of sight, which gives a new angular diameter of about 37' for the quasi-…
▽ More
We simulate the evolution of supernova remnant (SNR) W51C. The simulation shows the existence of a new northeast edge. We present magnetic field structure of the W51 complex (SNR W51C and two HII regions W51A/B) by employing the 11 cm survey data of Effelsberg. This new edge is identified and overlaps with W51A along the line of sight, which gives a new angular diameter of about 37' for the quasi-circular remnant. In addition, we assemble the OH spectral lines (1612/1665/1720 MHz) towards the complex by employing the newly released THOR (The HI OH Recombination line survey of Milky Way) data. We find that the known 1720 MHz OH maser in the W51B/C overlap area is located away from the detected 1612/1665MHz absorption region. The latter is sitting at the peak of the HII region G49.2-0.35 within W51B.
△ Less
Submitted 10 November, 2017; v1 submitted 12 October, 2017;
originally announced October 2017.
-
The Optical Luminosity Function of Gamma-ray Bursts deduced from ROTSE-III Observations
Authors:
X. H. Cui,
X. F. Wu,
J. J. Wei,
F. Yuan,
W. K. Zheng,
E. W. Liang,
C. W. Akerlof,
M. C. B. Ashley,
H. A. Flewelling,
E. Gogus,
T. Guver,
U. Kiziloglu,
T. A. McKay,
S. B. Pandey,
E. S. Rykoff,
W. Rujopakarn,
B. E. Schaefer,
J. C. Wheeler,
S. A. Yost
Abstract:
We present the optical luminosity function (LF) of gamma-ray bursts (GRBs) estimated from a uniform sample of 58 GRBs from observations with the Robotic Optical Transient Search Experiment III (ROTSE-III). Our GRB sample is divided into two sub-samples: detected afterglows (18 GRBs), and those with upper limits (40 GRBs). The $R$ band fluxes 100s after the onset of the burst for these two sub-samp…
▽ More
We present the optical luminosity function (LF) of gamma-ray bursts (GRBs) estimated from a uniform sample of 58 GRBs from observations with the Robotic Optical Transient Search Experiment III (ROTSE-III). Our GRB sample is divided into two sub-samples: detected afterglows (18 GRBs), and those with upper limits (40 GRBs). The $R$ band fluxes 100s after the onset of the burst for these two sub-samples are derived. The optical LFs at 100s are fitted by assuming that the co-moving GRB rate traces the star-formation rate. The detection function of ROTSE-III is taken into account during the fitting of the optical LFs by using Monte Carlo simulations. We find that the cumulative distribution of optical emission at 100s is well-described with an exponential rise and power-law decay (ERPLD), broken power-law (BPL), and Schechter LFs. A single power-law (SPL) LF, on the other hand, is ruled out with high confidence.
△ Less
Submitted 15 September, 2014;
originally announced September 2014.
-
X-ray and Optical Plateau Following the Main Bursts in Gamma-Ray Bursts and SNe II-P: A hint to the similar late injection behavior?
Authors:
X. H. Cui,
R. X. Xu
Abstract:
We analyze the emission plateaus in the X-ray afterglow lightcurves of gamma-ray bursts (GRBs) and in the optical lightcurves of Type II superpernovae (SNe IIP) in order to study whether they have similar late energy injection behaviors. We show that correlations of bolometric energies (or luminosities) between the prompt explosions and the plateaus for the two phenomena are similar. The Type II S…
▽ More
We analyze the emission plateaus in the X-ray afterglow lightcurves of gamma-ray bursts (GRBs) and in the optical lightcurves of Type II superpernovae (SNe IIP) in order to study whether they have similar late energy injection behaviors. We show that correlations of bolometric energies (or luminosities) between the prompt explosions and the plateaus for the two phenomena are similar. The Type II SNe are in the low energy end of the GRBs. The bolometric energies (or luminosities) in prompt phase E_{\rm expl} (or L_{\rm expl}) and in plateau phase E_{\rm plateau} (or L_{\rm plateau}) share relations of E_{\rm expl} \propto E_{\rm plateau}^{0.73\pm 0.14} and L_{\rm expl} \propto L_{\rm plateau}^{\sim 0.70}. These results may indicate a similar late energy injection behavior to reproduce the observed plateaus of the plateaus in the two phenomena.
△ Less
Submitted 11 March, 2013; v1 submitted 1 April, 2011;
originally announced April 2011.
-
Pulsar's kicks and Gamma-ray bursts
Authors:
X. H. Cui,
H. G. Wang,
R. X. Xu,
G. J. Qiao
Abstract:
The consistence of the distributions of pulsar's kick velocities from the model of GRB and from the pulsar observations is tested based on the supernova-GRB ($γ$-ray burst) association and under the assumption that the GRB asymmetric explosions produce pulsars. The deduced distribution of kick velocity from the model of GRB and the observed kick distribution of radio pulsars are checked by K-S t…
▽ More
The consistence of the distributions of pulsar's kick velocities from the model of GRB and from the pulsar observations is tested based on the supernova-GRB ($γ$-ray burst) association and under the assumption that the GRB asymmetric explosions produce pulsars. The deduced distribution of kick velocity from the model of GRB and the observed kick distribution of radio pulsars are checked by K-S test. These two distributions are found to come from a same parent population. This result may indicate that GRBs could be really related to supernova, and that the asymmetry of GRB associated with supernova would cause the kick of pulsars.
△ Less
Submitted 12 June, 2007; v1 submitted 31 January, 2007;
originally announced January 2007.
-
On a synthetical method to constrain 3-D structure of emission regions of pulsars
Authors:
H. G. Wang,
H. Zhang,
X. H. Cui,
G. J. Qiao,
K. J. Lee,
Y. Liu,
R. X. Xu
Abstract:
The main idea and procedure are presented for a new synthetical method on constraining 3-D structure of emission regions of pulsars. With this method the emission regions can be synthetically constrained by fitting multi-wavelength observations features, e.g. pulse widths, phase offsets between radio pulse profiles and high energy light curves, and radio polarization properties. The main techniq…
▽ More
The main idea and procedure are presented for a new synthetical method on constraining 3-D structure of emission regions of pulsars. With this method the emission regions can be synthetically constrained by fitting multi-wavelength observations features, e.g. pulse widths, phase offsets between radio pulse profiles and high energy light curves, and radio polarization properties. The main technique is based on numerically calculating the emission directions along any open field line via including aberration and retardation effects under a certain magnetic field configuration, e.g. vacuum dipole field. It can be easily extended by involving some further effects or magnetic field line configuration. The role of observed linear polarization properties in the method is discussed. The future application to known X/gamma-ray and radio pulsars will be helpful to test and improve current theoretical models on pulsar emission.
△ Less
Submitted 9 January, 2007;
originally announced January 2007.
-
To Constrain the Parameter Space of the ICS Model
Authors:
X. H. Cui,
K. J. Lee,
G. J. Qiao,
Y. L. Yue,
R. X. Xu,
H. G. Wang
Abstract:
The pulse profile of pulsar gives geometric information about pulsar's radiation model. After investigating the pulse profiles of PSR B1642-03 and PSR B0950+08, we calculate the ratios of beam width and the emission height between different frequencies. We find that the ratios are almost constants as inclination angle $α$ changes from $0^{\circ}$ to $90^{\circ}$. The ratios can be used to test p…
▽ More
The pulse profile of pulsar gives geometric information about pulsar's radiation model. After investigating the pulse profiles of PSR B1642-03 and PSR B0950+08, we calculate the ratios of beam width and the emission height between different frequencies. We find that the ratios are almost constants as inclination angle $α$ changes from $0^{\circ}$ to $90^{\circ}$. The ratios can be used to test pulsar's radiation model. In particular it can well constrain the parameter space of the inverse Compton scattering (ICS) model. These constrained parameters indicate some physical implication for the ICS model.
△ Less
Submitted 1 November, 2006;
originally announced November 2006.
-
Is PSR B0943+10 a low-mass quark star?
Authors:
Y. L. Yue,
X. H. Cui,
R. X. Xu
Abstract:
A recent X-ray observation has shown that the radio pulsar PSR B0943+10, with clear drifting subpulses, has a much smaller polar cap area than that of conventional pulsars with mass of $\sim\msun$ and radius of $\sim10$ km. Zhang et al. (2005) addressed then that this new result conflicts with the standard vacuum gap model. Nonetheless, the discrepancy could be explained if PSR B0943+10 is actua…
▽ More
A recent X-ray observation has shown that the radio pulsar PSR B0943+10, with clear drifting subpulses, has a much smaller polar cap area than that of conventional pulsars with mass of $\sim\msun$ and radius of $\sim10$ km. Zhang et al. (2005) addressed then that this new result conflicts with the standard vacuum gap model. Nonetheless, the discrepancy could be explained if PSR B0943+10 is actually a low-mass quark star. It is found that the potential drop in the open-field-line region of oblique pulsars (i.e., inclination angle $α\neq 0$) might be $\sim 10^2$ times that of aligned pulsars, and that PSR B0943+10 with $α= 12.4^{\rm o}$ could be well above the deathline. We thus conclude that the Ruderman-Sutherland-type vacuum gap model still works well for this pulsar if it is a bare quark star with a mass of $\sim 0.02M_\odot$ and a radius of $\sim 2.6$ km.
△ Less
Submitted 29 August, 2006; v1 submitted 17 March, 2006;
originally announced March 2006.
-
To Understand the UV-optical Excess of RX J1856.5-3754
Authors:
Y. L. Yue,
X. H. Cui,
R. X. Xu
Abstract:
The enigma source, RX J1856.5-3754, is one of the so-called dim thermal neutron stars. Two puzzles of RXJ1856.5-3754 exist: (1) the observational X-ray spectrum is completely featureless; (2) the UV-optical intensity is about seven times larger than that given by the continuation of the blackbody model yielded by the X-ray data. Both the puzzles would not exist anymore if RX J1856.5-3754 is a lo…
▽ More
The enigma source, RX J1856.5-3754, is one of the so-called dim thermal neutron stars. Two puzzles of RXJ1856.5-3754 exist: (1) the observational X-ray spectrum is completely featureless; (2) the UV-optical intensity is about seven times larger than that given by the continuation of the blackbody model yielded by the X-ray data. Both the puzzles would not exist anymore if RX J1856.5-3754 is a low mass bare strange quark star, which is in a propeller phase with a low accretion rate. A boundary layer of RX J1856.5-3754 is suggested and modelled, from which the UV-optical emission is radiated. Free-free absorption dominates the opacity of the boundary layer, which results in the opacity to be high in UV-optical but low in X-ray bands. The star's magnetic field, spin period, as well as the accretion rate are constrained by observations.
△ Less
Submitted 2 October, 2005;
originally announced October 2005.
-
Testing the Universal Structured Jet Models of Gamma-Ray Bursts by BATSE Observations
Authors:
X. H. Cui,
E. W. Liang,
R. J. Lu
Abstract:
Assuming that the observed gamma-ray burst (GRB) rate as a function of redshift is proportional to a corrected star formation rate, we derive the empirical distribution of the viewing angles of long BATSE GRBs, $P^{\rm em}(θ)$, and the distribution of these bursts in the plane of $θ$ against redshift, $P^{\rm em}(θ, z)$, by using a tight correlation between $E_γ$) and $E_{\rm p}^{'}$). Our resul…
▽ More
Assuming that the observed gamma-ray burst (GRB) rate as a function of redshift is proportional to a corrected star formation rate, we derive the empirical distribution of the viewing angles of long BATSE GRBs, $P^{\rm em}(θ)$, and the distribution of these bursts in the plane of $θ$ against redshift, $P^{\rm em}(θ, z)$, by using a tight correlation between $E_γ$) and $E_{\rm p}^{'}$). Our results show that $P^{\rm em}(θ)$ is well fitted by a log-normal distribution centering at $\log θ/{\rm rad}=-0.76$ with a width of $σ_{\log θ}=0.57$. We test different universal structured jet models by comparing model predictions with our empirical results. To make the comparisons reasonable, an "effective" threshold, which corresponds to the sample selection criteria of the long GRB sample, is used. We find that the predictions of a two-Gaussian jet model are roughly consistent with our empirical results. A brief discussion shows that cosmological effect on the $E_γ-E_{\rm p}^{'}$ relation does not significantly affect our results, but sample selection effects on this relationship might significantly influence our results.
△ Less
Submitted 17 April, 2005;
originally announced April 2005.
-
Gamma-Ray Burst Sequences in Hardness Ratio-Peak Energy Plane
Authors:
X. H. Cui,
E. W. Liang,
R. J. Lu
Abstract:
The narrowness of the distribution of the peak energy of $νF_ν$ spectrum of gamma-ray bursts (GRBs) and the unification of GRB population are great puzzles yet to be solved. We investigate the two puzzles based on the global spectral behaviors of different GRB population in the $HR-E_{\rm{p}}$ plane (HR the spectral hardness ratio) with BATSE and HETE-2 observations. It is found that long GRBs a…
▽ More
The narrowness of the distribution of the peak energy of $νF_ν$ spectrum of gamma-ray bursts (GRBs) and the unification of GRB population are great puzzles yet to be solved. We investigate the two puzzles based on the global spectral behaviors of different GRB population in the $HR-E_{\rm{p}}$ plane (HR the spectral hardness ratio) with BATSE and HETE-2 observations. It is found that long GRBs and XRFs observed by HETE-2 seem to follow the same sequence in the $HR-E_{\rm{p}}$ plane, with the XRFs at the low end of this sequence. The long and short GRBs observed by BATSE follow significantly different sequences in the $HR-E_{\rm p}$ plane, with most of the short GRBs having a larger hardness ratio than the long GRBs at a given $E_{\rm{p}}$. These results indicate that the global spectral behaviors of the long GRB sample and the XRF sample are similar, while that of short GRBs is different. The short GRBs seem to be a unique subclass of GRBs, and they are not the higher energy extension of the long GRBs (abridged).
△ Less
Submitted 17 April, 2005;
originally announced April 2005.