Search

The determination of Parton Distribution Functions from a finite set of data is a typical example of an inverse problem. Inverse problems are notoriously difficult to solve, in particular when a robust determination of the uncertainty in the result is needed. We present a Bayesian framework to deal with this problem and discuss first results from a closure test., Comment: 9 pages, 1 figure

The BESIII experiment at the Beijing Electron Positron collider, run as a ``hyperon factory'', recently reported the first measurement of the asymmetry parameter of the $\Lambda\to n\gamma$ decay and updated its branching fraction, which differs from the PDG average by 5.6 sigma. We highlight the impact of this new measurement on our understanding of the puzzle of weak radiative hyperon decays and pinpoint what needs to be done in the future., Comment: 5 pages, 2 figures; to appear in Science Bulletin as a Perspective

Weak radiative hyperon decays, important to test the strong interaction and relevant in searches for beyond the standard model physics, have remained puzzling both experimentally and theoretically for a long time. The recently updated branching fraction and first measurement of the asymmetry parameter of $\Lambda\to n\gamma$ by the BESIII Collaboration further exacerbate the issue, as none of the existing predictions can describe the data. We show in this letter that the covariant baryon chiral perturbation theory, with constraints from the latest measurements of hyperon non-leptonic decays, can well describe the BESIII data. The predicted branching fraction and asymmetry parameter for $\Xi^-\to\Sigma^-\gamma$ are also in agreement with the experimental data. We note that a more precise measurement of the asymmetry parameter, which is related with that of $\Sigma^+\to p\gamma$, is crucial to test Hara's theorem. We further predict the branching fraction and asymmetry parameter of $\Sigma^0\to n\gamma$, whose future measurement can serve as a highly nontrivial check on our understanding of weak radiative hyperon decay and on the covariant baryon chiral perturbation theory., Comment: 11 pages, 4 figures

In the past two decades, one of the most puzzling phenomena discovered in hadron physics is that a nominal hadronic state can actually correspond to two poles on the complex energy plane. This phenomenon was first noticed for the $\Lambda(1405)$, and then for $K_1(1270)$ and to a less extent for $D_0^*(2300)$. In this Letter, we show explicitly how the two-pole structures emerge from the underlying chiral dynamics describing the coupled-channel interactions between heavy matter particles and Nambu-Goldstone bosons. In particular, the fact that two poles appear between the two dominant coupled channels can be attributed to the particular form of the leading order chiral potentials of the Weinberg-Tomozawa form. Their lineshapes overlap with each other because the degeneracy of the two coupled channels is only broken by explicit chiral symmetry breaking of higher order. We predict that for light-quark~(pion) masses heavier than their physical values, the two-pole structures disappear, which can be easily verified by future lattice QCD simulations. Furthermore, we anticipate similar two-pole structures in other systems, such as the isopin $1/2$ $\bar{K}\Sigma_c-\pi\Xi'_c$ coupled channel, which await for experimental discoveries., Comment: 5 pages, 4 figures

The two gravitational form factors of the pion, $A^{\pi}(t)$ and $D^{\pi}(t)$, are computed as functions of the momentum transfer squared $t$ in the kinematic region $0\leq -t< 2~\text{GeV}^2$ on a lattice QCD ensemble with quark masses corresponding to a close-to-physical pion mass $m_{\pi}\approx 170~\text{MeV}$ and $N_f=2+1$ quark flavors. The flavor decomposition of these form factors into gluon, up/down light-quark, and strange quark contributions is presented in the $\overline{\text{MS}}$ scheme at energy scale $\mu=2~\text{GeV}$, with renormalization factors computed non-perturbatively via the RI-MOM scheme. Using monopole and (modified) $z$-expansion fits to the gravitational form factors, we obtain estimates for the pion momentum fraction and $D$-term that are consistent with the momentum fraction sum rule and the next-to-leading order chiral perturbation theory prediction for $D^{\pi}(0)$., Comment: 28 pages, 17 figures, 7 tables

The leading jet transport coefficients $\hat{q}$ or $\hat{e}_{2}$ encode transverse or longitudinal momentum broadening of a hard parton traversing a hot medium. Understanding their temperature dependence is key to appreciating the observed suppression of high-transverse momentum probes at RHIC or LHC collision energies. We present a first continuum extrapolated result of $\hat{q}$ computed on pure SU(3) lattices with non-trivial temperature dependence different from the weak-coupling expectation. We discuss our formalism and its challenges and status in view of obtaining $\hat{e}_{2}$ or of unquenching the calculation. We consider a hard quark subject to a single scattering on the plasma. The transport coefficients are factorized in terms of matrix elements given as integrals of non-local gauge-covariant gluon field-strength field-strength correlators. After the analytic continuation to the deep-Euclidean region, the hard scale permits to recast these as a series of local, gauge-invariant operators. The renormalized leading-twist term in this expansion is closely related to static quantities, and is computed on pure SU(3) lattices ($N_τ=4,~6,~8,$ and $10$) for a wide range of temperatures, ranging from 200MeV < T < 1GeV. Our estimate for the unquenched result in $2+1$-flavor QCD has very similar features., 6 pages, 1 figure, HardProbes2023, 26-31 March 2023, Aschaffenburg, Germany

The temperature dependence of the Landau gauge ghost propagator is investigated in pure SU(3) Yang-Mills theory with lattice QCD simulations. Its behavior around the confined-deconfined phase transition temperature, $T_c \sim 270$ MeV, is investigated. The simulations show that in the deconfined phase, the ghost propagator is enhanced for small momenta, $\lesssim 1$ GeV. Furthermore, the analysis of the spontaneous breaking of center symmetry on the ghost propagator is studied. Similarly as observed for the gluon propagator, the simulations result in a decoupling of the sectors where the phase of the Polyakov loop is either 0 or $\pm 2\pi/3$ sectors, with the latter remaining indistinguishable. The results point to the possible use of the ghost propagator as an "order parameter" for the confined-deconfined phase transition.

