Researcher profile

J. Chang

J. Chang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
37works
0followers
10topics
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

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

Magnetic Excitations of a Half-Filled Tl-based Cuprate

Strong electron correlations drive Mott insulator transitions. Yet, there exists no framework to classify Mott insulators by their degree of correlation. Cuprate superconductors, with their tunable doping and rich phase diagrams, offer a unique platform to investigate the evolution of these interactions. However, spectroscopic access to a clean half-filled Mott-insulating state is lacking in compounds with the highest superconducting onset temperature. To fill this gap, we introduce a pristine, half-filled thallium-based cuprate system, Tl$_2$Ba$_5$Cu$_4$O$_{x}$. Using high-resolution resonant inelastic x-ray scattering, we probe long-lived magnon excitations and uncover a pronounced kink in the magnon dispersion, marked by a simultaneous change in group velocity and lifetime broadening. Modeling the dispersion within a Hubbard-Heisenberg approach, we extract the interaction strength and compare it with other cuprate systems. Our results establish a cuprate universal relation between electron-electron interaction and magnon zone-boundary dispersion. Superconductivity seems to be optimal at intermediate correlation strength, suggesting an optimal balance between localization and itinerancy.

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

Orbital Separation of Charge Order and Superconductivity in La$_{2-x}$Sr$_{x}$CuO$_4$

We report a combined x-ray absorption (XAS) and oxygen $K$-edge resonant inelastic x-ray scattering (RIXS) study of charge order in underdoped La$_{2-x}$Sr$_x$CuO$_4$ ($x=0.125$) under uniaxial $c$-axis pressure. We find that compressive $c$-axis strain modifies the charge order only within the superconducting state, with a striking polarization dependence: suppression in the $d_{x^2-y^2}$ channel and enhancement in the $d_{z^2}$ channel. X-ray absorption spectra reveal concomitant strain-induced modifications of the oxygen pre-edge and upper Hubbard band, consistent with increased $d_{z^2}$ orbital admixture. Our results suggest that $c$-axis pressure drives an orbital separation between superconductivity, rooted in $d_{x^2-y^2}$ states, and charge order, which gradually shifts to the $d_{z^2}$ channel. This orbital separation reveals a way for superconductivity and charge order to coexist in the cuprates with minimal competition. Furthermore, it suggests that the multi-order phase diagram of La$_{2-x}$Sr$_{x}$CuO$_4$ cannot be realistically described within single band models usually used to describe cuprate physics.

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

Decoupling of Static and Dynamic Charge Correlations revealed by Uniaxial Strain in a Cuprate Superconductor

We use uniaxial strain in combination with ultra-high-resolution Resonant Inelastic X-ray Scattering (RIXS) at the oxygen K- and copper L3-edges to study the excitations stemming from the charge ordering wave vector in La1.875Sr0.125CuO4. By detwinning stripe ordering, we demonstrate that the optical phonon anomalies do not show any stripe anisotropy. The low-energy charge excitations also retain an in-plane four-fold symmetry. As such, we find that both phonon and charge excitations are decoupled entirely from the strength of static charge ordering. The almost isotropic character of charge excitations is indicative of a quantum critical behaviour and remains a possible source for the strange metal properties found in the normal state of cuprate superconductors.

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

A Study on Monte Carlo simulation of the radiation environment above GeV at the DAMPE orbit

The Dark Matter Particle Explorer (DAMPE) has been undergoing a stable on-orbit operation for more than 6 years and acquired observation of over 11 billion events. And a better understanding of the overall radiation environment on the DAMPE orbit is crucial for both simulation data production and flight data analysis. In this work, we study the radiation environment at the low Earth orbit and develop a simulation software package using the framework of ATMNC3, in which state-of-the-art full 3D models of the Earth&#39;s atmospheric and magnetic-field configurations is integrated. We consider in our Monte Carlo procedure event-by-event propagation of the cosmic rays in the geomagnetic field and their interaction with the Earth&#39;s atmosphere, focusing on the particles above GeV that are able to trigger the DAMPE data acquisition system. We compare the simulation results with the cosmic-ray electrons and positrons (CREs) flux measurements made by DAMPE. The overall agreement on both the spectral and angular distribution of the CREs flux demonstrates that our simulation is well established. Our software package could be of more general usage for the simulation of the radiation environment at the low Earth orbit of various altitudes.

