BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

J. H. Fan

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

ResearcherUnclaimed source record
astro-ph.HEastro-phastro-ph.COastro-ph.GAhep-phastro-ph.IMhep-ex

Catalog footprint

What is connected

34works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

34 published item(s)

preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Constraints on heavy decaying dark matter from 570 days of LHAASO observations

The Kilometer Square Array~(KM2A) of the Large High Altitude Air Shower Observatory (LHAASO) aims at surveying the northern gamma-ray sky at energies above 10 TeV with unprecedented sensitivity. Gamma-ray observations have long been one of the most powerful tools for dark matter searches, as e.g., high-energy gamma-rays could be produced by the decays of heavy dark matter particles. In this letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340~days of data from 1/2-KM2A and 230~days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after gamma/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest gamma-ray constraints on the lifetime of heavy dark matter particles with mass between 10^5 and 10^9~GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Discovery of a new $γ$-ray source LHAASO J0341+5258 with emission up to 200TeV

We report the discovery of a new unidentified extended $γ$-ray source in the Galactic plane named LHAASO J0341+5258 with a pre-trial significance of 8.2 standard deviations above 25 TeV. The best fit position is R.A.$=55.34^{\circ}\pm0.11^{\circ}$ and Dec$=52.97^{\circ}\pm0.07^{\circ}$. The angular size of LHAASO J0341+5258 is $0.29^\circ \pm 0.06^\circ_{stat} \pm0.02^\circ_{sys}$. The flux above 25 TeV is about $20\%$ of the flux of Crab Nebula. Although a power-law fit of the spectrum from 10 TeV to 200 TeV with the photon index $α=2.98 \pm 0.19_{stat} \pm 0.02_{sys}$ is not excluded, the LHAASO data together with the flux upper limit at 10 GeV set by the Fermi LAT observation, indicate a noticeable steepening of an initially hard power-law spectrum %($α\leq 1.75$) spectrum with a cutoff at $\approx 50$ TeV. We briefly discuss the origin of UHE gamma-rays. The lack of an energetic pulsar and a young SNR inside or in the vicinity of LHAASO J0341+5258 challenge, but do not exclude both the leptonic and hadronic scenarios of gamma-ray production.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Discovery of the Ultra-high energy gamma-ray source LHAASO J2108+5157

We report the discovery of a UHE gamma-ray source, LHAASO J2108+5157, by analyzing the LHAASO-KM2A data of 308.33 live days. Significant excess of gamma-ray induced showers is observed in both energy bands of 25-100 TeV and $\gt$100 TeV with 9.5 sigma and 8.5 sigma, respectively. This source is not significantly favored as an extensive source with the angular extension smaller than the point-spread function of KM2A. The measured energy spectrum from 20 to 200 TeV can be approximately described by a power-law function with an index of -2.83$\pm$ 0.18stat. A harder spectrum is demanded at lower energies considering the flux upper limit set by Fermi-LAT observations. The position of the gamma-ray emission is correlated with a giant molecular cloud, which favors a hadronic origin. No obvious counterparts have been found, deeper multiwavelength observations will help to shed new light on this intriguing UHE source.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Evidence for particle acceleration approaching PeV energies in the W51 complex

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Exploring Lorentz Invariance Violation from Ultra-high-energy Gamma Rays Observed by LHAASO

Recently the LHAASO Collaboration published the detection of 12 ultra-high-energy gamma-ray sources above 100 TeV, with the highest energy photon reaching 1.4 PeV. The first detection of PeV gamma rays from astrophysical sources may provide a very sensitive probe of the effect of the Lorentz invariance violation (LIV), which results in decay of high-energy gamma rays in the superluminal scenario and hence a sharp cutoff of the energy spectrum. Two highest energy sources are studied in this work. No signature of the existence of LIV is found in their energy spectra, and the lower limits on the LIV energy scale are derived. Our results show that the first-order LIV energy scale should be higher than about 10^5 times the Planck scale M_{pl} and that the second-order LIV scale is >10^{-3}M_{pl}. Both limits improve by at least one order of magnitude the previous results.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Possible 3-Year Quasi-Periodic Oscillation in γ-Ray Emission from the FSRQ S5 1044+71

