Pith. sign in

REVIEW 4 major objections 4 minor 1 cited by

Quantum correlations in a gravitational collapse simulation with SpheriCo.jl

T0 review · 4 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read SpheriCo.jl simulates stable semiclassical collapse and sees horizon-linking quantum correlation tongues.

desk verdict Solid code paper with an honest but not-yet-established physics hint: the across-horizon correlator tongues are shown only where the paper's own regulator condition is violated and fast modes have lost convergence. read the letter →

arxiv 2412.19722 v2 pith:JWWEZQCE submitted 2024-12-27 gr-qc hep-th

Benjamin Berczi , Magdalena Eriksson , Thanasis Giannakopoulos , Paul M. Saffin This is my paper · ORCID
classification gr-qchep-th PACS 04.25.dg04.62.+v04.70.Dy
keywords quantumfieldsincurvedspacetimegravitationalcollapseHawkingradiationnumericalrelativityPauli-Villarsregularizationcorrelationsapparenthorizonsummation-by-parts
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

SpheriCo.jl is a new open-access code that simulates the spherical gravitational collapse of a scalar field whose matter sector can be classical or a quantised field on a classical, dynamical background. The paper's central claim is that with summation-by-parts differentiation the semiclassical system stays stable long enough to compute field correlations after an apparent horizon forms. Running a supercritical pulse with Pauli-Villars-regularized quantum modes, the code sees two-point correlation tongues connecting points inside and outside the horizon, growing with time; subcritical runs show no persistent analogue. The authors present this as a hint of Hawking-pair correlations while explicitly noting that a more systematic parameter study is needed.

What carries the argument

The argument is carried by the equal-time two-point correlation function from equation (7.1), evaluated from the evolved quantum modes, together with the Pauli-Villars ghost-field sector that cancels ultraviolet divergences and the summation-by-parts finite-difference operators of [21] that keep the $1/r$ terms stable near the origin. The SBP operators are the practical enabler: they let the code evolve faster-oscillating modes and reach times where the horizon has formed, while the Pauli-Villars fields make the stress-energy expectation value finite on the lattice. The correlator is then sliced in time: tongues appear off the diagonal when one point sits inside the apparent-horizon radius and the other outside.

What would settle it

Re-run the supercritical correlator with the Pauli-Villars mass raised tenfold while keeping the same physical wavelength, and with more than 30 radial modes; if the horizon-linking tongues do not persist and converge, the reported correlation is a regularization or truncation artifact.

Watch

Extended reading notes

Core claim

The central discovery the authors report is that the equal-time two-point correlator $C(t; r_1, r_2)$ of the quantum scalar, computed in the semiclassical approximation on a dynamically collapsing spacetime, develops non-trivial 'tongues' away from the main diagonal that correlate a point inside the apparent horizon with one outside, and this happens only in supercritical evolutions. This is presented as evidence hinting at correlation between pairs of Hawking quanta. The correlation persists when backreaction is included, with the same qualitative shape and somewhat stronger oscillations around the horizon.

Load-bearing premise

The simulations use a Pauli-Villars ghost mass that is not much larger than the shortest wavelength resolved, so the correlation tongues could be a regularization artifact rather than a sign of Hawking radiation.

Editorial extensions

If this is right

  • Semiclassical collapse simulations can now run past apparent-horizon formation with stable evolution near the origin, allowing correlators to be evaluated after the horizon exists.
  • The classical module reproduces Choptuik critical-collapse behavior with second-order convergence, so the code can be used as a testbed for critical phenomena.
  • In supercritical runs the equal-time correlator develops off-diagonal tongues linking points inside and outside the apparent horizon, while subcritical runs relax back to a Minkowski-like profile.
  • When backreaction is included, the apparent-horizon location, area, and mass track the classical run within ten percent, and the correlation picture persists with stronger oscillations around the horizon.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the tongues are physical, equal-time correlation across the horizon could become a practical numerical witness of Hawking-pair formation, though the paper stops short of claiming entanglement.
  • The tongues' angle may depend on the gauge slicing; computing the correlator on the double-null surfaces the paper already provides could give a slicing-invariant diagnostic.
  • Because the SBP treatment is generic, the same pipeline could be applied to other spherically symmetric semiclassical settings, such as regular black holes or dynamical horizons, without new regularization work.
  • The reported failure of Z4 constraint damping for this variable set suggests that damping schemes should be re-tuned or extended to reduction constraints before being used in long runs.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 4 minor

