Search

We investigate an improved method to extract nucleon matrix elements from lattice 3-point functions using a generalized eigenvalue problem (GEVP) with nucleon and pion-nucleon interpolating fields. Our method avoids the computation of the costly three-point functions that have pion-nucleon interpolators at both source and sink. We demonstrate that excited state contamination from $N\pi$ is minimized in nucleon matrix elements of the scalar, vector, pseudoscalar, axial, and tensor currents and discuss our results based on a physical-point ensemble with a pion mass value of 131 MeV. We find that the GEVP is most significant for the isovector pseudoscalar and axial currents., Comment: 9 pages, 5 figures, 41th International Symposium on Lattice Field Theory - Lattice2024

We determine the nucleon axial, scalar and tensor charges and the nucleon $\sigma$-terms using twisted mass fermions. We employ three ensembles with approximately equal physical volume of about 5.5~fm, three values of the lattice spacing, approximately 0.06~fm, 0.07~fm and 0.08~fm, and with the mass of the degenerate up and down, strange and charm quarks tuned to approximately their physical values. We compute both isovector and isoscalar charges and $\sigma$-terms and their flavor decomposition including the disconnected contributions. We use the Akaike Information Criterion to evaluate systematic errors due to excited states and the continuum extrapolation. For the nucleon isovector axial charge we find $g_A^{u-d}=1.250(24)$, in agreement with the experimental value. Moreover, we extract the nucleon $\sigma$-terms and find for the light quark content $\sigma_{\pi N}=41.9(8.1)$~MeV and for the strange $\sigma_{s}=30(17)$~MeV., Comment: 18 pages, 17 figures

This work studies the initial-boundary value problem for both the linear Schr\"odinger equation and the cubic nonlinear Schr\"odinger equation on the half-space in higher dimensions ($n\ge 2$). First, the forced linear problem is solved on the half-space via the Fokas method and then using the obtained solution formula new and interesting linear estimates are derived with data and forcing in appropriate spaces. Second, the well-posedness of the nonlinear problem on the half-space is proved with initial data in Sobolev spaces $H^s(\mathbb{R}^n_+)$, with $s>\frac{n}{2}-1$, and boundary data in natural Bourgain spaces $\mathcal{B}^s$ that reflect the boundary regularity of the linear problem. The proof method consists of showing that the iteration map defined via the Fokas solution formula is a contraction by establishing sharper trilinear estimates. The presence of the boundary introduces solution spaces that involve temporal Bourgain spaces., Comment: 39 pages, 4 figures

We present a lattice calculation of the Hadronic Vacuum Polarization (HVP) contribution of the strange and charm quarks to the anomalous magnetic moment of the muon in isospin symmetric QCD. We employ the gauge configurations generated by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks at five lattice spacings and at values of the quark mass parameters that are close and/or include the isospin symmetric QCD point of interest. After computing the small corrections necessary to precisely match this point, and carrying out an extrapolation to the continuum limit based on the data at lattice spacings $a \simeq 0.049, 0.057, 0.068, 0.080$~fm and spatial lattice sizes up to $L \simeq 7.6$~fm, we obtain $a_\mu^{\rm HVP}(s) = (53.57 \pm 0.63) \times 10^{-10}$ and $a_\mu^{\rm HVP}(c) = (14.56 \pm 0.13) \times 10^{-10}$, for the quark-connected strange and charm contributions, respectively. Our findings agree well with the corresponding results by other lattice groups., Comment: 33 pages, 18 figures, 3 tables

We give a brief overview of recent progress in lattice QCD simulations that is enabling precision studies of the three-dimensional structure of the nucleon. We present results on nucleon charges and second Mellin moments of parton distribution functions and generalized parton distributions, highlighting results obtained at the continuum limit using only gauge ensembles simulated with physical quark masses. The tensor charge and transversity moments are determined with controlled systematics and used to extract nucleon transverse density distributions. We also discuss progress towards the direct evaluation of generalized parton distributions in lattice QCD and the impact they are having on phenomenology., Comment: 9 pages, 7 figures, Invited Talk at 7th International Workshop on "Transverse phenomena in hard processes and the transverse structure of the proton (Transversity2024)", 3-7 June 2024, Trieste, Italy