Variability is a typical observation feature of Fermi blazars, sometimes it shows quasi-periodic oscillation (QPO). In this work, we obtained 5-day binned light curves (with a time coverage of $\sim$ 12.9 yr) for S5 1044+71 based on Fermi LAT data, adopted five different methods: Date-compensated Discrete Fourier Transform (DCDFT), Jurkevich (JV), Lomb-Scargle Periodogram (LSP), a Fortran 90 program (REDFIT) and the Weighted Wavelet Z-transform (WWZ) to the $γ$-ray light curve, and found a possible QPO of 3.06 $\pm$ 0.43 yr at the significance level of $\sim3.6σ$. A binary black hole model including accretion model and dual-jets model is used to explain this quasi-periodic variability. We also estimated the Doppler factors and the apparent velocity for the two jet components. We speculate that this $γ$-ray quasi-periodic modulation suggest the presence of a binary supermassive black hole in S5 1044+71.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Comparison between TeV and non-TeV BL Lac Objects

BL Lacertae objects (BL Lacs) is the dominant population of TeV emitting blazars. In this work, we investigate whether there is any special observational properties for TeV sources. To do so, we will compare the observational properties of TeV detected BL Lacs (TeV BLs) and non-TeV detected BL Lac objects (non-TeV BLs). From the 3rd $Fermi$/LAT catalog (3FGL), we can get 662 BL Lacs, out of which, 47 are TeV BLs and 615 are non-TeV BLs. Their multi-wavelength flux densities ($F_{\rm R}$, $F_{\rm O}$, $F_{\rm X}$, $F_γ$), photon spectral indexes ($α_{\rm X}^{\rm ph}$, $α_γ^{\rm ph}$), and effective spectral indexes ($α_{\rm RO}$ and $α_{\rm OX}$) are compiled from the available literatures. Then the luminosities ($\log\,νL_{\rm R}$, $\log\,νL_{\rm O}$, $\log\,νL_{\rm X}$, $\log\,νL_γ$) are calculated. From comparisons, we found that TeV BLs are different from low-synchrotron-peaked BLs (LSP) and intermediate-synchrotron-peaked BLs (ISP), but TeV BLs show similar properties as high-synchrotron-peaked BLs (HSP). Therefore, we concentrated on comparison between TeV HSP BLs and non-TeV HSP BLs. Analysis results suggest that TeV HSP BLs and non-TeV HSP BLs show some differences in their $α_{\rm RO}$ and $α_γ^{\rm ph}$, while their other properties are quite similar.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiband optical variability of the blazar OJ 287 during its outbursts in 2015 -- 2016

We present recent optical photometric observations of the blazar OJ 287 taken during September 2015 -- May 2016. Our intense observations of the blazar started in November 2015 and continued until May 2016 and included detection of the large optical outburst in December 2016 that was predicted using the binary black hole model for OJ 287. For our observing campaign, we used a total of 9 ground based optical telescopes of which one is in Japan, one is in India, three are in Bulgaria, one is in Serbia, one is in Georgia, and two are in the USA. These observations were carried out in 102 nights with a total of ~ 1000 image frames in BVRI bands, though the majority were in the R band. We detected a second comparably strong flare in March 2016. In addition, we investigated multi-band flux variations, colour variations, and spectral changes in the blazar on diverse timescales as they are useful in understanding the emission mechanisms. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux, colour and spectral variability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Spectral Energy Distributions of Fermi Blazars

