Time-reversal symmetry-breaking charge order in a kagome superconductor.

The kagome lattice1, which is the most prominent structural motif in quantum physics, benefits from inherent non-trivial geometry so that it can host diverse quantum phases, ranging from spin-liquid phases, to topological matter, to intertwined orders2-8 and, most rarely, to unconventional superconductivity6,9. Recently, charge sensitive probes have indicated that the kagome superconductors AV3Sb5 (A = K, Rb, Cs)9-11 exhibit unconventional chiral charge order12-19, which is analogous to the long-sought-after quantum order in the Haldane model20 or Varma model21. However, direct evidence for the time-reversal symmetry breaking of the charge order remains elusive. Here we use muon spin relaxation to probe the kagome charge order and superconductivity in KV3Sb5. We observe a noticeable enhancement of the internal field width sensed by the muon ensemble, which takes place just below the charge ordering temperature and persists into the superconducting state. Notably, the muon spin relaxation rate below the charge ordering temperature is substantially enhanced by applying an external magnetic field. We further show the multigap nature of superconductivity in KV3Sb5 and that the [Formula: see text] ratio (where Tc is the superconducting transition temperature and λab is the magnetic penetration depth in the kagome plane) is comparable to those of unconventional high-temperature superconductors. Our results point to time-reversal symmetry-breaking charge order intertwining with unconventional superconductivity in the correlated kagome lattice.

Ultrahigh-energy photons up to 1.4 petaelectronvolts from 12 γ-ray Galactic sources.

The extension of the cosmic-ray spectrum beyond 1 petaelectronvolt (PeV; 1015 electronvolts) indicates the existence of the so-called PeVatrons-cosmic-ray factories that accelerate particles to PeV energies. We need to locate and identify such objects to find the origin of Galactic cosmic rays1. The principal signature of both electron and proton PeVatrons is ultrahigh-energy (exceeding 100 TeV) γ radiation. Evidence of the presence of a proton PeVatron has been found in the Galactic Centre, according to the detection of a hard-spectrum radiation extending to 0.04 PeV (ref. 2). Although γ-rays with energies slightly higher than 0.1 PeV have been reported from a few objects in the Galactic plane3-6, unbiased identification and in-depth exploration of PeVatrons requires detection of γ-rays with energies well above 0.1 PeV. Here we report the detection of more than 530 photons at energies above 100 teraelectronvolts and up to 1.4 PeV from 12 ultrahigh-energy γ-ray sources with a statistical significance greater than seven standard deviations. Despite having several potential counterparts in their proximity, including pulsar wind nebulae, supernova remnants and star-forming regions, the PeVatrons responsible for the ultrahigh-energy γ-rays have not yet been firmly localized and identified (except for the Crab Nebula), leaving open the origin of these extreme accelerators.

Designing the stripe-ordered cuprate phase diagram through uniaxial-stress.

The ability to efficiently control charge and spin in the cuprate high-temperature superconductors is crucial for fundamental research and underpins technological development. Here, we explore the tunability of magnetism, superconductivity, and crystal structure in the stripe phase of the cuprate La[Formula: see text]Ba[Formula: see text]CuO[Formula: see text], with [Formula: see text] = 0.115 and 0.135, by employing temperature-dependent (down to 400 mK) muon-spin rotation and AC susceptibility, as well as X-ray scattering experiments under compressive uniaxial stress in the CuO[Formula: see text] plane. A sixfold increase of the three-dimensional (3D) superconducting critical temperature [Formula: see text] and a full recovery of the 3D phase coherence is observed in both samples with the application of extremely low uniaxial stress of [Formula: see text]0.1 GPa. This finding demonstrates the removal of the well-known 1/8-anomaly of cuprates by uniaxial stress. On the other hand, the spin-stripe order temperature as well as the magnetic fraction at 400 mK show only a modest decrease under stress. Moreover, the onset temperatures of 3D superconductivity and spin-stripe order are very similar in the large stress regime. However, strain produces an inhomogeneous suppression of the spin-stripe order at elevated temperatures. Namely, a substantial decrease of the magnetic volume fraction and a full suppression of the low-temperature tetragonal structure is found under stress, which is a necessary condition for the development of the 3D superconducting phase with optimal [Formula: see text]. Our results evidence a remarkable cooperation between the long-range static spin-stripe order and the underlying crystalline order with the three-dimensional fully coherent superconductivity. Overall, these results suggest that the stripe- and the SC order may have a common physical mechanism.

Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A.