preprint2022arXiv

Magnetic Excitations in Strained Infinite-layer Nickelate PrNiO2

Strongly correlated materials often respond sensitively to the external perturbations. In the recently discovered superconducting infinite-layer nickelates, the superconducting transition temperature can be dramatically enhanced via only ~1% compressive strain-tuning enabled by substrate design. However, the root of such enhancement remains elusive. While the superconducting pairing mechanism is still not settled, magnetic Cooper pairing - similar to the cuprates has been proposed. Using resonant inelastic x-ray scattering, we investigate the magnetic excitations in infinite-layer PrNiO2 thin films for different strain conditions. The magnon bandwidth of PrNiO2 shows only marginal response to strain-tuning, in sharp contrast to the striking enhancement of the superconducting transition temperature Tc in the doped superconducting samples. These results suggest the enhancement of Tc is not mediated by spin excitations and thus provide important empirics for the understanding of superconductivity in infinite-layer nickelates.

preprint2022arXiv

Uniaxial Pressure Induced Stripe Order Rotation in La$_{1.88}$Sr$_{0.12}$CuO$_4$

Static stripe order is detrimental to superconductivity. Yet, it has been proposed that transverse stripe fluctuations may enhance the inter-stripe Josephson coupling and thus promote superconductivity. Direct experimental studies of stripe dynamics, however, remain difficult. From a strong-coupling perspective, transverse stripe fluctuations are realized in the form of dynamic &#34;kinks&#34; -- sideways shifting stripe sections. Here, we show how modest uniaxial pressure tuning reorganizes directional kink alignment. Our starting point is La$_{1.88}$Sr$_{0.12}$CuO$_4$, where transverse kink ordering results in a rotation of stripe order away from the crystal axis. Application of mild uniaxial pressure changes the ordering pattern and pins the stripe order to the crystal axis. This reordering occurs at a much weaker pressure than that to detwin the stripe domains and suggests a rather weak transverse stripe stiffness. Weak spatial stiffness and transverse quantum fluctuations are likely key prerequisites for stripes to coexist with superconductivity.

preprint2021arXiv

Evidence for a square-square vortex lattice transition in a high-$T_\textrm{c}$ cuprate superconductor

Using sound velocity and attenuation measurements in high magnetic fields, we identify a new transition in the vortex lattice state of La$_{2-x}$Sr$_{x}$CuO$_4$ (LSCO). The transition, observed in magnetic fields exceeding 35 T and temperatures far below zero field $T_c$, is detected in the compression modulus of the vortex lattice, at a doping level $x=p=0.17$. Our theoretical analysis based on Eilenberger theory of vortex lattice shows that the transition corresponds to the long-sought 45 degrees rotation of the square vortex lattice, predicted to occur in $d$-wave superconductors near a van Hove singularity.

preprint2021arXiv

Unconventional transverse transport above and below the magnetic transition temperature in Weyl semimetal EuCd$_2$As$_2$

As exemplified by the growing interest in the quantum anomalous Hall effect, the research on topology as an organizing principle of quantum matter is greatly enriched from the interplay with magnetism. In this vein, we present a combined electrical and thermoelectrical transport study on the magnetic Weyl semimetal EuCd$_2$As$_2$. Unconventional contribution to the anomalous Hall and anomalous Nernst effects were observed both above and below the magnetic transition temperature of EuCd$_2$As$_2$, indicating the existence of significant Berry curvature. EuCd$_2$As$_2$ represents a rare case in which this unconventional transverse transport emerges both above and below the magnetic transition temperature in the same material. The transport properties evolve with temperature and field in the antiferromagnetic phase in a different manner than in the paramagnetic phase, suggesting different mechanisms to their origin. Our results indicate EuCd$_2$As$_2$ is a fertile playground for investigating the interplay between magnetism and topology, and potentially a plethora of topologically nontrivial phases rooted in this interplay.

preprint2020arXiv

Cooperative Interactions Govern the Fermiology of the Polar Metal Ca$_3$Ru$_2$O$_7$

