Pith. sign in

REVIEW 3 major objections 6 minor 48 references

Inflation needs a patch about 2.5 times the Hubble scale

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · glm-5.2

2026-07-08 05:35 UTC pith:D2J6TMHC

load-bearing objection Solid numerical study generalizing the Goldwirth-Piran inflation threshold to asymptotically flat spacetimes with inhomogeneous curvature; the proper-volume unification is a genuine new observation but rests on a narrow parameter slice. the 3 major comments →

arxiv 2607.06441 v1 pith:D2J6TMHC submitted 2026-07-07 gr-qc

Starting inflation in asymptotically flat spacetimes

classification gr-qc PACS 98.80.Cq04.25.D-04.20.Ex
keywords inflationnumerical relativityinitial dataasymptotically flat spacetimeconstraint equationse-foldsproper radius thresholdspherical symmetry
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper asks whether a localized blob of inflaton field, sitting in an otherwise empty and asymptotically flat spacetime, can grow into a period of accelerated cosmic expansion. The authors construct spherically symmetric initial data in which the balance between intrinsic and extrinsic curvature is allowed to vary freely, removing the periodic-boundary conditions that prior simulations imposed and which may have biased those studies toward finding inflation. They then evolve dozens of such configurations through full nonlinear general relativity. The central result is that the proper physical size of the initial fluctuation—not the coordinate size, not the split between curvature types, and not the branch of the constraint solution—is the single quantity that determines how many e-folds of inflation result. Below a proper radius of roughly 2.5 times the inflationary length scale, the fluctuation collapses to a black hole rather than inflating; above it, inflation proceeds robustly and approaches the homogeneous limit. This generalizes the classic Goldwirth-Piran threshold from closed, periodic universes to open, asymptotically flat spacetimes with inhomogeneous curvature profiles.

Core claim

When initial data configurations with different mixes of intrinsic and extrinsic curvature, different branches of the constraint equations, and different coordinate radii are all mapped to a single physical measure—the proper radius of the initial fluctuation—they collapse onto approximately one universal curve relating proper size to the number of e-folds of inflation achieved. The threshold for successful inflation sits at a proper radius of about 2.5 times the inflationary length scale, below which black hole formation occurs instead. This universality holds despite the fact that pairs of solutions on different branches can have radically different geometric properties, including one case

What carries the argument

The paper introduces a parameter epsilon that controls what fraction of the inflaton's potential energy density sources the intrinsic curvature (via the conformal factor) versus the extrinsic curvature (via the mean curvature K). For a given fluctuation width, solutions to the Hamiltonian constraint exist only up to a critical epsilon, and below that critical value there are two distinct branches—a strong-field branch with larger spatial volume (sometimes featuring a throat or 'bag of gold' geometry) and a weak-field branch with smaller volume. The key diagnostic is a set of normalized contributions to the Hamiltonian constraint (omega_V for potential energy, omega_R for intrinsic curvature,

Load-bearing premise

The claim that proper volume is the universal driver of inflationary outcomes rests on simulations restricted to spherical symmetry, a simple quadratic inflaton potential, and zero initial conjugate momentum for the scalar field. The authors state they do not expect relaxing these restrictions to change the picture, but non-spherical perturbations, gravitational wave content, or kinetic inhomogeneities could in principle alter the threshold or break the universality of the e-

What would settle it

A 3+1D simulation with the same asymptotically flat boundary conditions but including non-spherical perturbations or nonzero initial momentum, in which two configurations with the same proper radius yield significantly different e-folds of inflation.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Periodic-box simulations of inflation with box sizes near the inflationary scale may systematically overestimate inflation's robustness, because periodicity forces a relationship between average energy density and average expansion that does not hold in genuinely open spacetimes.
  • The minimum proper size threshold of ~2.5 L_infl provides a concrete target for future 3+1D simulations: any computational domain smaller than about five times the inflationary scale risks biasing the outcome toward inflation regardless of the initial fluctuation profile.
  • The observation that strong-field and weak-field branch solutions with the same proper radius yield the same e-folds suggests a kind of geometric universality that could be exploited to reduce the parameter space of future inflation-robustness studies.
  • Black hole formation at the center of sub-threshold fluctuations in asymptotically flat spacetimes provides a channel for primordial black hole production that is absent in periodic simulations where the topology prevents genuine collapse.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. This paper studies whether a localized inflaton fluctuation can seed inflation in an asymptotically flat spacetime, generalizing the Goldwirth-Piran (GP) threshold result. The authors construct spherically symmetric initial data satisfying the Einstein constraint equations, parametrized by the fluctuation width σ and a parameter ε controlling the split between intrinsic and extrinsic curvature sources. They find two branches of solutions (strong-field and weak-field) and evolve them using BSSN in spherical symmetry. The central result is that the proper radius σ_prop ≈ 2.5 L_infl is the universal threshold for successful inflation, regardless of the intrinsic/extrinsic curvature split, and that the number of e-folds correlates primarily with the proper volume of the initial fluctuation. This extends the GP result to open spacetimes with inhomogeneous curvature profiles.

