Showing 1–50 of 68 results for author: Dai, B

The Nature of the High-energy Gamma-Ray Radiation Associated with the High-redshift Blazar B3 1343+451

Authors: Fan Wu, Wen Hu, Benzhong Dai

Abstract: High-redshift blazars are the most powerful extragalactic astrophysical sources ever detected in the high-energy gamma-ray band. In this study, we present a temporal and spectral analysis of the high-redshift blazar B3 1343+451 based on 14 years of Fermi-LAT observations, spanning from 2008 August 4 to 2022 June 6 (MJD 54686-59733). We extract a seven-day binned $γ$-ray light curve in the energy r… ▽ More High-redshift blazars are the most powerful extragalactic astrophysical sources ever detected in the high-energy gamma-ray band. In this study, we present a temporal and spectral analysis of the high-redshift blazar B3 1343+451 based on 14 years of Fermi-LAT observations, spanning from 2008 August 4 to 2022 June 6 (MJD 54686-59733). We extract a seven-day binned $γ$-ray light curve in the energy range 0.1--500 GeV and identify seven outburst periods with a peak flux of $>4.32\times10^{-7} \rm ph \cdot cm^{-2} \cdot s^{-1}$. The highest seven day flux (above 100 MeV) reaches $(8.06\pm0.56)\times10^{-7} \rm erg \ cm^{-2} \ s^{-1}$ on MJD = 56,177.16, which is 10 times higher than the flux in the quiescent period. To understand the properties of distant blazar jets, we employ a standard one-zone leptonic scenario and model the multiwavelength spectral energy distributions of one quiescent and seven flaring periods. We find that the $γ$-ray spectrum is better reproduced if we assume that the dissipation region of the jet, $R_{\rm diss}$, is located within the molecular torus, where infrared emission is the dominant external photon field. We infer that the jets in higher-redshift blazars have larger power and kinetic energy, where the kinetic energy is significantly greater than the radiation power, and the jet production efficiency suggests that we need to lower the accretion efficiency. These results imply that B3 1343+451 may have a standard thin disk surrounding its massive black hole, and the jets of B3 1343+451 may not be fully explained by the Blandford--Payne process. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 13 pages, 5 figures.Accepted for publication in APJ

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

The influence of the Sun and Moon on the observation of very high energy gamma-ray sources using EAS arrays

Authors: Tao Wen, Songzhan Chen, BenZhong Dai

Abstract: With great advance of ground-based extensive air shower array, such as LHAASO and HAWC, many very high energy (VHE) gamma-ray sources have been discovered and are been monitored regardless of the day and the night. Hence, the Sun and Moon would have some compact on the observation of gamma-ray sources, which have not been taken into account in previous analysis. In this paper, the influence of the… ▽ More With great advance of ground-based extensive air shower array, such as LHAASO and HAWC, many very high energy (VHE) gamma-ray sources have been discovered and are been monitored regardless of the day and the night. Hence, the Sun and Moon would have some compact on the observation of gamma-ray sources, which have not been taken into account in previous analysis. In this paper, the influence of the Sun and Moon on the observation of very high energy gamma-ray sources when they are near the line of sight of the Sun or Moon is estimated. The tracks of all the known VHE sources are scanned and several VHE sources are found to be very close to the line of sight of the Sun or Moon during some period. The absorption of very high energy gamma-ray by sunlight is estimated with detailed method and some usefully conclusions are achieved. The main influence is the block of the Sun and Moon on gamma-ray and their shadow on the cosmic ray background. The influence is investigated considering the detector angular resolution and some strategy on data analysis are proposed to avoid the underestimation of the gamma-ray emission. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 11 pages, 7 figures, accepted to RAA

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

Discovering the Mass-Scaled Damping Timescale from Microquasars to Blazars

Authors: Haoyang Zhang, Shenbang Yang, Benzhong Dai

Abstract: Studying the variability of the accretion disks of black holes and jets is important to identify their internal physical processes. In this letter, we obtain the characteristic damping timescale of 34 blazars and seven microquasars from the Fermi-Large Area Telescope and the XMM-Newton X-ray telescope, respectively. We found that the mass-scaled characteristic timescales, ranging from the microqua… ▽ More Studying the variability of the accretion disks of black holes and jets is important to identify their internal physical processes. In this letter, we obtain the characteristic damping timescale of 34 blazars and seven microquasars from the Fermi-Large Area Telescope and the XMM-Newton X-ray telescope, respectively. We found that the mass-scaled characteristic timescales, ranging from the microquasars of stellar-mass black holes to the blazars of supermassive black holes, exhibited a linear relationship with a slope of $\sim$0.57. Given the fact the damping timescales of the $γ$-ray in the blazars are associated with the jet, we propose that the timescales of the X-ray in these microquasars are also related with the jet. The mass-scaled damping timescale that we found was consistent with the radiation of the optical accretion disk. This can be attributed to the viscous timescale at the ultraviolet-emitting radii of the disk, which can affect the jet. Our study provides a new perspective on the origin of the region of radiation and the possible disk--jet connection based on time-domain analysis. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: 11 pages, 2 figures, 3 tables, Accepted for publication in ApJ Letters

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)

A comparative study of cosmological constraints from weak lensing using Convolutional Neural Networks

Authors: Divij Sharma, Biwei Dai, Uros Seljak

Abstract: Weak Lensing (WL) surveys are reaching unprecedented depths, enabling the investigation of very small angular scales. At these scales, nonlinear gravitational effects lead to higher-order correlations making the matter distribution highly non-Gaussian. Extracting this information using traditional statistics has proven difficult, and Machine Learning based summary statistics have emerged as a powe… ▽ More Weak Lensing (WL) surveys are reaching unprecedented depths, enabling the investigation of very small angular scales. At these scales, nonlinear gravitational effects lead to higher-order correlations making the matter distribution highly non-Gaussian. Extracting this information using traditional statistics has proven difficult, and Machine Learning based summary statistics have emerged as a powerful alternative. We explore the capabilities of a discriminative, Convolutional Neural Networks (CNN) based approach, focusing on parameter constraints in the ($Ω_m$, $σ_8$) cosmological parameter space. Leveraging novel training loss functions and network representations on WL mock datasets without baryons, we show that our models achieve $\sim 5$ times stronger constraints than the power spectrum, $\sim 3$ stronger constraints than peak counts, and $\sim 2$ stronger constraints than previous CNN-learned summary statistics and scattering transforms, for noise levels relevant to Rubin or Euclid. For WL convergence maps with baryonic physics, our models achieve $\sim 2.3$ times stronger constraining power than the power spectrum at these noise levels, also outperforming previous summary statistics. To further explore the possibilities of CNNs for this task, we also discuss transfer learning where we adapt pre-trained models, trained on different tasks or datasets, for cosmological inference, finding that these do not improve the performance. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 16 pages, 2 figures. Comments Welcome

A field-level emulator for modeling baryonic effects across hydrodynamic simulations

Authors: Divij Sharma, Biwei Dai, Francisco Villaescusa-Navarro, Uros Seljak