The antiferromagnetic Ruddlesden-Popper ruthenate Ca$_3$Ru$_2$O$_7$ is a model polar metal, combining inversion symmetry breaking with metallic conductivity; however, its low temperature ($T < 48$ K) crystal structure and Fermi surface topology remain ambiguous despite numerous measurements and theoretical studies. Here we perform both first principles calculations with static correlations and angle resolved photoelectron spectroscopy experiments to construct a complete model of Ca$_3$Ru$_2$O$_7$, reconciling inconsistencies among interpretations of electrical transport, thermopower measurements, and momentum- and energy-resolved band dispersions. The solution relies on treating the interplay among Coulomb repulsion, magnetic ordering, spin-orbit interactions, and the RuO$_6$ octahedral degrees-of-freedom on equal footing. For temperatures $30<T < 48$ K, we propose weak electron-electron interactions produce a symmetry-preserving metal-semimetal transition with Weyl nodes in proximity to the Fermi level, whereas a new orthorhombic $Pn2_1a$ structure emerges for $T<30$ K, exhibiting charge and spin density waves from enhanced Coulombic interactions.

preprint2020arXiv

Disentangling Intertwined Quantum States in a Prototypical Cuprate Superconductor

Spontaneous symmetry breaking constitutes a paradigmatic classification scheme of matter. However, broken symmetry also entails domain degeneracy that often impedes identification of novel low symmetry states. In quantum matter, this is additionally complicated by competing intertwined symmetry breaking orders. A prime example is that of unconventional superconductivity and density-wave orders in doped cuprates in which their respective symmetry relation remains a key question. Using uniaxial pressure as a domain-selective stimulus in combination with x-ray diffraction, we unambiguously reveal that the fundamental symmetry of the charge order in the prototypical cuprate La$_{1.88}$Sr$_{0.12}$CuO$_4$ is characterized by uniaxial stripes. We further demonstrate the direct competition of this stripe order with unconventional superconductivity via magnetic field tuning. The stripy nature of the charge-density-wave state established by our study is a prerequisite for the existence of a superconducting pair-density-wave -- a theoretical proposal that clarifies the interrelation of intertwined quantum phases in unconventional superconductors -- and paves the way for its high-temperature realization.

preprint2020arXiv

Gaia-DR2 extended kinematical maps. Part II: Dynamics in the Galactic disk explaining radial and vertical velocities

Context: In our Paper I, by using statistical deconvolution methods, extended kinematics maps of Gaia-DR2 data have been produced in a range of heliocentric distances that are a factor of two to three larger than those analyzed previously by the Gaia Collaboration with the same data. It added the range of Galactocentric distances between 13 kpc and 20 kpc to the previous maps. Aims: Here, we investigate the dynamical effects produced by different mechanisms that can explain the radial and vertical components of these extended kinematic maps, including a decomposition of bending and breathing of the vertical components. This paper as a whole tries to be a compendium of different dynamical mechanisms whose predictions can be compared to the kinematic maps. Methods: Using analytical methods or simulations, we are able to predict the main dynamical factors and compare them to the predictions of the extended kinematic maps of Gaia-DR2. Results: The gravitational influence of Galactic components that are different from the disk, such as the long bar or bulge, the spiral arms, or a tidal interaction with Sagittarius dwarf galaxy, may explain some features of the velocity maps, especially in the inner parts of the disk. However, they are not sufficient in explaining the most conspicuous gradients in the outer disk. Vertical motions might be dominated by external perturbations or mergers, although a minor component may be due to a warp whose amplitude evolves with time. Here, we show with two different methods, which analyze the dispersion of velocities, that the mass distribution of the disk is flared. Despite these partial explanations, the main observed features can only be explained in terms of out-of-equilibrium models, which are either due to external perturbers or to the fact that the disk has not had time to reach equilibrium since its formation.

preprint2020arXiv

High-Temperature Charge-Stripe Correlations in La$_{1.675}$Eu$_{0.2}$Sr$_{0.125}$CuO$_4$