We investigate contributions of excited states to nucleon matrix elements computed in lattice QCD by employing, in addition to the standard nucleon interpolating operator, pion-nucleon ($\pi$-$N$) operators. We solve a generalized eigenvalue problem (GEVP) to obtain an optimal interpolating operator that minimizes overlap with the $\pi$-$N$ states. We derive a variant of the standard application of the GEVP method, which allows for constructing 3-point correlation functions using the optimized interpolating operator without requiring the computationally demanding combination that includes $\pi$-$N$ operators in both sink and source. We extract nucleon matrix elements using two twisted mass fermion ensembles, one ensemble generated using pion mass of 346 MeV and one ensemble tuned to reproduce the physical value of the pion mass. Especially, we determine the isoscalar and isovector scalar, pseudoscalar, vector, axial, and tensor matrix elements. We include results obtained using a range of kinematic setups, including momentum in the sink. Our results using this variational approach are compared with previous results obtained using the same ensembles and multi-state fits without GEVP improvement. We find that for the physical mass point ensemble, the improvement, in terms of suppression of excited states using this method, is most significant for the case of the matrix elements of the isovector axial and pseudoscalar currents., Comment: 37 pages, 39 figures

For a large class of cohomology theories, we prove that refined unramified cohomology is canonically isomorphic to the hypercohomology of a natural truncated complex of Zariski sheaves. This generalizes a classical result of Bloch and Ogus and solves a conjecture of Kok and Zhou., Comment: 21 pages, final version, to appear in Advances in Mathematics, updated some references and added Corollaries 1.4, 5.1, and 5.2

We present the full decomposition of the momentum fraction carried by quarks and gluons in the pion and the kaon. We employ three gauge ensembles generated with $N_f=2+1+1$ Wilson twisted-mass clover-improved fermions at the physical quark masses. For both mesons we perform a continuum extrapolation directly at the physical pion mass, which allows us to determine for the first time the momentum decomposition at the physical point. We find that the total momentum fraction carried by quarks is 0.532(56) and 0.618(32) and by gluons 0.388(49) and 0.408(61) in the pion and in the kaon, respectively, in the $\overline{\mathrm{MS}}$ scheme and at the renormalization scale of 2 GeV. Having computed both the quark and gluon contributions in the continuum limit, we find that the momentum sum is 0.926(68) for the pion and 1.046(90) for the kaon, verifying the momentum sum rule., Comment: 14 pages, 10 figures

We study the $\bar b \bar b u d$ tetraquark with quantum numbers $I(J^P) = 0(1^+)$ as well as the $\bar b \bar b u s$ tetraquark with quantum numbers $J^P = 1^+$ using lattice QCD. We improve on existing work by including both local and scattering interpolating operators on both sides of the correlation functions and use symmetric correlation matrices. This allows not only a reliable determination of the energies of QCD-stable tetraquark ground states, but also of low-lying excited states, which are meson-meson scattering states. The latter is particularly important for future finite-volume scattering analyses. Here, we perform chiral and continuum extrapolations of just the ground-state energies, for which finite-volume effects are expected to be small. Our resulting tetraquark binding energies, $-100 \pm 10\:^{+36}_{-51}\:\:{\rm MeV}$ for $\bar b \bar b u d$ and $-30 \pm 3\:^{+11}_{-31}\:\:{\rm MeV}$ for $\bar b \bar b u s$, are consistent with other recent lattice-QCD predictions., Comment: 28 pages, 15 figures; version published in Phys. Rev. D

