Skip to content

Lattice Observables and Continuum Inference

A lattice observable claim is only as strong as its weakest inference step. A correlation function is not an energy, a fitted amplitude is not yet a matrix element, a bare matrix element is not renormalized, and a renormalized finite-spacing quantity is not a continuum prediction. This chapter connects those objects without collapsing their uncertainties or assumptions.

The recurring chain is

operators⟶C2,C3⟶levels and bare matrix elements⟶renormalization and scale⟶continuum observable.\text{operators} \longrightarrow C_2,C_3 \longrightarrow \text{levels and bare matrix elements} \longrightarrow \text{renormalization and scale} \longrightarrow \text{continuum observable}.

Every arrow has a failure test. The pages below let you enter at the missing arrow while keeping the complete chain visible.

QuestionPageRequired output
What do the measured Euclidean correlators encode before fitting?Euclidean Correlators and Spectral InformationNormalized finite-volume spectral sum, boundary images, overlaps, and noise limits
Can the operator basis isolate the desired level?Operator Bases, Effective Masses, and Excited-State ControlCorrelation matrix, conditioning analysis, GEVP or multi-state result, and competing-model tests
How is a bare matrix element obtained in the presence of disconnected terms?Three-Point Functions, Matrix Elements, and Disconnected ContributionsConnected/disconnected decomposition, excited-state control, stochastic validation, and bare result
What continuous spectral information is identifiable from Euclidean data?Spectral Reconstruction and Ill-Posed Euclidean Inverse ProblemsResolution function, regularization assumptions, mock recovery, and a bounded claim
How are lattice units converted without circularity?Scale Setting and Dimensionless RatiosIndependent reference, unit conversion, covariance, and held-out prediction
How does a bare operator reach a named continuum scheme and scale?Nonperturbative Renormalization, Mixing, and Step ScalingMixing matrix, renormalization condition, step-scaling chain, conversion, and window tests
Which lattice-side loop coefficients predict cutoff improvement?Lattice Perturbation Theory, Symanzik Analysis, and ImprovementLattice Feynman rules, Brillouin-zone integral, matched coefficient, allowed improvement basis, and residual power law
Does the complete multi-spacing analysis support a continuum claim?Lines of Constant Physics and Continuum ExtrapolationCorrelated fit with tuned inputs, volume treatment, alternative artifact models, and stability tests

Sampling algorithms and generic correlated-data methodology are developed in Sampling Algorithms for Lattice Fields and Statistical Inference and Error Budgets. Here they enter only through the properties needed to interpret a declared observable.

Choose interpolating operators Oi\mathcal O_i with fixed lattice quantum numbers and zero-temperature transfer evolution. This overview uses unit-normalized finite-volume states. For 0<t<T0<t<T, a two-point matrix has the low-temperature decomposition

Cij(t)=⟨Oi(t)Oj†(0)⟩=∑nZi(n)Zj(n)∗e−ΔEnt+∑nYi(n)Yj(n)∗e−ΔEn(T−t)+O(e−TΔE1),C_{ij}(t) =\langle\mathcal O_i(t)\mathcal O_j^\dagger(0)\rangle =\sum_n Z_i^{(n)}Z_j^{(n)*}e^{-\Delta E_nt} +\sum_n Y_i^{(n)}Y_j^{(n)*}e^{-\Delta E_n(T-t)} +O(e^{-T\Delta E_1}),

where ΔEn=En−E0\Delta E_n=E_n-E_0, Zi(n)=⟨0∣Oi∣n⟩Z_i^{(n)}=\langle0|\mathcal O_i|n\rangle, and Yi(n)=⟨n∣Oi∣0⟩Y_i^{(n)}=\langle n|\mathcal O_i|0\rangle. What matters is that one state normalization is used consistently in two- and three-point functions. At finite TT the backward image is part of the model. At finite temperature, additional transitions between nonvacuum states appear; this vacuum-dominance form must not be reused without a suppression argument.

Only for a self-adjoint or suitably self-conjugate basis does reflection relate YY to ZZ and reduce a diagonal correlator to a common-amplitude cosh. For a charged operator or a basis not closed under conjugation, backward propagation belongs to the reversed channel and can carry different overlaps; the exact thermal trace is the safe starting point. If relativistically normalized states are used instead, completeness supplies inverse state-normalization factors in both C2C_2 and C3C_3, and the reported matrix element must be converted with the same convention.

A three-point function with insertion JJ at time τ\tau and sink at TsT_s is

C3,ij(Ts,τ)=∑m,nZi(m)Zj(n)∗e−Em(Ts−τ)e−Enτ⟨m∣J∣n⟩,C_{3,ij}(T_s,\tau) =\sum_{m,n} Z_i^{(m)}Z_j^{(n)*} e^{-E_m(T_s-\tau)}e^{-E_n\tau} \langle m|J|n\rangle,

