---
title: Hopping Parameter Expansion in Lattice QCD
url: https://www.emergentmind.com/topics/hopping-parameter-expansion-hpe
type: topic
---

# Hopping Parameter Expansion in Lattice QCD

Searching arXiv for the cited HPE literature to ground the article in current records.
{"query":"hopping parameter expansion lattice QCD heavy quark finite density Wilson fermions arXiv", "max_results": 10}
{"query":"2311.11508 hopping parameter expansion", "max_results": 5}
Hopping parameter expansion (HPE) is an expansion in powers of a lattice hopping parameter around a limit in which propagation between neighboring sites is suppressed. In lattice QCD with Wilson fermions, it is the Taylor expansion of the fermion determinant, or of \(\ln \det M\), about the heavy-quark point \(\kappa=0\), where \(\kappa\) is small for heavy quarks. In that setting, each order of the expansion is a sum over closed lattice trajectories, so the determinant is reorganized into Wilson-loop and Polyakov-loop sectors [2112.06340]. Closely related expansions also appear in strong-coupling treatments of bosonic lattice models, transfer-matrix analyses of gauge theories with matter, and linked-cluster formulations of lattice field theory [1210.2322].

## 1. Formal definition and loop interpretation

For Wilson fermions, the lattice Dirac operator is commonly written as
\[
M=1-\kappa Q
\]
or equivalently
\[
M_{xy}(\kappa)=\delta_{xy}-\kappa B_{xy},
\]
with \(\kappa\) related naively to the bare quark mass by
\[
\kappa=\frac{1}{2m+8}
\]
in one standard convention, and by
\[
\kappa_f=\frac{1}{2am_f+8}
\]
when the lattice spacing is written explicitly [1408.3770], [2211.08631]. The HPE is then
\[
\log \det M=\mathrm{Tr}\log(1-\kappa Q)
=-\sum_{n=1}^{\infty}\frac{\kappa^n}{n}\,\mathrm{Tr}\,Q^n,
\]
or, in the Wilson-kernel notation,
\[
\ln \det M(\kappa_f)=-\sum_{n=1}^{\infty}\frac{\kappa_f^n}{n}\,\mathrm{Tr}[B^n].
\]
Because \(Q\) or \(B\) is a nearest-neighbor hopping matrix, \(\mathrm{Tr}\,Q^n\) and \(\mathrm{Tr}[B^n]\) are nonzero only for closed paths of length \(n\). The order in \(\kappa\) is therefore the loop length, and the determinant becomes a weighted sum of gauge-invariant closed trajectories [1503.08813], [2112.06340].

This loop interpretation is the structural reason HPE is useful. In heavy-quark lattice QCD, short contractible loops renormalize gauge-action-like terms, while temporally winding loops couple directly to thermal observables such as the Polyakov loop. The lowest nontrivial nonwinding term is the plaquette contribution at order \(\kappa^4\), whereas the leading thermal term is the straight Polyakov loop at order \(\kappa^{N_t}\) [2401.00400]. For example, in one normalization,
\[
W(4)\kappa^4=288\,\kappa^4\,\hat P,
\qquad
L(N_t,N_t)\kappa^{N_t}=\frac{12}{N_t}2^{N_t}\kappa^{N_t}\,\mathrm{Re}\,\hat\Omega,
\]
so already at leading order the expansion produces a shift of the plaquette coupling and a linear source for the Polyakov loop [2211.08631].

The same logic appears outside QCD. In the spin-1 Bose–Hubbard model, the strong-coupling expansion is an expansion in the hopping amplitude \(t\) around the atomic limit \(t=0\), with the Mott state as the unperturbed vacuum and defect energies computed order by order in \(t\) [1210.2322]. In \(1+1\)-dimensional gauge theory with fundamental matter, the transfer matrix is expanded in powers of
\[
K=\frac{1}{2+(ma)^2},
\]
so that powers of \(K\) count matter-hopping processes in the heavy-mass regime [1705.01549].

## 2. Temporal winding, chemical potential, and heavy-quark thermodynamics

At finite density, HPE acquires a particularly transparent structure because the chemical potential enters only through temporal hoppings. In the Wilson-fermion matrix, the temporal forward and backward terms carry factors \(e^{\pm \mu a}\), so a closed path with temporal winding number \(m\) acquires \(e^{\pm m\mu/T}\) [2311.11508]. The expansion coefficients can therefore be decomposed as
\[
D_n=\hat W(n)+\sum_{m=1}^\infty \hat L_m^+(N_t,n)e^{m\mu/T}
+\sum_{m=1}^\infty \hat L_m^-(N_t,n)e^{-m\mu/T},
\]
or equivalently into \(\cosh(m\mu/T)\) and \(\sinh(m\mu/T)\) pieces, with the latter generating the complex phase of the determinant [2311.11508].

