Showing 1–50 of 175 results for author: Urbach, C

The $ω$-meson from lattice QCD

Authors: Haobo Yan, Marco Garofalo, Maxim Mai, Ulf-G. Meißner, Carsten Urbach

Abstract: Many excited states in the hadron spectrum have large branching ratios to the three-hadron final states. Understanding such particles from first principles QCD requires input from lattice QCD with one-, two-, and three-meson interpolators as well as a reliable three-body formalism relating finite-volume spectra at unphysical pion masses to the scattering amplitudes at the physical point. In this w… ▽ More Many excited states in the hadron spectrum have large branching ratios to the three-hadron final states. Understanding such particles from first principles QCD requires input from lattice QCD with one-, two-, and three-meson interpolators as well as a reliable three-body formalism relating finite-volume spectra at unphysical pion masses to the scattering amplitudes at the physical point. In this work, we provide the first-ever calculation of the resonance parameters of the $ω$-meson from lattice QCD to the physical point, including an update of the formalism through matching to effective field theories. The main result of this pioneering study, the pole position of the $ω$-meson at $\sqrt{s_ω}= (778.0(11.2)-i\,3.0(5) )\,\mathrm{MeV}$, agrees reasonably well with experiment providing a pathway to study heavier three-hadron resonances. In addition we provide an estimate of the $ω-ρ$ mass difference as $29(15)\,\mathrm{MeV}$. △ Less

Submitted 23 July, 2024; originally announced July 2024.

Comments: 13 pages, 10 figures, 3 tables

Quark and gluon momentum fractions in the pion and in the kaon

Authors: Constantia Alexandrou, Simone Bacchio, Martha Constantinou, Joseph Delmar, Jacob Finkenrath, Bartosz Kostrzewa, Marcus Petschlies, Luis Alberto Rodriguez Chacon, Gregoris Spanoudes, Fernanda Steffens, Carsten Urbach, Urs wenger

Abstract: 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 de… ▽ More 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. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 14 pages, 10 figures

Towards determining the (2+1)-dimensional Quantum Electrodynamics running coupling with Monte Carlo and quantum computing methods

Authors: Arianna Crippa, Simone Romiti, Lena Funcke, Karl Jansen, Stefan Kühn, Paolo Stornati, Carsten Urbach

Abstract: In this paper, we examine a compact $U(1)$ lattice gauge theory in $(2+1)$ dimensions and present a strategy for studying the running coupling and extracting the non-perturbative $Λ$-parameter. To this end, we combine Monte Carlo simulations and quantum computing, where the former can be used to determine the numerical value of the lattice spacing $a$, and the latter allows for reaching the pertur… ▽ More In this paper, we examine a compact $U(1)$ lattice gauge theory in $(2+1)$ dimensions and present a strategy for studying the running coupling and extracting the non-perturbative $Λ$-parameter. To this end, we combine Monte Carlo simulations and quantum computing, where the former can be used to determine the numerical value of the lattice spacing $a$, and the latter allows for reaching the perturbative regime at very small values of the bare coupling and, correspondingly, small values of $a$. The methodology involves a series of sequential steps (i.e., the step scaling function) to bridge results from small lattice spacings to non-perturbative large-scale lattice calculations. Focusing on the pure gauge case, we demonstrate that these quantum circuits, adapted to gauge degrees of freedom, are able to capture the relevant physics by studying the expectation value of the plaquette operator, for matching with corresponding Monte Carlo simulations. We also present results for the static potential and static force, which can be related to the renormalized coupling. The procedure outlined in this work can be extended to Abelian and non-Abelian lattice gauge theories with matter fields and might provide a way towards studying lattice quantum chromodynamics utilizing both quantum and classical methods. △ Less

Submitted 12 June, 2024; v1 submitted 26 April, 2024; originally announced April 2024.

Comments: Edit and add references. Fix typos. Replace Figure 7. Add acknowledgements

Inclusive hadronic decay rate of the $τ$ lepton from lattice QCD: the $\bar u s$ flavour channel and the Cabibbo angle

Authors: C. Alexandrou, S. Bacchio, A. De Santis, A. Evangelista, J. Finkenrath, R. Frezzotti, G. Gagliardi, M. Garofalo, B. Kostrzewa, V. Lubicz, S. Romiti, F. Sanfilippo, S. Simula, N. Tantalo, C. Urbach, U. Wenger

Abstract: We present a lattice determination of the inclusive decay rate of the process $τ\mapsto X_{us} ν_τ$ in which the $τ$ 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 isospi… ▽ More We present a lattice determination of the inclusive decay rate of the process $τ\mapsto X_{us} ν_τ$ in which the $τ$ 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 $τ\mapsto X_{ud} ν_τ$ 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. △ Less

Submitted 10 June, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Comments: 9 pages, 5 figures, 1 table. Version accepted by PRL, expanded technical discussion moved to an appendix, results unchanged

Testing the $\mathrm{SU}(2)$ lattice Hamiltonian built from $S_3$ partitionings

Authors: Marco Garofalo, Tobias Hartung, Timo Jakobs, Karl Jansen, Johann Ostmeyer, Dominik Rolfes, Simone Romiti, Carsten Urbach

Abstract: We test a possible digitization of $\mathrm{SU}(2)$ lattice gauge theories based on partitionings of the sphere $S_3$. In our construction the link operators are unitary and diagonal, with eigenvalues determined by the vertices of the partitioning. The canonical momenta are finite difference operators approximating the Lie derivatives on the manifold. In this formalism we implement the standard Wi… ▽ More We test a possible digitization of $\mathrm{SU}(2)$ lattice gauge theories based on partitionings of the sphere $S_3$. In our construction the link operators are unitary and diagonal, with eigenvalues determined by the vertices of the partitioning. The canonical momenta are finite difference operators approximating the Lie derivatives on the manifold. In this formalism we implement the standard Wilson Hamiltonian. We show results for a 2-site Schwinger-type model in 1D and a single-plaquette system in 2D. Our calculations are performed on a classical computer, though in principle they can be implemented also on a quantum device. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Digitizing lattice gauge theories in the magnetic basis: reducing the breaking of the fundamental commutation relations

Authors: Simone Romiti, Carsten Urbach

Abstract: We present a digitization scheme for the lattice $\mathrm{SU}(2)$ gauge theory Hamiltonian in the $\mathit{magnetic}$ $\mathit{basis}$, where the gauge links are unitary and diagonal. The digitization is obtained from a particular partitioning of the $\mathrm{SU}(2)$ group manifold, with the canonical momenta constructed by an approximation of the Lie derivatives on this partitioning. This constru… ▽ More We present a digitization scheme for the lattice $\mathrm{SU}(2)$ gauge theory Hamiltonian in the $\mathit{magnetic}$ $\mathit{basis}$, where the gauge links are unitary and diagonal. The digitization is obtained from a particular partitioning of the $\mathrm{SU}(2)$ group manifold, with the canonical momenta constructed by an approximation of the Lie derivatives on this partitioning. This construction, analogous to a discrete Fourier transform, preserves the spectrum of the kinetic part of the Hamiltonian and the canonical commutation relations exactly on a subspace of the truncated Hilbert space, while the residual subspace can be projected above the cutoff of the theory. △ Less

Submitted 17 July, 2024; v1 submitted 20 November, 2023; originally announced November 2023.

Nucleon axial and pseudoscalar form factors using twisted-mass fermion ensembles at the physical point

Authors: Constantia Alexandrou, Simone Bacchio, Martha Constantinou, Jacob Finkenrath, Roberto Frezzotti, Bartosz Kostrzewa, Giannis Koutsou, Gregoris Spanoudes, Carsten Urbach

Abstract: 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_πL$>3.6. Convergence to the gro… ▽ More 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_π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_{πNN} \equiv \lim_{Q^2 \rightarrow -m_π^2} G_{πNN}(Q^2)$=13.25(67)(69) and for $G_P(0.88m_μ^2)\equiv g_P^*$=8.99(39)(49). △ Less

Submitted 9 January, 2024; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: accepted version, 32 pages, 39 figures, 14 tables

Pion Transition Form Factor from Twisted-Mass Lattice QCD and the Hadronic Light-by-Light $π^0$-pole Contribution to the Muon $g-2$

Authors: C. Alexandrou, S. Bacchio, G. Bergner, S. Burri, J. Finkenrath, A. Gasbarro, K. Hadjiyiannakou, K. Jansen, G. Kanwar, B. Kostrzewa, G. Koutsou, K. Ottnad, M. Petschlies, F. Pittler, F. Steffens, C. Urbach, U. Wenger