The structure of hadrons relevant for deep-inelastic scattering are completely characterised by the Compton amplitude. A direct calculation of the Compton amplitude in a lattice QCD setup provides a way to accessing the structure functions, circumventing the operator mixing and renormalisation issues of the standard operator product expansion approach. In this contribution, we focus on the QCDSF/UKQCD Collaboration's advances in calculating the forward Compton amplitude via an implementation of the second-order Feynman-Hellmann theorem. We highlight our progress in investigating the moments of nucleon structure functions., Presented at DIS2023: XXX International Workshop on Deep-Inelastic Scattering and Related Subjects, Michigan State University, USA, 27-31 March 2023. arXiv admin note: substantial text overlap with arXiv:2212.09197, arXiv:2110.01310

We investigate the finite-temperature quantum chromodynamics (QCD) on a rotating lattice with $N_f=2+1$ staggered fermions and the projective plane boundary condition. We observe a negative rotational rigidity (defined in the main text) and a negative quark spin susceptibility associated with the chiral vortical effect. In contrast to most of the effective model predictions, we find that the chiral condensate decreases and the Polyakov loop increases with imaginary rotation, implying a rotational catalysis of chiral symmetry breaking and confinement by real rotation. We determine the phase boundaries for both chiral and confinement-deconfinement phase transitions on the $\Omega_I$-$T$ plane, where $\Omega_I$ is the imaginary angular velocity., Comment: 17 pages, 7 figures