This decomposition is central in heavy-quark finite-temperature QCD. The leading winding contribution is the singly wound Polyakov-loop term, while higher winding sectors are parametrically smaller in the regime studied around the heavy-quark critical point [2112.06340]. In thermodynamic applications, the determinant can be reorganized by net winding number \(w\), yielding a grand potential of the form
\[
-\frac{\Omega_{\rm q}(T,\mu_{\rm q})}{T}
=
X_0+\sum_{w=1}^{\infty}
\left(e^{w\hat\mu_{\rm q}}+e^{-w\hat\mu_{\rm q}}\right)X_w,
\]
with \(\hat\mu_{\rm q}=\mu_{\rm q}/T\) and \(X_w\) built from connected Yang–Mills expectation values of Polyakov-loop operators [2508.09927].

A notable consequence is the leading-order prediction for quark-number susceptibility ratios in heavy-quark QCD. In the deconfined phase, the leading sector is \(w=1\), and the paper reports
\[
\frac{\chi_4}{\chi_2}=1
\qquad (T>T_{\rm c},\ \text{LO}),
\]
whereas in the confined phase, center symmetry removes \(w=1\) and \(w=2\), so the leading sector is \(w=3\), giving
\[
\frac{\chi_4}{\chi_2}=9
\qquad (T<T_{\rm c},\ \text{LO})
\]
[2508.09927]. The same HPE analysis yields an analytic formula for the quark excitation energy in the deconfined phase,
\[
E
=
m_{\rm pole}
-
T\ln\left[
\frac1{N_{\rm c}N_x^3}
\sum_{\mathbf s}
\langle \operatorname{tr}_{\rm c}\ell_{\rm P}(\mathbf s)\rangle
\right],
\]
and a baryonic decomposition in the confined phase in terms of winding-three Polyakov-loop cumulants [2508.09927].

## 3. Effective theories for the heavy-quark critical point

In finite-temperature heavy-quark QCD, HPE is not merely a formal determinant expansion; it is the basis of an effective theory in plaquette and Polyakov-loop variables. Keeping the leading HPE terms gives
\[
S_{g+\rm LO}
=
-6N_{\rm site}\beta^*\hat P
-\lambda N_s^3\,\mathrm{Re}\,\hat\Omega,
\]
with
\[
\beta^*=\beta+48N_{\rm f}\kappa^4,
\qquad
\lambda=2^{N_t+1}N_{\rm c}N_{\rm f}\kappa^{N_t}
\]
in the formulation used for two-flavor heavy-quark QCD [2401.00400]. This converts the fermion determinant into an effective external field for the deconfinement order parameter, and it enables efficient pseudo-heat-bath and over-relaxation updates instead of full dynamical-fermion simulation [2401.00400].

That strategy has been used to study the endpoint where the first-order deconfining transition turns into a crossover. In the \(N_t=4\) and \(N_t=6\) heavy-quark studies, leading-order terms are included in the simulation measure, next-to-leading-order terms are incorporated by reweighting, and finite-size scaling of the Binder cumulant of \(\mathrm{Re}\,\hat\Omega\) is used to extract the critical point [2211.08631], [2401.00400]. For \(N_f=2\), one cited NLO determination is
\[
\kappa_{\rm c,NLO}=0.0602(4)\quad (N_t=4),
\qquad
\kappa_{\rm c,NLO}=0.09003(19)\quad (N_t=6)
\]
[2311.11508].

A key HPE result at finite density is that the critical line can be expressed through an effective Polyakov-loop coupling \(\lambda^*\). For \(2+1\) flavors,
\[
2\sum_{n=N_t}^{n_{\max}}
L^0(N_t,n)\cosh\!\left(\frac{\mu}{T}\right)c_n\kappa_{\rm c,ud}^n
+
\sum_{n=N_t}^{n_{\max}}
L^0(N_t,n)\cosh\!\left(\frac{\mu}{T}\right)c_n\kappa_{\rm c,s}^n
=
\frac{\lambda_c^*}{N_t},
\]
so the same critical \(\lambda_c^*\) corresponds to smaller \(\kappa_c\) as \(\mu\) increases [2311.11508]. The reported physical conclusion is that the first-order phase-transition region in the heavy-quark region becomes narrower exponentially with increasing chemical potential, and that the critical hopping parameter decreases roughly exponentially with \(\mu/T\) [2311.11508]. In this regime, the sign problem remains mild because the critical \(\kappa_c\) moves deeper into the heavy-quark region, and the phase fluctuation is controlled by a small effective coupling [2311.11508].

