Pith. sign in

REVIEW 4 major objections 5 minor 6 cited by

Simulations indicate that the curvature perturbation from hybrid inflation can carry a global topological imprint only when the waterfall sector has a single field.

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 · deepseek-v4-flash

2026-08-04 05:52 UTC pith:Q5DHJO5Z

load-bearing objection Solid STOLAS power-spectrum work and a genuinely new ζ-topology diagnostic, but the abstract's 'sub-Hubble defect reconnection' claim is unsupported by any measured correlation length. the 4 major comments →

arxiv 2603.04850 v2 pith:Q5DHJO5Z submitted 2026-03-05 astro-ph.CO gr-qchep-phhep-th

STOchastic LAttice Simulation of hybrid inflation

classification astro-ph.CO gr-qchep-phhep-th PACS 98.80.Cq
keywords hybrid inflationstochastic inflationlattice simulationcurvature perturbationtopological defectsEuler characteristicprimordial black holeswaterfall transition
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 uses a lattice simulation of stochastic inflation to ask what the curvature perturbation looks like in multi-waterfall hybrid inflation, and whether topological defects leave a mark on it. The authors try to establish two things: that stochastic noise continuously reconnects the topological defects — domain walls (n=1), cosmic strings (n=2), monopoles (n=3) — as soon as they form, so their correlation length at the waterfall critical point ends up much smaller than the Hubble scale; and that only the single-waterfall n=1 case imprints a global, connected structure on the curvature perturbation, detected as a negative Euler characteristic. The simulation also validates the stochastic-δN algorithm by matching the analytic power-spectrum fitting formula, and reproduces the cubic-potential upper bound on δN that suppresses primordial black hole formation. If these claims are right, the curvature perturbation itself becomes a potential probe of topological-defect physics from the early universe, beyond the power spectrum.

Core claim

The paper's central claim is that topological defects produced at the waterfall transition of hybrid inflation are not left as a few large objects per Hubble patch; the stochastic white noise that continuously kicks the super-Hubble fields reconnects the defect network during inflation, fragmenting walls, strings, and monopoles into fine structures with correlation lengths much smaller than the Hubble scale at the critical point. This is shown through snapshots of the radial waterfall field below a threshold and through the time evolution of the Euler characteristic χ = V - E + F of the simulated grid: for n=1 and n=2, χ first dips negative at N ~ 4-5 (holes in domain walls, loops of cosmic

What carries the argument

The central object is STOLAS, a lattice implementation of the stochastic formalism of inflation in which each grid point evolves by a Langevin equation with a spatially correlated white-noise term representing quantum modes crossing the Hubble horizon, and the curvature perturbation ζ(x) is computed as the stochastic fluctuation δN(x) in the first-passage e-folding number to the end surface, averaged over 20 uncorrelated-noise realizations per point. The topological diagnostic is the Euler characteristic χ = V - E + F, counted over lattice cubes whose radial waterfall field is below hand-set thresholds ψ_r,th, with the convention that χ < 2 signals holes or global connected structures and χ

Load-bearing premise

The entire defect-reconnection conclusion rests on identifying 'defects' as the 256^3 grid cells whose radial waterfall field falls below a hand-set threshold (Eq. 5.1), so if that threshold or the lattice resolution is changed, the apparent fragmentation of defects into sub-Hubble pieces could disappear.

What would settle it

Rerun the n=1, n=2, and n=3 simulations with N_L = 512 and 1024 and with ψ_r,th varied by a factor of two. If the Euler characteristic no longer jumps to large positive values after N ~ 4-5, or if the defect correlation length stays comparable to the Hubble scale instead of shrinking, the stochastic-reconnection claim is refuted.

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

