Source author record

H. G. Wang

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

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

37 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

preprint2022arXiv

Circular polarization of fast radio bursts in the curvature radiation scenario

The curvature radiation is applied to the explain the circular polarization of FRBs. Significant circular polarization is reported in both apparently non-repeating and repeating FRBs. Curvature radiation can produce significant circular polarization at the wing of the radiation beam. In the curvature radiation scenario, in order to see significant circular polarization in FRBs (1) more energetic bursts, (2) burst with electrons having higher Lorentz factor, (3) a slowly rotating neutron star at the centre are required. Different rotational period of the central neutron star may explain why some FRBs have high circular polarization, while others don't. Considering possible difference in refractive index for the parallel and perpendicular component of electric field, the position angle may change rapidly over the narrow pulse window of the radiation beam. The position angle swing in FRBs may also be explained by this non-geometric origin, besides that of the rotating vector model.

preprint2022arXiv

The discovery of a rotating radio transient J1918$-$0449 with intriguing emission properties with the five hundred meter aperture spherical radio telescope

In this study, we report on a detailed single pulse analysis of the radio emission from a rotating radio transient (RRAT) J1918$-$0449 which is the first RRAT discovered with the five hundred meter aperture spherical radio telescope (FAST). The sensitive observations were carried out on 30 April 2021 using the FAST with a central frequency of 1250 MHz and a short time resolution of 49.152 $μ$s, which forms a reliable basis to probe single pulse emission properties in detail. The source was successively observed for around 2 hours. A total of 83 dispersed bursts with significance above 6$σ$ are detected over 1.8 hours. The source's DM and rotational period are determined to be 116.1$\pm$0.4 \pcm \ and 2479.21$\pm$0.03 ms, respectively. The share of registered pulses from the total number of observed period is 3.12\%. No underlying emission is detected in the averaged off pulse profile. For bursts with fluence larger than 10 Jy ms, the pulse energy follows a power-law distribution with an index of $-3.1\pm0.4$, suggesting the existence of bright pulse emission. We find that the distribution of time between subsequent pulses is consistent with a stationary Poisson process and find no evidence of clustering over the 1.8 h observations, giving a mean burst rate of one burst every 66 s. Close inspection of the detected bright pulses reveals that 21 pulses exhibit well-defined quasi-periodicities. The subpulse drifting is present in non-successive rotations with periodicity of $2.51\pm0.06$ periods. Finally, possible physical mechanisms are discussed.

preprint2021arXiv

Rotating vector model for magnetars

The modification of the rotating vector model in the case of magnetars are calculated. Magnetars may have twisted magnetic field compared with normal pulsars. The polarization position angle of magnetars will change in the case of a twisted magnetic field. For a twisted dipole field, we found that the position angle will change both vertically and horizontally. During the untwisting process of the magnetar magnetosphere, the modifications of the position angle will evolve with time monotonously. This may explain the evolution of the position angle in magnetar PSR J1622-4950 and XTE J1810-197. The relation between the emission point and the line of sight will also change. We suggest every magnetospheric models of magnetars also calculate the corresponding changes of position angle in their models. Order of magnitude estimation formula for doing this is given. This opens the possibility to extract the magnetic field geometry of magnetars from their radio polarization observations.

preprint2020arXiv

Giant micropulse emission in the Vela pulsar at C band

We present here the analysis of giant micropulses from the Vela pulsar. A total of 4187 giant micropulses with peak flux density $>$2.5 Jy were detected during almost 4 hours of observations carried out with the Yunnan 40-m radio telescope at 6800 MHz. Nine of the giant micropulses arrived approximately 3 to 4 ms earlier than the peak of average pulse profile, longer than that at lower frequencies. The remaining giant micropulses were clustered into three distributions which correspond to three main emission regions, including four occurring on the trailing edge of averaged profile.

preprint2020arXiv

Role of Fan Beam Model in Population Synthesis of Isolated Radio Pulsars

