Tunneling Potential Method Overview

• Tunneling potential method is a set of reformulations that calculate decay probabilities via auxiliary potential functions rather than traditional boundary-value problems.
• It applies to diverse areas including one-dimensional quantum barrier scattering, vacuum decay in scalar fields, and Q-ball theory, offering both numerical and analytical insights.
• The approach transforms complex variational problems into lower-dimensional functionals, leading to more stable computations and exact analytic benchmarks.

The tunneling potential method denotes a family of tunneling formalisms in which a decay probability, decay exponent, or transmitted amplitude is computed through an auxiliary potential-like object rather than exclusively through a direct coordinate-space boundary-value problem. In one-dimensional quantum mechanics, the phrase can refer to direct integration of the time-independent Schrödinger equation from the transmission region backward through an arbitrary barrier, followed by matching to incident and reflected waves (Rundquist, 2011). In relativistic vacuum decay, it refers more specifically to the introduction of a field-space function $V_t$ that connects the metastable and stable phases and turns the bounce calculation into a variational problem for a one-dimensional action functional in field space (Espinosa, 2018). Closely related constructions also appear in Q-ball theory, gauge-invariant false-vacuum decay, and several applied tunneling problems (Espinosa et al., 2023).

1. Terminological scope and defining constructions

A common misconception is that the tunneling potential method is a single universal algorithm. The literature instead uses the expression for several technically distinct constructions. In one usage, the relevant object is the transmitted-side boundary condition imposed on the wavefunction in a one-dimensional barrier problem. In another, the central object is an auxiliary function $V_t(\phi)$ or $\mathcal E(\varphi)$ defined in field space and optimized by an action principle. A plausible implication is that the phrase is best understood as a class of reformulations rather than a single formalism.

Setting	Auxiliary object	Representative source
One-dimensional barrier scattering	Backward-integrated $\Psi(x)$ with transmitted-side data	(Rundquist, 2011)
Vacuum decay of a scalar field	$V_t(\phi)$	(Espinosa, 2018)
Q-balls	$V_t(\varphi)$ or $\mathcal E(\varphi)$	(Espinosa et al., 2023)

For four-dimensional vacuum decay, the defining relation is

$$V_t(\phi)\equiv V(\phi)-\tfrac12\,\bigl(d\phi_b/dr\bigr)^2,$$

with $\phi_b(r)$ the $O(4)$-symmetric bounce. The associated endpoint data are

$V_t(\phi)$0

where $V_t(\phi)$1 is the false vacuum and $V_t(\phi)$2 the bounce exit point (Espinosa, 2018). For Q-balls, the analogous construction starts from a stationary complex scalar $V_t(\phi)$3 and defines

$V_t(\phi)$4

so that the radial equation becomes structurally identical to a three-dimensional bounce equation (Espinosa et al., 2023).

2. One-dimensional quantum-mechanical barrier formulations

In the direct-integration or shooting formulation of one-dimensional tunneling, the system is partitioned into Region I $V_t(\phi)$5, Region II $V_t(\phi)$6, and Region III $V_t(\phi)$7. Region I contains incident plus reflected waves, Region II contains an arbitrary barrier $V_t(\phi)$8, and Region III has constant $V_t(\phi)$9 with only a transmitted right-travelling wave. The method begins from the time-independent Schrödinger equation,

$\mathcal E(\varphi)$0

or equivalently

$\mathcal E(\varphi)$1

Rather than fixing the incident amplitude at the left, one sets $\mathcal E(\varphi)$2 in Region III and imposes at $\mathcal E(\varphi)$3

$\mathcal E(\varphi)$4

then integrates backward to $\mathcal E(\varphi)$5. Matching there to

$\mathcal E(\varphi)$6

yields

$\mathcal E(\varphi)$7

and hence

$\mathcal E(\varphi)$8

This construction automatically enforces a purely outgoing condition in Region III, typically needs fewer steps than a transfer-matrix discretization of comparable accuracy, yields the full complex $\mathcal E(\varphi)$9, and extends easily to arbitrary, even discontinuous, $\Psi(x)$0 (Rundquist, 2011).