Significance. The paper addresses a well-motivated question about the robustness of inflation to inhomogeneous initial conditions. The key advance over prior work is the use of asymptotically flat boundary conditions, which avoids the potential bias of periodic domains where the average expansion is tied to the average energy density. The finding that proper volume unifies the threshold across different curvature splits is a clean, falsifiable result. The numerical methods are well-established (BSSN in spherical symmetry with a reference-metric approach), and the authors provide convergence tests (Appendix B: second-order for initial data, fourth-order for evolution) and a careful shooting-method solution of the constraint equations (Appendix A). The discussion of how periodic boundary conditions may bias results towards inflation is a valuable contribution to the ongoing debate.

major comments (3)
  1. §IV.C, Fig. 9: The universality claim — that proper volume is the 'primary driver' of e-folds regardless of the curvature split — rests on the convergence of only three curves (ε/ε_crit = 0, 0.5, 1.0) for a single potential V(φ) = ½m²φ², with Π = 0, and in spherical symmetry. The authors acknowledge these restrictions in §V but offer no quantitative argument for why they should be innocuous. In particular, the ε-parametrization in Eq. (9) routes only V(φ) between intrinsic and extrinsic curvature sources, while gradient energy (D_iφ)² always contributes to K via Eq. (9b). This means the parameter space does not explore all possible distributions of energy among geometric terms. The claim would be strengthened by either (a) testing at least one alternative splitting convention to confirm the curves still collapse, or (b) softening the universality language to clearly scope it to the ε-spl
  2. §IV.C, Fig. 9: The threshold σ_prop ≈ 2.5 L_infl is identified visually from the collapse of three curves. No error bars or quantitative measure of the 'roughly single curve' collapse is provided. Given that this threshold is the paper's central quantitative result, a more systematic characterization — e.g., the scatter in σ_prop at fixed e-fold count across the three cases, or a fit with an uncertainty estimate — would make the claim more rigorous.
  3. §V, footnote 3: The authors note that the choice of initially negative K (expansion) is 'a necessary condition for successful inflation that cannot be avoided, and as far as we know it is not motivated by any first principles argument.' This is an important caveat that is somewhat buried in a footnote. Since the entire parameter space explored has K < 0 (or K = 0), the threshold result is conditional on this choice. The paper would benefit from a brief discussion in the main text of how this restriction affects the generality of the claimed threshold, and whether the proper-volume scaling would survive for initial data with mixed-sign K.
