Constrained instantons in scalar field theories

This paper revitalizes Affleck's 1981 concept of constrained instantons into a complete method for computing vacuum decay rates in scalar field theories lacking traditional saddle points, demonstrating its efficacy through numerical solutions in a massive scalar theory with negative quartic interaction that reveal a distinct two-branch structure separating constrained instantons from action minima.

The Gist

The Big Picture: The Unstable Hill

Imagine you are standing on top of a small, grassy hill. In the world of classical physics (the rules of everyday life), if you stand perfectly still, you stay there. You are stable.

But in the quantum world (the world of atoms and subatomic particles), things are different. You might suddenly "tunnel" through the hill and roll down to a much deeper valley on the other side. This is called vacuum decay. It's like the universe spontaneously deciding to fall into a new, lower-energy state.

Physicists usually calculate how likely this is to happen by finding a special "tunneling path" called an instanton. Think of an instanton as a perfect, magical map that shows exactly how the particle rolls over the hill.

The Problem: The Map Doesn't Exist

Here is the catch: In some specific theories (like the one studied in this paper), the hill is shaped in such a way that no perfect map exists.

Imagine the hill is made of slippery ice. If you try to find a path to roll down, you realize that no matter how you try to position yourself, you can always slide a little bit further down. There is no "saddle point" (a stable spot to pause and plan your route). Because the map doesn't exist, standard math breaks down, and physicists can't calculate how fast the universe might collapse.

The Solution: The "Constrained" Path

This is where the authors, Benjamin Elder, Kinga Gawrych, and Arttu Rajantie, step in with a clever trick called Constrained Instantons.

The Analogy: The Hiker and the Rope
Imagine you are a hiker trying to cross a slippery, shifting mountain. You can't find a stable path because the ground keeps sliding under your feet.

• The Old Way: You try to find a path on the open mountain. You fail because the ground is too unstable.
• The New Way (Constrained Instanton): You tie a rope to a specific point on the mountain. You tell yourself, "I will only look for paths that pass exactly through this point on the rope."

By forcing the path to go through a specific "constraint" (the rope), you stabilize the problem. Suddenly, a path appears! You can now calculate the tunneling rate based on this forced path.

The authors took this idea, which was first proposed in 1981, and turned it into a complete, modern toolkit. They didn't just guess; they built a rigorous mathematical method to find these "forced" paths and prove they are the right ones to use.

The Experiment: Two Different Ropes

To test their method, the authors applied it to a specific type of unstable field (a "scalar field" with a negative interaction). They tried tying the rope in two different ways:

1. The $\phi^3$ Constraint: Imagine tying the rope based on the cube of the field's height.
2. The $\phi^6$ Constraint: Imagine tying the rope based on the sixth power of the field's height.

They used powerful computers to solve the equations for these two scenarios.

The Discovery: The Two-Branch Mystery

When they looked at the results, they found something fascinating: Every time they set a constraint, they found two different solutions.

Think of it like this: You ask, "Show me a path that goes through point X."

• Solution A (The Instanton): This is the "tunneling" path. It has a specific property (a "negative mode") that means it represents a real, unstable transition. This is the path that leads to the universe decaying.
• Solution B (The Minimum): This is a "safe" path. It's the lowest energy state given that you are tied to the rope. It's stable and doesn't lead to tunneling.

The authors had to figure out which one was which. They did this by counting "negative modes" (a technical way of checking if the path is unstable).

• The Result: They found that for both types of ropes ($\phi^3$ and $\phi^6$), the "upper branch" of solutions were the real tunneling paths (instantons), and the "lower branch" were just stable, safe spots.

Why This Matters

1. Solving the Unsolvables: This paper provides a way to calculate decay rates in theories where standard methods fail completely. It's like finding a way to calculate the speed of a car that doesn't have an engine, by forcing it to roll down a specific ramp.
2. A New Toolkit: They proved that this "constrained" method works even when the solutions aren't just tiny tweaks of the old, easy solutions. It works for big, complex changes too.
3. Future Applications: While they used a simple "toy model" for this paper, the method can be applied to real-world problems, like:
   • The Higgs Field: Checking if our universe's vacuum is truly stable or if it could collapse in the distant future.
   • Baryon Number Violation: Understanding how matter might be created or destroyed in the early universe.

