Researcher profile

Y. Su

Y. Su contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

32 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

Evidence for the rare decay $B^{+} \to \bar Λp μ^{+} μ^{-}$

A search for the rare decay $B^{+} \to \bar Λp μ^{+} μ^{-}$ is performed using proton-proton collision data recorded by the LHCb experiment at a center-of-mass energy of $\sqrt{s}= 13$ TeV, corresponding to an integrated luminosity of 5.4 fb$^-1$. An excess of events is found with respect to the background-only expectation, with a signal significance of 3.5 standard deviations, in the low invariant-mass region of $m(\bar Λp)<2.8$ GeV/$c^2$. The branching fraction is measured to be ${\cal B}_{\rm low}({B^{+}}\to {\bar Λp μ^{+} μ^{-}})=\left(1.70 ^{+0.65}_{-0.56}(\rm stat) \pm 0.17(\rm syst) \pm 0.14(\rm ext)\right) \times 10^{-8}$, where the last uncertainty is due to external inputs on ${\cal B}(B^+\to{J/ψ\bar Λp})\times {\cal B}({J/ψ}\to{μ^+μ^-})$. With no significant signal observed in the high $m(\bar Λp)$ region above 2.8 GeV/$c^2$, an upper limit is set to be ${\cal B}_{\rm high}(B^{+}\to {\bar Λp μ^{+} μ^{-}})<2.8\,(3.7) \times 10^{-9}$ at the $90\%$ ($95\%$) confidence level.

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&#39; 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 transverse $Λ$ and $\barΛ$ hyperon polarization in $p$Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV

The transverse polarization of $Λ$ and $\barΛ$ hyperons is measured in $p$Pb collisions collected by the LHCb experiment at a nucleon-nucleon center-of-mass energy of $5.02$ TeV. The polarization is averaged over hyperon transverse momentum in the ranges $0.15 < p_{T} < 6.00 $ GeV/$c$, and Feynman-$x$ in the ranges $0.005 < x_{F} < 0.040$ (forward region) and $-0.10 < x_{F} < -0.01$ (backward region) defined relative to the proton beam direction. The transverse polarization is found to be compatible with zero for both $Λ$ and $\barΛ$ hyperons. The results are also measured as a function of $p_{T}$ and $x_{F}$ with no significant dependence on these variables observed. The results are compared with previous experimental measurements at different center-of-mass energies and collision environments.

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

Search for lepton-number-violating $B^-\to D^{(*)+}μ^-μ^-$ decays

A search is performed for lepton-number-violating $B^-\to D^{(*)+}μ^-μ^-$ decays, using data collected by the LHCb experiment in proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 5.4 fb$^{-1}$. No significant signal is observed, and upper limits are set on the branching fractions, ${\cal B}(B^-\to D^{+}μ^-μ^-) < 4.6 \times 10^{-8}$ and ${\cal B}(B^-\to D^{*+}μ^-μ^-) < 5.9 \times 10^{-8}$, at the 95% confidence level.

preprint2026arXiv

Study of charm mixing and CP violation with $D^0\to K^\pmπ^\mpπ^\pmπ^\mp$ decays

A study of charm mixing and CP violation in $D^0\to K^\pmπ^\mpπ^\pmπ^\mp$ decays is performed using data collected by the LHCb experiment in proton-proton collisions from 2015 to 2018, corresponding to an integrated luminosity of 6$\text{fb}^{-1}$. The ratio of promptly produced $D^0\to K^+π^- π^+π^-$ to $D^0\to K^-π^+ π^-π^+$ decay rates is measured as a function of $D^0$ decay time, both inclusive over phase space and in bins of phase space. Taking external inputs for the $D^0 -\overline{D}^0$ mixing parameters $x$ and $y$ allows constraints to be obtained on the hadronic parameters of the charm decay. When combined with previous measurements from charm-threshold experiments and at LHCb, improved knowledge is obtained for these parameters, which is valuable for studies of the angle $γ$ of the Unitarity Triangle. An alternative analysis is also performed, in which external inputs are taken for the hadronic parameters, and the mixing parameters are determined, including $Δx$ and $Δy$, which are nonzero in the presence of CP violation. It is found that $x=\left(0.85^{+0.15}_{-0.24}\right)\%$, $y=\left( 0.21^{+0.29}{-0.27} \right)\%$, $Δx=\left( -0.02\pm {0.04} \right)\% $ and $Δy=\left( 0.02^{+0.04}_{-0.03} \right)\%$. These results are consistent with previous measurements and the hypothesis of CP conservation.

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

