Source author record

R. Zhou

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

ResearcherUnclaimed source record

hep-lat hep-ph cond-mat.supr-con cond-mat.str-el cond-mat.mtrl-sci Artificial Intelligence cond-mat.mes-hall nucl-th physics.optics

Catalog footprint

What is connected

70works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

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

DESI Data Release 1: Stellar Catalogue

In this paper we present the stellar Value-Added Catalogue (VAC) based on the DESI Data Release 1. This VAC contains stellar parameter, abundance and radial velocity measurements for more than 4 million stars. It also contains, for the first time, measurements from individual epochs for more than a million stars with at least two observations. The main contribution to the catalogue comes from the bright program of the main survey, which includes $\sim $2.5 million stars, and the backup program, which includes $\sim $ 1 million stars. The combined magnitude range for the stars in the catalogue extends from ${\it Gaia}$ G $\sim 12$ to G $\sim 21$. For the magnitude range $17.5<G<21$ this catalogue represents a factor of 10 increase in the number of stars with radial velocity and abundance measurements compared to existing surveys. Despite DESI's resolution (R $\sim 2500-5000$), the median radial velocity uncertainty for stars in the catalogue is better than 1 km s$^{-1}$. The stellar parameters and abundances of stars in DESI are measured by two independent pipelines, and after applying a temperature-dependent calibration, [Fe/H] abundances of high signal-to-noise stars are accurate to better than $\sim$ 0.1 dex when compared to high-resolution surveys. The catalogue probes different Galactic components including a particularly large number of distant stars: tens of thousands of stars further than 10 kpc, and thousands further than 50 kpc. The catalogue also contains several thousand extremely metal-poor stars with ${\rm [Fe/H]}<-3$. The released sample of stars includes measurements for thousands of stars that are members of dwarf galaxies, open and globular clusters as well as members of several dozen stellar streams. The next public DESI data release is expected in less than two years and will contain three times as many stars as DR1.

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

Evidence for particle acceleration approaching PeV energies in the W51 complex

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2025arXiv

Counterpart identification and classification for eRASS1 and characterisation of the AGN content

[abridged] Accurately accounting for the AGN phase in galaxy evolution requires a large, clean AGN sample. This is now possible with SRG/eROSITA. The public Data Release 1 (DR1, Jan 31, 2024) includes 930,203 sources from the Western Galactic Hemisphere. The data enable the selection of a large AGN sample and the discovery of rare sources. However, scientific return depends on accurate characterisation of the X-ray emitters, requiring high-quality multiwavelength data. This paper presents the identification and classification of optical and infrared counterparts to eRASS1 sources using Gaia DR3, CatWISE2020, and Legacy Survey DR10 (LS10) with the Bayesian NWAY algorithm and trained priors. Sources were classified as Galactic or extragalactic via a Machine Learning model combining optical/IR and X-ray properties, trained on a reference sample. For extragalactic LS10 sources, photometric redshifts were computed using Circlez. Within the LS10 footprint, all 656,614 eROSITA/DR1 sources have at least one possible optical counterpart; about 570,000 are extragalactic and likely AGN. Half are new detections compared to AllWISE, Gaia, and Quaia AGN catalogues. Gaia and CatWISE2020 counterparts are less reliable, due to the surveys shallowness and the limited amount of features available to assess the probability of being an X-ray emitter. In the Galactic Plane, where the overdensity of stellar sources also increases the chance of associations, using conservative reliability cuts, we identify approximately 18,000 Gaia and 55,000 CatWISE2020 extragalactic sources. We release three high-quality counterpart catalogues, plus the training and validation sets, as a benchmark for the field. These datasets have many applications, but in particular empower researchers to build AGN samples tailored for completeness and purity, accelerating the hunt for the Universes most energetic engines.

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

Pressure-enhanced $f$-electron orbital weighting in UTe2 mapped by quantum interferometry

The phase landscape of UTe$_2$ features a remarkable diversity of superconducting phases under applied pressure and magnetic field. Recent quantum oscillation studies at ambient pressure have revealed the quasi-2D Fermi surface of this material. However, the pressure-dependence of the Fermi surface remains an open question. Here we track the evolution of the UTe$_2$ Fermi surface as a function of pressure up to 19.5 kbar by measuring quantum interference oscillations. We find that in sufficient magnetic field to suppress both superconductivity at low pressures and incommensurate antiferromagnetism at higher pressures, the quasi-2D Fermi surface found at ambient pressure smoothly connects to that at 19.5 kbar, with no signs of a reconstruction over this pressure interval. We observe a smooth increase in oscillatory frequency with increasing pressure, indicating that the warping of the cylindrical Fermi sheets continuously increases with pressure. By computing a tight-binding model, we show that this enhanced warping indicates increased $f$-orbital contribution at the Fermi level - up to and beyond the critical pressure at which superconductivity is truncated. These findings highlight the value of high-pressure quantum interference measurements as a new probe of the electronic structure in heavy fermion materials.

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.

preprint2025arXiv

The Compilation and Validation of the Spectroscopic Redshift Catalogs for the DESI-COSMOS and DESI-XMMLSS Fields

Over several dedicated programs that include targets beyond the main cosmological samples, the Dark Energy Spectroscopic Instrument (DESI) collected spectra for 304,970 unique objects in two fields centered on the COSMOS and XMM-LSS fields. In this work, we develop spectroscopic redshift robustness criteria for those spectra, validate these criteria using visual inspection, and provide two custom Value-Added Catalogs with our redshift characterizations. With these criteria, we reliably classify 212,935 galaxies below z < 1.6, 9,713 quasars and 35,222 stars. As a critical element in characterizing the selection function, we provide the description of 70 different algorithms that were used to select these targets from imaging data. To facilitate joint imaging/spectroscopic analyses, we provide row-matched photometry from the Dark Energy Camera, Hyper-Suprime Cam, and public COSMOS2020 photometric catalogs. Finally, we demonstrate example applications of these large catalogs to photometric redshift estimation, cluster finding, and completeness studies.

preprint2023arXiv

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more.

preprint2022arXiv

All-sky search for gravitational wave emission from scalar boson clouds around spinning black holes in LIGO O3 data

This paper describes the first all-sky search for long-duration, quasi-monochromatic gravitational-wave signals emitted by ultralight scalar boson clouds around spinning black holes using data from the third observing run of Advanced LIGO. We analyze the frequency range from 20~Hz to 610~Hz, over a small frequency derivative range around zero, and use multiple frequency resolutions to be robust towards possible signal frequency wanderings. Outliers from this search are followed up using two different methods, one more suitable for nearly monochromatic signals, and the other more robust towards frequency fluctuations. We do not find any evidence for such signals and set upper limits on the signal strain amplitude, the most stringent being $\approx10^{-25}$ at around 130~Hz. We interpret these upper limits as both an "exclusion region" in the boson mass/black hole mass plane and the maximum detectable distance for a given boson mass, based on an assumption of the age of the black hole/boson cloud system.

preprint2022arXiv

First joint observation by the underground gravitational-wave detector, KAGRA, with GEO600