The Bottom Line

The authors took a 40-year-old idea, refined it into a precise mathematical machine, and showed that even when nature refuses to give us a clear path, we can force a path into existence by adding a "constraint." By doing so, they can finally calculate how likely it is for the universe to undergo a dramatic, non-perturbative change.

They didn't finish the whole calculation (that's for a future paper), but they built the engine and proved it runs. Now, we can use this engine to explore the deepest, most unstable corners of quantum physics.

Technical Summary

1. Problem Statement

The paper addresses a fundamental limitation in non-perturbative quantum field theory (QFT): the computation of vacuum decay rates in theories where no standard instanton solutions exist.

• Context: Vacuum decay (tunneling from a false vacuum to a true vacuum) is typically calculated using the saddle-point approximation of the Euclidean path integral. The dominant contribution comes from "instantons"—localized saddle points of the action with a single negative mode.
• The Issue: In certain theories, such as massive $\phi^4$ scalar field theory with a negative quartic coupling (or electroweak theory with a Higgs field), the action lacks stationary points other than the false vacuum. This is due to scale invariance (in the massless limit) or the presence of a mass term that allows the action to be monotonically reduced by spatial rescaling. Consequently, standard instanton methods fail.
• Previous Work: Affleck (1981) proposed "constrained instantons" to solve this by fixing the "size" of the instanton via a constraint in the path integral. However, previous applications were limited to perturbative regimes where the solution is a small deformation of the massless instanton.
• Goal: To develop a fully non-perturbative formulation of the constrained instanton method to compute vacuum decay rates in theories where no instanton exists, specifically testing this on massive scalar field theory with negative self-coupling.

2. Methodology

The authors generalize Affleck's approach by defining a constraint functional $\xi[\phi]$ rather than relying on the concept of "size" (which is ill-defined in massive theories).

A. Theoretical Framework

1. Path Integral with Constraint: The path integral is restricted to field configurations satisfying $\xi[\phi] = \bar{\xi}$ using a delta functional.
2. Lagrange Multiplier: To find the saddle points (constrained instantons) numerically, they introduce a modified action:
   $$\tilde{S}_\kappa[\phi] = S[\phi] + \kappa \xi[\phi]$$
   where $\kappa$ is a Lagrange multiplier. The stationary points of $\tilde{S}_\kappa$ correspond to constrained instantons for a specific $\bar{\xi} = \xi[\phi]$.
3. Vacuum Decay Rate Formula: The decay rate $\Gamma$ is expressed as an integral over the constraint parameter $\bar{\xi}$:
   $$\Gamma \approx \int d\bar{\xi} \left( \frac{S[\phi_{\bar{\xi}}]}{2\pi} \right)^2 \left| \frac{\text{Det}_{\bar{\xi}} M_{\bar{\xi}}}{\text{Det} M_0} \right|^{-1/2} e^{-(S[\phi_{\bar{\xi}}] - S[\phi_0])}$$
   where $M$ is the fluctuation operator.

B. Numerical Implementation

• Model: A real scalar field in 4D Euclidean space with potential $V(\phi) = \frac{1}{2}m^2\phi^2 - \frac{\lambda}{4!}\phi^4$.
• Constraints Tested: Two distinct constraint operators were used to test robustness:
  1. Cubic Constraint: $O = \phi^3$ (scaling dimension $d=3$).
  2. Sextic Constraint: $O = \phi^6$ (scaling dimension $d=6$).
• Equation of Motion: The authors solve the resulting non-linear Ordinary Differential Equation (ODE) for the spherically symmetric profile $\phi(r)$ using the shooting method.
  $$\phi'' + \frac{3}{r}\phi' - m^2\phi + \frac{\lambda}{6}\phi^3 + \kappa \frac{\delta \xi}{\delta \phi} = 0$$
