Finite Volume Corrections Research Articles

Magnetic polarizability of octet baryons via lattice QCD

Drawing on recent advances in lattice QCD background-field techniques, the magnetic polarizability of octet baryons is calculated from the first principles of QCD. The results are presented in the context of new constituent quark model calculations providing a framework for understanding the lattice results and a direct comparison with simulation results at unphysical quark masses. Using smeared quark sources, low-lying Laplacian eigenmode projection and final-state Landau-mode projection, considerable attention is devoted to ensuring single-state isolation in the lattice correlation functions. We also introduce new weighting methods to reduce the sensitivity to correlation-function fits, averaging over many fits based on merit drawn from the full correlated χ2 of the fits. The techniques are implemented on the 323×64, 2+1-flavor dynamical-fermion lattices provided by the PACS-CS Collaboration following the introduction of uniform magnetic fields quantized to the lowest nontrivial values available. After some scaling of the constituent quark model parameters, we find the model captures the patterns observed in the lattice QCD results very well, providing important insights into the physics underpinning the magnetic polarizabilities. Finally, comparison with the most recent results from experiment proceeds through an effective field theory formalism which incorporates estimates of finite-volume corrections and small electroquenching corrections as the results are brought to the physical point. We find excellent agreement with experiment where available, including the proton and neutron polarizabilities. Published by the American Physical Society 2024

Discontinuous Galerkin schemes for hyperbolic systems in non-conservative variables: Quasi-conservative formulation with subcell finite volume corrections

We present a novel quasi-conservative arbitrary high order accurate ADER (Arbitrary-Derivative) discontinuous Galerkin method allowing to efficiently use a non-conservative form of the considered partial differential system, so that the governing equations can be solved directly in the most physically relevant set of variables. This is particularly interesting for multi-material flows with moving interfaces and steep, large magnitude contact discontinuities, as well as in presence of highly non-linear thermodynamics. However, the non-conservative formulation of course introduces a conservation error which would normally lead to a wrong approximation of shock waves. Hence, from the theoretical point of view, we give a formal definition of the conservation defect of non-conservative schemes and we analyze this defect providing a local quasi-conservation condition, which allows us to prove a modified Lax–Wendroff theorem. Within this formalism, we also reformulate classical results concerning smooth solutions, contact discontinuities and moving interfaces. Then, to deal with shock waves in practice, we exploit the framework of the so-called a posteriori subcell finite volume (FV) limiter, so that, in troubled cells appropriately detected, we can incorporate a local conservation correction. Our corrected FV update entirely removes the local conservation defect, allowing, at least formally, to fit in the hypotheses of the proposed modified Lax–Wendroff theorem. Here, the shock-triggered troubled cells are detected by combining physical admissibility criteria, a discrete maximum principle and a shock sensor inspired by Lagrangian hydrodynamics.To prove the capabilities of our novel approach, first, we show that we are able to recover the same results given by conservative schemes on classical benchmarks for the single-fluid Euler equations. We then conclude the presentation by showing the improved reliability of our scheme on the multi-fluid Euler system on examples like the interaction of a shock with a helium bubble for which we are able to avoid the development of any spurious oscillations.

Quantum electrodynamics on the lattice and numerical perturbative computation of g − 2

Abstract We compute the electron g factor to the $\mathcal {O}(\alpha ^5)$ order on the lattice in quenched quantum electrodynamics (QED). We first study finite volume (FV) corrections in various infrared regularization methods to discuss which regularization is optimal for our purpose. We find that in QEDL the FV correction to the effective mass can have different parametric dependences depending on the size of Euclidean time t and match the ‘naive on-shell result’ only at the very large t region, t ≫ L. We adopt finite photon mass regularization to suppress FV effects exponentially and also discuss our strategy for selecting simulation parameters and the order of extrapolations to efficiently obtain the g factor. We perform lattice simulation using small lattices to test the feasibility of our calculation strategy. This study can be regarded as an intermediate step toward giving the five-loop coefficient independently of preceding studies.