Abstract: 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}_{π^0γγ}(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 para… ▽ More 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}_{π^0γγ}(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_μ^{π^0\text{-pole}} = 56.7(3.2) \times 10^{-11}$. △ Less

Submitted 3 January, 2024; v1 submitted 23 August, 2023; originally announced August 2023.

Comments: 21 pages, 17 figures, 4 tables, updated to published version

Journal ref: Phys.Rev.D 108 (2023) 9, 094514

Exploring a new approach to Hadronic Parity Violation from Lattice QCD

Authors: Marcus Petschlies, Nikolas Schlage, Aniket Sen, Carsten Urbach

Abstract: The long-range, parity-odd nucleon interaction generated by single pion exchange is captured in the parity-odd pion-nucleon coupling $h^1_π$. Its calculation in lattice QCD requires the evaluation of 4-quark operator nucleon 3-point functions. We investigate a new numerical approach to compute $h^1_π$ based on nucleon matrix elements of parity-even 4-quark operators and related to the parity-viola… ▽ More The long-range, parity-odd nucleon interaction generated by single pion exchange is captured in the parity-odd pion-nucleon coupling $h^1_π$. Its calculation in lattice QCD requires the evaluation of 4-quark operator nucleon 3-point functions. We investigate a new numerical approach to compute $h^1_π$ based on nucleon matrix elements of parity-even 4-quark operators and related to the parity-violating electro-weak theory by PCAC and chiral perturbation theory. This study is performed with 2+1+1 dynamical flavors of twisted mass fermions at pion mass $m_π \approx 260\,\text{MeV}$ in a lattice box of $L \approx 3 \,\text{fm} $ and with a lattice spacing of $a \approx 0.091 \,\text{fm}$. From a calculation excluding fermion loop diagrams we find a bare coupling of $h^1_π= 8.08 \,(98) \cdot 10^{-7}$. △ Less

Submitted 28 November, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: 14 pages, 9 figures, 4 tables. Parity-violating Lagrangian Eq. (1), PCAC relation Eq. (2), details on 4-quark operators Eq. (4), summary of statistics Tab. I., summary of our results for the bare matrix elements weighted by the Wilson coefficients Tab. IV. and references to our preliminary results of this work added. Results of the paper are unaffected

Canonical Momenta in Digitized SU(2) Lattice Gauge Theory: Definition and Free Theory

Authors: Timo Jakobs, Marco Garofalo, Tobias Hartung, Karl Jansen, Johann Ostmeyer, Dominik Rolfes, Simone Romiti, Carsten Urbach

Abstract: Hamiltonian simulations of quantum systems require a finite-dimensional representation of the operators acting on the Hilbert space H. Here we give a prescription for gauge links and canonical momenta of an SU(2) gauge theory, such that the matrix representation of the former is diagonal in H. This is achieved by discretising the sphere $S_3$ isomorphic to SU(2) and the corresponding directional d… ▽ More Hamiltonian simulations of quantum systems require a finite-dimensional representation of the operators acting on the Hilbert space H. Here we give a prescription for gauge links and canonical momenta of an SU(2) gauge theory, such that the matrix representation of the former is diagonal in H. This is achieved by discretising the sphere $S_3$ isomorphic to SU(2) and the corresponding directional derivatives. We show that the fundamental commutation relations are fulfilled up to discretisation artefacts. Moreover, we directly construct the Casimir operator corresponding to the Laplace-Beltrami operator on $S_3$ and show that the spectrum of the free theory is reproduced again up to discretisation effects. Qualitatively, these results do not depend on the specific discretisation of SU(2), but the actual convergence rates do. △ Less

Submitted 28 July, 2023; v1 submitted 5 April, 2023; originally announced April 2023.

Comments: typos corrected, matches published version in EPJC

Hadronic Parity Violation from Twisted Mass Lattice QCD

Authors: Nikolas Schlage, Marcus Petschlies, Aniket Sen, Carsten Urbach

Abstract: We present results for an exploratory lattice calculation of the leading parity-violating pion-nucleon coupling $h_π^1$. Based on the PCAC relation we use a parity-conserving Lagrangian and focus on the techniques to determine the nucleon matrix elements of the effective 4-quark operators. For our study we employ an ensemble of twisted mass fermions with $260$ MeV pion mass. Barring mixing with lo… ▽ More We present results for an exploratory lattice calculation of the leading parity-violating pion-nucleon coupling $h_π^1$. Based on the PCAC relation we use a parity-conserving Lagrangian and focus on the techniques to determine the nucleon matrix elements of the effective 4-quark operators. For our study we employ an ensemble of twisted mass fermions with $260$ MeV pion mass. Barring mixing with lower-dimensional operators and renormalization at this stage, we discuss our estimate for $h_π^1$. △ Less

Submitted 19 May, 2023; v1 submitted 16 January, 2023; originally announced January 2023.

Comments: LATTICE2022 proceedings, revised coupling $h_π^1$ in Tabs. 3 and 4 and in corresponding Fig. 5

Lattice calculation of the R-ratio smeared with Gaussian kernel

Authors: Constantia Alexandrou, Simone Bacchio, Alessandro De Santis, Petros Dimopoulos, Jacob Finkenrath, Roberto Frezzotti, Giuseppe Gagliardi, Marco Garofalo, Kyriakos Hadjiyiannakou, Bartosz Kostrzewa, Karl Jansen, Vittorio Lubicz, Marcus Petschlies, Francesco Sanfilippo, Silvano Simula, Nazario Tantalo, Carsten Urbach, Urs Wenger

Abstract: The ratio $R(E)$ of the cross-sections for $e^+e^-\to$ hadrons and $e^+e^-\to μ^+μ^-$ 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 cal… ▽ More The ratio $R(E)$ of the cross-sections for $e^+e^-\to$ hadrons and $e^+e^-\to μ^+μ^-$ 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 $σ$ 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 $σ\sim 600$ MeV and central energies around the $ρ$-resonance peak. △ Less

Submitted 23 December, 2022; originally announced December 2022.

Comments: Contribution to the 39th International Symposium on Lattice Field Theory (Lattice2022), 8th-13th August, 2022, Bonn, Germany

Short \& intermediate distance HVP contributions to muon g-2: SM (lattice) prediction versus $e^+e^-$ annihilation data

Authors: Constantia Alexandrou, Simone Bacchio, Petros Dimopoulos, Jacob Finkenrath, Roberto Frezzotti, Giuseppe Gagliardi, Marco Garofalo, Kyriakos Hadjiyiannakou, Bartosz Kostrzewa, Karl Jansen, Vittorio Lubicz, Marcus Petschlies, Francesco Sanfilippo, Silvano Simula, Carsten Urbach, Urs Wenger

Abstract: 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 t… ▽ More 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_μ$. Results for $a_μ^{\rm LO-HVP,W}$ and $a_μ^{\rm LO-HVP,SD}$, besides representing a step forward to a complete lattice computation of $a_μ^{\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_μ^{\rm LO-HVP,W}$, already noted in 2020 by the BMW Collaboration, while showing no tension in $a_μ^{\rm LO-HVP,SD}$. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: Talk given at ICHEP 2022 (6-13 July 2022, Bologna - Italy) - Results here are almost final - Contribution accepted for publication on PoS

Pseudoscalar-pole contributions to the muon $g-2$ at the physical point

Authors: Sebastian Burri, Gurtej Kanwar, Constantia Alexandrou, Simone Bacchio, Georg Bergner, Jacob Finkenrath, Andrew Gasbarro, Kyriakos Hadjiyiannakou, Karl Jansen, Bartosz Kostrzewa, Giannis Koutsou, Konstantin Ottnad, Marcus Petschlies, Ferenc Pittler, Fernanda Steffens, Carsten Urbach, Urs Wenger

Abstract: 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 γ^* γ^*}$ for the neutral pseudoscalar mesons $P=π^0$ and $η$. The calculation is performed using twisted-mass lattice QCD with physical quark masses. On the lattice, we have… ▽ More 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 γ^* γ^*}$ for the neutral pseudoscalar mesons $P=π^0$ and $η$. 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 $η$-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 $η$ meson, we report on a preliminary determination from a single lattice spacing. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Hamiltonian limit of lattice QED in 2+1 dimensions

Authors: L. Funcke, C. F. Groß, K. Jansen, S. Kühn, S. Romiti, C. Urbach