If this is right

  • If the reconnection picture holds, defect networks in mild-waterfall hybrid inflation are much finer and more numerous at the end of inflation than the naive one-defect-per-Hubble-patch expectation, so any observable signature of these defects would come from many small objects rather than a few large ones.
  • The n=1 case singles itself out: the curvature perturbation has a global connected structure at thresholds around σζ, while n=2 and n=3 do not, making ζ a possible observable carrier of domain-wall-related information even if the walls themselves are unresolved.
  • The cubic-potential upper bound on δN means primordial black holes should rarely form in these models, while the nearly-uniform saturated regions in the density map may affect halo formation and induced gravitational waves.
  • The agreement between the STOLAS power spectra and the analytic fitting formula validates the stochastic-δN algorithm beyond the slow-roll approximation, explaining the earlier PDF discrepancy as a coarse-graining effect.
  • Because the curvature maps contain all spatial-correlation information, the same simulation pipeline can be used to compute bispectra and higher-order correlators that the stochastic-δN algorithm cannot yet provide.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • I would read the defect-fragmentation result cautiously until a resolution study is done: with N_L=256 and a hand-chosen threshold, the apparent shredding into sub-Hubble pieces could be a lattice artifact, and a convergence run at N_L=512 or 1024 would distinguish a physical reconnection process from discretization noise.
  • If the n=1 global structure in ζ is real, it should show up in the bispectrum as a distinctive squeezed-limit signature; computing the connected three-point function from the STOLAS maps would be a direct test of whether the topological imprint extends beyond the Euler characteristic.
  • The coarse-graining dependence of the PDF noted by the authors has a practical consequence: PBH abundance estimates from hybrid inflation must specify the smoothing scale explicitly, because the same model can show a heavy tail on one scale and a truncated tail on another.
  • The Euler-characteristic method could be extended to measure loop statistics of reconnected strings in the n=2 case and two-point clustering of monopoles in n=3, giving quantitative correlation lengths that are currently left for future work.

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

4 major / 5 minor

Summary. The authors extend their lattice stochastic-inflation code STOLAS to mild-waterfall hybrid inflation with n=1,2,3,15 waterfall fields and 'Quadratic'/'Cubic' inflaton potentials. They generate 256^3 maps of the waterfall fields and the curvature perturbation ζ, compare one-point PDFs and Fourier power spectra with the stochastic-δN fitting formula of Ref. [48], and use the Euler characteristic of thresholded regions to study the topology of waterfall-field defects and of ζ. They report that domain walls (n=1), cosmic strings (n=2), and monopoles (n=3) are reconnected into finer structures whose correlation lengths are much smaller than the Hubble scale, and that only n=1 leaves global topological structures in ζ. They also reproduce the Cubic upper bound of Ref. [54] and conclude that the narrow-sense stochastic-δN algorithm is validated. The paper itself notes that the PDF tail for δN≳1 disagrees with Ref. [54] and that a detailed correlation-length study is left to future work.

Significance. The paper makes a useful technical advance: it applies the public STOLAS code to multi-waterfall hybrid inflation, produces 256^3 real-space maps for six models, and introduces the Euler characteristic as a topological diagnostic for both waterfall fields and ζ. The consistency for the Quadratic n=1 power spectrum with the independent Monte-Carlo calculation of Ref. [57] is a genuine strength, and the reproduction of the Cubic upper bound is a nontrivial cross-check. If the claimed defect-reconnection phenomenon were established, it would change expectations for topological-defect production in hybrid inflation and could provide a new probe through the n=1 curvature perturbation. As it stands, however, the topological evidence is not yet convincing: the paper lacks a correlation-length measurement and resolution/threshold convergence studies, and the 'validation' of stochastic-δN partly relies on a two-parameter fit to the authors' own fitting formula. The significance of the paper is therefore conditional on the additional analysis outlined in the major comments.