The numerical implementation described in the same source rewrites the complex second-order ODE as four coupled first-order equations for $\Psi(x)$1, $\Psi(x)$2, $\Psi(x)$3, and $\Psi(x)$4, and advances them with a standard fourth-order Runge–Kutta scheme. The step size $\Psi(x)$5 must resolve the shortest local wavelength in oscillatory regions or the exponential decay length in forbidden regions; a rule of thumb is 20–50 points per local half-wavelength. Because the equation has no first-derivative term it is not stiff in the classical sense, although strong forbidden regions can generate overflow or underflow unless $\Psi(x)$6 is sufficiently small or periodic renormalization is used (Rundquist, 2011).

For continuous barriers, a distinct but nearby methodology uses a uniform asymptotic expansion based on Airy functions. There the forbidden-region solution is written in terms of $\Psi(x)$7 and $\Psi(x)$8, remains uniformly valid up to the turning points, and gives

$\Psi(x)$9

This approach captures the $V_t(\phi)$0 prefactor absent in naive WKB and the first beyond-WKB correction, while retaining the assumptions of a smooth barrier, linear turning points, and $V_t(\phi)$1 (Khorasani, 2011).

3. Field-space variational reformulation of vacuum decay

In the vacuum-decay formalism introduced for metastable scalar potentials, the bounce action is rewritten as a field-space functional of $V_t(\phi)$2. In flat space,

$V_t(\phi)$3

and extremizing $V_t(\phi)$4 reproduces the ordinary Euclidean bounce action. The corresponding Euler–Lagrange equation is

$V_t(\phi)$5

The Euclidean radius can then be reconstructed from

$V_t(\phi)$6

so the usual profile $V_t(\phi)$7 is not discarded but encoded indirectly in $V_t(\phi)$8 (Espinosa, 2018).

The same formalism admits an alternative derivation via canonical transformations. In that derivation one trades the canonical pair $V_t(\phi)$9 for $V_t(\varphi)$0, finds a generating function $V_t(\varphi)$1, and then Legendre-transforms to recover exactly the field-space action above. The same canonical logic extends to multi-field systems and to gravity, where the Hamiltonian constraint plays the role of the Euclidean Friedmann equation (Espinosa et al., 2022).

With gravity included, the functional becomes

$V_t(\varphi)$2

with

$V_t(\varphi)$3

The gravitational Euler–Lagrange equation is

$V_t(\varphi)$4

In the limit $V_t(\varphi)$5, the gravitational functional reduces to the flat-space one (Espinosa, 2024).

A central structural property of the method is that the physical $V_t(\varphi)$6 is a minimum of $V_t(\varphi)$7, not a saddle point. The second variation is positive definite, or equivalently the Hessian operator governing $V_t(\varphi)$8 fluctuations has a strictly positive spectrum. This is the feature on which the frequent description of the method as a downhill minimization rests (Espinosa, 2024).

4. Boundary conditions, limiting cases, and gauge-invariant implementation

The boundary data depend on the asymptotics of the false vacuum. For Minkowski or AdS false vacua,

$V_t(\varphi)$9

while at the endpoint on the true side of the barrier,

$\mathcal E(\varphi)$0

For a de Sitter false vacuum, there is always an initial segment $\mathcal E(\varphi)$1 on which $\mathcal E(\varphi)$2, and the actual bounce range is $\mathcal E(\varphi)$3 with the same endpoint condition $\mathcal E(\varphi)$4 at both ends of that range. When the range collapses to the barrier top $\mathcal E(\varphi)$5, the minimizing configuration is simply $\mathcal E(\varphi)$6 from $\mathcal E(\varphi)$7 to $\mathcal E(\varphi)$8, and the action reduces to the Hawking–Moss result

$\mathcal E(\varphi)$9