Abstract: The Hamiltonian limit of lattice gauge theories can be found by extrapolating the results of anisotropic lattice computations, i.e., computations using lattice actions with different temporal and spatial lattice spacings ($a_t\neq a_s$), to the limit of $a_t\to 0$. In this work, we present a study of this Hamiltonian limit for a Euclidean $U(1)$ gauge theory in 2+1 dimensions (QED3), regularized o… ▽ More The Hamiltonian limit of lattice gauge theories can be found by extrapolating the results of anisotropic lattice computations, i.e., computations using lattice actions with different temporal and spatial lattice spacings ($a_t\neq a_s$), to the limit of $a_t\to 0$. In this work, we present a study of this Hamiltonian limit for a Euclidean $U(1)$ gauge theory in 2+1 dimensions (QED3), regularized on a toroidal lattice. The limit is found using the renormalized anisotropy $ξ_R=a_t/a_s$, by sending $ξ_R \to 0$ while keeping the spatial lattice spacing constant. We compute $ξ_R$ in $3$ different ways: using both the ``normal'' and the ``sideways'' static quark potential, as well as the gradient flow evolution of gauge fields. The latter approach will be particularly relevant for future investigations of combining quantum computations with classical Monte Carlo computations, which requires the matching of lattice results obtained in the Hamiltonian and Lagrangian formalisms. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Report number: MIT-CTP/5480

Digitizing SU(2) Gauge Fields and What to Look Out for When Doing So

Authors: Tobias Hartung, Timo Jakobs, Karl Jansen, Johann Ostmeyer, Carsten Urbach

Abstract: With the long term perspective of using quantum computers and tensor networks for lattice gauge theory simulations, an efficient method of digitizing gauge group elements is needed. We thus present our results for a handful of discretization approaches for the non-trivial example of SU(2), such as its finite subgroups, as well as different classes of finite subsets. We focus our attention on a fre… ▽ More With the long term perspective of using quantum computers and tensor networks for lattice gauge theory simulations, an efficient method of digitizing gauge group elements is needed. We thus present our results for a handful of discretization approaches for the non-trivial example of SU(2), such as its finite subgroups, as well as different classes of finite subsets. We focus our attention on a freezing transition observed towards weak couplings. A generalized version of the Fibonacci spiral appears to be particularly efficient and close to optimal. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: Proceedings of the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, Rheinische Friedrich-Wilhelms-Universität Bonn, Bonn, Germany

Time windows of the muon HVP from twisted-mass lattice QCD

Authors: C. Alexandrou, S. Bacchio, P. Dimopoulos, J. Finkenrath, R. Frezzotti, G. Gagliardi, M. Garofalo, K. Hadjiyiannakou, B. Kostrzewa, K. Jansen, V. Lubicz, M. Petschlies, F. Sanfilippo, S. Simula, C. Urbach, U. Wenger