We use resonant inelastic x-ray scattering to investigate charge-stripe correlations in La$_{1.675}$Eu$_{0.2}$Sr$_{0.125}$CuO$_4$. By differentiating elastic from inelastic scattering, it is demonstrated that charge-stripe correlations precede both the structural low-temperature tetragonal phase and the transport-defined pseudogap onset. The scattering peak amplitude from charge stripes decays approximately as $T^{-2}$ towards our detection limit. The in-plane integrated intensity, however, remains roughly temperature independent. Therefore, although the incommensurability shows a remarkably large increase at high temperature, our results are interpreted via a single scattering constituent. In fact, direct comparison to other stripe-ordered compounds (La$_{1.875}$Ba$_{0.125}$CuO$_4$, La$_{1.475}$Nd$_{0.4}$Sr$_{0.125}$CuO$_4$ and La$_{1.875}$Sr$_{0.125}$CuO$_4$) suggests a roughly constant integrated scattering intensity across all these compounds. Our results therefore provide a unifying picture for the charge-stripe ordering in La-based cuprates. As charge correlations in La$_{1.675}$Eu$_{0.2}$Sr$_{0.125}$CuO$_4$ extend beyond the low-temperature tetragonal and pseudogap phase, their emergence heralds a spontaneous symmetry breaking in this compound.

preprint2020arXiv

Magnetism and anomalous transport in the Weyl semimetal PrAlGe: Possible route to axial gauge fields

In magnetic Weyl semimetals, where magnetism breaks time-reversal symmetry, large magnetically sensitive anomalous transport responses are anticipated that could be useful for topological spintronics. The identification of new magnetic Weyl semimetals is therefore in high demand, particularly since in these systems Weyl node configurations may be easily modified using magnetic fields. Here we explore experimentally the magnetic semimetal PrAlGe, and unveil a direct correspondence between easy-axis Pr ferromagnetism and anomalous Hall and Nernst effects. With sizes of both the anomalous Hall conductivity and Nernst effect in good quantitative agreement with first principles calculations, we identify PrAlGe as a system where magnetic fields can connect directly to Weyl nodes via the Pr magnetization. Furthermore, we find the predominantly easy-axis ferromagnetic ground state co-exists with a low density of nanoscale textured magnetic domain walls. We describe how such nanoscale magnetic textures could serve as a local platform for tunable axial gauge fields of Weyl fermions.

preprint2020arXiv

Mapping the Galactic disk with the LAMOST and Gaia Red clump sample. VI: An evidence for the long lived non-steady warp of non-gravitational scenarios

By combining LAMOST DR4 and Gaia DR2 common red clump stars with age and proper motion, we analyze the amplitude evolution of the stellar warp independently of any assumption with a simple model. The greatest height of the warp disk increases with Galactocentric distance in different populations and it is dependent on the age: the younger stellar populations exhibit stronger warp features than the old ones, accompanied with the warp amplitude $γ({\rm age})$ decreasing with age and its first derivative $\dotγ ({\rm age})$ is different from zero. The azimuth of line of nodes $ϕ_w$ is stable at $-$5 degree without clear time evolution, which perfectly confirms some of previous works. All these self-consistent evidences support that our Galactic warp should most likely be a long lived, but non-steady structure and not a transient one, which is supporting the warp is originated from gas infall onto the disk or other hypotheses that suppose that the warp mainly affects to the gas, and consequently younger populations tracing the gas are stronger than the older ones. In other words, the Galactic warp is induced by the non-gravitational interaction over the disk models.

preprint2020arXiv

Multiple- to single-gap superconductivity crossover in Nb$_5$Ir$_{3-x}$Pt$_x$O alloys

By using mostly the muon-spin rotation/relaxation ($μ$SR) technique, we investigate the superconductivity (SC) of Nb$_5$Ir$_{3-x}$Pt$_x$O ($x = 0$ and 1.6) alloys, with $T_c = 10.5$ K and 9.1 K, respectively. At a macroscopic level, their superconductivity was studied by electrical resistivity, magnetization, and specific-heat measurements. In both compounds, the electronic specific heat and the low-temperature superfluid density data suggest a nodeless SC. The superconducting gap value and the specific heat discontinuity at $T_c$ are larger than that expected from the Bardeen-Cooper-Schrieffer theory in the weak-coupling regime, indicating strong-coupling superconductivity in the Nb$_5$Ir$_{3-x}$Pt$_x$O family. In Nb$_5$Ir$_3$O, multigap SC is evidenced by the field dependence of the electronic specific heat coefficient and the superconducting Gaussian relaxation rate, as well as by the temperature dependence of the upper critical field. Pt substitution suppresses one of the gaps, and Nb$_5$Ir$_{1.4}$Pt$_{1.6}$O becomes a single-gap superconductor. By combining our extensive experimental results, we provide evidence for a multiple- to single-gap SC crossover in the Nb$_5$Ir$_{3-x}$Pt$_x$O family.