We present a lattice determination of the inclusive decay rate of the process $\tau\mapsto X_{us} \nu_\tau$ in which the $\tau$ lepton decays into a generic hadronic state $X_{us}$ with $\bar u s$ flavour quantum numbers. Our results have been obtained in $n_f=2+1+1$ iso-symmetric QCD with full non-perturbative accuracy, without any OPE approximation and, except for the presently missing long-distance isospin-breaking corrections, include a solid estimate of all sources of theoretical uncertainties. This has been possible by using the Hansen-Lupo-Tantalo method [1] that we have already successfully applied in [2] to compute the inclusive decay rate of the process $\tau\mapsto X_{ud} \nu_\tau$ in the $\bar u d$ flavour channel. By combining our first-principles theoretical results with the presently-available experimental data we extract the CKM matrix element $\vert V_{us}\vert$, the Cabibbo angle, with a $0.9$\% accuracy, dominated by the experimental error., Comment: 9 pages, 5 figures, 1 table. Version accepted by PRL, expanded technical discussion moved to an appendix, results unchanged

We compare two cycle class maps on torsion cycles and show that they agree up to a minus sign. The first one goes back to Bloch (1979), with recent generalizations to non-closed fields. The second is the \'etale motivic cycle class map $\alpha^{i}_{X}\colon \text{CH}^{i}(X)_{\mathbb{Z}_{\ell}}\to H^{2i}_{L}(X,\mathbb{Z}_{\ell}(i))$ restricted to torsion cycles., Comment: 17 pages

We present an update on our lattice calculations of the Mellin moments of PDFs and GPDs for the pion and kaon, using momentum-boosted meson states. In particular, we focus on the calculation of the scalar and tensor local operators, and the vector operator with up to three-covariant derivatives. The corresponding matrix elements allow us to extract the scalar and tensor charges, as well as $\langle x^n \rangle$ with $n=1,2,3$. In addition, we introduce momentum transfer between the initial and final meson state that leads to the scalar, vector, and tensor form factors, and the generalized form factors up to three covariant derivatives. The above results are obtained using two ensembles of maximally twisted mass fermions with clover improvement with two degenerate light, a strange, and a charm quark $(N_f=2+1+1)$ at lattice spacings of 0.093 fm and 0.081 fm. The pion mass of the ensembles is about 260 MeV. We study excited-states effects by analyzing four values of the source-sink time separation (1.12 - 1.67 fm). We also examine discretization and volume effects. The lattice data are renormalized non-perturbatively, and the results are presented in the MS scheme at a scale of 2 GeV., Comment: 11 pages, 7 figures, Contribution to The 40th International Symposium on Lattice Field Theory (Lattice 2023)

State-of-the-art lattice QCD simulations enable the evaluation of nucleon form factors and Mellin moments with controlled systematics, yielding results with unprecedented accuracy. At the same time, new theoretical approaches are allowing the direct computation of nucleon generalized parton distributions. We review recent lattice QCD results on these quantities that are paving the way for extracting a wealth of information on the 3-dimensional structure of the nucleon., Comment: Two figures updated and typos corrected. 8 pages and 7 figures. Invited talk at the 16th edition of the "International Conference on Meson-Nucleon Physics and the Structure of the Nucleon" (MENU 2023), October 16-20, 2023, Institute for Nuclear Physics of the Johannes Gutenberg University of Mainz, Germany. arXiv admin note: text overlap with arXiv:2309.05774

We investigate contributions of excited states to nucleon matrix elements by studying the two- and three-point functions using nucleon and pion-nucleon interpolating fields. This study is carried out using twisted mass fermion ensembles with pion masses 346 MeV and 131 MeV. We compare the results obtained using these two ensembles and show preliminary results for nucleon charges., Comment: 8 pages, 3 figures, 40th International Symposium on Lattice Field Theory - Lattice2023

