Skip to content

Fluctuation Operators and Determinant Ratios

The one-loop prefactor around a nontrivial saddle is the Gaussian measure of its nonzero fluctuations, normalized to the same measure in a reference sector and renormalized with the same boundary conditions and subtraction scheme. It is therefore a determinant ratio, not an isolated formal product of eigenvalues. This page derives that ratio, states when the Gel’fand–Yaglom method computes it, and evaluates the fluctuation operator of the quartic-double-well instanton.

Required background. Saddles, control parameters, and loop counting supplies the rescaling that isolates the Hessian and the criterion selecting a contributing saddle.

Helpful background. Heat kernels, zeta functions, and spectral determinants supplies regulator-independent definitions of determinant ratios and their ultraviolet expansion; second variation, Hessians, and Jacobi fields supplies the boundary-value problem defined by the quadratic action.

From the second variation to a determinant ratio

Section titled “From the second variation to a determinant ratio”

Let ϕσ\phi_\sigma be a Euclidean saddle of a regulated action S/g\mathcal S/g. For a bosonic fluctuation η\eta obeying the linearized boundary conditions,

S[ϕσ+g η]=Sσ+g2⟨η,Mση⟩+O(g3/2),\mathcal S[\phi_\sigma+\sqrt g\,\eta] =\mathcal S_\sigma +\frac g2\langle\eta,M_\sigma\eta\rangle +O(g^{3/2}),

where

Mσ(x,y)=δ2Sδϕ(x)δϕ(y)∣ϕσ.M_\sigma(x,y) =\left.\frac{\delta^2\mathcal S} {\delta\phi(x)\delta\phi(y)}\right|_{\phi_\sigma}.

In a finite spectral regulator with strictly positive eigenvalues λnσ\lambda_n^\sigma, the Gaussian integral is proportional to

∏n(λnσ)−1/2=(det⁡Mσ)−1/2.\prod_n(\lambda_n^\sigma)^{-1/2} =(\det M_\sigma)^{-1/2}.

The normalization of the functional measure cancels only in a ratio. Relative to a reference saddle with operator MrefM_{\rm ref},

Aσ(1)=(det⁡Mσdet⁡Mref)−1/2.\mathcal A_\sigma^{(1)} =\left(\frac{\det M_\sigma}{\det M_{\rm ref}}\right)^{-1/2}.

If MσM_\sigma has kk zero modes, they must be omitted and replaced by collective-coordinate integrals:

Aσ(1)=(det⁡′Mσdet⁡Mref)−1/2.\mathcal A_\sigma^{(1)} =\left(\frac{\det{}' M_\sigma}{\det M_{\rm ref}}\right)^{-1/2}.

The prime has operational content: it specifies a projection, reduced determinant, or limiting prescription. It must not mean “discard whichever eigenvalues make the answer inconvenient.” Negative eigenvalues likewise cannot be hidden inside an absolute value; their phase is fixed by the integration cycle.

For a positive elliptic operator with eigenvalues of mass dimension two, zeta regularization introduces a reference scale μ\mu and defines

ζM(s;μ)=∑n(λnμ2)−s,log⁡det⁡ ⁣(Mμ2)=−ζM′(0;μ).\zeta_M(s;\mu)=\sum_n\left(\frac{\lambda_n}{\mu^2}\right)^{-s}, \qquad \log\det\!\left(\frac{M}{\mu^2}\right)=-\zeta_M'(0;\mu).

Equivalently,

log⁡det⁡Mσdet⁡Mref=−∫0∞dtt Tr⁡ ⁣(e−tMσ−e−tMref),\log\frac{\det M_\sigma}{\det M_{\rm ref}} =-\int_0^\infty\frac{\mathrm dt}{t}\, \operatorname{Tr}\!\left(e^{-tM_\sigma}-e^{-tM_{\rm ref}}\right),

after the required small-tt subtractions. In quantum mechanics the difference is often finite directly. In QFT the small-tt terms are local ultraviolet divergences and must be canceled by counterterms fixed in the reference vacuum; the residual μ\mu dependence of the regulated determinant cancels against those counterterms and the running parameters through the retained order. A determinant value without its regulator, boundary conditions, zero-mode prescription, and renormalization scale is not a reproducible prefactor. Dunne 2008, §§2–4, printed pp. 3–8 reviews these complementary definitions and their domain of validity.

For scalar Sturm–Liouville operators on [a,b][a,b],

M=−d2dτ2+U(τ),M0=−d2dτ2+U0(τ),M=-\frac{d^2}{d\tau^2}+U(\tau), \qquad M_0=-\frac{d^2}{d\tau^2}+U_0(\tau),

with the same Dirichlet boundary conditions, solve

My=0,y(a)=0,y′(a)=1,My=0,\quad y(a)=0,\quad y'(a)=1,

and the analogous problem for y0y_0. When neither operator has a zero mode,

det⁡Mdet⁡M0=y(b)y0(b).\frac{\det M}{\det M_0}=\frac{y(b)}{y_0(b)}.

This formula is powerful because it replaces an infinite spectral product by an initial-value problem. Its hypotheses are equally important:

  • the operators have the same leading symbol and compatible boundary conditions;
  • the comparison uses the same interval and normalization;
  • a zero mode is treated by a reduced-determinant limit;
  • matrix-valued problems use a fundamental solution matrix and its determinant;
  • ultraviolet counterterms are still required when transverse momenta or higher-dimensional fields are present.

One reduced-determinant prescription is to shift M↦M+λM\mapsto M+\lambda, compute the ratio, and remove the simple zero:

det⁡′M=lim⁡λ→0det⁡(M+λ)λ.\det{}'M =\lim_{\lambda\to0} \frac{\det(M+\lambda)}{\lambda}.

On a finite interval an infinite-volume translation zero mode is usually lifted to a small eigenvalue. One must project the known mode or take the correlated large-interval and zero-eigenvalue limit; simply deleting the numerically smallest eigenvalue can remove a physical soft mode instead. Mariño 2015, §1.6, pp. 30–35 derives the reduced Gel’fand–Yaglom construction for instanton operators.

Pöschl–Teller fluctuations of the double well

Section titled “Pöschl–Teller fluctuations of the double well”

For the action

SE[x]=1g∫dτ [12x˙2+12(x2−1)2]S_E[x]=\frac1g\int \mathrm d\tau\, \left[\frac12\dot x^2+\frac12(x^2-1)^2\right]

and instanton xI(τ)=tanh⁡(τ−τ0)x_I(\tau)=\tanh(\tau-\tau_0), the Hessian is

MI=−d2dτ2+4−6 sech⁡2(τ−τ0).M_I=-\frac{d^2}{d\tau^2} +4-6\,\operatorname{sech}^2(\tau-\tau_0).

The vacuum operator at x=±1x=\pm1 is

M0=−d2dτ2+4.M_0=-\frac{d^2}{d\tau^2}+4.

This is the reflectionless l=2l=2 Pöschl–Teller problem. Its discrete spectrum consists of

λ0=0,λ1=3,\lambda_0=0,\qquad \lambda_1=3,

with continuum λ(k)=k2+4\lambda(k)=k^2+4. The normalized zero mode is proportional to sech⁡2(τ−τ0)\operatorname{sech}^2(\tau-\tau_0); there is no negative eigenvalue.

The spectral calculation can be made explicit. For the reflectionless family

Mℓ,m=−d2dτ2+m2−ℓ(ℓ+1)sech⁡2τ,M_{\ell,m}=-\frac{\mathrm d^2}{\mathrm d\tau^2} +m^2-\ell(\ell+1)\operatorname{sech}^2\tau,

the Jost factors give

det⁡Mℓ,mdet⁡M0,m=∏j=1ℓm−jm+j.\frac{\det M_{\ell,m}}{\det M_{0,m}} =\prod_{j=1}^{\ell}\frac{m-j}{m+j}.

At m=ℓ=2m=\ell=2, the factor m−2m-2 marks the translation eigenvalue. The actual eigenvalue in this family is m2−4=(m−2)(m+2)m^2-4=(m-2)(m+2), so dividing by it before taking m→2m\to2 gives

det⁡′MIdet⁡M0=lim⁡m→21m2−4∏j=12m−jm+j=22−12(2+1)2(2+2)2=148.\frac{\det{}'M_I}{\det M_0} =\lim_{m\to2} \frac{1}{m^2-4} \prod_{j=1}^{2}\frac{m-j}{m+j} =\frac{2^2-1^2}{(2+1)^2(2+2)^2} =\frac1{48}.

An initial-value calculation reaches the same number without multiplying the spectrum. Put the operators on [−L,L][-L,L] with Dirichlet boundaries. For MIyI=0M_Iy_I=0, yI(−L)=0y_I(-L)=0, and yI′(−L)=1y_I'(-L)=1, reduction of order from u0=sech⁡2τu_0=\operatorname{sech}^2\tau gives

yI(L)=2sech⁡4L I(L),I(L)=∫0Lcosh⁡4t dt=3L8+sinh⁡2L4+sinh⁡4L32.y_I(L)=2\operatorname{sech}^4L\,I(L), \qquad I(L)=\int_0^L\cosh^4t\,\mathrm dt =\frac{3L}{8}+\frac{\sinh2L}{4}+\frac{\sinh4L}{32}.

The reference solution is y0(L)=sinh⁡(4L)/2y_0(L)=\sinh(4L)/2, so Gel’fand–Yaglom gives

rL≡det⁡MI,Ldet⁡M0,L=4sech⁡4L I(L)sinh⁡4L=2e−4L[1+o(1)].r_L\equiv\frac{\det M_{I,L}}{\det M_{0,L}} =\frac{4\operatorname{sech}^4L\,I(L)}{\sinh4L} =2e^{-4L}\bigl[1+o(1)\bigr].

The coefficient of the lifted mode can also be read from the even shooting solution. Write u=u0+λu1+O(λ2)u=u_0+\lambda u_1+O(\lambda^2) with u(0)=1u(0)=1, u′(0)=0u'(0)=0, and MIu1=u0M_Iu_1=u_0. Variation of parameters with the second zero-energy solution v=u0Iv=u_0I gives

u1(L)=−e2L24[1+o(1)],u0(L)=4e−2L[1+o(1)].u_1(L)=-\frac{e^{2L}}{24}\bigl[1+o(1)\bigr], \qquad u_0(L)=4e^{-2L}\bigl[1+o(1)\bigr].

The Dirichlet condition u(L)=0u(L)=0 therefore lifts the translation mode to λ0(L)=96e−4L[1+o(1)]\lambda_0(L)=96e^{-4L}[1+o(1)]. Projecting precisely that mode yields

lim⁡L→∞rLλ0(L)=148,\lim_{L\to\infty}\frac{r_L}{\lambda_0(L)}=\frac1{48},

in agreement with the spectral calculation. This comparison also shows why deleting an arbitrary small numerical eigenvalue is not an acceptable reduced-determinant prescription.

Before nondimensionalization, det⁡′MI/det⁡M0\det{}'M_I/\det M_0 has dimensions of time squared because one inverse-time-squared eigenvalue has been removed. Its inverse square root supplies one inverse time. In the normalization above the Jacobian J=(SI/2πg)1/2J=(\mathcal S_I/2\pi g)^{1/2} is dimensionless, while J dτ0J\,\mathrm d\tau_0 carries one time; after extracting ∫dτ0=T\int\mathrm d\tau_0=T, the remaining fugacity has dimensions of inverse time. The value 1/481/48 is tied to the displayed operator normalization; rescaling τ\tau rescales a primed determinant differently from an unprimed one. Mariño 2015, §§1.5–1.8, pp. 26–42 derives both evaluations.

The complete one-instanton factor is not just 48\sqrt{48}. It has the form

κ=(SI2πg)1/2(det⁡′MIdet⁡M0)−1/2e−SI/g[1+O(g)],SI=43,\kappa =\left(\frac{\mathcal S_I}{2\pi g}\right)^{1/2} \left(\frac{\det{}'M_I}{\det M_0}\right)^{-1/2} e^{-\mathcal S_I/g} \bigl[1+O(g)\bigr], \qquad \mathcal S_I=\frac43,

up to the stated normalization of the Euclidean time kernel. The first factor is the dimensionless translation Jacobian, the second supplies the inverse-time dimension of the event fugacity, and the exponential is the classical weight. The separate integral over τ0\tau_0 produces the total-time factor. Their separation prevents a zero eigenvalue from being counted twice.

Shared calculation. The saddle-contribution anatomy shows where this determinant sits relative to the contour phase, collective measure, counterterms, and observable. Shared comparison. The canonical comparison table records the corresponding prescription and failure boundary for several standard saddles.

In a field theory, write the one-loop contribution schematically as

Γσ(1)=12log⁡det⁡′Mσ−12log⁡det⁡Mref+Sct(1)[ϕσ]−Sct(1)[ϕref].\Gamma_\sigma^{(1)} =\frac12\log\det{}'M_\sigma -\frac12\log\det M_{\rm ref} +S_{\rm ct}^{(1)}[\phi_\sigma] -S_{\rm ct}^{(1)}[\phi_{\rm ref}].

The same renormalization conditions must fix both sectors. A successful calculation passes four checks:

  1. all ultraviolet regulator dependence cancels through the claimed order;
  2. zero and negative mode counts agree between spectral and variational analyses;
  3. the result has the dimension required by the moduli measure and observable;
  4. an independent method, such as phase shifts versus Gel’fand–Yaglom, agrees.

The ordinary determinant expansion is not uniform when a nonzero eigenvalue becomes comparable to the interaction corrections. Near such a soft mode, integrate that coordinate non-Gaussianly or use a uniform approximation. It also fails if the reference operator uses different boundary data, if an infinite-volume limit is taken before projecting the zero mode, or if a QFT determinant is quoted before counterterm subtraction.

Taking the determinant of a differential expression without a domain. Boundary conditions are part of the operator. Changing Dirichlet to periodic data changes both the spectrum and the determinant.

Deleting the smallest numerical eigenvalue automatically. Finite-volume translation modes are small but nonzero, and other physical modes may also soften. Project the analytically identified zero mode and test stability as the volume changes.

Using an absolute determinant for a negative mode. This destroys the contour phase that distinguishes a decay saddle from a stable one. Continue the Gaussian along the prescribed steepest direction instead.

  1. Derive the double-well fluctuation potential 4−6 sech⁡2τ4-6\,\operatorname{sech}^2\tau.
Solution

For V(x)=12(x2−1)2V(x)=\tfrac12(x^2-1)^2,

V′′(x)=6x2−2.V''(x)=6x^2-2.

Substituting xI=tanh⁡τx_I=\tanh\tau and using tanh⁡2τ=1−sech⁡2τ\tanh^2\tau=1-\operatorname{sech}^2\tau gives

V′′(xI)=6(1−sech⁡2τ)−2=4−6 sech⁡2τ.V''(x_I)=6(1-\operatorname{sech}^2\tau)-2 =4-6\,\operatorname{sech}^2\tau.
  1. Why does det⁡′MI/det⁡M0\det{}'M_I/\det M_0 acquire dimensions even if the unprimed ratio is dimensionless?
Solution

Both unprimed determinants contain the same number of eigenvalues, so their dimensions cancel in a common regulator. Removing one eigenvalue from the numerator removes one factor of inverse time squared. The primed ratio therefore has dimensions of time squared. Its inverse square root has dimensions of inverse time and combines with the collective-coordinate integral to give the required dimension.

  1. Let Mϵ=M+ϵM_\epsilon=M+\epsilon and suppose MM has one normalized zero mode and otherwise positive spectrum. Show that det⁡Mϵ=ϵ det⁡′M+O(ϵ2)\det M_\epsilon=\epsilon\,\det{}'M+O(\epsilon^2).
Solution

In a spectral regulator,

det⁡Mϵ=ϵ∏n≥1(λn+ϵ)=ϵ(∏n≥1λn)[1+O(ϵ)].\det M_\epsilon =\epsilon\prod_{n\geq1}(\lambda_n+\epsilon) =\epsilon\left(\prod_{n\geq1}\lambda_n\right) \bigl[1+O(\epsilon)\bigr].

The product over nonzero eigenvalues is det⁡′M\det{}'M, which proves the stated limit. The same conclusion follows after a consistent zeta or determinant-ratio regularization.

Continue from spectra to the remaining saddle data

Section titled “Continue from spectra to the remaining saddle data”
  • Dunne, Gerald V. “Functional Determinants in Quantum Field Theory.” Journal of Physics A: Mathematical and Theoretical 41 (2008): 304006. DOI.
  • Mariño, Marcos. Instantons and Large N: An Introduction to Non-Perturbative Methods in Quantum Field Theory. Cambridge University Press, 2015, ch. 1, pp. 3–61. DOI.

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