On the basis of the PSRPOPPY software package, we developed an evolution model of population synthesis for isolated radio pulsars by incorporating the fan beam model, an alternative radio emission beam model to the widely used conal beam model. With proper choice of related models and parameters, we succesfully reproduced the observational distributions of Galactic longitude ($l$) and latitude ($b$), spin period, period derivative, dispersion measure and 1.4-GHz flux density for the 1214 isolated pulsars discovered by the Parkes multibeam and Swinburne pulsar surveys. The number of underlying population of radio-loud pulsars is found to be $2.27\times10^6$, which is larger than the underlying population of radio-loud pulsars simulated with conal beam model. The model is used to estimate the number of isolated radio pulsars detected by the Galactic-plane pulsar survey with Five hundred-metre Aperture Spherical Telescope (FAST). Approximately 2700 and 240 new isolated pulsars are expected to be discovered in the inner galactic plane ($20^{\circ} < l < 90^{\circ}$, $|b| < 10^{\circ}$) and the outer galactic plane ($150^{\circ} < l < 210^{\circ}$, $|b| < 10^{\circ}$), respectively. These numbers are smaller than those estimated by the evolution models with conal beam and the snapshot models.

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.

preprint2015arXiv

A Study of Multi-frequency Polarization Pulse Profiles of Millisecond Pulsars

We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $λ^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile.

preprint2015arXiv

Operation and performance of the ICARUS-T600 cryogenic plant at Gran Sasso underground Laboratory

ICARUS T600 liquid argon time projection chamber is the first large mass electronic detector of a new generation able to combine the imaging capabilities of the old bubble chambers with the excellent calorimetric energy measurement. After the three months demonstration run on surface in Pavia during 2001, the T600 cryogenic plant was significantly revised, in terms of reliability and safety, in view of its long-term operation in an underground environment. The T600 detector was activated in Hall B of the INFN Gran Sasso Laboratory during Spring 2010, where it was operated without interruption for about three years, taking data exposed to the CERN to Gran Sasso long baseline neutrino beam and cosmic rays. In this paper the T600 cryogenic plant is described in detail together with the commissioning procedures that lead to the successful operation of the detector shortly after the end of the filling with liquid Argon. Overall plant performance and stability during the long-term underground operation are discussed. Finally, the decommissioning procedures, carried out about six months after the end of the CNGS neutrino beam operation, are reported.

preprint2014arXiv

Frequency dependence of pulse width for 150 radio normal pulsars

The frequency dependence of pulse width is studied for 150 normal pulsars, mostly selected from the European Pulsar Network, for which the 10\% multifrequency pulse widths can be well fit with the Thorsett relationship W_{10}=Aν^μ+W_{\rm 10, min}. The relative fraction of pulse width change between 0.4~GHz and 4.85~GHz, η=(W_{4.85}-W_{0.4})/W_{0.4}, is calculated in terms of the best-fit relationship for each pulsar. It is found that 81 pulsars (54\%) have $η<-10$\% (group A), showing considerable profile narrowing at high frequencies, 40 pulsars (27\%) have $-$10\%$\leqη\leq 10$\% (group B), meaning a marginal change in pulse width, and 29 pulsars (19\%) have $η>10$\% (group C), showing a remarkable profile broadening at high frequencies. The fractions of the group-A and group-C pulsars suggest that the profile narrowing phenomenon at high frequencies is more common than the profile broadening phenomenon, but a large fraction of the group-B and group-C pulsars (a total of 46\%) is also revealed. The group-C pulsars, together with a portion of group-B pulsars with a slight pulse broadening, can hardly be explained using the conventional radius-to-frequency mapping, which only applies to the profile narrowing phenomenon. Based on a recent version of the fan beam model, a type of broadband emission model, we propose that the diverse frequency dependence of pulse width is a consequence of different types of distribution of emission spectra across the emission region. The geometrical effect predicting a link between the emission beam shrinkage and the spectrum steepening is tested but disfavored.

preprint2014arXiv

ICARUS at FNAL