We present the first determination of the energy dependence of the $B$-$\bar{D}$ and $B^*$-$\bar{D}$ isospin-0, $S$-wave scattering amplitudes both below and above the thresholds using lattice QCD, which allows us to investigate rigorously whether mixed bottom-charm $\bar{b}\bar{c}ud$ tetraquarks exist as bound states or resonances. The scattering phase shifts are obtained using L\"uscher's method from the energy spectra in two different volumes. To ensure that no relevant energy level is missed, we use large, symmetric $7 \times 7$ and $8 \times 8$ correlation matrices that include, at both source and sink, $B^{(*)}$-$\bar{D}$ scattering operators with the lowest three or four possible back-to-back momenta in addition to local $\bar{b}\bar{c}ud$ operators. We fit the energy dependence of the extracted scattering phase shifts using effective-range expansions. We observe sharp peaks in the $B^{(*)}$-$\bar{D}$ scattering rates close to the thresholds, which are associated with shallow bound states, either genuine or virtual, a few MeV or less below the $B^{(*)}$-$\bar{D}$ thresholds. In addition, we find hints for resonances with masses of order $100$ MeV above the thresholds and decay widths of order $200$ MeV., Comment: 5 pages + 10 pages of references and supplemental material. Added estimates of bound-state probabilities and improved calculation of bound-state mass uncertainties

In this paper, we present lattice QCD results for the $x$-dependence of the unpolarized gluon PDF for the proton. We use one ensemble of $N_f=2+1+1$ maximally twisted mass fermions with a clover improvement, and the Iwasaki improved gluon action. The quark masses are tuned to produce a pion with a mass of 260 MeV. The ensemble has a lattice spacing of $a=0.093$ fm and a spatial extent of 3 fm. We employ the pseudo-distribution approach, which relies on matrix elements of non-local operators that couple to momentum-boosted hadrons. In this work, we use five values of the momentum boost between 0 and 1.67 GeV. The gluon field strength tensors of the non-local operator are connected with straight Wilson lines of varying length $z$. The light-cone Ioffe time distribution (ITD) is extracted utilizing data with $z$ up to 0.56 fm and a quadratic parametrization in terms of the Ioffe time at fixed values of $z$. We explore systematic effects, such as the effect of the stout smearing for the gluon operator, excited states effects, and the dependence on the maximum value of $z$ entering the fits to obtain the gluon PDF. Also, for the first time, the mixing with the quark singlet PDFs is eliminated using matrix elements with non-local quark operators that were previously analyzed within the quasi-PDF framework on the same ensemble. Here, we expand the data set for the quark singlet and reanalyze within the pseudo-PDFs method eliminating the corresponding mixing in the gluon PDF.

We compute the nucleon axial and pseudoscalar form factors using three $N_f=$2+1+1 twisted mass fermion ensembles with all quark masses tuned to approximately their physical values. The values of the lattice spacings of these three physical point ensembles are 0.080 fm, 0.068 fm, and 0.057 fm, and spatial sizes 5.1 fm, 5.44 fm, and 5.47 fm, respectively, yielding $m_\pi L$>3.6. Convergence to the ground state matrix elements is assessed using multi-state fits. We study the momentum dependence of the three form factors and check the partially conserved axial-vector current (PCAC) hypothesis and the pion pole dominance (PPD). We show that in the continuum limit, the PCAC and PPD relations are satisfied. We also show that the Goldberger-Treimann relation is approximately fulfilled and determine the Goldberger-Treiman discrepancy. We find for the nucleon axial charge $g_A$=1.245(28)(14), for the axial radius $\langle r^2_A \rangle$=0.339(48)(06) fm$^2$, for the pion-nucleon coupling constant $g_{\pi NN} \equiv \lim_{Q^2 \rightarrow -m_\pi^2} G_{\pi NN}(Q^2)$=13.25(67)(69) and for $G_P(0.88m_{\mu}^2)\equiv g_P^*$=8.99(39)(49)., Comment: accepted version, 32 pages, 39 figures, 14 tables