(Abridged) In this paper, multi-wavelength data are compiled for a sample of 1425 Fermi blazars to calculate their spectral energy distributions (SEDs). A parabolic function, $\log(νF_ν) = P_1(\logν- P_2)^2 + P_3,$ is used for SED fitting. Synchrotron peak frequency ($\log ν_p$), spectral curvature ($P_1$), peak flux ($ν_{\rm p}F_{\rm ν_p}$), and integrated flux ($νF_ν$) are successfully obtained for 1392 blazars (461 flat spectrum radio quasars-FSRQs, 620 BL Lacs-BLs and 311 blazars of uncertain type-BCUs, 999 sources have known redshifts). Monochromatic luminosity at radio 1.4 GHz, optical R band, X-ray at 1 keV and $γ$-ray at 1 GeV, peak luminosity, integrated luminosity and effective spectral indexes of radio to optical ($α_{\rm RO}$), and optical to X-ray ($α_{\rm OX}$) are calculated. The "Bayesian classification" is employed to log$ν_{\rm p}$ in the rest frame for 999 blazars with available redshift and the results show that 3 components are enough to fit the $\logν_{\rm p}$ distribution, there is no ultra high peaked subclass. Based on the 3 components, the subclasses of blazars using the acronyms of Abdo et al. (2010a) are classified, and some mutual correlations are also studied. Conclusions are finally drawn as follows: (1) SEDs are successfully obtained for 1392 blazars. The fitted peak frequencies are compared with common sources from samples available (Sambruna et al. 1996, Nieppola et al. 2006, 2008, Abdo et al. 2010a). (2) Blazars are classified as low synchrotron peak sources (LSPs) if $\logν_{\rm p}$(Hz) $\leq 14.0$, intermediate synchrotron peak sources (ISPs) if $14.0 < \logν_{\rm p}$(Hz) $\leq 15.3$, and high synchrotron peak sources (HSPs) if $\logν_{\rm p}$(Hz) $> 15.3$. (3) $γ$-ray emissions are strongly correlated with radio emissions. (...)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Beaming Effect in Fermi Blazars

The γ-ray loud blazars (flat spectrum radio quasars--FSRQs and BL Lacertae objects-BLs) are very bright in the γ-ray bands, which is perhaps associated with a beaming effect. Therefore, one can expect that the γ-ray luminosity is correlated with the beaming factor. In this paper, we investigated the relation between the radio Doppler factors and the gamma-ray luminosities. Our analysis suggests that the γ-ray luminosity be strongly correlated with the factor of δ_R for the whole sample, FSRQs, and BLs. When the effect of a common redshift is excluded, the correlation still exists for the FSRQs sub-sample suggesting that the γ-rays are strongly beamed. However, the partial correlation analysis does not show a correlation for the small BL Lac sample.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Duration distributions for different softness groups of gamma-ray bursts

Gamma-ray bursts (GRBs) are divided into two classes according to their durations. We investigate if the softness of bursts plays a role in the conventional classification of the objects. We employ the BATSE (Burst and Transient Source Experiment) catalog and analyze the duration distributions of different groups of GRBs associated with distinct softness. Our analysis reveals that the conventional classification of GRBs with the duration of bursts is influenced by the softness of the objects. There exits a bimodality in the duration distribution of GRBs for each group of bursts and the time position of the dip in the bimodality histogram shifts with the softness parameter. Our findings suggest that the conventional classification scheme should be modified by separating the two well-known populations in different softness groups, which would be more reasonable than doing so with a single sample. According to the relation between the dip position and the softness parameter, we get an empirical function that can roughly set apart the short-hard and long-soft bursts: $SP = (0.100 \pm 0.028) T_{90}^{-(0.85 \pm 0.18)}$, where $SP$ is the softness parameter adopted in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the nuclear obscuration of H2O maser galaxies