minor comments (6)
  1. §II, Eq. (12): The Gaussian profile φ = φ₀ exp(-r²/σ²) has φ → 0 at large r, meaning the asymptotic field value is at the minimum of V(φ) = ½m²φ². This is a specific choice (the field sits at the vacuum rather than on the plateau asymptotically). This should be stated explicitly, as it affects the interpretation of the asymptotically flat boundary conditions.
  2. §III: The inflation criterion ω_V ≥ 0.9 Σ|ω_i| is stated, but the threshold value 0.9 is not justified. How sensitive are the results (particularly σ_min) to this choice? A brief comment on robustness would help.
  3. Fig. 1: The two panels are labeled only by σ values (0.77 L_infl and 10 L_infl) in the caption, but the axis labels and panel titles are not clearly distinguished. Adding explicit labels to each panel would improve readability.
  4. §IV.C: The statement 'for σ ⪅ 2 L_infl a black hole forms at the centre' uses the symbol ⪅ without definition. Standardize notation (≤ or < with approximate qualifier in text).
  5. Appendix D, Eq. (D1): The slicing condition is introduced without much explanation of the term √(24πρ). A brief comment on why this particular combination (K + √(24πρ)) is chosen, and how it relates to the FLRW value of K, would help readers reproduce the results.
  6. References: The paper by Corman & East (2022) [20] and Elley et al. (2024) [21] are cited in the introduction but their specific findings about kinetic inhomogeneities (Π ≠ 0) are not discussed in §V despite being directly relevant to the paper's limitations. A brief comparison would strengthen the discussion.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for a careful and constructive report. The recommendation of minor revision is appropriate, and we agree with all three major comments. We will (1) soften the universality language and scope it to the ε-splitting convention, (2) add a quantitative measure of the curve collapse in Fig. 9, and (3) promote the K < 0 caveat from footnote 3 to the main text. We provide point-by-point responses below.

read point-by-point responses
  1. Referee: §IV.C, Fig. 9: The universality claim rests on only three curves for a single potential with Π=0 in spherical symmetry, and the ε-parametrization does not explore all possible distributions of energy among geometric terms. The claim would be strengthened by (a) testing an alternative splitting convention or (b) softening the universality language.

    Authors: The referee is correct on both counts. We have not tested alternative splitting conventions, and the current ε-parametrization routes only V(φ) between intrinsic and extrinsic curvature sources while gradient energy always contributes to K via Eq. (9b). We agree that this limits the generality of the universality claim as stated. Rather than undertaking new simulations of an alternative splitting (which would be a substantial extension beyond the scope of a minor revision), we will adopt option (b): we will revise the language in §IV.C and §V to scope the universality claim explicitly to the ε-splitting convention employed, noting that gradient energy always sources the extrinsic curvature in our parametrization and that we have not explored alternative partitions. We will also note in §V that testing alternative splitting conventions is a natural direction for future work. revision: yes

  2. Referee: §IV.C, Fig. 9: The threshold σ_prop ≈ 2.5 L_infl is identified visually with no error bars or quantitative measure of the curve collapse. A more systematic characterization would make the central quantitative result more rigorous.

    Authors: We agree. The threshold is currently identified by eye, and a quantitative measure would strengthen the result. In the revised manuscript we will add a quantitative characterization of the curve collapse: specifically, we will report the scatter in σ_prop at a fixed e-fold count (e.g., at N_e = 30, the midpoint of the transition) across the three ε/ε_crit cases, and we will provide a fit to the threshold crossing with an uncertainty estimate derived from the spread between the three curves. This will be included as supplementary quantitative information alongside Fig. 9. revision: yes

  3. Referee: §V, footnote 3: The K < 0 restriction is an important caveat buried in a footnote. Since the entire parameter space has K ≤ 0, the threshold result is conditional on this choice. The paper would benefit from a brief main-text discussion of how this restriction affects generality and whether proper-volume scaling would survive for mixed-sign K.

    Authors: We agree that this caveat deserves more prominence. The restriction to initially negative K (expansion) is indeed a necessary condition for successful inflation in our setup, and the entire threshold result is conditional on it. We will promote this discussion from footnote 3 to the main text of §V, where we will explicitly state that the threshold σ_prop ≈ 2.5 L_infl applies only to initial data with K ≤ 0, and that we have not explored mixed-sign K profiles. We will note that for initial data with converging normal observers (K > 0), inflation cannot begin without a bounce mechanism, and that whether the proper-volume scaling would survive for more general K profiles remains an open question. revision: yes

Circularity Check

0 steps flagged

No significant circularity; self-citation [27] provides initial-data framework but central dynamical result is independently obtained

full rationale