We investigate the low-lying baryon spectrum using three $N_f$=2+1+1 ensembles simulated with physical values of the quark masses and lattice spacings of 0.080, 0.069, and 0.057 fm. The ensembles are generated using twisted mass clover-improved fermions and the Iwasaki gauge action. The spatial length is kept approximately the same at about 5.1 fm to 5.5 fm fulfilling the condition $m_\pi L$> 3.6. We investigate isospin splitting within isospin multiples and verify that for most cases the isospin splitting for these lattice spacing is consistent with zero. In the couple of cases, for which there is a non-zero value, in the continuum limit, the mass splitting goes to zero. The baryon masses are extrapolated to the continuum limit using the three $N_f$=2+1+1 ensembles and are compared to other recent lattice QCD results. For the strange and charm quark masses, we find, respectively, $m_s$(2 GeV)=99.2(2.7) MeV and $m_c$(3 GeV)=1.015(39) GeV. The values predicted for the masses of the doubly charmed $\Xi_{cc}^\star$, $\Omega_{cc}$ and $\Omega_{cc}^\star$ baryons are 3.676(55) GeV, 3.703(51) GeV and 3.803(50) GeV, respectively, and for the triply charmed $\Omega_{ccc}$ baryon is 4.785(71) GeV., Comment: accepted version, 15 pages, 18 figures, 10 tables

The neutral pion generates the leading pole contribution to the hadronic light-by-light tensor, which is given in terms of the nonperturbative transition form factor $\mathcal{F}_{\pi^0\gamma\gamma}(q_1^2,q_2^2)$. Here we present an ab-initio lattice calculation of this quantity in the continuum and at the physical point using twisted-mass lattice QCD. We report our results for the transition form factor parameterized using a model-independent conformal expansion valid for arbitrary space-like kinematics and compare it with experimental measurements of the single-virtual form factor, the two-photon decay width, and the slope parameter. We then use the transition form factors to compute the pion-pole contribution to the hadronic light-by-light scattering in the muon $g-2$, finding $a_\mu^{\pi^0\text{-pole}} = 56.7(3.2) \times 10^{-11}$., Comment: 21 pages, 17 figures, 4 tables, updated to published version

We present an investigation of pion-nucleon elastic scattering in the $I\,(J^P) = \frac{3}{2}\,(\frac{3}{2}^+)$ channel using lattice QCD with degenerate up and down, strange and charm quarks with masses tuned to their physical values. We use an ensemble of twisted mass fermions with box size $L = 5.1\,\mathrm{fm}$ and lattice spacing $a = 0.08\,\mathrm{fm}$ and we consider the $\pi N$ system in rest and moving frames up to total momentum $\vec{P}^2 = 3\,(2\pi/L)^2$ = 0.17 GeV$^2$. We take into account the finite volume symmetries and $S$- and $P$-wave mixing, and use the L\"uscher formalism to simultaneously constrain the $J = 1/2,\,\ell = 0$ and $J = 3/2,\,\ell = 1$ scattering amplitudes. We estimate the $\Delta$ resonance pole in the $P$-wave channel as well as the $S$-wave isospin-3/2 scattering length., Comment: 18 pages, 18 figures, 10 tables; v2 matches published version

Quantum computers offer an intriguing path for a paradigmatic change of computing in the natural sciences and beyond, with the potential for achieving a so-called quantum advantage, namely a significant (in some cases exponential) speed-up of numerical simulations. The rapid development of hardware devices with various realizations of qubits enables the execution of small scale but representative applications on quantum computers. In particular, the high-energy physics community plays a pivotal role in accessing the power of quantum computing, since the field is a driving source for challenging computational problems. This concerns, on the theoretical side, the exploration of models which are very hard or even impossible to address with classical techniques and, on the experimental side, the enormous data challenge of newly emerging experiments, such as the upgrade of the Large Hadron Collider. In this roadmap paper, led by CERN, DESY and IBM, we provide the status of high-energy physics quantum computations and give examples for theoretical and experimental target benchmark applications, which can be addressed in the near future. Having the IBM 100 x 100 challenge in mind, where possible, we also provide resource estimates for the examples given using error mitigated quantum computing.

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 momentum subtraction (RI/MOM) scheme 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: 20 pages, 14 figures