Abstract: We develop a new and simple method to model baryonic effects at the field level relevant for weak lensing analyses. We analyze thousands of state-of-the-art hydrodynamic simulations from the CAMELS project, each with different cosmology and strength of feedback, and we find that the cross-correlation coefficient between full hydrodynamic and N-body simulations is very close to 1 down to… ▽ More We develop a new and simple method to model baryonic effects at the field level relevant for weak lensing analyses. We analyze thousands of state-of-the-art hydrodynamic simulations from the CAMELS project, each with different cosmology and strength of feedback, and we find that the cross-correlation coefficient between full hydrodynamic and N-body simulations is very close to 1 down to $k\sim10~h{\rm Mpc}^{-1}$. This suggests that modeling baryonic effects at the field level down to these scales only requires N-body simulations plus a correction to the mode's amplitude given by: $\sqrt{P_{\rm hydro}(k)/P_{\rm nbody}(k)}$. In this paper, we build an emulator for this quantity, using Gaussian processes, that is flexible enough to reproduce results from thousands of hydrodynamic simulations that have different cosmologies, astrophysics, subgrid physics, volumes, resolutions, and at different redshifts. Our emulator is accurate at the percent level and exhibits a range of validation superior to previous studies. This method and our emulator enable field-level simulation-based inference analyses and accounting for baryonic effects in weak lensing analyses. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 12 pages, 9 figures. Comments welcome

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 detection of possible $γ$-ray quasi-periodic modulation with $\sim$600 days from the blazar S2 0109+22

Authors: Haoyang Zhang, Fan Wu, Benzhong Dai

Abstract: In this work, we analyzed the long term gamma-ray data by a Fermi Large Area Telescope (Fermi-LAT) of blazar S2 0109+22, ranging from 2008 to 2023. The quasi-periodic oscillations (QPOs) of blazars aided in investigating the physical properties of internal supermassive black holes, the nature of variability, and the underlying radiation mechanism. We employed four different methods--Weighted Wavel… ▽ More In this work, we analyzed the long term gamma-ray data by a Fermi Large Area Telescope (Fermi-LAT) of blazar S2 0109+22, ranging from 2008 to 2023. The quasi-periodic oscillations (QPOs) of blazars aided in investigating the physical properties of internal supermassive black holes, the nature of variability, and the underlying radiation mechanism. We employed four different methods--Weighted Wavelet Z-transform, Lomb-Scargle periodogram, REDFIT and phase folded light curve analysis, for searching QPO signals. Our analysis identified a possible QPO behavior with a periodicity of $\sim$600 days in November 2013 to January 2023 at a significance level of 3.5 $σ$. This QPO signal sustained $\sim$9 years, corresponding to 5.6 cycles, which was in good agreement with the previously observed of periodicity $\sim$657 days in radio. We explained this phenomenon based on the accretion model and the lighthouse effect, in a binary black hole system. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: 12 pages, 8 figures,3 tables,accepted for publication in PASP

Multiscale Flow for Robust and Optimal Cosmological Analysis

Authors: Biwei Dai, Uros Seljak

Abstract: We propose Multiscale Flow, a generative Normalizing Flow that creates samples and models the field-level likelihood of two-dimensional cosmological data such as weak lensing. Multiscale Flow uses hierarchical decomposition of cosmological fields via a wavelet basis, and then models different wavelet components separately as Normalizing Flows. The log-likelihood of the original cosmological field… ▽ More We propose Multiscale Flow, a generative Normalizing Flow that creates samples and models the field-level likelihood of two-dimensional cosmological data such as weak lensing. Multiscale Flow uses hierarchical decomposition of cosmological fields via a wavelet basis, and then models different wavelet components separately as Normalizing Flows. The log-likelihood of the original cosmological field can be recovered by summing over the log-likelihood of each wavelet term. This decomposition allows us to separate the information from different scales and identify distribution shifts in the data such as unknown scale-dependent systematics. The resulting likelihood analysis can not only identify these types of systematics, but can also be made optimal, in the sense that the Multiscale Flow can learn the full likelihood at the field without any dimensionality reduction. We apply Multiscale Flow to weak lensing mock datasets for cosmological inference, and show that it significantly outperforms traditional summary statistics such as power spectrum and peak counts, as well as novel Machine Learning based summary statistics such as scattering transform and convolutional neural networks. We further show that Multiscale Flow is able to identify distribution shifts not in the training data such as baryonic effects. Finally, we demonstrate that Multiscale Flow can be used to generate realistic samples of weak lensing data. △ Less

Submitted 14 February, 2024; v1 submitted 7 June, 2023; originally announced June 2023.

Comments: 14 pages, 7 figures. Comments welcome

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)

Search for X-ray Quasi-periodicity of six AGNs using Gaussian Process method

Authors: Haoyang Zhang, Shenbang Yang, Benzhong Dai

Abstract: The quasi-periodic oscillations (QPOs) found in active galactic nuclei (AGNs) are a very interesting observational phenomenon implying an unknown physical mechanism around supermassive black holes. Several AGNs have been found to have QPO phenomena in the X-ray energy band. Long-duration X-ray observations were collected and reduced for six AGNs with a suspected QPO. The Gaussian process (GP) mode… ▽ More The quasi-periodic oscillations (QPOs) found in active galactic nuclei (AGNs) are a very interesting observational phenomenon implying an unknown physical mechanism around supermassive black holes. Several AGNs have been found to have QPO phenomena in the X-ray energy band. Long-duration X-ray observations were collected and reduced for six AGNs with a suspected QPO. The Gaussian process (GP) model celerite was used to fit the light curves and to search for the quasi-periodicity behavior. The power spectral density and parameter posterior distributions of each light curve were calculated with the optimal model. Of the six AGNs, only RE J1034+396 was found to have a QPO of about 3600 s. The other five sources do not show QPO modulation behavior. We propose that a hot spot on the accretion disk is a possible physical mechanism resulting in this quasi-periodic behavior of AGNs. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 20 pages, 12 figures, accepted for publication in ApJ

The optical spectral features of 27 Fermi blazars

Authors: Bing-Kai Zhang, Wei-Feng Tang, Chun-Xiao Wang, Qi Wu, Min Jin, Ben-Zhong Dai, Feng-Rong Zhu

Abstract: A spectral variation accompanied with flux variability is a commonly-observed phenomenon for blazars. In order to further investigate the optical spectral feature of blazars, we have collected the long-term optical V and R band data of 27 blazars (14 BL Lacs and 13 FSRQs), and calculated their optical spectral indices. The results show that the spectral indices vary with respect to the brightness… ▽ More A spectral variation accompanied with flux variability is a commonly-observed phenomenon for blazars. In order to further investigate the optical spectral feature of blazars, we have collected the long-term optical V and R band data of 27 blazars (14 BL Lacs and 13 FSRQs), and calculated their optical spectral indices. The results show that the spectral indices vary with respect to the brightness for all of these blazars. In general, the optical spectrum progressively becomes flatter (or steeper), when the brightness increases. However the spectrum changes more and more slowly, until it tends to be stable. In other words, the source becomes bluer (or redder) and then gradually stabilizes when it brightens, which are briefly named the bluer-stable-when-brighter (BSWB) and redder-stable-when-brighter (RSWB) behaviors, respectively. Thirteen of the 14 BL Lacs show the BSWB behavior, with an exception of AO 0235+164. On the contrary, most of FSRQs (10 out of 13) exhibit the RSWB trend. It is confirmed that blazars follow the two universal optical spectral behaviors, namely, BSWB and RSWB. The model of two constant-spectral-index components can well explain the optical spectral features qualitatively and quantitatively. The results illustrate that the optical emission are mainly composed of two stable-color components, i.e., less variable thermal emission and high variable synchrotron radiation. And in most cases, the thermal component of BL Lacs is redder than that of synchrotron radiation, whereas FSRQs are the opposite. △ Less

Submitted 13 January, 2023; v1 submitted 23 December, 2022; originally announced December 2022.

Comments: 9 pages, 14 figures, accepted by MNRAS

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)

Translation and Rotation Equivariant Normalizing Flow (TRENF) for Optimal Cosmological Analysis

Authors: Biwei Dai, Uros Seljak

