Non-perturbative False Vacuum Decay Using Lattice Monte Carlo in Imaginary Time

This paper introduces a non-perturbative lattice Monte Carlo method in imaginary time, featuring a novel sampling technique to overcome false vacuum suppression and a new Fermi's Golden Rule-like formula, to accurately calculate quantum tunneling rates, as validated by reproducing one-dimensional Schrödinger equation results.

The Gist

The Big Picture: The "Sleeping Giant" Problem

Imagine the universe is a giant landscape. Usually, things settle into the lowest valley possible—that's the True Vacuum, the most stable, happy state for everything.

But sometimes, the universe gets stuck in a False Vacuum. Think of this as a ball sitting in a small, shallow dip on the side of a mountain. It looks stable; it's resting comfortably. But there is a huge mountain peak (a barrier) separating it from the deep valley below.

In the quantum world, even if the ball doesn't have enough energy to roll over the mountain, it can still tunnel through it. It's like the ball suddenly teleporting to the other side. Once it does, it rolls down into the deep valley, releasing a massive amount of energy. In physics, this is called False Vacuum Decay. If our universe is in a false vacuum, this could theoretically happen at any moment, rewriting the laws of physics instantly.

The Problem: Physicists want to know: How long will the ball stay in that shallow dip before it tunnels out?

• Old methods (Semi-classical) are like guessing the answer based on a rough sketch of the mountain. They work well if the mountain is simple, but they fail if the terrain is complex and "strongly coupled" (messy).
• The Goal: This paper introduces a new, super-precise way to calculate this tunneling time using a computer simulation, without relying on rough guesses.

---

The Challenge: The "Ghost" in the Machine

The authors tried to simulate this on a computer using a method called Lattice Monte Carlo. Imagine the computer is trying to walk through the mountain landscape, step by step, to see how often the ball falls into the deep valley.

But there are three huge problems:

1. The "Rare Event" Problem: The ball staying in the shallow dip is very common. The ball tunneling out is incredibly rare. In a standard simulation, the computer spends 99.99% of its time walking around the shallow dip and almost never sees the tunneling event. It's like trying to find a specific grain of sand on a beach by walking the beach once.
2. The "Ergodicity" Problem (The Maze): The computer gets stuck. It wanders around the shallow dip and can't figure out how to cross the mountain to the other side. It's like a mouse in a maze that only knows one room; it never finds the exit.
3. The "Time Travel" Problem: Real tunneling happens in real time. But the computer simulation works in Imaginary Time (a mathematical trick that turns time into a distance). It's like trying to calculate how fast a car drives by looking at a map of where it would be if it drove forever, rather than watching it drive. You have to translate the map back to reality, which is tricky.

---

The Solution: The "Magic Mirror" and the "Staircase"

The authors developed a clever three-part strategy to solve these problems.

1. The "Magic Mirror" (Implicit Decay Amplitude)

Instead of waiting for the ball to actually tunnel (which takes forever), they changed the rules of the game slightly.