The INFN and the ICARUS collaboration originally developed the technology of the LAr-TPC. Located the underground LNGS Hall-B, the ICARUS T600 detector has been performed over three years with remarkable detection efficiency featuring a smooth operation, high live time, and high reliability. About 3000 CNGS neutrino events have been collected and are being actively analyzed. ICARUS will now be moved to CERN for an extensive R&D program. The T600 detector will be overhauled and complemented with a similar T150 detector. These improvements are performed in collaboration with the LBNE experiment, of which several INFN Institutions are now members. As a novelty, a SC magnetic field of about 1 T will be introduced. During 2016 it is proposed to move the experiment to FNAL where short base line neutrino beams are available, complementing the approved MicroBooNe experiment which will start operation in 2014. The ICARUS detectors at FNAL will be an important addition since, in absence of anomalies, the signals of several detectors at different distances from the target should be a copy of each other for all experimental signatures. Due to the reduced mass, in MicroBooNE the anti-neutrino signal is too weak for a sensitive comparison. Hence, a definitive clarification of the LSND anomaly requires the exploration of the anti-neutrino signal provided by the much larger T600. The magnetic field will allow separating the anti-neutrino signal from the neutrino-induced background. It is proposed to expose the T600 at the Booster NuBeam at ~700 m from target; the T150 will be located at ~150 m. The T600 will also receive >10^4 nu_e events/year from the off-axis NUMI beam peaked around 1 GeV and exploitable to prepare for the LBNE experiment. The ICARUS teams are also interested in extending the participation to other short baseline neutrino activities collaborating with existing FNAL groups.

preprint2012arXiv

Search for "anomalies" from neutrino and anti-neutrino oscillations at Delta_m^2 ~ 1eV^2 with muon spectrometers and large LAr-TPC imaging detectors

This proposal describes an experimental search for sterile neutrinos beyond the Standard Model with a new CERN-SPS neutrino beam. The experiment is based on two identical LAr-TPC's followed by magnetized spectrometers, observing the electron and muon neutrino events at 1600 and 300 m from the proton target. This project will exploit the ICARUS T600, moved from LNGS to the CERN "Far" position. An additional 1/4 of the T600 detector will be constructed and located in the "Near" position. Two spectrometers will be placed downstream of the two LAr-TPC detectors to greatly complement the physics capabilities. Spectrometers will exploit a classical dipole magnetic field with iron slabs, and a new concept air-magnet, to perform charge identification and muon momentum measurements in a wide energy range over a large transverse area. In the two positions, the radial and energy spectra of the nu_e beam are practically identical. Comparing the two detectors, in absence of oscillations, all cross sections and experimental biases cancel out, and the two experimentally observed event distributions must be identical. Any difference of the event distributions at the locations of the two detectors might be attributed to the possible existence of ν-oscillations, presumably due to additional neutrinos with a mixing angle sin^2(2theta_new) and a larger mass difference Delta_m^2_new. The superior quality of the LAr imaging TPC, in particular its unique electron-pi_zero discrimination allows full rejection of backgrounds and offers a lossless nu_e detection capability. The determination of the muon charge with the spectrometers allows the full separation of nu_mu from anti-nu_mu and therefore controlling systematics from muon mis-identification largely at high momenta.

preprint2012arXiv

Search for anomalies in the neutrino sector with muon spectrometers and large LArTPC imaging detectors at CERN

A new experiment with an intense ~2 GeV neutrino beam at CERN SPS is proposed in order to definitely clarify the possible existence of additional neutrino states, as pointed out by neutrino calibration source experiments, reactor and accelerator experiments and measure the corresponding oscillation parameters. The experiment is based on two identical LAr-TPCs complemented by magnetized spectrometers detecting electron and muon neutrino events at Far and Near positions, 1600 m and 300 m from the proton target, respectively. The ICARUS T600 detector, the largest LAr-TPC ever built with a size of about 600 ton of imaging mass, now running in the LNGS underground laboratory, will be moved at the CERN Far position. An additional 1/4 of the T600 detector (T150) will be constructed and located in the Near position. Two large area spectrometers will be placed downstream of the two LAr-TPC detectors to perform charge identification and muon momentum measurements from sub-GeV to several GeV energy range, greatly complementing the physics capabilities. This experiment will offer remarkable discovery potentialities, collecting a very large number of unbiased events both in the neutrino and antineutrino channels, largely adequate to definitely settle the origin of the observed neutrino-related anomalies.

preprint2011arXiv

Constraint on parameters of the Inverse Compton Scattering model for radio pulsars

