---
title: BUQEYE Model for EFT Uncertainty Quantification
url: https://www.emergentmind.com/topics/buqeye-model
type: topic
---

# BUQEYE Model for EFT Uncertainty Quantification

Searching arXiv for BUQEYE-related papers to ground the article in the latest research.
BUQEYE, expanded by the collaboration as **Bayesian Uncertainty Quantification: Errors in Your EFT**, denotes a Bayesian framework for uncertainty quantification in effective field theory (EFT), especially for correlated truncation errors in low-energy nuclear physics, and more broadly a research program that also includes projection-based emulators for fast surrogate modeling [2212.04912]. In its truncation-error formulation, the central hypothesis is that dimensionless coefficient functions extracted from order-by-order EFT calculations can be treated as independent draws from a common Gaussian process (GP), enabling lower-order results to inform predictive distributions for omitted higher-order terms [2402.13165]. Subsequent studies have used this model to assess convergence patterns across modern nucleon-nucleon potentials, infer breakdown scales such as $\Lambda_b$, and construct statistically calibrated uncertainty bands for observables including Wolfenstein amplitudes [2508.17558; 2501.09231].

## 1. Definition and scope

In the nuclear-physics literature represented here, the BUQEYE model is a **correlated EFT truncation-error model** in which order-by-order observable calculations are recast in terms of dimensionless coefficients and those coefficients are modeled probabilistically by a GP [2402.13165]. The collaboration’s guide to projection-based emulators places this statistical work within a larger BUQEYE framework concerned with Bayesian uncertainty quantification and computational acceleration for many-body and scattering calculations [2212.04912].

The core application is the EFT expansion of an observable. One writes the prediction at truncation order $k$ as
$$
y_k(x)=y_{\rm ref}(x)\sum_{n=0}^k c_n(x)\,Q^n(x),
$$
with $y_{\rm ref}(x)$ a reference scale, $c_n(x)$ dimensionless coefficients, and $Q(x)$ the expansion parameter [2402.13165]. The truncation error is then
$$
\delta y_k(x)=y_{\rm ref}(x)\sum_{n=k+1}^\infty c_n(x)\,Q^n(x),
$$
so the statistical problem reduces to learning the size and correlation structure of the unknown higher-order coefficients from the known lower-order ones [2402.13165].

A related but distinct BUQEYE usage appears in the literature on Wolfenstein amplitudes, where the same GP machinery is applied to angle-correlated truncation uncertainties for the five independent on-shell nucleon-nucleon amplitude functions [2501.09231]. This suggests that “BUQEYE model” is best understood not as a single narrow algorithm, but as a family of closely related Bayesian GP constructions for EFT-order uncertainty propagation.

## 2. Bayesian hypothesis and mathematical structure