• They created a "mirror world" (a modified Hamiltonian) where the mountain barrier is slightly different.
• They derived a new formula (similar to Fermi's Golden Rule, a famous physics equation) that says: "You don't need to watch the ball fall. You just need to measure how much the ball 'wiggles' at the edge of the cliff."
• By measuring this "wiggle" (the decay amplitude) in the mirror world, they can mathematically calculate the tunneling rate without ever waiting for the rare event to happen.

2. The "Staircase" (Intermediate Ratios)

To solve the problem of the computer getting stuck in the shallow dip, they built a staircase.

• Imagine you want to walk from the shallow dip (Action A) to the deep valley (Action B), but the jump is too big.
• Instead of jumping, they built a series of small steps (Intermediate Actions) between A and B.
• They ran many small simulations, stepping from one rung to the next. By multiplying the results of these small steps, they could bridge the gap between the two worlds.
• Analogy: Instead of trying to jump a 10-foot gap, you build a ladder with 10 rungs. You climb one rung at a time, and the computer never gets "stuck" because the steps are small and easy.

3. The "Spectral Reconstruction" (Reading the Crystal Ball)

Since the computer works in "Imaginary Time," the data it gives back is a bit blurry. It's like looking at a reflection in a foggy mirror.

• The computer gives them a curve of data over time.
• They need to reverse-engineer this curve to find the "spectrum" (the specific energy levels involved in the tunneling).
• They used a Gaussian Fit (a bell curve) to guess the shape of this spectrum.
• The Catch: This is the weakest link. It's like trying to guess the exact shape of a cloud just by looking at its shadow. Sometimes the guess is off by a factor of two, but it's still the best we can do with current tools.

---

The Results: Did it Work?

The authors tested their new method on a simple 1D quantum system (a particle in a box with a bump in the middle).

• They compared their computer results against the "Exact" answer (solved using the standard Schrödinger equation, which is the gold standard for simple problems).
• The Verdict: Their method worked! They could calculate decay rates that were incredibly small (like 1 in a trillion).
• The Accuracy: The results were very close to the exact answer, usually within a factor of 2. The main source of error was the "foggy mirror" (spectral reconstruction), not the staircase method itself.

Why Does This Matter?

This paper is a proof-of-concept. It's like building a new type of engine in a garage to see if it runs.

• For now: It proves we can calculate these "impossible" tunneling rates without using rough approximations.
• For the future: The authors hope to use this method on full Quantum Field Theories (like the Standard Model of particle physics). If the universe is indeed in a false vacuum, this method could help us calculate exactly how long we have before the "ball" rolls down the mountain.

In a nutshell: They built a digital ladder to climb out of a simulation trap, used a magic mirror to predict a rare event, and successfully calculated how long a quantum system stays stuck in a "false" state. It's a major step toward understanding the ultimate fate of the universe.

Technical Summary

1. Problem Statement

The paper addresses the challenge of calculating false vacuum decay rates (quantum tunneling rates) in strongly coupled quantum field theories and quantum mechanical systems.

• Context: False vacuum decay occurs when a metastable state (false vacuum) tunnels through a potential barrier to a lower energy state (true vacuum). This is relevant for first-order phase transitions, baryogenesis, and the stability of the Standard Model electroweak vacuum.
• Limitations of Existing Methods:
  • Semi-classical methods (Callan-Coleman): Rely on tree-level effective actions and perturbative expansions, which fail for strongly coupled theories or when quantum effects create false vacua not apparent in the classical Lagrangian.
  • Lattice Monte Carlo Challenges:
    1. Exponential Suppression: Configurations corresponding to the false vacuum have much larger Euclidean actions than the ground state, leading to exponential suppression in the path integral.
    2. Ergodicity: Markov chains struggle to traverse the configuration space separated by high potential barriers.
    3. Time Domain Mismatch: Lattice simulations are performed in imaginary (Euclidean) time, whereas decay rates are inherently real-time phenomena. Direct analytic continuation is often ill-posed or impossible due to the non-exponential behavior of decay at very short and very long times.

2. Methodology

The authors propose a new non-perturbative framework combining Lattice Monte Carlo (MCMC) simulations in imaginary time with a novel spectral reconstruction technique.

A. Theoretical Framework: Implicit Decay Amplitude

Instead of simulating real-time evolution, the authors derive a formula for the decay rate $\Gamma$ analogous to Fermi's Golden Rule, but formulated entirely in terms of Euclidean-time operators.

• Hamiltonian Splitting: They define a "False Vacuum Hamiltonian" ($H_{FV}$) and a "True Vacuum Hamiltonian" ($H_{TV}$).
  • $H_{FV} = H + V_{bar, FV}$: Adds a potential barrier to prevent the system from entering the true vacuum region.
  • $H_{TV} = H + V_{bar, TV}$: Adds a barrier to prevent entry into the false vacuum region.
• The Formula: The decay rate is expressed as:
  $$\Gamma = 2\pi \langle FV | (H_{FV} - H) \delta(H_{TV} - E_{FV}) (H_{FV} - H) | FV \rangle$$
  Here, $|FV\rangle$ is the ground state of $H_{FV}$, and $E_{FV}$ is its energy. The term $(H_{FV} - H)$ acts as the transition operator (implicitly replacing the perturbation $H_{tr}$).
• Implicit Amplitude: The delta function $\delta(H_{TV} - E_{FV})$ enforces energy conservation. The matrix element is treated as an "implicit decay amplitude" that can be computed via Euclidean time evolution.

B. Spectral Reconstruction (The Inverse Problem)

To extract the decay rate from the Euclidean correlator, the authors define a quantity $Q(t)$:
$$Q(t) = \langle D | e^{-(H_{TV} - E_{FV})t} | D \rangle$$
where $|D\rangle = (H_{FV} - H)|FV\rangle$.

• The Challenge: Recovering the spectral density $\rho(E)$ (specifically at $E = E_{FV}$) from $Q(t)$ is an ill-posed inverse problem, especially with noisy Monte Carlo data.
• The Solution: The authors assume the spectral density $\rho(E)$ follows a Gaussian distribution centered at $E_{FV}$. They fit the calculated $Q(t)$ data to this Gaussian ansatz to extract the peak height, which corresponds to the decay rate. They argue this is valid for intermediate times where high-energy states are suppressed but the ground state hasn't completely dominated.

C. Sampling Strategy: Intermediate Ratios Method

To overcome the exponential suppression and ergodicity issues in the path integral:

• Problem: Directly sampling the ratio of path integrals for the numerator (involving $H_{FV}$) and denominator (involving $H$) is impossible because the ensembles are disjoint; configurations important for the numerator are exponentially suppressed in the denominator.
• Solution: The authors introduce a chain of intermediate actions ($S_L$ and $S_M$) that interpolate smoothly between the numerator and denominator actions.
  • They define families of actions $S_L(L)$ and $S_M(M)$ that ensure the action in the denominator is always less than or equal to the action in the numerator for any configuration ($S_D \leq S_N$).
  • They calculate the total ratio as a product of ratios of adjacent intermediate ensembles:
    $$\frac{Z_N}{Z_D} = \prod \frac{Z_{L_{j+1}}}{Z_{L_j}} \times \left( \prod \frac{Z_{M_{k+1}}}{Z_{M_k}} \right)^{-1}$$
  • This allows the use of standard MCMC to sample configurations that are relevant for both the initial and final states, effectively bridging the gap between the false vacuum and true vacuum regions.

3. Key Contributions

1. New Non-Perturbative Formula: Derivation of a decay rate formula based on an implicit decay amplitude that avoids semi-classical approximations and works for strongly coupled systems.
2. Imaginary Time to Real Rate: A robust method to extract real-time decay rates from Euclidean lattice data without requiring direct analytic continuation of the full time-dependent signal.
3. Intermediate Ratios Sampling: A novel sampling technique that solves the ergodicity and signal suppression problems by constructing a continuous path of intermediate actions, allowing the calculation of extremely small decay rates (down to $10^{-15}$ in their tests) without restricting the lattice volume.
4. Spectral Reconstruction Strategy: A practical approach using Gaussian fitting for spectral reconstruction, validated against exact solutions.

4. Results

The method was tested on a one-dimensional single-particle quantum system with a double-well potential (Equation 5 in the paper).

• Comparison: The authors compared three methods:
  1. Exact: Numerical solution of the Schrödinger equation.
  2. Implicit Amplitude (FGR): Using the derived formula with exact spectral data (no Monte Carlo noise).
  3. Lattice Monte Carlo: Using the full proposed method with MCMC and spectral reconstruction.
• Performance:
  • The Monte Carlo results successfully reproduced the exact decay rates across a wide range of parameters ($\alpha$ and $\beta$).
  • The method calculated decay rates as small as $\sim 10^{-15}$.
  • Accuracy: The Monte Carlo results differed from the exact results by a factor of less than 2.
• Error Analysis:
  • Statistical Error: Controlled via standard jackknife analysis.
  • Systematic Error: The dominant source of error was spectral reconstruction (the Gaussian ansatz). The authors note that while the error is currently a factor of ~2, it is systematically improvable by using more sophisticated reconstruction methods (e.g., Maximum Entropy, Bayesian reconstruction, or Hansen-Lupo-Tantalo methods).

5. Significance and Future Work

• Significance: This work provides the first fully non-perturbative lattice method capable of calculating false vacuum decay rates for strongly coupled systems without relying on semi-classical approximations. It overcomes the fundamental barrier of simulating rare tunneling events in Euclidean time.
• Future Directions:
  • Field Theory: Generalizing the method to full Quantum Field Theories (e.g., QCD or Electroweak theory). The authors argue that spectral reconstruction might actually be easier in field theory due to the higher density of states, which naturally leads to a sharper energy peak after Euclidean evolution.
  • Improved Reconstruction: Implementing advanced spectral reconstruction algorithms to reduce the systematic error from the current factor-of-two level.
  • Finite Temperature: Extending the formalism to finite temperature systems.

In summary, Jin and Swaim have established a viable pathway for computing non-perturbative tunneling rates on the lattice, solving the critical issues of ergodicity and real-time extraction, and demonstrating its efficacy on a benchmark quantum mechanical system.