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.

China

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

United States

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

Germany

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

The ALICE Collaboration at the CERN LHC has measured the inclusive production cross section of isolated photons at midrapidity as a function of the photon transverse momentum ( $$p_{\textrm{T}}^{\gamma }$$ ), in Pb–Pb collisions in different centrality intervals, and in pp collisions, at centre-of-momentum energy per nucleon pair of $$\sqrt{s_{\textrm{NN}}}~=~5.02$$ TeV. The photon transverse momentum range is between 10–14 and 40–140 GeV/ $$c$$ , depending on the collision system and on the Pb–Pb centrality class. The result extends to lower $$p_{\textrm{T}}^{\gamma }$$ than previously published results by the ATLAS and CMS experiments at the same collision energy. The covered pseudorapidity range is $$|\eta ^{\gamma } | <0.67$$ . The isolation selection is based on a charged particle isolation momentum threshold $$p_{\textrm{T}}^\mathrm{iso,~ch} = 1.5$$ GeV/ $$c$$ within a cone of radii $$R=0.2$$ and 0.4. The nuclear modification factor is calculated and found to be consistent with unity in all centrality classes, and also consistent with the HG-PYTHIA model, which describes the event selection and geometry biases that affect the centrality determination in peripheral Pb–Pb collisions. The measurement is compared to next-to-leading order perturbative QCD calculations and to the measurements of isolated photons and Z $$^{0}$$ bosons from the CMS experiment, which are all found to be in agreement.

France

S. Acharya

India

A. Agarwal

Switzerland

G. Rinella

Italy

L. Aglietta

M. Agnello

N. Agrawal

Z. Ahammed

S. Ahmad

South Korea

S. Ahn

Slovakia

I. Ahuja

A. Akindinov

V. Akishina

M. Al-Turany

D. Aleksandrov

B. Alessandro

H. Alfanda

Mexico

R. Alfaro Molina

B. Ali

A. Alici

N. Alizadehvandchali

A. Alkin

Norway

J. Alme

G. Alocco

T. Alt

A. Altamura

I. Altsybeev

J. Alvarado

C. Alvarez

M. Anaam

Romania

C. Andrei

N. Andreou

A. Andronic

E. Andronov

V. Anguelov

F. Antinori

P. Antonioli

N. Apadula

L. Aphecetche

H. Appelshäuser

C. Arata

S. Arcelli

R. Arnaldi

Brazil

J. Arneiro

I. Arsene

M. Arslandok

A. Augustinus

R. Averbeck

D. Averyanov

M. Azmi

Japan

H. Baba

A. Badalà

J. Bae

Y. Bae

Y. Baek

X. Bai

R. Bailhache

Y. Bailung

R. Bala

A. Balbino

A. Baldisseri

Poland

B. Balis

Z. Banoo

V. Barbasova

F. Barile

L. Barioglio

Greece

M. Barlou

B. Barman

Hungary

G. Barnaföldi

L. Barnby

E. Barreau

V. Barret

L. Barreto

C. Bartels

K. Barth

E. Bartsch

N. Bastid

Sweden

S. Basu

G. Batigne

D. Battistini

B. Batyunya

D. Bauri

Peru

J. Alba

Denmark

I. Bearden

C. Beattie

P. Becht

D. Behera

I. Belikov

A. Hechavarria

F. Bellini

R. Bellwied

S. Belokurova

L. Beltran

Y. Beltran

G. Bencedi

A. Bensaoula

S. Beole

Y. Berdnikov

A. Berdnikova

L. Bergmann

M. Besoiu

Measurement of the inclusive isolated-photon production cross section in pp and Pb–Pb collisions at $$\mathbf {\sqrt{\textit{s}_{NN }} = 5.02}$$ TeV

The JETSCAPE Collaboration reports a new determination of the jet transport parameter <math xmlns="http://www.w3.org/1998/Math/MathML"><mover accent="true"><mi>q</mi><mo>̂</mo></mover></math> in the quark-gluon plasma (QGP) using Bayesian inference, incorporating all available inclusive hadron and jet yield suppression data measured in heavy-ion collisions at the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC). This multi-observable analysis extends the previously published JETSCAPE Bayesian inference determination of <math xmlns="http://www.w3.org/1998/Math/MathML"><mover accent="true"><mi>q</mi><mo>̂</mo></mover></math>, which was based solely on a selection of inclusive hadron suppression data. jetscape is a modular framework incorporating detailed dynamical models of QGP formation and evolution, and jet propagation and interaction in the QGP. Virtuality-dependent partonic energy loss in the QGP is modeled as a thermalized weakly coupled plasma, with parameters determined from Bayesian calibration using soft-sector observables. This Bayesian calibration of <math xmlns="http://www.w3.org/1998/Math/MathML"><mover accent="true"><mi>q</mi><mo>̂</mo></mover></math> utilizes active learning, a machine-learning approach, for efficient exploitation of computing resources. The experimental data included in this analysis span a broad range in collision energy and centrality, and in transverse momentum. In order to explore the systematic dependence of the extracted parameter posterior distributions, several different calibrations are reported, based on combined jet and hadron data; on jet or hadron data separately; and on restricted kinematic or centrality ranges of the jet and hadron data. Tension is observed in comparison of these variations, providing new insights into the physics of jet transport in the QGP and its theoretical formulation.

R. Ehlers

Y. Chen

J. Mulligan

Y. Ji

Canada

A. Kumar

S. Mak

P.M. Jacobs

A. Majumder

A. Angerami

R. Arora

S.A. Bass

R. Datta

L. Du

H. Elfner

R.J. Fries

C. Gale

Y. He

B.V. Jacak

S. Jeon

F. Jonas

L. Kasper

M. Kordell

R. Kunnawalkam-Elayavalli

J. Latessa

Y.J. Lee

R. Lemmon

M. Luzum

A. Mankolli

C. Martin

H. Mehryar

T. Mengel

C. Nattrass

J. Norman

C. Parker

J.F. Paquet

J.H. Putschke

H. Roch

G. Roland

B. Schenke

L. Schwiebert

A. Sengupta

C. Shen

M. Singh

C. Sirimanna

D. Soeder

R.A. Soltz

Finland

I. Soudi

Bayesian inference analysis of jet quenching using inclusive jet and hadron suppression measurements