• Negative Mode Counting: A critical step is identifying which solutions are true instantons (contributing to decay) and which are just minima.
  • They utilize the relation between the constrained fluctuation determinant and the unconstrained modified Hessian $\tilde{M}_\kappa$.
  • Using the Cauchy Interlacing Theorem, they determine that if the unconstrained operator $\tilde{M}_\kappa$ has one negative mode, the constrained operator has one negative mode if the projection factor $\nu(\bar{\xi}) > 0$, and zero if $\nu(\bar{\xi}) < 0$.

3. Key Results

A. Two-Branch Structure

For both constraint types ($\phi^3$ and $\phi^6$), the numerical solutions exhibit a two-branch structure when plotting the Action $S$ versus the constraint value $\bar{\xi}$:

1. Upper Branch: Corresponds to small absolute values of the Lagrange multiplier $|\kappa|$.
   • These solutions possess exactly one negative mode.
   • They are identified as the constrained instantons.
   • As $\bar{\xi} \to 0$ (for $\phi^3$) or $\bar{\xi} \to \infty$ (for $\phi^6$), the action approaches the massless instanton action ($16\pi^2/\lambda$).
2. Lower Branch: Corresponds to large $|\kappa|$.
   • These solutions have zero negative modes.
   • They are identified as minima of the action subject to the constraint, not instantons. They do not contribute to the tunneling rate.

B. Behavior of the Action

• $\phi^3$ Constraint: The constraint $\bar{\xi}$ has a maximum value $\bar{\xi}_{max}$. The two branches meet at a cusp at this maximum. The action on the instanton branch is always higher than the massless instanton action but approaches it as $\kappa \to 0$.
• $\phi^6$ Constraint: The constraint $\bar{\xi}$ has a minimum value $\bar{\xi}_{min}$. The branches meet at a cusp at this minimum. The action decreases rapidly as $\kappa$ approaches its upper limit.
• Consistency: The results for both constraints are consistent with the theoretical expectation that the relevant solutions correspond to "small" instanton sizes (where the mass term is negligible compared to the quartic term).

C. Validation

The authors performed rigorous consistency checks:

• Lagrange Multiplier Relations: Verified that $\bar{\xi} = d\tilde{S}/d\kappa$ and $\kappa = -dS/d\bar{\xi}$ numerically.
• Scaling Identities: Checked integral identities derived from scaling transformations of the modified action.
• Numerical Stability: Confirmed that results are independent of the simulation box size ($R_{min}, R_{max}$) and working precision, with errors following expected exponential/power-law decay patterns.

4. Significance and Contributions

1. Non-Perturbative Formulation: The paper moves beyond the perturbative regime of Affleck's original work. It provides a complete, non-perturbative method to find constrained instantons even when they are not small perturbations of massless solutions.
2. Resolution of the "No Instanton" Problem: It demonstrates a viable mathematical and numerical pathway to compute vacuum decay rates in theories (like massive $\phi^4$ with negative coupling) where standard instanton calculus fails due to the absence of stationary points.
3. Identification of Physical Solutions: The method of using the projection factor $\nu(\bar{\xi})$ to distinguish between constrained instantons (negative mode) and constrained minima (no negative mode) is a crucial technical contribution. It resolves the ambiguity of the two-branch solution structure.
4. Future Applications: The authors argue that this method is directly applicable to more complex and physically relevant scenarios, such as:
   • Electroweak Vacuum Metastability: Where the Higgs potential behaves similarly to the model studied.
   • Baryon Number Violation: Calculating rates in the Standard Model where instantons are forbidden by the Higgs mechanism.
   • Non-renormalizable Operators: Studying the impact of higher-dimensional operators on tunneling rates.

Conclusion

The authors successfully developed and validated a robust numerical framework for calculating vacuum decay rates using constrained instantons. By applying this to a massive scalar field theory with negative quartic coupling, they confirmed the existence of a two-branch solution structure and successfully isolated the physical instanton branch. This work lays the groundwork for computing non-perturbative rates in theories previously considered intractable due to the lack of standard instanton solutions.