QCD thermodynamics and neutral pion in a uniform magnetic field: Finite volume effects

We address finite volume effects of lattice QCD calculations in background magnetic fields. Using chiral perturbation theory at next-to-leading order, volume effects are calculated for thermodynamic quantities: the chiral condensate, pressure anisotropy, and magnetization. The neutral pion effective action in a finite volume is additionally derived. For these charge neutral observables, volume and source averaging are shown to capitalize on magnetic periodicity, which is the remnant translational invariance of the finite-volume theory. For a fixed magnetic field strength, certain volume and source averaged quantities are independent of the size of the lattice transverse to the magnetic field. Despite this simplifying feature, finite volume corrections to the magnetic field dependence of the chiral condensate and neutral pion magnetic polarizability can be non-negligible. The pressure anisotropy at fixed magnetic flux, moreover, appears acutely sensitive to the lattice volume.

Isospin-breaking corrections to light-meson leptonic decays from lattice simulations at physical quark masses

The decreasing uncertainties in theoretical predictions and experimental measurements of several hadronic observables related to weak processes, which in many cases are now smaller than O(1%), require theoretical calculations to include subleading corrections that were neglected so far. Precise determinations of leptonic and semi-leptonic decay rates, including QED and strong isospin-breaking effects, can play a central role in solving the current tensions in the first-row unitarity of the CKM matrix. In this work we present the first RBC/UKQCD lattice calculation of the isospin-breaking corrections to the ratio of leptonic decay rates of kaons and pions into muons and neutrinos. The calculation is performed at fixed lattice spacing (a−1 ≃ 1.730 GeV) on a 483× 96 volume with Nf = 2 + 1 dynamical quarks close to the physical point and domain wall fermions in the Möbius formulation are employed. Long-distance QED interactions are included according to the QEDL prescription and the crucial role of finite-volume electromagnetic corrections in the determination of leptonic decay rates, which produce a large systematic uncertainty, is extensively discussed. Finally, we study the different sources of uncertainty on |Vus|/|Vud| and observe that, if finite-volume systematics can be reduced, the error from isospin-breaking corrections is potentially sub-dominant in the final precision of the ratio of the CKM matrix elements.

Window observable for the hadronic vacuum polarization contribution to the muon g−2 from lattice QCD

Euclidean time windows in the integral representation of the hadronic vacuum polarization contribution to the muon $g-2$ serve to test the consistency of lattice calculations and may help in tracing the origins of a potential tension between lattice and data-driven evaluations. In this paper, we present results for the intermediate time window observable computed using O($a$) improved Wilson fermions at six values of the lattice spacings below 0.1\,fm and pion masses down to the physical value. Using two different sets of improvement coefficients in the definitions of the local and conserved vector currents, we perform a detailed scaling study which results in a fully controlled extrapolation to the continuum limit without any additional treatment of the data, except for the inclusion of finite-volume corrections. To determine the latter, we use a combination of the method of Hansen and Patella and the Meyer-Lellouch-L\"uscher procedure employing the Gounaris-Sakurai parameterization for the pion form factor. We correct our results for isospin-breaking effects via the perturbative expansion of QCD+QED around the isosymmetric theory. Our result at the physical point is $a_\mu^{\mathrm{win}}=(237.30\pm0.79_{\rm stat}\pm1.22_{\rm syst})\times10^{-10}$, where the systematic error includes an estimate of the uncertainty due to the quenched charm quark in our calculation. Our result displays a tension of 3.9$\sigma$ with a recent evaluation of $a_\mu^{\mathrm{win}}$ based on the data-driven method.

A unified formulation of one-loop tensor integrals for finite volume effects

A unified formulation of one-loop tensor integrals is proposed for systematical calculations of finite volume corrections. It is shown that decomposition of the one-loop tensor integrals into a series of tensors accompanied by tensor coefficients is feasible, if a unit space-like four vector nμ, originating from the discretization effects at finite volume, is introduced. A generic formula has been derived for numerical computations of all the involved tensor coefficients. For the vanishing external three-momenta, we also investigate the feasibility of the conventional Passarino-Veltmann reduction of the tensor integrals in a finite volume. Our formulation can be easily used to realize the automation of the calculations of finite volume corrections to any interesting quantities at one-loop level. Besides, it provides finite volume result in a unique and concise form, which is suited for, e.g., carrying out precision determination of physical observable from modern lattice QCD data.