Summary. This paper presents SpheriCo.jl, an open-access Julia code for spherically symmetric gravitational collapse of a scalar field in classical and semiclassical (quantum-field-on-classical-background) settings, with the quantum stress tensor regularized by five Pauli-Villars ghost fields and WKB-derived counterterms for the cosmological constant and Planck mass. The numerical method combines summation-by-parts operators for the 1/r terms at the origin, a Z4-type constraint-damping option, Kreiss-Oliger dissipation, and an infalling-outer-boundary option, and it is validated extensively: second-order convergence for smooth and noisy data in super- and subcritical runs, reproduction of Choptuik-type echoing and universality, convergence to analytic Minkowski mode solutions, and a Minkowski-consistency test for backreaction. The central physics result is the equal-time two-point correlator (7.1) in collapse runs: supercritical runs develop off-diagonal 'tongues' correlating points inside and outside the apparent horizon (Fig. 16), which the authors present as a hint of Hawking-pair correlations, while subcritical runs show no persistent analogue. Sections 7.2 and 8 report that the correlator and horizon quantities are qualitatively similar with backreaction, within parameter ranges the authors themselves flag as limited.

Significance. If the correlator tongues are physical, this would be a notable advance: a four-dimensional, spherically symmetric numerical computation of across-horizon quantum correlations in a dynamically collapsing spacetime, going beyond earlier two-dimensional semiclassical collapse studies and complementing analogue-horizon observations. The manuscript has real strengths independent of the headline claim: the code ships with a public repository and archived data; the validation is machine-checkable and broad (analytic Minkowski modes, Choptuik universality, robust-stability tests); the WKB derivation of the Pauli-Villars counterterms (Appendix C) is a genuine methodological contribution; and the authors are unusually candid about their diagnostics, including the failure of their constraint-damping implementation (Section 5.3). The skeptic's concern lands: the paper's own Section 7.1 admits that the runs producing the tongues violate the PV-decoupling condition MPV ≫ 2π/λ̃, and Section 6.1 (Fig. 12) shows the fastest ghost modes losing second-order convergence within the displayed time range.

major comments (4)
  1. [§7.1 (final paragraph), §4 (third paragraph), Eq. (7.1)] The central result — the across-horizon correlation tongues of Fig. 16 — is computed with MPV = 1, a value that violates the paper's own decoupling condition MPV ≫ 2π/λ̃ ≈ 3.6 stated in Section 4 and explicitly acknowledged in Section 7.1 ('we were not able to respect the condition ... Here we have 2π/λ̃ ≃ 3.6'). Because the correlator (7.1) is an alternating sum over the physical field and the five ghost fields with masses (3.12) of order MPV, the ghost sectors are not suppressed at the wavelengths of the physical pulse. No computation of C(t; r1, r2) for larger MPV — even in the no-backreaction case where the Section 6.2 backreaction instability does not apply — and no estimate of the ghost-sector fraction of the correlator is provided. The abstract's claim that the results 'hint at a non-trivial correlation across the horizon of Hawking quanta' therefore rests on an untested assumption about regulator robustness rather than on the paper's usual standard of evidence.
  2. [§6.1, Fig. 12; §7.1, Fig. 16] The fastest modes of the heaviest ghost field (l = 60, n = 5, µ5 = 2MPV) lose second-order convergence for t ≳ 7.9 (Fig. 12), which is inside the time range where the tongues are displayed in Fig. 16 (last panel, t = 7.96875), and the correlator runs use lmax = 90, beyond the l = 60 at which convergence was tested. No resolution study or kmax/lmax-study of the correlator itself is presented: the subcritical control in Fig. 16 only shows that the tonguelike feature is associated with horizon formation, not that the supercritical feature is converged. Since (7.1) sums over all modes up to the truncation, oscillatory errors in high-l, high-mass ghost modes at t ≳ 7.9 could contribute to the observed off-diagonal structure, so the late-time growth of the tongues is not currently established as a physical effect.
  3. [§7.1 (gauge paragraph), §8] The foliation dependence of the claimed signal is asserted but not tested. The text states that 'we expect the details of the correlation tongues to depend on the choice of lapse function, but the qualitative structure of correlations outside and within the horizon to remain' and immediately afterwards that 'It is possible that our picture is greatly affected by the collapse of the lapse inside the apparent horizon', while Section 8 defers a wider parameter scan to future work. Because the slicing determines which points are related as 'inside' and 'outside' at equal time, the physical interpretation of the tongues requires at least one independent gauge computation (e.g., harmonic or polar slicing), which is not provided. In addition, the abstract's phrase 'correlation across the horizon of Hawking quanta' is stronger than the body's hedged statement (Section 7, opening) that the correlation 'may be interpreted as correlating pairs of Hawking quanta'; if a gauge test is not feasible, the abstract and Section 8 should be re-scoped accordingly.
  4. [§3.2, Eqs. (3.13)–(3.15) vs Appendix C, Eqs. (C.11)–(C.15)] The counterterm expressions in the main text and in the WKB derivation do not agree: Eqs. (3.13)–(3.14) quote log arguments 39/1216 and 24/33 with coefficient 1/(12(2π)²) on the Gab term, whereas Eqs. (C.11)–(C.15) give 39/216 and 28/36 with coefficient 1/(24(2π)²). The effective Planck mass (3.15) constructed from these counterterms is used in the backreaction runs of Section 7.2 (Figs. 17–18), so the inconsistency leaves the reader unable to determine which values were implemented. Please reconcile the two presentations and state explicitly which counterterms were used in production.
