Logo
ScienceStack
Table of Contents
Fetching ...
Sign In

LHAASO observation of Mrk 421 during 2021 March - 2024 March: a comprehensive VHE catalog of multi-timescale outbursts and its time average behavior

The LHAASO Collaboration, Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, J. Blunier, A. V. Bukevich, C. M. Cai, Y. Y. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, E. S. Chen, G. H. Chen, H. K. Chen, L. F. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen, S. H. Chen, S. Z. Chen, T. L. Chen, X. B. Chen, X. J. Chen, X. P. Chen, Y. Chen, N. Cheng, Q. Y. Cheng, Y. D. Cheng, M. Y. Cui, S. W. Cui, X. H. Cui, Y. D. Cui, B. Z. Dai, H. L. Dai, Z. G. Dai, Danzengluobu, Y. X. Diao, A. J. Dong, X. Q. Dong, K. K. Duan, J. H. Fan, Y. Z. Fan, J. Fang, J. H. Fang, K. Fang, C. F. Feng, H. Feng, L. Feng, S. H. Feng, X. T. Feng, Y. Feng, Y. L. Feng, S. Gabici, B. Gao, Q. Gao, W. Gao, W. K. Gao, M. M. Ge, T. T. Ge, L. S. Geng, G. Giacinti, G. H. Gong, Q. B. Gou, M. H. Gu, F. L. Guo, J. Guo, K. J. Guo, X. L. Guo, Y. Q. Guo, Y. Y. Guo, R. P. Han, O. A. Hannuksela, M. Hasan, H. H. He, H. N. He, J. Y. He, X. Y. He, Y. He, S. Hernández-Cadena, B. W. Hou, C. Hou, X. Hou, H. B. Hu, S. C. Hu, C. Huang, D. H. Huang, J. J. Huang, X. L. Huang, X. T. Huang, X. Y. Huang, Y. Huang, Y. Y. Huang, A. Inventar, X. L. Ji, H. Y. Jia, K. Jia, H. B. Jiang, K. Jiang, X. W. Jiang, Z. J. Jiang, M. Jin, S. Kaci, M. M. Kang, I. Karpikov, D. Khangulyan, D. Kuleshov, K. Kurinov, Cheng Li, Cong Li, D. Li, F. Li, H. B. Li, H. C. Li, Jian Li, Jie Li, K. Li, L. Li, R. L. Li, S. D. Li, T. Y. Li, W. L. Li, X. R. Li, Y. Li, Zhe Li, Zhuo Li, E. W. Liang, Y. F. Liang, S. J. Lin, B. Liu, C. Liu, D. Liu, D. B. Liu, H. Liu, J. Liu, J. L. Liu, J. R. Liu, M. Y. Liu, R. Y. Liu, S. M. Liu, W. Liu, X. Liu, Y. Liu, Y. Liu, Y. N. Liu, Y. Q. Lou, Q. Luo, Y. Luo, H. K. Lv, B. Q. Ma, L. L. Ma, X. H. Ma, I. O. Maliy, J. R. Mao, Z. Min, W. Mitthumsiri, Y. Mizuno, G. B. Mou, A. Neronov, K. C. Y. Ng, M. Y. Ni, L. Nie, L. J. Ou, Z. W. Ou, P. Pattarakijwanich, Z. Y. Pei, D. Y. Peng, J. C. Qi, M. Y. Qi, J. J. Qin, D. Qu, A. Raza, C. Y. Ren, D. Ruffolo, A. Sáiz, D. Savchenko, D. Semikoz, L. Shao, O. Shchegolev, Y. Z. Shen, X. D. Sheng, Z. D. Shi, F. W. Shu, H. C. Song, Yu. V. Stenkin, Y. Su, D. X. Sun, H. Sun, J. X. Sun, Q. N. Sun, X. N. Sun, Z. B. Sun, N. H. Tabasam, J. Takata, P. H. T. Tam, H. B. Tan, Q. W. Tang, R. Tang, Z. B. Tang, W. W. Tian, C. N. Tong, L. H. Wan, C. Wang, D. H. Wang, G. W. Wang, H. G. Wang, J. C. Wang, K. Wang, Kai Wang, Kai Wang, L. P. Wang, L. Y. Wang, L. Y. Wang, R. Wang, W. Wang, X. G. Wang, X. J. Wang, X. Y. Wang, Y. Wang, Y. D. Wang, Z. H. Wang, Z. X. Wang, Zheng Wang, D. M. Wei, J. J. Wei, Y. J. Wei, T. Wen, S. S. Weng, C. Y. Wu, H. R. Wu, Q. W. Wu, S. Wu, X. F. Wu, Y. S. Wu, S. Q. Xi, J. Xia, J. J. Xia, G. M. Xiang, D. X. Xiao, G. Xiao, Y. F. Xiao, Y. L. Xin, H. D. Xing, Y. Xing, D. R. Xiong, B. N. Xu, C. Y. Xu, D. L. Xu, R. F. Xu, R. X. Xu, S. S. Xu, W. L. Xu, L. Xue, D. H. Yan, T. Yan, C. W. Yang, C. Y. Yang, F. F. Yang, L. L. Yang, M. J. Yang, R. Z. Yang, W. X. Yang, Z. H. Yang, Z. G. Yao, X. A. Ye, L. Q. Yin, N. Yin, X. H. You, Z. Y. You, Q. Yuan, H. Yue, H. D. Zeng, T. X. Zeng, W. Zeng, X. T. Zeng, M. Zha, B. B. Zhang, B. T. Zhang, C. Zhang, H. Zhang, H. M. Zhang, H. Y. Zhang, J. L. Zhang, J. Y. Zhang, Li Zhang, P. F. Zhang, R. Zhang, S. R. Zhang, S. S. Zhang, S. Y. Zhang, W. Zhang, W. Y. Zhang, X. Zhang, X. P. Zhang, Yi Zhang, Yong Zhang, Z. P. Zhang, J. Zhao, L. Zhao, L. Z. Zhao, S. P. Zhao, X. H. Zhao, Z. H. Zhao, F. Zheng, T. C. Zheng, B. Zhou, H. Zhou, J. N. Zhou, M. Zhou, P. Zhou, R. Zhou, X. X. Zhou, X. X. Zhou, B. Y. Zhu, C. G. Zhu, F. R. Zhu, H. Zhu, K. J. Zhu, Y. C. Zou, X. Zuo