π−π scattering, QED, and finite-volume quantization

Using the Coulomb gauge formulation of QED we present a lattice QCD procedure to calculate the $\pi^+\pi^+$ scattering phase shift including the effects of the Coulomb potential which appears in this formulation. The approach described here incorporates the effects of relativity and avoids finite-volume corrections that vanish as a power of the volume in which the lattice calculation is performed. This is the first step in developing a complete lattice QCD calculation of the electromagnetic and isospin-breaking light-quark mass contributions to $\varepsilon'$, the parameter describing direct CP violating effects in $K_L\to\pi\pi$ decay.

Finite volume effects in the McLerran–Venugopalan initial condition for the JIMWLK equation

We revisit the numerical construction of the initial condition for the dipole amplitude from the McLerran–Venugopalan model in the context of the JIMWLK evolution equation. We observe large finite volume effects induced by the Poisson equation formulated on a torus. We show that the situation can be partially cured by introducing an infrared regularization. We propose a procedure which has negligible finite volume corrections. The control of the finite volume and finite lattice spacings effects is crucial when considering the numerical solutions of the JIMWLK evolution equation with the collinear improvement.

Correlations and Critical Behavior in Lattice Gluodynamics

In the Landau-gauge lattice gluodynamics we find that, both in the SU(2) and SU(3) theory, a correlation of the Polyakov loop with the asymmetry of the A2 gluon condensate as well as with the longitudinal propagator makes it possible to determine the critical behavior of these quantities. We discuss finitevolume corrections and reveal that they can be reduced by the use of regression analysis. We also analyze the temperature dependence of low-momenta propagators in different Polyakov-loop sectors.

Model-independent framework for determining finite-volume effects of spatially nonlocal operators

We present a model-independent framework to determine finite-volume corrections of matrix elements of spatially-separated current-current operators. We define these matrix elements in terms of Compton-like amplitudes, i.e. amplitudes coupling single-particle states via two current insertions. We show that the infrared behavior of these matrix elements is dominated by the single-particle pole, which is approximated by the elastic form factors of the lowest-lying hadron. Therefore, given lattice data on the relevant elastic form factors, the finite-volume effects can be estimated non-perturbatively and without recourse to effective field theories. For illustration purposes, we investigate the implications of the proposed formalism for a class of scalar theories in two and four dimensions.

Nucleon momentum fraction, helicity and transversity from 2+1-flavor lattice QCD

A detailed analysis of the systematic uncertainties in the calculation of the isovector momentum fraction, 〈x〉u − d, helicity moment, 〈x〉Δu − Δd, and the transversity moment, 〈x〉δu − δd, of the nucleon is presented using high-statistics data on seven ensembles of gauge configurations generated by the JLab/W&M/LANL/MIT collaborations using 2 + 1-flavors of dynamical Wilson-clover quarks. The much higher statistics have facilitated better control over all systematics compared to previous lattice calculations. The least understood systematic — excited-state contamination — is quantified by studying the variation of the results as a function of different estimates of the mass gap of the first excited state, obtained from two- and three-point correlation functions, and as a function of the pion mass Mπ. The final results are obtained using a simultaneous fit in the lattice spacing a, pion mass Mπ and the finite volume parameter MπL keeping leading order corrections. The data show no significant dependence on the lattice spacing and some evidence for finite-volume corrections. Our final results, in the overline{mathrm{MS}} scheme at 2 GeV, are 〈x〉u − d = 0.155(17)(20), 〈x〉Δu − Δd = 0.183(14)(20) and 〈x〉δu − δd = 0.220(18)(20), where the first error is the overall analysis uncertainty assuming excited-state contributions have been removed, and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with phenomenological global fit values.