We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing run from April 7 to 20, 2020. We present the results of the joint analysis of the GEO--KAGRA data for transient gravitational-wave signals, including the coalescence of neutron-star binaries and generic unmodeled transients. We also perform dedicated searches for binary coalescence signals and generic transients associated with gamma-ray burst events observed during the joint run. No gravitational-wave events were identified. We evaluate the minimum detectable amplitude for various types of transient signals and the spacetime volume for which the network is sensitive to binary neutron-star coalescences. We also place lower limits on the distances to the gamma-ray bursts analysed based on the non-detection of an associated gravitational-wave signal for several signal models, including binary coalescences. These analyses demonstrate the feasibility and utility of KAGRA as a member of the global gravitational-wave detector network.

preprint2022arXiv

GWTC-2.1: Deep Extended Catalog of Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run

The second Gravitational-Wave Transient Catalog reported on 39 compact binary coalescences observed by the Advanced LIGO and Advanced Virgo detectors between 1 April 2019 15:00 UTC and 1 October 2019 15:00 UTC. We present GWTC-2.1, which reports on a deeper list of candidate events observed over the same period. We analyze the final version of the strain data over this period with improved calibration and better subtraction of excess noise, which has been publicly released. We employ three matched-filter search pipelines for candidate identification, and estimate the astrophysical probability for each candidate event. While GWTC-2 used a false alarm rate threshold of 2 per year, we include in GWTC-2.1, 1201 candidates that pass a false alarm rate threshold of 2 per day. We calculate the source properties of a subset of 44 high-significance candidates that have an astrophysical probability greater than 0.5. Of these candidates, 36 have been reported in GWTC-2. If the 8 additional high-significance candidates presented here are astrophysical, the mass range of events that are unambiguously identified as binary black holes (both objects $\geq 3M_\odot$) is increased compared to GWTC-2, with total masses from $\sim 14 M_\odot$ for GW190924_021846 to $\sim 182 M_\odot$ for GW190426_190642. The primary components of two new candidate events (GW190403_051519 and GW190426_190642) fall in the mass gap predicted by pair instability supernova theory. We also expand the population of binaries with significantly asymmetric mass ratios reported in GWTC-2 by an additional two events (the mass ratio is less than $0.65$ and $0.44$ at $90\%$ probability for GW190403_051519 and GW190917_114630 respectively), and find that 2 of the 8 new events have effective inspiral spins $χ_\mathrm{eff} > 0$ (at $90\%$ credibility), while no binary is consistent with $χ_\mathrm{eff} < 0$ at the same significance.

preprint2022arXiv

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets.

preprint2022arXiv

Possible star-of-David pattern charge density wave with additional modulation in the kagome superconductor CsV$_3$Sb$_5$

$A$V$_3$Sb$_5$ ($A$ = K, Rb, Cs) is a novel kagome superconductor coexisting with the charge density wave (CDW) order. Identifying the structure of the CDW order is crucial for understanding the exotic normal state and superconductivity in this system. Here, we report $^{51}$V nuclear magnetic resonance (NMR) and $^{121/123}$Sb nuclear quadrupole resonance (NQR) studies on kagome-metal CsV$_3$Sb$_5$. Below the CDW transition temperature $T_\textrm{CDW} \sim$ 98 K, an abrupt change of spectra was observed, indicating that the transition is of the first order. By further analysing the spectra, we find that the CDW order is commensurate. And most remarkably, the obtained experimental results suggest that the charge modulation of the CDW order is of star-of-David pattern and accompanied by an additional charge modulation in bulk below $T^* \sim$ 40 K. Our results revealing the unconventional CDW order provide new insights into $A$V$_3$Sb$_5$.

preprint2022arXiv

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO--Virgo data

We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo run in the detector frequency band $[10,2000]\rm~Hz$ have been used. No significant detection was found and 95$\%$ confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about $7.6\times 10^{-26}$ at $\simeq 142\rm~Hz$. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass -- boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC.

preprint2022arXiv

Search for continuous gravitational waves from 20 accreting millisecond X-ray pulsars in O3 LIGO data

Results are presented of searches for continuous gravitational waves from 20 accreting millisecond X-ray pulsars with accurately measured spin frequencies and orbital parameters, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. The search algorithm uses a hidden Markov model, where the transition probabilities allow the frequency to wander according to an unbiased random walk, while the $\mathcal{J}$-statistic maximum-likelihood matched filter tracks the binary orbital phase. Three narrow sub-bands are searched for each target, centered on harmonics of the measured spin frequency. The search yields 16 candidates, consistent with a false alarm probability of 30% per sub-band and target searched. These candidates, along with one candidate from an additional target-of-opportunity search done for SAX J1808.4$-$3658, which was in outburst during one month of the observing run, cannot be confidently associated with a known noise source. Additional follow-up does not provide convincing evidence that any are a true astrophysical signal. When all candidates are assumed non-astrophysical, upper limits are set on the maximum wave strain detectable at 95% confidence, $h_0^{95\%}$. The strictest constraint is $h_0^{95\%} = 4.7\times 10^{-26}$ from IGR J17062$-$6143. Constraints on the detectable wave strain from each target lead to constraints on neutron star ellipticity and $r$-mode amplitude, the strictest of which are $ε^{95\%} = 3.1\times 10^{-7}$ and $α^{95\%} = 1.8\times 10^{-5}$ respectively. This analysis is the most comprehensive and sensitive search of continuous gravitational waves from accreting millisecond X-ray pulsars to date.

preprint2022arXiv

Search of the Early O3 LIGO Data for Continuous Gravitational Waves from the Cassiopeia A and Vela Jr. Supernova Remnants

We present directed searches for continuous gravitational waves from the neutron stars in the Cassiopeia A (Cas A) and Vela Jr. supernova remnants. We carry out the searches in the LIGO data from the first six months of the third Advanced LIGO and Virgo observing run, using the Weave semi-coherent method, which sums matched-filter detection-statistic values over many time segments spanning the observation period. No gravitational wave signal is detected in the search band of 20--976 Hz for assumed source ages greater than 300 years for Cas A and greater than 700 years for Vela Jr. Estimates from simulated continuous wave signals indicate we achieve the most sensitive results to date across the explored parameter space volume, probing to strain magnitudes as low as ~$6.3\times10^{-26}$ for Cas A and ~$5.6\times10^{-26}$ for Vela Jr. at frequencies near 166 Hz at 95% efficiency.

preprint2022arXiv

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

preprint2021arXiv

All-sky search for long-duration gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

After the detection of gravitational waves from compact binary coalescences, the search for transient gravitational-wave signals with less well-defined waveforms for which matched filtering is not well-suited is one of the frontiers for gravitational-wave astronomy. Broadly classified into "short" $ \lesssim 1~$\,s and "long" $ \gtrsim 1~$\,s duration signals, these signals are expected from a variety of astrophysical processes, including non-axisymmetric deformations in magnetars or eccentric binary black hole coalescences. In this work, we present a search for long-duration gravitational-wave transients from Advanced LIGO and Advanced Virgo's third observing run from April 2019 to March 2020. For this search, we use minimal assumptions for the sky location, event time, waveform morphology, and duration of the source. The search covers the range of $2~\text{--}~ 500$~s in duration and a frequency band of $24 - 2048$ Hz. We find no significant triggers within this parameter space; we report sensitivity limits on the signal strength of gravitational waves characterized by the root-sum-square amplitude $h_{\mathrm{rss}}$ as a function of waveform morphology. These $h_{\mathrm{rss}}$ limits improve upon the results from the second observing run by an average factor of 1.8.

preprint2021arXiv

All-sky search for short gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

