---
title: Classical Holstein-Spin-Fermion Model
url: https://www.emergentmind.com/topics/classical-holstein-spin-fermion-model
type: topic
---

# Classical Holstein-Spin-Fermion Model

Searching arXiv for the specified papers and closely related Holstein / spin-fermion formulations.
Searching for `2509.05876 classical Holstein spin fermion model`.
The classical Holstein-spin-fermion model is a sign-problem-free classical auxiliary-field formulation that combines a Holstein-type lattice coupling with a spin-fermion interaction, while retaining quantum itinerant electrons. In the formulation used to study its phase transition, the model is obtained from a Hubbard-Holstein starting point by a mean-field approximation, and its defining feature is that both the lattice displacement field and the local spin field are treated classically. In the adiabatic limit, the lattice field is static rather than a quantum phonon operator, and the local spins are classical unit vectors. The resulting model is designed to retain competing microscopic interactions—specifically antiferromagnetic and charge-order tendencies—while making large-scale Monte Carlo sampling feasible [2509.05876].

## 1. Definition and Hamiltonian

The starting point is a Hubbard-Holstein-type Hamiltonian written as
\[
H=H_0+H_{\mathrm{e\text{-}ph}}+H_{\mathrm{e\text{-}e}}+H_{\mathrm{lattice}},
\]
with
\[
H_0=-t\sum_{\langle i, j \rangle,\sigma} c_{i,\sigma}^{\dagger} c_{j,\sigma}^{\phantom\dagger} - \mu\sum_{i,\sigma} \hat{n}_{i,\sigma},
\]
\[
H_{\mathrm{e\text{-}ph}}=g\sum_{i,\sigma} \hat{x}_i \left(\hat{n}_{i,\sigma}-\frac{1}{2}\right),
\]
\[
H_{\mathrm{e\text{-}e}}=U\sum_{i} \hat{n}_{i,\uparrow} \hat{n}_{i,\downarrow},
\]
and
\[
H_{\mathrm{lattice}}=\sum_i \left[\frac{1}{2M}\left(\frac{d\hat{x}_i}{dt}\right)^2 + \frac{K^2}{2} \hat{x}_i^2\right].
\]
Here \(t\) is the nearest-neighbor hopping, \(\mu\) the chemical potential, \(g\) the electron-phonon coupling strength, \(U\) the onsite Coulomb repulsion, \(M\) the lattice mass, \(K\) the elastic constant, \(\hat{x}_i\) the lattice displacement operator, and \(\hat{n}_{i,\sigma}\) the fermion density operator. In the reported simulations, \(t=1\) and \(K=1\), and \(\mu\) is adjusted to keep the electron density at 1 [2509.05876].

To eliminate the sign problem of the full interacting fermion model, the mean-field approximation replaces the Hubbard interaction by a spin-fermion coupling and yields the classical Holstein-spin-fermion Hamiltonian
\[
H=H_0+H_{\mathrm{e\text{-}ph}}+H_{\mathrm{sf}}+H_{\mathrm{lattice}},
\]
with
\[
H_{\mathrm{sf}}=J\sum_i \mathbf{S}_i \cdot \hat{\mathbf{s}}_i,
\]
where \(\mathbf{S}_i\) is a classical local spin vector with fixed amplitude 1, \(\hat{\mathbf{s}}_i\) is the quantum electron spin operator, and \(J\) is the spin-fermion coupling. The benchmark study sets \(J=2\) [2509.05876].

In the adiabatic limit \(M\to\infty\), the lattice kinetic term drops out and the lattice sector becomes
\[
H_{\mathrm{lattice}}=\sum_i \frac{K^2}{2}x_i^2,
\]
with \(x_i\) treated as a classical variable. The Monte Carlo auxiliary fields are then
\[
u=\{x_i, S_{i,x}, S_{i,y}, S_{i,z}\}.
\]
This formulation isolates a hybrid problem: electrons remain quantum mechanical, but the spin and lattice backgrounds are classical [2509.05876].

## 2. Classical character and relation to the Holstein family