major comments (4)
  1. [Abstract; Secs. 5.1, 5.2, 6] The headline statement that topological defects are 'reconnected ... into finer structures ... correlation lengths much smaller than the Hubble scale' is not supported by the quantitative analysis. The diagnostic is the Euler characteristic of voxels with ψ_r<ψ_r,th (Eq. (5.1), thresholds set by hand) on a 256^3 lattice; no two-point function or correlation length is computed, and Sec. 6 explicitly leaves this to future work. With Eq. (2.4), the physical grid spacing is (2π e^N)/(256 H); at N=5 it is ≈3.6/H and at N=7 ≈27/H, so sub-Hubble structures at late times are not resolvable. The rise of χ to large positive values (Fig. 5) is the expected signature of unresolved point-like clusters; Sec. 5.2 itself calls the late-time objects unresolved. Without N_L/threshold convergence and a direct correlation-length measurement, the reconnection claim can be a discretization artifact. It should
  2. [Sec. 4, Eq. (3.6), Table 2] The conclusion that the narrow-sense stochastic-δN algorithm is 'validated by STOLAS' is stronger than the evidence. The comparison fits Eq. (3.6) with N_water and P_ζ^(peak) as free parameters per model (Table 2); a two-parameter fit cannot independently test the algorithm's predicted amplitude or normalization. The independent calculation in Ref. [57] for Quadratic n=1 is a genuine supporting check, but for a validation claim the authors should either compute P_ζ^(peak) a priori or provide parameter uncertainties and a goodness-of-fit statistic. Also, each model uses one lattice realization (Fig. 3 has no error bars). Please quantify the fit and soften the wording.
  3. [Sec. 4, Fig. 2; Sec. 6] The authors acknowledge that the PDF at δN≳1 contradicts Ref. [54] and attribute this to the lattice-scale coarse-graining, but they do not demonstrate this: no resolution study or simulation to the end of inflation is presented, and they only 'expect' to reproduce the earlier tail. Because the heavy tail is central to PBH formation and to the claimed consistency of the stochastic-δN framework, this unresolved discrepancy should be presented as an open issue. A convergence test varying N_L and the number of noise realizations (currently 20) is needed before the PDF comparison can be called 'broadly consistent.'
  4. [Sec. 5.3, Fig. 6] The n=1 'global structures' conclusion rests on χ/N becoming negative for ζ_th∼σ_ζ in one 64^3 subregion, and the authors themselves state that they do not conclude whether it is due to the domain wall or the PDF shape. No error bars, bootstrap, or dependence on the chosen subregion and threshold are shown. As presented, a single negative χ in one realization is not a statistically supported global structure. Please add a significance estimate or, failing that, present this as a tentative hint rather than a main conclusion of the abstract.
minor comments (5)
  1. [Eq. (2.12)] As typeset the two expectation values are identical; if the second is the simulation-box average, please use an overbar or distinct notation.
  2. [Appendix A] The text gives χ=8−12−6=2 for a cube; it should be 8−12+6=2.
  3. [Fig. 5] The caption says orange dots show negative values, but the ordinate is logarithmic; clarify whether |χ| is plotted or how negative values are represented.
  4. [Sec. 2.2] The choice of 20 noise realizations per grid point is not justified; a convergence check in this number would strengthen the coarse-graining procedure.
  5. [Sec. 5.1] The statement that the structure 'hardly depends' on the potential is based on a single realization per model; please state this limitation. Also, the animation link should be a persistent URL.

Circularity Check

1 steps flagged

Power-spectrum 'validation' of stochastic-δN is a two-parameter fit to the authors' own formula; the topological analysis is independent, though its headline correlation-length claim is under-supported rather than circular.

specific steps
  1. fitted input called prediction [Sec. 4, Fig. 3 and Table 2; Eq. (3.6)]
    "T able 2: Simulation results on the average and variance of N, and the fitting parameters for the analytic formula (3.6). ... The black dashed line shows the analytic formula (3.6) with fitting parameters (N water,P (peak) ζ ), whose values are listed in Table 2. ... One finds that the fitting formula works well."

    The comparison that underlies the claim that the results are 'broadly consistent with the stochastic-δN algorithm' uses Eq. (3.6), the stochastic-δN fitting formula, with its two free constants N_water and P_peak^(ζ) adjusted to the STOLAS spectrum. The normalization and peak position are therefore not predicted by the algorithm but imposed by the fit; what remains is only the Gaussian-in-N_k shape. Hence the consistency check is partly a fitted-input re-description rather than an independent confirmation of the algorithm's predictions.

full rationale