again up to the declared state normalization. A plateau approximation keeps only m=n=0m=n=0; its error is controlled by gaps times both τ\tau and Ts−τT_s-\tau. Increasing only one side cannot suppress both contaminations.

The chapter’s first four pages determine what may be inferred from these kernels. The next four decide whether the inferred bare quantities support a matched continuum observable. Gattringer and Lang 2010, chs. 4–6 give a standard treatment of the lattice spectral, operator, and renormalization ingredients connected here.

The observable chain and its independent checks

Section titled “The observable chain and its independent checks”
StageDefinition and finite-regulator dataEstimator or matching stepCorrelations and cutoff controlIndependent check and stop condition
Action and lattice perturbation theoryAction, measure, Fourier convention, vertex normalization, Brillouin zone, and symmetry factorNamed loop coefficient in a stated lattice schemeCommon ensembles and input parameters retained through matchingFree rule or Ward identity; stop if a coefficient lacks its measure, domain, or normalization
Correlator and state isolationOperator quantum numbers, boundary rule, Cij(t)C_{ij}(t), and covarianceGEVP or correlated multi-state fit for EnE_n and Zi(n)Z_i^{(n)}Basis rank, time window, images, and excited states varied jointlyHermiticity, positivity where applicable, and an exact fixture; stop if the level label is unstable
Bare insertionConnected and disconnected contractions for ⟨m∣J0∣n⟩\langle m\lvert J_0\rvert n\ranglePlateau, summation, and multi-state estimators with stochastic-noise treatmentSource–sink separations, shared configurations, and stochastic replicas remain correlatedExact small system and estimator agreement; stop if disconnected bias or excited states are unresolved
Renormalization and mixingComplete operator basis and JR(μ)=Z(μ,a)J0J_R(\mu)=Z(\mu,a)J_0 in a named schemeNonperturbative condition, step scaling, and continuum conversionMatrix-valued mixing and scale dependence propagated with the bare resultWard identity, scheme round trip, and step-scaling closure; stop if the basis or scale window is incomplete
Scale settingDimensionless reference quantity and declared physical inputCorrelated conversion from lattice unitsTuning, scale, and predicted observables keep their shared covarianceSecond reference and held-out round trip; stop if the claimed prediction was also an input
Improvement and continuum limitSymanzik-allowed action and operator basis, line of constant physics, QR(a,L)Q_R(a,L)Correlated finite-volume and continuum fit with a declared artifact ansatzMatching, renormalization, scale, volume, action choice, and fit alternatives propagated onceAlternative action, artifact model, and leave-one-spacing-out tests; stop if constant physics or the asymptotic window is unsupported
Reported continuum quantityQR(0,∞)Q_R(0,\infty) with scheme, scale, units, and covarianceRounding from the combined uncertainty, not from fit precision aloneStatistical and systematic components retain correlations and avoid double countingHeld-out observable or cross-formulation comparison; report only digits supported by the complete uncertainty

The repeated fields make omissions visible as the calculation moves from a regulated action to a continuum observable.

The dependency-free Chapter 2 benchmark reproduces eight bounded checks used across the chapter: the periodic free-scalar correlator and midpoint estimator, GEVP rank versus scalar time-signal identifiability, exact stochastic-trace enumeration, a spectral null direction and smeared broad peak, correlated scale propagation, the two-dimensional lattice tadpole, noncommuting mixing matrices, and continuum-fit rank and conditioning failures. Run it from the repository root with node scripts/benchmark-lattice-observables-continuum-inference.mjs; all eight groups must report PASS.

The validated 25 August 2026 source has SHA-256 d76924d8c724cb7c9307fb38265336495fa3ff1a98103b3fa4fa734262647c2a. The linked repository path may advance, so compare that digest before treating later output as the fixture documented here. The checks test the declared algebra and numerical constructions; they are not physical measurements or evidence for lattice actions outside their stated fixtures.

The chapter overview has no hard prerequisite. For the first four pages, you should be able to manipulate transfer spectral sums and covariance matrices. Repair the first capability with Euclidean Correlators and Schwinger Functions and Spectral Decomposition of Two-Point Functions. Repair lattice normalization with Lattice Momentum, Propagators, and Cutoff Dispersion.

For the final four pages, you should already distinguish bare and renormalized operators, scheme from observable, and a tuned trajectory from a bare-coupling scan. The repair routes are Local and Composite Operator Insertions, Nonperturbative Renormalization Schemes and Step Scaling, and Bare Parameters, Tuning Conditions, and Continuum Targets.

