We discuss prompt production of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi></math> mesons in proton-proton collisions at the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msqrt><mi>s</mi></msqrt><mo>=</mo><mrow><mn>68.5</mn><mtext> </mtext><mtext> </mtext><mi>GeV</mi></mrow></math> energy within the nonrelativistic quantum chromodynamics <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>k</mi><mi>t</mi></msub></math>-factorization approach using different unintegrated gluon distribution functions (UGDFs). We include both direct color-singlet production (<math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>g</mi><mi>g</mi><mo stretchy="false">→</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>g</mi></math>) as well as a feed-down from <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>χ</mi><mi>c</mi></msub><mo stretchy="false">→</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>γ</mi></math> and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>ψ</mi><mo>′</mo></msup><mo stretchy="false">→</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>X</mi></math> decays. The production of the decaying mesons <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>χ</mi><mi>c</mi></msub></math> or <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>ψ</mi><mo>′</mo></msup></math> is also calculated within the nonrelativistic quantum chromodynamics <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>k</mi><mi>t</mi></msub></math>-factorization approach. Differential distributions in the rapidity and transverse momentum of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi></math> are calculated and compared with experimental data of LHCb collaboration. The LHCb fixed target data allow us to test UGDFs in a new range of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>x</mi></math>.

Poland

Anna Cisek

Antoni Szczurek

Prompt inclusive production of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi></math>, <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>ψ</mi><mo>′</mo></msup></math>, and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>χ</mi><mi>c</mi></msub></math> mesons at <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msqrt><mi>s</mi></msqrt><mo>=</mo><mn>68.5</mn><mtext> </mtext><mtext> </mtext><mi>GeV</mi></math> energy within the NRQCD <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>k</mi><mi>t</mi></msub></math>-factorization approach

Based on <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mo stretchy="false">(</mo><mn>2712.4</mn><mo>±</mo><mn>14.3</mn><mo stretchy="false">)</mo><mo>×</mo><msup><mn>10</mn><mn>6</mn></msup><mtext> </mtext><mtext> </mtext><mi>ψ</mi><mo stretchy="false">(</mo><mn>3686</mn><mo stretchy="false">)</mo></math> events collected by the BESIII detector operating at the BEPCII storage ring, we perform a search for the baryon- and lepton-number violating decay <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi><mo stretchy="false">→</mo><mi>p</mi><msup><mi>e</mi><mo>−</mo></msup><mo>+</mo><mi mathvariant="normal">c</mi><mo>.</mo><mi mathvariant="normal">c</mi><mo>.</mo></math> via <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ψ</mi><mo stretchy="false">(</mo><mn>3686</mn><mo stretchy="false">)</mo><mo stretchy="false">→</mo><msup><mi>π</mi><mo>+</mo></msup><msup><mi>π</mi><mo>−</mo></msup><mi>J</mi><mo>/</mo><mi>ψ</mi></math>. No significant signal is found. An upper limit on the branching fraction of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi mathvariant="script">B</mi><mo stretchy="false">(</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mo stretchy="false">→</mo><mi>p</mi><msup><mi>e</mi><mo>−</mo></msup><mo>+</mo><mi mathvariant="normal">c</mi><mo>.</mo><mi mathvariant="normal">c</mi><mo>.</mo><mo stretchy="false">)</mo></math> is found to be <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mn>3.1</mn><mo>×</mo><msup><mn>10</mn><mrow><mo>−</mo><mn>8</mn></mrow></msup></math> at 90% confidence level.

China

M. Ablikim

Russia

M.N. Achasov

Sweden

P. Adlarson

X.C. Ai

Germany

R. Aliberti

Italy

A. Amoroso

Q. An

Y. Bai

JINR

O. Bakina

Y. Ban

H.R. Bao

V. Batozskaya

Mongolia

K. Begzsuren

N. Berger

M. Berlowski

M. Bertani

D. Bettoni

F. Bianchi

E. Bianco

A. Bortone

I. Boyko

United States

R.A. Briere

A. Brueggemann

H. Cai

M.H. Cai

X. Cai

A. Calcaterra

G.F. Cao

N. Cao

Türkiye

S.A. Cetin

X.Y. Chai

J.F. Chang

G.R. Che

Y.Z. Che

G. Chelkov

C. Chen

C.H. Chen

Chao Chen

G. Chen