We show that for any integer $n\geq2$ there is a smooth complex projective variety $X$ of dimension $5$ whose third Griffiths group $\text{Griff}^{3}(X)$ contains infinitely many torsion elements of order $n$. This generalises a recent theorem of Schreieder who proved the result for $n=2$., Comment: 25 pages

The ratio $R(E)$ of the cross-sections for $e^+e^-\to$ hadrons and $e^+e^-\to \mu^+\mu^-$ is a valuable energy-dependent probe of the hadronic sector of the Standard Model. Moreover, the experimental measurements of $R(E)$ are the inputs of the dispersive calculations of the leading hadronic vacuum polarization contribution to the muon $g-2$ and these are in significant tension with direct lattice calculations and with the muon $g-2$ experiment. In this talk we discuss the results of our first-principles lattice study of $R(E)$. By using a recently proposed method for extracting smeared spectral densities from Euclidean lattice correlators, we have calculated $R(E)$ convoluted with Gaussian kernels of different widths $\sigma$ and central energies up to $2.5$ GeV. Our theoretical results have been compared with the KNT19 [1] compilation of experimental results smeared with the same Gaussian kernels and a tension (about three standard deviations) has been observed for $\sigma\sim 600$ MeV and central energies around the $\rho$-resonance peak., Comment: Contribution to the 39th International Symposium on Lattice Field Theory (Lattice2022), 8th-13th August, 2022, Bonn, Germany

We present results of the x-dependence of the unpolarized gluon PDF for the proton. We use an $N_f=2+1+1$ ensemble of maximally twisted mass fermions with clover improvement and the Iwasaki improved gluon action. The quark masses are tuned so that the pion mass is 260 MeV. We use a $32^3\times64$ lattice size with a lattice spacing $a=0.093$ fm giving a spatial extent of 3 fm. We employ the pseudo-distribution approach and obtain the light-cone Ioffe time distribution (ITD) combining data for nucleon momentum boosts up to 1.67 GeV and Wilson line lengths, $z$, up to 0.56 fm. We explore systematic effects such as the dependence on the maximum value of $z$ entering the fits to obtain the gluon PDF.

Pseudoscalar-pole diagrams are an important component of estimates of the hadronic light-by-light (HLbL) contribution to the muon $g-2$. We report on our computation of the transition form factors $\mathcal{F}_{P \rightarrow \gamma^* \gamma^*}$ for the neutral pseudoscalar mesons $P=\pi^0$ and $\eta$. The calculation is performed using twisted-mass lattice QCD with physical quark masses. On the lattice, we have access to a broad range of (space-like) photon four-momenta and therefore produce form factor data complementary to the experimentally accessible single-virtual direction, which directly leads to an estimate of the pion- and $\eta$-pole components of the muon $g-2$. For the pion, our result for the $g-2$ contribution in the continuum is comparable with previous lattice and data-driven determinations, with combined relative uncertainties below $10\%$. For the $\eta$ meson, we report on a preliminary determination from a single lattice spacing.

We present new lattice results of the ETM Collaboration, obtained from extensive simulations of lattice QCD with dynamical up, down, strange and charm quarks at physical mass values, different volumes and lattice spacings, concerning the SM prediction for the so-called intermediate window (W) and short-distance (SD) contributions to the leading order hadronic vacuum polarization (LO-HVP) term of the muon anomalous magnetic moment, $a_\mu$. Results for $a_{\mu}^{\rm LO-HVP,W}$ and $a_{\mu}^{\rm LO-HVP,SD}$, besides representing a step forward to a complete lattice computation of $a_{\mu}^{\rm LO-HVP}$ and a useful benchmark among lattice groups, are compared here with their dispersive counterparts based on experimental data for $e^+e^-$ into hadrons. The comparison confirms the tension in $a_{\mu}^{\rm LO-HVP,W}$, already noted in 2020 by the BMW Collaboration, while showing no tension in $a_{\mu}^{\rm LO-HVP,SD}$., Comment: Talk given at ICHEP 2022 (6-13 July 2022, Bologna - Italy) - Results here are almost final - Contribution accepted for publication on PoS