Coherent photoproduction of $ρ^0, ω$ and excited vector mesons in ultraperipheral PbPb collisions

The invariant-mass distribution for the coherent photoproduction of dipions in ultraperipheral PbPb collisions is measured using data, corresponding to an integrated luminosity of $224.6 \pm 9.6 μ$b$^{-1}$, collected by the LHCb experiment in 2018 at a nucleon-nucleon centre-of-mass energy $\sqrt{s_{\rm NN}}=5.02$ TeV. In the mass range from 400 to 1200 MeV, the results are consistent with previous experiments, with the spectrum dominated by the $ρ^0$ meson, which interferes with a nonresonant component, together with a smaller $ω$ meson contribution. In an extended mass range up to 2300 MeV, models previously used do not fit the data and a consistent description requires the introduction of two resonances at masses of $1350\pm20$ MeV and $1790\pm20$ MeV with widths of about 300 MeV. The cross-section for each meson is measured differentially in twelve bins of rapidity from 2.05 to 4.90. The $ρ^0$ cross-section increases with rapidity from about 400 to 600 mb and is measured with a typical precision of 8\%, while the cross-section times branching fraction for the $ω,ρ^\prime$ and $ρ^{\prime\prime}$, with the statistical precision of the data, do not have a pronounced rapidity dependence and are between 0.5 and 1.5 mb, with uncertainties up to 30\%. A large nuclear suppression is observed for the $ρ^0$ meson compared to expectations based on photoproduction on the proton that use the impulse approximation. Significant suppression is also observed compared to that predicted by elastic scattering described in the Glauber approach, or with the addition of inelastic scattering in a Gribov--Glauber model.

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.

preprint2023arXiv

Longitudinal and transverse spin transfer to $Λ$ and $\overlineΛ$ hyperons in polarized $p$+$p$ collisions at $\sqrt{s} = 200$ GeV

The longitudinal and transverse spin transfers to $Λ$ ($\overlineΛ$) hyperons in polarized proton-proton collisions are expected to be sensitive to the helicity and transversity distributions, respectively, of (anti-)strange quarks in the proton, and to the corresponding polarized fragmentation functions. We report improved measurements of the longitudinal spin transfer coefficient, $D_{LL}$, and the transverse spin transfer coefficient, $D_{TT}$, to $Λ$ and $\overlineΛ$ in polarized proton-proton collisions at $\sqrt{s}$ = 200 GeV by the STAR experiment at RHIC. The data set includes longitudinally polarized proton-proton collisions with an integrated luminosity of 52 pb$^{-1}$, and transversely polarized proton-proton collisions with a similar integrated luminosity. Both data sets have about twice the statistics of previous results and cover a kinematic range of $|η_{Λ(\overlineΛ)}|$ $<$ 1.2 and transverse momentum $p_{T,{Λ(\overlineΛ)}}$ up to 8 GeV/$c$. We also report the first measurements of the hyperon spin transfer coefficients $D_{LL}$ and $D_{TT}$ as a function of the fractional jet momentum $z$ carried by the hyperon, which can provide more direct constraints on the polarized fragmentation functions.

preprint2022arXiv

AGN feedback duty cycle in Planck SZ selected clusters using Chandra observations

We present a systematic study of X-ray cavities using archival Chandra observations of nearby galaxy clusters selected by their Sunyaev-Zel&#39;dovich (SZ) signature in the Planck survey, which provides a nearly unbiased mass-selected sample to explore the entire AGN feedback duty cycle. Based on X-ray image analysis, we report that 30 of the 164 clusters show X-ray cavities, which corresponds to a detection fraction of 18%. After correcting for spatial resolution to match the high-$z$ SPT-SZ sample, the detection fraction decreases to 9%, consistent with the high-z sample, hinting that the AGN feedback has not evolved across almost 8 Gyrs. Our finding agrees with the lack of evolution of cool-core clusters fraction. We calculate the cavity power, P_{\rm cav}, and find that most systems of our sample have enough AGN heating to offset the radiative losses of the intracluster medium.