We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at 3.67±0.05±0.15 PeV. Below the knee, the spectral index is found to be -2.7413±0.0004±0.0050, while above the knee, it is -3.128±0.005±0.027, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -0.1200±0.0003±0.0341. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components.

First Study of Antihyperon-Nucleon Scattering Λ[over ¯]p→Λ[over ¯]p and Measurement of Λp→Λp Cross Section.

Using (10.087±0.044)×10^{9} J/ψ events collected with the BESIII detector at the BEPCII storage ring, the processes Λp→Λp and Λ[over ¯]p→Λ[over ¯]p are studied, where the Λ/Λ[over ¯] baryons are produced in the process J/ψ→ΛΛ[over ¯] and the protons are the hydrogen nuclei in the cooling oil of the beam pipe. Clear signals are observed for the two reactions. The cross sections in -0.9≤cosθ_{Λ/Λ[over ¯]}≤0.9 are measured to be σ(Λp→Λp)=(12.2±1.6_{stat}±1.1_{syst}) and σ(Λ[over ¯]p→Λ[over ¯]p)=(17.5±2.1_{stat}±1.6_{syst}) mb at the Λ/Λ[over ¯] momentum of 1.074 GeV/c within a range of ±0.017 GeV/c, where the θ_{Λ/Λ[over ¯]} are the scattering angles of the Λ/Λ[over ¯] in the Λp/Λ[over ¯]p rest frames. Furthermore, the differential cross sections of the two reactions are also measured, where there is a slight tendency of forward scattering for Λp→Λp, and a strong forward peak for Λ[over ¯]p→Λ[over ¯]p. We present an approach to extract the total elastic cross sections by extrapolation. The study of Λ[over ¯]p→Λ[over ¯]p represents the first study of antihyperon-nucleon scattering, and these new measurements will serve as important inputs for the theoretical understanding of the (anti)hyperon-nucleon interaction.

Determination of Spin-Parity Quantum Numbers of X(2370) as 0

Based on (10087±44)×10^{6} J/ψ events collected with the BESIII detector, a partial wave analysis of the decay J/ψ→γK_{S}^{0}K_{S}^{0}η^{'} is performed. The mass and width of the X(2370) are measured to be 2395±11(stat)_{-94}^{+26}(syst) MeV/c^{2} and 188_{-17}^{+18}(stat)_{-33}^{+124}(syst) MeV, respectively. The corresponding product branching fraction is B[J/ψ→γX(2370)]×B[X(2370)→f_{0}(980)η^{'}]×B[f_{0}(980)→K_{S}^{0}K_{S}^{0}]=(1.31±0.22(stat)_{-0.84}^{+2.85}(syst))×10^{-5}. The statistical significance of the X(2370) is greater than 11.7σ and the spin parity is determined to be 0^{-+} for the first time. The measured mass and spin parity of the X(2370) are consistent with the predictions of the lightest pseudoscalar glueball.

First Measurement of the Decay Asymmetry in the Pure W-Boson-Exchange Decay Λ_{c}

Based on 4.4 fb^{-1} of e^{+}e^{-} annihilation data collected at the center-of-mass energies between 4.60 and 4.70 GeV with the BESIII detector at the BEPCII collider, the pure W-boson-exchange decay Λ_{c}^{+}→Ξ^{0}K^{+} is studied with a full angular analysis. The corresponding decay asymmetry is measured for the first time to be α_{Ξ^{0}K^{+}}=0.01±0.16(stat)±0.03(syst). This result reflects the noninterference effect between the S- and P-wave amplitudes. The phase shift between S- and P-wave amplitudes has two solutions, which are δ_{p}-δ_{s}=-1.55±0.25(stat)±0.05(syst) rad or 1.59±0.25(stat)±0.05(syst) rad.

Observation of Structures in the Processes e