Abstract: Our universe is homogeneous and isotropic, and its perturbations obey translation and rotation symmetry. In this work we develop Translation and Rotation Equivariant Normalizing Flow (TRENF), a generative Normalizing Flow (NF) model which explicitly incorporates these symmetries, defining the data likelihood via a sequence of Fourier space-based convolutions and pixel-wise nonlinear transforms. TR… ▽ More Our universe is homogeneous and isotropic, and its perturbations obey translation and rotation symmetry. In this work we develop Translation and Rotation Equivariant Normalizing Flow (TRENF), a generative Normalizing Flow (NF) model which explicitly incorporates these symmetries, defining the data likelihood via a sequence of Fourier space-based convolutions and pixel-wise nonlinear transforms. TRENF gives direct access to the high dimensional data likelihood p(x|y) as a function of the labels y, such as cosmological parameters. In contrast to traditional analyses based on summary statistics, the NF approach has no loss of information since it preserves the full dimensionality of the data. On Gaussian random fields, the TRENF likelihood agrees well with the analytical expression and saturates the Fisher information content in the labels y. On nonlinear cosmological overdensity fields from N-body simulations, TRENF leads to significant improvements in constraining power over the standard power spectrum summary statistic. TRENF is also a generative model of the data, and we show that TRENF samples agree well with the N-body simulations it trained on, and that the inverse mapping of the data agrees well with a Gaussian white noise both visually and on various summary statistics: when this is perfectly achieved the resulting p(x|y) likelihood analysis becomes optimal. Finally, we develop a generalization of this model that can handle effects that break the symmetry of the data, such as the survey mask, which enables likelihood analysis on data without periodic boundaries. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: 11 pages, 10 figures. Submitted to MNRAS. Comments welcome

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

SN 2015bq: A Luminous Type Ia Supernova with Early Flux Excess

Authors: Liping Li, Jujia Zhang, Benzhong Dai, Wenxiong Li, Xiaofeng Wang, Qian Zhai, Jinming Bai

Abstract: We present optical and ultraviolet (UV) observations of a luminous type Ia supernova (SN Ia) SN 2015bq characterized by the early flux excess. This SN reaches a B-band absolute magnitude at $M_B = -19.68 \pm 0.41$ mag and a peak bolometric luminosity at $L = (1.75 \pm 0.37) \times 10^{43}$ erg s$^{-1}$, with a relatively small post-maximum decline rate [$Δm_{15}(B) = 0.82 \pm 0.05$ mag]. The flux… ▽ More We present optical and ultraviolet (UV) observations of a luminous type Ia supernova (SN Ia) SN 2015bq characterized by the early flux excess. This SN reaches a B-band absolute magnitude at $M_B = -19.68 \pm 0.41$ mag and a peak bolometric luminosity at $L = (1.75 \pm 0.37) \times 10^{43}$ erg s$^{-1}$, with a relatively small post-maximum decline rate [$Δm_{15}(B) = 0.82 \pm 0.05$ mag]. The flux excess observed in the light curves of SN 2015bq a few days after the explosion, especially seen in the UV bands, might be due to the radioactive decay of $^{56}$Ni mixed into the surface. The radiation from the decay of the surface $^{56}$Ni heats the outer layer of this SN. It produces blue $U-B$ color followed by monotonically reddening in the early phase, dominated iron-group lines, and weak intermediate-mass elements absorption features in the early spectra. The scenario of enhanced $^{56}$Ni in the surface is consistent with a large amount of $^{56}$Ni ($M_{ \rm ^{56}{\rm Ni}}$ = 0.97 $\pm 0.20$ $M_{\odot}$) synthesized during the explosion. The properties of SN 2015bq are found to locate between SN 1991T and SN 1999aa, suggesting the latter two subclasses of SNe Ia may have a common origin. △ Less

Submitted 17 October, 2021; originally announced October 2021.

Comments: 18 pages, 12 figures, accepted for publication in ApJ

Kinetic powers of the relativistic jets in Mrk 421 and Mrk 501

Authors: Xiao-Chun Deng, Wen Hu, Fang-Wu Lu, Ben-Zhong Dai

Abstract: Using the standard one-zone synchrotron self-Compton (SSC) model and the Markov chain Monte Carlo (MCMC) technique, we systematically analyze the quasi-simultaneous multi-wavelength (MWL) spectral energy distributions (SEDs) of Mrk 421 and Mrk 501 during states of relatively low activity. With this model in place, a semi-analytical method is developed to examine the uncertainty in jet power estima… ▽ More Using the standard one-zone synchrotron self-Compton (SSC) model and the Markov chain Monte Carlo (MCMC) technique, we systematically analyze the quasi-simultaneous multi-wavelength (MWL) spectral energy distributions (SEDs) of Mrk 421 and Mrk 501 during states of relatively low activity. With this model in place, a semi-analytical method is developed to examine the uncertainty in jet power estimation caused by degeneracy of the radiative models. The semi-analytical method, in combination with the MCMC technique, allows us to explore the jet properties over a wide range of the variability timescale. Our results seem to support: (1) In both Mrk 421 and Mrk 501 the jets are powered by rapidly rotating black holes (BHs). The BH spin in Mrk 501 could be lower than that in Mrk 421 or possibility they are equal under the assumption of Blandford-Znajek mechanism. (2)The energy losses, which could be used to form the large-scale radio structure, are important for reconciling the differences of the kinetic power derived from the observations of the large-scale structure and the SED fitting results. Moreover, the jet energy losses in the propagation are more significant for Mrk 501 than for Mrk 421. △ Less

Submitted 18 April, 2021; originally announced April 2021.

Comments: Accepted for publication in MNRAS; 11 pages, 5 figures and 2 tables

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)

Long-term multi-wavelength variations of Fermi blazar 3C 279

Authors: Bing-Kai Zhang, Min Jin, Xiao-Yun Zhao, Li Zhang, Ben-Zhong Dai

Abstract: The long-term optical, X-ray and $γ$-ray data of blazar 3C 279 have been compiled from $Swift$-XRT, $RXTE$ PCA, $Fermi$-LAT, SMARTS and literature. The source exhibits strong variability on long time scales. Since 1980s to now, the optical $R$ band light curve spans above 32 yr, and a possible 5.6-yr-long quasi-periodic variation component has been found in it. The optical spectral behavior has be… ▽ More The long-term optical, X-ray and $γ$-ray data of blazar 3C 279 have been compiled from $Swift$-XRT, $RXTE$ PCA, $Fermi$-LAT, SMARTS and literature. The source exhibits strong variability on long time scales. Since 1980s to now, the optical $R$ band light curve spans above 32 yr, and a possible 5.6-yr-long quasi-periodic variation component has been found in it. The optical spectral behavior has been investigated. In the optical band, the mean spectral index is -1.71. The source exhibits an obvious special spectral behavior. In the low state, the source shows a clear bluer-when-brighter behavior in a sense that the optical spectrum turns harder (flatter) when the brightness increases. While in the high state, the optical spectrum is stable, that means the source spectral index does not vary with the brightness. The correlation analysis has been performed among optical, X-ray and $γ$-ray energy bands. The result indicates that the variations of $γ$-ray and X-ray bands are well correlated without time delay on the time scale of days, and their variations exhibit weak correlations with those of optical band. The variations, especial outbursts, are simultaneous, but the magnitude of variations is disproportionate. The detailed analysis reveals that the main outbursts exhibit strong correlations in different $γ$-ray, X-ray and optical bands. △ Less

Submitted 20 March, 2021; originally announced March 2021.

Comments: 17 pages, 13 figures. Accepted by Research in Astron. Astrophys.(RAA)

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)

MADLens, a python package for fast and differentiable non-Gaussian lensing simulations

Authors: Vanessa Böhm, Yu Feng, Max E. Lee, Biwei Dai