minor comments (4)
  1. [§2.2, Eq. (2.17)] The domain is written as '0 ≥ t ≥ tmax, 0 ≥ r ≥ rmax'; the inequalities are reversed and should read 0 ≤ t ≤ tmax and 0 ≤ r ≤ rmax.
  2. [§7.1, Eq. (7.1)] In the definition of C(t1, t2; r1, r2), the second field factor is written as 'Φ̂(t, r2)' but should be 'Φ̂(t2, r2)'; as written the expression mixes t and t2 in the same formula.
  3. [§6.2, Fig. 13] The three MPV panels use very different vertical scales (MPV = 0.1: up to 2e-5; MPV = 1: up to 1; MPV = 2: up to 6), which makes the several-orders-of-magnitude difference between the cases difficult to read; a common logarithmic scale would be clearer, and the caption should state that the MPV = 2 runs crash during the evolution.
  4. [§5.3, §7.1] Since the constraint-damping option is reported in Section 5.3 to be unstable or to increase constraint violation for many parameter choices, please state explicitly whether damping was enabled (and with which κ1, κ2) in the Section 7 correlator runs, so that the reader can assess the quality of the background geometry used in the mode evolution.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the correlator is computed directly from evolved mode functions; self-cited Pauli-Villars scheme is a method, not a fitted target.

full rationale

The central claim is the equal-time two-point correlator (7.1), which is defined directly as a mode sum over the evolved physical and Pauli-Villars mode functions and is evaluated without fitting any parameter to the target signal. The mode functions are evolved from Minkowski initial data (3.21), and the quantum sector is validated against the analytic Minkowski solutions (6.1) with second-order convergence. The counterterms used in the backreaction sector are derived in appendix C via a WKB calculation on an FRW background, not imported from the claimed result. Self-citations [17,18] supply the Pauli-Villars regularization method and the empirical lmax/kmax ratio, but these are methodological inputs rather than a re-statement of the correlator result; the correlation tongues are not chosen or fitted by that scheme. The paper explicitly flags robustness limitations in section 7.1, admitting that MPV = 1 fails the condition MPV ≫ 2π/λ ≈ 3.6, and section 6.1 shows the fastest modes (l = 60, n = 5) lose second-order convergence after t ≈ 7.9, which overlaps with the displayed correlator times. These are correctness and convergence concerns, not circularity: they weaken confidence in the physical interpretation but do not make the prediction equivalent to its input by construction. No self-definitional step, fitted input renamed as prediction, or load-bearing self-citation chain is present.

Assumptions & free parameters 5 free parameters · 7 assumptions · 1 invented entities

The central claim rests on the numerical scheme's ability to approximate the renormalized quantum stress tensor with a truncated mode sum and a non-decoupled Pauli-Villars regulator, and on the gauge choice for the correlator. Free parameters are the regulator mass, mode truncations, momentum spacing, filter and dissipation settings; the five PV ghost fields are invented entities with no independent observational handle. The principal axioms are the validity of the semiclassical approximation, the sufficiency of the five-ghost spectrum, the WKB counterterm computation, the Minkowski definition of the vacuum, the extrapolated SBP operators, the innocuousness of boundary and filter treatments, and the assumed gauge-independence of the correlation structure.