Abstract: We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_μ^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_μ^{\rm SD}$ and $a_μ^{\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-m… ▽ More We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_μ^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_μ^{\rm SD}$ and $a_μ^{\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_μ^{\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_μ^{\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\,σ$. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 10 pages, 4 figures, 1 table, LATTICE 2022

Probing the energy-smeared R-ratio on the lattice

Authors: Constantia Alexandrou, Simone Bacchio, Alessandro De Santis, Petros Dimopoulos, Jacob Finkenrath, Roberto Frezzotti, Giuseppe Gagliardi, Marco Garofalo, Kyriakos Hadjiyiannakou, Bartosz Kostrzewa, Karl Jansen, Vittorio Lubicz, Marcus Petschlies, Francesco Sanfilippo, Silvano Simula, Nazario Tantalo, Carsten Urbach, Urs Wenger

Abstract: 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.… ▽ More 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 $ρ$-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. △ Less

Submitted 19 May, 2023; v1 submitted 16 December, 2022; originally announced December 2022.

Comments: Version accepted for publication on PRL. Results unchanged

Disconnected contribution to the LO HVP term of muon g-2 from ETMC

Authors: Constantia Alexandrou, Simone Bacchio, Petros Dimopoulos, Jacob Finkenrath, Roberto Frezzotti, Giuseppe Gagliardi, Marco Garofalo, Kyriakos Hadjiyiannakou, Bartosz Kostrzewa, Karl Jansen, Vittorio Lubicz, Marcus Petschlies, Francesco Sanfilippo, Silvano Simula, Carsten Urbach, Urs Wenger

Abstract: 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… ▽ More 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. △ Less

Submitted 14 December, 2022; originally announced December 2022.

Comments: 10 pages, 4 figures, 7 tables, LATTICE 2022

The $η\rightarrow γ^* γ^*$ transition form factor and the hadronic light-by-light $η$-pole contribution to the muon $g-2$ from lattice QCD

Authors: Constantia Alexandrou, Simone Bacchio, Sebastian Burri, Jacob Finkenrath, Andrew Gasbarro, Kyriakos Hadjiyiannakou, Karl Jansen, Gurtej Kanwar, Bartosz Kostrzewa, Konstantin Ottnad, Marcus Petschlies, Ferenc Pittler, Carsten Urbach, Urs Wenger

Abstract: We calculate the double-virtual $η\rightarrow γ^* γ^*$ transition form factor $\mathcal{F}_{η\to γ^* γ^*}(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 $η$-p… ▽ More We calculate the double-virtual $η\rightarrow γ^* γ^*$ transition form factor $\mathcal{F}_{η\to γ^* γ^*}(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 $η$-pole contribution to the hadronic light-by-light scattering in the anomalous magnetic moment $a_μ= (g-2)/2$ of the muon. From the form factor calculation we extract the partial decay width $Γ(η\rightarrow γγ) = 323(85)_\text{stat}(22)_\text{syst}$ eV and the slope parameter $b_η=1.19(36)_\text{stat}(16)_\text{syst}$ GeV${}^{-2}$. For the $η$-pole contribution to $a_μ$ we obtain $a_μ^{η-\text{pole}} = 13.2(5.2)_\text{stat}(1.3)_\text{syst} \cdot 10^{-11}$. △ Less

Submitted 21 December, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

Comments: 12 pages, 8 figures; updated to published version

Journal ref: Phys. Rev. D 108 (2023) 054509

Twisted mass ensemble generation on GPU machines

Authors: Bartosz Kostrzewa, Simone Bacchio, Jacob Finkenrath, Marco Garofalo, Ferenc Pittler, Simone Romiti, Carsten Urbach

Abstract: We present how we ported the Hybrid Monte Carlo implementation in the tmLQCD software suite to GPUs through offloading its most expensive parts to the QUDA library. We discuss our motivations and some of the technical challenges that we encountered as we added the required functionality to both tmLQCD and QUDA. We further present some performance details, focussing in particular on the usage of QU… ▽ More We present how we ported the Hybrid Monte Carlo implementation in the tmLQCD software suite to GPUs through offloading its most expensive parts to the QUDA library. We discuss our motivations and some of the technical challenges that we encountered as we added the required functionality to both tmLQCD and QUDA. We further present some performance details, focussing in particular on the usage of QUDA's multigrid solver for poorly conditioned light quark monomials as well as the multi-shift solver for the non-degenerate strange and charm sector in $N_f=2+1+1$ simulations using twisted mass clover fermions, comparing the efficiency of state-of-the-art simulations on CPU and GPU machines. We also take a look at the performance-portability question through preliminary tests of our HMC on a machine based on AMD's MI250 GPU, finding good performance after a very minor additional porting effort. Finally, we conclude that we should be able to achieve GPU utilisation factors acceptable for the current generation of (pre-)exascale supercomputers with subtantial efficiency improvements and real time speedups compared to just running on CPUs. At the same time, we find that future challenges will require different approaches and, most importantly, a very significant investment of personnel for software development. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: 10 pages, 7 figures, Contribution to the 39th International Symposium on Lattice Field Theory (Lattice2022), 8-13 August, 2022, Bonn, Germany

Three-body resonances in the $\varphi^4$ theory

Authors: Marco Garofalo, Maxim Mai, Fernando Romero-López, Akaki Rusetsky, Carsten Urbach

Abstract: We study the properties of three-body resonances using a lattice complex scalar $\varphi^4$ theory with two scalars, with parameters chosen such that one heavy particle can decay into three light ones. We determine the two- and three-body spectra for several lattice volumes using variational techniques, and then analyze them with two versions of the three-particle finite-volume formalism: the Rela… ▽ More We study the properties of three-body resonances using a lattice complex scalar $\varphi^4$ theory with two scalars, with parameters chosen such that one heavy particle can decay into three light ones. We determine the two- and three-body spectra for several lattice volumes using variational techniques, and then analyze them with two versions of the three-particle finite-volume formalism: the Relativistic Field Theory approach and the Finite-Volume Unitarity approach. We find that both methods provide an equivalent description of the energy levels, and we are able to fit the spectra using simple parametrizations of the scattering quantities. By solving the integral equations of the corresponding three-particle formalisms, we determine the pole position of the resonance in the complex energy-plane and thereby its mass and width. We find very good agreement between the two methods at different values of the coupling of the theory. △ Less

Submitted 7 March, 2023; v1 submitted 10 November, 2022; originally announced November 2022.

Comments: 26 pages, 6 figures, 4 tables

Report number: MIT-CTP/5487

Journal ref: JHEP 02 (2023) 252

Defining Canonical Momenta for Discretised SU(2) Gauge Fields

Authors: Marco Garofalo, Tobias Hartung, Karl Jansen, Johann Ostmeyer, Simone Romiti, Carsten Urbach

Abstract: In this proceeding contribution we discuss how to define canonical momenta for SU(N) lattice gauge theories in the Hamiltonian formalism in a basis where the gauge field operators are diagonal. For an explicit discretisation of SU(2) we construct the momenta and check the violation of the fundamental commutation relations. In this proceeding contribution we discuss how to define canonical momenta for SU(N) lattice gauge theories in the Hamiltonian formalism in a basis where the gauge field operators are diagonal. For an explicit discretisation of SU(2) we construct the momenta and check the violation of the fundamental commutation relations. △ Less

Submitted 27 October, 2022; originally announced October 2022.

Comments: Proceedings of the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, Rheinische Friedrich-Wilhelms-Universität Bonn, Bonn, Germany

Nucleon transverse quark spin densities

Authors: Constantia Alexandrou, Simone Bacchio, Martha Constantinou, Petros Dimopoulos, Jacob Finkenrath, Roberto Frezzotti, Kyriakos Hadjiyiannakou, Karl Jansen, Bartosz Kostrzewa, Giannis Koutsou, Gregoris Spanoudes, Carsten Urbach

Abstract: We present a calculation of the Mellin moments of the nucleon transverse quark spin densities extracted from the unpolarized and transversity generalized form factors. We use three $N_F=2+1+1$ ensembles of twisted mass fermions with quark masses tuned to their physical values and lattice spacings $a\sim 0.08$~fm, $a\sim 0.07$~fm and $a\sim 0.06$~fm and extrapolate the form factors to the continuum… ▽ More We present a calculation of the Mellin moments of the nucleon transverse quark spin densities extracted from the unpolarized and transversity generalized form factors. We use three $N_F=2+1+1$ ensembles of twisted mass fermions with quark masses tuned to their physical values and lattice spacings $a\sim 0.08$~fm, $a\sim 0.07$~fm and $a\sim 0.06$~fm and extrapolate the form factors to the continuum limit. Besides isovector densities we also include results for the tensor charge for each quark flavor using the ensemble with $a\sim 0.08$~fm for which we include the disconnected contributions. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: 9 pages, 6 figures, talk presented at the 39th International Symposium on Lattice Field Theory, LATTICE2022 8th-13th August, 2022, University of Bonn, Germany

Lattice calculation of the short and intermediate time-distance hadronic vacuum polarization contributions to the muon magnetic moment using twisted-mass fermions

Authors: C. Alexandrou, S. Bacchio, P. Dimopoulos, J. Finkenrath, R. Frezzotti, G. Gagliardi, M. Garofalo, K. Hadjiyiannakou, B. Kostrzewa, K. Jansen, V. Lubicz, M. Petschlies, F. Sanfilippo, S. Simula, C. Urbach, U. Wenger

Abstract: We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_μ^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_μ^{\rm SD}$ and $a_μ^{\rm W}$, defined by the RBC/UKQCD Collaboration [1]. We employ gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with… ▽ More We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_μ^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_μ^{\rm SD}$ and $a_μ^{\rm W}$, defined by the RBC/UKQCD Collaboration [1]. We employ gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks with masses of all the dynamical quark flavors 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_μ^{\rm SD}({\rm ETMC}) = 69.27\,(34) \cdot 10^{-10}$, which is consistent with the recent dispersive value of $a_μ^{\rm SD}(e^+ e^-) = 68.4\,(5) \cdot 10^{-10}$ [2]. In the case of the intermediate window we get the value $a_μ^{\rm W}({\rm ETMC}) = 236.3\,(1.3) \cdot 10^{-10}$, which is consistent with the result $a_μ^{\rm W}({\rm BMW}) = 236.7\,(1.4) \cdot 10^{-10}$ [3] by the BMW collaboration as well as with the recent determination by the CLS/Mainz group of $a_μ^{\rm W}({\rm CLS}) = 237.30\,(1.46) \cdot 10^{-10}$ [4]. However, it is larger than the dispersive result of $a_μ^{\rm W}(e^+ e^-) = 229.4\,(1.4) \cdot 10^{-10}$ [2] by approximately $3.6$ standard deviations. The tension increases to approximately $4.5$ standard deviations if we average our ETMC result with those by BMW and CLS/Mainz. Our accurate lattice results in the short and intermediate windows point to a possible deviation of the $e^+ e^-$ cross section data with respect to Standard Model predictions in the low and intermediate energy regions, but not in the high energy region. △ Less

Submitted 29 November, 2022; v1 submitted 30 June, 2022; originally announced June 2022.

Comments: 78 pages, 22 figures, 14 tables. Analysis improved with more data and fits, presentation reorganized, more material in appendices, fixed typos

Towards a theory of hadron resonances

Authors: Maxim Mai, Ulf-G. Meißner, Carsten Urbach

Abstract: In this review, we present the current state of the art of our understanding of the spectrum of excited strongly interacting particles and discuss methods that allow for a systematic and model-independent calculation of the hadron spectrum. These are lattice QCD and effective field theories. Synergies between both approaches can be exploited allowing for deeper understanding of the hadron spectrum… ▽ More In this review, we present the current state of the art of our understanding of the spectrum of excited strongly interacting particles and discuss methods that allow for a systematic and model-independent calculation of the hadron spectrum. These are lattice QCD and effective field theories. Synergies between both approaches can be exploited allowing for deeper understanding of the hadron spectrum. Results based on effective field theories and hadron-hadron scattering in lattice QCD or combinations thereof are presented and discussed. We also show that the often used Breit-Wigner parameterization is at odds with chiral symmetry and should not be used in case of strongly coupled channels. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: 75 pages, 21 figures, commissioned article for Physics Reports

Journal ref: Phys. Rept. 1001 (2023) 1-66

First moments of the nucleon transverse quark spin densities using lattice QCD

Authors: C. Alexandrou, S. Bacchio, M. Constantinou, P. Dimopoulos, J. Finkenrath, R. Frezzotti, K. Hadjiyiannakou, K. Jansen, B. Kostrzewa, G. Koutsou, G. Spanoudes, C. Urbach

Abstract: We present a calculation of the Mellin moments of the transverse quark spin densities in the nucleon using lattice QCD. The densities are extracted from the unpolarized and transversity generalized form factors extrapolated to the continuum limit using three $N_f=2+1+1$ twisted mass fermion gauge ensembles simulated with physical quark masses and spanning three lattice spacings. The first moment o… ▽ More We present a calculation of the Mellin moments of the transverse quark spin densities in the nucleon using lattice QCD. The densities are extracted from the unpolarized and transversity generalized form factors extrapolated to the continuum limit using three $N_f=2+1+1$ twisted mass fermion gauge ensembles simulated with physical quark masses and spanning three lattice spacings. The first moment of transversely polarized quarks in an unpolarized nucleon shows an interesting distortion, which can be traced back to the sharp falloff of the transversity generalized form factor $\bar{B}_{Tn0}(t)$. The isovector tensor anomalous magnetic moment is determined to be $κ_T=1.051(94)$, which confirms a negative and large Boer-Mulders function, $h_1^{\perp}$, in the nucleon. △ Less

Submitted 20 February, 2022; originally announced February 2022.

PLANETAMOS, A Physics Show Musical (Phyusical)

Authors: Lara Becker, Erik Busley, Jakob Dietl, Herbi K. Dreiner, Till Fohrmann, Kathrin Grunthal, Jana Heysel, Finn Jaekel, Kristoffer Kerkhof, Michael Kortmann, Barbara Leibrock, Viola Middelhauve, Steffi Moll, David Ohse, Johann Ostmeyer, Laura Rodríguez Gómez, Christoph Schürmann, Anne Stockhausen, Joshua Streichhahn, Carsten Urbach, Heinrich von Campe, Alexandra Wald, Laura Weber, Inga Woeste

Abstract: We present a physics show musical with live physics experiments and live performed songs with live orchestral accompaniment. The musical was first put on stage in German at the Physikalische Institut, University of Bonn, on the 24th of March, 2019. Here we present the original German script as well as an English translation, including a translation of the songs. We also give brief descriptions of… ▽ More We present a physics show musical with live physics experiments and live performed songs with live orchestral accompaniment. The musical was first put on stage in German at the Physikalische Institut, University of Bonn, on the 24th of March, 2019. Here we present the original German script as well as an English translation, including a translation of the songs. We also give brief descriptions of all the experiments employed, with photos as well as how we present the experiments. The musical tells the story of a couple, Life and Death, nicknamed Vita and Mortis, who want to buy a planet where they can settle down and raise a family. They arrive at the store Planetamos, run by the manager Luna Callisto and the scientist Jupi Mercury. The latter two present to Vita and Mortis what is required for life on a planet, as well as all of the extra features for sale, to make it that very special place for their family. Through Vita and Mortis the audience learns about the physical requirements for life on a planet orbiting a star. The show involves 14 live physics experiments. The four actors sing 7 songs, which are to well-known tunes with new lyrics, in German. During the show the singers are accompanied by a small live orchestra, consisting of a flute, an oboe, a piano, a violin, a cello, a harp and percussion. The show was developed by the University of Bonn physics students, who have all been very active in the University of Bonn Physics Show (Physikshow der Universitaet Bonn) the past few years. To the best of our knowledge this is the first physics show musical with live experiments and live song and orchestra. △ Less

Submitted 26 January, 2022; originally announced January 2022.

Comments: 86 pages, 18 figures

Digitising SU(2) Gauge Fields and the Freezing Transition

Authors: Tobias Hartung, Timo Jakobs, Karl Jansen, Johann Ostmeyer, Carsten Urbach

Abstract: Efficient discretisations of gauge groups are crucial with the long term perspective of using tensor networks or quantum computers for lattice gauge theory simulations. For any Lie group other than U$(1)$, however, there is no class of asymptotically dense discrete subgroups. Therefore, discretisations limited to subgroups are bound to lead to a freezing of Monte Carlo simulations at weak coupling… ▽ More Efficient discretisations of gauge groups are crucial with the long term perspective of using tensor networks or quantum computers for lattice gauge theory simulations. For any Lie group other than U$(1)$, however, there is no class of asymptotically dense discrete subgroups. Therefore, discretisations limited to subgroups are bound to lead to a freezing of Monte Carlo simulations at weak couplings, necessitating alternative partitionings without a group structure. In this work we provide a comprehensive analysis of this freezing for all discrete subgroups of SU$(2)$ and different classes of asymptotically dense subsets. We find that an appropriate choice of the subset allows unfrozen simulations for arbitrary couplings, though one has to be careful with varying weights of unevenly distributed points. A generalised version of the Fibonacci spiral appears to be particularly efficient and close to optimal. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Twisted mass gauge ensembles at physical values of the light, strange and charm quark masses

Authors: Jacob Finkenrath, Constantia Alexandrou, Simone Bacchio, Martha Constantinou, Petros Dimopoulos, Roberto Frezzotti, Karl Jansen, Bartosz Kostrzewa, Giannis Koutsou, Giancarlo Rossi, Carsten Urbach, Urs Wenger

Abstract: Lattice QCD simulations directly at physical masses of dynamical light, strange and charm quarks are highly desirable especially to remove systematic errors due to chiral extrapolations. However such simulations are still challenging. We discuss the adaption of efficient algorithms, like multi-grid methods or higher order integrators, within the molecular dynamic steps of the Hybrid Monte Carlo al… ▽ More Lattice QCD simulations directly at physical masses of dynamical light, strange and charm quarks are highly desirable especially to remove systematic errors due to chiral extrapolations. However such simulations are still challenging. We discuss the adaption of efficient algorithms, like multi-grid methods or higher order integrators, within the molecular dynamic steps of the Hybrid Monte Carlo algorithm, that are enabling simulations of a new set of gauge ensembles by the Extended Twisted Mass collaboration (ETMC). We present the status of the on-going ETMC simulation effort that aims to enabling studies of finite size and discretization effects. We work within the twisted mass discretization which is free of odd-discretization effects at maximal twist and present our tuning procedure. △ Less

Submitted 7 January, 2022; originally announced January 2022.

Comments: 12 pages, 12 figures, The 38th International Symposium on Lattice Field Theory (LATTICE2021), 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

Pion-pole contribution to HLbL from twisted mass lattice QCD at the physical point

Authors: Sebastian Burri, Constantia Alexandrou, Simone Bacchio, Georg Bergner, Jacob Finkenrath, Andrew Gasbarro, Kyriakos Hadjiyiannakou, Karl Jansen, Bartosz Kostrzewa, Giannis Koutsou, Konstantin Ottnad, Marcus Petschlies, Ferenc Pittler, Fernanda Steffens, Carsten Urbach, Urs Wenger

Abstract: We report on our computation of the pion transition form factor ${\cal F}_{P\rightarrow γ^*γ^*}$ from twisted mass lattice QCD in order to determine the numerically dominant light pseudoscalar pole contribution in the hadronic light-by-light scattering contribution to the anomalous magnetic moment of the muon $a_μ=(g-2)_μ$. The pion transition form factor is computed directly at the physical point… ▽ More We report on our computation of the pion transition form factor ${\cal F}_{P\rightarrow γ^*γ^*}$ from twisted mass lattice QCD in order to determine the numerically dominant light pseudoscalar pole contribution in the hadronic light-by-light scattering contribution to the anomalous magnetic moment of the muon $a_μ=(g-2)_μ$. The pion transition form factor is computed directly at the physical point. We present first results for our estimate of the pion-pole contribution with kinematic setup for the pion at rest. △ Less

Submitted 2 May, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: 9 pages, 7 figures; Proceedings of the 38th International Symposium on Lattice Field Theory - Lattice 2021, Zoom/Gather@Massachusetts Institute of Technology

Scattering from generalised lattice $φ^4$ theory

Authors: Marco Garofalo, Fernando Romero-López, Akaki Rusetsky, Carsten Urbach

Abstract: We investigate numerically different techniques to extract scattering amplitudes from the Euclidean Lattice φ4 theory with two fields, having different masses. We present an exploratory study of the recently proposed method by Bruno and Hansen for extracting the scattering length from a four-point function (cf 10.1007/JHEP06(2021)043), and a study of the two- and three-particle quantization condit… ▽ More We investigate numerically different techniques to extract scattering amplitudes from the Euclidean Lattice φ4 theory with two fields, having different masses. We present an exploratory study of the recently proposed method by Bruno and Hansen for extracting the scattering length from a four-point function (cf 10.1007/JHEP06(2021)043), and a study of the two- and three-particle quantization condition. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Comments: Contribution to the 38th International Symposium on Lattice Field Theory (LATTICE21), 26-30 July 2021, Zoom/Gather@Massachusetts Institute of Technology. This contribution overlap with arXiv:2107.04853

Report number: MIT-CTP/5366

Gradient flow scale-setting with $N_f=2+1+1$ Wilson-clover twisted-mass fermions

Authors: Constantia Alexandrou, Simone Bacchio, Georg Bergner, Petros Dimopoulos, Jacob Finkenrath, Roberto Frezzotti, Marco Garofalo, Bartosz Kostrzewa, Giannis Koutsou, Peter Labus, Francesco Sanfilippo, Silvano Simula, Martin Ueding, Carsten Urbach, Urs Wenger

Abstract: We present a determination of the gradient flow scales $w_0$, $\sqrt{t_0}$ and $t_0/w_0$ in isosymmetric QCD, making use of the gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f=2+1+1$ flavours of Wilson-clover twisted-mass quarks including configurations close to the physical point for all dynamical flavours. The simulations are carried out at three values of th… ▽ More We present a determination of the gradient flow scales $w_0$, $\sqrt{t_0}$ and $t_0/w_0$ in isosymmetric QCD, making use of the gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f=2+1+1$ flavours of Wilson-clover twisted-mass quarks including configurations close to the physical point for all dynamical flavours. The simulations are carried out at three values of the lattice spacing and the scale is set through the PDG value of the pion decay constant, yielding $w_0=0.17383(63)$ fm, $\sqrt{t_0}=0.14436(61)$ fm and $t_0/w_0=0.11969(62)$ fm. Finally, fixing the kaon mass to its isosymmetric value, we determine the ratio of the kaon and pion leptonic decay constants to be $f_K/f_π=1.1995(44)$. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: arXiv admin note: text overlap with arXiv:2104.06747

Hadronic Parity Violation from 4-quark Interactions

Authors: Aniket Sen, Marcus Petschlies, Nikolas Schlage, Carsten Urbach

Abstract: We present an exploratory investigation of the parity odd $ΔI = 1$ pion-nucleon coupling $h_π^1$ from lattice QCD. Based on the PCAC relation, we study the parity-conserving effective Hamiltonian and extract the coupling by determining the nucleon mass splitting arising from the effective 4-quark interactions using the Feynman-Hellmann theorem. We present preliminary results of the mass shift for… ▽ More We present an exploratory investigation of the parity odd $ΔI = 1$ pion-nucleon coupling $h_π^1$ from lattice QCD. Based on the PCAC relation, we study the parity-conserving effective Hamiltonian and extract the coupling by determining the nucleon mass splitting arising from the effective 4-quark interactions using the Feynman-Hellmann theorem. We present preliminary results of the mass shift for a $32^3 \times 64$ ensemble of $N_f = 2 + 1 + 1$ twisted mass fermions at pion mass $260$ MeV and lattice spacing $a = 0.097$ fm. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: LATTICE2021 proceedings

Lattice QCD results for the topological up-quark mass contribution: too small to rescue the $m_u=0$ solution to the strong CP problem

Authors: Constantia Alexandrou, Jacob Finkenrath, Lena Funcke, Karl Jansen, Bartosz Kostrzewa, Ferenc Pittler, Carsten Urbach

Abstract: A vanishing Yukawa coupling of the up quark could in principle solve the strong CP problem. To render this solution consistent with current algebra results, the up quark must receive an alternative mass contribution that conserves CP symmetry. Such a contribution could be provided by QCD through non-perturbative topological effects, including instantons. In this talk, we present the first direct l… ▽ More A vanishing Yukawa coupling of the up quark could in principle solve the strong CP problem. To render this solution consistent with current algebra results, the up quark must receive an alternative mass contribution that conserves CP symmetry. Such a contribution could be provided by QCD through non-perturbative topological effects, including instantons. In this talk, we present the first direct lattice computation of this topological mass contribution, using gauge configurations generated by the Extended Twisted Mass collaboration. We use the Iwasaki gauge action, Wilson twisted mass fermions at maximal twist, and dynamical up, down, strange and charm quarks. Our result for the topological mass contribution is an order of magnitude too small to account for the phenomenologically required up-quark mass. This rules out the "massless" up-quark solution to the strong CP problem, in accordance with previous results relying on $χ$PT fits to lattice data. The talk is based on [Alexandrou et al., PRL 125, 232001 (2020)], where more details can be found. △ Less

Submitted 30 November, 2021; v1 submitted 30 October, 2021; originally announced November 2021.

Comments: 12 pages, 3 figures, proceedings of the 38th International Symposium on Lattice Field Theory, 26th-30th July 2021, Zoom/Gather@Massachusetts Institute of Technology

Report number: MIT-CTP/5347

The Semimetal-Antiferromagnetic Mott Insulator Quantum Phase Transition of the Hubbard Model on the Honeycomb Lattice

Authors: Johann Ostmeyer, Evan Berkowitz, Stefan Krieg, Timo A. Lähde, Thomas Luu, Carsten Urbach

Abstract: The Hubbard model on the honeycomb lattice undergoes a quantum phase transition from a semimetallic to a Mott insulating phase and from a disordered to an anti-ferromagnetically phase. We show that these transitions occur simultaneously and we calculate the critical coupling $U_c=3.835(14)$ as well as the critical exponents $ν=1.181(43)$ and $β=0.898(37)$ which are expected to fall into the… ▽ More The Hubbard model on the honeycomb lattice undergoes a quantum phase transition from a semimetallic to a Mott insulating phase and from a disordered to an anti-ferromagnetically phase. We show that these transitions occur simultaneously and we calculate the critical coupling $U_c=3.835(14)$ as well as the critical exponents $ν=1.181(43)$ and $β=0.898(37)$ which are expected to fall into the $SU(2)$ Gross-Neveu universality class. For this we employ Hybrid Monte Carlo simulations, extrapolate the single particle gap and the spin structure factors to the thermodynamic and continuous time limits, and perform a data collapse fit. We also determine the zero temperature values of single particle gap and staggered magnetisation on both sides of the phase transition. △ Less

Submitted 28 October, 2021; originally announced October 2021.

Comments: LATTICE2021 proceedings, more details in arXiv:2005.11112 and arXiv:2105.06936

The Hubbard model with fermionic tensor networks

Authors: Manuel Schneider, Johann Ostmeyer, Karl Jansen, Thomas Luu, Carsten Urbach

Abstract: Many electromagnetic properties of graphene can be described by the Hubbard model on a honeycomb lattice. However, this system suffers strongly from the sign problem if a chemical potential is included. Tensor network methods are not affected by this problem. We use the imaginary time evolution of a fermionic projected entangled pair state, which allows to simulate both parity sectors independentl… ▽ More Many electromagnetic properties of graphene can be described by the Hubbard model on a honeycomb lattice. However, this system suffers strongly from the sign problem if a chemical potential is included. Tensor network methods are not affected by this problem. We use the imaginary time evolution of a fermionic projected entangled pair state, which allows to simulate both parity sectors independently. Incorporating the fermionic nature on the level of the tensor network allows to fix the particle number to be either even or odd. This way we can access the states at half filling and with one additional electron. We calculate the energy and other observables of both states, which was not possible before with Monte Carlo methods. △ Less

Submitted 23 December, 2021; v1 submitted 28 October, 2021; originally announced October 2021.

Comments: Proceedings of the 38th International Symposium on Lattice Field Theory (Lattice 2021)

Determination of the light, strange and charm quark masses using twisted mass fermions

Authors: C. Alexandrou, S. Bacchio, G. Bergner, M. Constantinou, M. Di Carlo, P. Dimopoulos, J. Finkenrath, E. Fiorenza, R. Frezzotti, M. Garofalo, K. Hadjiyiannakou, B. Kostrzewa, G. Koutsou, K. Jansen, V. Lubicz, M. Mangin-Brinet, F. Manigrasso, G. Martinelli, F. Pittler, G. C. Rossi, F. Sanfilippo, S. Simula, C. Tarantino, A. Todaro, C. Urbach , et al. (1 additional authors not shown)

Abstract: We present results for the light, strange and charm quark masses using $N_f=2+1+1$ twisted mass fermion ensembles at three values of the lattice spacing, including two ensembles simulated with the physical value of the pion mass. The analysis is done both in the meson and baryon sectors. The difference in the mean values found in the two sectors is included as part of the systematic error. The pre… ▽ More We present results for the light, strange and charm quark masses using $N_f=2+1+1$ twisted mass fermion ensembles at three values of the lattice spacing, including two ensembles simulated with the physical value of the pion mass. The analysis is done both in the meson and baryon sectors. The difference in the mean values found in the two sectors is included as part of the systematic error. The presentation is based on the work of Ref. [1], where more details can be found. △ Less

Submitted 9 October, 2021; originally announced October 2021.

Comments: Talk presented at the 38th International Symposium on Lattice Field Theory, Lattice 2021, Zoom/Gather@Massachusetts Institute of Technology, 26-30 Jul. 2021; 8 pages, 5 figures

Quark and gluon momentum fractions in the pion from $N_f=2+1+1$ lattice QCD

Authors: Constantia Alexandrou, Simone Bacchio, Georg Bergner, Jacob Finkenrath, Andrew Gasbarro, Kyriakos Hadjiyiannakou, Karl Jansen, Bartosz Kostrzewa, Konstantin Ottnad, Marcus Petschlies, Ferenc Pittler, Fernanda Steffens, Carsten Urbach, Urs Wenger

Abstract: We perform the first full decomposition of the pion momentum into its gluon and quark contributions. We employ an ensemble generated by the Extended Twisted Mass Collaboration with $N_f=2 + 1 +1$ Wilson twisted mass clover fermions at maximal twist tuned to reproduce the physical pion mass. We present our results in the $\overline{\mathrm{MS}}$ scheme at $2\gev$. We find $\avgx_{u+d}=0.601(28)$,… ▽ More We perform the first full decomposition of the pion momentum into its gluon and quark contributions. We employ an ensemble generated by the Extended Twisted Mass Collaboration with $N_f=2 + 1 +1$ Wilson twisted mass clover fermions at maximal twist tuned to reproduce the physical pion mass. We present our results in the $\overline{\mathrm{MS}}$ scheme at $2\gev$. We find $\avgx_{u+d}=0.601(28)$, $\avgx_s=0.059(13)$, $\avgx_c=0.019(05)$, and $\avgx_g=0.52(11)$ for the separate contributions, respectively, whose sum saturates the momentum sum rule. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Testing a new method for scattering in finite volume in the $φ^4$ theory

Authors: Marco Garofalo, Fernando Romero-López, Akaki Rusetsky, Carsten Urbach

Abstract: We test an alternative proposal by Bruno and Hansen [1] to extract the scattering length from lattice simulations in a finite volume. For this, we use a scalar $φ^4$ theory with two mass nondegenerate particles and explore various strategies to implement this new method. We find that the results are comparable to those obtained from the Lüscher method, with somewhat smaller statistical uncertainti… ▽ More We test an alternative proposal by Bruno and Hansen [1] to extract the scattering length from lattice simulations in a finite volume. For this, we use a scalar $φ^4$ theory with two mass nondegenerate particles and explore various strategies to implement this new method. We find that the results are comparable to those obtained from the Lüscher method, with somewhat smaller statistical uncertainties at larger volumes. △ Less

Submitted 10 July, 2021; originally announced July 2021.

Comments: 4 pages, 3 figures

Simulating both parity sectors of the Hubbard Model with Tensor Networks

Authors: Manuel Schneider, Johann Ostmeyer, Karl Jansen, Thomas Luu, Carsten Urbach

Abstract: Tensor networks are a powerful tool to simulate a variety of different physical models, including those that suffer from the sign problem in Monte Carlo simulations. The Hubbard model on the honeycomb lattice with non-zero chemical potential is one such problem. Our method is based on projected entangled pair states (PEPS) using imaginary time evolution. We demonstrate that it provides accurate es… ▽ More Tensor networks are a powerful tool to simulate a variety of different physical models, including those that suffer from the sign problem in Monte Carlo simulations. The Hubbard model on the honeycomb lattice with non-zero chemical potential is one such problem. Our method is based on projected entangled pair states (PEPS) using imaginary time evolution. We demonstrate that it provides accurate estimators for the ground state of the model, including cases where Monte Carlo simulations fail miserably. In particular it shows near to optimal, that is linear, scaling in lattice size. We also present a novel approach to directly simulate the subspace with an odd number of fermions. It allows to independently determine the ground state in both sectors. Without a chemical potential this corresponds to half filling and the lowest energy state with one additional electron or hole. We identify several stability issues, such as degenerate ground states and large single particle gaps, and provide possible fixes. △ Less

Submitted 2 July, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Comments: 20 pages, 20 figures

ACM Class: J.2

Journal ref: Phys. Rev. B 104 (2021) 155118

Beer Mats make bad Frisbees

Authors: Johann Ostmeyer, Christoph Schürmann, Carsten Urbach

Abstract: In this article we show why flying and rotating beer mats, CDs, or other flat disks will eventually flip in the air and end up flying with backspin, thus, making them unusable as frisbees. The crucial effect responsible for the flipping is found to be the lift attacking not in the center of mass but slightly offset to the forward edge. This induces a torque leading to a precession towards backspin… ▽ More In this article we show why flying and rotating beer mats, CDs, or other flat disks will eventually flip in the air and end up flying with backspin, thus, making them unusable as frisbees. The crucial effect responsible for the flipping is found to be the lift attacking not in the center of mass but slightly offset to the forward edge. This induces a torque leading to a precession towards backspin orientation. An effective theory is developed providing an approximate solution for the disk's trajectory with a minimal set of parameters. Our theoretical results are confronted with experimental results obtained using a beer mat shooting apparatus and a high speed camera. Very good agreement is found. △ Less

Submitted 28 July, 2021; v1 submitted 15 June, 2021; originally announced June 2021.

Comments: 4 videos in ancillary files

Journal ref: Eur. Phys. J. Plus 136, 769 (2021)

The Antiferromagnetic Character of the Quantum Phase Transition in the Hubbard Model on the Honeycomb Lattice

Authors: Johann Ostmeyer, Evan Berkowitz, Stefan Krieg, Timo A. Lähde, Thomas Luu, Carsten Urbach

Abstract: We provide a unified, comprehensive treatment of all operators that contribute to the anti-ferromagnetic, ferromagnetic, and charge-density-wave structure factors and order parameters of the hexagonal Hubbard Model. We use the Hybrid Monte Carlo algorithm to perform a systematic, carefully controlled analysis in the temporal Trotter error and of the thermodynamic limit. We expect our findings to i… ▽ More We provide a unified, comprehensive treatment of all operators that contribute to the anti-ferromagnetic, ferromagnetic, and charge-density-wave structure factors and order parameters of the hexagonal Hubbard Model. We use the Hybrid Monte Carlo algorithm to perform a systematic, carefully controlled analysis in the temporal Trotter error and of the thermodynamic limit. We expect our findings to improve the consistency of Monte Carlo determinations of critical exponents. We perform a data collapse analysis and determine the critical exponent $β=0.898(37)$ for the semimetal-Mott insulator transition in the hexagonal Hubbard Model. Our methods are applicable to a wide range of lattice theories of strongly correlated electrons. △ Less

Submitted 28 October, 2021; v1 submitted 14 May, 2021; originally announced May 2021.

Journal ref: Phys. Rev. B 104, 155142 (2021)

Quark masses using twisted mass fermion gauge ensembles

Authors: C. Alexandrou, S. Bacchio, G. Bergner, M. Constantinou, M. Di Carlo, P. Dimopoulos, J. Finkenrath, E. Fiorenza, R. Frezzotti, M. Garofalo, K. Hadjiyiannakou, B. Kostrzewa, G. Koutsou, K. Jansen, V. Lubicz, M. Mangin-Brinet, F. Manigrasso, G. Martinelli, E. Papadiofantous, F. Pittler, G. C. Rossi, F. Sanfilippo, S. Simula, C. Tarantino, A. Todaro , et al. (2 additional authors not shown)

Abstract: We present a calculation of the up, down, strange and charm quark masses performed within the lattice QCD framework. We use the twisted mass fermion action and carry out simulations that include in the sea two light mass-degenerate quarks, as well as the strange and charm quarks. In the analysis we use gauge ensembles simulated at three values of the lattice spacing and with light quarks that corr… ▽ More We present a calculation of the up, down, strange and charm quark masses performed within the lattice QCD framework. We use the twisted mass fermion action and carry out simulations that include in the sea two light mass-degenerate quarks, as well as the strange and charm quarks. In the analysis we use gauge ensembles simulated at three values of the lattice spacing and with light quarks that correspond to pion masses in the range from 350 MeV to the physical value, while the strange and charm quark masses are tuned approximately to their physical values. We use several quantities to set the scale in order to check for finite lattice spacing effects and in the continuum limit we get compatible results. The quark mass renormalization is carried out non-perturbatively using the RI'-MOM method converted into the $\overline{\rm MS}$ scheme. For the determination of the quark masses we use physical observables from both the meson and the baryon sectors, obtaining $m_{ud} = 3.636(66)(^{+60}_{-57})$~MeV and $m_s = 98.7(2.4)(^{+4.0}_{-3.2})$~MeV in the $\overline{\rm MS}(2\,{\rm GeV})$ scheme and $m_c = 1036(17)(^{+15}_{-8})$~MeV in the $\overline{\rm MS}(3\,{\rm GeV})$ scheme, where the first errors are statistical and the second ones are combinations of systematic errors. For the quark mass ratios we get $m_s / m_{ud} = 27.17(32)(^{+56}_{-38})$ and $m_c / m_s = 11.48(12)(^{+25}_{-19})$. △ Less

Submitted 11 October, 2021; v1 submitted 27 April, 2021; originally announced April 2021.

Comments: 37 pages, 23 figures, 24 tables. One reference added

Ratio of kaon and pion leptonic decay constants with $N_f = 2 + 1 + 1$ Wilson-clover twisted-mass fermions

Authors: C. Alexandrou, S. Bacchio, G. Bergner, P. Dimopoulos, J. Finkenrath, R. Frezzotti, M. Garofalo, B. Kostrzewa, G. Koutsou, P. Labus, F. Sanfilippo, S. Simula, M. Ueding, C. Urbach, U. Wenger

Abstract: We present a determination of the ratio of kaon and pion leptonic decay constants in isosymmetric QCD (isoQCD), $f_K / f_π$, making use of the gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks, including configurations close to the physical point for all dynamical flavors. The simulations are carried out a… ▽ More We present a determination of the ratio of kaon and pion leptonic decay constants in isosymmetric QCD (isoQCD), $f_K / f_π$, making use of the gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks, including configurations close to the physical point for all dynamical flavors. The simulations are carried out at three values of the lattice spacing ranging from $\sim 0.068$ to $\sim 0.092$ fm with linear lattice size up to $L \sim 5.5$~fm. The scale is set by the PDG value of the pion decay constant, $f_π^{isoQCD} = 130.4~(2)$ MeV, at the isoQCD pion point, $M_π^{isoQCD} = 135.0~(2)$ MeV, obtaining for the gradient-flow (GF) scales the values $w_0 = 0.17383~(63)$ fm, $\sqrt{t_0} = 0.14436~(61)$ fm and $t_0 / w_0 = 0.11969~(62)$ fm. The data are analyzed within the framework of SU(2) Chiral Perturbation Theory (ChPT) without resorting to the use of renormalized quark masses. At the isoQCD kaon point $M_K^{isoQCD} = 494.2~(4)$ MeV we get $(f_K / f_π)^{isoQCD} = 1.1995~(44)$, where the error includes both statistical and systematic uncertainties. Implications for the Cabibbo-Kobayashi-Maskawa (CKM) matrix element $|V_{us}|$ and for the first-row CKM unitarity are discussed. △ Less

Submitted 3 October, 2021; v1 submitted 14 April, 2021; originally announced April 2021.

Comments: 68 pages, 14 figures, 12 tables. Version to appear in PRD

Nucleon axial and pseudoscalar form factors from lattice QCD at the physical point

Authors: C. Alexandrou, S. Bacchio, M. Constantinou, P. Dimopoulos, J. Finkenrath, K. Hadjiyiannakou, K. Jansen, G. Koutsou, B. Kostrzewa, T. Leontiou, C. Urbach

Abstract: We compute the nucleon axial and induced pseudoscalar form factors using three ensembles of gauge configurations, generated with dynamical light quarks with mass tuned to approximately their physical value. One of the ensembles also includes the strange and charm quarks with their mass close to physical. The latter ensemble has large statistics and finer lattice spacing and it is used to obtain fi… ▽ More We compute the nucleon axial and induced pseudoscalar form factors using three ensembles of gauge configurations, generated with dynamical light quarks with mass tuned to approximately their physical value. One of the ensembles also includes the strange and charm quarks with their mass close to physical. The latter ensemble has large statistics and finer lattice spacing and it is used to obtain final results, while the other two are used for assessing volume effects. The pseudoscalar form factor is also computed using these ensembles. We examine the momentum dependence of these form factors as well as relations based on pion pole dominance and the partially conserved axial-vector current hypothesis. △ Less

Submitted 26 November, 2020; originally announced November 2020.

Comments: 25 pages and 27 figures

Journal ref: Phys. Rev. D 103, 034509 (2021)

Relativistic $N$-particle energy shift in finite volume

Authors: Fernando Romero-López, Akaki Rusetsky, Nikolas Schlage, Carsten Urbach

Abstract: We present a general method for deriving the energy shift of an interacting system of $N$ spinless particles in a finite volume. To this end, we use the nonrelativistic effective field theory (NREFT), and match the pertinent low-energy constants to the scattering amplitudes. Relativistic corrections are explicitly included up to a given order in the $1/L$ expansion. We apply this method to obtain… ▽ More We present a general method for deriving the energy shift of an interacting system of $N$ spinless particles in a finite volume. To this end, we use the nonrelativistic effective field theory (NREFT), and match the pertinent low-energy constants to the scattering amplitudes. Relativistic corrections are explicitly included up to a given order in the $1/L$ expansion. We apply this method to obtain the ground state of $N$ particles, and the first excited state of two and three particles to order $L^{-6}$ in terms of the threshold parameters of the two- and three-particle relativistic scattering amplitudes. We use these expressions to analyze the $N$-particle ground state energy shift in the complex $\varphi^4$ theory. △ Less

Submitted 1 March, 2021; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 52 pages, 12 figures, 8 tables. Typos in Appendix A removed, numerical values for the integrals added. Results of the paper are unaffected

Scattering of two and three physical pions at maximal isospin from lattice QCD

Authors: Matthias Fischer, Bartosz Kostrzewa, Liuming Liu, Fernando Romero-López, Martin Ueding, Carsten Urbach

Abstract: We present the first direct $N_f=2$ lattice QCD computation of two- and three-$π^+$ scattering quantities that includes an ensemble at the physical point. We study the quark mass dependence of the two-pion phase shift, and the three-particle interaction parameters. We also compare to phenomenology and chiral perturbation theory (ChPT). In the two-particle sector, we observe good agreement to the p… ▽ More We present the first direct $N_f=2$ lattice QCD computation of two- and three-$π^+$ scattering quantities that includes an ensemble at the physical point. We study the quark mass dependence of the two-pion phase shift, and the three-particle interaction parameters. We also compare to phenomenology and chiral perturbation theory (ChPT). In the two-particle sector, we observe good agreement to the phenomenological fits in $s$- and $d$-wave, and obtain $M_πa_0 = -0.0481(86)$ at the physical point from a direct computation. In the three-particle sector, we observe reasonable agreement at threshold to the leading order chiral expansion, i.e.\@ a mildly attractive three-particle contact term. In contrast, we observe that the energy-dependent part of the three-particle quasilocal scattering quantity is not well described by leading order ChPT. △ Less

Submitted 25 May, 2021; v1 submitted 7 August, 2020; originally announced August 2020.

Comments: Published version

The $ρ$-resonance with physical pion mass from $N_f=2$ lattice QCD

Authors: Matthias Fischer, Bartosz Kostrzewa, Maxim Mai, Marcus Petschlies, Ferenc Pittler, Martin Ueding, Carsten Urbach, Markus Werner

Abstract: We present the first-ever lattice computation of pi pi-scattering in the I=1 channel with Nf=2 dynamical quark flavours obtained including an ensemble with physical value of the pion mass. Employing a global fit to data at three values of the pion mass, we determine the universal parameters of the rho-resonance. We carefully investigate systematic uncertainties by determining energy eigenvalues us… ▽ More We present the first-ever lattice computation of pi pi-scattering in the I=1 channel with Nf=2 dynamical quark flavours obtained including an ensemble with physical value of the pion mass. Employing a global fit to data at three values of the pion mass, we determine the universal parameters of the rho-resonance. We carefully investigate systematic uncertainties by determining energy eigenvalues using different methods and by comparing inverse amplitude method and Breit-Wigner type parametrizations. Overall, we find mass 786(20) MeV and width 180(6) MeV, including statistical and systematic uncertainties. In stark disagreement with the previous Nf=2 extrapolations from higher than physical pion mass results, our mass value is in good agreement with experiment, while the width is slightly too high. △ Less

Submitted 24 June, 2020; originally announced June 2020.

Comments: raw data available in data repository: https://github.com/urbach/datarhonf2

The Semimetal-Mott Insulator Quantum Phase Transition of the Hubbard Model on the Honeycomb Lattice

Authors: Johann Ostmeyer, Evan Berkowitz, Stefan Krieg, Timo A. Lähde, Thomas Luu, Carsten Urbach

Abstract: We take advantage of recent improvements in the grand canonical Hybrid Monte Carlo algorithm, to perform a precision study of the single-particle gap in the hexagonal Hubbard model, with on-site electron-electron interactions. After carefully controlled analyses of the Trotter error, the thermodynamic limit, and finite-size scaling with inverse temperature, we find a critical coupling of… ▽ More We take advantage of recent improvements in the grand canonical Hybrid Monte Carlo algorithm, to perform a precision study of the single-particle gap in the hexagonal Hubbard model, with on-site electron-electron interactions. After carefully controlled analyses of the Trotter error, the thermodynamic limit, and finite-size scaling with inverse temperature, we find a critical coupling of $U_c/κ=3.834(14)$ and the critical exponent $zν=1.185(43)$. Under the assumption that this corresponds to the expected anti-ferromagnetic Mott transition, we are also able to provide a preliminary estimate $β=1.095(37)$ for the critical exponent of the order parameter. We consider our findings in view of the $SU(2)$ Gross-Neveu, or chiral Heisenberg, universality class. We also discuss the computational scaling of the Hybrid Monte Carlo algorithm, and possible extensions of our work to carbon nanotubes, fullerenes, and topological insulators. △ Less

Submitted 7 December, 2020; v1 submitted 22 May, 2020; originally announced May 2020.

Journal ref: Phys. Rev. B 102, 245105 (2020)