Abstract: We present MADLens a python package for producing non-Gaussian lensing convergence maps at arbitrary source redshifts with unprecedented precision. MADLens is designed to achieve high accuracy while keeping computational costs as low as possible. A MADLens simulation with only $256^3$ particles produces convergence maps whose power agree with theoretical lensing power spectra up to $L{=}10000$ wit… ▽ More We present MADLens a python package for producing non-Gaussian lensing convergence maps at arbitrary source redshifts with unprecedented precision. MADLens is designed to achieve high accuracy while keeping computational costs as low as possible. A MADLens simulation with only $256^3$ particles produces convergence maps whose power agree with theoretical lensing power spectra up to $L{=}10000$ within the accuracy limits of HaloFit. This is made possible by a combination of a highly parallelizable particle-mesh algorithm, a sub-evolution scheme in the lensing projection, and a machine-learning inspired sharpening step. Further, MADLens is fully differentiable with respect to the initial conditions of the underlying particle-mesh simulations and a number of cosmological parameters. These properties allow MADLens to be used as a forward model in Bayesian inference algorithms that require optimization or derivative-aided sampling. Another use case for MADLens is the production of large, high resolution simulation sets as they are required for training novel deep-learning-based lensing analysis tools. We make the MADLens package publicly available under a Creative Commons License (https://github.com/VMBoehm/MADLens). △ Less

Submitted 16 December, 2020; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: 10 pages, 15 figures. Update matches version submitted to journal. Acknowledgments added and typos fixed

Gaussian Process Modeling $\it{Fermi}$-LAT $γ$-ray Blazar Variability: A Sample of Blazars with $γ$-ray Quasi-periodicities

Authors: Shenbang Yang, Dahai Yan, Pengfei Zhang, Benzhong Dai, Li Zhang

Abstract: Blazar variability may be driven by stochastic processes. On the other hand, quasi-periodic oscillation (QPO) behaviors are recently reported to be detected in $\it{Fermi}$-LAT data of blazars. However, the significances of these QPO signals given by traditional Fourier-like methods are still questioned. We analyze $γ$-ray light curves of the QPO blazars with two Gaussian process methods, CARMA an… ▽ More Blazar variability may be driven by stochastic processes. On the other hand, quasi-periodic oscillation (QPO) behaviors are recently reported to be detected in $\it{Fermi}$-LAT data of blazars. However, the significances of these QPO signals given by traditional Fourier-like methods are still questioned. We analyze $γ$-ray light curves of the QPO blazars with two Gaussian process methods, CARMA and $\it{celerite}$, to examine the appropriateness of Gaussian processes for characterizing $γ$-ray light curves of blazars and the existence of the reported QPOs. We collect a sample of 27 blazars with possible $γ$-ray periodicity and generate their $\sim11$ years $\it{Fermi}$-LAT light curves. We apply the Gaussian process models to the $γ$-ray light curves, and build their intrinsic power spectral densities (PSDs). The results show that in general the $γ$-ray light curves can be characterized by CARMA and $\textit{celerite}$ models, indicating that $γ$-ray variabilities of blazars are essentially Gaussian processes. The resulting PSDs are generally the red noise shapes with slopes between $-0.6$ and $-1.7$. Possible evidence for the $γ$-ray QPOs in PKS 0537$-$441 and PG 1553$+$113 are found in the Gaussian process modelings. △ Less

Submitted 19 November, 2020; originally announced November 2020.

Comments: 21 pages, 9 figures, 3 tables. Accepted for publication in ApJ

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)

Learning effective physical laws for generating cosmological hydrodynamics with Lagrangian Deep Learning

Authors: Biwei Dai, Uros Seljak

Abstract: The goal of generative models is to learn the intricate relations between the data to create new simulated data, but current approaches fail in very high dimensions. When the true data generating process is based on physical processes these impose symmetries and constraints, and the generative model can be created by learning an effective description of the underlying physics, which enables scalin… ▽ More The goal of generative models is to learn the intricate relations between the data to create new simulated data, but current approaches fail in very high dimensions. When the true data generating process is based on physical processes these impose symmetries and constraints, and the generative model can be created by learning an effective description of the underlying physics, which enables scaling of the generative model to very high dimensions. In this work we propose Lagrangian Deep Learning (LDL) for this purpose, applying it to learn outputs of cosmological hydrodynamical simulations. The model uses layers of Lagrangian displacements of particles describing the observables to learn the effective physical laws. The displacements are modeled as the gradient of an effective potential, which explicitly satisfies the translational and rotational invariance. The total number of learned parameters is only of order 10, and they can be viewed as effective theory parameters. We combine N-body solver FastPM with LDL and apply them to a wide range of cosmological outputs, from the dark matter to the stellar maps, gas density and temperature. The computational cost of LDL is nearly four orders of magnitude lower than the full hydrodynamical simulations, yet it outperforms it at the same resolution. We achieve this with only of order 10 layers from the initial conditions to the final output, in contrast to typical cosmological simulations with thousands of time steps. This opens up the possibility of analyzing cosmological observations entirely within this framework, without the need for large dark-matter simulations. △ Less

Submitted 6 October, 2020; originally announced October 2020.

Comments: 10 pages, 6 figures

Compact Radio Emission from Nearby Galaxies with Mid-infrared Nuclear Outbursts

Authors: B. B. Dai, X. W. Shu, N. Jiang, L. M. Dou, D. Z. Liu, C. W. Yang, F. B. Zhang, T. G. Wang

Abstract: We present 5.5 GHz observations with the VLA of a sample of nearby galaxies with energetic nuclear outbursts at mid-infrared (MIR) bands. These observations reach a uniform depth down to a median rms of ~10 uJy, representing one of most sensitive searches for radio emission associated with nuclear transients. We detect radio emission in 12 out of 16 galaxies at a level of >5sigma, corresponding to… ▽ More We present 5.5 GHz observations with the VLA of a sample of nearby galaxies with energetic nuclear outbursts at mid-infrared (MIR) bands. These observations reach a uniform depth down to a median rms of ~10 uJy, representing one of most sensitive searches for radio emission associated with nuclear transients. We detect radio emission in 12 out of 16 galaxies at a level of >5sigma, corresponding to a detection rate of 75%. Such a high detection is remarkably different from previous similar searches in stellar tidal disruption events. The radio emission is compact and not resolved for the majority of sources on scales of ~<0.5" (<0.9 kpc at z<0.1). We find the possibility of the star-formation contributing to the radio emission is low, but an AGN origin remains a plausible scenario, especially for sources that show evidence of AGN activity in their optical spectra. If the detections could represent radio emission associated with nuclear transient phenomenon such as jet or outflow, we use the blast wave model by analogy with the GRB afterglows to describe the evolution of radio light curves. In this context, the observations are consistent with a decelerating jet with an energy of ~10^{51-52} erg viewed at 30\degree-60\degree off-axis at later times, suggesting that powerful jets may be ubiquitous among MIR-burst galaxies. Future continuous monitoring observations will be crucial to decipher the origin of radio emission through detections of potential flux and spectral evolution. Our results highlight the importance of radio observations to constrain the nature of nuclear MIR outbursts in galaxies. △ Less

Submitted 3 June, 2020; originally announced June 2020.

Comments: 8 pages, 4 figures, 1 table, Accepted for publication in ApJ Letters

On the Injection of Relativistic Electrons in the Jet of 3C 279

Authors: Wen Hu, Dahai Yan, Benzhong Dai, Wei Zeng, Qianglin Hu