The paper's headline defect-reconnection claim is not circular: it is an interpretation of Euler-characteristic counts on thresholded waterfall-field maps (Eq. 5.1, Sec. 5.2), with limitations acknowledged in Sec. 6 ('We leave a detailed study about the correlation length ... for future work'). That claim may be resolution-limited or threshold-dependent, but that is a correctness concern, not a reduction of the result to its inputs. The Cubic upper bound is also observed in the STOLAS PDFs, not merely cited from Ref. [54]. The main circularity-adjacent step is the power-spectrum validation: the stochastic-δN consistency is established by fitting the two parameters of the authors' own formula (Eq. 3.6, Ref. [48]) to the STOLAS spectrum, and the cited prior confirmations (Refs. [52,54,57]) are from the same author group. Since the central new topological analysis is independent and externally checkable, the overall circularity burden is modest.

Axiom & Free-Parameter Ledger

4 free parameters · 6 axioms · 0 invented entities

The model parameters (Lambda, M, phi_c, mu_1, mu_2, mu_3) are fixed by CMB constraints from prior literature; the free parameters listed here are the ones the paper itself fits or chooses by hand to reach its claims. The simulation relies on the standard stochastic-inflation and deltaN assumptions, plus an ad hoc threshold-based identification of defects and a numerical averaging procedure. No new entities are introduced.

free parameters (4)
  • Nwater (waterfall e-fold scale in power spectrum fit) = Quadratic n=1: 18.4; see Table 2 for all models
    Free parameter in the fitting formula (3.6) adjusted to match the STOLAS power spectrum.
  • P_zeta^(peak) (peak amplitude of power spectrum fit) = Quadratic n=1: 6.27e-2; see Table 2 for all models
    Free normalization in Eq. (3.6) fitted to each simulation.
  • psi_r,th (waterfall defect threshold) = 0.05 sigma_psi (n=1), 0.1 sigma_psi (n=2), 0.5 sigma_psi (n=3), 2.5 sigma_psi (n=15)
    Hand-chosen thresholds in Eq. (5.1) defining what counts as a topological-defect site; interpretation of defect structure depends on this choice.
  • sigma (coarse-graining ratio) = 1/16
    Simulation setup Eq. (2.3) fixing the noise correlation scale; results may depend on this choice, and no convergence study over sigma is given.
axioms (6)
  • domain assumption Stochastic Langevin equations (2.1) with noise amplitude P^{1/2}=H/(2pi) delta_ij and neglect of momentum and inflaton noise describe the super-Hubble dynamics.
    Standard stochastic-inflation approximation; adopted without derivation in Sec. 2.1.
  • domain assumption Curvature perturbation zeta equals the fluctuation in the first-passage e-fold number N (deltaN formalism).
    Invoked in Sec. 2.2 as the basis for computing zeta from the lattice.
  • domain assumption The surface eta_r=-2 approximates the end of inflation / uniform-density hypersurface.
    Defined in Sec. 3 (Eq. 3.3) and used as the stopping condition for first-passage times; affects all zeta values.
  • ad hoc to paper Continuing with uncorrelated noise for 20 realizations per grid point yields the lattice-scale coarse-grained zeta_c.
    Procedure in Sec. 2.2 (Eq. 2.12) introduced for this paper; the number 20 is a numerical choice with no stated convergence test.
  • ad hoc to paper Topological-defect content is captured by the threshold psi_r < psi_r,th and the Mathematica vertex/edge/face Euler-characteristic convention.
    Thresholds chosen by hand (Eq. 5.1); counting convention described in Appendix A; quantitative conclusions about defect sizes depend on these choices.
  • domain assumption In the Cubic model, parameter choices are retained only if they satisfy CMB constraints (3.4) after solving the background dynamics.
    A posteriori selection of parameters, Sec. 3; a selection but standard for model-building against CMB.

pith-pipeline@v1.3.0-alltime-deepseek · 15542 in / 14356 out tokens · 135569 ms · 2026-08-04T05:52:37.313120+00:00 · methodology

0 comments
read the original abstract