free parameters (5)
  • Pauli-Villars mass MPV = 1
    Chosen empirically: larger values make the vacuum stress tensor converge too slowly for the affordable number of modes (figure 14), and the paper uses MPV = 1 for the correlator runs even though it violates the condition MPV >> 2*pi/lambda.
  • Mode truncation (kmax, lmax) = 30 and 90, with lmax/kmax = 3
    Truncation of the mode expansion; the ratio 3 was found empirically to make the vacuum stress tensor vanish over a sufficiently large r-domain (section 4.3). No convergence study of the correlator with respect to these truncations is provided.
  • Momentum spacing dk = pi/30
    Momentum spacing for the discretized mode sum.
  • Filter radius rcut and steepness p = rcut = 20, p = 1
    Smooth step function that suppresses the quantum stress tensor near rmax; chosen to guarantee asymptotic flatness in the simulation.
  • Kreiss-Oliger dissipation coefficient sigma = 0.02
    Artificial dissipation coefficient chosen empirically for stability.
assumptions (7)
  • domain assumption The semiclassical approximation (classical geometry sourced by the expectation value of the quantum stress tensor) is valid throughout the collapse.
    Invoked in section 3; assumes curvature stays well below the Planck scale and the mass is much larger than M_Pl.
  • domain assumption Five Pauli-Villars ghost fields with masses MPV, MPV, sqrt(3) MPV, sqrt(3) MPV, 2 MPV fully regularize the stress-energy tensor.
    Inherited from [18] and used in section 3.2; the counterterm calculation in appendix C assumes this spectrum.
  • standard math The truncated WKB expansion (zeroth plus second adiabatic order) is sufficient to extract the Lambda and delta M_Pl^2 counterterms.
    Used in appendix C to derive equation (3.13); higher-order terms are assumed negligible at the considered order.
  • domain assumption The vacuum state is defined by the Minkowski mode functions of the asymptotic past.
    Stated in section 3; the initial data (3.21) use these modes. Any ambiguity of the vacuum in a dynamical spacetime affects the correlator.
  • domain assumption The summation-by-parts operators of [21], designed for the flat-space wave equation, remain second-order accurate and stable when extrapolated to the curved, non-linear evolution system.
    Section 4.2 states 'we heavily extrapolate their use past their original purpose'; the classical and quantum systems couple them to metric functions.
  • domain assumption The radiative boundary conditions at rmax and the filter function F(r) do not contaminate the causal domain r < rcausal.
    Sections 2.2 and 4.3; results are presented only in the causally disconnected region, but the filter also multiplies the stress tensor and cosmological constant on the whole domain.
  • ad hoc to paper The qualitative structure of the across-horizon correlations is independent of the 1+log gauge choice.
    Section 7.1 states that 'we expect the details of the correlation tongues to depend on the choice of lapse function, but the qualitative structure to remain' and suggests testing harmonic gauge in future work; this is an untested assumption on which the physical interpretation rests.
invented entities (1)
  • Five Pauli-Villars ghost fields
    purpose: Cancel the ultraviolet divergences of the physical scalar field's stress-energy tensor
    Their masses (3.12) are chosen for regularization, not tied to any observed particle; with MPV = 1 they are not decoupled from the dynamics, so they could influence the computed correlator.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Quantum correlations in a gravitational collapse simulation with SpheriCo.jl." pith.science (2026). https://pith.science/paper/JWWEZQCE

@misc{pith2026241219722,
  author       = {Pith},
  title        = {Pith review of: Quantum correlations in a gravitational collapse simulation with SpheriCo.jl},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JWWEZQCE}},
  note         = {Machine review of arXiv:2412.19722}
}
read the original abstract

We report on work using a newly developed code, SpheriCo.jl, that computes the gravitational collapse of a spherical scalar field, where the scalar can be either a classical field, or a quantum field operator. By utilising summation-by-parts methods for the numerical derivatives we are able to simulate the collapse longer than was possible previously due to enhanced numerical stability. We present a suite of tests for the code that tests its accuracy and stability, both for the classical and quantum fields. We are able to observe critical behavior of gravitational collapse for the classical setup, in agreement with expected results. The code is also used to compute two-point correlation functions, with results that hint at a non-trivial correlation across the horizon of Hawking quanta.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. Quantum fields in boson star spacetime

    gr-qc 2026-01 conditional novelty 6.0 of 10

    In boson star spacetimes, the renormalized quantum stress tensor has mostly positive energy density and negative radial pressure that grow with curvature, rivaling the classical stress in the most compact solutions.