Abstract: The acceleration of electrons in 3C 279 is investigated through analyzing the injected electron energy distribution (EED) in a time-dependent synchrotron self-Compton + external Compton emission model. In this model, it is assumed that relativistic electrons are continuously injected into the emission region, and the injected EED [$Q_e^\prime(γ^\prime)$] follows a single power-law form with low- a… ▽ More The acceleration of electrons in 3C 279 is investigated through analyzing the injected electron energy distribution (EED) in a time-dependent synchrotron self-Compton + external Compton emission model. In this model, it is assumed that relativistic electrons are continuously injected into the emission region, and the injected EED [$Q_e^\prime(γ^\prime)$] follows a single power-law form with low- and high-energy cutoffs $\rm γ_{min}'$ and $\rm γ_{max}'$, respectively, and the spectral index $n$, i.e, $Q_e^\prime(γ^\prime)\proptoγ^{\prime-n}$. This model is applied to 14 quasi-simultaneous spectral energy distributions (SEDs) of 3C 279. The Markov Chain Monte Carlo fitting technique is performed to obtain the best-fitting parameters and the uncertainties on the parameters. The results show that the injected EED is well constrained in each state. The value of $n$ is in the range of 2.5 to 3.8, which is larger than that expected by the classic non-relativistic shock acceleration. However, the large value of $n$ can be explained by the relativistic oblique shock acceleration. The flaring activity seems to be related to an increased acceleration efficiency, reflected in an increased $γ'_{\rm min}$ and electron injection power. △ Less

Submitted 24 January, 2020; originally announced January 2020.

Comments: 18 pages, 12 figures, 3 tables, accepted for publication in MNRAS

Statistical analysis on X-ray flares from the nucleus and HST-1 knot in the M87 jet

Authors: Shenbang Yang, Dahai Yan, Benzhong Dai, Pengfei Zhang, Qianqian Zhu, Jiancheng Wang, Li Zhang

Abstract: The statistical properties of X-ray flares from two separate locations (nucleus and HST-1) in the M87 jet are investigated to reveal the physical origin of the flares. We analyse the archival \textit{Chandra} data for M87, and identify 14 flares in the nucleus and 9 flares in HST-1. The peak intensity ($I_{\rm{P}}$) and the flaring duration time ($T_{\rm{fl}}$) for each flare are obtained. It is f… ▽ More The statistical properties of X-ray flares from two separate locations (nucleus and HST-1) in the M87 jet are investigated to reveal the physical origin of the flares. We analyse the archival \textit{Chandra} data for M87, and identify 14 flares in the nucleus and 9 flares in HST-1. The peak intensity ($I_{\rm{P}}$) and the flaring duration time ($T_{\rm{fl}}$) for each flare are obtained. It is found that the distributions of both $I_{\rm{P}}$ and $T_{\rm{fl}}$ for the nucleus obey a power-law form with a similar index. A similar result is also obtained for HST-1, and no significant inconsistency between the nucleus and HST-1 is found for the indices. Similar to solar X-ray flares, the power-law distributions of the flare event parameters can be well explained by a self-organized criticality (SOC) system, which are triggered by magnetic reconnection. Our results suggest that the flares from nucleus and HST-1 are possibly triggered by magnetic reconnection process. The consistent indices for the distributions of $I_{\rm{P}}$ and $T_{\rm{fl}}$ in the CORE and HST-1 indicate that the dimensions of the energy dissipation of the magnetic reconnection is identical in the two regions. A strong correlation between the flares in the two regions also suggests a similar physical origin for the flares. △ Less

Submitted 16 August, 2019; originally announced August 2019.

Comments: 10 pages, 6 figures, accepted for publication in MNRAS

High mass and halo resolution from fast low resolution simulations

Authors: Biwei Dai, Yu Feng, Uros Seljak, Sukhdeep Singh

Abstract: Generating mocks for future sky surveys requires large volumes and high resolutions, which is computationally expensive even for fast simulations. In this work we try to develop numerical schemes to calibrate various halo and matter statistics in fast low resolution simulations compared to high resolution N-body and hydrodynamic simulations. For the halos, we improve the initial condition accuracy… ▽ More Generating mocks for future sky surveys requires large volumes and high resolutions, which is computationally expensive even for fast simulations. In this work we try to develop numerical schemes to calibrate various halo and matter statistics in fast low resolution simulations compared to high resolution N-body and hydrodynamic simulations. For the halos, we improve the initial condition accuracy and develop a halo finder "relaxed-FOF", where we allow different linking length for different halo mass and velocity dispersions. We show that our relaxed-FoF halo finder improves the common statistics, such as halo bias, halo mass function, halo auto power spectrum in real space and in redshift space, cross correlation coefficient with the reference halo catalog, and halo-matter cross power spectrum. We also incorporate the potential gradient descent (PGD) method into fast simulations to improve the matter distribution at nonlinear scale. By building a lightcone output, we show that the PGD method significantly improves the weak lensing convergence tomographic power spectrum. With these improvements FastPM is comparable to the high resolution full N-body simulation of the same mass resolution, with two orders of magnitude fewer time steps. These techniques can be used to improve the halo and matter statistics of FastPM simulations for mock catalogs of future surveys such as DESI and LSST. △ Less

Submitted 14 August, 2019; originally announced August 2019.

Comments: 19 pages, 13 figures

Search for gamma-ray emission from the Sun during solar minimum with the ARGO-YBJ experiment

Authors: B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. D'Amone, Danzengluobu, I. De Mitri, B. D'Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Zhaoyang Feng, Zhenyong Feng, W. Gao, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu, M. Iacovacci , et al. (48 additional authors not shown)

Abstract: The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar mini… ▽ More The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Statistical analysis on XMM-Newton X-ray flares of Mrk 421: distributions of peak flux and flaring time duration

Authors: Dahai Yan, Shenbang Yang, Pengfei Zhang, Benzhong Dai, Jiancheng Wang, Li Zhang

Abstract: The energy dissipation mechanism in blazar jet is unknown. Blazar's flares could provide insights on this problem. Here we report statistical results of XMM-Newton X-ray flares of Mrk 421. We analyze all public XMM-Newton X-ray observations for Mrk 421, and construct the light curves. Through fitting light curves, we obtain the parameters of flare-profiles, such as peak flux ($F_{\rm p}$) and flar… ▽ More The energy dissipation mechanism in blazar jet is unknown. Blazar's flares could provide insights on this problem. Here we report statistical results of XMM-Newton X-ray flares of Mrk 421. We analyze all public XMM-Newton X-ray observations for Mrk 421, and construct the light curves. Through fitting light curves, we obtain the parameters of flare-profiles, such as peak flux ($F_{\rm p}$) and flaring time duration ($T_{\rm fl}$). It is found that both the distributions of $F_{\rm p}$ and $T_{\rm fl}$ obey a power-law form, with the same index of $α_{\rm F}=α_{\rm T}\approx1$. The statistical properties are consistent with the predictions by a self-organized criticality (SOC) system with energy dissipation in one-dimensional space. This is similar to solar flare, but with different space dimensions of the energy dissipation domain. This suggests that X-ray flaers of Mrk 421 are possibly driven by a magnetic reconnection mechanism. Moreover, in the analysis, we find that variability on timescale of $\sim1000\ $s frequently appears. Such rapid variability indicates a magnetic field of $\geq 2.1δ_{\rm D}^{-1/3}$ G ($δ_{\rm D}$ is the Doppler factor) in emission region. △ Less

Submitted 17 August, 2018; originally announced August 2018.

Comments: Accepted by ApJ

Galactic Cosmic-Ray Anisotropy in the Northern hemisphere from the ARGO-YBJ Experiment during 2008-2012

Authors: B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. D'Amone, Danzengluobu, I. De Mitri, B. D'Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Z. Y. Feng, W. Gao, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu, M. Iacovacci, R. Iuppa , et al. (47 additional authors not shown)

Abstract: This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maxi… ▽ More This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV. △ Less

Submitted 23 May, 2018; originally announced May 2018.

A gradient based method for modeling baryons and matter in halos of fast simulations

Authors: Biwei Dai, Yu Feng, Uros Seljak