We present measurements of the Born cross sections for the processes e^{+}e^{-}→ωχ_{c1} and ωχ_{c2} at center-of-mass energies sqrt[s] from 4.308 to 4.951 GeV. The measurements are performed with data samples corresponding to an integrated luminosity of 11.0 fb^{-1} collected with the BESIII detector operating at the Beijing Electron Positron Collider storage ring. Assuming the e^{+}e^{-}→ωχ_{c2} signals come from a single resonance, the mass and width are determined to be M=(4413.6±9.0±0.8) MeV/c^{2} and Γ=(110.5±15.0±2.9) MeV, respectively, which is consistent with the parameters of the well-established resonance ψ(4415). In addition, we also use one single resonance to describe the e^{+}e^{-}→ωχ_{c1} line shape and determine the mass and width to be M=(4544.2±18.7±1.7) MeV/c^{2} and Γ=(116.1±33.5±1.7) MeV, respectively. The structure of this line shape, observed for the first time, requires further understanding.

First Observation of a Three-Resonance Structure in e

We report the measurement of the inclusive cross sections for e^{+}e^{-}→nOCH (where nOCH denotes non-open charm hadrons) with improved precision at center-of-mass (c.m.) energies from 3.645 to 3.871 GeV. We observe three resonances: R(3760), R(3780), and R(3810) with significances of 8.1σ, 13.7σ, and 8.8σ, respectively. The R(3810) state is observed for the first time, while the R(3760) and R(3780) states are observed for the first time in the nOCH cross sections. Two sets of resonance parameters describe the energy-dependent line shape of the cross sections well. In set I [set II], the R(3810) state has mass (3805.7±1.1±2.7) [(3805.7±1.1±2.7)] MeV/c^{2}, total width (11.6±2.9±1.9) [(11.5±2.8±1.9)] MeV, and an electronic width multiplied by the nOCH decay branching fraction of (10.9±3.8±2.5) [(11.0±3.4±2.5)] eV. In addition, we measure the branching fractions B[R(3760)→nOCH]=(25.2±16.1±30.4)%[(6.4±4.8±7.7)%] and B[R(3780)→nOCH]=(12.3±6.6±8.3)%[(10.4±4.8±7.0)%] for the first time. The R(3760) state can be interpreted as an open-charm (OC) molecular state, but containing a simple four-quark state component. The R(3810) state can be interpreted as a hadrocharmonium state.

Study of the f_{0}(980) and f_{0}(500) Scalar Mesons through the Decay D_{s}

Using e^{+}e^{-} collision data corresponding to an integrated luminosity of 7.33 fb^{-1} recorded by the BESIII detector at center-of-mass energies between 4.128 and 4.226 GeV, we present an analysis of the decay D_{s}^{+}→π^{+}π^{-}e^{+}ν_{e}, where the D_{s}^{+} is produced via the process e^{+}e^{-}→D_{s}^{*±}D_{s}^{∓}. We observe the f_{0}(980) in the π^{+}π^{-} system and the branching fraction of the decay D_{s}^{+}→f_{0}(980)e^{+}ν_{e} with f_{0}(980)→π^{+}π^{-} measured to be (1.72±0.13_{stat}±0.10_{syst})×10^{-3}, where the uncertainties are statistical and systematic, respectively. The dynamics of the D_{s}^{+}→f_{0}(980)e^{+}ν_{e} decay are studied with the simple pole parametrization of the hadronic form factor and the Flatté formula describing the f_{0}(980) in the differential decay rate, and the product of the form factor f_{+}^{f_{0}}(0) and the c→s Cabibbo-Kobayashi-Maskawa matrix element |V_{cs}| is determined for the first time to be f_{+}^{f_{0}}(0)|V_{cs}|=0.504±0.017_{stat}±0.035_{syst}. Furthermore, the decay D_{s}^{+}→f_{0}(500)e^{+}ν_{e} is searched for the first time but no signal is found. The upper limit on the branching fraction of D_{s}^{+}→f_{0}(500)e^{+}ν_{e}, f_{0}(500)→π^{+}π^{-} decay is set to be 3.3×10^{-4} at 90% confidence level.