## 4. Higher orders, convergence, and algorithmic reformulations

A recurrent theme in the HPE literature is the distinction between formal convergence and practical truncation error. A detailed worst-case analysis with all gauge links set to unity computes HPE coefficients to more than \(100\)th order and finds that the expansion converges up to
\[
\kappa=\frac18,
\]
the free Wilson-fermion chiral-limit value [2112.06340]. This result addresses a common objection: the issue is not that HPE has an exceptionally small convergence radius, but that the order required for accurate truncation grows with \(N_t\) and with the critical \(\kappa_c\) [2112.06340].

The same study shows that higher-order Polyakov-type loop terms are strongly correlated with the ordinary Polyakov loop on a configuration-by-configuration basis,
\[
L(N_t,n)\approx L^0(N_t,n)c_n\,\mathrm{Re}\,\hat\Omega,
\]
and reports that this linear relation holds up to the \(20\)th order for the measured ensembles [2112.06340]. Finite-density work extends this observation to the complex phase,
\[
\mathrm{Arg}\,\hat L_1^+(N_t,n)\approx \mathrm{Arg}\,\hat\Omega,
\]
which permits a “high-order from low-order” strategy in which higher-order effects are absorbed into an effective Polyakov-loop coupling rather than evaluated operator by operator [2311.11508]. For \(N_t=6\), the critical line is reported to stabilize for \(n_{\max}\gtrsim 10\), while explicit calculations were carried up to \(n=22\) [2311.11508].

These empirical correlations underpin the assessment of truncation errors in critical-point studies. The heavy-quark phase-structure analysis states that LO is fairly accurate for \(N_t=4\), NLO is fairly accurate for \(N_t=6\), and still higher orders are needed for larger \(N_t\) [2211.08631]. On finer lattices, the \(N_t=6\) Binder-cumulant study reports
\[
\lambda_c=0.000818(10),
\qquad
\kappa_c^{\rm NLO}=0.09003(19),
\qquad
\kappa_c=0.08781(17)
\]
after incorporating yet higher-order HPE contributions, showing that higher-order corrections remain quantitatively relevant even when LO+NLO captures the basic physics [2401.00400].

A separate line of work reformulates HPE to all orders while preserving the full Yang–Mills action exactly. One formulation is the direct \(\kappa\)-expansion,
\[
\det(1-\kappa Q)=
\exp\left[
-\sum_{n=1}^{\infty}\frac{\kappa^n}{n}\mathrm{Tr}\,Q^n
\right],
\]
and a second is the \(\kappa_s\)-expansion,
\[
\det M=
\det(1-R)\,
\exp\left[
-\sum_{n=1}^{\infty}
\frac{\kappa_s^n}{n}
\mathrm{Tr}\left(\frac{1}{1-R}S\right)^n
\right],
\]
which treats temporal hopping and its \(\mu\)-dependence analytically and exactly while expanding only spatial hopping [1408.3770], [1503.08813]. Simulations with complex Langevin show that at \(\mu=0.7\) both expansions approach full-QCD results at sufficiently high order, while at \(\mu=1.1\) the direct \(\kappa\)-expansion breaks down and the \(\kappa_s\)-expansion still converges, already agreeing well with full QCD by about order \(\kappa^{10}\) [1408.3770].

On the computational side, HPE has also been used as a UV filter in Hybrid Monte Carlo. For two degenerate Wilson fermions, the leading nonvanishing contribution \(\mathrm{Tr}\,H^4\) is a plaquette term and gives
\[
\Delta\beta=96\kappa^4,
\]
while numerical tests report a speed-up of roughly a factor \(2\) for \(\kappa^2\)-filtering and \(3\) for \(\kappa^4\)-filtering [1805.03560]. More recently, explicit high-order evaluation of \(\kappa^8\), \(\kappa^{10}\), and \(\kappa^{12}\) terms in \(\mathrm{Tr}\ln M\) has been made practical with trie-based algorithms, with computational costs of approximately \(20\), \(460\), and \(8900\) times that of a single staple evaluation, respectively [2606.31492].

## 5. Uses beyond heavy-quark lattice QCD

Although HPE is most closely associated with heavy-quark lattice QCD, the same expansion principle appears in other gauge and many-body systems. In the slave-fermion \(t\)-\(J\) model for underdoped cuprates, fermionic holons are integrated out by a hopping expansion in the holon hopping amplitude. The effective small parameter is stated to be \(t\times \delta\), and the leading induced term is
\[
\tilde A_{\rm hop}
=
\delta\left(\frac{c_3}{2}\right)^2
\sum_{x,\mu}
|\bar z_{x+\mu}z_x|^2
+
\delta\left(\frac{c_3}{2}\right)^4
\sum_{x,\mu<\nu}\prod_{\rm plaq.}(\bar z z)+\cdots.
\]
This generates an effective bosonic lattice gauge theory with emergent link fields \(U_{x\mu}\), \(V_{x\mu}\), and \(M^\star_{x\mu}\), used to analyze antiferromagnetic, metal–insulator, and superconducting transitions [1007.4273].