Abstract: Fast N-body PM simulations with a small number of time steps such as FastPM or COLA have been remarkably successful in modeling the galaxy statistics, but their lack of small scale force resolution and long time steps cannot give accurate halo matter profiles or matter power spectrum. High resolution N-body simulations can improve on this, but lack baryonic effects, which can only be properly incl… ▽ More Fast N-body PM simulations with a small number of time steps such as FastPM or COLA have been remarkably successful in modeling the galaxy statistics, but their lack of small scale force resolution and long time steps cannot give accurate halo matter profiles or matter power spectrum. High resolution N-body simulations can improve on this, but lack baryonic effects, which can only be properly included in hydro simulations. Here we present a scheme to calibrate the fast simulations to mimic the precision of the hydrodynamic simulations or high resolution N-body simulations. The scheme is based on a gradient descent of either effective gravitational potential, which mimics the short range force, or of effective enthalpy, which mimics gas hydrodynamics and feedback. The scheme is fast and differentiable, and can be incorporated as a post-processing step into any simulation. It gives very good results for the matter power spectrum for several of the baryonic feedback and dark matter simulations, and also gives improved dark matter halo profiles. The scheme is even able to find the large subhalos, and increase the correlation coefficient between the fast simulations and the high resolution N-body or hydro simulations. It can also be used to add baryonic effects to the high resolution N-body simulations. While the method has free parameters that can be calibrated on various simulations, they can also be viewed as astrophysical nuisance parameters describing baryonic effects that can be marginalized over during the data analysis. In this view these parameters can be viewed as an efficient parametrization of baryonic effects. △ Less

Submitted 2 April, 2018; originally announced April 2018.

Comments: 20 pages, 11 figures

Around The Way: Testing $Λ$CDM with Milky Way Stellar Stream Constraints

Authors: Biwei Dai, Brant E. Robertson, Piero Madau

Abstract: Recent analyses of the Pal 5 and GD-1 tidal streams suggest that the inner dark matter halo of the Milky Way is close to spherical, in tension with predictions from collisionless N-body simulations of cosmological structure formation. We use the Eris simulation to test whether the combination of dissipative physics and hierarchical structure formation can produce Milky Way-like galaxies whose dark… ▽ More Recent analyses of the Pal 5 and GD-1 tidal streams suggest that the inner dark matter halo of the Milky Way is close to spherical, in tension with predictions from collisionless N-body simulations of cosmological structure formation. We use the Eris simulation to test whether the combination of dissipative physics and hierarchical structure formation can produce Milky Way-like galaxies whose dark matter halos match the tidal stream constraints from the GD-1 and Pal 5 clusters. We use a dynamical model of the simulated Eris galaxy to generate many realizations of the GD-1 and Pal 5 tidal streams, marginalize over observational uncertainties in the cluster galactocentric positions and velocities, and compare with the observational constraints. We find that the total density and potential of Eris contributed by baryons and dark matter satisfies constraints from the existing Milky Way stellar stream data, as the baryons both round and redistribute the dark matter during the dissipative formation of the galaxy, and provide a centrally-concentrated mass distribution that rounds the inner potential. The Eris dark matter halo or a spherical Navarro-Frenk-White dark matter work comparably well in modeling the stream data. In contrast, the equivalent dark matter-only ErisDark simulation produces a prolate halo that cannot reproduce the observed stream data. The on-going Gaia mission will provide decisive tests of the consistency between $Λ$CDM and Milky Way streams, and should distinguish between models like Eris and more spherical halos. △ Less

Submitted 2 April, 2018; originally announced April 2018.

Comments: Accepted by ApJ, 19 pages, 8 figures

Observation of the thunderstorm-related ground cosmic ray flux variations by ARGO-YBJ

Authors: B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. D Amone, Danzeng Luobu, I. De Mitri, B. D Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Zhaoyang Feng, Zhenyong Feng, W. Gao, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu, M. Iacovacci , et al. (48 additional authors not shown)

Abstract: A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be s… ▽ More A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be strongly dependent upon the intensity and polarity of the electric field measured during the course of 15 thunderstorms. In negative electric fields (i.e. accelerating negative charges downwards), the counting rates increase with increasing electric field strength. In positive fields, the rates decrease with field intensity until a certain value of the field EFmin (whose value depends on the event multiplicity), above which the rates begin increasing. By using Monte Carlo simulations, we found that this peculiar behavior can be well described by the presence of an electric field in a layer of thickness of a few hundred meters in the atmosphere above the detector, which accelerates/decelerates the secondary shower particles of opposite charge, modifying the number of particles with energy exceeding the detector threshold. These results, for the first time, give a consistent explanation for the origin of the variation of the electron/positron flux observed for decades by high altitude cosmic ray detectors during thunderstorms. △ Less

Submitted 13 February, 2018; v1 submitted 4 December, 2017; originally announced December 2017.

Comments: 17 pages, 11 figures. The paper has been accepted by PRD

Journal ref: Phys. Rev. D 97, 042001 (2018)

The role of jet properties in the spectral evolution of the powerful blazars

Authors: Wen Hu, Wei Zeng, Ben-Zhong Dai

Abstract: In the work, we explore the role of jet properties in the spectral evolution for a sample of Fermi-LAT bright blazars composed primarily by flat spectrum radio quasars (FSRQs). We introduce a near-equipartition log-parabola (NELP) model to fit the quasi-simultaneous multi-waveband spectral energy distributions (SEDs). The Markov Chain Monte Carlo (MCMC) method is employed to determine the best val… ▽ More In the work, we explore the role of jet properties in the spectral evolution for a sample of Fermi-LAT bright blazars composed primarily by flat spectrum radio quasars (FSRQs). We introduce a near-equipartition log-parabola (NELP) model to fit the quasi-simultaneous multi-waveband spectral energy distributions (SEDs). The Markov Chain Monte Carlo (MCMC) method is employed to determine the best values of spectral parameters and its errors. The correlations of synchrotron peak luminosity $L_{pk}^*$ and its peak frequency $ν_{pk}^*$, SSC-dominant factor $ξ_s$ and $ν_{pk}^*$, and their implications as to the spectral evolution were studied. The statistical analysis and the comparison with theoretical approximation indicate that: (1)the spectral evolution in $L_{pk}^*-ν_{pk}^*$ plane may not only ascribed to variation of the characteristic energy $γ_{br}^\prime$ of the employed electron energy distribution (EED), but that other quantities must be varying as well.(2) for the spectral evolution in $ξ_s-ν_{pk}^*$ plane the magnetic field $B^\prime$ may play the dominant role, and $γ_{br}^\prime$ may be considered to be the subdominant role. (3) Variation of the beaming factor $δ_b$ may be mainly responsible for the absence of a possible correlation, since $B^\prime$ is strongly anti-correlated with the size of radiation zone $R_b^\prime$. The relation is agreement with that derived from the analytical theory, and should be self-consistent with the underling acceleration mechanism of the employed EED. By assuming a conical jet with opening angle $θ_{op}\sim1/Γ_b$, we find that the $γ$-ray emission site being located at a distance $\sim10^{17}-10^{19}$ cm. △ Less

Submitted 9 February, 2018; v1 submitted 15 November, 2017; originally announced November 2017.

Comments: 15 pages, 26 figures and 3 tables

EAS age determination from the study of the lateral distribution of charged particles near the shower axis with the ARGO-YBJ experiment

Authors: ARGO-YBJ Collaboration, :, B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. D'Amone, Danzengluobu, I. De Mitri, B. D'Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Zhaoyang Feng, Zhenyong Feng, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu , et al. (48 additional authors not shown)