Check. Show that the two-state effective energy from C(t)=Ae−E0t+Be−E1tC(t)=Ae^{-E_0t}+Be^{-E_1t} approaches E0E_0 with a correction proportional to e−(E1−E0)te^{-(E_1-E_0)t} when A,B>0A,B>0.

Derive. Starting from the three-point spectral sum, identify the leading source-side and sink-side excited-state corrections to a ground-state matrix element.

Compute. For a 2×22\times2 mixing matrix ZZ, propagate the covariance of a bare vector v0v_0 and of the uncertain entries of ZZ into vR=Zv0v_R=Zv_0.

Diagnose. A study uses three lattice spacings, tunes the mass separately at each spacing, sets the scale with the same mass, and quotes that mass as its principal prediction. Explain which part is circular and propose a held-out observable.

Synthesize. Design a full calculation record for one continuum matrix element: operators, state isolation, insertion estimator, renormalization, scale, volume study, constant-physics conditions, artifact models, and two independent validation routes.

Solutions and checkpoints

Writing C(t)=Ae−E0t[1+re−Δt]C(t)=Ae^{-E_0t}[1+r e^{-\Delta t}] gives

Eeff(t)=log⁡C(t)C(t+a)1a=E0+1alog⁡1+re−Δt1+re−Δ(t+a),E_{\mathrm{eff}}(t)=\log\frac{C(t)}{C(t+a)}\frac1a =E_0+\frac1a\log\frac{1+r e^{-\Delta t}} {1+r e^{-\Delta(t+a)}},

whose leading correction is r(1−e−aΔ)e−Δt/ar(1-e^{-a\Delta})e^{-\Delta t}/a.

For the three-point function, terms (m,n)=(1,0)(m,n)=(1,0) and (0,1)(0,1) scale as e−Δm(Ts−τ)e^{-\Delta_m(T_s-\tau)} and e−Δnτe^{-\Delta_n\tau} relative to the ground term. Both distances must be varied.

For covariance propagation, linearize vR=Zv0v_R=Zv_0. The Jacobians are ∂vR,i/∂v0,j=Zij\partial v_{R,i}/\partial v_{0,j}=Z_{ij} and ∂vR,i/∂Zkl=δikv0,l\partial v_{R,i}/\partial Z_{kl}=\delta_{ik}v_{0,l}. Apply them to the joint covariance, retaining cross-covariances if ZZ and v0v_0 share ensembles.

The mass agreement is guaranteed because it tunes and sets the scale. A different mass ratio, decay constant, or matrix element not used in either step is a valid held-out prediction.

For the synthesis, a complete record names the action and tuned trajectory; source and sink operators; correlator normalization and covariance; state-isolation model and competing windows; every connected and disconnected contraction; the bare insertion and its estimator tests; the complete mixing basis, scheme, scale, and matching matrix; the independent scale input; the volume sequence; and the correlated continuum models. One validation route should close a regulated calculation internally—for example, an exact small-system correlator or Ward identity. A second should test the final inference externally—for example, a held-out observable, an alternative discretization with the same continuum target, or a scheme round trip. The record must also state a stop condition for unresolved state contamination, incomplete mixing, scale circularity, volume effects, or unsupported cutoff power counting.

After completing the chapter, you should be able to trace a lattice result from measured correlation functions to a continuum observable and identify exactly where state, model, renormalization, scale, or cutoff assumptions enter. Lattice Gauge Theory supplies gauge-specific actions and observables; Sampling Algorithms for Lattice Fields and Statistical Inference and Error Budgets supply the ensemble-generation and full statistical machinery used to execute the chain.

  • Gattringer, Christof, and Christian B. Lang. Quantum Chromodynamics on the Lattice: An Introductory Presentation. Springer, 2010, chs. 4–6. doi:10.1007/978-3-642-01850-3.
  • OpenAI Codex for QFT.org. “Lattice Observables and Continuum Inference Benchmark.” JavaScript source, validated 25 August 2026. SHA-256 d76924d8c724cb7c9307fb38265336495fa3ff1a98103b3fa4fa734262647c2a. JavaScript source on GitHub.
  • Beane, Silas R., William Detmold, Kostas Orginos, and Martin J. Savage. “Nuclear Physics from Lattice QCD.” Progress in Particle and Nuclear Physics 66, no. 1 (2011): 1–40. doi:10.1016/j.ppnp.2010.08.002.
  • Lüscher, Martin. “Advanced Lattice QCD.” In Les Houches 1997: Probing the Standard Model of Particle Interactions, 1998, pp. 229–280. arXiv:hep-lat/9802029.
  • Montvay, István, and Gernot Münster. Quantum Fields on a Lattice. Cambridge University Press, 1994. doi:10.1017/CBO9780511470783.

Original QFT.org content:CC BY 4.0, unless an item supplies different terms. Third-party material retains its own terms.