BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Z. G. Dai

Z. G. Dai contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.COastro-ph.SRastro-ph.IMastro-ph.GAhep-exhep-phastro-ph.EPphysics.geo-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

33 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

Flaring-associated Complex Dynamics in Two M-dwarfs Revealed by Fast, Time-resolved Spectroscopy

Habitability of an exoplanet is believed to be profoundly affected by activities of the host stars, although the related coronal mass ejections (CMEs) are still rarely detected in solar-like and late-type stars. We here report an observational study on flares of two M-dwarfs triggered by the high-cadence survey performed by the Ground Wide-angle Camera system. In both events, the fast, time-resolved spectroscopy enables us to identify symmetric broad H$α$ emission with not only a nearly zero bulk velocity, but also a large projected maximum velocity as high as $\sim700-800\ \mathrm{km\ s^{-1}}$. This broadening could be resulted from either Stark (pressure) effect or a flaring-associated CME at stellar limb. In the context of the CME scenario, the CME mass is estimated to be $\sim4\times10^{18}$ g and $2\times10^{19}$ g. In addition, our spectral analysis reveals a temporal variation of the line center of the narrow H$α$ emission in both events. The variation amplitudes are at tens of $\mathrm{km\ s^{-1}}$, which could be ascribed to the chromospheric evaporation in one event, and to a binary scenario in the other one. With the total flaring energy determined from our photometric monitor, we show a reinforced trend in which larger the flaring energy, higher the CME mass is.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quasi-periodic Eruptions from Helium Envelope of Hydrogen-deficient Stars Stripped by Supermassive Black Holes

Quasi-periodic eruptions (QPEs), which are a new kind of X-ray bursts with a recurrence time of several hours, have been detected from supermassive black holes (SMBHs) in galactic nuclei. Recently, the two QPEs discovered by the \textit{eROSITA} show asymmetric light curves with a fast rise and a slow decline. Current models cannot explain the observational characteristics of QPEs. Here we show that QPEs can be generated from the Roche lobe overflows at each periapsis passage of an evolved star orbiting an SMBH. The properties of the companion stars are constrained via analytic estimations. We find that hydrogen-deficient post-AGB stars are promising candidates. Modules for Experiments in Stellar Astrophysics (MESA) stellar evolution code is used to construct the hydrogen-deficient stars which can fulfill the requirements, as obtained through analytical estimates, to produce the properties of QPEs, including the fast-rise and slow-decay light curves, periods, energetics, and rates. Furthermore, the extreme mass ratio $\sim 10^5$ between the SMBH and the donor will lead to a phenomenon called extreme mass-ratio inspiral (EMRI), producing millihertz gravitational waves. These QPEs would be detected as EMRI sources with electromagnetic counterparts for space-based GW detectors, such as Laser Interferometer Space Antenna (LISA) and Tianqin. They would provide a new way to measure the Hubble constant and further test the Hubble constant tension.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Standardized long gamma-ray bursts as a cosmic distance indicator

Gamma-ray bursts (GRBs) are the most luminous explosions and can be detectable out to the edge of Universe. It has long been thought they can extend the Hubble diagram to very high redshifts. Several correlations between temporal or spectral properties and GRB luminosities have been proposed to make GRBs cosmological tools. However, those correlations cannot be properly standardized. In this paper, we select a long GRB sample with X-ray plateau phases produced by electromagnetic dipole emissions from central new-born magnetars. A tight correlation is found between the plateau luminosity and the end time of the plateau in X-ray afterglows out to the redshift $z=5.91$. We standardize these long GRBs X-ray light curves to a universal behavior by this correlation for the first time, with a luminosity dispersion of 0.5 dex. The derived distance-redshift relation of GRBs is in agreement with the standard $Λ$CDM model both at low and high redshifts. The evidence of accelerating universe from this GRB sample is $3σ$, which is the highest statistical significance from GRBs to date.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Magnetar-Asteroid Impact Model for FRB 200428 Associated with an X-ray Burst from SGR 1935+2154