Abstract: The ARGO-YBJ experiment, a full coverage extensive air shower (EAS) detector located at high altitude (4300 m a.s.l.) in Tibet, China, has smoothly taken data, with very high stability, since November 2007 to the beginning of 2013. The array consisted of a carpet of about 7000 m$^2$ Resistive Plate Chambers (RPCs) operated in streamer mode and equipped with both digital and analog readout, providi… ▽ More The ARGO-YBJ experiment, a full coverage extensive air shower (EAS) detector located at high altitude (4300 m a.s.l.) in Tibet, China, has smoothly taken data, with very high stability, since November 2007 to the beginning of 2013. The array consisted of a carpet of about 7000 m$^2$ Resistive Plate Chambers (RPCs) operated in streamer mode and equipped with both digital and analog readout, providing the measurement of particle densities up to few particles per cm$^2$. The unique detector features (full coverage, readout granularity, wide dynamic range, etc) and location (very high altitude) allowed a detailed study of the lateral density profile of charged particles at ground very close to the shower axis and its description by a proper lateral distribution function (LDF). In particular, the information collected in the first 10 m from the shower axis have been shown to provide a very effective tool for the determination of the shower development stage ("age") in the energy range 50 TeV - 10 PeV. The sensitivity of the age parameter to the mass composition of primary Cosmic Rays is also discussed. △ Less

Submitted 5 July, 2017; originally announced July 2017.

Journal ref: Astroparticle Physics 93 (2017) 46-55

Search for Gamma Ray Bursts with the ARGO-YBJ Detector in Shower Mode

Authors: B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. D Amone, Danzeng Luobu, I. De Mitri, B. D Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Zhaoyang Feng, Zhenyong Feng, W. Gao, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu, M. Iacovacci , et al. (47 additional authors not shown)

Abstract: The ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013… ▽ More The ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2. △ Less

Submitted 23 March, 2017; originally announced March 2017.

Comments: 24pages and 12 figures

Journal ref: The Astrophysical Journal, 2017, 842(31)

Detection of thermal neutrons with the PRISMA-YBJ array in Extensive Air Showers selected by the ARGO-YBJ experiment

Authors: B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. D'Amone, Danzengluobu, I. De Mitri, B. D'Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Zhaoyang Feng, Zhenyong Feng, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu, M. Iacovacci, R. Iuppa , et al. (57 additional authors not shown)

Abstract: We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory i… ▽ More We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory in Yangbajing (Tibet, 4300 m a.s.l.). Due to the tight correlation between the air shower hadrons and thermal neutrons, this technique can be envisaged as a simple way to estimate the number of high energy hadrons in EAS. Coincident events generated by primary cosmic rays of energies greater than 100 TeV have been selected and analyzed. The EN-detectors have been used to record simultaneously thermal neutrons and the air shower electromagnetic component. The density distributions of both components and the total number of thermal neutrons have been measured. The correlation of these data with the measurements carried out by ARGO-YBJ confirms the excellent performance of the EN-detector. △ Less

Submitted 17 May, 2016; v1 submitted 4 December, 2015; originally announced December 2015.

Comments: 35 pages, 16 figures

4.5 years multi-wavelength observations of Mrk 421 during the ARGO-YBJ and Fermi common operation time

Authors: B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. Damone, Danzengluobu, I. De Mitri, B. D Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Zhaoyang Feng, Zhenyong Feng, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu, M. Iacovacci, R. Iuppa , et al. (46 additional authors not shown)

Abstract: We report on the extensive multi-wavelength observations of the blazar Markarian 421 (Mrk 421) covering radio to gamma-rays, during the 4.5 year period of ARGO-YBJ and Fermi common operation time, from August 2008 to February 2013. In particular, thanks to the ARGO-YBJ and Fermi data, the whole energy range from 100 MeV to 10 TeV is covered without any gap. In the observation period, Mrk 421 showe… ▽ More We report on the extensive multi-wavelength observations of the blazar Markarian 421 (Mrk 421) covering radio to gamma-rays, during the 4.5 year period of ARGO-YBJ and Fermi common operation time, from August 2008 to February 2013. In particular, thanks to the ARGO-YBJ and Fermi data, the whole energy range from 100 MeV to 10 TeV is covered without any gap. In the observation period, Mrk 421 showed both low and high activity states at all wavebands. The correlations among flux variations in different wavebands were analyzed. Seven large flares, including five X-ray flares and two GeV gamma-ray flares with variable durations (3-58 days), and one X-ray outburst phase were identified and used to investigate the variation of the spectral energy distribution with respect to a relative quiescent phase. During the outburst phase and the seven flaring episodes, the peak energy in X-rays is observed to increase from sub-keV to few keV. The TeV gamma-ray flux increases up to 0.9-7.2 times the flux of the Crab Nebula. The behavior of GeV gamma-rays is found to vary depending on the flare, a feature that leads us to classify flares into three groups according to the GeV flux variation. Finally, the one-zone synchrotron self-Compton model was adopted to describe the emission spectra. Two out of three groups can be satisfactorily described using injected electrons with a power-law spectral index around 2.2, as expected from relativistic diffuse shock acceleration, whereas the remaining group requires a harder injected spectrum. The underlying physical mechanisms responsible for different groups may be related to the acceleration process or to the environment properties. △ Less

Submitted 21 November, 2015; originally announced November 2015.

Comments: 17 pages, 9 figures, 5 tables, Accepted for publication in ApJS

Journal ref: ApJS, 222, 6 (2016)

Study of the diffuse gamma-ray emission from the Galactic plane with ARGO-YBJ

Authors: B. Bartoli, P. Bernardini, X. J. Bi, P. Branchini, A. Budano, P. Camarri, Z. Cao, R. Cardarelli, S. Catalanotti, S. Z. Chen, T. L. Chen, P. Creti, S. W. Cui, B. Z. Dai, A. D'Amone, Danzengluobu, I. De Mitri, B. D'Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Zhaoyang Feng, Zhenyong Feng, Q. B. Gou, Y. Q. Guo , et al. (56 additional authors not shown)

Abstract: The events recorded by ARGO-YBJ in more than five years of data collection have been analyzed to determine the diffuse gamma-ray emission in the Galactic plane at Galactic longitudes 25° < l < 100° and Galactic latitudes . The energy range covered by this analysis, from ~350 GeV to ~2 TeV, allows the connection of the region explored by Fermi with the multi-TeV measurements carried out by Milagro.… ▽ More The events recorded by ARGO-YBJ in more than five years of data collection have been analyzed to determine the diffuse gamma-ray emission in the Galactic plane at Galactic longitudes 25° < l < 100° and Galactic latitudes . The energy range covered by this analysis, from ~350 GeV to ~2 TeV, allows the connection of the region explored by Fermi with the multi-TeV measurements carried out by Milagro. Our analysis has been focused on two selected regions of the Galactic plane, i.e., 40° < l < 100° and 65° < l < 85° (the Cygnus region), where Milagro observed an excess with respect to the predictions of current models. Great care has been taken in order to mask the most intense gamma-ray sources, including the TeV counterpart of the Cygnus cocoon recently identified by ARGO-YBJ, and to remove residual contributions. The ARGO-YBJ results do not show any excess at sub-TeV energies corresponding to the excess found by Milagro, and are consistent with the predictions of the Fermi model for the diffuse Galactic emission. From the measured energy distribution we derive spectral indices and the differential flux at 1 TeV of the diffuse gamma-ray emission in the sky regions investigated. △ Less

Submitted 24 July, 2015; originally announced July 2015.

Comments: 11 pages, 6 figures, published in APJ

Journal ref: The Astrophysical Journal, 806:20 (11pp), 2015 June 10

The analog Resistive Plate Chamber detector of the ARGO-YBJ experiment

