Tunneling from an oscillating initial state in quantum… — Plain-Language Explanation

Imagine you are in a deep valley (a "metastable well") surrounded by a high mountain pass (a "barrier"). In the world of classical physics, if you don't have enough energy to climb over the mountain, you are stuck there forever. But in the quantum world, particles have a weird superpower: they can "tunnel" through the mountain, appearing on the other side even without climbing over it.

This paper is about figuring out exactly how fast a particle escapes this valley, but with a twist: the particle isn't just sitting still at the bottom of the valley. It's oscillating, bouncing back and forth like a ball in a bowl.

Here is the breakdown of their discovery using simple analogies:

1. The Problem: A Bouncing Ball vs. a Still Ball

Usually, scientists calculate tunneling for a particle that is perfectly still at the bottom of the valley (the "ground state"). It's like a ball sitting quietly; it leaks out very slowly and steadily.

But in many real-world situations (like in superconducting circuits or the early universe), the particle is moving. It's oscillating back and forth. The authors asked: Does the fact that the particle is moving change how it escapes?

2. The Solution: Breaking the Motion into "Resonant States"

To solve this, the authors used a mathematical trick. Imagine the bouncing particle is actually a choir of many different singers, each singing a specific note (a "resonant state").

Some notes are low and slow; others are high and fast.
Each note has its own specific "leakiness" (how easily it tunnels through the mountain).
Because the particle is a mix of all these notes, they interfere with each other.

The authors derived a master formula (Equation 18) that adds up all these individual notes. It tells you not just the average rate of escape, but the exact probability of the particle escaping at any specific moment in time.

3. The Big Surprise: The "Burst" Effect

The most exciting finding is what happens when the particle is oscillating coherently (moving in a smooth, rhythmic pattern).

The Old View: You might expect the particle to leak out at a steady, slow drip, like water leaking from a bucket.
The New View: The paper shows the particle doesn't leak steadily. Instead, it leaks in sudden, sharp bursts.

The Analogy: Think of a person trying to sneak out of a guarded house through a narrow, dark tunnel.

If they just stand in the hallway, they might slip out slowly.
But if they are running back and forth, they only have a chance to slip through the tunnel when they are closest to the door.
Every time they bounce off the wall and rush toward the tunnel entrance, there is a tiny window of opportunity where the "quantum magic" works best.

The authors found that the particle escapes almost entirely during these brief moments when it is closest to the barrier. For the rest of the time, it is effectively trapped. This creates a "spiky" pattern of escape rather than a smooth curve.

4. The "Saddle-Point" Shortcut

Calculating this for every single moment is incredibly hard. The authors used a method called the "saddle-point approximation."

The Metaphor: Imagine a hiker trying to cross a mountain range. Instead of checking every single path, they realize the hiker will almost certainly take the one specific pass that is the lowest point.
In their math, they found that the "escape" happens almost exclusively at one specific point in the particle's oscillation cycle (the classical turning point). They calculated the exact width and height of these escape "bursts" using this shortcut.

5. What They Tested

They didn't just do math on paper; they ran computer simulations to prove it works.

They simulated a particle in a valley with a barrier.
They compared their new formula against the raw computer simulation.
The Result: The formula matched the simulation perfectly. It correctly predicted the "spiky" bursts of escape and the exact timing of when the particle would leak out.

6. Why It Matters (According to the Paper)

The paper notes that this is crucial for understanding:

Superconducting Circuits: Specifically, Josephson junctions where current flows. The decay rate depends on whether the system is in a quiet state or an excited, oscillating state.
Cosmology: The early universe might have had fields (like axion dark matter) that were oscillating. If these fields were trying to "tunnel" to a lower energy state (creating bubbles of a new universe), this paper suggests they would do so in rhythmic bursts rather than a steady stream.

Summary

The paper provides a new, precise recipe for calculating how a moving, oscillating quantum particle escapes a trap. It reveals that instead of leaking out slowly and evenly, the particle waits until it is closest to the exit, then "pops" out in a rapid, rhythmic burst. This happens because the different "notes" of the particle's motion interfere with each other to create these precise moments of opportunity.

Technical Summary: Tunneling from an Oscillating Initial State in Quantum Mechanics

Problem Statement
The paper addresses the calculation of quantum tunneling decay rates for general initial states localized in a metastable potential well, specifically focusing on cases where the initial state is not the perturbative vacuum (ground state). While tunneling from the ground state is well-understood via WKB or semiclassical path integrals, many physically relevant scenarios involve excited or time-dependent initial states. Examples include current-biased Josephson junctions starting in excited states and cosmological models involving oscillating scalar fields (e.g., axion dark matter) or fields evolving along slow-roll trajectories. Existing frameworks for non-vacuum tunneling are often incomplete; in quantum field theory, different approaches yield conflicting results even at the level of the exponential exponent, and previous quantum mechanical treatments often rely on fitting coefficients to numerical data rather than providing closed-form analytic expressions.