Abstract

The Large High Altitude Air Shower Observatory (LHAASO) monitors sources within its field of view for up to 7 hours daily, achieving a duty cycle exceeding 98% and an annual point-source sensitivity of 1.5% Crab Units (CU) in the very high energy (VHE) band. This unbiased sky-survey mode facilitates systematic monitoring and investigation of outburst phenomena. In this paper, we present results from an unprecedented three-year monitoring campaign (March 2021--March 2024) of Mrk421 using LHAASO, spanning energies from 0.4 TeV to 20 TeV. We find that the blazar stayed in a quiescent state in 2021 and became active starting in 2022 with a total of 23 VHE outburst events identified, where the highest observed daily significance reaches $20\,σ$ with a flux equivalent to approximately 3.3~CU. LHAASO's continuous monitoring suggests the flaring occupancy of Mrk~421 to be around 14%. During long-term monitoring, multiwavelength (MWL) variability and correlation analyses are conducted using complementary data from Fermi-LAT, MAXI-GSC, Swift-XRT, and ZTF. A significant correlation ($>3\,σ$) is observed between X-ray and VHE bands with no detectable time lag, while the correlation between GeV and TeV bands is weaker. The flux distribution of the TeV emission during the quiescent state is different from that in the active state, implying the existence of two modes of energy dissipation in the blazar jet. Using simultaneous MWL data, we also analyzed both the long-term and outburst-period SEDs, and discussed the possible origin of the outburst events.

LHAASO observation of Mrk 421 during 2021 March - 2024 March: a comprehensive VHE catalog of multi-timescale outbursts and its time average behavior

Abstract