This paper presents the results of a search for generic short-duration gravitational-wave transients in data from the third observing run of Advanced LIGO and Advanced Virgo. Transients with durations of milliseconds to a few seconds in the 24--4096 Hz frequency band are targeted by the search, with no assumptions made regarding the incoming signal direction, polarization or morphology. Gravitational waves from compact binary coalescences that have been identified by other targeted analyses are detected, but no statistically significant evidence for other gravitational wave bursts is found. Sensitivities to a variety of signals are presented. These include updated upper limits on the source rate-density as a function of the characteristic frequency of the signal, which are roughly an order of magnitude better than previous upper limits. This search is sensitive to sources radiating as little as $\sim$10$^{-10} M_{\odot} c^2$ in gravitational waves at $\sim$70 Hz from a distance of 10~kpc, with 50\% detection efficiency at a false alarm rate of one per century. The sensitivity of this search to two plausible astrophysical sources is estimated: neutron star f-modes, which may be excited by pulsar glitches, as well as selected core-collapse supernova models.

preprint2021arXiv

All-sky, all-frequency directional search for persistent gravitational-waves from Advanced LIGO's and Advanced Virgo's first three observing runs

We present the first results from an all-sky all-frequency (ASAF) search for an anisotropic stochastic gravitational-wave background using the data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. Upper limit maps on broadband anisotropies of a persistent stochastic background were published for all observing runs of the LIGO-Virgo detectors. However, a broadband analysis is likely to miss narrowband signals as the signal-to-noise ratio of a narrowband signal can be significantly reduced when combined with detector output from other frequencies. Data folding and the computationally efficient analysis pipeline, {\tt PyStoch}, enable us to perform the radiometer map-making at every frequency bin. We perform the search at 3072 {\tt{HEALPix}} equal area pixels uniformly tiling the sky and in every frequency bin of width $1/32$~Hz in the range $20-1726$~Hz, except for bins that are likely to contain instrumental artefacts and hence are notched. We do not find any statistically significant evidence for the existence of narrowband gravitational-wave signals in the analyzed frequency bins. Therefore, we place $95\%$ confidence upper limits on the gravitational-wave strain for each pixel-frequency pair, the limits are in the range $(0.030 - 9.6) \times10^{-24}$. In addition, we outline a method to identify candidate pixel-frequency pairs that could be followed up by a more sensitive (and potentially computationally expensive) search, e.g., a matched-filtering-based analysis, to look for fainter nearly monochromatic coherent signals. The ASAF analysis is inherently independent of models describing any spectral or spatial distribution of power. We demonstrate that the ASAF results can be appropriately combined over frequencies and sky directions to successfully recover the broadband directional and isotropic results.

preprint2021arXiv

Higher-order valley vortices enabled by synchronized rotation in a photonic crystal

Synchronized rotation of unit cells in a periodic structure provides a novel design perspective for manipulation of band topology. We then design a two-dimensional version of higher-order topological insulators (HOTI), by such rotation in a triangular photonic lattice with $\mathcal{C}_3$ symmetry. This HOTI supports the hallmark zero-dimensional corner states and simultaneously the one-dimensional edge states. We also find that our photonic corner states carry chiral orbital angular momenta locked by valleys, whose wavefunctions are featured by the phase vortex (singularity) positioned at the maximal Wyckoff points. Moreover, when excited by a fired source with various frequencies, the valley topological states of both one-dimensional edges and zero-dimensional corners emerge simultaneously. Extendable to higher or synthetic dimensions, our work provides access to a chiral vortex platform for HOTI realisations in the THz photonic system.

preprint2021arXiv

Microscopic evidence for the intra-unit-cell electronic nematicity inside the pseudogap phase in YBa$_2$Cu$_4$O$_8$

Understanding the nature of the mysterious pseudogap phenomenon is one of the most important issues associated with cuprate high-$T_c$ superconductors. Here, we report $^{17}$O nuclear magnetic resonance (NMR) studies on two planar oxygen sites in stoichiometric cuprate YBa$_2$Cu$_4$O$_8$ to investigate the symmetry breaking inside the pseudogap phase. We observe that the Knight shifts of the two oxygen sites are identical at high temperatures but different below $T_{\rm nem} \sim$ 185 K, which is close to the pseudogap temperature $T^{\ast}$. Our result provides a microscopic evidence for intra-unit-cell electronic nematicity. The difference in quadrupole resonance frequency between the two oxygen sites is unchanged below $T_{\rm nem}$, which suggests that the observed nematicity does not directly stem from the local charge density modulation. Furthermore, a short-range charge density wave (CDW) order is observed below $T \simeq$ 150 K. The additional broadening in the $^{17}$O-NMR spectra because of this CDW order is determined to be inequivalent for the two oxygen sites, which is similar to that observed in case of nematicity. These results suggest a possible connection between nematicity, CDW order, and pseudogap.

preprint2021arXiv

Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift During the LIGO-Virgo Run O3b

We search for gravitational-wave signals associated with gamma-ray bursts detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (1 November 2019 15:00 UTC-27 March 2020 17:00 UTC).We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 gamma-ray bursts and an analysis to target binary mergers with at least one neutron star as short gamma-ray burst progenitors for 17 events. We find no significant evidence for gravitational-wave signals associated with any of these gamma-ray bursts. A weighted binomial test of the combined results finds no evidence for sub-threshold gravitational wave signals associated with this GRB ensemble either. We use several source types and signal morphologies during the searches, resulting in lower bounds on the estimated distance to each gamma-ray burst. Finally, we constrain the population of low luminosity short gamma-ray bursts using results from the first to the third observing runs of Advanced LIGO and Advanced Virgo. The resulting population is in accordance with the local binary neutron star merger rate.

preprint2021arXiv

Search for subsolar-mass binaries in the first half of Advanced LIGO and Virgo's third observing run