Methodology
The authors develop a framework based on the expansion of the initial wavefunction in a basis of resonant (Gamow-Siegert) states. Unlike standard Hermitian eigenstates, these states satisfy outgoing boundary conditions, possess complex energies $\epsilon_n = E_n - i\hbar\Gamma_n/2$ , and describe the decay of probability from the metastable well.

Resonant State Expansion: The initial state $\psi(x)$ is decomposed into a sum of resonant states $\psi_n(x)$ with coefficients $c_n$ . Because the resonant Hamiltonian is non-Hermitian, the states are not strictly orthogonal or normalizable in the standard $L^2$ sense. The authors define a physically meaningful normalization within the well and construct an overlap matrix $S_{mn}$ to determine the coefficients $c_n$ via matrix inversion.
Probability Current Calculation: The time-dependent probability current $j(x, t)$ $j (x, t)$ at the barrier exit is derived by evolving the resonant components. The resulting expression (Eq. 18) contains two distinct terms:
- A time-averaged decay term proportional to $|c_n|^2 \Gamma_n e^{-\Gamma_n t}$ .
- An interference term arising from the cross-terms between different resonant states, which introduces oscillations in the current.
Semiclassical Approximation (WKB): The authors compute all ingredients ( $E_n$ $E_{n}$ , $\Gamma_n$ $Γ_{n}$ , $c_n$ $c_{n}$ ) analytically using the WKB approximation to first subleading order.
- Energies and Widths: Complex energies are derived by imposing purely outgoing boundary conditions on WKB wavefunctions, yielding the standard Breit-Wigner form for the width $\Gamma_n \propto e^{-2S_n/\hbar}$ .
- Coherent States: Specializing to a coherent initial state (a superposition of harmonic oscillator eigenstates), the authors apply a saddle-point approximation to the sum over resonant states. This assumes the dominant contributions come from high quantum numbers ( $n \gg 1$ ).

Key Results

Closed-Form Expression for Current: The primary result is Eq. (18), a closed-form expression for the probability current through the barrier. This formula captures both the exponential decay of individual resonances and the interference effects between them.
Time-Structured Tunneling: For a coherent initial state, the tunneling current is not a smooth exponential decay. Instead, it is concentrated in narrow bursts (spikes) occurring near the classical turning point closest to the barrier.
- The time window of effective tunneling is $\Delta t \sim (\omega \sqrt{n_0})^{-1}$ , which is much shorter than the oscillation period $2\pi/\omega$ .
- The total probability leaked per oscillation cycle is given by Eq. (42), which depends exponentially on the semiclassical action evaluated at the dominant resonance energy.
Role of Interference: The authors demonstrate that the "spiky" structure is a direct consequence of phase coherence among the resonant states. If the phases of the expansion coefficients are randomized, the interference terms average out, and the current reverts to a smooth decay tracking the average rate $\bar{\Gamma} = \sum |c_n|^2 \Gamma_n$ .
Numerical Verification: The analytic results are validated against numerical solutions of the Schrödinger equation using a Complex Absorbing Potential (CAP) to simulate outgoing boundary conditions.
- For low-energy coherent states ( $|\alpha|=1.1$ ), the fully analytic WKB prediction matches the numerical current and cumulative probability with high accuracy.
- For higher-energy states ( $|\alpha|=2$ ), where energies approach the barrier top, the WKB approximation for the decay width $\Gamma_n$ shows larger deviations, though the resonant-state decomposition formula itself remains accurate when exact numerical complex energies are used.

Significance and Claims
The paper claims to provide a fully analytic, closed-form expression for the time-dependent tunneling current of non-vacuum states, avoiding the need for numerical fitting of coefficients found in previous works (e.g., [10, 11]).

Contrast with Existing Methods: The authors note that their result contrasts favorably with path integral methods in quantum field theory, which currently disagree on the exponent, and with previous quantum mechanical approaches that required numerical input.
Physical Insight: The work offers a precise quantum-mechanical realization of the intuitive picture that a particle in a well "attempts" to tunnel once per bounce off the barrier, with the tunneling probability exponentially enhanced at the classical turning point.
Limitations: The authors explicitly state that the resonant-state expansion is not a complete spectral decomposition; it neglects the non-resonant continuum. Consequently, the formula is accurate only after a brief initial transient period (of order one oscillation) and before very late times where power-law tails dominate. Furthermore, the WKB approximation for decay widths becomes less reliable as the energy approaches the top of the barrier.

The paper concludes by identifying the extension of this framework to Quantum Field Theory (specifically for oscillating scalar fields and bubble nucleation) as the most important future direction, noting that the path integral formulation in [12] may serve as an entry point.

Tunneling from an oscillating initial state in quantum mechanics