We revisit the simplest Bars-Yankielowicz (BY) model (the $\psi\eta$ model), starting from a model with an additional Dirac pair of fermions in the fundamental representation, together with a complex color-singlet scalar $\phi$ coupled to them through a Yukawa interaction. This model possesses a color-flavor-locked 1-form ${\mathbbm Z}_N$ symmetry, due to the intersection of the color $SU(N)$ and two nonanomalous $U(1)$ groups. In the bulk, the model reduces to the $\psi\eta$ model studied earlier when $\phi$ acquires a nonzero vacuum expectation value and the extra fermions pair up, get massive and decouple (thus we will call our extended theory as the ``X-ray model"), while it provides a regularization of the $\mathbbm Z_2$ fluxes needed to study the $\mathbbm Z_2$ anomaly. The anomalies involving the 1-form ${\mathbbm Z}_N$ symmetry reduce, for $N$ even, exactly to the mixed ${\mathbbm Z}_2$ anomaly found earlier in the $\psi\eta$ model. The present work is a first significant step to clarify the meaning of the mixed ${\mathbbm Z}_2-[{\mathbbm Z}_N^{(1)}]^2$ anomaly found in the $\psi\eta$ and in other BY and Georgi-Glashow type $SU(N)$ models with even $N$., Comment: 26 pages

Puzzles in the determination of the hadronic-vacuum-polarization contribution currently impede a conclusive interpretation of the precision measurement of the anomalous magnetic moment of the muon at the Fermilab experiment. One such puzzle concerns tensions between evaluations in lattice QCD and using $e^+e^-\to\text{hadrons}$ cross-section data. In lattice QCD, the dominant isospin-symmetric part and isospin-breaking (IB) corrections are calculated separately, with very different systematic effects. Identifying these two pieces in a data-driven approach provides an opportunity to compare them individually and trace back the source of the discrepancy. Here, we estimate the IB component of the lattice-QCD calculations from phenomenology, based on a comprehensive study of exclusive contributions that can be enhanced via infrared singularities, threshold effects, or hadronic resonances, including, for the first time, in the $e^+e^-\to3\pi$ channel. We observe sizable cancellations among different channels, with a sum that even suggests a slightly larger result for the QED correction than obtained in lattice QCD. We conclude that the tensions between lattice QCD and $e^+e^-$ data therefore cannot be explained by the IB contributions in the lattice-QCD calculations., Comment: 7 pages, 2 figures

We present the results of a survey meant to assess the opinion of the high-energy physics theory (HET) community on the January 7th postdoc acceptance deadline - specifically, whether there is a preference to shift the deadline to later in January or February. This survey, which served for information-gathering purpose only, is part of a community conversation on the optimal timing of an acceptance deadline and whether the community would be better served by a later date. In addition, we present an analysis of data from the postdoc Rumor Mill, which gives a picture of the current hiring landscape in the field. We emphasize the importance of preserving a universal deadline, and the current results of our survey show broad support for a shift to a later date. A link to the survey, frequently asked questions, a running list of supporters, and next steps can be found on our companion web page., 13 pages, 8 figures and 2 tables webpage: https://het-postdoc-accord.github.io/information/

We investigate the mechanics of the proton by examining the flavor-decomposed proton cosmological constants and generalized vector form factors. The interplay of up, down, and strange quarks within the proton is explored, shedding light on its internal structure. The contributions of strange quarks play a crucial role in the $D$-term and cosmological constants. We find that the flavor blindness of the isovector $D$-term form factor is only valid in flavor SU(3) symmetry., 6 pages, 2 figures

We apply dilaton chiral perturbation theory (dChPT) at next-to-leading order to lattice data from the LatKMI collaboration for the eight-flavor SU(3) gauge theory. In previous work, we found that leading-order dChPT does not account for these data, but that a model extension of leading-order dChPT with a varying mass anomalous dimension describes these data well. Here we calculate the next-to-leading order corrections for the pion mass and decay constant. We focus on these quantities, as data for the dilaton mass are of poorer quality. The application of next-to-leading order dChPT is difficult because of the large number of new low-energy constants, and the the results of our fits turn out to be inconclusive. They suggest -- yet cannot firmly establish -- that the LatKMI mass range might be outside the scope of dChPT., RevTex, 28 pages

Assuming that the $P_c(4312)^+$ is a $I(J^P)=\frac12(\frac32^-)$ compact pentaquark, we study the mass spectrum of its S-wave hidden-charm partner states in a color-magnetic interaction model. Combining the information from their decays obtained in a simple rearrangement scheme, one finds that the quantum numbers of $P_c(4457)^+$, $ P_c(4440)^+$, and $P_c(4337)^+$ can be assigned to be $I(J^P)=\frac12(\frac32^-)$, $\frac12(\frac12^-)$, and $\frac12(\frac12^-)$, respectively, while both $P_{cs}(4338)^0$ and $P_{cs}(4459)^0$ can be interpreted as $I(J^P)=0(\frac12^-)$ $udsc\bar{c}$ compact states. Based on the numerical results, we also find narrow pentaquarks in $ssnc\bar{c}$ ($n=u,d$) and $sssc\bar{c}$ systems. The decay properties of the studied pentaquarks and the searching channels for them can be tested in future experiments., 16 pages, 14 tables, 4 figures

Pulsar Timing Arrays (PTAs) have reported evidence for a stochastic gravitational wave (GW) background at nHz frequencies, possibly originating in the early Universe. We show that the spectral shape of the low-frequency (causality) tail of GW signals sourced at temperatures around $T\gtrsim 1$ GeV is distinctively affected by confinement of strong interactions (QCD), due to the corresponding sharp decrease in the number of relativistic species. Bayesian analyses in the NANOGrav 15 years and the previous International PTA datasets reveal a significant improvement in the fit with respect to cubic power-law spectra, previously employed for the causality tail. This suggests that the inclusion of Standard Model effects on GWs can have a potentially decisive impact on model selection., 5+11 pages, 7 figures. v2: Analysis with NANOGrav 15 yrs added, minor corresponding changes in abstract and text, references added

We present a detailed exploration of certain outstanding features of the transversely-projected three-gluon vertex, using the corresponding Schwinger-Dyson equation in conjunction with key results obtained from quenched lattice simulations. The main goal of this study is the scrutiny of the approximate property denominated ``planar degeneracy'', unveiled when the Bose symmetry of the vertex is properly exploited. The planar degeneracy leads to a particularly simple parametrization of the vertex, reducing its kinematic dependence to essentially a single variable. Our analysis, carried out in the absence of dynamical quarks, reveals that the planar degeneracy is particularly accurate for the description of the form factor associated with the classical tensor, for a wide array of arbitrary kinematic configurations. Instead, the remaining three form factors display considerable violations of this property. In addition, and in close connection with the previous point, we demonstrate the numerical dominance of the classical form factor over all others, except in the vicinity of the soft-gluon kinematics. The final upshot of these considerations is the emergence of a very compact description for the three-gluon vertex in general kinematics, which may simplify significantly nonperturbative applications involving this vertex., 39 pages, 10 figures

We present the first data-driven result for $a_\mu^{\rm win,lqc}$, the isospin-limit light-quark connected component of the intermediate-window Hadronic-Vacuum-Polarization contribution to the muon anomalous magnetic moment. Our result, $(198.8\pm 1.1)\times 10^{-10}$, is in significant tension with eight recent mutually compatible high-precision lattice-QCD determinations, and provides enhanced evidence for a puzzling discrepancy between lattice and data-driven determinations of the intermediate window quantity, one driven largely by a difference in the light-quark connected component., Comment: 6 pages, 1 figure

Using the PDG 22 compilation of the $e^+e^-\to$ Hadrons $\oplus$ the recent CMD3 data for the pion form factor and the value of gluon condensate $$ from heavy quarkonia, we extract the value of the four-quark condensate : $\rho\alpha_s^2= (5.98\pm 0.64)\times 10^{-4}$ GeV$^6$ and the dimension eight condensate: $d_8= (4.3\pm 3.0)\times 10^{-2}$ GeV$^8$from the ratio ${\cal R}_{10}$ of Laplace sum rules to order $\alpha_s^4$. We show the inconsistency in using at the same time the standard SVZ value of the gluon and the vacuum saturation of the four-quark condensates. Using the previous values of the four-quark and $d_8$ condensates, we re-extract $$ from ${\cal R}_{10}$ to be: $(6.12\pm 0.61)\times 10^{-2}$ Gev$^4$ in perfect agreement with the one from heavy quarkonia. We also use the lowest $\tau$-like decay moment ${\cal R}_\tau^{ee}$ to extract the value of the QCD coupling $\alpha_s(M^2_\tau)$=0.3385(145)[resp. 0.3262(86)] (mean of fixed order (FO) and Contour Improved (CI) PT series) to order $\alpha_s^4$ [resp. $\alpha_s^5$] and the standard OPE. The corresponding value of the sum of the non-perturbative contribution is: $\delta_{NP}(M_\tau)=(3.74\pm 0.40)\times 10^{-2}$. Reciprocally, using $\alpha_s(M_\tau)$, $$ and $d_8$ as inputs, we test the stability of the value of the four-quark condensate obtained from the lowest $\tau$-like moment. We complete our analysis by updating our previous determinations of the lowest order hadronic vacuum polarization contributions to the lepton anomalies and to $\alpha(M^2_Z)$. We obtain in Table 2 : $a_\mu\vert^{hvp}_{l.o}= (7036.5\pm 38.9)\times10^{-11}, a_\tau\vert^{hvp}_{l.o}= (3494.8\pm 24.7)\times10^{-9} $. and $\alpha(M^2_Z)=(2766.3\pm 4.5)\times 10^{-5}$. This new value of $a_\mu$ leads to: $\Delta a_\mu\equiv a_\mu^{exp}-a_\mu^{th} = (142\pm 42_{th}\pm 41_{exp})\times 10^{-11}$., Comment: Latex files, 37 pages, 19 figures, 3 Tables. Inclusion of PT a_s^4 term for the determination of the QCD condensates. Some references added. Complementary discussions on the determination of \alpha_s from $\tau$-like decay and on the muon anomaly. Revised version submitted for publication to the journal

We propose a new method to calculate parton distribution functions (PDFs) from correlations of boosted quarks and gluons in the Coulomb gauge. Compared to the widely used quasi-PDFs defined from gauge-invariant Wilson-line operators, such correlations offer advantages including accessibility to larger off-axis momenta, absence of linear power divergence, and enhanced long-range precision. We verify the validity of this method at next-to-leading order in perturbation theory and use it to calculate the pion valence quark PDF on a lattice with spacing $a=0.06$ fm and valence pion mass $m_\pi=300$ MeV. Our result agrees with that from the gauge-invariant quasi-PDF at similar precision, achieved with only half the computational cost through a large off-axis momentum $|\vec{p}|\sim2.2$ GeV. This opens the door to a more efficient way to calculate parton physics on the lattice., Comment: 5 pages, 4 figures, and the appendix

By extracting the transition amplitudes, we give the first lattice QCD prediction of the two-body decay partial widths of the $1^{-+}$ charmoniumlike hybrid $\eta_{c1}$. Given the calculated mass value $m_{\eta_{c1}}=4.329(36)$ GeV, the $\eta_{c1}$ decay is dominated by the open charm modes $D_1\bar{D}$, $D^*\bar{D}$ and $D^*\bar{D}^*$ with partial widths of $258(133)$ MeV, $88(18)$ MeV and $150(118)$ MeV, respectively. The coupling of $\eta_{c1}$ to $\chi_{c1}$ plus a flavor singlet pseudoscalar is not small, but $\chi_{c1}\eta$ decay is suppressed by the small $\eta-\eta'$ mixing angle. The partial width of $\eta_{c1}\to \eta_c\eta'$ is estimated to be around 1 MeV. We suggest experiments to search for $\eta_{c1}$ in the $P$-wave $D^*\bar{D}$ and $D^*\bar{D}^*$ systems. Especially, the polarization of $D^*\bar{D}^*$ can be used to distinguish the $1^{-+}$ product (total spin $S=1$) from $1^{--}$ products ($S=0$)., Comment: 16 pages, 8 figures

The free scalar field is only conformally invariant when non-minimally coupled to gravity. In flat space this amounts to amending, or improving, the energy momentum tensor. A no-go theorem prohibits the improvement for Goldstone bosons, originating from global internal spontaneous symmetry breaking. It is shown that the no-go theorem can be circumvented in the presence of a dilaton. The latter is a (pseudo) Goldstone boson originating from spontaneous conformal symmetry breaking in a theory with an infrared fixed point. Specifically, the tracelessness of the energy momentum tensor is demonstrated for a generic $d$-dimensional curved space. Additionally, the Goldstone gravitational form factors are shown to obey conformality constraints in the soft limit. The crucial point is that the remainder term of the soft theorem is non-zero due to the presence of the dilaton pole. For Goldstone systems with a trivial infrared fixed point the leading order analysis of this paper ought to be sufficient. Loop effects govern the improvement term outside the fixed point and are scheme-dependent as briefly discussed towards the end of the paper., 14 pages + refs, 1 figure

Non-perturbative results in QCD-like theories can be derived employing positivity of the Euclidean path integral measure, as pioneered by Weingarten, Vafa, Witten. We show that positivity of the measure can be generalized to parity-invariant theories with Yukawa couplings to fundamental scalars, provided the fermions are Dirac and carry a real representation of the gauge group. This result allows us to extend the proof of parity and vector-like flavor symmetries conservation to such theories, as well as to derive exact inequalities among hadrons' masses., 20 pages plus appendix

We determine masses of light quarks ($m_u$,$m_d$,$m_s$) using Borel-Laplace sum rules and renormalization group summed perturbation theory (RGSPT) from the divergence of the axial vector current. The RGSPT significantly reduces the scale dependence of the finite order perturbative series for the renormalization group (RG) invariant quantities such as spectral function, the second derivative of the polarization function of the pseudoscalar current correlator, and its Borel transformation. In addition, the convergence of the spectral function is significantly improved by summing all running logarithms and kinematical $\pi^2$-terms. Using RGSPT, we find $m_s(2 GeV)=104.34_{-4.21}^{+4.23}\hs\MeV$, and $m_d(2 GeV)=4.21_{-0.45}^{+0.48}\hs\MeV$ leading to $m_u(2 GeV)=2.00_{-0.40}^{+0.33}\hs\MeV$., Comment: 17 pages and 8 figures. Figures updated, new references added, and minor text corrections

Analytic continuation of the perturbative series from spacelike to timelike regions is performed using renormalization group summed perturbation theory (RGSPT). This method provides an all-order summation of kinematic ``$\pi^2$-terms'' accessible from a given order of a perturbative series. The impact of the summation of these terms is studied for Higgs boson decay and electromagnetic R-ratio in the perturbative QCD. Results obtained using RGSPT have improved convergence behavior in addition to significantly reduced renormalization scale dependence compared to fixed-order perturbation theory (FOPT). The higher-order behavior using the Pad\'e approximant is also studied for processes considered., Comment: 25 pages and 10 figures. References added, figures are updated, and minor text corrections

The shear viscosity of strongly interacting dense heavy-quarks plasma is evaluated analytically using a methodology valid for strongly-correlated nonequilibrium dense matter. The shear viscosity turns out to be directly proportional to the zero-frequency limit of the spectral function. By evaluating the latter using lattice QCD data from the literature, the vanishing of interaction-dominated viscosity contribution of heavy quarks to the QGP plasma is demonstrated., arXiv admin note: substantial text overlap with arXiv:2306.05771

We study a bottom-up, holographic description of a field theory yielding the spontaneous breaking of an approximate SO(5) global symmetry to its SO(4) subgroup. The weakly-coupled, six-dimensional gravity dual has regular geometry. One of the dimensions is compactified on a circle that shrinks smoothly to zero size at a finite value of the holographic direction, hence introducing a physical scale in a way that mimics the effect of confinement in the dual four-dimensional field theory. We study the spectrum of small fluctuations of the bulk fields carrying SO(5) quantum numbers, which can be interpreted as spin-0 and spin-1 bound states in the dual field theory. This work supplements an earlier publication, focused only on the SO(5) singlet states. We explore the parameter space of the theory, paying particular attention to composite states that have the right quantum numbers to be identified as pseudo-Nambu-Goldstone Bosons (PNGBs). We find that in this model the PNGBs are generally heavy, with masses of the same order as other bound states, indicating the presence of a sizeable amount of explicit symmetry breaking in the field theory side. But we also find a qualitatively new, unexpected result. When the dimension of the field-theory operator inducing SO(5) breaking is close to half of the space-time dimensionality, there exists a region of parameter space in which the PNGBs and the lightest scalar are both parametrically light in comparison to all other bound states of the field theory. Although this region is known to yield metastable classical backgrounds, this finding might be relevant to model building in the composite Higgs context., 26 pages, 6 figures; Version accepted for publication

We delve deeper into the potential composition of dark matter as stable scalar glueballs from a confining dark $SU(N)$ gauge theory, focusing on $N=\{3,4,5\}$. To predict the relic abundance of glueballs for the various gauge groups and scenarios of thermalization of the dark gluon gas, we employ a thermal effective theory that accounts for the strong-coupling dynamics in agreement with lattice simulations. We compare our methodology with previous works and find that our approach is more comprehensive and reliable. The results are encouraging and show that glueballs can account for the totality of dark matter in many unconstrained scenarios with a phase transition scale $20~{\rm MeV}\lesssim\Lambda\lesssim10^{10}~{\rm GeV}$, thus opening the possibility of exciting future studies., Comment: 12 pages, 6 figures; update and generalization of arXiv:2207.13716

This is the second of two companion papers in which we continue to develop the construction of the doubly heavy pentaquark systems using the gauge/string duality. In this paper, we propose a stringy description of the $Q\bar Qqqq$ system in the case of two light flavors. Our goal is to explore the lower-lying Born-Oppenheimer potentials as a function of the separation distance between the heavy quark-antiquark pair. The analysis shows that the ground state Born-Oppenheimer potential is described in terms of both hadro-quarkonia and hadronic molecules. Meanwhile a standard pentaquark configuration, which describes a genuine five-quark interaction, makes the dominant contribution to a higher lying potential. This configuration has an antiquark-diquark-diquark structure $\bar Q[qq][Qq]$ for separations larger than $0.1\,\text{fm}$. The latter enables us to establish a relation among the masses of hadrons in the heavy quark limit. To describe the structure of the potentials more clear, we define some critical separations that are related to the processes of string reconnection, breaking and junction annihilation. Additionally, we consider the generalized baryon vertices, where more than three strings can meet, and explore their implications for the pentaquark systems., 29 pages, many figures. arXiv admin note: text overlap with arXiv:2211.12305; v2: minor corrections

We present a detailed study of the finite momentum dynamics of the $O(4)$ critical point of QCD, which lies in the dynamic universality class of Model G. The critical scaling of the model is analyzed in multiple dynamical channels. For instance, the finite momentum analysis allows us to precisely extract the pion dispersion curve below the critical point. The pion velocity is in striking agreement with the predictions relation and static universality. The pion damping rate and velocity are both consistent with the dynamical critical exponent $\zeta = 3/2$ of Model G. Similarly, although the critical amplitude for the diffusion coefficient of the conserved $O(4)$ charges is small, it is clearly visible both in the restored phase and with finite explicit symmetry breaking, and its dynamical scaling is again consistent with $\zeta=3/2$. We determine a new set of universal dynamical critical amplitude ratios relating the diffusion coefficient to a suitably defined order parameter relaxation time. We also show that in a finite volume simulation, the chiral condensate diffuses on the coset manifold in a manner consistent with dynamical scaling, and with a diffusion coefficient that is determined by the transport coefficients of hydrodynamic pions. Finally, the amplitude ratios (together with other non-universal amplitudes also reported here) compile all relevant information for further studies of Model G both in and out of equilibrium., Comment: 34 pages, 13 figures

The possibility that gauge theories with chiral symmetry breaking exhibit an infrared fixed point is explored. First, explicit symmetry breaking by a quark mass term and renormalisation group arguments, with reasonable assumptions, imply that the quark mass anomalous dimension at the infrared fixed point is $\gamma_* = 1$. Second, the value $\gamma_* = 1$ is compatible with the long-distance interpretation of the adjoint scalar correlation function and leading order chiral perturbation theory. Third, the value $\gamma_* =1$ gives rise to a consistent matching of the trace anomaly and the Feynman-Hellmann theorem. This suggests the following picture for the conformal window. Close to its upper boundary $\gamma_*$ is zero and grows as the number of fermion is reduced until $\gamma_*=1$ is reached where chiral symmetry breaking sets in, marking its lower boundary. Below, the strongly coupled gauge theory is infrared dual to the free theory of pions. The dilaton sector of the scenario will be addressed in a companion paper., Comment: 12 pages+refs, 3 figures

We compute the structure-dependent axial and vector form factors for the radiative leptonic decays $D_s\to \ell\nu_\ell\gamma$, where $\ell$ is a charged lepton, as functions of the energy of the photon in the rest frame of the $D_s$ meson. The computation is performed using gauge-field configurations with 2+1+1 sea-quark flavours generated by the European Twisted Mass Collaboration and the results have been extrapolated to the continuum limit. For the vector form factor we observe a very significant partial cancellation between the contributions from the emission of the photon from the strange quark and that from the charm quark. The results for the form factors are used to test the reliability of various Anz\"atze based on single-pole dominance and its extensions, and we present a simple parametrization of the form factors which fits our data very well and which can be used in future phenomenological analyses. Using the form factors we compute the differential decay rate and the branching ratio for the process $D_s\to e\nu_e\gamma$ as a function of the lower cut-off on the photon energy. With a cut-off of 10 MeV for example, we find a branching ratio of Br$(E_\gamma>10\,\mathrm{MeV})=4.4(3)\times 10^{-6}$ which, unlike some model calculations, is consistent with the upper bound from the BESIII experiment Br$(E_\gamma>10\,\mathrm{MeV}), Comment: 31 pages, 14 figures

This work presents the first lattice-QCD calculation of the twist-3 axial quark GPDs for the proton using the large-momentum effective theory approach. We calculate matrix elements with momentum-boosted proton states and a non-local axial operator. The calculation is performed using one ensemble of two degenerate light, a strange and a charm quark ($N_f=2+1+1$) of maximally twisted mass fermions with a clover term. The ensemble has a volume $32^3\times64$, lattice spacing 0.0934 fm, and corresponds to a pion mass of 260 MeV. The matrix elements are calculated for three values of the proton momentum, namely 0.83, 1.25, and 1.67 GeV. The light-cone GPDs are defined in the symmetric frame, which we implement here with a (negative) 4-momentum transfer squared of 0.69, 1.38, and 2.76 GeV$^2$, all at zero skewness. We also conduct several consistency checks, including assessing the local limit of the twist-3 GPDs and examining the Burkhardt-Cottingham-type as well as Efremov-Teryaev-Leader-type sum rules., 22 pages, 17 figures

The Refined Gribov-Zwanziger (RGZ) action in the Landau gauge accounts for the existence of infinitesimal Gribov copies as well as the dynamical formation of condensates in the infrared of Euclidean Yang-Mills theories. We couple scalar fields to the RGZ action and compute the one-loop scalar propagator in the adjoint representation of the gauge group. We compare our findings with existing lattice data. The fate of BRST symmetry in this model is discussed, and we provide a comparison to a previous proposal for a non-minimal coupling between matter and the RGZ action. We find good agreement with the lattice data of the scalar propagator for the values of the mass parameters that fit the RGZ gluon propagator to the lattice. This suggests that the non-perturbative information carried by the gluon propagator in the RGZ framework provides a suitable mechanism to reproduce the behavior of correlation functions of colored matter fields in the infrared., 18 pages + refs.; 6 figures

We provide evidence for partial deconfinement -- the deconfinement of a SU($M$) subgroup of the SU($N$) gauge group -- by using lattice Monte Carlo simulations. We take matrix models as concrete examples. By appropriately fixing the gauge, we observe that the $M\times M$ submatrices deconfine. This gives direct evidence for partial deconfinement at strong coupling. We discuss the applications to QCD and holography., Comment: 49 pages, 24 figures, 3 tables. Version accepted for publication on JHEP

The gauge-invariant two-point function of the Higgs field at the same spacetime point can make a natural gauge-invariant order parameter for spontaneous gauge symmetry breaking. However, this composite operator is ultraviolet divergent and is not well defined. We propose using a gradient flow to cure the divergence from putting the fields at the same spacetime point. As a first step, we compute it for the Abelian Higgs model with a positive mass squared at the one-loop order in the continuum theory using the saddle-point method to estimate the finite part. The order parameter consistently goes to zero in the infrared limit of the infinite flow time., Comment: 37 pages, 2 figures, minor changes, published version on EPJC

The complex Langevin (CL) method is a promising approach to overcome the sign problem that occurs in real-time formulations of quantum field theories. Using the Schwinger-Keldysh formalism, we study SU($N_c$) gauge theories with CL. We observe that current stabilization techniques are insufficient to obtain correct results. Therefore, we revise the discretization of the CL equations on complex time contours, find a time reflection symmetric formulation and introduce a novel anisotropic kernel that enables CL simulations on discretized complex time paths. Applying it to SU(2) Yang-Mills theory in 3+1 dimensions, we obtain unprecedentedly stable results that we validate using additional observables and that can be systematically improved. For the first time, we are able to simulate non-Abelian gauge theory on time contours whose real-time extent exceeds its inverse temperature. Thus, our approach may pave the way towards an ab-initio real-time framework of QCD in and out of equilibrium with a potentially large impact on the phenomenology of heavy-ion collisions., 28 pages + references and appendices, 7 figures, 1 table; v2: revised, data and plot files attached, published version

The theory of the strong force, quantum chromodynamics, describes the proton in terms of quarks and gluons. The proton is a state of two up quarks and one down quark bound by gluons, but quantum theory predicts that in addition there is an infinite number of quark-antiquark pairs. Both light and heavy quarks, whose mass is respectively smaller or bigger than the mass of the proton, are revealed inside the proton in high-energy collisions. However, it is unclear whether heavy quarks also exist as a part of the proton wavefunction, which is determined by non-perturbative dynamics and accordingly unknown: so-called intrinsic heavy quarks. It has been argued for a long time that the proton could have a sizable intrinsic component of the lightest heavy quark, the charm quark. Innumerable efforts to establish intrinsic charm in the proton have remained inconclusive. Here we provide evidence for intrinsic charm by exploiting a high-precision determination of the quark-gluon content of the nucleon based on machine learning and a large experimental dataset. We disentangle the intrinsic charm component from charm-anticharm pairs arising from high-energy radiation. We establish the existence of intrinsic charm at the 3-standard-deviation level, with a momentum distribution in remarkable agreement with model predictions. We confirm these findings by comparing to very recent data on Z-boson production with charm jets from the LHCb experiment., Comment: 10 pages, 4 figures + Supplementary Information, published version available in https://www.nature.com/articles/s41586-022-04998-2. v2: updated comparison with the MBC model in Fig. 1 right (thanks to Tim Hobbs for the suggestion)

Gauge theories form the foundation of modern physics, with applications ranging from elementary particle physics and early-universe cosmology to condensed matter systems. We perform quantum simulations of the unitary dynamics of a U(1) symmetric gauge field theory and demonstrate emergent irreversible behavior. The highly constrained gauge theory dynamics is encoded in a one-dimensional Bose--Hubbard simulator, which couples fermionic matter fields through dynamical gauge fields. We investigate global quantum quenches and the equilibration to a steady state well approximated by a thermal ensemble. Our work may enable the investigation of elusive phenomena, such as Schwinger pair production and string-breaking, and paves the way for simulating more complex higher-dimensional gauge theories on quantum synthetic matter devices., Accepted version, 6+6 pages, 4+5 figures

A partial width formula is proposed using the analytical extension of the wave function in momentum space. The distinction of the Riemann sheets is explained from the perspective of the Schrodinger equation. The analytical form in coordinate space and the partial width are derived subsequently. Then a coupled-channel analysis is performed to investigate the open-charm branching ratios of the $P_c$ states, involving the contact interactions and one-pion-exchange potential with the three-body effects. The low energy constants are fitted using the experimental masses and widths as input. The $P_c(4312)$ is found to decay mainly to $\Lambda_c\bar{D}^*$, while the branching ratios of the $P_c(4440)$ and $P_c(4457)$ in different channels are comparable. Under the reasonable assumption that the off-diagonal contact interactions are small, the $J^P$ quantum numbers of the $P_c(4440)$ and the $P_c(4457)$ prefer $\frac{1}{2}^-$ and $\frac{3}{2}^-$ respectively. Three additional $P_c$ states at 4380 MeV, 4504 MeV and 4516 MeV, together with their branching ratios, are predicted. A deduction of the revised one-pion-exchange potential involving the on-shell three-body intermediate states is performed., Comment: 16 pages, 5 figures

Off-central heavy-ion collisions are known to feature magnetic fields with magnitudes and characteristic gradients corresponding to the scale of the strong interactions. In this work, we employ equilibrium lattice simulations of the underlying theory, QCD, involving similar inhomogeneous magnetic field profiles to achieve a better understanding of this system. We simulate three flavors of dynamical staggered quarks with physical masses at a range of magnetic fields and temperatures, and extrapolate the results to the continuum limit. Analyzing the impact of the field on the quark condensate and the Polyakov loop, we find non-trivial spatial features that render the QCD medium qualitatively different as in the homogeneous setup, especially at temperatures around the transition. In addition, we construct leading-order chiral perturbation theory for the inhomogeneous background and compare its prediction to our lattice results at low temperature. Our findings will be useful to benchmark effective theories and low-energy models of QCD for a better description of peripheral heavy-ion collisions., 24 pages, 15 figures

A novel approach is proposed to link the charge symmetry breaking (CSB) nuclear interaction and the low-energy constants in quantum chromodynamics (QCD) for the first time by matching the CSB effect in nuclear matter. Resulting CSB interaction is applied to study the Okamoto-Nolen-Schiffer anomaly on the energy differences of mirror nuclei by taking $ {}^{17} \mathrm{F} $-$ {}^{17} \mathrm{O} $, $ {}^{15} \mathrm{O} $-$ {}^{15} \mathrm{N} $, $ {}^{41} \mathrm{Sc} $-$ {}^{41} \mathrm{Ca} $, and $ {}^{39} \mathrm{Ca} $-$ {}^{39} \mathrm{K} $ as typical examples. The magnitude and sign of the QCD-based CSB interactions are found to resolve the anomaly successfully within theoretical uncertainties., 5 pages, 2 figures, 6 tables in the main text; 5 pages, 2 figures, 2 tables in the supplemental material

We compute the sphaleron rate in quenched QCD for a temperature $T \simeq 1.24~T_c$ from the inversion of the Euclidean lattice time correlator of the topological charge density. We explore and compare two different strategies: one follows a new approach proposed in this study and consists in extracting the rate from finite lattice spacing correlators, and then in taking the continuum limit at fixed smoothing radius followed by a zero-smoothing extrapolation; the other follows the traditional approach of extracting the rate after performing such double extrapolation directly on the correlator. In both cases the rate is obtained from a recently-proposed modification of the standard Backus-Gilbert procedure. The two strategies lead to compatible estimates within errors, which are then compared to previous results in the literature at the same or similar temperatures; the new strategy permits to obtain improved results, in terms of statistical and systematic uncertainties., 11 pages, 12 eps figures

We investigate the $Z_b$, $Z_c$ and $Z_{cs}$ states within the chiral effective field theory framework and the $S$-wave single channel molecule picture. With the complex scaling method, we accurately solve the Schr\"odinger equation in momentum space. Our analysis reveals that the $Z_b(10610)$, $Z_b(10650)$, $Z_c(3900)$ and $Z_c(4020)$ states are the resonances composed of the $S-$wave $(B\bar{B}^{*}+B^{*}\bar{B})/\sqrt{2}$, $B^{*}\bar{B}^*$, $(D\bar{D}^{*}+D^{*}\bar{D})/\sqrt{2}$ and $D^{*}\bar{D}^*$, respectively. Furthermore, although the $Z_{cs}(3985)$ and $Z_{cs}(4000)$ states exhibit a significant difference in width, these two resonances may originate from the same channel, the $S-$wave $(D_{s}\bar{D}^{*}+D_{s}^{*}\bar{D})/\sqrt{2}$. Additionally, we find two resonances in the $S-$wave $D_s^*\bar{D}^*$ channel, corresponding to the $Z_{cs}(4123)$ and $Z_{cs}(4220)$ states that await experimental confirmation., Comment: 10 pages, 5 figures, 4 tables

We present a nonperturbative recipe for directly computing the $S$-matrix in strongly-coupled QFTs. The method makes use of spectral data obtained in a Hamiltonian framework and can be applied to a wide range of theories, including potentially QCD. We demonstrate the utility of this prescription in the specific example of the 2+1d $O(N)$ model at large $N$, using energy eigenstates computed with Hamiltonian truncation to reproduce the full $2 \to 2$ scattering amplitude for arbitrary (complex) center-of-mass energy., 28 pages with appendices, 12 figures

We address the nonperturbative calculation of the inclusive decay rate of semileptonic $B_{(s)}$-meson decays from lattice QCD. Precise Standard-Model predictions are key ingredients in searches for new physics, and this type of computation may eventually provide new insight into the long-standing tension between the inclusive and exclusive determinations of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements $|V_{cb}|$ and $|V_{ub}|$. We present results from a pilot lattice computation for $B_s \rightarrow X_c\, l \nu_l$, where the initial $b$ quark described by the relativistic-heavy-quark (RHQ) formalism on the lattice and the other valence quarks discretised with domain-wall fermions are simulated approximately at their physical quark masses. We compare two different methods for computing the decay rate from lattice data of Euclidean $n$-point functions, namely Chebyshev and Backus-Gilbert approaches. We further study how much the ground-state meson dominates the inclusive decay rate and indicate our strategy towards a computation with a more comprehensive systematic error budget., Comment: 43 pages, 18 figures

Using proton-proton collisions at centre-of-mass energies $7$, $8$, and $13$ TeV, with a total integrated luminosity of $9\,\text{fb}^{-1}$, the LHCb collaboration has performed amplitude analyses of the $B^+\to D^+D^-K^+$, $B^+\to D^- D_s^+ \pi^+$ and $B^0\to \bar{D}^0 D_s^+ \pi^-$ decays, observing that new $T_{cs}$ and $T_{c\bar s}$ resonances are required in order to explain the experimental data. These signals could be the first observation of tetraquark candidates that do not contain a heavy quark-antiquark pair; in fact, they consist of four different flavours of quarks, one of which is a doubly charged open-charm state. We present herein an analysis of the $T_{cs}$ and $T_{c\bar s}$ states, which is an extension of our recently published study of similar $T_{cc}^+$ exotic candidates. Our theoretical framework is a constituent-quark-model-based coupled-channels calculation of $qq^\prime \bar s \bar c$ and $cq\bar s\bar q^{\prime}$ tetraquark sectors for $T_{cs}$ and $T_{c\bar s}$ structures, respectively. We explore the nature, and pole position, of the singularities that appear in the scattering matrix with spin-parity quantum numbers: $J^P=0^\pm$, $1^\mp$, and $2^\pm$. The constituent quark model has been widely used in the heavy quark sector, and thus all model parameters are already constrained from previous works. This makes our predictions robust and parameter-free. We find many singularities in the solution of various scattering-matrix problems which are either virtual states or resonances, but not bound states. Some of them fit well with the experimental observations of the spin-parity, mass and width of $T_{cs}$ and $T_{c\bar s}$ candidates, and thus tentative assignments are made; however, with caution, because the experimental Breit-Wigner parameters are related to the pole characteristics., Comment: 10 pages, 4 tables

Staple-shaped Wilson line operators are necessary for the study of transverse momentum-dependent parton distribution functions (TMDPDFs) in lattice QCD and beyond. In this work, we study the renormalization of such operators in the general case of an asymmetric staple. We analyze the mixing pattern of these operators using their symmetry properties, where we find that the possible mixing is restricted within groups of four operators. We then present numerical results using the regularization independent scheme, RI/MOM, to study the importance of mixing using one operator in particular, the $\gamma_0$ operator. Based on these results, we consider the short distance ratio (SDR) scheme, which is desirable in the absence of mixing. Finally, we investigate a variant of the RI/MOM scheme, where the renormalization factors are computed at short distances., Comment: 18 pages, 14 figures

Quantum simulation has become a promising avenue of research that allows one to simulate and gain insight into the models of High Energy Physics whose experimental realizations are either complicated or inaccessible with current technology. We demonstrate the quantum simulation of such a model, a quantum mechanical system with spatial noncommutativity, which is inspired by the works in Noncommutative Geometry and Noncommutative Field theory for a universal quantum computer. We use the novel group theoretical formalism to map the Hamiltonian of such a noncommutative quantum system into the ordinary quantum mechanical Hamiltonian and then carry out the quantum simulation using the Trotter-Suzuki product formula. Furthermore, we distinguish the impact of the noncommutativity parameter on the quantum simulation, especially on the Trotter error, and point out how its sizable value affects the simulation., 15 pages plus references, Fig. 2 redrawn and Fig. 5 repositioned, typos in Eq.(26) corrected

For the first time, we present a lattice QCD determination of Mellin moments of unpolarized generalized parton distributions (GPDs) of the proton from an analysis of the quasi-GPD matrix elements within the short-distance factorization framework. We perform our calculation on an $N_f$=2+1+1 twisted mass fermions ensemble with a clover improvement at lattice spacing $a=0.093$ fm and a pion mass of $m_\pi=260$ MeV. Focusing on the zero-skewness case, the iso-vector and iso-scalar quasi-GPDs are calculated from the $\gamma_0$ definition, as well as a recently proposed Lorentz-invariant definition. We utilize data on both symmetric and asymmetric kinematic frames, which allows us to obtain the Mellin moments for several values of the momentum transfer, $-t$, in the range 0.17 to $2.77~\rm{GeV}^2$. We use the ratio scheme for GPDs, i.e. renormalization group invariant ratios with leading-twist factorization formula and perturbatively calculated matching coefficients up to the next-next-to-leading order (NNLO) to extract Mellin moments of GPDs, which are consistent with renormalization-group improved results. We compare our determination from quasi-GPDs with the results extracted using standard calculations of Mellin moments of local operators, specifically those related to the electromagnetic and gravitational form factors. We estimated the moments of GPDs up to the fifth ones for the first time. By extrapolating the Mellin moments to $-t=0$, we obtained the quark charges, momentum fraction, as well as the angular momentum contributions to the proton spin. The impact parameter space interpretation of the GPD moments is discussed, which provides insights into the spatial distribution of unpolarized quarks and their correlations in the transverse plane of an unpolarized or transversely polarized proton., Comment: 28 pages, 22 figures