Finite volume corrections to forward Compton scattering off the nucleon

We calculate the spin-averaged amplitude for doubly virtual forward Compton scattering off nucleons in the framework of manifestly Lorentz invariant baryon chiral perturbation theory at complete one-loop order $O(p^4)$. The calculations are carried out both in the infinite and in a finite volume. The obtained results allow for a detailed estimation the finite-volume corrections to the amplitude which can be extracted on the lattice using the background field technique.

Finite-volume effects in long-distance processes with massless leptonic propagators

In Ref. [1], a method was proposed to calculate QED corrections to hadronic self energies from lattice QCD without power-law finite-volume errors. In this paper, we extend the method to processes which occur at second-order in the weak interaction and in which there is a massless (or almost massless) leptonic propagator. We demonstrate that, in spite of the presence of the propagator of an almost massless electron, such an infinite-volume reconstruction procedure can be used to obtain the amplitude for the rare kaon decay $K^+\to\pi^+\nu\bar\nu$ from a lattice quantum chromodynamics computation with only exponentially small finite-volume corrections.

P-wave two-particle bound and scattering states in a finite volume including QED

The mass shifts for two-fermion bound and scattering P-wave states subject to the long-range interactions due to QED in the non-relativistic regime are derived. Introducing a short range force coupling the spinless fermions to one unit of angular momentum in the framework of pionless EFT, we first calculate both perturbatively and non-perturbatively the Coulomb corrections to fermion–fermion scattering in the continuum and infinite volume context. Motivated by the research on particle–antiparticle bound states, we extend the results to fermions of identical mass and opposite charge. Second, we transpose the system onto a cubic box with periodic boundary conditions and we calculate the finite volume corrections to the energy of the lowest bound and unbound T_1^{-} eigenstates. In particular, power law corrections proportional to the fine structure constant and resembling the recent results for S-wave states are found. Higher order contributions in alpha are neglected, since the gapped nature of the momentum operator in the finite-volume environnement allows for a perturbative treatment of the QED interactions.

Application of effective field theory to finite-volume effects in aμHVP

One of the more important systematic effects affecting lattice computations of the hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon, $a_\mu^{\rm HVP}$, is the distortion due to a finite spatial volume. In order to reach sub-percent precision, these effects need to be reliably estimated and corrected for, and one of the methods that has been employed for doing this is finite-volume chiral perturbation theory. In this paper, we argue that finite-volume corrections to $a_\mu^{\rm HVP}$ can, in principle, be calculated at any given order in chiral perturbation theory. More precisely, once all low-energy constants needed to define the Effective Field Theory representation of $a_\mu^{\rm HVP}$ in infinite volume are known to a given order, also the finite-volume corrections can be predicted to that order in the chiral expansion.

Finite-volume and thermal effects in the leading-HVP contribution to muonic (g \u2212 2)

The leading finite-volume and thermal effects, arising in numerical lattice QCD calculations of {a}_{mu}^{mathrm{HVP},mathrm{LO}}equiv {left(g-2right)}_{mu}^{mathrm{HVP},mathrm{LO}}/2 , are determined to all orders with respect to the interactions of a generic, relativistic effective field theory of pions. In contrast to earlier work [1] based in the finite-volume Hamiltonian, the results presented here are derived by formally summing all Feynman diagrams contributing to the Euclidean electromagnetic-current two-point function, with any number of internal pion loops and interaction vertices. As was already found in ref. [1], the leading finite-volume corrections to {a}_{mu}^{mathrm{HVP},mathrm{LO}} scale as exp[−mL] where m is the pion mass and L is the length of the three periodic spatial directions. In this work we additionally control the two sub-leading exponentials, scaling as exp[− sqrt{2} mL] and exp[− sqrt{3} mL]. As with the leading term, the coefficient of these is given by the forward Compton amplitude of the pion, meaning that all details of the effective theory drop out of the final result. Thermal effects are additionally considered, and found to be sub-percent-level for typical lattice calculations. All finite-volume corrections are presented both for {a}_{mu}^{mathrm{HVP},mathrm{LO}} and for each time slice of the two-point function, with the latter expected to be particularly useful in correcting small to intermediate current separations, for which the series of exponentials exhibits good convergence.