The Large High Altitude Air Shower Observatory (LHAASO) monitors sources within its field of view for up to 7 hours daily, achieving a duty cycle exceeding 98% and an annual point-source sensitivity of 1.5% Crab Units (CU) in the very high energy (VHE) band. This unbiased sky-survey mode facilitates systematic monitoring and investigation of outburst phenomena. In this paper, we present results from an unprecedented three-year monitoring campaign (March 2021--March 2024) of Mrk421 using LHAASO, spanning energies from 0.4 TeV to 20 TeV. We find that the blazar stayed in a quiescent state in 2021 and became active starting in 2022 with a total of 23 VHE outburst events identified, where the highest observed daily significance reaches with a flux equivalent to approximately 3.3~CU. LHAASO's continuous monitoring suggests the flaring occupancy of Mrk~421 to be around 14%. During long-term monitoring, multiwavelength (MWL) variability and correlation analyses are conducted using complementary data from Fermi-LAT, MAXI-GSC, Swift-XRT, and ZTF. A significant correlation () is observed between X-ray and VHE bands with no detectable time lag, while the correlation between GeV and TeV bands is weaker. The flux distribution of the TeV emission during the quiescent state is different from that in the active state, implying the existence of two modes of energy dissipation in the blazar jet. Using simultaneous MWL data, we also analyzed both the long-term and outburst-period SEDs, and discussed the possible origin of the outburst events.
Paper Structure (30 sections, 3 equations, 25 figures, 8 tables)

This paper contains 30 sections, 3 equations, 25 figures, 8 tables.

Table of Contents

  1. Introduction
  2. LHAASO Observation and Data Analysis
  3. TeV observations from WCDA
  4. TeV gamma-ray data from KM2A
  5. Multi-wavelength Data Analysis
  6. HE observations from Fermi-LAT
  7. X-ray observation from MAXI-GSC
  8. X-ray data from Swift-XRT
  9. Optical data from ZTF
  10. Results
  11. Light Curves
  12. Long-term behavior
  13. Flux distribution analysis
  14. Fractional variability analysis
  15. Correlation analysis
  16. ...and 15 more sections

Figures (25)

  • Figure 1: Daily Light Curves from the Mrk 421 direction in different Energy Bands from 2021 March 8 to 2024 March 8. The panels from top to bottom refer to LHAASO-WCDA data ($N_{hit}> 100$), LHAASO-KM2A data (6 - 40 TeV), Fermi-LAT data (0.1 - 100 GeV), MAXI-GSC data (2 - 20 keV), Swift-XRT data (0.3 - 10 keV ) and optical data. In the top panel, the gray horizontal line represents zero flux level, and the pink line indicates the baseline signal (see Section \ref{['sec:dutycycle']}). The red lines show the distinct $N_{s}$ states identified between change points through Bayesian blocks analysis with a 5% false positive probability and the gray shaded areas indicate the marked outburst periods.
  • Figure 2: The left panel shows the significance sky map of Mrk 421 observed by WCDA from March 2021 to March 2024 with a maximum significance of approximately $215\ \sigma$. The right panel displays the significance sky map obtained by KM2A during the same period, with a peak significance of about $14\ \sigma$.
  • Figure 3: Histograms of the daily signal events $N_{s}$ from Mrk 421. The Left and right panels present the 3 years data and the 2021 data, respectively; in each plot, the red, blue and green corresponds to the log-normal and Gaussian and combination of Gaussian and log-normal fit.
  • Figure 4: Histograms of the MWL fluxes are shown, with (a) to (f) representing the distribution and model fits for LHAASO-KM2A, Fermi-LAT, MAXI-GSC, Swift-XRT, ZTF (g-band), and ZTF (r-band), respectively. In each plot, the red line indicates the Log-normal fit, the blue line indicates the Gaussian fit and the green line represents the combined Gaussian and Log-normal fit.
  • Figure 5: Left panel: Multi-instrument fractional variability $F_\mathrm{{var}}$ as a function of the energy over the entire observation period. Right panel: The variation of $F_\mathrm{{var}}$ in LHAASO-WCDA over , with each vertical dashed line representing a time range.
  • ...and 20 more figures