Observation of Significant Flavor-SU(3) Breaking in the Kaon Wave Function at 12

We present cross sections for the reaction e^{+}e^{-}→K_{S}^{0}K_{L}^{0} at center-of-mass energies ranging from 3.51 to 4.95 GeV using data samples collected in the BESIII experiment, corresponding to a total integrated luminosity of 26.5 fb^{-1}. The ratio of neutral-to-charged kaon form factors at large momentum transfers (12

Observation of D

We perform for the first time an amplitude analysis of the decay D^{+}→K_{S}^{0}π^{+}η and report the observation of the decay D^{+}→K_{S}^{0}a_{0}(980)^{+} using 2.93 fb^{-1} of e^{+}e^{-} collision data taken at a center-of-mass energy of 3.773 GeV with the BESIII detector. As the only W-annihilation-free decay among D to a_{0}(980) pseudoscalar, D^{+}→K_{S}^{0}a_{0}(980)^{+} is the ideal decay in extracting the contributions of the W-emission amplitudes involving a_{0}(980) and to study the final-state interactions. The absolute branching fraction of D^{+}→K_{S}^{0}π^{+}η is measured to be (1.27±0.04_{stat}±0.03_{syst})%. The branching fractions of intermediate processes D^{+}→K_{S}^{0}a_{0}(980)^{+} with a_{0}(980)^{+}→π^{+}η and D^{+}→π^{+}K[over ¯]_{0}^{*}(1430)^{0} with K[over ¯]_{0}^{*}(1430)^{0}→K_{S}^{0}η are measured to be (1.33±0.05_{stat}±0.04_{syst})% and (0.14±0.03_{stat}±0.01_{syst})%, respectively.

Determination of the Σ

Based on data samples collected with the BESIII detector at the BEPCII collider, the process e^{+}e^{-}→Σ^{+}Σ[over ¯]^{-} is studied at center-of-mass energies sqrt[s]=2.3960, 2.6454, and 2.9000 GeV. Using a fully differential angular description of the final state particles, both the relative magnitude and phase information of the Σ^{+} electromagnetic form factors in the timelike region are extracted. The relative phase between the electric and magnetic form factors is determined to be sinΔΦ=-0.67±0.29(stat)±0.18(syst) at sqrt[s]=2.3960 GeV, ΔΦ=55°±19°(stat)±14°(syst) at sqrt[s]=2.6454 GeV, and 78°±22°(stat)±9°(syst) at sqrt[s]=2.9000 GeV. For the first time, the phase of the hyperon electromagnetic form factors is explored in a wide range of four-momentum transfer. The evolution of the phase along with four-momentum transfer is an important input for understanding its asymptotic behavior and the dynamics of baryons.

Observation of D_{s}

By analyzing 7.33 fb^{-1} of e^{+}e^{-} annihilation data collected at center-of-mass energies between 4.128 and 4.226 GeV with the BESIII detector, we report the observation of the semileptonic decay D_{s}^{+}→η^{'}µ^{+}ν_{µ}, with a statistical significance larger than 10σ, and the measurements of the D_{s}^{+}→ηµ^{+}ν_{µ} and D_{s}^{+}→η^{'}µ^{+}ν_{µ} decay dynamics for the first time. The branching fractions of D_{s}^{+}→ηµ^{+}ν_{µ} and D_{s}^{+}→η^{'}µ^{+}ν_{µ} are determined to be (2.235±0.051_{stat}±0.052_{syst})% and (0.801±0.055_{stat}±0.028_{syst})%, respectively, with precision improved by factors of 6.0 and 6.6 compared to the previous best measurements. Combined with the results for the decays D_{s}^{+}→ηe^{+}ν_{e} and D_{s}^{+}→η^{'}e^{+}ν_{e}, the ratios of the decay widths are examined both inclusively and in several ℓ^{+}ν_{ℓ} four-momentum transfer ranges. No evidence for lepton flavor universality violation is found within the current statistics. The products of the hadronic form factors f_{+,0}^{η^{(')}}(0) and the c→s Cabibbo-Kobayashi-Maskawa matrix element |V_{cs}| are determined. The results based on the two-parameter series expansion are f_{+,0}^{η}(0)|V_{cs}|=0.452±0.010_{stat}±0.007_{syst} and f_{+,0}^{η^{'}}(0)|V_{cs}|=0.504±0.037_{stat}±0.012_{syst}, which help to constrain present models on f_{+,0}^{η^{(')}}(0). The forward-backward asymmetries are determined to be ⟨A_{FB}^{η}⟩=-0.059±0.031_{stat}±0.005_{syst} and ⟨A_{FB}^{η^{'}}⟩=-0.064±0.079_{stat}±0.006_{syst} for the first time, which are consistent with the theoretical calculation.

Coupled-Channel Analysis of the χ_{c1}(3872) Line Shape with BESIII Data.

We perform a study of the χ_{c1}(3872) line shape using the data samples of e^{+}e^{-}→γχ_{c1}(3872), χ_{c1}(3872)→D^{0}D[over ¯]^{0}π^{0}, and π^{+}π^{-}J/ψ collected with the BESIII detector. The effects of the coupled channels and the off-shell D^{*0} are included in the parametrization of the line shape. The line shape mass parameter is obtained to be M_{X}=(3871.63±0.13_{-0.05}^{+0.06}) MeV. Two poles are found on the first and second Riemann sheets corresponding to the D^{*0}D[over ¯]^{0} branch cut. The pole location on the first sheet is much closer to the D^{*0}D[over ¯]^{0} threshold than the other, and is determined to be 7.04±0.15_{-0.08}^{+0.07} MeV above the D^{0}D[over ¯]^{0}π^{0} threshold with an imaginary part -0.19±0.08_{-0.19}^{+0.14} MeV.

Observation of the Anomalous Shape of X(1840) in J/ψ→γ3(π

Using a sample of (10087±44)×10^{6} J/ψ events, which is about 45 times larger than that was previously analyzed, a further investigation on the J/ψ→γ3(π^{+}π^{-}) decay is performed. A significant distortion at 1.84 GeV/c^{2} in the line shape of the 3(π^{+}π^{-}) invariant mass spectrum is observed for the first time, which could be resolved by two overlapping resonant structures, X(1840) and X(1880). The new state X(1880) is observed with a statistical significance larger than 10σ. The mass and width of X(1880) are determined to be 1882.1±1.7±0.7 MeV/c^{2} and 30.7±5.5±2.4 MeV, respectively, which indicates the existence of a pp[over ¯] bound state.

Investigation of the ΔI=1/2 Rule and Test of CP Symmetry through the Measurement of Decay Asymmetry Parameters in Ξ

Using (10087±44)×10^{6} J/ψ events collected with the BESIII detector, numerous Ξ^{-} and Λ decay asymmetry parameters are simultaneously determined from the process J/ψ→Ξ^{-}Ξ[over ¯]^{+}→Λ(pπ^{-})π^{-}Λ[over ¯](n[over ¯]π^{0})π^{+} and its charge-conjugate channel. The precisions of α_{Λ0} for Λ→nπ^{0} and α[over ¯]_{Λ0} for Λ[over ¯]→n[over ¯]π^{0} compared to world averages are improved by factors of 4 and 1.7, respectively. The ratio of decay asymmetry parameters of Λ→nπ^{0} to that of Λ→pπ^{-}, ⟨α_{Λ0}⟩/⟨α_{Λ-}⟩, is determined to be 0.873±0.012_{-0.010}^{+0.011}, where the first and the second uncertainties are statistical and systematic, respectively. The ratio is smaller than unity more than 5σ, which signifies the existence of the ΔI=3/2 transition in Λ for the first time. Besides, we test for CP symmetry in Ξ^{-}→Λπ^{-} and in Λ→nπ^{0} with the best precision to date.

Experimental Limits on Solar Reflected Dark Matter with a New Approach on Accelerated-Dark-Matter-Electron Analysis in Semiconductors.

Recently a dark matter-electron (DM-electron) paradigm has drawn much attention. Models beyond the standard halo model describing DM accelerated by high energy celestial bodies are under intense examination as well. In this Letter, a velocity components analysis (VCA) method dedicated to swift analysis of accelerated DM-electron interactions via semiconductor detectors is proposed and the first HPGe detector-based accelerated DM-electron analysis is realized. Utilizing the method, the first germanium based constraint on sub-GeV solar reflected DM-electron interaction is presented with the 205.4 kg·day dataset from the CDEX-10 experiment. In the heavy mediator scenario, our result excels in the mass range of 5-15 keV/c^{2}, achieving a 3 orders of magnitude improvement comparing with previous semiconductor experiments. In the light mediator scenario, the strongest laboratory constraint for DM lighter than 0.1 MeV/c^{2} is presented. The result proves the feasibility and demonstrates the vast potential of the VCA technique in future accelerated DM-electron analyses with semiconductor detectors.

Transcriptomic and clinical heterogeneity of metastatic disease timing within metastatic castration-sensitive prostate cancer.

BACKGROUND: Metastatic castration-sensitive prostate cancer (mCSPC) is commonly classified into high- and low-volume subgroups which have demonstrated differential biology, prognosis, and response to therapy. Timing of metastasis has similarly demonstrated differences in clinical outcomes; however, less is known about any underlying biologic differences between these disease states. Herein, we aim to compare transcriptomic differences between synchronous and metachronous mCSPC and identify any differential responses to therapy. PATIENTS AND METHODS: We performed an international multi-institutional retrospective review of men with mCSPC who completed RNA expression profiling evaluation of their primary tumor. Patients were stratified according to disease timing (synchronous versus metachronous). The primary endpoint was to identify differences in transcriptomic profiles between disease timing. The median transcriptomic scores between groups were compared with the Mann-Whitney U test. Secondary analyses included determining clinical and transcriptomic variables associated with overall survival (OS) from the time of metastasis. Survival analysis was carried out with the Kaplan-Meier method and multivariable Cox regression. RESULTS: A total of 252 patients were included with a median follow-up of 39.6 months. Patients with synchronous disease experienced worse 5-year OS (39% versus 79%; P < 0.01) and demonstrated lower median androgen receptor (AR) activity (11.78 versus 12.64; P < 0.01) and hallmark androgen response (HAR; 3.15 versus 3.32; P < 0.01). Multivariable Cox regression identified only high-volume disease [hazard ratio (HR) = 4.97, 95% confidence interval (CI) 2.71-9.10; P < 0.01] and HAR score (HR = 0.51, 95% CI 0.28-0.88; P = 0.02) significantly associated with OS. Finally, patients with synchronous (HR = 0.47, 95% CI 0.30-0.72; P < 0.01) but not metachronous (HR = 1.37, 95% CI 0.50-3.92; P = 0.56) disease were found to have better OS with AR and non-AR combination therapy as compared with monotherapy (P value for interaction = 0.05). CONCLUSIONS: We have demonstrated a potential biologic difference between metastatic timing of mCSPC. Specifically, for patients with low-volume disease, those with metachronous low-volume disease have a more hormone-dependent transcriptional profile and exhibit a better prognosis than synchronous low-volume disease.

Measurement of Ultra-High-Energy Diffuse Gamma-Ray Emission of the Galactic Plane from 10 TeV to 1 PeV with LHAASO-KM2A.

The diffuse Galactic γ-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this Letter, we report the measurements of diffuse γ rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer array of the Large High Altitude Air Shower Observatory (LHAASO). Diffuse emissions from the inner (15°10 TeV). The energy spectrum in the inner Galaxy regions can be described by a power-law function with an index of -2.99±0.04, which is different from the curved spectrum as expected from hadronic interactions between locally measured cosmic rays and the line-of-sight integrated gas content. Furthermore, the measured flux is higher by a factor of ∼3 than the prediction. A similar spectrum with an index of -2.99±0.07 is found in the outer Galaxy region, and the absolute flux for 10≲E≲60 TeV is again higher than the prediction for hadronic cosmic ray interactions. The latitude distributions of the diffuse emission are consistent with the gas distribution, while the longitude distributions show clear deviation from the gas distribution. The LHAASO measurements imply that either additional emission sources exist or cosmic ray intensities have spatial variations.