To shed light onto the circumnuclear environment of 22 GHz H2O maser galaxies, we have analyzed some of their multi-wavelength properties, including the far infrared luminosity (FIR), the luminosity of the [O III]\lambda5007 emission line, the nuclear X-ray luminosity, and the equivalent width of the neutral iron Ka emission line (EW (Ka)). Our statistical analysis includes a total of 85 sources, most of them harboring an active galactic nucleus (AGN). There are strong anti-correlations between EW (Ka) and two "optical thickness parameters", i.e. the ratios of the X-ray luminosity versus the presumably more isotropically radiated [O III] and far infrared (FIR) luminosities. Based on these anti-correlations, a set of quantitative criteria, EW (Ka) > 300eV, L_{2-10keV} < 2L_[O III] and L_{FIR} > 600L_{2-10keV} can be established for Compton-thick nuclear regions. 18 H2O maser galaxies belong to this category. There are no obvious correlations between the EW (Ka), the [O III] luminosity and the isotropic H2O maser luminosity. When comparing samples of Seyfert 2s with and without detected H2O maser lines, there seem to exist differences in EW (Ka) and the fraction of Compton-thick nuclei. This should be studied further. For AGN masers alone, there is no obvious correlation between FIR and H2O maser luminosities. However, including masers associated with star forming regions, a linear correlation is revealed. Overall, the extragalactic FIR-H2O data agree with the corresponding relation for Galactic maser sources, extrapolated by several orders of magnitude to higher luminosities.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The optical microvariability and spectral changes of the BL Lacertae object S5 0716+714

We monitored the BL Lac object S5 0716+714 in the optical band during October 2008, December 2008 and February 2009 with a best temporal resolution of about 5 minutes in the BVRI bands. Four fast flares were observed with amplitudes ranging from 0.3 to 0.75 mag. The source remained active during the whole monitoring campaign, showing microvariability in all days except for one. The overall variability amplitudes are delta B ~ 0.89 mag, delta V ~ 0.80 mag, delta R ~ 0.73 mag and delta I ~0.51 mag. Typical timescales of microvariability range from 2 to 8 hours. The overall V - R color index ranges from 0.37 to 0.59. Strong bluer- when-brighter chromatism was found on internight timescales. However, different spectral behavior was found on intranight timescales. A possible time lag of ~ 11 mins between B and I bands was found on one night. The shock-in-jet model and geometric effects can be applied to explain the source's intranight behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Optical Photometrical Observations and Variability for Quasar 4C 29.45

We reported the result of long term optical variability of the blazar 4C 29.45 (QSO 1156+295, Ton 599), carried out optical photometric observations in Johnson V, Cousins RI passbands during April 1997 to March 2002 using the 1.56 meter telescope of the Shanghai Astronomical Observatory (SHAO) at Sheshan, China, compiled the post-1974 optical photometric data of the blazar by combining our new observations with the published optical data, and found maximum variations in different passbands: $Δ$U = 4.41 mag, $Δ$B = 5.55 mag, $Δ$V = 4.53 mag, $Δ$R = 5.80 mag, and $Δ$I = 5.34 mag. The average color indices are: U$-$B = $-0.54\pm$0.18 mag, B$-$V = 0.56$\pm$0.21 mag, B$-$R = 0.93$\pm$0.18 mag, B$-$I = 1.51$\pm$0.24 mag, V$-$R = 0.44$\pm$0.15 mag and V$-$I = 1.03$\pm$0.23 mag. The post-1974 data give us an excellent opportunity to search for the existence of possible periodicity in the light curve. In search for periodicity in the R passband light curve, we performed Jurkevich test and power spectral (Fourier) analysis methods, and CLEANest algorithms to remove false signals. We found possible periods of 3.55 and 1.58 years. The possible mechanism for the periodic variability was discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The binary black hole scenario for the BL Lacertae object AO 0235+16

Recent analysis of the long term radio lightcurve of the extremely variable BL Lacertae object AO 0235+16 by Raiteri et al. (2001) have revealed the presence of recurrent outbursts in this source with a period of ~ 5.7 +/- 0.5 yr. Periodicity analysis of the optical lightcurve also show evidence for a shorter period. Here we discuss whether such a behavior can be explained by a binary black hole model where the accretion disk of one of the supermassive black holes is precessing due to the tidal effects of the companion. We estimate the mass of the accreting hole and analyze under what constraints onto the secondary mass and the orbital parameters of the system it is possible to provide a viable interpretation of the available multiwavelength data.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Optical observations of BL Lacertae from 1997 to 1999