Authors: B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. D'Amone, Danzengluobu, I. De Mitri, B. D'Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Zhaoyang Feng, Zhenyong Feng, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu, M. Iacovacci, R. Iuppa , et al. (46 additional authors not shown)

Abstract: The ARGO-YBJ experiment has been in stable data taking from November 2007 till February 2013 at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l.). The detector consists of a single layer of Resistive Plate Chambers (RPCs) ( about 6700 m^2}) operated in streamer mode. The signal pick-up is obtained by means of strips facing one side of the gas volume. The digital readout of the signals, while a… ▽ More The ARGO-YBJ experiment has been in stable data taking from November 2007 till February 2013 at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l.). The detector consists of a single layer of Resistive Plate Chambers (RPCs) ( about 6700 m^2}) operated in streamer mode. The signal pick-up is obtained by means of strips facing one side of the gas volume. The digital readout of the signals, while allows a high space-time resolution in the shower front reconstruction, limits the measurable energy to a few hundred TeV. In order to fully investigate the 1-10 PeV region, an analog readout has been implemented by instrumenting each RPC with two large size electrodes facing the other side of the gas volume. Since December 2009 the RPC charge readout has been in operation on the entire central carpet (about 5800 m^2). In this configuration the detector is able to measure the particle density at the core position where it ranges from tens to many thousands of particles per m^2. Thus ARGO-YBJ provides a highly detailed image of the charge component at the core of air showers. In this paper we describe the analog readout of RPCs in ARGO-YBJ and discuss both the performance of the system and the physical impact on the EAS measurements. △ Less

Submitted 7 April, 2015; originally announced April 2015.

Long-term optical-infrared color variability of blazars

Authors: Bing-Kai Zhang, Xiao-Shan Zhou, Xiao-Yun Zhao, Ben-Zhong Dai

Abstract: The long-term optical and infrared color variability of blazars has been investigated with the SMARTS monitoring data. The sample in this study consists of 49 flat spectrum radio quasars and 22 BL Lac objects. The fractional variability amplitudes of each source have been calculated in both optical R band and infrared J band. Overall, the variability amplitudes of FSRQs are larger than those of BL… ▽ More The long-term optical and infrared color variability of blazars has been investigated with the SMARTS monitoring data. The sample in this study consists of 49 flat spectrum radio quasars and 22 BL Lac objects. The fractional variability amplitudes of each source have been calculated in both optical R band and infrared J band. Overall, the variability amplitudes of FSRQs are larger than those of BL Lac objects. The results also suggest that the variability amplitude of most FSRQs is larger at lower energy band (J band) than at higher one (R band), while the variability amplitude of BL Lacs are larger at higher energy band. Two types of blazars both display color variation along with the variability in brightness. However, they show different variation behaviors in general. With the whole data set, 35 FSRQs exhibit redder-when-brighter trends, and 11 FSRQs exhibit opposite trends; 11 BL Lacs follow bluer-when-brighter trends, and 7 BL Lacs follow opposite trends. The examination in detail shows that there are 10 blazars showing redder-when-brighter trend in their low state, and bluer-when-brighter or stable-when-brighter trends in their high state. Some more complicated color behaviors have also been detected in several blazars. The non-thermal jet emission and the thermal emission from the accretion disc are employed to explain the observed color behaviors. △ Less

Submitted 3 April, 2015; originally announced April 2015.

Comments: 11 pages, 5 figures. Accepted by RAA (Res. Astron. Astrophys.)

The cosmic ray proton plus helium energy spectrum measured by the ARGO-YBJ experiment in the energy range 3-300 TeV

Authors: The ARGO-YBJ Collaboration, :, B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. D'Amone, Danzengluobu, I. De Mitri, B. D'Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Zhaoyang Feng, Zhenyong Feng, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu , et al. (49 additional authors not shown)

Abstract: The ARGO-YBJ experiment is a full-coverage air shower detector located at the Yangbajing Cosmic Ray Observatory (Tibet, People's Republic of China, 4300 m a.s.l.). The high altitude, combined with the full-coverage technique, allows the detection of extensive air showers in a wide energy range and offer the possibility of measuring the cosmic ray proton plus helium spectrum down to the TeV region,… ▽ More The ARGO-YBJ experiment is a full-coverage air shower detector located at the Yangbajing Cosmic Ray Observatory (Tibet, People's Republic of China, 4300 m a.s.l.). The high altitude, combined with the full-coverage technique, allows the detection of extensive air showers in a wide energy range and offer the possibility of measuring the cosmic ray proton plus helium spectrum down to the TeV region, where direct balloon/space-borne measurements are available. The detector has been in stable data taking in its full configuration from November 2007 to February 2013. In this paper the measurement of the cosmic ray proton plus helium energy spectrum is presented in the region 3-300 TeV by analyzing the full collected data sample. The resulting spectral index is $γ= -2.64 \pm 0.01$. These results demonstrate the possibility of performing an accurate measurement of the spectrum of light elements with a ground based air shower detector. △ Less

Submitted 24 March, 2015; originally announced March 2015.

Comments: 18 pages, 8 figures, preprint submitted to Phys. Rev. D

Journal ref: Phys. Rev. D 91, 112017 (2015)

Search for GeV Gamma Ray Bursts with the ARGO-YBJ Detector: Summary of Eight Years of Observations

Authors: B. Bartoli, P. Bernardini, X. J. Bi, P. Branchini, A. Budano, P. Camarri, Z. Cao, R. Cardarelli, S. Catalanotti, S. Z. Chen, T. L. Chen, P. Creti, S. W. Cui, B. Z. Dai, A. D'Amone, Danzengluobu, I. De Mitri, B. D'Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Zhaoyang Feng, Zhenyong Feng, Q. B. Gou, Y. Q. Guo , et al. (56 additional authors not shown)

Abstract: The search for Gamma Ray Burst (GRB) emission in the energy range 1-100 GeV in coincidence with the satellite detection has been carried out using the Astrophysical Radiation with Ground-based Observatory at YangBaJing (ARGO-YBJ) experiment. The high altitude location (4300 m a.s.l.), the large active surface ($\sim$ 6700 m$^2$ of Resistive Plate Chambers), the wide field of view ($\sim 2~$sr, lim… ▽ More The search for Gamma Ray Burst (GRB) emission in the energy range 1-100 GeV in coincidence with the satellite detection has been carried out using the Astrophysical Radiation with Ground-based Observatory at YangBaJing (ARGO-YBJ) experiment. The high altitude location (4300 m a.s.l.), the large active surface ($\sim$ 6700 m$^2$ of Resistive Plate Chambers), the wide field of view ($\sim 2~$sr, limited only by the atmospheric absorption) and the high duty cycle ($>$ 86 %) make the ARGO-YBJ experiment particularly suitable to detect short and unexpected events like GRBs. With the scaler mode technique, i.e., counting all the particles hitting the detector with no measurement of the primary energy and arrival direction, the minimum threshold of $\sim$ 1 GeV can be reached, overlapping the direct measurements carried out by satellites. During the experiment lifetime, from December 17, 2004 to February 7, 2013, a total of 206 GRBs occurring within the ARGO-YBJ field of view (zenith angle $θ$ $\le$ 45$^{\circ}$) have been analyzed. This is the largest sample of GRBs investigated with a ground-based detector. Two lightcurve models have been assumed and since in both cases no significant excess has been found, the corresponding fluence upper limits in the 1-100 GeV energy region have been derived, with values as low as 10$^{-5}~$erg cm$^{-2}$. The analysis of a subset of 24 GRBs with known redshift has been used to constrain the fluence extrapolation to the GeV region together with possible cutoffs under different assumptions on the spectrum. △ Less

Submitted 19 February, 2015; originally announced February 2015.

Comments: 43 pages, 6 figures, published in APJ

Journal ref: APJ, 794:82 (12pp), 2014 October 10