The paper's central claim — that a proper size of ~2.5 L_infl is required for inflation and that proper volume is the unifying driver of e-folds regardless of the intrinsic/extrinsic curvature split — emerges from numerical dynamical evolutions (Figs. 7-9), not from the initial-data construction. The parameter ε parametrizes the split between intrinsic and extrinsic curvature sources in the Hamiltonian constraint (Eq. 9), but nothing in this parametrization forces the e-folds vs. σ_prop curves to collapse onto a single track; that is an empirical observation from the simulation output. The self-citation [27] (Baumgarte, Clough, Giblin — two of three current authors) is invoked for the mathematical framework of the constraint equations: the existence of strong-field/weak-field branches and the maximum intrinsic curvature limit. However, [27] itself generalizes results from [31, 32] (Pfeiffer & York; Baumgarte, Murchadha & Pfeiffer), which are by partially different authors and concern well-posedness properties of elliptic PDEs. The self-citation provides the setup, not the conclusion. The dynamical results (which configurations inflate, how many e-folds result, the threshold value) are obtained by evolving the initial data through Einstein's equations and are not determined by the constraint-solving procedure. The proper radius σ_prop = ∫ψ² dr (Eq. 17) depends on the conformal factor ψ solved from the constraints, but the claim that e-folds align when plotted against σ_prop is a finding about the dynamics, not a tautology of the initial data construction. No step in the derivation chain reduces to its own inputs by construction.

Axiom & Free-Parameter Ledger

3 free parameters · 5 axioms · 0 invented entities

No new physical entities, particles, forces, or dimensions are introduced. The paper works within standard GR plus a canonical scalar field. The free parameters (σ, ϵ, ϕ_0) are initial-data choices, not fitted constants. The initially negative K is a necessary condition for expansion that the authors explicitly flag as a choice without first-principles motivation.

free parameters (3)
  • σ (fluctuation width) = varied over range ~0.77 to ~10 L_infl
    The coordinate width of the Gaussian inflaton profile, varied to find the minimum for successful inflation.
  • ϵ (curvature split fraction) = varied 0 to ϵ_crit
    Parametrizes the proportion of potential energy density sourcing intrinsic vs. extrinsic curvature in the initial data.
  • ϕ_0 (field amplitude) = chosen for 60 e-folds in FLRW
    Fixed so that a homogeneous universe with φ=ϕ_0 would undergo 60 e-folds of inflation.
axioms (5)
  • domain assumption Spherical symmetry of the spacetime and initial data
    Invoked throughout; simplifies the problem to 1+1D and eliminates gravitational wave degrees of freedom. Acknowledged as a limitation in Sec. V.
  • domain assumption Quadratic inflaton potential V(φ) = ½m²φ²
    Stated in Eq. (10); a standard large-field inflation model. Model dependence is discussed but not explored with other potentials.
  • domain assumption Zero initial conjugate momentum Π = 0
    Stated in Eq. (12); imposes moment of time symmetry on the scalar field. Acknowledged as a simplification in Sec. V.
  • domain assumption Initially negative mean curvature K (expansion)
    Discussed in footnote 3 of Sec. V: the negative K solution is a choice necessary for inflation; the alternative (convergence) will always see inflation fail.
  • standard math BSSN formulation correctly evolves Einstein's equations
    Standard numerical relativity formulation, referenced via [33-35]. Validated by convergence tests in Appendix B.

pith-pipeline@v1.1.0-glm · 21169 in / 2762 out tokens · 212131 ms · 2026-07-08T05:35:22.907359+00:00 · methodology

0 comments
read the original abstract

A key question in early universe cosmology is whether inflation can be successfully seeded by a generic, localised fluctuation in the inflaton field that probes the inflationary part of the potential. Past simulations have mainly considered periodic spacetimes representing either a closed universe of a specific size or a typical patch of a larger one, and as a result have needed to impose restrictive conditions on the extrinsic and intrinsic curvatures, which are arguably not generic. In this work we consider initial fluctuations in asymptotically flat spacetimes, allowing more general profiles of the intrinsic and extrinsic curvature. Our findings confirm and generalise the result of Goldwirth and Piran that a fluctuation with proper size several times the inflationary scale $(G\rho_{\rm infl})^{-1/2}$ is required for successful inflation. We also discuss inherent restrictions on the initial data, and how imposing a periodic length close to the inflationary scale may bias results.

Figures

Figures reproduced from arXiv: 2607.06441 by Katy Clough, Sam E. Brady, Thomas W. Baumgarte.

