Gaussian fluctuations in the tunneling probability of a closed universe

Abstract: We consider the quantum creation of a closed universe within the Euclidean path-integral formalism. An analytical expression for the tunneling probability is derived, including both the exponential suppression and the exact Gaussian prefactor due to quadratic fluctuations around the instanton. The calculation is performed in a fixed-interval minisuperspace formulation, where the Hamiltonian constraint is imposed at the level of the classical instanton, while the full lapse integration is not included beyond the leading semiclassical approximation. The result provides a transparent and self-consistent semiclassical estimate of the nucleation rate, refining previous analyses with the inclusion of Gaussian fluctuations.

• The paper presents a semiclassical derivation of an exact analytic Gaussian prefactor that refines the standard tunneling probability in a closed FLRW universe.
• It employs both an approximate harmonic method and an exact approach via the Gel'fand-Yaglom theorem to compute the fluctuation determinant.
• The refined calculation highlights how Gaussian fluctuations significantly affect cosmic nucleation rates and the interpretation of initial conditions.

Analytical Evaluation of Gaussian Fluctuations in the Tunneling Probability of a Closed Universe

Introduction

The paper "Gaussian fluctuations in the tunneling probability of a closed universe" (2605.16150) provides a rigorous semiclassical analysis of the quantum nucleation of a closed Friedmann-Lemaître-Robertson-Walker (FLRW) universe within the Euclidean path-integral formalism. Specifically, it derives a fully analytical expression for the universe creation tunneling probability, not only recovering the standard exponential suppression governed by the instanton action, but also calculating—exactly—the Gaussian prefactor arising from quadratic fluctuations around the instanton. This constitutes a refinement over previous treatments which did not fully address the explicit fluctuation determinant in the closed-FLRW cosmological setting.

Theoretical Framework and Methodology

The quantum nucleation of a closed universe is analyzed as a tunneling event from "nothing," corresponding to a transition from a classically forbidden region (vanishing scale factor) via an instanton trajectory in Euclidean spacetime. The standard minisuperspace action is considered for the scale factor $a(t)$, neglecting inhomogeneities and anisotropies, and reduced to a gauge-fixed path integral. The cosmological constant $\Lambda$ plays a pivotal role, setting both the effective action and the kinematics of the instanton solution.

The Euclidean path-integral amplitude is dominated by a one-dimensional instanton, with the tunneling amplitude given by

$$P_T = \left|\left\langle a_F=\sqrt{3}, t_F | a_I=0, t_I=0 \right\rangle\right|^2,$$

where the endpoints correspond to the forbidden and allowed regions of the scale factor. After Wick rotation and extremization, the instanton profile

$$a_c(\tau) = \sqrt{3} \sin \left(\frac{\tau}{\sqrt{3}}\right),$$

with total Euclidean time interval $\tau_F = \frac{\pi}{2}\sqrt{3}$, provides an exact stationary solution. The leading order action is computed explicitly, yielding the familiar exponential suppression: $$P_T = \exp \left(-\frac{3 \pi c^3}{G \Lambda \hbar}\right).$$ This is in concordance with the Vilenkin tunneling proposal, which contrasts with the Hartle-Hawking prescription via the sign of the exponential.

Calculation of the Gaussian Prefactor

The main technical advance lies in the analytic computation of the Gaussian prefactor, corresponding to the determinant of a nontrivial differential operator arising from quadratic fluctuations about the instanton. The path integral is expanded to leading order in $\delta a(\tau)$, leading to the evaluation of

$$F_E = \mathcal{C} \, \left[\det \left(\frac{1}{\hbar} \frac{\delta^2 S_E}{\delta a(\tau) \delta a(\tau')}\right)\right]^{-1/2}.$$

Two approaches are developed:

1. Approximate (harmonic) method: The potential is linearized near the top of the barrier, allowing the system to be mapped to an inverted harmonic oscillator. The prefactor is then computable in closed form:

$$F_E^{(\text{approx})} = \sqrt{ \frac{\sqrt{3} A}{2\pi \hbar \sinh(\sqrt{3}\pi/2)} }.$$

This leads to

$$P_T^{(\text{approx})} \simeq 0.018 \left( \frac{3\pi c^3}{G\Lambda \hbar}\right) \exp\left( -\frac{3\pi c^3}{G\Lambda\hbar} \right).$$

2. Exact method: The full fluctuation operator is transformed, via an isospectral map, to a Schrödinger-like operator with a nontrivial potential. Application of the Gel'fand-Yaglom theorem allows the determinant ratio (and thus the prefactor) to be reduced to ratios of special functions:

$\Lambda$0

The resulting tunneling probability is

$\Lambda$1

This constitutes a significant correction in the prefactor as compared to the harmonic estimate, increasing the normalization by an order of magnitude but—in the semiclassical limit—leaving the exponential suppression dominant.

Comparison to Previous Work and Consistency

While the exponential scaling of the tunneling probability is a standard result, the explicit analytic evaluation of the fluctuation determinant for the closed FLRW case had not been carried out in the literature. Notably, previous works addressed related conceptual issues, such as time-reparametrization invariance and the conformal factor problem, but either avoided the explicit prefactor or relied on approximations or numerical treatments [pagels, vilenkin, linde, halliwell2019, feldbrugge2017, jia2023].

The current work is careful to situate its calculation within an explicit gauge-fixed minisuperspace context without imposing the full Wheeler-DeWitt constraint at the quantum (fluctuation) level, acknowledging the limitations of the semiclassical approach and the potential for further corrections via fully nonperturbative treatments or Lorentzian contour formulations.

Implications and Future Directions

The results yield an explicit, compact formula for the nucleation rate of a closed universe under the Vilenkin proposal, now including the proper semiclassical normalization. This analytic control is directly relevant for theoretical investigations into initial conditions for cosmic inflation and the probabilistic interpretation of universe creation measures within semiclassical gravity.

The main theoretical implication is that, while the inclusion of the Gaussian prefactor modifies only the normalization and not the exponential hierarchy dictated by the action, precise calculations of creation rates and their dependence on $\Lambda$2 are now possible—including for model comparison between different wave function proposals (Vilenkin vs. Hartle-Hawking). Furthermore, the methodology sets a standard for similar fluctuation analyses in other quantum gravity tunneling contexts, particularly where leading-order approximations are unreliable or insufficient.

Potential future developments include:

• Numerical evaluation of the full path integral or direct solution of the Wheeler-DeWitt equation to investigate non-Gaussian and nonperturbative corrections.
• Extension of analytic techniques to more general minisuperspace models incorporating scalar fields, anisotropies, or alternative gravity theories.
• Comparative studies between Euclidean and Lorentzian quantum cosmological path integrals, especially in view of resurgence ideas and Picard-Lefschetz theory.
• Applications to the statistical prediction of inflationary initial conditions and their dependence on precise measure assignments.

Conclusion

The paper (2605.16150) presents a complete semiclassical path-integral treatment of the quantum creation rate of a closed universe, culminating in an analytic expression for the tunneling probability that includes both exponential suppression and the exact quadratic fluctuation prefactor. The analytic methodology and result fill a longstanding gap in quantum cosmology, enhancing the precision and clarity of competing wave function proposals for the universe's quantum origin. While the semiclassical result is only valid in the regime $\Lambda$3, it provides a critical analytic baseline against which future fully quantum or numerical treatments may be compared.