We investigate the spatial profile of the curvature perturbation generated in multi-waterfall hybrid inflation models, which are known to produce various topological defects. Using the lattice simulation code STOchastic LAttice Simulation, based on the stochastic formalism of inflation, we analyse six cases by varying the number of waterfall fields $n$ and the functional form of the inflaton potential (``Quadratic'' and ``Cubic'' cases). Our statistical analysis shows that the probability density functions (PDFs) and power spectra are broadly consistent with the so-called stochastic-$\delta N$ algorithm. The ``Cubic'' case also exhibits a characteristic upper bound in the PDF, as discovered in our previous work, that suppresses \acl{PBH} formation while potentially affecting halo formation. Furthermore, we employ the Euler characteristic as a topological diagnostic tool to identify the structures of the waterfall fields as well as the curvature perturbation. We find that the topological defects, such as domain walls ($n=1$), cosmic strings ($n=2$), and monopoles ($n=3$), are reconnected during inflation into finer structures by the stochastic noise, making their correlation lengths much smaller than the Hubble scale at the critical point of the waterfall phase transition counterintuitively. The Euler characteristic also implies global structures of the curvature perturbation for $n=1$, though we do not conclude if they are due to the domain wall, because neither the strings ($n=2$) nor monopoles ($n=3$) leave such structures. The global structures of the curvature perturbation will provide a novel probe for the physics of the early universe.

discussion (0)

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

Forward citations

Cited by 6 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Nonlinear Lattice Framework for Inflation: Bridging stochastic inflation and the $\delta{N}$ formalism

    gr-qc 2026-04 unverdicted novelty 8.0

    A shear-free lattice method bridges stochastic inflation and δN formalism by enabling fully nonlinear calculations of curvature perturbations in single-field models with ultra-slow-roll phases.

  2. Stochastic scalar-tensor inflation and beyond

    gr-qc 2026-06 unverdicted novelty 7.0

    Extends the stochastic inflation formalism to a wide class of scalar-tensor theories by mapping EFT of dark energy coefficients to stochastic equations of motion.

  3. A consistent formulation of stochastic inflation I: Non-Markovian effects and issues beyond linear perturbations

    astro-ph.CO 2026-05 unverdicted novelty 7.0

    The conventional truncation in stochastic inflation is inconsistent because quadratic-noise contributions are the same perturbative order as the deterministic non-Markovian corrections.

  4. Quantum algorithm for solving high-dimensional linear stochastic differential equations via amplitude encoding of the noise term

    quant-ph 2026-04 unverdicted novelty 7.0

    Quantum algorithms achieve polylog(N) complexity for high-dimensional linear SDEs by amplitude-encoding the solution and noise via Dyson series or Euler-Maruyama approximations plus quantum linear systems solvers.

  5. Nonperturbative stochastic inflation in perturbative dynamical background

    astro-ph.CO 2026-04 unverdicted novelty 6.0

    Derives stochastic equations from Schwinger-Keldysh formalism that include quantum diffusion and classical metric perturbations for non-perturbative ultra-slow-roll inflation, validated on Starobinsky and critical Hig...

  6. Nonlinear Lattice Framework for Inflation: Bridging stochastic inflation and the $\delta{N}$ formalism

    gr-qc 2026-04 conditional novelty 6.0

    A shear-free locally FLRW lattice framework for single-field inflation captures spatially varying expansion, curvature corrections, and nonlinear δN observables at a fraction of the cost of full numerical relativity.

Reference graph

Works this paper leans on