preprint2020arXiv

Oxide Fermi liquid universality revealed by electron spectroscopy

We present a combined soft x-ray and high-resolution vacuum-ultraviolet angle-resolved photoemission spectroscopy study of the electron-overdoped cuprate Pr$_{1.3-x}$La$_{0.7}$Ce$_{x}$CuO$_4$ (PLCCO). Demonstration of its highly two-dimensional band structure enabled precise determination of the in-plane self-energy dominated by electron-electron scattering. Through analysis of this self-energy and the Fermi-liquid cut-off energy scale, we find -- in contrast to hole-doped cuprates -- a momentum isotropic and comparatively weak electron correlation in PLCCO. Yet, the self-energies extracted from multiple oxide systems combine to demonstrate a logarithmic divergent relation between the quasiparticle scattering rate and mass. This constitutes a spectroscopic version of the Kadowaki-Woods relation with an important merit -- the demonstration of Fermi liquid quasiparticle lifetime and mass being set by a single energy scale.

preprint2020arXiv

Spatially Inhomogeneous Competition between Superconductivity and the Charge Density Wave in YBa$_2$Cu$_3$O$_{6.67}$

The charge density wave in the high-temperature superconductor YBa$_2$Cu$_3$O$_{7-x}$ (YBCO) is now known to have two different ordering tendencies differentiated by their $c$-axis correlations. These correspond to ferro- (F-CDW) and antiferro- (AF-CDW) couplings between CDW in neighbouring CuO$_2$ bilayers. This discovery has prompted a number of fundamental questions. For example, how does superconductivity adjust to two competing orders and are either of these orders responsible for the electronic reconstruction? Here we use high-energy x-ray diffraction to study YBa$_2$Cu$_3$O$_{6.67}$ as a function of magnetic field and temperature. We show that regions of the sample with F-CDW correlations suppress superconductivity more strongly than those with AF-CDW correlations. This implies that an inhomogeneous superconducting state exists, in which some regions show a weak or fragile form of superconductivity. By comparison of F-CDW and AF-CDW correlation lengths, it is furthermore concluded that F-CDW ordering is sufficiently long-range to modify the electronic structure. Our study thus suggests that F-CDW correlations have an important impact on superconducting and normal state properties of underdoped YBCO.

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.

preprint2019arXiv

Mapping the Galactic disk with the LAMOST and Gaia Red clump sample. II. 3D asymmetrical kinematics of mono-age populations in the disk between 6$-14$ \,kpc

We perform analysis of the three-dimensional kinematics of Milky Way disk stars in mono-age populations. We focus on stars between Galactocentric distances of $R=6$ and 14 \,kpc, selected from the combined LAMOST DR4 red clump giant stars and Gaia DR2 proper motion catalogue. We confirm the 3D asymmetrical motions of recent works and provide time tagging of the Galactic outer disk asymmetrical motions near the anticenter direction out to Galactocentric distances of 14\,kpc. Radial Galactocentric motions reach values up to 10 km s$^{-1}$, depending on the age of the population, and present a north-south asymmetry in the region corresponding to density and velocity substructures that were sensitive to the perturbations in the early 6 \,Gyr. After that time, the disk stars in this asymmetrical structure have become kinematically hotter, and are thus not sensitive to perturbations, and we find the structure is a relatively younger population. With quantitative analysis, we find stars both above and below the plane at $R\gtrsim 9$ kpc that exhibit bending mode motions of which the sensitive duration is around 8 \,Gyr. We speculate that the in-plane asymmetries might not be mainly caused by a fast rotating bar, intrinsically elliptical outer disk, secular expansion of the disk, or streams. Spiral arm dynamics, out-of-equilibrium models, minor mergers or others are important contributors. Vertical motions might be dominated by bending and breathing modes induced by complicated inner or external perturbers. It is likely that many of these mechanisms are coupled together.