H.S. Chen

H.Y. Chen

M.L. Chen

S.J. Chen

S.L. Chen

S.M. Chen

T. Chen

X.R. Chen

X.T. Chen

Y.B. Chen

Y.Q. Chen

Z.J. Chen

South Korea

S.K. Choi

X. Chu

G. Cibinetto

F. Cossio

J.J. Cui

H.L. Dai

J.P. Dai

A. Dbeyssi

R.E. de Boer

D. Dedovich

C.Q. Deng

Z.Y. Deng

A. Denig

I. Denysenko

M. Destefanis

F. De Mori

B. Ding

X.X. Ding

Y. Ding

Y.X. Ding

J. Dong

L.Y. Dong

M.Y. Dong

X. Dong

M.C. Du

S.X. Du

Y.Y. Duan

Z.H. Duan

P. Egorov

G.F. Fan

J.J. Fan

Y.H. Fan

J. Fang

S.S. Fang

W.X. Fang

Y.Q. Fang

R. Farinelli

L. Fava

F. Feldbauer

G. Felici

C.Q. Feng

J.H. Feng

Y.T. Feng

M. Fritsch

C.D. Fu

J.L. Fu

Y.W. Fu

Search for the baryon and lepton number violating decay <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi><mo stretchy="false">→</mo><mi>p</mi><msup><mi>e</mi><mo>−</mo></msup><mo>+</mo><mi mathvariant="normal">c</mi><mo>.</mo><mi mathvariant="normal">c</mi><mo>.</mo></math>

In this work, the electromagnetic structure of the hidden-charm pentaquark $$P_c(4337)$$ is investigated within the diquark–diquark–antiquark model using the QCD light-cone sum rule approach. The magnetic moments of the $$ P_c(4337) $$ state are calculated for the spin-parity assignments $$ J^P = \frac{1}{2}^- $$ and $$\frac{3}{2}^-.$$ The results are found to be $$ \mu _{P_c} = 1.76 \pm 0.44~\mu _N $$ for the $$ \frac{1}{2}^- $$ case and $$ \mu _{P_c} = -1.38 \pm 0.35~\mu _N $$ for the $$ \frac{3}{2}^- $$ scenario. These findings offer important insights into the internal quark-gluon structure and electromagnetic features of this multiquark system. Beyond their theoretical relevance, the results serve as essential benchmarks for future experimental studies aimed at determining the quantum numbers and underlying configuration of the $$ P_c(4337) .$$ Additionally, the electric quadrupole and magnetic octupole moments of the spin- $$\frac{3}{2}$$ state are extracted, indicating a non-spherical charge distribution for this exotic pentaquark.

Ulaş Özdem

Probing the electromagnetic structure of the $$P_c(4337)^+$$ pentaquark: insights from a diquark–diquark–antiquark picture for $$J^P = \frac{1}{2}^-$$ and $$\frac{3}{2}^-$$ states

Exclusive $$J/\psi $$ production is investigated in proton–proton collisions employing GPD approach with GK model. Three sets gluon density are used to calculate exclusive $$J/\psi $$ production. The survival factors and equivalent photon approximation are applied to predict the exclusive $$J/\psi $$ photoproduction in proton–proton collisions. The GPD method prediction gives a good agreement with the experimental data at LHCb. The exclusive $$J/\psi $$ production is given in proton–proton collisions at RHIC and NICA. These predictions of GPD approach can be employed to estimate the exclusive $$J/\psi $$ cross sections in future proton–proton collisions.

Ya-Ping Xie

S. Goloskokov

Exclusive $$J/\psi $$ production in proton–proton collisions adopting GPD approach