The same framework also accommodates pseudo-bounces and bubbles of nothing by relaxing the endpoint condition at $$V_t(\phi)\equiv V(\phi)-\tfrac12\,\bigl(d\phi_b/dr\bigr)^2,$$0; in these cases the Euclidean approach often requires boundary terms, while the $$V_t(\phi)\equiv V(\phi)-\tfrac12\,\bigl(d\phi_b/dr\bigr)^2,$$1-action gives the correct $$V_t(\phi)\equiv V(\phi)-\tfrac12\,\bigl(d\phi_b/dr\bigr)^2,$$2 without adding any extra boundary term (Espinosa, 2024).

Near-degenerate vacua recover the thin-wall limit. In the original field-space reformulation, the compact action density reduces to Coleman’s thin-wall expression,

$$V_t(\phi)\equiv V(\phi)-\tfrac12\,\bigl(d\phi_b/dr\bigr)^2,$$3

so the field-space formalism functions as a generalization of the thin-wall action to arbitrary potentials (Espinosa, 2018). A concrete analytic realization appears in a piecewise linear/quadratic landscape potential, where the bubble radius $$V_t(\phi)\equiv V(\phi)-\tfrac12\,\bigl(d\phi_b/dr\bigr)^2,$$4 obeys

$$V_t(\phi)\equiv V(\phi)-\tfrac12\,\bigl(d\phi_b/dr\bigr)^2,$$5

and the tunneling exponent is

$$V_t(\phi)\equiv V(\phi)-\tfrac12\,\bigl(d\phi_b/dr\bigr)^2,$$6

In the large-radius regime $$V_t(\phi)\equiv V(\phi)-\tfrac12\,\bigl(d\phi_b/dr\bigr)^2,$$7, the result reproduces the Coleman–de Luccia thin-wall expression; in the small-radius regime $$V_t(\phi)\equiv V(\phi)-\tfrac12\,\bigl(d\phi_b/dr\bigr)^2,$$8, it reproduces the Duncan–Jensen linear-potential limit. Across all $$V_t(\phi)\equiv V(\phi)-\tfrac12\,\bigl(d\phi_b/dr\bigr)^2,$$9, one of the two closed-form approximations is within 10% of the exact $\phi_b(r)$0 and within 50% in $\phi_b(r)$1 (Dutta et al., 2011).

Gauge invariance requires additional care. In the Abelian Higgs model, a naive implementation of the tunneling-potential approach inherits the gauge dependence of the bounce formalism. A gauge-invariant procedure uses the counting $\phi_b(r)$2, $\phi_b(r)$3, expands $\phi_b(r)$4 and $\phi_b(r)$5, computes the leading action from $\phi_b(r)$6, and adds the subleading corrections from $\phi_b(r)$7 and $\phi_b(r)$8 evaluated on the leading $\phi_b(r)$9. With a quartic ansatz for $O(4)$0, the method is typically 10–100 times faster than shooting the bounce ODE, and the residual gauge dependence is attributed solely to the polynomial approximation; with the refinements described in the source, it is typically at the $O(4)$1–$O(4)$2 level for $O(4)$3 (Arunasalam et al., 2021).

5. Exact solutions, multi-field embeddings, and Q-balls

The tunneling-potential formalism is unusually effective at generating exact solutions. One result is that any single-field pair $O(4)$4 can be embedded in an $O(4)$5-field potential by choosing a path $O(4)$6 in field space and extending the potential transversely so as to satisfy the multi-field equations. Along the path, the ODE for $O(4)$7 is identical to the single-field equation, while the transverse force is fixed by

$O(4)$8

This construction yields explicit analytic two-field and three-field examples, including a quarter-circle trajectory, a catenary, and a helix, and the bounce action matches the single-field value. In one thin-wall example, the exact action is $O(4)$9 and FindBounce returns $V_t(\phi)$00 with default settings, making these solutions useful as benchmarks for multi-field numerical codes (Espinosa et al., 2023).

The same paper gives a gravitating extension with

$V_t(\phi)$01

and an analytic example based on

$V_t(\phi)$02