In \(1+1\)-dimensional \(SU(N)\) lattice gauge theory with matter, HPE is performed on the transfer matrix in powers of
\[
K=\frac{1}{2+(ma)^2}.
\]
The ground state is found to contain local mesons at \(O(K^2)\), nearest-neighbor mesons crossing a cut at \(O(K^3)\), and longer strings at higher orders [1705.01549]. The entanglement analysis based on that HPE shows that Shannon-sector entropy and color entanglement appear at \(O(K^3)\), while the Bell-pair contribution first appears at \(O(K^6)\) in the wavefunction, giving an entropy contribution
\[
S_{EE}^{\rm Bell}
=
\left(1-\log K^{12}\right)K^{12}
+\mathcal{O}(K^{14})
\]
[1705.01549].

In the spin-1 Bose–Hubbard model with antiferromagnetic interaction, the strong-coupling expansion in the hopping amplitude \(t\) is carried through third order to determine Mott-state and defect energies. The paper concludes that the Mott insulator phase is considerably more stable against the superfluid phase when filling with an even number of bosons than when filling with an odd number of bosons, reflecting the role of on-site singlet formation [1210.2322].

At a more formal level, the Functional Renormalization Group has been used to reinterpret the hopping expansion as a linked-cluster expansion for the Legendre effective action,
\[
\Gamma_{\kappa}[\phi]
=
\Gamma_0[\phi]
+
\sum_{l\ge 2}\kappa^l\Gamma_l[\phi],
\]
with the critical hopping parameter identified with the finite radius of convergence of the susceptibilities and with the unstable manifold of a Gaussian or non-Gaussian fixed point of the FRG flow [1812.02251]. This suggests a direct bridge between HPE, linked-cluster graph rules, and nonperturbative renormalization-group resummations.

## 6. Scope, limitations, and recurring misconceptions

A persistent misconception is that HPE is only a low-order heavy-quark approximation with no systematic extension. The all-orders \(\kappa\)- and \(\kappa_s\)-expansions, the explicit high-order convergence studies, and the algorithmic evaluation of \(\kappa^8\), \(\kappa^{10}\), and \(\kappa^{12}\) terms show that the method can be systematically improved to high order [1408.3770], [2112.06340], [2606.31492]. A second misconception is that the existence of many higher-order loops automatically destroys predictivity; the heavy-quark critical-point literature instead reports strong configuration-by-configuration correlations between higher-order Polyakov-type terms and the ordinary Polyakov loop, which permit an effective one-coupling description over the regime studied [2311.11508].

The principal limitation is unchanged across applications: HPE is controlled only when the hopping parameter is sufficiently small. In lattice QCD, this means the heavy-quark region; once \(\kappa\) is no longer small, neglected higher-order and non-Polyakov-loop structures become important [2311.11508]. The required truncation order also grows with \(N_t\), so low-order results that are reliable at \(N_t=4\) do not automatically remain reliable on finer lattices [2211.08631]. In finite density, higher winding sectors \(m>1\) can eventually matter because of the \(e^{m\mu/T}\) enhancement, even though they are numerically small in the parameter regions emphasized in the heavy-quark critical-point studies [2112.06340].

Another recurring issue concerns the sign problem. The finite-density heavy-quark studies do not claim that HPE removes the sign problem universally; rather, they report that along the heavy-quark critical line the sign problem does not become serious even when the density increases, because \(\kappa_c\) decreases exponentially and the phase fluctuation remains small [2311.11508]. This is a regime statement, not a general theorem.

Taken together, the literature presents HPE as a family of controlled expansions around a static or atomic limit, organized by closed trajectories in the hopping matrix. In heavy-quark QCD it reduces the fermion determinant to an effective plaquette–Polyakov-loop theory, makes the \(\mu\)-dependence explicit through temporal winding, and supports critical-point, thermodynamic, and algorithmic analyses [2311.11508], [2508.09927]. In broader lattice field theory and condensed-matter contexts, it serves the same structural purpose: integrating out mobile degrees of freedom in powers of their hopping and replacing them with effective local or gauge-invariant operators [1007.4273], [1705.01549], [1210.2322].

Source: https://www.emergentmind.com/topics/hopping-parameter-expansion-hpe