Skip to content

Derrick Scaling, Virial Tests, and Nontopological Stability

Derrick scaling tests whether a proposed static localized field can be stationary under a uniform change of size. If every positive contribution to the energy changes with the same sign, no stationary finite-size lump exists under the theorem’s hypotheses. An evasion must supply a term with the opposite scaling, impose a conserved-charge constraint, change the field content or geometry, or leave the static ansatz through genuine time dependence. Naming the stabilizer is part of the result.

Required background. Variational field equations and conserved currents supplies the stationary-action argument, and interactions, potentials, and stability supplies the bounded-energy assumptions. Helpful background. EFT power counting is needed when higher-derivative stabilization is used near a cutoff.

Consider nn real scalar fields on flat RD\mathbb R^D with static energy

E[ϕ]=E2+E0=dDx[12Gab(ϕ)iϕaiϕb+V(ϕ)],E[\phi]=E_2+E_0 =\int\mathrm d^D x \left[ \frac12 G_{ab}(\phi)\, \partial_i\phi^a\partial_i\phi^b+V(\phi) \right],

where GabG_{ab} is positive definite, V0V\geq0 after subtracting the vacuum energy, and the fields approach a vacuum fast enough for all integrals and integrations by parts to exist. Define the size variation

ϕλ(x)=ϕ(λx),λ>0.\phi_\lambda(\mathbf x)=\phi(\lambda\mathbf x), \qquad \lambda>0 .

Changing variables to y=λx\mathbf y=\lambda\mathbf x gives

E(λ)=λ2DE2+λDE0.E(\lambda) =\lambda^{2-D}E_2+\lambda^{-D}E_0 .

Any regular static solution is stationary under this admissible variation, so

dEdλλ=1=(2D)E2DE0=0.\left.\frac{\mathrm dE}{\mathrm d\lambda}\right|_{\lambda=1} =(2-D)E_2-DE_0=0 .

This is the virial identity. For D3D\geq3, both coefficients are negative and a nontrivial solution with E2,E00E_2,E_0\geq0 is impossible. In D=2D=2, stationarity requires E0=0E_0=0 and leaves the two-derivative energy scale invariant; size is not fixed. In D=1D=1, the relation E2=E0E_2=E_0 permits a kink. This is the core of Derrick 1964, pp. 1252–1254; Coleman 1985, § 2, pp. 194–195 gives the field-theory interpretation.

More generally, if a positive energy term contains a total of kk spatial derivatives and its field amplitudes are not rescaled, then

Ek[ϕλ]=λkDEk[ϕ].E_k[\phi_\lambda]=\lambda^{k-D}E_k[\phi].

For E=kEkE=\sum_k E_k, the scale residual is

RDdE(λ)dλ1=k(kD)Ek.\mathcal R_{\mathrm D} \equiv \left.\frac{\mathrm dE(\lambda)}{\mathrm d\lambda}\right|_1 =\sum_k(k-D)E_k .

A numerical profile claimed to be static should satisfy RD=0\mathcal R_{\mathrm D}=0 within a tolerance controlled separately from the field-equation residual. A small equation residual on a finite grid can coexist with a large virial error if the box, tails, or boundary data are wrong.

Hypotheses that must travel with the obstruction

Section titled “Hypotheses that must travel with the obstruction”

The conclusion is not “solitons do not exist above one spatial dimension.” It concerns a particular variational problem:

  • the configuration is time independent and localized on unbounded flat space;
  • the scale variation stays within the admissible field space and preserves the declared boundary sector;
  • the energy consists of the stated positive terms and has no explicit position-dependent scale;
  • the fields are regular enough that E(λ)E(\lambda) is differentiable at λ=1\lambda=1;
  • no conserved quantity is being held fixed by an additional Lagrange multiplier;
  • boundary, curvature, background fields, gauge constraints, and higher derivatives have not been omitted.