Scaling of the formation probabilities and universal boundary entropies in the quantum XY spin chain

We calculate exactly the probability to find the ground state of the XY chain in a given spin configuration in the transverse σz-basis. By determining finite-volume corrections to the probabilities for a wide variety of configurations, we obtain the universal boundary entropy at the critical point. The latter is a benchmark of the underlying boundary conformal field theory characterizing each quantum state. To determine the scaling of the probabilities, we prove a theorem that expresses, in a factorized form, the eigenvalues of a sub-matrix of a circulant matrix as functions of the eigenvalues of the original matrix. Finally, the boundary entropies are computed by exploiting a generalization of the Euler–MacLaurin formula to non-differentiable functions. It is shown that, in some cases, the spin configuration can flow to a linear superposition of Cardy states. Our methods and tools are rather generic and can be applied to all the periodic quantum chains which map to free-fermionic Hamiltonians.

Constraints on charm-anticharm asymmetry in the nucleon from lattice QCD

We present the first lattice QCD calculation of the charm quark contribution to the nucleon electromagnetic form factors GE,Mc(Q2) in the momentum transfer range 0≤Q2≤1.4GeV2. The quark mass dependence, finite lattice spacing and volume corrections are taken into account simultaneously based on the calculation on three gauge ensembles including one at the physical pion mass. The nonzero value of the charm magnetic moment μMc=−0.00127(38)stat(5)sys, as well as the Pauli form factor, reflects a nontrivial role of the charm sea in the nucleon spin structure. The nonzero GEc(Q2) indicates the existence of a nonvanishing asymmetric charm-anticharm sea in the nucleon. Performing a nonperturbative analysis based on holographic QCD and the generalized Veneziano model, we study the constraints on the [c(x)−c¯(x)] distribution from the lattice QCD results presented here. Our results provide complementary information and motivation for more detailed studies of physical observables that are sensitive to intrinsic charm and for future global analyses of parton distributions including asymmetric charm-anticharm distribution.

Consistency of hadronic vacuum polarization between lattice QCD and the R ratio

There are emerging tensions for theory results of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment both within recent lattice QCD calculations and between some lattice QCD calculations and R-ratio results. In this paper we work towards scrutinizing critical aspects of these calculations. We focus in particular on a precise calculation of Euclidean position-space windows defined by RBC/UKQCD that are ideal quantities for cross-checks within the lattice community and with R-ratio results. We perform a lattice QCD calculation using physical up, down, strange, and charm sea quark gauge ensembles generated in the staggered formalism by the MILC collaboration. We study the continuum limit using inverse lattice spacings from $a^{-1}\approx 1.6$ GeV to $3.5$ GeV, identical to recent studies by FNAL/HPQCD/MILC and Aubin et al. and similar to the recent study of BMW. Our calculation exhibits a tension for the particularly interesting window result of $a_\mu^{\rm ud, conn.,isospin, W}$ from $0.4$ fm to $1.0$ fm with previous results obtained with a different discretization of the vector current on the same gauge configurations. Our results may indicate a difficulty related to estimating uncertainties of the continuum extrapolation that deserves further attention. In this work we also provide results for $a_\mu^{\rm ud,conn.,isospin}$, $a_\mu^{\rm s,conn.,isospin}$, $a_\mu^{\rm SIB,conn.}$ for the total contribution and a large set of windows. For the total contribution, we find $a_\mu^{\rm HVP~LO}=714(27)(13) 10^{-10}$, $a_\mu^{\rm ud,conn.,isospin}=657(26)(12) 10^{-10}$, $a_\mu^{\rm s,conn.,isospin}=52.83(22)(65) 10^{-10}$, and $a_\mu^{\rm SIB,conn.}=9.0(0.8)(1.2) 10^{-10}$, where the first uncertainty is statistical and the second systematic. We also comment on finite-volume corrections for the strong-isospin-breaking corrections.