Initiated by the recent observation of a flattened line shape of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>I</mi><msup><mi>J</mi><mrow><mi>P</mi><mi>C</mi></mrow></msup><mo>=</mo><mn>0</mn><msup><mn>0</mn><mrow><mo>−</mo><mo>+</mo></mrow></msup></math> around 1.4–1.5 GeV in the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>K</mi><mi>S</mi></msub><msub><mi>K</mi><mi>S</mi></msub><msup><mi>π</mi><mn>0</mn></msup></math> invariant mass spectrum by BESIII, we make a systematic partial wave analysis of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi><mo stretchy="false">→</mo><mi>γ</mi><msub><mi>η</mi><mi>X</mi></msub><mo stretchy="false">→</mo><mi>γ</mi><mi>K</mi><mover accent="true"><mi>K</mi><mo stretchy="false">¯</mo></mover><mi>π</mi></math> based on an isobaric approach. We demonstrate that in the scenario of the first radial excitations of the isoscalar pseudoscalar from the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>K</mi><mover accent="true"><mi>K</mi><mo stretchy="false">¯</mo></mover><mi>π</mi></math> threshold to about 1.6 GeV the non-trivial <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>K</mi><mi>S</mi></msub><msub><mi>K</mi><mi>S</mi></msub><msup><mi>π</mi><mn>0</mn></msup></math> invariant mass spectrum can be explained by the coupled-channel effects with the presence of the triangle singularity mechanism. It shows that a combined fit of the Dalitz plots, three-body and two-body spectra can be achieved, which suggests that the one-state solution around 1.4–1.5 GeV proposed before still holds well. In particular, we show that the coupled-channel effects between the three important quasi-two-body decay channels, <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mrow><msup><mrow><mi>K</mi></mrow><mrow><mo>*</mo></mrow></msup><mover accent="true"><mrow><mi>K</mi></mrow><mrow><mo stretchy="false">¯</mo></mrow></mover><mo>+</mo><mi mathvariant="normal">c</mi><mo>.</mo><mi mathvariant="normal">c</mi><mo>.</mo></mrow></math>, <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mrow><mi>κ</mi><mover accent="true"><mrow><mi>K</mi></mrow><mrow><mo stretchy="false">¯</mo></mrow></mover><mo>+</mo><mi mathvariant="normal">c</mi><mo>.</mo><mi mathvariant="normal">c</mi><mo>.</mo></mrow></math> and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>a</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mn>980</mn><mo stretchy="false">)</mo><mi>π</mi></math>, can be well described by taking into account the one-loop corrections in the isobaric approach. This is because the isoscalar pseudoscalar states are coupled to the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>K</mi><mo>*</mo></msup><mover accent="true"><mi>K</mi><mo stretchy="false">¯</mo></mover><mo>+</mo><mi mathvariant="normal">c</mi><mo>.</mo><mi mathvariant="normal">c</mi><mo>.</mo></math> and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>a</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mn>980</mn><mo stretchy="false">)</mo><mi>π</mi></math> (<math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mrow><mi>κ</mi><mover accent="true"><mrow><mi>K</mi></mrow><mrow><mo stretchy="false">¯</mo></mrow></mover><mo>+</mo><mi mathvariant="normal">c</mi><mo>.</mo><mi mathvariant="normal">c</mi><mo>.</mo></mrow></math>) channels via the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>P</mi></math> and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>S</mi></math> waves, respectively. As a consequence, the coupled-channel effects can be largely absorbed into the redefinition of the tree-level effective couplings with the transition amplitudes computed to the order of one-loop corrections. Then, the coupled-channel effects can be estimated by the contributions from the one-loop rescattering amplitudes in comparison with the tree-level ones, where we find that the rescattering contributions from the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>P</mi></math> wave into the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>S</mi></math> wave, or vice versa, are apparently suppressed in the kinematic region near threshold.

Yin Cheng

Lin Qiu

Qiang Zhao

Understanding the nontrivial isoscalar pseudoscalar structures in the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mrow><msub><mrow><mi>K</mi></mrow><mrow><mi>S</mi></mrow></msub><msub><mrow><mi>K</mi></mrow><mrow><mi>S</mi></mrow></msub><msup><mrow><mi>π</mi></mrow><mrow><mn>0</mn></mrow></msup></mrow></math> spectra in the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi></math> radiative decay

A key search mode for axionlike particles (ALPs) that couple to charm quarks is <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi><mo stretchy="false">→</mo><mi>γ</mi><mi>a</mi></math>. Here we calculate the form factor that allows the rate of this process to be determined using lattice QCD for the first time. Our calculations use the relativistic highly improved staggered quark (HISQ) action for the valence charm quarks on gluon field configurations generated by the MILC collaboration that include <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>u</mi></math>, <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>d</mi></math>, <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>s</mi></math> and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>c</mi></math> HISQ quarks in the sea at four values of the lattice spacing and both unphysical and physical sea quark masses. We determine the form factor as a function of ALP mass with an uncertainty of less than 2% across our full range of ALP masses from zero up to 95% of the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi></math> mass. This represents a substantial improvement in accuracy of the theoretical picture of this decay compared to the previously used tree-level and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi mathvariant="script">O</mi><mo stretchy="false">(</mo><msub><mi>α</mi><mi>s</mi></msub><mo stretchy="false">)</mo></math> perturbation theory. We use our form factor to determine constraints on ALP masses and couplings to charm quarks and photons in several different scenarios using recent experimental data from BESIII. Our calculation paves the way for further lattice QCD input on new physics constraints from radiative decays.