The model is “classical” in two distinct senses. First, the lattice field is classical: in the adiabatic limit, \(x_i\) is not a dynamical quantum phonon operator but a classical displacement. Second, the local spins are classical: \(\mathbf{S}_i\) are classical unit vectors rather than quantum spin operators. This distinguishes the model from the full Hubbard-Holstein problem, where both the electron-electron and electron-phonon sectors are quantum, and from more conventional spin-fermion models without an explicit Holstein lattice sector [2509.05876].

Its Holstein component is inherited from the standard local density-displacement structure. In the quantum Holstein model, the Hamiltonian is written as
\[
\hat{H} = \hat{H}_0 +\sum_i \left( \frac{1}{2M} P_i^2 + \frac{K}{2} Q_i^2 \right) + g \sum_i Q_i \hat{\rho}_i ,
\]
with
\[
\hat{H}_0 = -t \sum_{i\sigma} \left( c^\dagger_{i\sigma} c_{i+1\sigma} + c^\dagger_{i+1\sigma} c_{i\sigma} \right),
\qquad
\hat{\rho}_i= \sum_{\sigma} \left( \hat{n}_{i\sigma} - \frac{1}{2} \right).
\]
After integrating out the phonons, the interaction becomes retarded and purely fermionic. That construction is explicitly quantum and retains full imaginary-time phonon dynamics, so it is not a classical Holstein-spin-fermion model in the static-background sense [1609.02026].

A related but distinct generalization appears in trapped-ion realizations of spin-Holstein physics. There, an effective spin model coupled to phonons can be mapped by a Jordan–Wigner transformation to a fermionic Holstein-like problem. The analogy is direct—\(\sigma_i^z\) becomes fermion density, spin-flip terms become hopping, and spin-phonon coupling becomes electron-phonon coupling—but the effective model contains long-range \(1/|i-j|^3\) exchange, dispersive collective phonons, and nonlocal couplings \(g_{in}\). It is therefore described as a generalized Holstein model rather than the textbook local one [2108.13730].

## 3. Competing interactions and phase structure

The central physical content of the classical Holstein-spin-fermion model is the competition between spin-driven antiferromagnetism and lattice-driven charge ordering. The spin-spin or AFM-favoring tendencies originate from the spin-fermion term, while the electron-lattice coupling favors a charge-density-wave state. This competition drives an AFM-to-CDW transition [2509.05876].

At \(\beta=10\), the reported numerical regime is divided as follows: AFM fluctuations occur for \(g<1.5\), while CDW fluctuations occur for \(1.5\le g<1.8\). Accordingly, \(g<1.5\) is referred to as the AFM region and \(g>1.5\) as the CDW region [2509.05876]. The microscopic interpretation is explicit. In the AFM region, the electron-phonon term slightly raises the energy because it induces spin fluctuations away from the AFM ground state. In the CDW region, the coupling to the displacement field strongly lowers the energy by stabilizing a checkerboard charge modulation. The CDW picture is described in terms of atoms displaced by \(\pm a\), such that the density becomes 0 on one sublattice and 2 on the other, lowering the energy by
\[
N g a.
\]
This places the transition in the broader class of interacting-order competition problems [2509.05876].

The broader Holstein literature provides a useful comparison. In one-dimensional half-filled quantum Holstein models, the Peierls transition is accompanied by phonon softening at \(q=\pi\), and in the adiabatic regime this is consistent with a soft-mode transition. The ordered state is a Peierls insulating state with long-range charge-density-wave order at half filling [1609.02026]. A plausible implication is that the classical Holstein-spin-fermion model retains the charge-ordering logic of the Holstein problem while replacing quantum phonon dynamics by static distortions and adding an explicit spin background.

The same charge-order tendency also appears in multicomponent generalizations. In the SU(\(N\)) Holstein model at half filling, the low-temperature phase is an insulating CDW in which empty sites alternate with sites with \(N\) particles; for \(N=3\), the pattern is alternating empty and triply occupied sites. That ordered phase is treated as breaking a discrete \(\mathbb Z_2\) symmetry and is analyzed using 2D Ising finite-size scaling [2606.25240]. This suggests that classical or semiclassical Holstein-type charge order remains robust when the fermionic sector is generalized beyond the usual spin-\(1/2\) case.

## 4. Monte Carlo formulation and self-learning acceleration