If any item fails, the scaling calculation must be redone rather than cited away. The theorem is a diagnostic: it says exactly what kind of new contribution is needed. Manton and Sutcliffe 2004, § 4.2, pp. 82–87 discuss the general scaling method and its model dependence.

Let U(x)SU(2)U(\mathbf x)\in SU(2) in three spatial dimensions and define Li=UiUL_i=U^\dagger\partial_iU. A standard static Skyrme energy has two- and four-derivative pieces,

E2=fπ216d3x[tr(LiLi)],E_2 =\frac{f_\pi^2}{16} \int\mathrm d^3x\, \big[-\operatorname{tr}(L_iL_i)\big], E4=132e2d3x[tr([Li,Lj][Li,Lj])],E_4 =\frac{1}{32e^2} \int\mathrm d^3x\, \big[-\operatorname{tr}([L_i,L_j][L_i,L_j])\big],

with positive integrands for anti-Hermitian LiL_i. A nonnegative potential contribution E0E_0 may also be present. Under Uλ(x)=U(λx)U_\lambda(\mathbf x)=U(\lambda\mathbf x),

E(λ)=λ1E2+λE4+λ3E0,E(\lambda)=\lambda^{-1}E_2+\lambda E_4+\lambda^{-3}E_0 ,

so a stationary solution must satisfy

E2+E43E0=0.-E_2+E_4-3E_0=0 .

Without E4E_4, shrinking lowers the two-derivative energy and no finite size is selected. The four-derivative term grows under shrinking and can balance E2E_2 and E0E_0. This is a real evasion because it changes the energy functional; topology alone did not do the work. The original higher-derivative construction is due to Skyrme 1961, pp. 127–138, while the modern scaling analysis is summarized in Manton and Sutcliffe 2004, §§ 4.2 and 9.1, pp. 82–87 and 349–356.

If the four-derivative operator belongs to an EFT with cutoff Λ\Lambda, the solution’s inverse size must remain parametrically below Λ\Lambda, and omitted operators must be smaller on the profile. A formal virial balance at R1ΛR^{-1}\sim\Lambda is not a controlled EFT prediction.

Gauge fields. Gauge potentials carry their own scaling required by covariance, and magnetic-field energy can balance Higgs-gradient and potential terms. The smooth monopole is the canonical example. One must scale the full gauge–Higgs ansatz; applying the scalar formula only to the Higgs profile gives the wrong conclusion.

Fixed Noether charge. A Q-ball has Φ(t,x)=eiωtf(x)\Phi(t,\mathbf x)=e^{i\omega t}f(\mathbf x). It is stationary only in energy density, not as a field. The relevant variational problem minimizes EE at fixed QQ, or equivalently extremizes EωQE-\omega Q. The time-dependent phase contributes a term absent from the static scalar theorem.

Periodic or rotating motion. Oscillons and rotating solitons can evade a static obstruction dynamically. Their claim is then a lifetime, radiation rate, or orbital stability statement, not a static minimum.

Boundaries and curvature. A finite box, compact space, impurity, or curved metric introduces a length and changes the scaling variation. A profile stabilized by its container is not automatically a soliton of the infinite-volume theory.

Higher derivatives or noncanonical kinetics. These alter the exponents and may stabilize a size, but positivity, hyperbolicity, additional modes, and EFT control must be checked independently.

The shared stability taxonomy shows why passing the virial test is necessary but never sufficient: existence, the Hessian spectrum, nonlinear evolution, and quantum corrections remain separate. The family-by-family scaling mechanisms, fluctuation claims, and failure boundaries are collected in the soliton boundary and stability comparison.

The second scale derivative gives information only along the dilation direction. For the two-plus-four-derivative example in D=3D=3,

d2Edλ21=2E2+12E0>0,\left.\frac{\mathrm d^2E}{\mathrm d\lambda^2}\right|_1 =2E_2+12E_0>0 ,

