Hubbard–Stratonovich Fields and Collective Channels

A Hubbard–Stratonovich field is an exact Gaussian rewriting of an interaction, not a new approximation. It can expose density, spin, pairing, or order-parameter fluctuations and make selected resummations transparent. Approximation begins only when the auxiliary-field integral is truncated, sampled with a constraint, or evaluated about a saddle.

Required background. Coherent-State Path Integrals for Many-Body Systems fixes the regulated functional measure, and Gaussian Fields and Sources supplies Gaussian completion of the square.

Helpful background. Degrees of Freedom, Symmetry, and the Local Operator Expansion helps distinguish a true emergent mode from a convenient auxiliary variable; Emergent Variables and Reorganized Effective Descriptions supplies the broader reorganization principle.

The Gaussian identity

For a real positive constant $g$ and a real variable $O$,

$$e^{\frac g2O^2} =\frac1{\sqrt{2\pi g}} \int_{-\infty}^{\infty}\mathrm d\sigma\, e^{-\frac{\sigma^2}{2g}+\sigma O}.$$

Completing the square, $-\sigma^2/(2g)+\sigma O=-(\sigma-gO)^2/(2g)+gO^2/2$, proves the identity. Functional versions apply it at each spacetime point, with a regulator making the product finite.

The sign matters. If the quadratic form has the opposite sign, a convergent representation may require an imaginary coupling or a rotated contour. Writing a formally real auxiliary field without stating the contour can turn an identity into a divergent integral.

For a complex pairing operator $P=\psi_\downarrow\psi_\uparrow$ and $g>0$,

$$e^{gP^\dagger P} \propto \int\mathrm d\Delta^*\mathrm d\Delta\, \exp\!\left[-\frac{|\Delta|^2}{g} +\Delta^*P+P^\dagger\Delta\right].$$

Integrating out $\Delta$ recovers the original attraction. This round trip fixes every sign and normalization.

Density, spin, and pairing channels

Fermionic anticommutation permits algebraic rearrangements. For a one-band site,

$$n_\uparrow n_\downarrow =\frac14n^2-\frac14(n_\uparrow-n_\downarrow)^2, \qquad n=n_\uparrow+n_\downarrow.$$

One may therefore decouple the Hubbard interaction in charge or spin channels. An attractive interaction can instead be decoupled in the pairing channel. If the auxiliary integral is performed exactly, all valid decouplings give the same partition function and observables.

The apparent choice becomes consequential after truncation. A saddle in the pairing channel privileges anomalous order; a density saddle privileges Hartree structure. Fierz-equivalent decompositions need not yield identical mean-field theories because the omitted fluctuations differ. This is channel ambiguity, not evidence that the exact identities disagree.

Hirsch’s discrete transformation provides a useful lattice example. For a repulsive Hubbard time slice,

$$e^{-\epsilon U(n_\uparrow-\frac12)(n_\downarrow-\frac12)} =\frac12e^{-\epsilon U/4} \sum_{s=\pm1}e^{\lambda s(n_\uparrow-n_\downarrow)}, \qquad \cosh\lambda=e^{\epsilon U/2}.$$

Because $n_\sigma$ has eigenvalues zero or one, the identity is checked on four local occupation states. It is exact for the slice; Trotter error enters elsewhere Hirsch 1983, pp. 4059–4061.

From auxiliary field to collective propagator

After decoupling a fermion interaction, integrate out the fermions formally:

$$Z=\int\mathcal D\sigma\,e^{-S_{\rm eff}[\sigma]}, \qquad S_{\rm eff}[\sigma] =\frac12\sigma g^{-1}\sigma -\operatorname{Tr}\ln G^{-1}[\sigma].$$

A stationary field $\bar\sigma$ satisfies $\delta S_{\rm eff}/\delta\sigma=0$. Expanding $\sigma=\bar\sigma+\delta\sigma$ gives

$$S_{\rm eff}=S_{\rm eff}[\bar\sigma] +\frac12\delta\sigma\,D^{-1}\,\delta\sigma+\cdots, \qquad D^{-1}=g^{-1}-\Pi,$$

where $\Pi$ is a channel-dependent polarization kernel. Zeros of $D^{-1}$ can signal collective modes or instabilities, but only after analytic continuation, symmetry, and damping are checked.

The saddle is controlled when a parameter suppresses fluctuations—for example large flavor number, weak coupling in a suitable regime, high coordination, or proximity to an upper critical dimension. Merely obtaining a sharp stationary point numerically is not a proof of control.

Exact rewrite versus physical degree of freedom

An auxiliary field becomes an efficient low-energy degree of freedom when its propagator develops a pole or long correlation length below the cutoff and its couplings admit a controlled expansion. Before that, it is a redundant integration variable. Rescaling $\sigma\mapsto A\sigma$ changes its residue and Yukawa coupling without changing observables, so an auxiliary-field expectation value is not automatically normalized as a measured order parameter.

For gauge-charged pairing fields, $\langle\Delta\rangle$ is gauge dependent. Physical claims require gauge-invariant response, stiffness, or a properly fixed and interpreted formulation. The Gaussian identity alone establishes none of these.

Common pitfalls

Calling the transformation an approximation. The regulated Gaussian identity is exact. The saddle, loop truncation, contour deformation, or sampling constraint is the approximation.

Ignoring the contour. Repulsive and attractive channels require different signs. A divergent real Gaussian cannot be justified by notation.

Selecting a channel after seeing the answer. Equivalent exact decouplings can produce inequivalent truncations. Compare channels or justify the control parameter before making a phase claim.

Exercises

Integrate out the pairing field

Complete the square in the complex Gaussian and recover $e^{gP^\dagger P}$.

Solution

Write

$$-\frac{|\Delta|^2}{g}+\Delta^*P+P^\dagger\Delta =-\frac{|\Delta-gP|^2}{g}+gP^\dagger P.$$

The shifted Gaussian contributes only its normalization, leaving $e^{gP^\dagger P}$. This also fixes the relative signs of the Yukawa terms.

Verify the discrete Hubbard identity

Check the four occupation states and derive $\cosh\lambda=e^{\epsilon U/2}$.

Solution

For $(n_\uparrow,n_\downarrow)=(0,0)$ or $(1,1)$, the left side is $e^{-\epsilon U/4}$ and the sum gives the same value. For $(1,0)$ or $(0,1)$, the left side is $e^{+\epsilon U/4}$, while the right side is $e^{-\epsilon U/4}\cosh\lambda$. Equality requires $\cosh\lambda=e^{\epsilon U/2}$.

Continue

Nonrelativistic Power Counting and Universality determines whether an auxiliary-field expansion is controlled. Irreducible Vertices and Bethe–Salpeter Equations gives the two-particle kernel whose poles can reappear as collective auxiliary-field modes. Two-Channel Resonance Models shows when a molecular field carries physical low-energy resonance data.

References

• Hirsch, J. E. “Discrete Hubbard–Stratonovich Transformation for Fermion Lattice Models.” Physical Review B 28 (1983): 4059–4061. DOI.

Further reading

• Hubbard, John. “Calculation of Partition Functions.” Physical Review Letters 3 (1959): 77–78. DOI.
• Stratonovich, Ruslan L. “On a Method of Calculating Quantum Distribution Functions.” Soviet Physics Doklady 2 (1958): 416–419. Bibliographic record.