Very recently, an extremely bright fast radio burst (FRB) 200428 with two sub-millisecond pulses was discovered to come from the direction of the Galactic magnetar SGR 1935+2154, and an X-ray burst (XRB) counterpart was detected simultaneously. These observations favor magnetar-based interior-driven models. In this Letter, we propose a different model for FRB 200428 associated with an XRB from SGR 1935+2154, in which a magnetar with high proper velocity encounters an asteroid of mass $\sim10^{20}\,$g. This infalling asteroid in the stellar gravitational field is first possibly disrupted tidally into a great number of fragments at radius $\sim {\rm a\,\,few}$ times $10^{10}\,$cm, and then slowed around the Alfv$\acute{\rm e}$n radius by an ultra-strong magnetic field and in the meantime two major fragments of mass $\sim 10^{17}\,$g that cross magnetic field lines produce two pulses of FRB 200428. The whole asteroid is eventually accreted onto the poles along magnetic field lines, impacting the stellar surface, creating a photon-e$^\pm$ pair fireball trapped initially in the stellar magnetosphere, and further leading to an XRB. We show that this gravitationally-powered model can interpret all of the observed features self-consistently.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fast radio bursts from activities of neutron stars newborn in BNS mergers: offset, birth rate and observational properties

Young neutron stars (NSs) born in core-collapse explosions are promising candidates for the central engines of fast radio bursts (FRBs), since the first localized repeating burst FRB 121102 happens in a star forming dwarf galaxy, which is similar to the host galaxies of superluminous supernovae (SLSNe) and long gamma-ray bursts (LGRBs). However, FRB 180924 and FRB 190523 are localized to massive galaxies with low rates of star formation, compared with the host of FRB 121102. Meanwhile, the offsets between the bursts and host centers are about 4 kpc and 29 kpc for FRB 180924 and FRB 190523, respectively. These properties of hosts are similar to short gamma-ray bursts \textbf{(SGRBs)}, which are produced by mergers of binary neutron star (BNS) or neutron star-black hole (NS-BH). Therefore, the NSs powering FRBs may be formed in BNS mergers. In this paper, we study the BNS merger rates, merger times, and predict their most likely merger locations for different types of host galaxies using population synthesis method. We find that the BNS merger channel is consistent with the recently reported offsets of FRB 180924 and FRB 190523. The offset distribution of short GRBs is well reproduced by population synthesis using galaxy model which is similar to GRB hosts. The event rate of FRBs (including non-repeating and repeating), is larger than those of BNS merger and short GRBs, which requires a large fraction of observed FRBs emitting several bursts. Using curvature radiation by bunches in NS magnetospheres, we also predict the observational properties of FRBs from BNS mergers, including the dispersion measure, and rotation measure. At late times ($t\geq1$yr), the contribution to dispersion measure and rotation measure from BNS merger ejecta could be neglected.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Galactic and cosmological fast radio bursts as scaled-up solar radio bursts

Fast radio bursts (FRBs) are bright milliseconds radio transients with large dispersion measures. Recently, FRB 200428 was detected in temporal coincidence with a hard X-ray flare from the Galactic magnetar SGR 1935+2154, which supports that at least some FRBs are from magnetar activity. Interestingly, a portion of X-ray flares from magnetar XTE J1810-197 and the Sun are also accompanied by radio bursts. Many features of Galactic FRB 200428 and cosmological FRBs resemble solar radio bursts. However, a common physical origin among FRBs, magnetar radio pulses and solar radio bursts has not yet been established. Here we report a universal correlation between X-ray luminosity and radio luminosity over twenty orders of magnitude among solar type III radio bursts, XTE J1810-197 and Galactic FRB 200428. This universal correlation reveals that the energetic electrons which produce the X-ray flares and those which cause radio emissions have a common origin, which can give stringent limits on the generation process of radio bursts. Moreover, we find similar occurrence frequency distributions of energy, duration and waiting time for solar radio bursts, SGR 1935+2154 and repeating FRB 121102, which also support the tight correlation and the X-ray flares temporally associated with radio bursts. All of these distributions can be understood by avalanche models of self-organized criticality systems. The universal correlation and statistical similarities indicate that the Galactic FRB 200428 and FRBs seen at cosmological distances can be treated as scaled-up solar radio bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the Distance of SGR 1935+2154 Associated with FRB 200428 and Hosted in SNR G57.2+0.8