The inverse Compton scattering (ICS) model can explain various pulse profile shapes and diversity of pulse profile evolution based on the mechanism that the radio emission is generated through inverse Compton scattering between secondary relativistic particles and radio waves from polar gap avalanches. In this paper, we study the parameter space of ICS model for 15 pulsars, which share the common pulse profile evolution phenomena that the pulse profiles are narrower at higher observing frequencies. Two key parameters, the initial Lorentz factor and the energy loss factor of secondary particles are constrained using the least square fitting method, where we fit the theoretical curve of the "beam-frequency mapping" of the ICS model to the observed pulse widths at multiple frequencies. The uncertainty of the inclination and viewing angles are taken into account in the fitting process. It is found that the initial Lorentz factor is larger than 4000, and the energy loss factor is between 20 and 560. The Lorentz factor is consistent with the prediction of the inner vacuum gap model. Such high energy loss factors suggest significant energy loss for secondary particles at altitudes of a few tens to hundreds of kilometers.

preprint2011arXiv

Long term monitoring of mode switching for PSR B0329+54

The mode switching phenomenon of PSR B0329+54 is investigated based on the long-term monitoring from September 2003 to April 2009 made with the Urumqi 25m radio telescope at 1540 MHz. At that frequency, the change of relative intensity between the leading and trailing components is the predominant feature of mode switching. The intensity ratios between the leading and trailing components are measured for the individual profiles averaged over a few minutes. It is found that the ratios follow normal distributions, where the abnormal mode has a wider typical width than the normal mode, indicating that the abnormal mode is less stable than the normal mode. Our data show that 84.9% of the time for PSR B0329+54 was in the normal mode and 15.1% was in the abnormal mode. From the two passages of eight-day quasi-continuous observations in 2004, and supplemented by the daily data observed with 15 m telescope at 610 MHz at Jodrell Bank Observatory, the intrinsic distributions of mode timescales are constrained with the Bayesian inference method. It is found that the gamma distribution with the shape parameter slightly smaller than 1 is favored over the normal, lognormal and Pareto distributions. The optimal scale parameters of the gamma distribution is 31.5 minutes for the abnormal mode and 154 minutes for the normal mode. The shape parameters have very similar values, i.e. 0.75^{+0.22}_{-0.17} for the normal mode and 0.84^{+0.28}_{-0.22} for the abnormal mode, indicating the physical mechanisms in both modes may be the same. No long-term modulation of the relative intensity ratios was found for both the modes, suggesting that the mode switching was stable. The intrinsic timescale distributions, for the first time constrained for this pulsar, provide valuable information to understand the physics of mode switching.

preprint2010arXiv

Lower bounds of altitudes for pulsar $γ$-ray radiation

Determining radiation location observationally plays a very important role in testing the pulsar radiation models. One-photon pair production in the strong magnetic field, $γ-e^{+}e^{1}$, is one of the important physical processes in pulsar radiation mechanisms. Photons near pulsar surface with sufficient energy will be absorbed in the magnetosphere and the absorption optical depth for these GeV $γ$-ray photons is usually large. In this paper, we include the aberrational, rotational and general relativistic effects and calculate the $γ$-B optical depth for $γ$-ray photons. Then we use the derived optical depth to determine the radiation altitude lower bounds for photons with given energies. As a case study, we calculate the lower bounds of radiation altitudes of Crab pulsar for photons with energy from 5 GeV to 1 TeV.

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.

preprint2007arXiv

An annular gap acceleration model for $γ$-ray emission of pulsars

If the binding energy of the pulsar's surface is not so high (the case of a neutron star), both the negative and positive charges will flow out freely from the surface of the star. The annular free flow model for $γ$-ray emission of pulsars is suggested in this paper. It is emphasized that: (1). Two kinds of acceleration regions (annular and core) need to be taken into account. The annular acceleration region is defined by the magnetic field lines that cross the null charge surface within the light cylinder. (2). If the potential drop in the annular region of a pulsar is high enough (normally the cases of young pulsars), charges in both the annular and the core regions could be accelerated and produce primary gamma-rays. Secondary pairs are generated in both regions and stream outwards to power the broadband radiations. (3). The potential drop in the annular region grows more rapidly than that in the core region. The annular acceleration process is a key point to produce wide emission beams as observed. (4). The advantages of both the polar cap and outer gap models are retained in this model. The geometric properties of the $γ$-ray emission from the annular flow is analogous to that presented in a previous work by Qiao et al., which match the observations well. (5). Since charges with different signs leave the pulsar through the annular and the core regions, respectively, the current closure problem can be partially solved.

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.