We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_{\mu}^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_{\mu}^{\rm SD}$ and $a_{\mu}^{\rm W}$. We employ gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavours of Wilson-clover twisted-mass quarks with masses of all the dynamical quark flavours tuned close to their physical values. The simulations are carried out at three values of the lattice spacing equal to $\simeq 0.057, 0.068$ and $0.080$ fm with spatial lattice sizes up to $L \simeq 7.6$~fm. For the short distance window we obtain $a_\mu^{\rm SD} = 69.27\,(34) \cdot 10^{-10}$, in agreement with the dispersive determination based on experimental $e^+ e^-$ data. For the intermediate window we get instead $a_\mu^{\rm W} = 236.3\,(1.3) \cdot 10^{-10}$, which is consistent with recent determinations by other lattice collaborations, but disagrees with the dispersive determination at the level of $3.6\,\sigma$., Comment: 10 pages, 4 figures, 1 table, LATTICE 2022

We present a first-principles lattice QCD investigation of the $R$-ratio between the $e^+e^-$ cross-section into hadrons and that into muons. By using the method of Ref.[1], that allows to extract smeared spectral densities from Euclidean correlators, we compute the $R$-ratio convoluted with Gaussian smearing kernels of widths of about $600$ MeV and central energies from $220$ MeV up to $2.5$ GeV. Our theoretical results are compared with the corresponding quantities obtained by smearing the KNT19 compilation [2] of $R$-ratio experimental measurements with the same kernels and, by centring the Gaussians in the region around the $\rho$-resonance peak, a tension of about three standard deviations is observed. From the phenomenological perspective, we have not included yet in our calculation QED and strong isospin-breaking corrections and this might affect the observed tension. From the methodological perspective, our calculation demonstrates that it is possible to study the $R$-ratio in Gaussian energy bins on the lattice at the level of accuracy required in order to perform precision tests of the Standard Model., Comment: Version accepted for publication on PRL. Results unchanged

We present preliminary results for the renormalization functions (RFs) of a number of quark and gluon operators studied in lattice QCD using a gauge-invariant renormalization scheme (GIRS). GIRS is a variant of the coordinate-space renormalization prescription, in which Green's functions of gauge-invariant operators are calculated in position space. A novel aspect is that summations over different time slices of the positions of the operators are employed in order to reduce the statistical noise in lattice simulations. We test the reliability of this scheme by calculating RFs for the vector one-derivative quark bilinear operator, which enters the average momentum fraction of the nucleon. We use $N_f=4$ degenerate twisted mass clover-improved fermion ensembles of different volumes and lattice spacings. We also present first results of applying GIRS when operator mixing occurs: the mixing coefficients of the gluon and quark singlet energy-momentum tensor operators are evaluated by imposing appropriate renormalization conditions on the lattice., Comment: 10 pages, 5 figures, 3 tables, Talk at the 39th International Symposium on Lattice Field Theory, LATTICE2022, 8th-13th August, 2022, Rheinische Friedrich-Wilhelms-Universit\"at Bonn, Bonn, Germany

The Majda-Biello system models the interaction of Rossby waves. It consists of two coupled KdV equations one of which has a parameter $\alpha$ as coefficient of its dispersion. This work studies this system on the half line with Robin, Neumann, and Dirichlet boundary data. It shows that for $0<\alpha<1$ or $1<\alpha<4$ all these problems are well-posed for initial data in Sobolev spaces $H^s$, $s\ge 0$. For $\alpha=1$ or $\alpha>4$ well-posedness holds for Dirichlet data if $s>-3/4$, while for Neumann and Robin data it depends on the sign of the parameters involved in the data. For $\alpha=4$ well-posedness of all problems holds for $s\ge 3/4$. The Robin and Neumann boundary data are in $H^{s/3}$ while the Dirichlet boundary data are in $H^{(s+1)/3}$. These are consistent with the time regularity of the Cauchy problem for the corresponding linear system. The proof is based on linear estimates in Bourgain spaces derived by utilizing the Fokas solution formula for the forced linear system, and appropriate bilinear estimates suggested by the coupled nonlinearities. These show that the iteration map defined via the Fokas formula is a contraction in appropriate solution spaces. All the well-posedness results obtained here are optimal.