Owing to the detection of an extremely bright fast radio burst (FRB) 200428 associated with a hard X-ray counterpart from the magnetar soft gamma-ray repeater (SGR) 1935+2154, the distance of SGR 1935+2154 potentially hosted in the supernova remnant (SNR) G57.2+0.8 can be revisited. Under the assumption that the SGR and the SNR are physically related, in this Letter, by investigating the dispersion measure (DM) of the FRB contributed by the foreground medium of our Galaxy and the local environments and combining with other observational constraints, we find that the distance of SGR 1935+2154 turns out to be $9.0\pm2.5\,$kpc and the SNR radius falls into $10$ to $18\,$pc since the local DM contribution is as low as $0-18\,$pc cm$^{-3}$. These results are basically consistent with the previous studies. In addition, an estimate for the Faraday rotation measure of the SGR and SNR is also carried out.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Periodic Fast Radio Bursts as a Probe of Extragalactic Asteroid Belts

The periodic activity of repeating fast radio burst (FRB) 180916.J0158+65 was recently reported by the CHIME/FRB Collaboration team. 28 bursts from this source not only show a $\sim16$-day period with an active phase of $\sim 4.0$ days but also exhibit a broken power law in differential energy distribution. In this paper, we suggest that FRB 180916.J0158+65-like periodic FRBs would provide a unique probe of extragalactic asteroid belts (EABs), based on our previously-proposed pulsar-EAB impact model, in which repeating FRBs arise from an old-aged, slowly-spinning, moderately-magnetized pulsar traveling through an EAB around another stellar-mass object. These two objects form a binary and thus the observed period is in fact the orbital period. We show that this model can be used to well interpret all of the observed data of FRB 180916.J0158+65. Furthermore, we constrain the EAB&#39;s physical properties and find that (1) the outer radius of the EAB is at least an order of magnitude smaller than that of its analogue in the solar system, (2) the differential size distribution of the EAB&#39;s asteroids at small diameters (large diameters) is shallower (steeper) than that of solar-system small objects, and (3) the two belts have a comparable mass.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reply to &#34;Rapid $^{14}$C excursion at 3372-3371 BCE not observed at two different locations&#34;

The nuclide $^{14}$C can be produced in the atmosphere by high energy particles and $γ$-rays from high-energy phenomena. Through the carbon cycle, some of $^{14}$CO$_2$ produced in the atmosphere can be retained in annual tree rings. Four events of rapid increase of the $^{14}$C content occurred in AD 775, AD 994, BC 660 and BC 3371 were found. Recently, the data of Jull et al. (2020) was inconsistent with our records around BC 3371. We measured our sample again and found the $^{14}$C records are consistent with the value in Wang et al. (2017). Therefore, our $^{14}$C records are robust. The inconsistency may be caused by the difference of calendar ages for the wood samples, or the physical origin of the event. First, crossdating on ring width can be performed only between trees whose growth has the same environmental conditions. Because the master tree-ring for dendrochronology is lack for Chinese trees. The master tree-ring from California has to be used. Therefore, the calendar ages derived from dendrochronology may be not precise. Second, the $^{14}$C even may be not global. One evidence is the variation of $^{14}$C content around AD 1006. The $^{14}$C contents of Californian trees increase 12\textperthousand~ in two years, while Japanese trees show no $^{14}$C increase.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The braking index of PSR B0540-69 and the associated pulsar wind nebula emission after spin-down rate transition

In Dec. 2011 PSR B0540-69 experienced a spin-down rate transition (SRT), after which the spin-down power of the pulsar increased by ~36%. About 1000 days after the SRT, the X-ray luminosity of the associated pulsar wind nebula (PWN) was found to brighten by 32+/-8%. After the SRT, the braking index n of PSR B0540-69 changes from n=2.12 to n=0.03 and then keeps this value for about five years before rising to n=0.9 in the following years. We find that most of the current models have difficulties in explaining the measured braking index. One exceptive model of the braking index evolution is the increasing dipole magnetic field of PSR B0540-69. We suggest that the field increase may result from some instabilities within the pulsar core that enhance the poloidal component at the price of toroidal component of the magnetic field. The increasing dipole magnetic field will result in the X-ray brightening of the PWN. We fit the PWN X-ray light curve by two models: one assumes a constant magnetic field within the PWN during the brightening and the other assumes an enhanced magnetic field proportional to the energy density of the PWN. It appears that the two models fit the data well, though the later model seems to fit the data a bit better. This provides marginal observational evidence that magnetic field in the PWN is generated by the termination shock. Future high-quality and high-cadence data are required to draw a solid conclusion.

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
preprint2020arXiv

Understanding FRB 200428 in the synchrotron maser shock model: consistency and possible challenge