We present the optical (V, R, and I) photometry for BL Lacertae, which was observed from 1997 through 1999, with the 1.56-m telescope at the Shanghai astronomical observatory (SHAO). After the 1997 outburst, it dimmed to a low state and then brightened again. During the period JD 2450701 to JD 2450701.5, variations of 0.40mag, 0.27mag, and 0.21mag over a time scale of 100 minutes were found for V, R, and I bands, suggesting that the variations were decreasing with wavelength. The correlation between V, R, and I is also analyzed using the DCF (Discrete Correlation Function) method. This shows that the variability in the V, R, and I bands are correlated with no time delay longer than 0.2 day.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The Radio and Gamma-Ray Luminosities of Blazars

Based on the $γ$-ray data of blazars in the third EGRET catalog and radio data at 5 GHz, we studied the correlation between the radio and $γ$-ray luminosities using two statistical methods. The first method was the partial correlation analysis method, which indicates that there exist correlations between the radio and $γ$-ray luminosities in both high and low states as well as in the average case. The second method involved a comparison of expected $γ$-ray luminosity distribution with the observed data using the Kolmogorov-- Smirnov (KS) test. In the second method, we assumed that there is a correlation between the radio and $γ$-ray luminosities and that the $γ$-ray luminosity function is proportional to the radio luminosity function. The KS test indicates that the expected gamma-ray luminosity distributions are consistent with the observed data in a reasonable parameter range. Finally, we used different $γ$-ray luminosity functions to estimate the possible 'observed' $γ$-ray luminosity distributions by GLAST.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Long-term Variability Properties and Periodicity Analysis for Blazars

In this paper, the compiled long-term optical and infrared measurements of some blazars are used to analyze the variation properties and the optical data are used to search for periodicity evidence in the lightcurve by means of the Jurkevich technique and the discrete correlation function (DCF) method. Following periods are found: 4.52-year for 3C 66A; 1.56 and 2.95 years for AO 0235+164; 14.4, 18.6 years for PKS 0735+178; 17.85 and 24.7 years for PKS 0754+100; 5.53 and 11.75 for OJ 287. 4.45, and 6.89 years for PKS 1215; 9 and 14.84 years for PKS 1219+285; 2.0, 13.5 and 22.5 for 3C273; 7.1 year for 3C279; 6.07 for PKS 1308+326; 3.0 and 16.5 years for PKS 1418+546; 2.0 and 9.35 years for PKS 1514-241; 18.18 for PKS 1807+698; 4.16 and 7.0 for 2155-304; 14 and 20 years for BL Lacertae. Some explanations have been discussed.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Infrared Variation of Blazars

In this paper, the historical infrared (JHK) data compiled from the published literature are presented here for 30 blazars. Maximum near-IR variations are found and compared with the optical ones. For the color-magnitude relation, some objects show that the color index increases with magnitude while 1253-055 shows complex behaviour, which perhaps suggests that the emission mechanism consists of, at least, two parts. The spectral indexes are in the range of $α_{IR}$ = 0.77 to 2.37.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

BL Lac Objects and Relativistic Beaming Model

In this paper: 1. The assumption of the dependence of Doppler factor on the emission frequency (Fan et al. 1993) is used to explain the observational differences between the radio-selected BL Lac objects(RBLs) and the X-Ray-Selected BL Lac objects (XBLs): a) Hubble relation; b) different multiwavelength correlations; c) different regions in the effective spectral index diagram; d) different polarization. The results suggest that RBLs and XBLs are the same. 2. From the analysis of the relation between infrared magnitude and redshift, it is proposed that the parent population of BL Lac objects should be FRI radio galaxies and FRII(G) radio galaxies showing the optical spectra of a galaxy. 3. From the superluminal motion, the assumption (Fan et al. 1993) is confirmed. 4. Based on the relation between polarization and Doppler factor (Fan et al. 1997, A&A, 327, 947), it is proposed that the ratio of the beamed luminosity to the unbeamed luminosity in the source frame of OVVs/HPQs is smaller than that of BL Lac objects.

0 citations0 reviewsNo code linkedOpen access