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.

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.

[Preparation of a dual-specific antibody targeting human CD123 and exploration of its anti-acute myeloid leukemia effects].

Objective: To construct a novel dual-specific antibody targeting human CD123 (CD123 DuAb) and study its effects in acute myeloid leukemia (AML) . Methods: Based on the variable region of the CD123 monoclonal antibody independently developed at our institution, the CD123 DuAb expression plasmid was constructed by molecular cloning and transfected into ExpiCHO-S cells to prepare the antibody protein. Through a series of in vitro experiments, its activation and proliferation effect on T cells, as well as the effect of promoting T-cell killing of AML cells, were verified. Results: ① A novel CD123 DuAb plasmid targeting CD123 was successfully constructed and expressed in the Expi-CHO eukaryotic system. ②The CD123 DuAb could bind both CD3 on T cells and CD123 on CD123(+) tumor cells. ③When T cells were co-cultured with MV4-11 cells with addition of the CD123 DuAb at a concentration of 1 nmol/L, the positive expression rates of CD69 and CD25 on T cells were 68.0% and 44.3%, respectively, which were significantly higher than those of the control group (P<0.05). ④Co-culture with CD123 DuAb at 1 nmol/L promoted T-cell proliferation, and the absolute T-cell count increased from 5×10(5)/ml to 3.2×10(6)/ml on day 9, and CFSE fluorescence intensity decreased significantly. ⑤ With the increase in CD123 DuAb concentration in the culture system, T-cell exhaustion and apoptosis increased. When the CD123 DuAb was added at a concentration of 1 nmol/L to the culture system, the proportion of CD8(+) PD-1(+) LAG-3(+) T cells was 10.90%, and the proportion of propidium iodide (PI) (-) Annexin Ⅴ(+) T cells and PI(+) Annexin Ⅴ(+) T cells was 18.27% and 11.43%, respectively, which were significantly higher than those in the control group (P<0.05). ⑥ The CD123 DuAb significantly activated T cells, and the activation intensity was positively correlated with its concentration. The expression rate of CD107a on T cells reached 16.05% with 1 nmol/L CD123 DuAb, which was significantly higher than that of the control group (P<0.05). ⑦The CD123 DuAb promoted cytokine secretion by T cells at a concentration of 1 nmol/L, and the concentration of IFN-γ and TNF-α in the supernatant of the co-culture system reached 193.8 pg/ml and 169.8 pg/ml, respectively, which was significantly higher than that of the control group (P<0.05). ⑧When CD123 DuAb was added at a concentration of 1 nmol/L to the co-culture system of T cells and CD123(+) tumor cells, the killing intensity of T cells significantly increased, and the residual rates of CD123(+) MV4-11 cells, CD123(+) Molm13 cells, and CD123(+) THP-1 cells were 7.4%, 6.7%, and 14.6% on day 3, respectively, which were significantly lower than those in the control group (P<0.05) . Conclusion: In this study, a novel CD123 DuAb was constructed and expressed. In vitro experiments verified that the DuAb binds to CD123(+) tumor cells and T cells simultaneously, promotes T-cell activation and proliferation, and facilitates their anti-leukemia effect, which provides a basis for further clinical research.

[Investigation of the immune profile of multiple myeloma patients achieving long-term survival after autologous stem cell transplantation].

Objective: To identify the characteristics of the bone marrow immune microenvironment associated with long-term survival in multiple myeloma (MM) patients. Methods: In the follow-up cohort of patients with newly diagnosed MM and who received "novel agent induction therapy and subsequent autologous stem cell transplantation and immunomodulator maintenance therapy" in the First Affiliated Hospital of Sun Yat-sen University, a cross-sectional study was carried out between August 2019 and May 2020. Using NanoString technology, the RNA expression of 770 bone marrow immune-related markers was compared between 16 patients who had progression-free survival ≥5 years and 5 patients with progressive disease. Among the 16 patients who achieved long-term survival, 9 achieved persistent minimal residual disease (MRD) negative while the other 7 had persistent positive MRD. The functional scores of each kind of immune cells were calculated based on the expression level of characteristic genes, so as to indirectly obtained the proportion of each immune cell subset. The Mann-Whitney U test and the Kruskal Wallis test were used for statistical analysis. Results: The proportion of neutrophils was significantly higher in long-surviving MM patients than in patients with progressive disease [functional scores, 13.61 (13.33, 14.25) vs. 12.93 (12.58, 13.38); Z=2.31, P=0.021]. Among long-surviving patients, those who were MRD-positive had a significantly greater number of mast cells compared with those who were MRD-negative [functional scores, 7.09 (6.49, 8.57) vs. 6.03 (5.18, 6.69); H=2.18, P=0.029]. Compared with patients with progressive disease, four genes (CTSG, IFIT2, S100B, and CHIT1) were significantly downregulated and six (C4B, TNFRSF17, CD70, IRF4, C2, and GAGE1) were upregulated in long-surviving patients. Among long-surviving patients, only gene CMA1 was significantly upgraded, 10 genes (ISG15, OAS3, MX1, IFIT2, DDX58, SIGLEC1, CXCL10, IL1RN, SERPING and TNFSF10) were significantly downregulated in the MRD-positive group compared with that in the MRD-negative group, the first 5 of which are related to the interferon response pathway. Conclusions: The increased neutrophil and mast cell numbers may be related to long-term survival in MM. Interferon signaling activation may be a key bone marrow immune profiling feature for MRD-negative, long-surviving patients with MM.

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

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.

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.

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.

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.

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.

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.

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.

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 a Vector Charmoniumlike State at 4.7 GeV/c