Reference graph

Works this paper leans on

66 extracted references · 34 canonical work pages · cited by 1 Pith paper

  1. [1]

    Christodoulou, The Problem of a Selfgravitating Scalar Field , Commun

    D. Christodoulou, The Problem of a Selfgravitating Scalar Field , Commun. Math. Phys. 105 (1986) 337

    work page 1986

  2. [2]

    Christodoulou, Global Existence of Generalized Solutions of the Spherically Symmetric Einstein Scalar Equations in the Large , Commun

    D. Christodoulou, Global Existence of Generalized Solutions of the Spherically Symmetric Einstein Scalar Equations in the Large , Commun. Math. Phys. 106 (1986) 587

    work page 1986

  3. [3]

    Christodoulou, A Mathematical Theory of Gravitational Collapse , Commun

    D. Christodoulou, A Mathematical Theory of Gravitational Collapse , Commun. Math. Phys. 109 (1987) 613

    work page 1987

  4. [4]

    Christodoulou, The Structure and Uniqueness of Generalized Solutions of the Spherically Symmetric Einstein Scalar Equations , Commun

    D. Christodoulou, The Structure and Uniqueness of Generalized Solutions of the Spherically Symmetric Einstein Scalar Equations , Commun. Math. Phys. 109 (1987) 591

    work page 1987

  5. [5]

    Alcubierre, Introduction to 3+1 Numerical Relativity , International Series of Monographs on Physics, OUP Oxford (2008)

    M. Alcubierre, Introduction to 3+1 Numerical Relativity , International Series of Monographs on Physics, OUP Oxford (2008)

    work page 2008

  6. [6]

    Choptuik, Universality and scaling in gravitational collapse of a massless scalar field , Phys

    M.W. Choptuik, Universality and scaling in gravitational collapse of a massless scalar field , Phys. Rev. Lett. 70 (1993) 9

    work page 1993

  7. [7]

    Choptuik scaling in null coordinates

    D. Garfinkle, Choptuik scaling in null coordinates , Phys. Rev. D 51 (1995) 5558 [gr-qc/9412008]

    work page Pith review arXiv 1995

  8. [8]

    The spherically symmetric collapse of a massless scalar field

    R.S. Hamade and J.M. Stewart, The Spherically symmetric collapse of a massless scalar field, Class. Quant. Grav. 13 (1996) 497 [ gr-qc/9506044]

    work page Pith review arXiv 1996