United Kingdom

Brian Colquhoun

Christine T.H. Davies

G. Peter Lepage

Sophie Renner

Constraints on axionlike particles using lattice QCD calculations of the rate for <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi><mo stretchy="false">→</mo><mi>γ</mi><mi>a</mi></math>

We investigate the existence of bound states and resonances in the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ϕ</mi><mi>N</mi><mi>N</mi><mo>,</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>N</mi><mi>N</mi><mo>,</mo><msub><mi>η</mi><mi>c</mi></msub><mi>N</mi><mi>N</mi></math> systems using HAL QCD interactions for <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ϕ</mi><mi>N</mi><mo>,</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>N</mi></math>, and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>η</mi><mi>c</mi></msub><mi>N</mi></math>. We employ the Gaussian expansion method to solve the complex-scaled Schrödinger equation and find no resonances or bound states in the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>N</mi><mi>N</mi></math> and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>η</mi><mi>c</mi></msub><mi>N</mi><mi>N</mi></math> systems. We estimate the interaction between charmonium and nuclei, concluding that the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi></math> or <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>η</mi><mi>c</mi></msub></math> is likely to bind with <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mrow><mmultiscripts><mrow><mi mathvariant="normal">H</mi></mrow><mprescripts></mprescripts><none></none><mrow><mn>3</mn></mrow></mmultiscripts></mrow></math>, <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mrow><mmultiscripts><mrow><mi>He</mi></mrow><mprescripts></mprescripts><none></none><mrow><mn>3</mn></mrow></mmultiscripts></mrow></math>, <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mrow><mmultiscripts><mrow><mi>He</mi></mrow><mprescripts></mprescripts><none></none><mrow><mn>4</mn></mrow></mmultiscripts></mrow></math>, and heavier nuclei. For the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ϕ</mi><mi>N</mi><mi>N</mi></math> system, the lattice QCD <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mrow><mi>ϕ</mi><mi>N</mi><mo stretchy="false">(</mo><mrow><msub><mrow><mmultiscripts><mrow><mi>S</mi></mrow><mprescripts></mprescripts><none></none><mrow><mn>2</mn></mrow></mmultiscripts></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mrow></math> interaction is absent. We combine the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ϕ</mi><mi>p</mi></math> correlation function analysis and HAL QCD results in model A. We assume the spin-spin interactions for <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>N</mi></math> and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ϕ</mi><mi>N</mi></math> systems are inversely proportional to their masses in model B. Model A predicts a stronger <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ϕ</mi><mi>N</mi><mo stretchy="false">(</mo><msub><mrow><mmultiscripts><mrow><mi>S</mi></mrow><mprescripts></mprescripts><none></none><mrow><mn>2</mn></mrow></mmultiscripts></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msub><mo stretchy="false">)</mo></math> interaction and permits a two-body bound state, whereas model B suggests the interaction is attractive but too weak to form a bound state. Both models predict bound states for the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>I</mi><mo stretchy="false">(</mo><msup><mi>J</mi><mi>P</mi></msup><mo stretchy="false">)</mo><mo>=</mo><mn>0</mn><mo stretchy="false">(</mo><msup><mn>0</mn><mo>−</mo></msup><mo stretchy="false">)</mo></math> and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mn>0</mn><mo stretchy="false">(</mo><msup><mn>1</mn><mo>−</mo></msup><mo stretchy="false">)</mo></math> <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ϕ</mi><mi>N</mi><mi>N</mi></math> systems. In model A, these states are deeply bound with binding energies exceeding 15 MeV and remain existent when considering parameter uncertainties. In contrast, these states are very loosely bound in model B, with binding energies below 1 MeV and an existent probability of about 60% when parameter uncertainties are considered. In both models, there exist very loosely bound <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>I</mi><mo stretchy="false">(</mo><msup><mi>J</mi><mi>P</mi></msup><mo stretchy="false">)</mo><mo>=</mo><mn>0</mn><mo stretchy="false">(</mo><msup><mn>2</mn><mo>−</mo></msup><mo stretchy="false">)</mo></math> three-body states which resemble a <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ϕ</mi></math>-d atom with the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ϕ</mi></math> meson surrounding the deuteron, but their existence is sensitive to parameter uncertainties. No bound states or resonances are found in the isovector <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>I</mi><mo stretchy="false">(</mo><msup><mi>J</mi><mi>P</mi></msup><mo stretchy="false">)</mo><mo>=</mo><mn>1</mn><mo stretchy="false">(</mo><msup><mn>1</mn><mo>−</mo></msup><mo stretchy="false">)</mo></math> <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ϕ</mi><mi>N</mi><mi>N</mi></math> system.

Liang-Zhen Wen

Yao Ma

Lu Meng

Shi-Lin Zhu

<math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>ϕ</mi><mi>N</mi><mi>N</mi></math>, <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>N</mi><mi>N</mi></math>, <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>η</mi><mi>c</mi></msub><mi>N</mi><mi>N</mi></math> systems based on HAL QCD interactions

We calculate the cross sections for the processes <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>e</mi><mo>+</mo></msup><msup><mi>e</mi><mo>−</mo></msup><mo stretchy="false">→</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>η</mi></math> and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>e</mi><mo>+</mo></msup><msup><mi>e</mi><mo>−</mo></msup><mo stretchy="false">→</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><msup><mi>η</mi><mo>′</mo></msup></math> at various c.m. energies <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msqrt><mi>s</mi></msqrt></math>. We first predict these cross sections by combining nonrelativistic QCD with light cone factorization. The predicted cross sections are on the order of several femtobarns for <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>e</mi><mo>+</mo></msup><msup><mi>e</mi><mo>−</mo></msup><mo stretchy="false">→</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><msup><mi>η</mi><mo>′</mo></msup></math> and less than 1 fb for <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>e</mi><mo>+</mo></msup><msup><mi>e</mi><mo>−</mo></msup><mo stretchy="false">→</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>η</mi></math>, which are significantly smaller than the experimental measurements. It is anticipated that the cross sections are dominated by resonant contributions when <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msqrt><mi>s</mi></msqrt></math> is close to the resonance mass. In this study, we employ the vector meson dominance model to predict these resonant contributions. The effective coupling constants between the photon and the resonance, as well as between the resonance and <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>η</mi></math> are extracted from the data either provided by the latest Particle Data Group (PDG) or predicted by theoretical calculations. Taking the predictions from the factorization calculation as the continuum contribution, we predict the cross section of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>e</mi><mo>+</mo></msup><msup><mi>e</mi><mo>−</mo></msup><mo stretchy="false">→</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>η</mi></math> through a coherent sum of contributions from various resonances and the continuum. The relative phase angles between these contributions are determined through a least-<math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>χ</mi><mn>2</mn></msup></math> fit to the experimental data. We then compare our theoretical predictions with the experimental data. Additionally, we find our theoretical prediction for the cross section of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>e</mi><mo>+</mo></msup><msup><mi>e</mi><mo>−</mo></msup><mo stretchy="false">→</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>η</mi></math> is significantly larger than those for <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>e</mi><mo>+</mo></msup><msup><mi>e</mi><mo>−</mo></msup><mo stretchy="false">→</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><msup><mi>η</mi><mo>′</mo></msup></math> measured by the BESIII Collaboration.

Chang-Man Gan

Jun-Kang He

Wen-Long Sang

Min-Zhen Zhou

Revisiting the line shape of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>e</mi><mo>+</mo></msup><msup><mi>e</mi><mo>−</mo></msup><mo stretchy="false">→</mo><mi>J</mi><mo>/</mo><mi>ψ</mi><mi>η</mi></math> cross section

Within the framework of nonrelativistic-QCD factorization, we perform a comprehensive investigation of the color singlet (CS) contribution to prompt <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>J</mi><mo>/</mo><mi>ψ</mi></math> pair production at the CERN Large Hadron Collider at next-to-leading order (NLO) in <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>α</mi><mi>s</mi></msub></math>. Specifically, we compare our NLO predictions with measurements from the LHCb, CMS, and ATLAS collaborations. We find that the CS contribution itself can well describe the LHCb data in most of the experimental bins, except in the regions where the fixed-order calculation is spoiled by the emission of soft or hard-collinear gluons. In the CMS and ATLAS cases, however, the CS predictions greatly undershoot both the measured total and differential cross sections, despite sizable <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>K</mi></math> factors, of 2–3, for the total cross sections.

Zhi-Guo He

Xiao-Bo Jin

Bernd A. Kniehl

Next-to-leading-order QCD corrections to double prompt <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mrow><mrow><mi>J</mi></mrow><mo>/</mo><mi>ψ</mi></mrow></math> hadroproduction

We consider the beam-charge asymmetry of the <math altimg='si1.svg'><mi>J</mi><mo stretchy='false'>/</mo><mi>&#968;</mi></math> production cross sections in <math altimg='si2.svg'><msup><mrow><mi>&#960;</mi></mrow><mrow><mo linebreak='badbreak' linebreakstyle='after'>&#8722;</mo></mrow></msup></math>- versus <math altimg='si3.svg'><msup><mrow><mi>&#960;</mi></mrow><mrow><mo linebreak='badbreak' linebreakstyle='after'>+</mo></mrow></msup></math>-induced reactions on proton or nuclear targets. We show that the <math altimg='si1.svg'><mi>J</mi><mo stretchy='false'>/</mo><mi>&#968;</mi></math> production cross section difference between <math altimg='si2.svg'><msup><mrow><mi>&#960;</mi></mrow><mrow><mo linebreak='badbreak' linebreakstyle='after'>&#8722;</mo></mrow></msup></math> and <math altimg='si3.svg'><msup><mrow><mi>&#960;</mi></mrow><mrow><mo linebreak='badbreak' linebreakstyle='after'>+</mo></mrow></msup></math> beams impinging on a proton target has a positive sign with a magnitude proportional to the product of pion's valence quark distribution, <math altimg='si16.svg'><msub><mrow><mi>V</mi></mrow><mrow><mi>&#960;</mi></mrow></msub></math>, and proton's up and down valence quark distribution difference, <math altimg='si5.svg'><msup><mrow><mi>u</mi></mrow><mrow><mi>V</mi></mrow></msup><mo linebreak='goodbreak' linebreakstyle='after'>&#8722;</mo><msup><mrow><mi>d</mi></mrow><mrow><mi>V</mi></mrow></msup></math>. The existing <math altimg='si1.svg'><mi>J</mi><mo stretchy='false'>/</mo><mi>&#968;</mi></math> production data for <math altimg='si6.svg'><msup><mrow><mi>&#960;</mi></mrow><mrow><mo linebreak='badbreak' linebreakstyle='after'>&#8722;</mo></mrow></msup><mo linebreak='goodbreak' linebreakstyle='after'>+</mo><mi>p</mi></math> and <math altimg='si7.svg'><msup><mrow><mi>&#960;</mi></mrow><mrow><mo linebreak='badbreak' linebreakstyle='after'>+</mo></mrow></msup><mo linebreak='goodbreak' linebreakstyle='after'>+</mo><mi>p</mi></math> at 39.5 and 200 GeV/c are consistent with the expected positive beam-charge asymmetry. The magnitude of the asymmetry is compared with calculations performed within two theoretical frameworks, the Color Evaporation Model (CEM) and the Non-Relativistic QCD (NRQCD) formalism. We also examine the beam-charge dependence for pion-induced <math altimg='si1.svg'><mi>J</mi><mo stretchy='false'>/</mo><mi>&#968;</mi></math> production cross sections measured on the neutron-rich platinum target, and find good agreement between the data and theory for both the negative sign and the magnitude of the beam-charge asymmetry. The comparison between data and theoretical calculations for both proton and platinum targets suggests that the beam-charge asymmetry in pion-induced <math altimg='si1.svg'><mi>J</mi><mo stretchy='false'>/</mo><mi>&#968;</mi></math> production is a viable method of accessing the valence quark distribution of the pion.

Taiwan

Wen-Chen Chang

Japan

Marco Meyer-Conde

Jen-Chieh Peng

France

Stephane Platchkov

Takahiro Sawada

Probing pion valence quark distribution with beam-charge asymmetry of pion-induced <italic>J</italic>/ <italic>&#968;</italic> production