The defining BUQEYE hypothesis is that once the observable has been normalized by a reference scale and the powers of the EFT expansion parameter have been factored out, the coefficients behave like independent draws from the same GP:
$$
c_n(x)\overset{\rm iid}{\sim}\GP\bigl(\mu(x)=0,\;\sigma^2\,r(x,x';L)\bigr).
$$
In the formulation used for nucleon-nucleon observables, the correlation kernel is the squared-exponential kernel,
$$
K(x,x')=\sigma^2\exp\!\bigl[-\tfrac12\,(x-x')^T\,L^{-1}\,(x-x')\bigr],
$$
where $\sigma^2$ is the marginal variance and $L$ is a diagonal matrix of correlation-length scales such as $\ell_E$ in energy and $\ell_\theta$ in angle [2402.13165].

The same structure appears in the Wolfenstein-amplitude analysis, specialized to a one-dimensional input variable $q$:
$$
c_n(q)\sim\mathcal{GP}\bigl(0,\;\bar c^2\,r(q,q';\ell_q)\bigr),
$$
with
$$
r(q,q';\ell_q)=\exp\!\bigl[-\,\tfrac{(q-q')^2}{2\,\ell_q^2}\bigr].
$$
There, the hyperparameters are the signal variance $\bar c^2$, the correlation length $\ell_q$, and the EFT breakdown scale $\Lambda_b$ entering $Q$ [2501.09231].

The observable expansion is paired with a specific expansion parameter. One common choice is
$$
Q(p,m_\pi;\Lambda_b)=\frac{p+m_{\rm eff}}{\Lambda_b+m_{\rm eff}},
$$
with $p$ a characteristic momentum, $m_{\rm eff}$ a soft scale, and $\Lambda_b$ the EFT breakdown scale [2402.13165]. In the Wolfenstein-amplitude study, the expansion parameter is instead written
$$
Q=\frac{\max(p,q)+m_\pi}{\Lambda_b+m_\pi},
$$
with the same role assigned to $\Lambda_b$ [2501.09231].

Because the omitted terms are sums of GP-distributed coefficient functions weighted by powers of $Q$, BUQEYE induces a GP—or, after marginalizing over variance parameters, a Student-$t$ process—for the truncation error itself [2508.17558]. In the 2024 assessment paper, the predictive covariance for the missing terms at order $k$ is
$$
K^{(k)}(x,x')
=y_{\rm ref}(x)\,y_{\rm ref}(x')\,
\frac{[\,Q(x)\,Q(x')\,]^{\,k+1}}{1-Q(x)\,Q(x')}\,
r(x,x';L),
$$
which directly encodes correlated uncertainty bands across kinematics [2402.13165].

## 3. Inference, diagnostics, and workflow

BUQEYE proceeds by extracting empirical coefficients from order-by-order calculations. Given predictions $\{y_n(x)\}_{n=0}^k$, one computes
$$
c_n(x)=
\frac{1}{y_{\rm ref}(x)\,Q^n(x)}
\Bigl[y_n(x)-y_{n-1}(x)\Bigr]
\quad(n\ge1),
$$
with $c_0(x)=y_0(x)/y_{\rm ref}(x)$ [2402.13165]. These extracted coefficient curves are then used to fit GP hyperparameters by conditioning the prior on a training set and scoring candidate length scales via the GP marginal likelihood, or its Student-$t$ analog after integrating out $\sigma^2$ [2402.13165].

The published workflow is explicit. It begins with data preparation on a dense grid, choice of $y_{\rm ref}(x)$ and $Q$, coefficient extraction, and a train–test split. Hyperparameter inference is followed by a stationarity check using validation diagnostics. If the diagnostics pass, one performs scale inference over $(\Lambda_b,m_{\rm eff},L)$, marginalizes to obtain posteriors for the physical scales, and then refits the GP to compute final truncation-error bands [2402.13165].

Three diagnostics recur across the BUQEYE literature. The first is a **Mahalanobis-distance diagnostic**, in which
$$
\Delta^2=[c_n(x_{\rm val})]^T\,K_{\rm val}^{-1}\,[c_n(x_{\rm val})]
$$
should follow $\chi^2_M$ under the GP model, with outliers indicating mis-sized variance or length scale [2402.13165]. The second is the **pivoted-Cholesky diagnostic**, in which
$$
D_{\rm PC}=G^{-1}[c_n(x_{\rm val})]
$$
should have standard normal components if the covariance model is correct; funneling or trumpeting patterns reveal under- or overestimated length scales [2402.13165]. The third is the **credible-interval coverage** or “weather” plot, which compares empirical coverage to nominal probability and is also used in later convergence assessments [2402.13165; 2508.17558].

A concise summary of the standard BUQEYE analysis pipeline is given below.

| Step | Operation | Purpose |
|---|---|---|
| 1 | Compute $y_n(x)$ on a grid | Assemble order-by-order EFT data |
| 2 | Choose $y_{\rm ref}(x)$ and $Q$ | Define dimensionless coefficients |
| 3 | Form $c_n(x)$ | Extract empirical convergence pattern |
| 4 | Split train/validation points | Enable GP fitting and diagnostics |
| 5 | Fit $(L,\sigma^2)$ | Learn smoothness and variance |
| 6 | Run Mahalanobis, PC, and coverage diagnostics | Test model adequacy |
| 7 | Infer $(\Lambda_b,m_{\rm eff},L)$ | Calibrate physical scales |
| 8 | Compute predictive truncation bands | Propagate correlated uncertainty |

This workflow is reproducible through publicly available Jupyter notebooks, a point emphasized in both the 2024 assessment paper and the 2025 multi-potential study [2402.13165; 2508.17558].

## 4. Convergence patterns, stationarity, and model adequacy

A major theme of the BUQEYE program is that the GP prior is not automatically valid for every EFT interaction or every parametrization. The 2024 assessment of semi-local momentum-space nucleon-nucleon potentials concluded that the BUQEYE model is “generally applicable” to the potential investigated, but also found that the assumption of GP stationarity across lab energy and scattering angle is “not generally met,” requiring adjustments in future work [2402.13165].

The later multi-potential study sharpened this point. It found that some “soft” regulator choices exhibit irregular convergence, specifically a mismatch between the sizes of coefficients at even and odd orders [2508.17558]. In that setting, a simple stationary GP prior across $(p,\theta)$ is inadequate. The same paper reported that the GP length scale in the angular dimension, $\ell_\theta$, is approximately inversely proportional to the relative momentum, so the correlation structure is nonstationary in angle even when momentum smoothness is roughly stationary [2508.17558].

The proposed remediation is either a momentum-dependent warp of the angular coordinate,
$$
\theta\mapsto\tilde\theta=\theta\;\bigl(\tfrac{p}{p_{\max}}\bigr)^B,
\quad B\approx1,
$$
or a nonstationary kernel in which $\ell_\theta$ becomes a function of $p$ [2508.17558]. After this modification, diagnostics show significant improvement, and the BUQEYE model is validated for the retained hard potentials [2508.17558].

The literature therefore treats BUQEYE not as a fixed black-box prior, but as a statistically testable model class. Good practice begins by plotting the coefficient functions themselves to inspect whether they have roughly the same amplitude and wiggliness across orders, which the 2024 paper describes as naturalness and stationarity [2402.13165]. A plausible implication is that BUQEYE’s strongest contributions are methodological as much as inferential: it supplies an explicit criterion for when EFT convergence patterns can be regarded as probabilistically regular and when they cannot.

## 5. Applications to nucleon-nucleon observables and Wolfenstein amplitudes

The BUQEYE model has been applied to multiple nucleon-nucleon datasets and observable classes. The 2024 assessment paper considered the semi-local momentum-space implementation of the chiral EFT expansion of the nucleon-nucleon potential and concluded that the model enables statistically principled estimates of the impact of higher EFT orders on observables [2402.13165]. It also demonstrated posterior inference for the expansion parameter $Q$, the breakdown scale $\Lambda_b$, the soft scale $m_{\rm eff}$, and GP hyperparameters [2402.13165].

The 2025 Wolfenstein-amplitude analysis extends the same GP machinery to the five independent operator coefficients of the on-shell nucleon-nucleon amplitude. There, all five amplitudes $A$, $C$, $M$, $G$, and $H$, with real and imaginary parts treated separately, are modeled by separate scalar GPs; there is no multi-output covariance in that work [2501.09231]. At laboratory energies $E_{\rm lab}\in\{25,50,100,150,200\}$ MeV, one-dimensional GPs are built in the input variable $q=2p\sin\frac{\theta}{2}$, using training and testing point sets across the $q$-range [2501.09231].

The study reports that, with the inferred breakdown scales, the EFT truncation uncertainties cover both higher-order results and empirical Wolfenstein amplitudes well for all orders other than the leading order [2501.09231]. It further states that typical MAP breakdown scales are $\Lambda_b\simeq750\text{–}800$ MeV and typical correlation lengths are $\ell_q\simeq120$ MeV [2501.09231].

The later cross-potential convergence study broadens the application set to six common nucleon-nucleon scattering observables as functions of relative momentum and scattering angle [2508.17558]. After excluding soft regulators with severe rms splitting and applying angular warping or equivalent nonstationary modeling, diagnostics pass for the remaining six potentials: SMS 450, 500, 550 MeV; SCS 0.9, 1.0 fm; and EMN 500 MeV [2508.17558].

## 6. Hyperparameter inference, breakdown scales, and broader BUQEYE methods

A central scientific use of the BUQEYE model is the inference of EFT scales from convergence data. Because the extracted coefficients depend on the choice of $Q$, and $Q$ depends on $(\Lambda_b,m_{\rm eff})$, these physical quantities become inferential targets rather than fixed external inputs [2402.13165]. The 2024 paper recommends computing a joint likelihood over $(\Lambda_b,m_{\rm eff},L)$, combining it with priors, and marginalizing over $L$ to obtain $\pi(\Lambda_b,m_{\rm eff})$ [2402.13165].

In the multi-potential analysis, with $m_{\rm eff}$ fixed at the physical pion mass, the preferred breakdown scales for various N$^{3}$LO interactions lie in the range $\Lambda_b=600$--$750$ MeV [2508.17558]. The same study reports that statistically consistent distributions for $\Lambda_b$ across orders are only found for the SMS potential at a regulator scale of 450 or 500 MeV [2508.17558]. When $\Lambda_b$ and $m_{\rm eff}$ are inferred simultaneously, joint MAPs lie in $\Lambda_b\approx650\mbox{–}800$ MeV and $m_{\rm eff}\approx60\mbox{–}200$ MeV depending on the potential, but the authors caution that posterior ellipses do not collapse as expected, indicating an identifiability tension between the two parameters [2508.17558].

The BUQEYE collaboration’s program also extends beyond truncation errors. Its guide to projection-based emulators defines emulators as fast surrogate models that replace expensive, high-fidelity nuclear-physics calculations by solving a much smaller reduced problem built from representative snapshots [2212.04912]. The reduction is based on variational and Galerkin projection methods together with offline–online decomposition. For bound states, inserting $\psi=Xc$ into the Rayleigh–Ritz functional yields the reduced generalized eigenproblem
$$
X^\dagger H(\theta)X\,c=E(\theta)\,X^\dagger X\,c,
$$
or, after orthonormalization, $\hat H(\theta)c=E(\theta)c$ with $\hat H(\theta)=X^\dagger H(\theta)X$ [2212.04912]. This broader BUQEYE activity is methodologically adjacent to the truncation-error model: both are aimed at making Bayesian uncertainty quantification operational in computational nuclear theory.

## 7. Significance, limitations, and recurring misconceptions

The principal significance of the BUQEYE model is that it converts EFT truncation uncertainty from an informal order-estimate into a statistically calibrated predictive distribution with explicit correlation structure across kinematic variables [2402.13165]. Because the same GP prior governs all orders, lower-order information constrains omitted higher-order contributions, and credible intervals can be checked against withheld calculations through Mahalanobis, pivoted-Cholesky, and coverage diagnostics [2402.13165; 2508.17558].

A common misconception is that BUQEYE assumes all EFT calculations automatically satisfy a universal stationary GP prior. The published assessments argue otherwise. They show that the model’s success depends on a well-chosen expansion parameter $Q$, on regulator choices that produce regular convergence patterns, and on stationarity assumptions that may need to be relaxed or remapped, especially in the angular variable [2402.13165; 2508.17558]. Another misconception is that the breakdown scale is simply fixed once and for all. In the BUQEYE framework, $\Lambda_b$ is often treated as a hyperparameter to be inferred from the coefficient data, with order dependence and posterior stability used as diagnostics of model adequacy [2501.09231; 2508.17558].

The available literature also delineates clear limitations. Leading-order uncertainty bands can be too small in some channels, as seen for forward-angle tensor and spin-orbit Wolfenstein amplitudes [2501.09231]. Soft regulator scales can generate even–odd coefficient splitting severe enough that the single-GP prior should not be used without modification [2508.17558]. And simultaneous inference of $\Lambda_b$ and $m_{\rm eff}$ can exhibit identifiability problems [2508.17558].

Within those limits, BUQEYE provides a unified statistical language for EFT convergence assessment, correlated truncation-error propagation, and, through related projection-based emulators, fast many-query computation in nuclear theory [2212.04912]. This suggests that its lasting contribution is not merely a set of uncertainty bands for specific observables, but an integrated methodology in which EFT expansion structure, Bayesian inference, and model diagnostics are treated as parts of the same quantitative framework.

Source: https://www.emergentmind.com/topics/buqeye-model