We report on a search for compact binary coalescences where at least one binary component has a mass between 0.2 $M_\odot$ and 1.0 $M_\odot$ in Advanced LIGO and Advanced Virgo data collected between 1 April 2019 1500 UTC and 1 October 2019 1500 UTC. We extend previous analyses in two main ways: we include data from the Virgo detector and we allow for more unequal mass systems, with mass ratio $q \geq 0.1$. We do not report any gravitational-wave candidates. The most significant trigger has a false alarm rate of 0.14 $\mathrm{yr}^{-1}$. This implies an upper limit on the merger rate of subsolar binaries in the range $[220-24200] \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$, depending on the chirp mass of the binary. We use this upper limit to derive astrophysical constraints on two phenomenological models that could produce subsolar-mass compact objects. One is an isotropic distribution of equal-mass primordial black holes. Using this model, we find that the fraction of dark matter in primordial black holes is $f_\mathrm{PBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 6\%$. The other is a dissipative dark matter model, in which fermionic dark matter can collapse and form black holes. The upper limit on the fraction of dark matter black holes depends on the minimum mass of the black holes that can be formed: the most constraining result is obtained at $M_\mathrm{min}=1 M_\odot$, where $f_\mathrm{DBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 0.003\%$. These are the tightest limits on spinning subsolar-mass binaries to date.

preprint2020arXiv

Exotic Landau Diamagnetism and Weyl-Fermions Excitations in TaAs Revealed by $^{75}$As NMR and NQR

The electronic and superconducting properties associated with the topologically non-trivial bands in Weyl semimetals have recently attracted much attention. We report the microscopic properties of the type-I Weyl semimetal TaAs measured by $^{75}$As nuclear magnetic (quadrupole) resonance under zero and elevated magnetic fields over a wide temperature range up to 500 K. The magnetic susceptibility measured by the Knight shift $K$ is found to be negative at low magnetic fields and have a strong field ($B$) dependence as ln$B$ at $T$ = 1.56 K. Such nonlinear field-dependent magnetization can be well accounted for by Landau diamagnetism arising from the 3D linearly dispersed bands, and thus is a fingerprint of topological semimetals. We further study the low-energy excitations by the spin-lattice relaxation rate 1/$T_{1}$. At zero field and 30 K $\leq T\leq$ 250 K, 1/$T_{1}T$ shows a $T^{2}$ variation due to Weyl nodes excitations. At $B \sim$ 13 T, $1/T_1T$ exhibits the same $T$-dependence but with a smaller value, scaling with $K^2\propto T^2$, which indicates that the Korringa relation also holds for a Weyl semimetal. Analysis of the Korringa ratio reveals that the energy range of the linear bands is about 250 K in TaAs.

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.

preprint2016arXiv

Cu-NMR study of oxygen disorder in ortho-II YBa2Cu3Oy

We show that 63Cu NMR spectra place strong constraints on both the nature and the concentration of oxygen defects in ortho-II YBa2Cu3Oy. Systematic deviation from ideal ortho-II order is revealed by the presence of inequivalent Cu sites in either full or empty chains. The results can be explained by two kinds of defects: oxygen clustering into additional chains, or fragments thereof, most likely present at all concentrations (6.4<y<6.6) and oxygen vacancies randomly distributed in the full chains for y<6.50 only. Furthermore, the remarkable reproducibility of the spectra in different samples with optimal ortho-II order (y=6.55) shows that chain-oxygen disorder, known to limit electronic coherence, is ineluctable because it is inherent to these compounds.

preprint2016arXiv

Electromagnetic effects on the light pseudoscalar mesons and determination of $m_u/m_d$

The MILC Collaboration has completed production running of electromagnetic effects on light mesons using asqtad improved staggered quarks. In these calculations, we use quenched photons in the noncompact formalism. We study four lattice spacings from $\approx\!0.12\:$fm to $\approx\!0.045\:$fm. To study finite-volume effects, we used six spatial lattice sizes $L/a=12$, 16, 20, 28, 40, and 48, at $a\!\approx\!0.12\:$fm. We update our preliminary values for the correction to Dashen's theorem ($ε$) and the quark-mass ratio $m_u/m_d$.

preprint2016arXiv

Orbital Order and Spin Nematicity in the Tetragonal Phase of Electron-doped Iron-Pnictides NaFe$_{1-x}$Co$_{x}$As

In copper-oxide and iron-based high temperature (high-$T_{\rm c}$) superconductors, many physical properties exhibit in-plane anisotropy, which is believed to be caused by a rotational symmetry-breaking nematic order, whose origin and its relationship to superconductivity remain elusive. In many iron-pnictides, a tetragonal-to-orthorhombic structural transition temperature $T_{\rm s}$ coincides with the magnetic transition temperature $T_{\rm N}$, making the orbital and spin degrees of freedom highly entangled. NaFeAs is a system where $T_{\rm s}$ = 54 K is well separated from $T_{\rm N}$ = 42 K, which helps simplify the experimental situation. Here we report nuclear magnetic resonance (NMR) measurements on NaFe$_{1-x}$Co$_x$As (0 $\leq x \leq$ 0.042) that revealed orbital and spin nematicity occurring at a temperature $T^{\rm *}$ far above $T_{\rm s}$ in the tetragonal phase. We show that the NMR spectra splitting and its evolution can be explained by an incommensurate orbital order that sets in below $T^{\rm *}$ and becomes commensurate below $T_{\rm s}$, which brings about the observed spin nematicity.

preprint2015arXiv

Decay constants $f_B$ and $f_{B_s}$ from HISQ simulations

We give a progress report on a project aimed at a high-precision calculation of the decay constants $f_B$ and $f_{B_s}$ from simulations with HISQ heavy and light valence and sea quarks. Calculations are carried out with several heavy valence-quark masses on ensembles with 2+1+1 flavors of HISQ sea quarks at five lattice spacings and several light sea-quark mass ratios $m_{ud}/m_s$, including approximately physical sea-quark masses. This range of parameters provides excellent control of the continuum limit and of heavy-quark discretization errors. We present a preliminary error budget with projected uncertainties of 2.2~MeV and 1.5~MeV for $f_B$ and $f_{B_s}$, respectively.

preprint2015arXiv

Electromagnetic effects on the light hadron spectrum

For some time, the MILC Collaboration has been studying electromagnetic effects on light mesons. These calculations use fully dynamical QCD, but only quenched photons, which suffices to NLO in XPT. That is, the sea quarks are electrically neutral, while the valence quarks carry charge. For the photons we use the non-compact formalism. We have new results with lattice spacing as small as 0.045 fm and a large range of volumes. We consider how well chiral perturbation theory describes these results and the implications for light quark masses

preprint2015arXiv

JUNO Conceptual Design Report

The Jiangmen Underground Neutrino Observatory (JUNO) is proposed to determine the neutrino mass hierarchy using an underground liquid scintillator detector. It is located 53 km away from both Yangjiang and Taishan Nuclear Power Plants in Guangdong, China. The experimental hall, spanning more than 50 meters, is under a granite mountain of over 700 m overburden. Within six years of running, the detection of reactor antineutrinos can resolve the neutrino mass hierarchy at a confidence level of 3-4$σ$, and determine neutrino oscillation parameters $\sin^2θ_{12}$, $Δm^2_{21}$, and $|Δm^2_{ee}|$ to an accuracy of better than 1%. The JUNO detector can be also used to study terrestrial and extra-terrestrial neutrinos and new physics beyond the Standard Model. The central detector contains 20,000 tons liquid scintillator with an acrylic sphere of 35 m in diameter. $\sim$17,000 508-mm diameter PMTs with high quantum efficiency provide $\sim$75% optical coverage. The current choice of the liquid scintillator is: linear alkyl benzene (LAB) as the solvent, plus PPO as the scintillation fluor and a wavelength-shifter (Bis-MSB). The number of detected photoelectrons per MeV is larger than 1,100 and the energy resolution is expected to be 3% at 1 MeV. The calibration system is designed to deploy multiple sources to cover the entire energy range of reactor antineutrinos, and to achieve a full-volume position coverage inside the detector. The veto system is used for muon detection, muon induced background study and reduction. It consists of a Water Cherenkov detector and a Top Tracker system. The readout system, the detector control system and the offline system insure efficient and stable data acquisition and processing.

preprint2015arXiv

New Superconductivity Dome in LaFeAsO$_{1-x}$F$_{x}$ Accompanied by Structural Transition

High temperature superconductivity is often found in the vicinity of antiferromagnetism. This is also true in LaFeAsO$_{1-x}$F$_{x}$ ($x \leq$ 0.2) and many other iron-based superconductors, which leads to proposals that superconductivity is mediated by fluctuations associated with the nearby magnetism. Here we report the discovery of a new superconductivity dome without low-energy magnetic fluctuations in LaFeAsO$_{1-x}$F$_{x}$ with 0.25$\leq x \leq$0.75, where the maximal critical temperature $T_c$ at $x_{opt}$ = 0.5$\sim$0.55 is even higher than that at $x \leq$ 0.2. By nuclear magnetic resonance and Transmission Electron Microscopy, we show that a C4 rotation symmetry-breaking structural transition takes place for $x>$ 0.5 above $T_c$. Our results point to a new paradigm of high temperature superconductivity.

preprint2014arXiv

$B\toπ\ellν$ semileptonic form factors from unquenched lattice QCD and determination of $|V_{ub}|$

We compute the $B\toπ\ellν$ semileptonic form factors and update the determination of the CKM matrix element $|V_{ub}|$. We use the MILC asqtad ensembles with $N_f=2+1$ sea quarks at four different lattice spacings in the range $a \approx 0.045$~fm to $0.12$~fm. The lattice form factors are extrapolated to the continuum limit using SU(2) staggered chiral perturbation theory in the hard pion limit, followed by an extrapolation in $q^2$ to the full kinematic range using a functional $z$-parameterization. The extrapolation is combined with the experimental measurements of the partial branching fraction to extract $|V_{ub}|$. Our preliminary result is $|V_{ub}|=(3.72\pm 0.14)\times 10^{-3}$, where the error reflects both the lattice and experimental uncertainties, which are now on par with each other.

preprint2014arXiv

Charmed and light pseudoscalar meson decay constants from four-flavor lattice QCD with physical light quarks

We compute the leptonic decay constants $f_{D^+}$, $f_{D_s}$, and $f_{K^+}$, and the quark-mass ratios $m_c/m_s$ and $m_s/m_l$ in unquenched lattice QCD using the experimentally determined value of $f_{π^+}$ for normalization. We use the MILC highly improved staggered quark (HISQ) ensembles with four dynamical quark flavors---up, down, strange, and charm---and with both physical and unphysical values of the light sea-quark masses. The use of physical pions removes the need for a chiral extrapolation, thereby eliminating a significant source of uncertainty in previous calculations. Four different lattice spacings ranging from $a\approx 0.06$ fm to $0.15$ fm are included in the analysis to control the extrapolation to the continuum limit. Our primary results are $f_{D^+} = 212.6(0.4)({}^{+1.0}_{-1.2})\ \mathrm{MeV}$, $f_{D_s} = 249.0(0.3)({}^{+1.1}_{-1.5})\ \mathrm{MeV}$, and $f_{D_s}/f_{D^+} = 1.1712(10)({}^{+29}_{-32})$, where the errors are statistical and total systematic, respectively. The errors on our results for the charm decay constants and their ratio are approximately two to four times smaller than those of the most precise previous lattice calculations. We also obtain $f_{K^+}/f_{π^+} = 1.1956(10)({}^{+26}_{-18})$, updating our previous result, and determine the quark-mass ratios $m_s/m_l = 27.35(5)({}^{+10}_{-7})$ and $m_c/m_s = 11.747(19)({}^{+59}_{-43})$. When combined with experimental measurements of the decay rates, our results lead to precise determinations of the CKM matrix elements $|V_{us}| = 0.22487(51) (29)(20)(5)$, $|V_{cd}|=0.217(1) (5)(1)$ and $|V_{cs}|= 1.010(5)(18)(6)$, where the errors are from this calculation of the decay constants, the uncertainty in the experimental decay rates, structure-dependent electromagnetic corrections, and, in the case of $|V_{us}|$, the uncertainty in $|V_{ud}|$, respectively.

preprint2014arXiv

Charmed and light pseudoscalar meson decay constants from HISQ simulations

We compute the leptonic decay constants $f_{D^+}$, $f_{D_s}$, and $f_{K^+}$, and the quark-mass ratios $m_c/m_s$ and $m_s/m_l$ in unquenched lattice QCD. We use the MILC highly improved staggered quark (HISQ) ensembles with four dynamical quark flavors. Our primary results are $f_{D^+} = 212.6(0.4)({}^{+1.0}_{-1.2})\ \mathrm{MeV}$, $f_{D_s} = 249.0(0.3)({}^{+1.1}_{-1.5})\ \mathrm{MeV}$, and $f_{D_s}/f_{D^+} = 1.1712(10)({}^{+29}_{-32})$, where the errors are statistical and total systematic, respectively. We also obtain $f_{K^+}/f_{π^+} = 1.1956(10)({}^{+26}_{-18})$, updating our previous result, and determine the quark-mass ratios $m_s/m_l = 27.35(5)({}^{+10}_{-7})$ and $m_c/m_s = 11.747(19)({}^{+59}_{-43})$. When combined with experimental measurements of the decay rates, our results lead to precise determinations of the CKM matrix elements $|V_{us}| = 0.22487(51) (29)(20)(5)$, $|V_{cd}|=0.217(1) (5)(1)$ and $|V_{cs}|= 1.010(5)(18)(6)$, where the errors are from this calculation of the decay constants, the uncertainty in the experimental decay rates, structure-dependent electromagnetic corrections, and, in the case of $|V_{us}|$, the uncertainty in $|V_{ud}|$, respectively.

preprint2014arXiv

Gradient Flow Analysis on MILC HISQ Ensembles

We report on a preliminary scale determination with gradient-flow techniques on the $N_f = 2 + 1 + 1$ HISQ ensembles generated by the MILC collaboration. The ensembles include four lattice spacings, ranging from 0.15 to 0.06 fm, and both physical and unphysical values of the quark masses. The scales $\sqrt{t_0}/a$ and $w_0/a$ are computed using Symanzik flow and the cloverleaf definition of $\langle E \rangle$ on each ensemble. Then both scales and the meson masses $aM_π$ and $aM_K$ are adjusted for mistunings in the charm mass. Using a combination of continuum chiral perturbation theory and a Taylor series ansatz in the lattice spacing, the results are simultaneously extrapolated to the continuum and interpolated to physical quark masses. Our preliminary results are $\sqrt{t_0} = 0.1422(7)$fm and $w_0 = 0.1732(10)$fm. We also find the continuum mass-dependence of $w_0$.

preprint2014arXiv

The $D_s$, $D^+$, $B_s$ and $B$ decay constants from $2+1$ flavor lattice QCD

We present a study of the $D$ and $B$ leptonic decay constants on the MILC $N_f=2+1$ asqtad gauge ensembles using asqtad-improved staggered light quarks and clover heavy quarks in the Fermilab interpretation. Our previous analysis \cite{Bazavov:2011aa} computed the decay constants at lattice spacings $a \approx 0.14, 0.11$ and $0.083$ fm. We have extended the simulations to finer $a \approx 0.058$ and $0.043$ fm lattice spacings, and have also increased statistics; this allows us to address many important sources of uncertainty. Technical advances include a two-step two-point fit procedure, better tuning of the heavy quark masses and a better determination of the axial-vector current matching. The present analysis remains blinded, so here we focus on the improvements and their predicted impact on the error budget compared to the prior analysis.

preprint2013arXiv

Charmed and strange pseudoscalar meson decay constants from HISQ simulations

We update our determinations of $f_{D^+}$, $f_{D_s}$, $f_K$, and quark mass ratios from simulations with four flavors of HISQ dynamical quarks. The availability of ensembles with light quarks near their physical mass means that we can extract physical results with only small corrections for valence- and sea-quark mass mistunings instead of a chiral extrapolation. The adjusted valence-quark masses and lattice spacings may be determined from an ensemble-by-ensemble analysis, and the results for the quark mass ratios then extrapolated to the continuum limit. Our central values of the charmed meson decay constants, however, come from an alternative analysis, which uses staggered chiral perturbation theory for the heavy-light mesons, and allows us to incorporate data at unphysical quark masses where statistical errors are often smaller. A jackknife analysis propagated through all of these steps takes account of the correlations among all the quantities used in the analysis. Systematic errors from the finite spatial size and EM effects are estimated by varying the parameters in the analysis, and systematic errors from the assumptions in the continuum extrapolation are estimated from the spread of values from different extrapolations.

preprint2013arXiv

Correlation between Internal States and Strength in Bulk Metallic Glass

The internal states or local structures of bulk metallic glass (BMGs) can be well reflected from the changes of density, structural relaxation as well as the elastic constants. With the increasing free volume (FV) content, more local atomic clusters are capable of simultaneous plastic shear at different sites in the metallic glasses, inducing large plasticity. In this work, we report a close correlation between the internal states and strength in a BMG and discover that the yield strength can be changed by varying of the casting current, revealing that the yielding strength of BMGs is not only intrinsically associated with the glass transition, but also with the internal states, such as free volume and elastic properties. Such results may have some implications for understanding the correlations between the internal states and mechanical properties of BMGs.

preprint2013arXiv

Electromagnetic contributions to pseudoscalar masses

We report on the calculation by the MILC Collaboration of the electromagnetic effects on kaon and pion masses. These masses are computed in QCD with dynamical (asqtad staggered) quarks plus quenched photons at three lattice spacings varying from 0.12 to 0.06 fm. The masses are fit to staggered chiral perturbation theory with NLO electromagnetic terms, as well as analytic terms at higher order. We extrapolate the results to physical light-quark masses and to the continuum limit. At the current stage of the analysis, most, but not all, of the systematic errors have been estimated. The main goal is the comparison of kaon electromagnetic splittings to those of the pion, i.e., an evaluation of the corrections to "Dashen's theorem." This in turn will allow us to significantly reduce the systematic errors in our determination of m_u/m_d.

preprint2013arXiv

Leptonic decay-constant ratio f_{K^+}/f_{pi^+} from lattice QCD with physical light quarks

A calculation of the ratio of leptonic decay constants f_{K^+}/f_{π^+} makes possible a precise determination of the ratio of CKM matrix elements |V_{us}|/|V_{ud}| in the Standard Model, and places a stringent constraint on the scale of new physics that would lead to deviations from unitarity in the first row of the CKM matrix. We compute f_{K^+}/f_{π^+} numerically in unquenched lattice QCD using gauge-field ensembles recently generated that include four flavors of dynamical quarks: up, down, strange, and charm. We analyze data at four lattice spacings a ~ 0.06, 0.09, 0.12, and 0.15 fm with simulated pion masses down to the physical value 135 MeV. We obtain f_{K^+}/f_{π^+} = 1.1947(26)(37), where the errors are statistical and total systematic, respectively. This is our first physics result from our N_f = 2+1+1 ensembles, and the first calculation of f_{K^+}/f_{π^+} from lattice-QCD simulations at the physical point. Our result is the most precise lattice-QCD determination of f_{K^+}/f_{π^+}, with an error comparable to the current world average. When combined with experimental measurements of the leptonic branching fractions, it leads to a precise determination of |V_{us}|/|V_{ud}| = 0.2309(9)(4) where the errors are theoretical and experimental, respectively.

preprint2013arXiv

Quantum Criticality in Electron-doped BaFe_{2-x}Ni_xAs_2

A quantum critical point (QCP) is a point in a system's phase diagram at which an order is completely suppressed at absolute zero temperature (T). The presence of a quantum critical point manifests itself in the finite-T physical properties, and often gives rise to new states of matter. Superconductivity in the cuprates and in heavy fermion materials is believed by many to be mediated by fluctuations associated with a quantum critical point. In the recently-discovered iron-pnictide high temperature superconductors, it is unknown whether a QCP exists or not in a carrier-doped system. Here we report transport and nuclear magnetic resonance (NMR) measurements on BaFe_{2-x}Ni_xAs_2 (0 =< x =< 0.17). We find two critical points at x_{c1} = 0.10 and x_{c2} = 0.14. The electrical resistivity follows ρ= ρ_0 + A*T^n, with n = 1 around x_{c1} and another minimal n = 1.1 at x_{c2}. By NMR measurements, we identity x_{c1} to be a magnetic QCP and suggest that x_{c2} is a new type of QCP associated with a nematic structural phase transition. Our results suggest that the superconductivity in carrier-doped pnictides is closely linked to the quantum criticality.

preprint2013arXiv

Symanzik flow on HISQ ensembles

We report on a scale determination with gradient-flow techniques on the $N_f = 2 + 1 + 1$ HISQ ensembles generated by the MILC collaboration. The lattice scale $w_0/a$, originally proposed by the BMW collaboration, is computed using Symanzik flow at four lattice spacings ranging from 0.15 to 0.06 fm. With a Taylor series ansatz, the results are simultaneously extrapolated to the continuum and interpolated to physical quark masses. We give a preliminary determination of the scale $w_0$ in physical units, along with associated systematic errors, and compare with results from other groups. We also present a first estimate of autocorrelation lengths as a function of flowtime for these ensembles.

preprint2013arXiv

Update on the 2+1+1 flavor QCD equation of state with HISQ

We present recent results on the QCD equation of state with 2+1+1 flavors of highly improved staggered quarks (HISQ). We focus on three sets of ensembles with temporal extent 6, 8 and 10, that reach up to temperatures of 967, 725 and 580 MeV, respectively. The strange and charm quark masses are tuned to the physical values and the light quarks mass is set to one fifth of the strange. This corresponds to a Goldstone pion of about 300 MeV.

preprint2012arXiv

Kaon semileptonic decay form factors with HISQ valence quarks

We report on the status of our kaon semileptonic form factor calculations using the highly-improved staggered quark (HISQ) formulation to simulate the valence fermions. We present results for the form factor f_+^{K π}(0) on the asqtad N_f=2+1 MILC configurations, discuss the chiral-continuum extrapolation, and give a preliminary estimate of the total error. We also present a more preliminary set of results for the same form factor but with the sea quarks also simulated with the HISQ action; these results include data at the physical light quark masses. The improvements that we expect to achieve with the use of the HISQ configurations and simulations at the physical quark masses are briefly discussed.

preprint2012arXiv

Microscopic coexistence of antiferromagnetic order and superconductivity in Ba$_{0.77}$K$_{0.23}$Fe$_{2}$As$_{2}$

We report ^{75}As nuclear magnetic resonance studies on an underdoped single-crystal Ba0.77K0.23Fe2As2 with Tc = 16.5 K. Below T_N = 46 K, the NMR peaks for H//c split and those for H //a shift to higher frequencies, which indicates that an internal magnetic field along the c axis develops below T_N. The spin-lattice relaxation rate 1/T_1 measured at the shifted peak with H//a, which experiences the internal field, shows a distinct decrease below Tc(μ_0H = 12 T) = 16 K, following a T^3 relation at low temperatures. Our results show unambiguously that antiferromagnetic order and superconductivity coexist microscopically. The unusual superconducting state with the coexisting magnetism is highlighted.

preprint2012arXiv

Status of the MILC calculation of electromagnetic contributions to pseudoscalar masses

We calculate pseudoscalar masses on gauge configurations containing the effects of 2+1 flavors of dynamical asqtad quarks and quenched electromagnetism. The lattice spacings vary from 0.12 to 0.06 fm. The masses are fit with staggered chiral perturbation theory including NLO electromagnetic terms. We attempt to extract the fit parameters for the electromagnetic contributions, while taking into account the finite volume effects, and extrapolate them to the physical limit.

preprint2011arXiv

$K$ to $ππ$ Decay amplitudes from Lattice QCD

We report a direct lattice calculation of the $K$ to $ππ$ decay matrix elements for both the $ΔI=1/2$ and 3/2 amplitudes $A_0$ and $A_2$ on 2+1 flavor, domain wall fermion, $16^3\times32\times16$ lattices. This is a complete calculation in which all contractions for the required ten, four-quark operators are evaluated, including the disconnected graphs in which no quark line connects the initial kaon and final two-pion states. These lattice operators are non-perturbatively renormalized using the Rome-Southampton method and the quadratic divergences are studied and removed. This is an important but notoriously difficult calculation, requiring high statistics on a large volume. In this paper we take a major step towards the computation of the physical $K\toππ$ amplitudes by performing a complete calculation at unphysical kinematics with pions of mass 422\,MeV at rest in the kaon rest frame. With this simplification we are able to resolve Re$(A_0)$ from zero for the first time, with a 25% statistical error and can develop and evaluate methods for computing the complete, complex amplitude $A_0$, a calculation central to understanding the $Δ=1/2$ rule and testing the standard model of CP violation in the kaon system.

preprint2011arXiv

Anytime Heuristic Search

We describe how to convert the heuristic search algorithm A* into an anytime algorithm that finds a sequence of improved solutions and eventually converges to an optimal solution. The approach we adopt uses weighted heuristic search to find an approximate solution quickly, and then continues the weighted search to find improved solutions as well as to improve a bound on the suboptimality of the current solution. When the time available to solve a search problem is limited or uncertain, this creates an anytime heuristic search algorithm that allows a flexible tradeoff between search time and solution quality. We analyze the properties of the resulting Anytime A* algorithm, and consider its performance in three domains; sliding-tile puzzles, STRIPS planning, and multiple sequence alignment. To illustrate the generality of this approach, we also describe how to transform the memory-efficient search algorithm Recursive Best-First Search (RBFS) into an anytime algorithm.

preprint2011arXiv

Continuum Limit Physics from 2+1 Flavor Domain Wall QCD

We present physical results obtained from simulations using 2+1 flavors of domain wall quarks and the Iwasaki gauge action at two values of the lattice spacing $a$, ($a^{-1}$=\,1.73\,(3)\,GeV and $a^{-1}$=\,2.28\,(3)\,GeV). On the coarser lattice, with $24^3\times 64\times 16$ points, the analysis of ref.[1] is extended to approximately twice the number of configurations. The ensembles on the finer $32^3\times 64\times 16$ lattice are new. We explain how we use lattice data obtained at several values of the lattice spacing and for a range of quark masses in combined continuum-chiral fits in order to obtain results in the continuum limit and at physical quark masses. We implement this procedure at two lattice spacings, with unitary pion masses in the approximate range 290--420\,MeV (225--420\,MeV for partially quenched pions). We use the masses of the $π$ and $K$ mesons and the $Ω$ baryon to determine the physical quark masses and the values of the lattice spacing. While our data are consistent with the predictions of NLO SU(2) chiral perturbation theory, they are also consistent with a simple analytic ansatz leading to an inherent uncertainty in how best to perform the chiral extrapolation that we are reluctant to reduce with model-dependent assumptions about higher order corrections. Our main results include $f_π=124(2)_{\rm stat}(5)_{\rm syst}$\,MeV, $f_K/f_π=1.204(7)(25)$ where $f_K$ is the kaon decay constant, $m_s^{\bar{\textrm{MS}}}(2\,\textrm{GeV})=(96.2\pm 2.7)$\,MeV and $m_{ud}^{\bar{\textrm{MS}}}(2\,\textrm{GeV})=(3.59\pm 0.21)$\,MeV\, ($m_s/m_{ud}=26.8\pm 1.4$) where $m_s$ and $m_{ud}$ are the mass of the strange-quark and the average of the up and down quark masses respectively, $[Σ^{\msbar}(2 {\rm GeV})]^{1/3} = 256(6)\; {\rm MeV}$, where $Σ$ is the chiral condensate, the Sommer scale $r_0=0.487(9)$\,fm and $r_1=0.333(9)$\,fm.

preprint2010arXiv

Continuum Limit of $B_K$ from 2+1 Flavor Domain Wall QCD

We determine the neutral kaon mixing matrix element $B_K$ in the continuum limit with 2+1 flavors of domain wall fermions, using the Iwasaki gauge action at two different lattice spacings. These lattice fermions have near exact chiral symmetry and therefore avoid artificial lattice operator mixing. We introduce a significant improvement to the conventional NPR method in which the bare matrix elements are renormalized non-perturbatively in the RI-MOM scheme and are then converted into the MSbar scheme using continuum perturbation theory. In addition to RI-MOM, we introduce and implement four non-exceptional intermediate momentum schemes that suppress infrared non-perturbative uncertainties in the renormalization procedure. We compute the conversion factors relating the matrix elements in this family of RI-SMOM schemes and MSbar at one-loop order. Comparison of the results obtained using these different intermediate schemes allows for a more reliable estimate of the unknown higher-order contributions and hence for a correspondingly more robust estimate of the systematic error. We also apply a recently proposed approach in which twisted boundary conditions are used to control the Symanzik expansion for off-shell vertex functions leading to a better control of the renormalization in the continuum limit. We control chiral extrapolation errors by considering both the NLO SU(2) chiral effective theory, and an analytic mass expansion. We obtain $B_K^{\msbar}(3 GeV) = 0.529(5)_{stat}(15)_χ(2)_{FV}(11)_{NPR}$. This corresponds to $\hat{B}_K = 0.749(7)_{stat}(21)_χ(3)_{FV}(15)_{NPR}$. Adding all sources of error in quadrature we obtain $\hat{B}_K = 0.749(27)_{combined}$, with an overall combined error of 3.6%.

preprint2010arXiv

Electromagnetic mass splittings of the low lying hadrons and quark masses from 2+1 flavor lattice QCD+QED

Results are presented for the electromagnetic mass splittings of the low lying hadrons. These are used to determine the non-degenerate light quark masses. It is found that m_u=2.24(10)(34), m_d=4.65(15)(32), and $m_s=97.6(2.9)(5.5)$ MeV (MSbar scheme, 2 GeV scale). The first error is statistical and the second systematic. We find the lowest order electromagnetic splitting (m_pi+-m_pi0)_QED=3.38(23) MeV, the splittings including next-to-leading order, (m_pi+-m_pi0)_QED=4.50(23) MeV, (m_K+-m_K0)_QED=1.87(10) MeV, and the m_u != m_d contribution to the kaon mass difference, (m_K+-m_K0)_(m_u-m_d)=-5.840(96) MeV. All errors are statistical only, and the next-to-leading order pion splitting is only approximate; it does not contain all next-to-leading order contributions. We also computed the proton-neutron mass difference, including for the first time, QED interactions in a realistic 2+1 flavor calculation. We find $(m_p-m_n)_{\rm QED}=0.383(68)$ MeV, (m_p-m_n)_(m_u-m_d)=-2.51(14) MeV, and the total m_p-m_n=-2.13(16)(70) MeV, where the first error is statistical, and the second, part of the systematic error. We use domain wall fermions and the Iwasaki gauge action (gauge coupling beta=2.13). We use two lattice sizes, 16^3 and 24^3, to address finite volume effects. Non-compact QED is treated in the quenched approximation. We present new results for the electromagnetic low energy constants in SU(3) and SU(2) partially-quenched chiral perturbation theory to the next-to-leading order, obtained from fits to our data. Detailed analysis of systematic errors in our results and methods for improving them are discussed. Finally, new analytic results for SU(2)_L x SU(2)_R-plus-kaon chiral perturbation theory, including the one-loop logs proportional to alpha_em*m, are given.

preprint2010arXiv

The $η$ and $η^\prime$ mesons from Lattice QCD

The large mass of the ninth pseudoscalar meson, the $η^\prime$, is believed to arise from the combined effects of the axial anomaly and the gauge field topology present in QCD. We report a realistic, 2+1 flavor, lattice QCD calculation of the $η$ and $η^\prime$ masses and mixing which confirms this picture. The physical eigenstates show small octet-singlet mixing with a mixing angle of $θ= -14.1(2.8)^\circ$. Extrapolation to physical light quark mass gives, with statistical errors only, $m_η=573(6)$ MeV and $m_{η^\prime}=947(142)$ MeV, consistent with the experimental values of 548 MeV and 958 MeV.

preprint2009arXiv

Chiral magnetic effect in 2+1 flavor QCD+QED

The exciting possibility of direct observation of QCD instantons in heavy-ion collisions has recently been proposed by Kharzeev. The underlying phenomenon, known as the chiral magnetic effect, may have been observed recently at RHIC, and a first principles calculation is needed to confirm and understand the results. The chiral magnetic effect is thought to be visible in the symmetric phase, at temperatures above the QCD critical temperature, and in the presence of an external magnetic field. We report on first 2+1 flavor, domain wall fermion, QCD+QED dynamical simulations above the critical temperature, in a fixed topological sector(s), which are used to study the electric charge separation produced by the effect.

Institution

Affiliation not imported yet

This author record came from a source that does not expose affiliation metadata. Once the author claims the profile or we enrich the record from another provider, this section will link to the concrete institution.

Topic footprint

Fields this researcher appears in

hep-lat hep-ph cond-mat.supr-con cond-mat.str-el cond-mat.mtrl-sci Artificial Intelligence cond-mat.mes-hall nucl-th physics.optics

Source provenance

Where this author record came from

arxivconfidence 95%

external id: arxiv:2408.09905:author:320:r-zhou

Imported May 21, 2026Synced May 21, 2026

arxivconfidence 95%

external id: arxiv:2504.00601:author:320:r-zhou

Imported May 21, 2026Synced May 21, 2026

19 works

D. Bastieri

Researcher

D. Bastieri contributes to research discovery and scholarly infrastructure.

Open to collaborate

19 works

F. Aharonian

Researcher

F. Aharonian contributes to research discovery and scholarly infrastructure.

Open to collaborate

19 works

J. Chang

Researcher

J. Chang contributes to research discovery and scholarly infrastructure.

Open to collaborate

19 works

X. J. Bi

Researcher

X. J. Bi contributes to research discovery and scholarly infrastructure.

Open to collaborate

R. Zhou

What is connected

Connect this record

See the researcher in context

Building this map preview

70 published item(s)

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

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

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

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

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

DESI Data Release 1: Stellar Catalogue

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

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

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

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

Evidence for particle acceleration approaching PeV energies in the W51 complex

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

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

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

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

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

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

Counterpart identification and classification for eRASS1 and characterisation of the AGN content

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

Pressure-enhanced $f$-electron orbital weighting in UTe2 mapped by quantum interferometry

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

The Compilation and Validation of the Spectroscopic Redshift Catalogs for the DESI-COSMOS and DESI-XMMLSS Fields

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

All-sky search for gravitational wave emission from scalar boson clouds around spinning black holes in LIGO O3 data

First joint observation by the underground gravitational-wave detector, KAGRA, with GEO600

GWTC-2.1: Deep Extended Catalog of Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Possible star-of-David pattern charge density wave with additional modulation in the kagome superconductor CsV$_3$Sb$_5$

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO--Virgo data

Search for continuous gravitational waves from 20 accreting millisecond X-ray pulsars in O3 LIGO data

Search of the Early O3 LIGO Data for Continuous Gravitational Waves from the Cassiopeia A and Vela Jr. Supernova Remnants

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

All-sky search for long-duration gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

All-sky search for short gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

All-sky, all-frequency directional search for persistent gravitational-waves from Advanced LIGO's and Advanced Virgo's first three observing runs

Higher-order valley vortices enabled by synchronized rotation in a photonic crystal

Microscopic evidence for the intra-unit-cell electronic nematicity inside the pseudogap phase in YBa$_2$Cu$_4$O$_8$

Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift During the LIGO-Virgo Run O3b

Search for subsolar-mass binaries in the first half of Advanced LIGO and Virgo's third observing run

Exotic Landau Diamagnetism and Weyl-Fermions Excitations in TaAs Revealed by $^{75}$As NMR and NQR

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

Cu-NMR study of oxygen disorder in ortho-II YBa2Cu3Oy

Electromagnetic effects on the light pseudoscalar mesons and determination of $m_u/m_d$

Orbital Order and Spin Nematicity in the Tetragonal Phase of Electron-doped Iron-Pnictides NaFe$_{1-x}$Co$_{x}$As

Decay constants $f_B$ and $f_{B_s}$ from HISQ simulations

Electromagnetic effects on the light hadron spectrum

JUNO Conceptual Design Report

New Superconductivity Dome in LaFeAsO$_{1-x}$F$_{x}$ Accompanied by Structural Transition

$B\toπ\ellν$ semileptonic form factors from unquenched lattice QCD and determination of $|V_{ub}|$

Charmed and light pseudoscalar meson decay constants from four-flavor lattice QCD with physical light quarks

Charmed and light pseudoscalar meson decay constants from HISQ simulations

Gradient Flow Analysis on MILC HISQ Ensembles

The $D_s$, $D^+$, $B_s$ and $B$ decay constants from $2+1$ flavor lattice QCD

Charmed and strange pseudoscalar meson decay constants from HISQ simulations

Correlation between Internal States and Strength in Bulk Metallic Glass

Electromagnetic contributions to pseudoscalar masses

Leptonic decay-constant ratio f_{K^+}/f_{pi^+} from lattice QCD with physical light quarks

Quantum Criticality in Electron-doped BaFe_{2-x}Ni_xAs_2

Symanzik flow on HISQ ensembles

Update on the 2+1+1 flavor QCD equation of state with HISQ

Kaon semileptonic decay form factors with HISQ valence quarks

Microscopic coexistence of antiferromagnetic order and superconductivity in Ba$_{0.77}$K$_{0.23}$Fe$_{2}$As$_{2}$

Status of the MILC calculation of electromagnetic contributions to pseudoscalar masses

$K$ to $ππ$ Decay amplitudes from Lattice QCD

Anytime Heuristic Search

Continuum Limit Physics from 2+1 Flavor Domain Wall QCD

Continuum Limit of $B_K$ from 2+1 Flavor Domain Wall QCD

Electromagnetic mass splittings of the low lying hadrons and quark masses from 2+1 flavor lattice QCD+QED

The $η$ and $η^\prime$ mesons from Lattice QCD

Chiral magnetic effect in 2+1 flavor QCD+QED