65 extracted references · 42 linked inside Pith · cited by 5 Pith papers

  1. [1]

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

  2. [2]

    Sato,First-order phase transition of a vacuum and the expansion of the Universe,Mon

    K. Sato,First-order phase transition of a vacuum and the expansion of the Universe,Mon. Not. Roy. Astron. Soc.195(1981) 467

  3. [3]

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

  4. [4]

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

  5. [5]

    Albrecht and P

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

  6. [6]

    A. D. Linde,Chaotic Inflation,Phys. Lett. B129(1983) 177. [7]Planckcollaboration, Y. Akrami et al.,Planck 2018 results. X. Constraints on inflation, Astron. Astrophys.641(2020) A10 [1807.06211]

  7. [8]

    Hawking,Gravitationally collapsed objects of very low mass,Mon

    S. Hawking,Gravitationally collapsed objects of very low mass,Mon. Not. Roy. Astron. Soc. 152(1971) 75

  8. [9]

    B. J. Carr and S. W. Hawking,Black holes in the early Universe,Mon. Not. Roy. Astron. Soc. 168(1974) 399

  9. [10]

    B. J. Carr,The Primordial black hole mass spectrum,Astrophys. J.201(1975) 1

  10. [11]

    B. Carr, K. Kohri, Y. Sendouda and J. Yokoyama,Constraints on primordial black holes,Rept. Prog. Phys.84(2021) 116902 [2002.12778]

  11. [12]

    Escriv` a, F

    A. Escriv` a, F. Kuhnel and Y. Tada,Primordial Black Holes,2211.05767

  12. [13]

    Yoo,The Basics of Primordial Black Hole Formation and Abundance Estimation, Galaxies10(2022) 112 [2211.13512]

    C.-M. Yoo,The Basics of Primordial Black Hole Formation and Abundance Estimation, Galaxies10(2022) 112 [2211.13512]

  13. [14]

    B. Carr, S. Clesse, J. Garcia-Bellido, M. Hawkins and F. Kuhnel,Observational evidence for primordial black holes: A positivist perspective,Phys. Rept.1054(2024) 1 [2306.03903]

  14. [15]

    A. D. Linde,Hybrid inflation,Phys. Rev. D49(1994) 748 [astro-ph/9307002]

  15. [16]

    Garcia-Bellido, A

    J. Garcia-Bellido, A. D. Linde and D. Wands,Density perturbations and black hole formation in hybrid inflation,Phys. Rev. D54(1996) 6040 [astro-ph/9605094]

  16. [17]

    Clesse,Hybrid inflation along waterfall trajectories,Phys

    S. Clesse,Hybrid inflation along waterfall trajectories,Phys. Rev. D83(2011) 063518 [1006.4522]. – 17 –

  17. [18]

    D. H. Lyth,Contribution of the hybrid inflation waterfall to the primordial curvature perturbation,JCAP07(2011) 035 [1012.4617]

  18. [19]

    Kodama, K

    H. Kodama, K. Kohri and K. Nakayama,On the waterfall behavior in hybrid inflation,Prog. Theor. Phys.126(2011) 331 [1102.5612]

  19. [20]

    Mulryne, S

    D. Mulryne, S. Orani and A. Rajantie,Non-Gaussianity from the hybrid potential,Phys. Rev. D84(2011) 123527 [1107.4739]

  20. [21]

    Bugaev and P

    E. Bugaev and P. Klimai,Curvature perturbation spectra from waterfall transition, black hole constraints and non-Gaussianity,JCAP11(2011) 028 [1107.3754]

  21. [22]

    Bugaev and P

    E. Bugaev and P. Klimai,Formation of primordial black holes from non-Gaussian perturbations produced in a waterfall transition,Phys. Rev. D85(2012) 103504 [1112.5601]

  22. [23]

    D. H. Lyth,The hybrid inflation waterfall and the primordial curvature perturbation,JCAP05 (2012) 022 [1201.4312]

  23. [24]

    Clesse and B

    S. Clesse and B. Garbrecht,Slow Roll during the Waterfall Regime: The Small Coupling Window for SUSY Hybrid Inflation,Phys. Rev. D86(2012) 023525 [1204.3540]

  24. [25]

    A. H. Guth and E. I. Sfakianakis,Density Perturbations in Hybrid Inflation Using a Free Field Theory Time-Delay Approach,1210.8128

  25. [26]

    Clesse and J

    S. Clesse and J. Garc ´ ıa-Bellido,Massive Primordial Black Holes from Hybrid Inflation as Dark Matter and the seeds of Galaxies,Phys. Rev. D92(2015) 023524 [1501.07565]

  26. [27]

    A. A. Starobinsky,Dynamics of Phase Transition in the New Inflationary Universe Scenario and Generation of Perturbations,Phys. Lett. B117(1982) 175

  27. [28]

    A. A. Starobinsky,STOCHASTIC DE SITTER (INFLATIONARY) STAGE IN THE EARLY UNIVERSE,Lect. Notes Phys.246(1986) 107

  28. [29]

    Nambu and M

    Y. Nambu and M. Sasaki,Stochastic Stage of an Inflationary Universe Model,Phys. Lett. B 205(1988) 441

  29. [30]

    Nambu and M

    Y. Nambu and M. Sasaki,Stochastic Approach to Chaotic Inflation and the Distribution of Universes,Phys. Lett. B219(1989) 240

  30. [31]

    H. E. Kandrup,STOCHASTIC INFLATION AS A TIME DEPENDENT RANDOM WALK, Phys. Rev. D39(1989) 2245

  31. [32]

    Nakao, Y

    K.-i. Nakao, Y. Nambu and M. Sasaki,Stochastic Dynamics of New Inflation,Prog. Theor. Phys.80(1988) 1041

  32. [33]

    Nambu,Stochastic Dynamics of an Inflationary Model and Initial Distribution of Universes, Prog

    Y. Nambu,Stochastic Dynamics of an Inflationary Model and Initial Distribution of Universes, Prog. Theor. Phys.81(1989) 1037

  33. [34]

    Mollerach, S

    S. Mollerach, S. Matarrese, A. Ortolan and F. Lucchin,Stochastic inflation in a simple two field model,Phys. Rev. D44(1991) 1670

  34. [35]

    A. D. Linde, D. A. Linde and A. Mezhlumian,From the Big Bang theory to the theory of a stationary universe,Phys. Rev. D49(1994) 1783 [gr-qc/9306035]

  35. [36]

    A. A. Starobinsky and J. Yokoyama,Equilibrium state of a selfinteracting scalar field in the De Sitter background,Phys. Rev. D50(1994) 6357 [astro-ph/9407016]

  36. [37]

    A. A. Starobinsky,Multicomponent de Sitter (Inflationary) Stages and the Generation of Perturbations,JETP Lett.42(1985) 152

  37. [38]

    D. S. Salopek and J. R. Bond,Nonlinear Evolution of Long Wavelength Metric Fluctuations in Inflationary mMdels,Phys. Rev. D42(1990) 3936

  38. [39]

    Sasaki and E

    M. Sasaki and E. D. Stewart,A General analytic formula for the spectral index of the density perturbations produced during inflation,Prog. Theor. Phys.95(1996) 71 [astro-ph/9507001]. – 18 –

  39. [40]

    Sasaki and T

    M. Sasaki and T. Tanaka,Superhorizon scale dynamics of multiscalar inflation,Prog. Theor. Phys.99(1998) 763 [gr-qc/9801017]

  40. [41]

    Wands, K

    D. Wands, K. A. Malik, D. H. Lyth and A. R. Liddle,A New approach to the evolution of cosmological perturbations on large scales,Phys. Rev. D62(2000) 043527 [astro-ph/0003278]

  41. [42]

    D. H. Lyth, K. A. Malik and M. Sasaki,A General proof of the conservation of the curvature perturbation,JCAP05(2005) 004 [astro-ph/0411220]

  42. [43]

    D. H. Lyth and Y. Rodriguez,The Inflationary prediction for primordial non-Gaussianity, Phys. Rev. Lett.95(2005) 121302 [astro-ph/0504045]

  43. [44]

    Fujita, M

    T. Fujita, M. Kawasaki, Y. Tada and T. Takesako,A new algorithm for calculating the curvature perturbations in stochastic inflation,JCAP12(2013) 036 [1308.4754]

  44. [45]

    Fujita, M

    T. Fujita, M. Kawasaki and Y. Tada,Non-perturbative approach for curvature perturbations in stochasticδNformalism,JCAP10(2014) 030 [1405.2187]

  45. [46]

    Vennin and A

    V. Vennin and A. A. Starobinsky,Correlation Functions in Stochastic Inflation,Eur. Phys. J. C75(2015) 413 [1506.04732]

  46. [47]

    Kawasaki and Y

    M. Kawasaki and Y. Tada,Can massive primordial black holes be produced in mild waterfall hybrid inflation?,JCAP08(2016) 041 [1512.03515]

  47. [48]

    Tada and M

    Y. Tada and M. Yamada,Primordial black hole formation in hybrid inflation,Phys. Rev. D 107(2023) 123539 [2304.01249]

  48. [49]

    Mizuguchi, T

    Y. Mizuguchi, T. Murata and Y. Tada,STOLAS: STOchastic LAttice Simulation of cosmic inflation,JCAP12(2024) 050 [2405.10692]

  49. [50]

    Pinol, S

    L. Pinol, S. Renaux-Petel and Y. Tada,A manifestly covariant theory of multifield stochastic inflation in phase space: solving the discretisation ambiguity in stochastic inflation,JCAP04 (2021) 048 [2008.07497]

  50. [51]

    I. F. Halpern, M. P. Hertzberg, M. A. Joss and E. I. Sfakianakis,A Density Spike on Astrophysical Scales from an N-Field Waterfall Transition,Phys. Lett. B748(2015) 132 [1410.1878]

  51. [52]

    Tada and M

    Y. Tada and M. Yamada,Stochastic dynamics of multi-waterfall hybrid inflation and formation of primordial black holes,JCAP11(2023) 089 [2306.07324]

  52. [53]

    Tada and M

    Y. Tada and M. Yamada,Multifield stochastic dynamics in GUT hybrid inflation without monopole problem and with gravitational wave signatures of GUT Higgs representation,Phys. Lett. B855(2024) 138854 [2405.08859]

  53. [54]

    Murata and Y

    T. Murata and Y. Tada,Stochastic tail of the curvature perturbation in hybrid inflation,Phys. Rev. D113(2026) 023552 [2507.22439]. [55]Planckcollaboration, N. Aghanim et al.,Planck 2018 results. VI. Cosmological parameters, Astron. Astrophys.641(2020) A6 [1807.06209]

  54. [56]

    Clesse, B

    S. Clesse, B. Garbrecht and Y. Zhu,Non-Gaussianities and Curvature Perturbations from Hybrid Inflation,Phys. Rev. D89(2014) 063519 [1304.7042]

  55. [57]

    Miyamoto and Y

    K. Miyamoto and Y. Tada,Calculating the power spectrum in stochastic inflation by Monte Carlo simulation and least squares curve fitting,2508.17654

  56. [58]

    R. L. Davis,COSMIC TEXTURE,Phys. Rev. D35(1987) 3705

  57. [59]

    R. L. Davis,Cosmic Texture and the Microwave Background,Phys. Rev. D36(1987) 997

  58. [60]

    Turok,Global Texture as the Origin of Cosmic Structure,Phys

    N. Turok,Global Texture as the Origin of Cosmic Structure,Phys. Rev. Lett.63(1989) 2625

  59. [61]

    K. N. Ananda, C. Clarkson and D. Wands,The Cosmological gravitational wave background from primordial density perturbations,Phys. Rev. D75(2007) 123518 [gr-qc/0612013]. – 19 –

  60. [62]

    Baumann, P

    D. Baumann, P. J. Steinhardt, K. Takahashi and K. Ichiki,Gravitational Wave Spectrum Induced by Primordial Scalar Perturbations,Phys. Rev. D76(2007) 084019 [hep-th/0703290]

  61. [63]

    Saito and J

    R. Saito and J. Yokoyama,Gravitational wave background as a probe of the primordial black hole abundance,Phys. Rev. Lett.102(2009) 161101 [0812.4339]

  62. [64]

    J. R. Espinosa, D. Racco and A. Riotto,A Cosmological Signature of the SM Higgs Instability: Gravitational Waves,JCAP09(2018) 012 [1804.07732]

  63. [65]

    Kohri and T

    K. Kohri and T. Terada,Semianalytic calculation of gravitational wave spectrum nonlinearly induced from primordial curvature perturbations,Phys. Rev. D97(2018) 123532 [1804.08577]

  64. [66]

    Dom` enech,Scalar Induced Gravitational Waves Review,Universe7(2021) 398 [2109.01398]

    G. Dom` enech,Scalar Induced Gravitational Waves Review,Universe7(2021) 398 [2109.01398]

  65. [67]

    Saikawa,A review of gravitational waves from cosmic domain walls,Universe3(2017) 40 [1703.02576]

    K. Saikawa,A review of gravitational waves from cosmic domain walls,Universe3(2017) 40 [1703.02576]. – 20 –