This permits closed-form expressions for $V_t(\phi)$03, $V_t(\phi)$04, and hence $V_t(\phi)$05, with the Minkowski limit reducing to an elementary combination of $V_t(\phi)$06, $V_t(\phi)$07, and $V_t(\phi)$08 (Espinosa et al., 2023).

For Q-balls, the structural identity with three-dimensional bounce solutions allows a direct translation of the field-space strategy. Besides the definition

$V_t(\phi)$09

the problem can be reformulated in terms of a monotonic “mechanical energy” function

$V_t(\phi)$10

which satisfies

$V_t(\phi)$11

The corresponding boundary conditions are

$V_t(\phi)$12

Low-order polynomial ansätze $V_t(\phi)$13 then approximate small-$V_t(\phi)$14 and intermediate-$V_t(\phi)$15 Q-balls with percent-level accuracy, while exactly solvable examples follow by choosing $V_t(\phi)$16 directly and reconstructing $V_t(\phi)$17 from it (Espinosa et al., 2023).

6. Applied and adjacent usages

Outside false-vacuum decay, the term also appears in more application-specific settings. In nitride-based resonance-tunneling structures, an analytical method combines the effective-mass Schrödinger equation with the dielectric-continuum Poisson equation and iterates a self-consistent potential profile

$V_t(\phi)$18

At the zeroth iteration each layer has a linear slope $V_t(\phi)$19 and the Schrödinger equation is solved with Airy functions; after Poisson integration and inclusion of the Hedin–Lundquist exchange-correlation term, the potential is linearized on a fine mesh and the cycle is repeated until

$V_t(\phi)$20

For the experimentally realized three-well QCD cascade cited in the source, the iteration converges in approximately 25 steps, predicts $V_t(\phi)$21 meV, and agrees with the experimental $V_t(\phi)$22 meV to within 2.5% (Boyko, 2018).

For explicitly time-dependent barriers, a different construction introduces a localized negative imaginary absorber $V_t(\phi)$23 into the Schrödinger equation,

$V_t(\phi)$24

with $V_t(\phi)$25 built from Gamow-like outgoing states. In the current-biased Josephson-junction application, this leads to a TDIP of the form

$V_t(\phi)$26

and a survival probability $V_t(\phi)$27 with instantaneous decay rate

$V_t(\phi)$28

The resulting switching-current distributions reproduce the zero-field switching peak, microwave-induced resonance side peaks, multi-peak structures under chirp, dynamical bifurcation, and real-time oscillatory detection of a weak field (Andersen et al., 2013).

In field emission from curved nanotips, the tunneling potential is the local normal-direction barrier including curvature corrections. For sharply curved emitters the corrected form is

$V_t(\phi)$29

with the same coefficients $V_t(\phi)$30 and $V_t(\phi)$31 for hemiellipsoid, hyperboloid, and sphere. Numerical validation on a rounded conical tip shows the corrected expansion reproduces the exact barrier to better than 1–2% out to $V_t(\phi)$32 nm, and the empirical regime of validity quoted in the source is $V_t(\phi)$33 nm with local field strengths of $V_t(\phi)$34 V/nm or higher (Biswas et al., 2017).

A further adjacent usage appears in singular-interaction models, where tunneling is encoded not by a spatial barrier profile but by exponentially small off-diagonal matrix elements of a finite-dimensional principal matrix $V_t(\phi)$35 in Krein’s formula. In the one-dimensional double-$V_t(\phi)$36 problem,

$V_t(\phi)$37

and degenerate perturbation theory gives

$V_t(\phi)$38

This is not the same formalism as the field-space $V_t(\phi)$39 construction, but it shares the same basic strategy of trading a tunneling boundary-value problem for a more compact auxiliary object (Erman et al., 2019).

Across these usages, the common thread is methodological rather than notational: the tunneling problem is reformulated so that the decisive quantity is obtained from a potential-like auxiliary structure with simpler boundary conditions, a lower-dimensional variational principle, or a more stable numerical scheme. In vacuum decay, that reformulation is particularly sharp, since the tunneling potential $V_t(\phi)$40 replaces the Euclidean bounce by a true minimization problem in field space (Espinosa, 2024).