after using the virial identity. The solution is stable against infinitesimal uniform rescaling. It may still have a negative fluctuation with a different shape, so the full Hessian must be studied:

Hab(x,y)=δ2Eδϕa(x)δϕb(y)ϕsol.\mathcal H_{ab}(\mathbf x,\mathbf y) =\left. \frac{\delta^2E}{\delta\phi^a(\mathbf x)\delta\phi^b(\mathbf y)} \right|_{\phi_{\mathrm{sol}}}.

Zero eigenvalues generated by exact symmetries are collective coordinates, not instabilities. A negative eigenvalue is an unstable direction. Positive spectrum apart from normalizable symmetry zero modes establishes linear stability, not a general theorem of nonlinear or quantum stability.

For an analytic ansatz or numerical profile:

  1. list every contribution to the conserved energy, including boundary and constraint terms;
  2. declare how every field scales so that gauge covariance and boundary data are preserved;
  3. compute each exponent and the residual RD\mathcal R_{\mathrm D};
  4. check the field equations and RD\mathcal R_{\mathrm D} independently;
  5. identify the pair of terms that prevents shrinking and spreading;
  6. verify that the stabilizing term is inside its domain of validity;
  7. only then compute the full fluctuation spectrum.

Stop if the scale variation changes charge, leaves the admissible field space, or probes an EFT cutoff. In those cases the naive residual is not a valid test.

Treating a vanishing virial residual as a solution. It is one integrated consequence of the field equations. Many nonsolutions can satisfy it accidentally.

Applying the scalar theorem to a gauge–Higgs system. Gauge covariance fixes how the gauge potential scales. Omitting its energy or rescaling it inconsistently invalidates the identity.

Using an uncontrolled higher-derivative term as a stabilizer. A new term can balance the scaling while every still-higher operator is equally important. The solution is predictive only if its gradients remain below the cutoff and the truncation is demonstrably ordered.

  1. For E(λ)=λ2DE2+λDE0E(\lambda)=\lambda^{2-D}E_2+\lambda^{-D}E_0, derive the stationarity condition and classify D=1,2,3D=1,2,3 when both contributions are nonnegative.
Solution

Differentiation at λ=1\lambda=1 gives (2D)E2DE0=0(2-D)E_2-DE_0=0. For D=1D=1, this is E2=E0E_2=E_0 and permits a balance. For D=2D=2, it requires E0=0E_0=0, leaving scale invariance of E2E_2. For D=3D=3, it reads E23E0=0-E_2-3E_0=0, so only the vacuum with both terms zero is possible under the hypotheses.

  1. In D=3D=3, suppose E=E2+E4E=E_2+E_4 with positive two- and four-derivative terms. Find the stationary relation and show that the dilation mode is locally stable.
Solution

The scaled energy is E(λ)=λ1E2+λE4E(\lambda)=\lambda^{-1}E_2+\lambda E_4. Stationarity gives E4=E2E_4=E_2. The second derivative is 2E2>02E_2>0. This proves positivity only in the uniform-scale direction; the remaining fluctuation spectrum is not determined.

Lumps, Textures, and Skyrmions applies the test to scale-invariant and higher-derivative models. Q-Balls, Oscillons, and Sphalerons compares constrained, dynamical, and unstable-saddle evasions.

  • Coleman, Sidney. Aspects of Symmetry: Selected Erice Lectures. Cambridge University Press, 1985, ch. 6, pp. 185–264. DOI.
  • Derrick, G. H. “Comments on Nonlinear Wave Equations as Models for Elementary Particles.” Journal of Mathematical Physics 5 (1964): 1252–1254. DOI.
  • Manton, Nicholas, and Paul Sutcliffe. Topological Solitons. Cambridge Monographs on Mathematical Physics. Cambridge University Press, 2004, §§ 4.2 and 9.1. DOI.
  • Skyrme, T. H. R. “A Non-Linear Field Theory.” Proceedings of the Royal Society of London A 260 (1961): 127–138. DOI.