Using data samples with an integrated luminosity of 5.85 fb^{-1} collected at center-of-mass energies from 4.61 to 4.95 GeV with the BESIII detector operating at the BEPCII storage ring, we measure the cross section for the process e^{+}e^{-}→K^{+}K^{-}J/ψ. A new resonance with a mass of M=4708_{-15}^{+17}±21 MeV/c^{2} and a width of Γ=126_{-23}^{+27}±30 MeV is observed in the energy-dependent line shape of the e^{+}e^{-}→K^{+}K^{-}J/ψ cross section with a significance over 5σ. The K^{+}J/ψ system is also investigated to search for charged charmoniumlike states, but no significant Z_{cs}^{+} states are observed. Upper limits on the Born cross sections for e^{+}e^{-}→K^{-}Z_{cs}(3985)^{+}/K^{-}Z_{cs}(4000)^{+}+c.c. with Z_{cs}(3985)^{±}/Z_{cs}(4000)^{±}→K^{±}J/ψ are reported at 90% confidence levels. The ratio of branching fractions {[B(Z_{cs}(3985)^{+}→K^{+}J/ψ)]/B[Z_{cs}(3985)^{+}→(D[over ¯]^{0}D_{s}^{*+}+D[over ¯]^{*0}D_{s}^{+})]} is measured to be less than 0.03 at 90% confidence level.

Measurement of Energy-Dependent Pair-Production Cross Section and Electromagnetic Form Factors of a Charmed Baryon.

We study the process e^{+}e^{-}→Λ_{c}^{+}Λ[over ¯]_{c}^{-} at twelve center-of-mass energies from 4.6119 to 4.9509 GeV using data samples collected by the BESIII detector at the BEPCII collider. The Born cross sections and effective form factors (|G_{eff}|) are determined with unprecedented precision after combining the single and double-tag methods based on the decay process Λ_{c}^{+}→pK^{-}π^{+}. Flat cross sections around 4.63 GeV are obtained and no indication of the resonant structure Y(4630), as reported by Belle, is found. In addition, no oscillatory behavior is discerned in the |G_{eff}| energy dependence of Λ_{c}^{+}, in contrast to what is seen for the proton and neutron cases. Analyzing the cross section together with the polar-angle distribution of the Λ_{c}^{+} baryon at each energy point, the moduli of electric and magnetic form factors (|G_{E}| and |G_{M}|) are extracted and separated. For the first time, the energy dependence of the form factor ratio |G_{E}/G_{M}| is observed, which can be well described by an oscillatory function.

Test of CP Symmetry in Hyperon to Neutron Decays.

The quantum entangled J/ψ→Σ^{+}Σ[over ¯]^{-} pairs from (1.0087±0.0044)×10^{10} J/ψ events taken by the BESIII detector are used to study the nonleptonic two-body weak decays Σ^{+}→nπ^{+} and Σ[over ¯]^{-}→n[over ¯]π^{-}. The CP-odd weak decay parameters of the decays Σ^{+}→nπ^{+} (α_{+}) and Σ[over ¯]^{-}→n[over ¯]π^{-} (α[over ¯]_{-}) are determined to be 0.0481±0.0031_{stat}±0.0019_{syst} and -0.0565±0.0047_{stat}±0.0022_{syst}, respectively. The decay parameter α[over ¯]_{-} is measured for the first time, and the accuracy of α_{+} is improved by a factor of 4 compared to the previous results. The simultaneously determined decay parameters allow the first precision CP symmetry test for any hyperon decay with a neutron in the final state with the measurement of A_{CP}=(α_{+}+α[over ¯]_{-})/(α_{+}-α[over ¯]_{-})=-0.080±0.052_{stat}±0.028_{syst}. Assuming CP conservation, the average decay parameter is determined as ⟨α_{+}⟩=(α_{+}-α[over ¯]_{-})/2=-0.0506±0.0026_{stat}±0.0019_{syst}, while the ratios α_{+}/α_{0} and α[over ¯]_{-}/α[over ¯]_{0} are -0.0490±0.0032_{stat}±0.0021_{syst} and -0.0571±0.0053_{stat}±0.0032_{syst}, where α_{0} and α[over ¯]_{0} are the decay parameters of the decays Σ^{+}→pπ^{0} and Σ[over ¯]^{-}→p[over ¯]π^{0}, respectively.

Search for Λ[over ¯]-Λ Baryon-Number-Violating Oscillations in the Decay J/ψ→pK

We report on the first search for Λ[over ¯]-Λ oscillations in the decay J/ψ→pK^{-}Λ[over ¯]+c.c. by analyzing 1.31×10^{9} J/ψ events accumulated with the BESIII detector at the BEPCII collider. The J/ψ events are produced using e^{+}e^{-} collisions at a center of mass energy sqrt[s]=3.097 GeV. No evidence for hyperon oscillations is observed. The upper limit for the oscillation rate of Λ[over ¯] to Λ hyperons is determined to be P(Λ)=[B(J/ψ→pK^{-}Λ+c.c.)/B(J/ψ→pK^{-}Λ[over ¯]+c.c.)]<4.4×10^{-6} corresponding to an oscillation parameter δm_{ΛΛ[over ¯]} of less than 3.8×10^{-18} GeV at the 90% confidence level.

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.