For a proposed field update \(u\to u'\), the Monte Carlo acceptance probability is
\[
p=\min(1,p_u),\qquad p_u=e^{-\beta F(u^\prime)+\beta F(u)}.
\]
The free energy is written as
\[
F(u)=-\sum_i\ln(1+e^{-\beta E_i})/\beta+E_{\mathrm{lattice}},
\]
with
\[
E_{\mathrm{lattice}}=\sum_i \frac{K^2}{2}x_i^2 - g\sum_i x_i.
\]
Here \(E_i\) are the electronic eigenenergies for a given auxiliary-field configuration, obtained by exact diagonalization with LAPACK’s `zgeev` routine. This exact-diagonalization step is the computational bottleneck: it scales as \(\mathcal{O}(N^3)\), and because it must be repeated through the update sequence, a standard QMC epoch scales as \(\mathcal{O}(N^4)\) [2509.05876].

The self-learning quantum Monte Carlo strategy replaces the exact free-energy evaluation by a machine-learning surrogate \(F_{\mathrm{ML}}(u)\). Trial updates are accepted with
\[
p_{\mathrm{trial}}=\min(1,e^{-\beta F_{\mathrm{ML}}(u^\prime)+\beta F_{\mathrm{ML}}(u)}).
\]
After a full epoch, a cumulative correction restores consistency through
\[
p_c=\min\!\left(1, e^{-\beta (F(u^\prime)-F(u))} e^{\beta \sum_i \Delta F_i}\right),
\]
where \(\Delta F_i\) are the predicted free-energy changes accumulated over accepted trial updates. This reduces the per-epoch cost to \(\mathcal{O}(N^3)\) [2509.05876].

Two surrogate classes are compared. The linear-regression model takes the form
\[
F_{\mathrm{LR}}= E_{\mathrm{s\text{-}s}}+E_{\mathrm{l\text{-}l}}+E_{\mathrm{s\text{-}l}}+W(x)+c_0,
\]
where the effective interactions are included up to the 7th nearest neighbor and the model has 24 terms excluding the bias. The onsite term is
\[
W(x)=\sum_i c_1\left[\frac{x_{{\bf r}_i}^4}{4}-\frac{\overline{x^2}x_{{\bf r}_i}^2}{2}\right]+c_2 x_{{\bf r}_i}^2+c_3|x_{{\bf r}_i}|.
\]
The neural-network model supplements this with two-point correlation functions and the linear-regression feature set. At \(g=1.5\) and \(N=10\), the neural network reduces the MSE to about \(0.014\), compared with \(0.038\) for linear regression [2509.05876].

The simulations also employ annealing, cooling the system gradually from \(T=t\) to \(T=0.1t\) with decay rate 0.9 and 2000 epochs per temperature. In this setting, both linear regression and neural-network surrogates reproduce the AFM-to-CDW transition, while the neural network provides improved free-energy prediction and better cumulative acceptance [2509.05876].

## 5. Diagnostics, sampling efficiency, and finite-size effects

Phase identification is performed with magnetic and charge susceptibilities at momentum point \(M\),
\[
\chi_m({\bf q})=\frac{1}{N} \sum_{i,{\bf d}} \langle \hat{\mathbf s}_{{\bf r}_i}\hat{\mathbf s}_{{\bf r}_i+{\bf d}} \rangle e^{i{\bf q}\cdot {\bf d}},
\]
and
\[
\chi_c({\bf q})=\frac{1}{N} \sum_{i,{\bf d}} \langle \hat{n}_{{\bf r}_i}\hat{n}_{{\bf r}_i+{\bf d}} \rangle e^{i{\bf q}\cdot {\bf d}}.
\]
Using these observables, the self-learning results from both surrogate classes agree well with standard QMC and recover AFM fluctuations at small \(g\) and CDW fluctuations at large \(g\) [2509.05876].

The detailed energy decomposition from the linear-regression model clarifies the transition. The total free energy changes smoothly across the crossover region; spin-spin interaction energies decrease rapidly in the AFM region; lattice-lattice terms become important and drop in the CDW region; and spin-lattice terms are smaller and change sign across the transition. A specific reported feature is that the second-nearest-neighbor spin-spin interaction vanishes in the CDW region, while nearest-neighbor spin correlations remain finite. The CDW state also suppresses some longer-range lattice correlations [2509.05876].

The principal numerical limitation is a degradation of sampling efficiency near the AFM-CDW transition and on larger lattices. The cumulative update probability \(p_c\) is reasonably good in the AFM and CDW regimes far from criticality, but near \(g=1.5\) it becomes much smaller. For the linear-regression model at \(g=1.5\), \(p_c\) decreases from about 0.6 to 0.1 as \(N_x\) grows from 10 to 30. The neural network improves this behavior, but \(p_c\) still declines with size [2509.05876].

The stated interpretation links Monte Carlo efficiency directly to surrogate accuracy. Near the AFM-CDW transition, simultaneous spin and lattice fluctuations increase the MSE and RMSE of the free-energy predictor. Since the acceptance depends on \(\beta F\), the effective control parameter is \(\beta\times\mathrm{RMSE}\). The finite-size effect is described in terms of the decreasing energy gap between the ground state and low-lying excited states: when that gap becomes comparable to or smaller than the machine-learning error, the sampled state is no longer close to a pure ground state but becomes a mixture of ground and excited states. The qualitative criterion reported is
\[
\mathrm{RMSE}<\frac{1}{N^\alpha},
\]
for some model-dependent exponent \(\alpha\) [2509.05876].

A related Monte Carlo issue appears in determinant QMC for the Holstein model, where the configuration space develops extremely long autocorrelation times despite being sign-problem free. A Wang–Landau flat-histogram sampling in configuration-weight space was introduced to reduce those autocorrelation times, with the key tradeoff that lower-probability configurations are sampled more often at some loss of sampling efficiency per step [2107.14454]. This suggests that the classical Holstein-spin-fermion model belongs to a broader class of Holstein-type systems where the dominant obstacle is mixing time rather than the sign problem itself.

## 6. Related formulations and generalizations

Several neighboring models clarify what is specific to the classical Holstein-spin-fermion construction. The trapped-ion spin-Holstein model provides a generalized Holstein simulator in which effective spins, phonons, and laser-induced couplings emulate electron-phonon physics. After Jordan–Wigner transformation, it becomes a fermionic Holstein-like model with long-range hopping or exchange, dispersive phonons, and nonlocal couplings. Its phase diagram contains CDW order in the stiff regime, superconducting or pairing correlations in the soft regime, phase separation at strong coupling, and a mixed pCDW phase at quarter filling [2108.13730]. These ingredients are absent from the classical Holstein-spin-fermion benchmark, whose purpose is instead to isolate competing AFM and CDW tendencies in a sign-problem-free setting.

A different generalization is the spin-dependent electron-transfer model motivated by radical-pair chemistry. There, the Holstein electron-transfer framework is extended by explicit electron spin, magnetic-field orientation, and hyperfine coupling to nuclear spins. The total Hamiltonian is \(H_{\rm tot}=H_r+H_s\), and the key physical conclusion is that the triplet-state reaction rate depends on the direction of the magnetic field, whereas the singlet-state rate does not [1110.1918]. This is a spinful Holstein generalization, but it is not a lattice spin-fermion model and does not involve classical local-spin backgrounds.

The SU(\(N\)) Holstein model extends the fermionic sector from two species to \(N\) flavors coupled identically to a single local phonon mode per site. At half filling, defined by \(\rho=N/2\), it exhibits an insulating CDW phase in which empty sites alternate with fully occupied sites containing all \(N\) fermions. For \(N=3\), the reported critical temperature can be as high as twice the maximum attainable for \(N=2\), and the large-\(N\) trend approaches \(T_c\to 0.63\,t\) at fixed representative parameters [2606.25240]. This broader family emphasizes that the charge-order tendency of Holstein systems can be strengthened by increasing the number of fermionic components.

The quantum Holstein literature also contributes formal tools relevant to classical or semiclassical formulations. In continuous-time interaction-expansion Monte Carlo, bosonic observables can be reconstructed from fermionic correlation functions using generating functionals, and improved estimators can be built directly from the vertex distribution for total energy, phonon propagators, and fidelity susceptibility [1609.02026]. A plausible implication is that classical Holstein-spin-fermion studies can draw on this formalism for diagnostics even when the bosonic sector is treated statically rather than quantum mechanically.

Source: https://www.emergentmind.com/topics/classical-holstein-spin-fermion-model