preprint2022arXiv

Nonlinear co-generation of graphene plasmons for optoelectronic logic operations

Surface plasmons in graphene provide a compelling strategy for advanced photonic technologies thanks to their tight confinement, fast response and tunability. Recent advances in the field of all optical generation of graphene plasmons in planar waveguides offer a promising method for high speed signal processing in nanoscale integrated optoelectronic devices. Here, we use two counter propagating frequency combs with temporally synchronized pulses to demonstrate deterministic all optical generation and electrical control of multiple plasmon polaritons, excited via difference frequency generation (DFG). Electrical tuning of a hybrid graphene fibre device offers a precise control over the DFG phase matching, leading to tunable responses of the graphene plasmons at different frequencies across a broadband (0 - 50 THz) and provides a powerful tool for high speed logic operations. Our results offer insights for plasmonics on hybrid photonic devices based on layered materials and pave the way to high speed integrated optoelectronic computing circuits.

preprint2022arXiv

Pinch points and half-moons in dipolar-octupolar Nd$_2$Hf$_2$O$_7$

While it is established that the pinch point scattering pattern in spin ice arises from an emergent coulomb phase associated with magnetic moment that is divergence-free, more complex Hamiltonians can introduce a divergence-full part. If these two parts remain decoupled, they give rise to the co-existence of distinct features. Here we show that the moment in ${\rm Nd_2Hf_2O_7}$ forms a static long-range ordered ground state, a flat, gapped pinch point excitation and dispersive excitations. These results confirm recent theories which predict that the dispersive modes, which arise from the divergence-full moment, host a pinch point pattern of their own, observed experimentally as `half-moons&#39;.

preprint2022arXiv

The cycle of metals in the infalling elliptical galaxy NGC 1404

Hot atmospheres pervading galaxy clusters, groups, and early-type galaxies are rich in metals, produced during epochs and diffused via processes that are still to be determined. While this enrichment has been routinely investigated in clusters, metals in lower mass systems are more challenging to probe with standard X-ray exposures and spectroscopy. In this paper, we focus on very deep XMM-Newton ($\sim$350 ks) observations of NGC 1404, a massive elliptical galaxy experiencing ram-pressure stripping of its hot atmosphere while infalling toward the centre of the Fornax cluster, with the aim to derive abundances through its hot gas extent. Importantly, we report the existence of a new fitting bias - the &#34;double Fe bias&#34; - leading to an underestimate of the Fe abundance when two thermal components cannot realistically model the complex temperature structure present in the outer atmosphere of the galaxy. Contrasting with the &#39;&#39;metal conundrum&#39;&#39; seen in clusters, the Fe and Mg masses of NGC 1404 are measured 1-2 orders of magnitude below what stars and supernovae could have reasonably produced and released. In addition, we note the remarkable Solar abundance ratios of the galaxy&#39;s halo, different from its stellar counterpart but similar to the chemical composition of the ICM of rich clusters. Completing the clusters regime, all these findings provide additional support toward a scenario of early enrichment, at play over two orders of magnitude in mass. A few peculiar and intriguing features, such as a possible double metal peak as well as an apparent ring of enhanced Si near the galaxy core, are also discussed.

preprint2022arXiv

The Infrared Database of Extragalactic Observables from Spitzer. -- II. The Database & The Diagnostic Power of Crystalline Silicate Features in Galaxy Spectra

We present the Infrared Database of Extragalactic Observables from Spitzer (IDEOS), a homogeneous, publicly available, database of 77 fitted mid-infrared observables in the 5.4-36um range, comprising measurements for 3335 galaxies observed in the low-resolution staring mode of the Infrared Spectrometer onboard the Spitzer Space Telescope. Among the included observables are PAH fluxes and their equivalent widths, the strength of the 9.8um silicate feature, emission line fluxes, solid-state features, rest frame continuum fluxes, synthetic photometry, and a mid-infrared spectral classification. The IDEOS spectra were selected from the Cornell Atlas of Spitzer-IRS Sources. To our surprise we have detected at a >95% confidence level crystalline silicates in the spectra of 786 IDEOS galaxies. The detections range from single band detections to detections of all fitted crystalline bands (16, 19, 23, 28 and 33um). We find the strength of the crystalline silicate bands to correlate with the amorphous silicate strength, and the change from an emission to an absorption feature to occur at higher obscuration as the wavelength of the crystalline silicate band is longer. These observed characteristics are consistent with an origin for the amorphous and crystalline silicate features in a centrally heated dust geometry, either an edge-on disk or a cocoon. We find the 23 and 33um crystalline silicate bands to be well-suited to classify the obscuration level of galactic nuclei, even in the presence of strong circumnuclear star formation. Based on our detection statistics, we conclude that crystalline silicates are a common component of the interstellar medium of galactic nuclei.

preprint2021arXiv

Distinct magnetic ground states of $R_2$ZnIrO$_6$ ($R$ = La and Nd) determined by neutron powder diffraction

Double perovskite iridates $A_2$ZnIrO$_6$ ($A$ = alkaline or lanthanide) show complex magnetic behaviors ranging from weak ferromagnetism to successive antiferromagnetic transitions. Here we report the static ($dc$) and dynamic ($ac$) magnetic susceptibility, and neutron powder diffraction measurements for $A$ = La and Nd compounds to elucidate the magnetic ground state. Below 10~K, the $A$ = La compound is best described as canted iridium moments in an antiferromagnet arrangement with a propagation vector \textbf{k} = 0 and a net ferromagnetic component along the $c$-axis. On the other hand, Nd$_2$ZnIrO$_6$ is described well as an antiferromagnet with a propagation vector \textbf{k} = (1/2~1/2~0) below $T_\mathrm{N} \sim$ 17 K. Scattering from both the Nd and Ir magnetic sublattices were required to describe the data and both were found to lie almost completely within the $ab$-plane. $Dc$ susceptibility revealed a bifurcation between the zero-field-cooled and field-cooled curves below $\sim$13 K in Nd$_2$ZnIrO$_6$. A glassy state was ruled out by $ac$ susceptibility but detailed magnetic isotherms revealed the opening of the loop below 13~K. These results suggest a delicate balance exists between the Dzyaloshinskii-Moriya, crystal field schemes, and $d$-$f$ interaction in this series of compounds.

preprint2020arXiv

Cation controlled wetting properties of vermiculite membranes and its potential for fouling resistant oil-water separation

The surface free energy is one of the most fundamental properties of solids, hence, manipulating the surface energy and thereby the wetting properties of solids, has tremendous potential for various physical, chemical, biological as well as industrial processes. Typically, this is achieved by either chemical modification or by controlling the hierarchical structures of surfaces. Here we report a phenomenon whereby the wetting properties of vermiculite laminates are controlled by the hydrated cations on the surface and in the interlamellar space. We find that by exploiting this mechanism, vermiculite laminates can be tuned from superhydrophillic to hydrophobic simply by exchanging the cations; hydrophilicity decreases with increasing cation hydration free energy, except for lithium. Lithium, which has a higher hydration free energy than potassium, is found to provide a superhydrophilic surface due to its anomalous hydrated structure at the vermiculite surface. Building on these findings, we demonstrate the potential application of superhydrophilic lithium exchanged vermiculite as a thin coating layer on microfiltration membranes to resist fouling, and thus, we address a major challenge for oil-water separation technology.

preprint2019arXiv

Neutron scattering study of commensurate magnetic ordering in single crystal CeSb$_2$

Temperature and field-dependent magnetization $M(H,T)$ measurements and neutron scattering study of a single crystal CeSb$_2$ are presented. Several anomalies in the magnetization curves have been confirmed at low magnetic field, i.e., 15.6 K, 12 K, and 9.8 K. These three transitions are all metamagnetic transitions (MMT), which shift to lower temperatures as the magnetic field increases. The anomaly at 15.6 K has been suggested as paramagnetic (PM) to ferromagnetic (FM) phase transition. The anomaly located at around 12 K is antiferromagnetic-like transition, and this turning point will clearly split into two when the magnetic field $H\geq0.2$ T. Neutron scattering study reveals that the low temperature ground state of CeSb$_2$ orders antiferromagnetically with commensurate propagation wave vectors $\textbf{k}=(-1,\pm1/6,0)$ and $\textbf{k}=(\pm1/6,-1,0)$, with Néel temperature $T_N\sim9.8$ K. This transition is of first-order, as shown in the hysteresis loop observed by the field cooled cooling (FCC) and field cooled warming (FCW) processes.