Show all 66 references
  1. [9]

    Russo, L

    J.G. Russo, L. Susskind and L. Thorlacius, The Endpoint of Hawking radiation , Phys. Rev. D 46 (1992) 3444 [ hep-th/9206070]

    1992 arXiv

  2. [10]

    Strominger and L

    A. Strominger and L. Thorlacius, Universality and scaling at the onset of quantum black hole formation, Phys. Rev. Lett. 72 (1994) 1584 [ hep-th/9312017]

    1994 arXiv

  3. [11]

    Piran and A

    T. Piran and A. Strominger, Numerical analysis of black hole evaporation , Phys. Rev. D 48 (1993) 4729 [ hep-th/9304148]

    1993 arXiv

  4. [12]

    Lowe, Semiclassical approach to black hole evaporation , Phys

    D.A. Lowe, Semiclassical approach to black hole evaporation , Phys. Rev. D 47 (1993) 2446 [hep-th/9209008]

    1993 arXiv

  5. [13]

    Levi and A

    A. Levi and A. Ori, Pragmatic mode-sum regularization method for semiclassical black-hole spacetimes, Phys. Rev. D 91 (2015) 104028 [ 1503.02810]. – 41 –

    2015 arXiv

  6. [14]

    Levi and A

    A. Levi and A. Ori, Mode-sum regularization of ϕ2 in the angular-splitting method , Phys. Rev. D 94 (2016) 044054 [ 1606.08451]

    2016 arXiv

  7. [15]

    Levi and A

    A. Levi and A. Ori, Versatile method for renormalized stress-energy computation in black-hole spacetimes, Phys. Rev. Lett. 117 (2016) 231101 [ 1608.03806]

    2016 arXiv

  8. [16]

    A. Levi, E. Eilon, A. Ori and M. van de Meent, Renormalized stress-energy tensor of an evaporating spinning black hole , Phys. Rev. Lett. 118 (2017) 141102 [ 1610.04848]

    2017 arXiv

  9. [17]

    Berczi, P.M

    B. Berczi, P.M. Saffin and S.-Y. Zhou, Gravitational collapse with quantum fields , Phys. Rev. D 104 (2021) L041703 [ 2010.10142]

    2021 arXiv

  10. [18]

    Berczi, P.M

    B. Berczi, P.M. Saffin and S.-Y. Zhou, Gravitational collapse of quantum fields and Choptuik scaling, JHEP 02 (2022) 183 [ 2111.11400]

    2022 arXiv

  11. [19]

    Guenther, C

    J.N. Guenther, C. Hoelbling and L. Varnhorst, Semiclassical gravitational collapse of a radially symmetric massless scalar quantum field , Phys. Rev. D 105 (2022) 105010 [2010.13215]

    2022 arXiv

  12. [20]

    Hoelbling, J.N

    C. Hoelbling, J.N. Guenther and L. Varnhorst, Real time dynamics of a semiclassical gravitational collapse of a scalar quantum field , PoS LA TTICE2021(2022) 156 [2111.15562]

    2022 arXiv

  13. [21]

    Gundlach, J.M

    C. Gundlach, J.M. Martin-Garcia and D. Garfinkle, Summation by parts methods for spherical harmonic decompositions of the wave equation in any dimensions , Class. Quant. Grav. 30 (2013) 145003 [ 1010.2427]

    2013 arXiv

  14. [22]

    Hawking, Particle Creation by Black Holes , Commun

    S.W. Hawking, Particle Creation by Black Holes , Commun. Math. Phys. 43 (1975) 199

    1975

  15. [23]

    Schutzhold and W.G

    R. Schutzhold and W.G. Unruh, On Quantum Correlations across the Black Hole Horizon , Phys. Rev. D 81 (2010) 124033 [ 1002.1844]

    2010 arXiv

  16. [24]

    Fontana and M

    M. Fontana and M. Rinaldi, Stress-energy tensor correlations across regular black holes horizons, Phys. Rev. D 108 (2023) 125003 [ 2302.08804]

    2023 arXiv

  17. [25]

    Carusotto, S

    I. Carusotto, S. Fagnocchi, A. Recati, R. Balbinot and A. Fabbri, Numerical observation of Hawking radiation from acoustic black holes in atomic BECs , New J. Phys. 10 (2008) 103001 [0803.0507]

    2008 arXiv

  18. [26]

    Balbinot and A

    R. Balbinot and A. Fabbri, Quantum correlations across the horizon in acoustic and gravitational black holes , Phys. Rev. D 105 (2022) 045010 [ 2107.00702]

    2022 arXiv

  19. [27]

    Spherico.jl

    “Spherico.jl.” https://github.com/ThanasisGiannakopoulos/SpheriCo.jl

  20. [28]

    Bezanson, A

    J. Bezanson, A. Edelman, S. Karpinski and V.B. Shah, Julia: A fresh approach to numerical computing, SIAM Review 59 (2017) 65 [ https://doi.org/10.1137/141000671]

    2017 doi

  21. [29]

    Alcubierre and J.A

    M. Alcubierre and J.A. Gonzalez, Regularization of spherically symmetric evolution codes in numerical relativity, Comput. Phys. Commun. 167 (2005) 76 [ gr-qc/0401113]

    2005 arXiv

  22. [30]

    Pretorius, Evolution of binary black hole spacetimes , Phys

    F. Pretorius, Evolution of binary black hole spacetimes , Phys. Rev. Lett. 95 (2005) 121101 [gr-qc/0507014]

    2005 arXiv

  23. [31]

    Gundlach, J.M

    C. Gundlach, J.M. Martin-Garcia, G. Calabrese and I. Hinder, Constraint damping in the Z4 formulation and harmonic gauge , Class. Quant. Grav. 22 (2005) 3767 [ gr-qc/0504114]

    2005 arXiv

  24. [32]

    Bernuzzi and D

    S. Bernuzzi and D. Hilditch, Constraint violation in free evolution schemes: Comparing BSSNOK with a conformal decomposition of Z4 , Phys. Rev. D 81 (2010) 084003 [0912.2920]

    2010 arXiv

  25. [33]

    C. Bona, T. Ledvinka, C. Palenzuela and M. Zacek, General covariant evolution formalism – 42 – for numerical relativity , Phys. Rev. D 67 (2003) 104005 [ gr-qc/0302083]

    2003 arXiv

  26. [34]

    Bona and C

    C. Bona and C. Palenzuela, Dynamical shift conditions for the Z4 and BSSN hyperbolic formalisms, Phys. Rev. D 69 (2004) 104003 [ gr-qc/0401019]

    2004 arXiv

  27. [35]

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

    1995 arXiv

  28. [36]

    Anninos, K

    P. Anninos, K. Camarda, J. Masso, E. Seidel, W.-M. Suen and J. Towns, Three-dimensional numerical relativity: The Evolution of black holes , Phys. Rev. D 52 (1995) 2059 [gr-qc/9503025]

    1995 arXiv

  29. [37]

    Cardoso, D

    V. Cardoso, D. Hilditch, K. Marouda, J. Nat´ ario and U. Sperhake, Curvature and dynamical spacetimes: can we peer into the quantum regime? , Class. Quant. Grav. 40 (2023) 065008 [2209.12589]

    2023 arXiv

  30. [38]

    Pareschi and G

    L. Pareschi and G. Russo, Implicit-explicit runge-kutta schemes for stiff systems of differential equations, in Recent Trends in Numerical Analysis, (USA), p. 269–288, Nova Science Publishers, Inc. (2000)

    2000

  31. [39]

    Evans, An Approach for Calculating Axisymmetric Gravitational Collapse , phd thesis, University of Texas at Austin, 1984

    C.R. Evans, An Approach for Calculating Axisymmetric Gravitational Collapse , phd thesis, University of Texas at Austin, 1984. Available at https://users.physics.unc.edu/~evans/cr_evans_thesis_1984.pdf

    1984

  32. [40]

    Su´ arez Fern´ andez, R

    I. Su´ arez Fern´ andez, R. Vicente and D. Hilditch,Semilinear wave model for critical collapse , Phys. Rev. D 103 (2021) 044016 [ 2007.13764]

    2021 arXiv

  33. [41]

    Gustafsson, H.-O

    B. Gustafsson, H.-O. Kreiss and J. Oliger, Time dependent problems and difference methods , Wiley, New York (1995)

    1995

  34. [42]

    Kreiss and J

    H. Kreiss and J. Oliger, Methods for the approximate solution of time dependent problems , GARP Publications Series, SELBSTVERL.) FEBR (1973)

    1973

  35. [43]

    Interpolations.jl

    “Interpolations.jl.” https://github.com/JuliaMath/Interpolations.jl

  36. [44]

    Rinne, Type II critical collapse on a single fixed grid: a gauge-driven ingoing boundary method, Gen

    O. Rinne, Type II critical collapse on a single fixed grid: a gauge-driven ingoing boundary method, Gen. Rel. Grav. 52 (2020) 117 [ 2008.12726]

    2020 arXiv

  37. [45]

    Lindblom, M.A

    L. Lindblom, M.A. Scheel, L.E. Kidder, R. Owen and O. Rinne, A New generalized harmonic evolution system , Class. Quant. Grav. 23 (2006) S447 [ gr-qc/0512093]

    2006 arXiv

  38. [46]

    D. Cors, S. Renkhoff, H.R. R¨ uter, D. Hilditch and B. Br¨ ugmann,Formulation improvements for critical collapse simulations , Phys. Rev. D 108 (2023) 124021 [ 2308.01812]

    2023 arXiv

  39. [47]

    Alcubierre et al., Toward standard testbeds for numerical relativity , Class

    M. Alcubierre et al., Toward standard testbeds for numerical relativity , Class. Quant. Grav. 21 (2004) 589 [ gr-qc/0305023]

    2004 arXiv

  40. [48]

    Calabrese, I

    G. Calabrese, I. Hinder and S. Husa, Numerical stability for finite difference approximations of Einstein ’s equations, J. Comput. Phys. 218 (2006) 607 [ gr-qc/0503056]

    2006 arXiv

  41. [49]

    Witek, D

    H. Witek, D. Hilditch and U. Sperhake, Stability of the puncture method with a generalized baumgarte-shapiro-shibata-nakamura formulation, Phys. Rev. D 83 (2011) 104041

    2011

  42. [50]

    Cao and D

    Z. Cao and D. Hilditch, Numerical stability of the Z4c formulation of general relativity , Phys. Rev. D 85 (2012) 124032 [ 1111.2177]

    2012 arXiv

  43. [51]

    Babiuc et al., Implementation of standard testbeds for numerical relativity , Class

    M.C. Babiuc et al., Implementation of standard testbeds for numerical relativity , Class. Quant. Grav. 25 (2008) 125012 [ 0709.3559]

    2008 arXiv

  44. [52]

    Giannakopoulos, Characteristic formulations of general relativity and applications , Ph.D

    T. Giannakopoulos, Characteristic formulations of general relativity and applications , Ph.D. – 43 – thesis, Instituto Superior T´ ecnico, Lisbon, 8, 2023.2308.16001

    2023

  45. [53]

    Hinder, Well-posed formulations and stable finite differencing schemes for numerical relativity, Ph.D

    I. Hinder, Well-posed formulations and stable finite differencing schemes for numerical relativity, Ph.D. thesis, U. Southampton (main), 2005. 10.5258/SOTON/393564

    2005 doi

  46. [54]

    Gundlach and J.M

    C. Gundlach and J.M. Martin-Garcia, Critical phenomena in gravitational collapse , Living Rev. Rel. 10 (2007) 5 [ 0711.4620]

    2007 arXiv

  47. [55]

    Baumgarte, Aspherical deformations of the Choptuik spacetime , Phys

    T.W. Baumgarte, Aspherical deformations of the Choptuik spacetime , Phys. Rev. D 98 (2018) 084012 [ 1807.10342]

    2018 arXiv

  48. [56]

    Marouda, D

    K. Marouda, D. Cors, H.R. R¨ uter, F. Atteneder and D. Hilditch, Twist-free axisymmetric critical collapse of a complex scalar field , Phys. Rev. D 109 (2024) 124042 [ 2402.06724]

    2024 arXiv

  49. [57]

    Jimenez-Vazquez and M

    E. Jimenez-Vazquez and M. Alcubierre, Critical gravitational collapse of a massive complex scalar field, Phys. Rev. D 106 (2022) 044071 [ 2206.01389]

    2022 arXiv

  50. [58]

    Simulations from SpheriCo.jl

    “Simulations from SpheriCo.jl.” https://youtube.com/playlist?list=PLyJ_ PP8qOGF8pj5YvAKoQHOc9wp6fSUNk&si=ZDa-8cMWwsSRHhBV

  51. [59]

    Steinhauer, Observation of quantum Hawking radiation and its entanglement in an analogue black hole , Nature Phys

    J. Steinhauer, Observation of quantum Hawking radiation and its entanglement in an analogue black hole , Nature Phys. 12 (2016) 959 [ 1510.00621]

    2016 arXiv

  52. [60]

    Jacquet, M

    M.J. Jacquet, M. Joly, L. Giacomelli, F. Claude, Q. Glorieux, A. Bramati et al., Analogue quantum simulation of the Hawking effect in a polariton superfluid , Eur. Phys. J. D 76 (2022) 152 [ 2201.02038]

    2022 arXiv

  53. [61]

    Berczi, M

    B. Berczi, M. Eriksson, T. Giannakopoulos and P. Saffin, Quantum correlations in a gravitational collapse simulation with SpheriCo.jl: data , 12, 2024. https://doi.org/10.5281/zenodo.14365566

    2024 doi

  54. [62]

    Pla and E

    S. Pla and E. Winstanley, Equivalence of the adiabatic expansion and Hadamard renormalization for a charged scalar field , Phys. Rev. D 107 (2023) 025004 [ 2209.01079]

    2023 arXiv

  55. [63]

    Baacke, K

    J. Baacke, K. Heitmann and C. Patzold, Renormalization of nonequilibrium dynamics in FR W cosmology, Phys. Rev. D 56 (1997) 6556 [ hep-ph/9706274]

    1997 arXiv

  56. [64]

    Baacke and C

    J. Baacke and C. Patzold, Out-of-equilibrium evolution of scalar fields in FR W cosmology: Renormalization and numerical simulations , Phys. Rev. D 61 (2000) 024016 [hep-ph/9906417]

    2000 arXiv

  57. [65]

    Markkanen and A

    T. Markkanen and A. Tranberg, A Simple Method for One-Loop Renormalization in Curved Space-Time, JCAP 08 (2013) 045 [ 1303.0180]

    2013 arXiv

  58. [66]

    Paz and F.D

    J.P. Paz and F.D. Mazzitelli, Renormalized Evolution Equations for the Back Reaction Problem With a Selfinteracting Scalar Field , Phys. Rev. D 37 (1988) 2170. A Hyperbolicity of the classical system The hyperbolicity of a partial differential equation (PDE) system is connecte...

    1988

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.