Recently, the discovery of Galactic FRB 200428 associated with a X-ray burst (XRB) of SGR 1935+2154 has built a bridge between FRBs and magnetar activities. In this paper, we assume that the XRB occurs in the magnetar magnetosphere. We show that the observational properties of FRB 200428 and the associated XRB are consistent with the predictions of synchrotron maser emission at ultrarelativistic magnetized shocks, including radiation efficiency, similar energy occurrence frequency distributions, and simultaneous arrive times. It requires that the upstream medium is a mildly relativistic baryonic shell ejected by a previous flare. The energy injection by flares responsible for the radio bursts will produce a magnetar wind nebula, which has been used to explain the persistent radio source associated FRB 121102. We find that the radio continuum around SGR 1935+2154 can be well understood in the magnetar wind nebula model, by assuming the same energy injection rate $\dot{E} \propto t^{-1.37}$ as FRB 121102. The required baryonic mass is also estimated form the observations of FRB 121102 by GBT and FAST. By assuming the same radiation efficiency $η\sim 10^{-5}$, the total baryonic mass ejected from the central magnetar is about 0.005 solar mass. This value is much larger than the typical mass of a magnetar outer crust, but is comparable to the total mass of a magnetar crust.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A bright electromagnetic counterpart to extreme mass ratio inspirals

The extreme mass ratio inspiral (EMRI), defined as a stellar-mass compact object inspiraling into a supermassive black hole (SMBH), has been widely argued to be a low-frequency gravitational wave (GW) source. EMRIs providing accurate measurements of black hole mass and spin, are one of the primary interests for Laser Interferometer Space Antenna (LISA). However, it is usually believed that there are no electromagnetic (EM) counterparts to EMRIs. Here we show a new formation channel of EMRIs with tidal disruption flares as EM counterparts. In this scenario, flares can be produced from the tidal stripping of the helium (He) envelope of a massive star by an SMBH. The left compact core of the massive star will evolve into an EMRI. We find that, under certain initial eccentricity and semimajor axis, the GW frequency of the inspiral can enter LISA band within 10 $\sim$ 20 years, which makes the tidal disruption flare an EM precursor to EMRI. Although the event rate is just $2\times 10^{-4}~\rm Gpc^{-3}yr^{-1}$, this association can not only improve the localization accuracy of LISA and help to find the host galaxy of EMRI, but also serve as a new GW standard siren for cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Gamma-ray Emission of 60Fe and 26Al Radioactivities in our Galaxy

The isotopes $^{60}$Fe and $^{26}$Al originate from massive stars and their supernovae, reflecting ongoing nucleosynthesis in the Galaxy. We studied the gamma-ray emission from these isotopes at characteristic energies 1173, 1332, and 1809 keV with over 15 years of SPI data, finding a line flux in $^{60}$Fe combined lines of $(0.31\pm 0.06) \times 10^{-3}$ ph cm$^{-2}$ s$^{-1}$ and the $^{26}$Al line flux of $(16.8\pm 0.7) \times 10^{-4}$ ph cm$^{-2}$ s$^{-1}$ above the background and continuum emission for the whole sky. Based on the exponential-disk grid maps, we characterise the emission extent of $^{26}$Al to find scale parameters $R_0 =7.0^{+1.5}_{-1.0}$ kpc and $z_0=0.8^{+0.3}_{-0.2}$ kpc, however the $^{60}$Fe lines are too weak to spatially constrain the emission. Based on a point source model test across the Galactic plane, the $^{60}$Fe emission would not be consistent with a single strong point source in the Galactic center or somewhere else, providing a hint for a diffuse nature. We carried out comparisons of emission morphology maps using different candidate-source tracers for both $^{26}$Al and $^{60}$Fe emissions, and suggests that the $^{60}$Fe emission is more likely to be concentrated towards the Galactic plane. We determine the $^{60}$Fe/$^{26}$Al $γ$-ray flux ratio at $(18.4\pm4.2)\,\%$ , when using a parameterized spatial morphology model. Across the range of plausible morphologies, it appears possible that $^{26}$Al and $^{60}$Fe are distributed differently in the Galaxy. Using the best fitting maps for each of the elements, we constrain flux ratios in the range 0.2--0.4. We discuss its implications for massive star models and their nucleosynthesis.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

eXTP -- enhanced X-ray Timing and Polarimetry Mission

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

0 citations0 reviewsNo code linkedOpen access