Figure 1
Figure 1. Figure 1: FIG. 1. The initial Misner-Sharp mass ( [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The areal radius [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The areal radius [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Normalised contributions to the Hamiltonian con [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Extrapolated [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Same as Figs [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. The scalar field [PITH_FULL_IMAGE:figures/full_fig_p010_10.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Local values of the Hamiltonian constraint (top [PITH_FULL_IMAGE:figures/full_fig_p012_12.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

48 extracted references · 48 canonical work pages · 26 internal anchors

  1. [1]

    (A1a) forms a non-linear, one-dimensional elliptic equation forψthat we solve subject to the boundary con- ditions∂ rψ= 0 atr= 0, which is required by regularity, andψ→1 asr→ ∞

    Hamiltonian constraint Eq. (A1a) forms a non-linear, one-dimensional elliptic equation forψthat we solve subject to the boundary con- ditions∂ rψ= 0 atr= 0, which is required by regularity, andψ→1 asr→ ∞. One possible approach for solving this equation is to linearize and iterate. Specifically, we could write ψn+1 =ψ n +u ,(A2) 11 whereuis the correction ...

  2. [2]

    As described above, we iterate between solving Eq

    Momentum constraint The extrinsic curvatureK ij can be reconstructed from the mean curvatureK, the (flat) conformally related spa- tial metric ¯γij, the conformal factorψ, and the confor- mally related traceless part ¯Aij. As described above, we iterate between solving Eq. (A1b) for the mean curvature Kand Eq. (A1c) for ¯Aij. We solve the momentum constra...

  3. [3]

    A. H. Guth, The Inflationary Universe: A Possible Solu- tion to the Horizon and Flatness Problems, Phys. Rev. D23, 347 (1981)

  4. [4]

    A. A. Starobinsky, A New Type of Isotropic Cosmological Models Without Singularity, Phys. Lett. B91, 99 (1980)

  5. [5]

    A. D. Linde, A New Inflationary Universe Scenario: A Possible Solution of the Horizon, Flatness, Homogeneity, Isotropy and Primordial Monopole Problems, Phys. Lett. B108, 389 (1982)

  6. [6]

    Albrecht and P

    A. Albrecht and P. J. Steinhardt, Cosmology for Grand Unified Theories with Radiatively Induced Symmetry Breaking, Phys. Rev. Lett.48, 1220 (1982)

  7. [7]

    A. D. Linde, Chaotic Inflation, Phys. Lett. B129, 177 (1983)

  8. [8]

    Initial Conditions for Inflation - A Short Review

    R. Brandenberger, Initial conditions for inflation — A short review, Int. J. Mod. Phys. D26, 1740002 (2016), arXiv:1601.01918 [hep-th]

  9. [9]

    D. S. Goldwirth and T. Piran, Inhomogeneity and the Onset of Inflation, Phys. Rev. Lett.64, 2852 (1990)

  10. [10]

    D. S. Goldwirth, On inhomogeneous initial conditions for inflation, Phys. Rev. D43, 3204 (1991)

  11. [11]

    D. S. Goldwirth and T. Piran, Spherical Inhomogeneous Cosmologies and Inflation. 1. Numerical Methods, Phys. Rev. D40, 3263 (1989)

  12. [12]

    D. S. Goldwirth and T. Piran, Initial conditions for in- flation, Phys. Rept.214, 223 (1992)

  13. [13]

    Kurki-Suonio, R

    H. Kurki-Suonio, R. A. Matzner, J. Centrella, and J. R. Wilson, Inflation From Inhomogeneous Initial Data in a One-dimensional Back Reacting Cosmology, Phys. Rev. D35, 435 (1987)

  14. [14]

    Laguna, H

    P. Laguna, H. Kurki-Suonio, and R. A. Matzner, Inho- mogeneous inflation: The Initial value problem, Phys. Rev. D44, 3077 (1991)

  15. [15]

    Inhomogeneous Inflation: Numerical Evolution

    H. Kurki-Suonio, P. Laguna, and R. A. Matzner, Inho- mogeneous inflation: Numerical evolution, Phys. Rev. D 48, 3611 (1993), arXiv:astro-ph/9306009

  16. [16]

    Shibata, K.-i

    M. Shibata, K.-i. Nakao, T. Nakamura, and K.-i. Maeda, Dynamical evolution of gravitational waves in the asymp- totically de Sitter space-time, Phys. Rev. D50, 708 (1994)

  17. [17]

    W. E. East, M. Kleban, A. Linde, and L. Senatore, Be- ginning inflation in an inhomogeneous universe, JCAP 14 09, 010, arXiv:1511.05143 [hep-th]

  18. [18]

    Robustness of Inflation to Inhomogeneous Initial Conditions

    K. Clough, E. A. Lim, B. S. DiNunno, W. Fischler, R. Flauger, and S. Paban, Robustness of Inflation to Inhomogeneous Initial Conditions, JCAP09, 025, arXiv:1608.04408 [hep-th]

  19. [19]

    Robustness of Inflation to Large Tensor Perturbations

    K. Clough, R. Flauger, and E. A. Lim, Robustness of Inflation to Large Tensor Perturbations, JCAP05, 065, arXiv:1712.07352 [hep-th]

  20. [20]

    J. C. Aurrekoetxea, K. Clough, R. Flauger, and E. A. Lim, The Effects of Potential Shape on Inhomoge- neous Inflation, JCAP05, 030, arXiv:1910.12547 [astro- ph.CO]

  21. [21]

    Inhomogeneous initial conditions for inflation: A wibbly-wobbly timey-wimey path to salvation

    C. Joana and S. Clesse, Inhomogeneous preinflation across Hubble scales in full general relativity, Phys. Rev. D103, 083501 (2021), arXiv:2011.12190 [astro-ph.CO]

  22. [22]

    Starting inflation from inhomogeneous initial conditions with momentum

    M. Corman and W. E. East, Starting inflation from in- homogeneous initial conditions with momentum, JCAP 10, 046, arXiv:2212.04479 [gr-qc]

  23. [23]

    Robustness of inflation to kinetic inhomogeneities

    M. Elley, J. C. Aurrekoetxea, K. Clough, R. Flauger, P. Giannadakis, and E. A. Lim, Robustness of inflation to kinetic inhomogeneities, (2024), arXiv:2405.03490 [astro- ph.CO]

  24. [24]

    Beginning inflation in conformally curved spacetimes

    C. Joana, Beginning inflation in conformally curved spacetimes, Phys. Rev. D110, 063534 (2024), arXiv:2406.00811 [astro-ph.CO]

  25. [25]

    Initial conditions problem in cosmological inflation revisited

    D. Garfinkle, A. Ijjas, and P. J. Steinhardt, Initial con- ditions problem in cosmological inflation revisited, Phys. Lett. B843, 138028 (2023), arXiv:2304.12150 [gr-qc]

  26. [26]

    Smoothing and flattening the universe through slow contraction versus inflation

    A. Ijjas, P. J. Steinhardt, D. Garfinkle, and W. G. Cook, Smoothing and flattening the universe through slow contraction versus inflation, JCAP07, 077, arXiv:2404.00867 [gr-qc]

  27. [27]

    Giannadakis, M

    P. Giannadakis, M. Elley, R. Flauger, and E. A. Lim, A critical value of the inflationary tensor-to-scalar ratio from inhomogeneous inflation, JCAP06, 093, arXiv:2512.13673 [astro-ph.CO]

  28. [28]

    Gravitational waves and large field inflation

    A. Linde, Gravitational waves and large field inflation, JCAP02, 006, arXiv:1612.00020 [astro-ph.CO]

  29. [29]

    T. W. Baumgarte, K. Clough, and J. T. Giblin, Jr., Restrictions on initial conditions in cosmological scenar- ios and implications for simulations of primordial black holes and inflation, Phys. Rev. D112, 123528 (2025), arXiv:2509.26470 [gr-qc]

  30. [30]

    S. K. Blau, E. I. Guendelman, and A. H. Guth, The Dy- namics of False Vacuum Bubbles, Phys. Rev. D35, 1747 (1987)

  31. [31]

    Farhi and A

    E. Farhi and A. H. Guth, An Obstacle to Creating a Uni- verse in the Laboratory, Phys. Lett. B183, 149 (1987)

  32. [32]

    Farhi, A

    E. Farhi, A. H. Guth, and J. Guven, Is It Possible to Create a Universe in the Laboratory by Quantum Tun- neling?, Nucl. Phys. B339, 417 (1990)

  33. [33]

    H. P. Pfeiffer and J. W. York, Jr., Uniqueness and non- uniqueness in the Einstein constraints, Phys. Rev. Lett. 95, 091101 (2005), arXiv:gr-qc/0504142

  34. [34]

    T. W. Baumgarte, N. O. Murchadha, and H. P. Pfeif- fer, The Einstein constraints: Uniqueness and non- uniqueness in the conformal thin sandwich approach, Phys. Rev. D75, 044009 (2007), arXiv:gr-qc/0610120

  35. [35]

    Nakamura, K

    T. Nakamura, K. Oohara, and Y. Kojima, General Rela- tivistic Collapse to Black Holes and Gravitational Waves from Black Holes, Prog. Theor. Phys. Suppl.90, 1 (1987)

  36. [36]

    Shibata and T

    M. Shibata and T. Nakamura, Evolution of three- dimensional gravitational waves: Harmonic slicing case, Phys. Rev. D52, 5428 (1995)

  37. [37]

    T. W. Baumgarte and S. L. Shapiro, On the numerical integration of Einstein’s field equations, Phys. Rev. D59, 024007 (1998), arXiv:gr-qc/9810065

  38. [38]

    Deriving formulations for numerical computation of binary neutron stars in quasicircular orbits

    M. Shibata, K. Uryu, and J. L. Friedman, Deriv- ing formulations for numerical computation of binary neutron stars in quasicircular orbits, Phys. Rev. D 70, 044044 (2004), [Erratum: Phys.Rev.D 70, 129901 (2004)], arXiv:gr-qc/0407036

  39. [39]

    A constrained scheme for Einstein equations based on Dirac gauge and spherical coordinates

    S. Bonazzola, E. Gourgoulhon, P. Grandclement, and J. Novak, A Constrained scheme for Einstein equations based on Dirac gauge and spherical coordinates, Phys. Rev. D70, 104007 (2004), arXiv:gr-qc/0307082

  40. [40]

    J. D. Brown, Covariant formulations of BSSN and the standard gauge, Phys. Rev. D79, 104029 (2009), arXiv:0902.3652 [gr-qc]

  41. [41]

    Gourgoulhon, 3+1 Formalism in General Relativity, Lecture Notes in Physics (Springer, 2012)

    E. Gourgoulhon, 3+1 Formalism in General Relativity, Lecture Notes in Physics (Springer, 2012)

  42. [42]

    T. W. Baumgarte, P. J. Montero, I. Cordero-Carrion, and E. Muller, Numerical Relativity in Spherical Polar Coor- dinates: Evolution Calculations with the BSSN Formu- lation, Phys. Rev. D87, 044026 (2013), arXiv:1211.6632 [gr-qc]

  43. [43]

    T. W. Baumgarte, P. J. Montero, and E. M¨ uller, Nu- merical Relativity in Spherical Polar Coordinates: Off- center Simulations, Phys. Rev. D91, 064035 (2015), arXiv:1501.05259 [gr-qc]

  44. [44]

    C. Bona, J. Masso, E. Seidel, and J. Stela, A New for- malism for numerical relativity, Phys. Rev. Lett.75, 600 (1995), arXiv:gr-qc/9412071

  45. [45]

    M. Kopp, S. Hofmann, and J. Weller, Separate Uni- verses Do Not Constrain Primordial Black Hole Forma- tion, Phys. Rev. D83, 124025 (2011), arXiv:1012.4369 [astro-ph.CO]

  46. [46]

    T. W. Baumgarte, Aspherical deformations of the Choptuik spacetime, Phys. Rev. D98, 084012 (2018), arXiv:1807.10342 [gr-qc]

  47. [47]

    J. C. Aurrekoetxea, K. Clough, and E. A. Lim, Cosmol- ogy using numerical relativity, Living Rev. Rel.28, 5 (2025), arXiv:2409.01939 [gr-qc]

  48. [48]

    Growth of a Black Hole in a Scalar Field Cosmology

    J. Doherty, M. Gracia-Linares, and P. Laguna, Growth of a black hole in a scalar field cosmology, Phys. Rev. D 112, 044029 (2025), arXiv:2504.20272 [gr-qc]