We present a lattice determination of the disconnected contributions to the leading-order hadronic vacuum polarization (HVP) to the muon anomalous magnetic moment in the so-called short and intermediate time-distance windows. We employ gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavours of Wilson twisted-mass clover-improved quarks with masses approximately tuned to their physical value. We take the continuum limit employing three lattice spacings at about 0.08, 0.07 and 0.06 fm., Comment: 10 pages, 4 figures, 7 tables, LATTICE 2022

We calculate the double-virtual $\eta \rightarrow \gamma^* \gamma^*$ transition form factor $\mathcal{F}_{\eta \to \gamma^* \gamma^*}(q_1^2,q_2^2)$ from first principles using a lattice QCD simulation with $N_f=2+1+1$ quark flavors at the physical pion mass and at one lattice spacing and volume. The kinematic range covered by our calculation is complementary to the one accessible from experiment and is relevant for the $\eta$-pole contribution to the hadronic light-by-light scattering in the anomalous magnetic moment $a_\mu = (g-2)/2$ of the muon. From the form factor calculation we extract the partial decay width $\Gamma(\eta \rightarrow \gamma \gamma) = 323(85)_\text{stat}(22)_\text{syst}$ eV and the slope parameter $b_\eta=1.19(36)_\text{stat}(16)_\text{syst}$ GeV${}^{-2}$. For the $\eta$-pole contribution to $a_\mu$ we obtain $a_\mu^{\eta-\text{pole}} = 13.2(5.2)_\text{stat}(1.3)_\text{syst} \cdot 10^{-11}$., Comment: 12 pages, 8 figures; updated to published version

We present our recent results on antiheavy-antiheavy-light-light tetraquark systems using lattice QCD. Our study of the $ \bar{b}\bar{b}us $ four-quark system with quantum numbers $ J^P=1^+ $ and the $ \bar{b}\bar{c}ud $ four-quark systems with $ I(J^P)=0(0^+) $ and $ I(J^P)=0(1^+) $ utilizes scattering operators at the sink to improve the extraction of the low-lying energy levels. We found a bound state for $ \bar{b}\bar{b}us $ with $ E_{\textrm{bind},\bar{b}\bar{b}us} = (-86 \pm 22 \pm 10)\,\textrm{MeV} $, but no indication for a bound state in both $ \bar{b}\bar{c}ud $ channels. Moreover, we show preliminary results for $ \bar{b}\bar{b}ud $ with $ I(J^P)=0(1^+) $, where we used scattering operators both at the sink and the source. We found a bound state and determined its infinite-volume binding energy with a scattering analysis, resulting in $ E_{\textrm{bind},\bar{b}\bar{b}ud} =(-103 \pm 8 )\,\textrm{MeV} $., Comment: 9 pages, 3 figures, talk given at "The 39th International Symposium on Lattice Field Theory", 08th-13th August 2022, Bonn, Germany

We present first results of a recently started lattice QCD investigation of antiheavy-antiheavy-light-light tetraquark systems including scattering interpolating operators in correlation functions both at the source and at the sink. In particular, we discuss the importance of such scattering interpolating operators for a precise computation of the low-lying energy levels. We focus on the $\bar b \bar b u d$ four-quark system with quantum numbers $I(J^P) = 0(1^+)$, which has a ground state below the lowest meson-meson threshold. We carry out a scattering analysis using L\"uscher's method to extrapolate the binding energy of the corresponding QCD-stable tetraquark to infinite spatial volume. Our calculation uses clover $u$, $d$ valence quarks and NRQCD $b$ valence quarks on gauge-link ensembles with HISQ sea quarks that were generated by the MILC collaboration., Comment: 6 pages, 2 figures, talk given at "The 39th International Symposium on Lattice Field Theory", 08th-13th August 2022, Bonn, Germany