Dawn and Twilight Time in Quantum Tunneling

Metastable decay exhibits a familiar exponential regime bracketed by early-time deviations and late-time power-law tails. We adopt the real-time, flux-based definition of the decay rate in the spirit of Andreassen et al.\ direct method and present a complete analysis of one-dimensional quantum-mechanical resonance models. We show that the kernel admits a universal pole–plus–branch decomposition and use it to define two computable time scales: a dawn time, when a single resonant contribution starts dominating and exponential decay sets in, and a twilight time, when the branch-cut tail overtakes exponential decay. The latter can be expressed in closed form via the Lambert $W$ function, making its parametric dependence manifest without fitting. For square, modified square, and Pöschl–Teller barriers we obtain simple thick-barrier formulas, clarify the relation $ΓT = T_{\text{trans}}$ between the decay rate $Γ$, oscillation period $T$, and transmission probability $T_{\text{trans}}$, and indicate how our spectral picture can be naturally extended to quantum field theoretic vacuum decay.

💡 Research Summary

The paper presents a comprehensive real‑time analysis of metastable decay in one‑dimensional quantum mechanics, focusing on the onset and termination of the familiar exponential regime. Starting from the continuity equation, the authors define the decay rate Γ as the outward probability flux from a metastable region L (Eq. 3). They then introduce the propagator (kernel) K(t,x;0,y), which encodes the full dynamics independent of the initial state. By performing a contour integral in the complex energy plane (Fig. 2), the kernel is split into two distinct contributions: a sum over isolated complex poles (resonances) and an integral along the branch cut (continuum).

Each pole corresponds to a quasi‑bound state with complex energy Eₙ = Eₙ⁽ᴿ⁾ – iΓₙ/2, giving rise to exponential terms exp(–Γₙ t). The branch‑cut contribution yields a power‑law tail ∝ t^{–γ}, where the exponent γ depends on the barrier shape (γ = 3/2 for a square barrier, γ = 2 for a modified square barrier, and γ = 1/2 for a Pöschl–Teller barrier).

Two characteristic times are defined. The “dawn time” t◒ marks when the lowest‑lying resonance (n = 0) begins to dominate over all higher resonances. Its value depends on the initial wave packet through overlap coefficients c(p) and is estimated by comparing the exponential contributions of different poles (Eqs. 23‑24). The “twilight time” t◓ is the moment when the exponential decay is overtaken by the universal power‑law tail. By equating the magnitude of the leading exponential term with the branch‑cut term, the authors obtain a transcendental equation that can be solved analytically using the Lambert W function. The final closed‑form expression (Eq. 15) is
t◓ = –2γ W_{–1}!\Bigl( –ℏ E₀ G F /