pith. machine review for the scientific record. sign in

arxiv: 2604.25635 · v2 · submitted 2026-04-28 · 🧮 math.DS · gr-qc· physics.comp-ph

Recognition: 2 theorem links

· Lean Theorem

Numerical Investigations of Stable Dynamics in the Presence of Ghosts

Jax Wysong , Samara Overvaag , Hyun Lim , Jung-Han Kimn
Authors on Pith no claims yet

Pith reviewed 2026-05-13 07:58 UTC · model grok-4.3

classification 🧮 math.DS gr-qcphysics.comp-ph
keywords ghost degrees of freedomnonlinear field dynamicsscalar fieldsnumerical simulationsmetastable statesoscillonsstability
0
0 comments X

The pith

Ghost-normal scalar field systems can maintain dynamically bounded evolution for extended times when initial data favors ultraviolet modes and small amplitudes.

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

The paper performs numerical simulations of two coupled scalar fields with opposite kinetic signs in flat spacetime, one normal and one ghost. It finds that these systems avoid immediate runaway and instead show long-lived bounded behavior whose duration depends on the initial frequency content and size. Ultraviolet-heavy small-amplitude data stays stable much longer than infrared or large-amplitude data because instability requires nonlinear energy exchange between the sectors. Specific potentials, such as a lifted sixth-order interaction, further create temporary metastable regimes by supporting oscillon-like objects that slow the growth.

Core claim

Ghost-normal systems can exhibit long-lived, dynamically bounded evolution over extended time intervals, with stability strongly controlled by spectral content and amplitude. Ultraviolet-dominated and small-amplitude configurations remain stable significantly longer than infrared-dominated or large-amplitude data, indicating that instability is mediated by nonlinear spectral energy transfer rather than instantaneous runaway. Nonlinear self-interactions play a dual role: while they can accelerate energy exchange between sectors, certain potentials generate transient metastable regimes that partially suppress ghost-induced growth.

What carries the argument

Spacetime finite-element discretization of two coupled scalar fields carrying opposite-sign kinetic terms, used to evolve broad classes of initial data in one and two spatial dimensions.

If this is right

  • Instability develops through gradual nonlinear transfer of energy between the normal and ghost sectors rather than linear exponential growth.
  • Initial data dominated by high-frequency modes and small overall amplitude postpones the onset of instability.
  • Potentials that admit oscillon-like structures can temporarily suppress ghost-induced growth and create metastable intervals.
  • The overall dynamical outcome depends on the balance of dispersion, nonlinearity strength, and phase structure rather than the mere presence of a ghost.

Where Pith is reading between the lines

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

  • If the discretization is faithful, classical field theories containing ghosts need not collapse on short timescales and may admit long-lived metastable configurations.
  • The same numerical strategy could be used to test whether analogous metastability appears in higher dimensions or with gravitational coupling.
  • Observing persistent oscillon-like objects in ghost models would connect these results to soliton stability questions in other nonlinear theories.

Load-bearing premise

The finite-element grid reproduces the continuous equations without adding artificial damping or dispersion that would hide genuine ghost-driven growth.

What would settle it

A sufficiently long simulation starting from ultraviolet-dominated small-amplitude data that eventually shows unbounded growth or loss of boundedness would falsify the reported stability hierarchy.

Figures

Figures reproduced from arXiv: 2604.25635 by Hyun Lim, Jax Wysong, Jung-Han Kimn, Samara Overvaag.

Figure 1
Figure 1. Figure 1: FIG. 1: Representative evolution of the coupled view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Evolution with Gaussian packet initial data. The view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Here we illustrate the evolution of the system view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Phase-correlated plane waves initialization case. view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Evolution of the view at source ↗
Figure 6
Figure 6. Figure 6: figure 6 view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Evolution of the view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Evolution of the view at source ↗
Figure 10
Figure 10. Figure 10: shows a non-monotonic dependence of life￾time on amplitude. Three regimes can be identified. In the small amplitude (A << Ac) case, dynamics remain perturbative. Nonlinear self-interactions are weak, and lifetime decreases gradually with increasing amplitude due to enhanced quartic coupling. At near critical am￾plitude (A ≃ Ac), the effective potential flattens near its minimum as the amplitude approaches… view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: Comparison of spacetime surface plots for the view at source ↗
read the original abstract

We explore the nonlinear dynamics of classical field theories containing ghost degrees of freedom, focusing on two coupled scalar fields with opposite kinetic terms in (1+1) and (2+1) dimensional Minkowski spacetime. Using a spacetime finite element formulation, we perform a systematic numerical study across a broad class of initial data. We find that ghost-normal systems can exhibit long-lived, dynamically bounded evolution over extended time intervals, with stability strongly controlled by spectral content and amplitude. Ultraviolet-dominated and small-amplitude configurations remain stable significantly longer than infrared-dominated or large-amplitude data, indicating that instability is mediated by nonlinear spectral energy transfer rather than instantaneous runaway. Nonlinear self-interactions play a dual role: while they can accelerate energy exchange between sectors, certain potentials, including a lifted $\phi^6$ interaction supporting oscillon-like structures, generate transient metastable regimes that partially suppress ghost-induced growth. Our results demonstrate that the dynamical consequences of ghost modes in classical field theory depend sensitively on dispersion, nonlinearity, and phase structure, revealing a richer metastability landscape than commonly assumed.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript reports a numerical study of two coupled scalar fields with opposite-sign kinetic terms in (1+1) and (2+1)D Minkowski spacetime, using a spacetime finite element formulation. It claims that ghost-normal systems can exhibit long-lived bounded evolution over extended times, with stability strongly controlled by spectral content and amplitude: UV-dominated and small-amplitude initial data remain stable significantly longer than IR-dominated or large-amplitude data. Instability is attributed to nonlinear spectral energy transfer rather than instantaneous runaway, and certain potentials (including a lifted φ⁶ interaction) are shown to generate transient metastable regimes that partially suppress ghost-induced growth.

Significance. If the reported metastability accurately reflects the continuous dynamics, the work would demonstrate that ghost instabilities in classical field theories are not invariably immediate but can be delayed or modulated by dispersion relations, nonlinearity, and phase structure, revealing a richer metastability landscape. The systematic scan across initial-data classes is a positive feature. However, the numerical character of the study means its significance hinges on validation of the discretization.

major comments (1)
  1. [Abstract and numerical results] Abstract and numerical results: the central claim that UV-dominated configurations remain stable significantly longer depends on the spacetime finite element discretization faithfully reproducing the continuous dynamics. For systems with opposite-sign kinetic terms (non-hyperbolic), standard FEM stabilization can preferentially damp the high-frequency modes that mediate the nonlinear energy transfer described in the abstract. No mesh-convergence data, a-posteriori error estimates, or cross-validation against alternative discretizations are reported, so the observed metastability windows could be numerical artifacts rather than properties of the continuous theory.
minor comments (2)
  1. The abstract would benefit from a brief statement of the specific potentials and initial-data families used, to make the trends immediately quantifiable.
  2. Figure captions should explicitly state the simulation duration and mesh parameters for each panel to allow readers to assess the reported stability intervals.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and the constructive comment on numerical validation. We address the concern directly below and will strengthen the manuscript accordingly.

read point-by-point responses

  1. Referee: [Abstract and numerical results] Abstract and numerical results: the central claim that UV-dominated configurations remain stable significantly longer depends on the spacetime finite element discretization faithfully reproducing the continuous dynamics. For systems with opposite-sign kinetic terms (non-hyperbolic), standard FEM stabilization can preferentially damp the high-frequency modes that mediate the nonlinear energy transfer described in the abstract. No mesh-convergence data, a-posteriori error estimates, or cross-validation against alternative discretizations are reported, so the observed metastability windows could be numerical artifacts rather than properties of the continuous theory.

    Authors: We agree that explicit convergence diagnostics are required to substantiate the claims for a non-hyperbolic system. In the revised version we will add mesh-refinement studies in both (1+1)D and (2+1)D, demonstrating that the reported metastability lifetimes for UV-dominated, low-amplitude data converge under successive h-refinement. We will also include a-posteriori residual-based error estimates for representative runs. Our spacetime Galerkin formulation employs no artificial viscosity or high-mode damping; the weak form is consistent and the time integrator is implicit and energy-consistent. To address possible scheme dependence we will include a short cross-check against a second-order centered finite-difference discretization on a subset of initial data, confirming that the qualitative distinction between UV- and IR-dominated regimes persists. These additions will be placed in a new subsection of the numerical-methods section together with the existing spectral diagnostics. revision: yes

Circularity Check

0 steps flagged

No significant circularity in numerical study

full rationale

The paper performs direct numerical simulations of ghost-normal scalar field systems using a spacetime finite element discretization in (1+1) and (2+1) dimensions. All reported stability windows, spectral dependence, and metastability observations are generated by evolving initial data under the discretized equations of motion; no parameters are fitted to subsets of results and then invoked as predictions, no self-definitional equations appear, and no load-bearing uniqueness theorems or ansatzes are imported via self-citation. The central claims therefore rest on computational output rather than any derivation that reduces to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that the chosen numerical scheme accurately captures the continuous ghost dynamics and that the observed bounded intervals are not artifacts of finite simulation time or discretization.

axioms (1)
  • domain assumption The spacetime finite element method converges to the true solution of the continuous field equations as mesh size and time step approach zero.
    Invoked implicitly when interpreting simulation outputs as physical stability statements.

pith-pipeline@v0.9.0 · 5492 in / 1272 out tokens · 35417 ms · 2026-05-13T07:58:04.184444+00:00 · methodology

Review history (2 revisions) →

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

65 extracted references · 65 canonical work pages · 1 internal anchor

  1. [1]

    We follow the dis- cretization scheme in [23, 24]

    Space-Time Finite Element Method A space-time FEM uses continuous approximation functions in both space and time. We follow the dis- cretization scheme in [23, 24]. In this scheme, space and time are discretized together for the entire domain using a finite element space which does not discriminate be- tween space and time basis functions. In this way, we...

    work page

  2. [2]

    Letu=∂ tϕandv=∂ tχ

    Setting Up The Weak Form Since the original system is second order in time, it is necessary for the weak formulation to introduce auxiliary variables such that our system can be reduced to one that is first order in time. Letu=∂ tϕandv=∂ tχ. From Eqns. 2 and 3, we 4 have ∂u ∂t − ∇2ϕ+m 2 ϕϕ+ ∂V ∂ϕ = 0,(16) ∂ϕ ∂t −u= 0,(17) ∂v ∂t − ∇2χ+m 2 χχ+γ ∂V ∂χ = 0,(1...

    work page

  3. [3]

    Nonlinear Solver: PETSc SNES In order to computationally handle the nonlinearity of the PDE system, we utilize PETSc’s SNES library. The library contains methods, like Newton’s method with line search, for solving nonlinear equations of the form F(U) = 0, whereFis the nonlinear differential operator, andUis the solution vector. The two main user-created c...

    work page

  4. [4]

    R k ←F(U k)(FormResidual()) 2.If∥R k∥ ≤rtol,stop (converged)

    work page

  5. [5]

    J k ← ∂F ∂U (Uk)(FormJacobian()) 4.SolveJ k δUk =−R k (KSP/PC linear solve, PETSc) 5.Choose stepα∈(0,1] (line search, PETSc)

    work page

  6. [6]

    Initial Data In this section, we describe different initial conditions (IC) that are motivated by different physical implica- 5 tions

    U k+1 ←U k +α δU k (PETSc) TABLE I B. Initial Data In this section, we describe different initial conditions (IC) that are motivated by different physical implica- 5 tions. We first start with simple plane wave and Gaussian packet IC which are similarly explored in [21]. Then, we expand the discussion and explore more interesting sce- narios

    work page

  7. [7]

    where cϕ =± s k2 ϕ +m 2 ϕ k2 ϕ , cχ =± s k2χ +m 2χ k2χ

    Plane waves The plane wave initial conditions are given by ϕ0(x) =A ϕ sin kϕ(x−x ϕ) , ˙ϕ0(x) =−c ϕ kϕ Aϕ cos kϕ(x−x ϕ) , χ0(x) =A χ sin kχ(x−x χ) , ˙χ0(x) =−c χ kχ Aχ cos kχ(x−x χ) . where cϕ =± s k2 ϕ +m 2 ϕ k2 ϕ , cχ =± s k2χ +m 2χ k2χ . (21) We fixx ϕ = 0,x χ =L/3,A ϕ =A χ =A,k ϕ =k χ/2 = k, and choose the two plane waves to move in opposite directions...

    work page

  8. [8]

    The initial positions are described byx ϕ = 0.3Landx χ = 0.7Land the amplitudes and lengths are kept consistent throughout,A≡A ϕ =A χ andℓ≡ℓ ϕ =ℓ χ

    Gaussian packets The Gaussian wave packet initial data are defined by ϕ0(x) =A ϕ exp −(x−x ϕ)2 2ℓ2 ϕ ! , ˙ϕ0(x) =A ϕ cϕ(x−x ϕ) ℓ2 ϕ exp −(x−x ϕ)2 2ℓ2 ϕ ! , χ0(x) =A χ exp −(x−x χ)2 2ℓ2χ , ˙χ0(x) =A χ cχ(x−x χ) ℓ2χ exp −(x−x χ)2 2ℓ2χ , wherec ϕ andc χ are defined oppositely (1 and -1) so that the waves travel in opposite directions or defined equivalently ...

    work page

  9. [9]

    integrable turbulence

    Colored-noise spectra (IR/UV-tilted) Choose integersn 1 ≤n≤n 2, phasesθ n ∈[0,2π), and a spectral tiltn s ∈R. Define ϕ(x,0) =C ϕ n2X n=n1 k ns/2 n cos knx+θ n ,(22) ˙ϕ(x,0) =− n2X n=n1 sn ωϕ(kn) kn ∂x h Cϕ k ns/2 n cos knx+θ n i , (23) wheres n ∈ {+1,−1}selects co-/counter-propagating content. Set Cϕ = Aq 1 2 Pn2 n=n1 k ns n ,(24) so that RMSx[ϕ] =A(for r...

    work page

  10. [10]

    This is the cleanest way to control initial cross- correlation between sectors

    Phase-correlated two-field plane waves For a single carrierk, ϕ(x,0) =Acos k(x−x ϕ) ,(25) ˙ϕ(x,0) =− ωϕ(k) k ∂xϕ(x,0),(26) χ(x,0) =A rcos k(x−x χ) + ∆ϕ ,(27) ˙χ(x,0) =−σ ωχ(k) k ∂xχ(x,0),(28) with amplitude ratior >0, relative phase ∆ϕ∈[0, π], andσ∈ {+1,−1}for co-/counter-propagation. This is the cleanest way to control initial cross- correlation between ...

    work page

  11. [11]

    Oscillons are long-lived, localized, nearly periodic lumps that appear in many scalar theories [44–46]

    Oscillon-like, time-symmetric seeds ϕ(x,0) =Asech x−x 0 σ , ˙ϕ(x,0) = 0,(29) χ(x,0) =A rsech x−x 0 σ cos ∆ϕ,˙χ(x,0) = 0.(30) To add a carrier, we may multiply each profile by cos(k0(x−x 0)). Oscillons are long-lived, localized, nearly periodic lumps that appear in many scalar theories [44–46]. Seed- ing an oscillon-like profile probes whether the ghostly ...

    work page

  12. [12]

    Each sim- ulation was run until a time step numerically diverged

    Plane Wave Initial Conditions The plane wave tests sweep over the field initial am- plitudeAand wavenumber ratioC=kL/2π. Each sim- ulation was run until a time step numerically diverged. Thus, we definet long-lived as the final converged time. We see from table II that systems are much longer lived if the initial amplitude is decreased. TABLE II: Plane-wa...

    work page

  13. [13]

    Three widths were considered, and for each one, the amplitude was varied

    Gaussian packets Gaussian packet initial conditions have tunable ampli- tudeAand widthℓ. Three widths were considered, and for each one, the amplitude was varied. Note that since k= 1/(4ℓ), and we have definedC=kL/2π, we tune Crather thanℓdirectly. The packet centers were fixed atx ϕ = 0.3Landx χ = 0.7L. We look at scenarios with co-/counter-propagating p...

    work page

  14. [14]

    We sets n =−1 for consistent counter propagation

    Colored-Noise Spectra Colored-noise initial data were generated with different spectral tiltsn s andn 2 values. We sets n =−1 for consistent counter propagation. From table VIII, we find that a more negative tilt (ns <0) leads to earlier instabilities. Recall that a more negative tilt concentrates power in the IR range, and a more positive tilt biases tow...

    work page

  15. [15]

    We keep initial amplitudeA= 1 and wave number ratioC=kL/2π= 1.0 for all runs

    Phase-Correlated two-field Plane Waves For the phase-correlated plane wave initial data, we vary the relative phase ∆ϕ, the amplitude ratiorbetween the two fields, and the co-/counter-propagation which is controlled byσ. We keep initial amplitudeA= 1 and wave number ratioC=kL/2π= 1.0 for all runs. From tables X and XI, we find that the long lived-ness, tl...

    work page

  16. [16]

    Tests were performed both with (XIII) and without (XII) an added spatial carrier (k0L/2π= 1)

    Oscillon-Like, Time-Symmetric Seeds These tests explore the oscillon-like seeds characterized by widthσand amplitudeA. Tests were performed both with (XIII) and without (XII) an added spatial carrier (k0L/2π= 1). These results show us that instability is reached sooner as the width,σ, increases. FIG. 4: Phase-correlated plane waves initialization case. Th...

    work page

  17. [17]

    J. M. Cline, S. Jeon, and G. D. Moore, The phantom menaced: Constraints on low-energy effective ghosts, Phys. Rev. D70, 043543 (2004), arXiv:hep-ph/0311312

    work page Pith review arXiv 2004

  18. [18]

    R. P. Woodard, Avoiding dark energy with 1/rmodi- fications of gravity, Lect. Notes Phys.720, 403 (2007), arXiv:astro-ph/0601672

    work page Pith review arXiv 2007

  19. [19]

    R. R. Caldwell, A phantom menace? cosmological consequences of a dark energy component with super- negative equation of state, Phys. Lett. B545, 23 (2002), arXiv:astro-ph/9908168

    work page Pith review arXiv 2002

  20. [20]

    K. S. Stelle, Renormalization of higher derivative quan- tum gravity, Phys. Rev. D16, 953 (1977)

    work page 1977

  21. [21]

    K. S. Stelle, Classical gravity with higher derivatives, Gen. Rel. Grav.9, 353 (1978)

    work page 1978

  22. [22]

    J. Z. Simon, Higher-derivative lagrangians, nonlocality, problems and solutions, Phys. Rev. D41, 3720 (1990)

    work page 1990

  23. [23]

    Ostrogradsky, M´ emoires sur les ´ equations diff´ erentielles, relatives au probl` eme des isop´ erim` etres, Mem

    M. Ostrogradsky, M´ emoires sur les ´ equations diff´ erentielles, relatives au probl` eme des isop´ erim` etres, Mem. Acad. St. Petersbourg6, 385 (1850)

    work page

  24. [24]

    The galileon as a local modification of gravity

    A. Nicolis, R. Rattazzi, and E. Trincherini, The galileon as a local modification of gravity, Phys. Rev. D79, 064036 (2009), arXiv:0811.2197

    work page Pith review arXiv 2009

  25. [25]

    Degenerate higher derivative theories beyond Horndeski: evading the Ostrogradski instability

    D. Langlois and K. Noui, Degenerate higher derivative theories beyond horndeski: evading the ostrogradski in- stability, JCAP1602, 034, arXiv:1510.06930

    work page Pith review arXiv

  26. [26]

    C. P. Burgess, Introduction to effective field theory, Ann. Rev. Nucl. Part. Sci.57, 329 (2007), arXiv:hep- th/0701053

    work page arXiv 2007

  27. [27]

    J. F. Donoghue, General relativity as an effective field theory: The leading quantum corrections, Phys. Rev. D 50, 3874 (1994), arXiv:gr-qc/9405057

    work page Pith review arXiv 1994

  28. [28]

    Can dark energy evolve to the Phantom?

    A. Vikman, Can dark energy evolve to the phantom?, Phys. Rev. D71, 023515 (2005), arXiv:astro-ph/0407107

    work page internal anchor Pith review arXiv 2005

  29. [29]

    Creminelli, G

    P. Creminelli, G. D’Amico, J. Norena, L. Senatore, and F. Vernizzi, The effective theory of quintessence: thew < −1 side unveiled, JCAP0902, 018, arXiv:0811.0827

    work page arXiv

  30. [30]

    Deffayet, S

    C. Deffayet, S. Mukohyama, and A. Vikman, Ghosts without runaway instabilities, Phys. Rev. Lett.128, 041301 (2022)

    work page 2022

  31. [31]

    Deffayet, A

    C. Deffayet, A. Held, S. Mukohyama, and A. Vikman, Global and local stability for ghosts coupled to positive energy degrees of freedom, Journal of Cosmology and As- troparticle Physics2023(11), 031

    work page

  32. [32]

    Errasti D´ ıez, J

    V. Errasti D´ ıez, J. Gaset Rif` a, and G. Staudt, Foundations of ghost stability, Fortschritte der Physik73, 2400268 (2025), https://onlinelibrary.wiley.com/doi/pdf/10.1002/prop.202400268

    work page doi:10.1002/prop.202400268 2025

  33. [33]

    D. R. Noakes, The initial value formulation of higher derivative gravity, Journal of Mathematical Physics24, 1846 (1983)

    work page 1983

  34. [34]

    Held and H

    A. Held and H. Lim, Nonlinear dynamics of quadratic gravity in spherical symmetry, Phys. Rev. D104, 084075 (2021)

    work page 2021

  35. [35]

    Held and H

    A. Held and H. Lim, Nonlinear evolution of quadratic gravity in 3 + 1 dimensions, Phys. Rev. D108, 104025 (2023)

    work page 2023

  36. [36]

    Held and H

    A. Held and H. Lim, Black-hole binaries and waveforms in quadratic gravity (2025), arXiv:2503.13428 [gr-qc]

    work page arXiv 2025

  37. [37]

    Deffayet, A

    C. Deffayet, A. Held, S. Mukohyama, and A. Vikman, Ghostly interactions in (1 + 1)-dimensional classical field theory, Phys. Rev. D112, 065011 (2025)

    work page 2025

  38. [38]

    Balay, S

    S. Balay, S. Abhyankar, M. F. Adams, S. Benson, J. Brown, P. Brune, K. Buschelman, E. M. Constan- tinescu, L. Dalcin, A. Dener, V. Eijkhout, J. Faibus- sowitsch, W. D. Gropp, V. Hapla, T. Isaac, P. Jo- livet, D. Karpeev, D. Kaushik, M. G. Knepley, F. Kong, S. Kruger, D. A. May, L. C. McInnes, R. T. Mills, L. Mitchell, T. Munson, J. E. Roman, K. Rupp, P. S...

    work page 2025

  39. [39]

    Anderson and J.-H

    M. Anderson and J.-H. Kimn, A numerical approach to space-time finite elements for the wave equation, Journal of Computational Physics226, 466 (2007)

    work page 2007

  40. [40]

    D. A. French and T. E. Peterson, A continuous space- time finite element method for the wave equation, Math. Comput.65, 491–506 (1996)

    work page 1996

  41. [41]

    Cherubini and S

    C. Cherubini and S. Filippi, Using femlab for gravita- tional problems: numerical simulations for all, Journal- Korean Physical Society49(2005)

    work page 2005

  42. [42]

    Cherubini, F

    C. Cherubini, F. Federici, S. Succi, and M. P. Tosi, Ex- cised acoustic black holes: The scattering problem in the time domain, Phys. Rev. D72, 084016 (2005)

    work page 2005

  43. [43]

    Z. Cao, P. Fu, L.-W. Ji, and Y. Xia, Binary black hole simulation with an adaptive finite element method ii: Ap- plication of local discontinuous galerkin method to ein- stein equations (2018), arXiv:1805.10640 [gr-qc]

    work page arXiv 2018

  44. [44]

    H. Lim, M. Anderson, and J.-H. Kimn, Exploring critical collapse in the semilinear wave equation using space-time finite elements (2020), arXiv:1410.0708 [gr-qc]

    work page arXiv 2020

  45. [45]

    Kim, On the Numerical Computation for Solv- ing the Two-Dimensional Parabolic Equations by Space- Time Finite Element Method, JSME International Jour- nal Series B44, 434 (2001)

    C.-K. Kim, On the Numerical Computation for Solv- ing the Two-Dimensional Parabolic Equations by Space- Time Finite Element Method, JSME International Jour- nal Series B44, 434 (2001)

    work page 2001

  46. [46]

    Guddati and J

    M. Guddati and J. Tassoulas, Space-time finite elements for the analysis of transient wave propagation in un- bounded layered media, International Journal of Solids and Structures36, 4699 (1999)

    work page 1999

  47. [47]

    Dyniewicz, Space-time finite element approach to gen- eral description of a moving inertial load, Finite Elem

    B. Dyniewicz, Space-time finite element approach to gen- eral description of a moving inertial load, Finite Elem. Anal. Des.62, 8–17 (2012)

    work page 2012

  48. [48]

    Zank, Efficient direct space-time finite element solvers for the wave equation in second-order formulation, J

    M. Zank, Efficient direct space-time finite element solvers for the wave equation in second-order formulation, J. Sci. Comput.105, 10.1007/s10915-025-03038-1 (2025)

    work page doi:10.1007/s10915-025-03038-1 2025

  49. [49]

    Schwing, I

    A. Schwing, I. Nompelis, and G. V. Candler, Im- plementation of adaptive mesh refinement in an im- plicit unstructured finite-volume flow solver, in21st AIAA Computational Fluid Dynamics Conference, https://arc.aiaa.org/doi/pdf/10.2514/6.2013-2446

    work page doi:10.2514/6.2013-2446 2013

  50. [50]

    X. Xu, Y. Chen, Z. Han, and F. Zhou, A total energy- based adaptive mesh refinement technique for the simu- lation of compressible flow, Frontiers in Energy Research Volume 11 - 2023, 10.3389/fenrg.2023.1203801 (2023). [35]Adaptive Mesh Refinement for DDES Sim- ulation on Transonic Compressor Cas- cade With Unstructured Mesh, Turbo Expo, Vol. Volume 2C: T...

    work page doi:10.3389/fenrg.2023.1203801 2023

  51. [51]

    G. L. Bryan, M. L. Norman, B. W. O’Shea, T. Abel, J. H. Wise, M. J. Turk, D. R. Reynolds, D. C. Collins, P. Wang, S. W. Skillman, B. Smith, R. P. Harkness, J. Bordner, J.-h. Kim, M. Kuhlen, H. Xu, N. Goldbaum, C. Hummels, A. G. Kritsuk, E. Tasker, S. Skory, C. M. Simpson, O. Hahn, J. S. Oishi, G. C. So, F. Zhao, R. Cen, and Y. Li, Enzo: An adaptive mesh r...

    work page 2014

  52. [52]

    G. Xu, Hydrodynamic and n-body schemes on an un- structured, adaptive mesh with applications to cosmo- logical simulations, Monthly Notices of the Royal Astro- nomical Society288, 903–919 (1997)

    work page 1997

  53. [53]

    Mazurek, M

    E. Evans, S. Iyer, E. Schnetter, W.-M. Suen, J. Tao, R. Wolfmeyer, and H.-M. Zhang, Computa- tional relativistic astrophysics with adaptive mesh refine- ment: Testbeds, Physical Review D71, 10.1103/phys- revd.71.081301 (2005)

    work page doi:10.1103/phys- 2005

  54. [54]

    Nazarenko,Wave Turbulence, Lecture Notes in Physics, Vol

    S. Nazarenko,Wave Turbulence, Lecture Notes in Physics, Vol. 825 (Springer, 2011)

    work page 2011

  55. [55]

    Walczak, S

    P. Walczak, S. Randoux, and P. Suret, Optical rogue waves in integrable turbulence, Phys. Rev. Lett.114, 143903 (2015)

    work page 2015

  56. [56]

    Randoux and P

    S. Randoux and P. Suret, Integrable turbulence: a review of recent results in optics (2016), preprint

    work page 2016

  57. [57]

    Kofman, A

    L. Kofman, A. Linde, and A. Starobinsky, Towards the theory of reheating after inflation, Phys. Rev. D56, 3258 (1997), arXiv:hep-ph/9704452

    work page arXiv 1997

  58. [58]

    P. B. Greene, L. Kofman, A. Linde, and A. A. Starobin- sky, Structure of resonance in preheating after inflation, Phys. Rev. D56, 6175 (1997), arXiv:hep-ph/9705347

    work page arXiv 1997

  59. [59]

    Gleiser and D

    M. Gleiser and D. Sicilia, General theory of oscillon dy- namics, Phys. Rev. D80, 125037 (2009)

    work page 2009

  60. [60]

    M. A. Amin and D. Shirokoff, Flat-top oscillons in an expanding universe, Phys. Rev. D81, 085045 (2010)

    work page 2010

  61. [61]

    Hindmarsh and P

    M. Hindmarsh and P. Salmi, Oscillons and domain walls, Phys. Rev. D77, 105025 (2008). 16 Appendix A: Element Stiffness Matrix Calculations An important step in any FEM is determining the el- ement stiffness matrices. These matrices are found by solving the weak form of the equations of motion, Eq

    work page 2008

  62. [62]

    The weak equations are displayed in equation 20, but in equations A1, A2, A3, and A4 we display them in their (1 + 1) form

    Solving the space-time FEM for the (1 + 1) case re- sults in 4x4 matrices, the (2 + 1) case finds 8x8 matrices, and the 3+1 case has 16x16 matrices. The weak equations are displayed in equation 20, but in equations A1, A2, A3, and A4 we display them in their (1 + 1) form. K1 = Z X,T utΨ +ϕ xΨx +m 2 ϕϕΨ +V ϕΨ dxdt(A1) K2 = Z X,T (ϕtΨ−uΨ)dxdt(A2) G1 = Z X,T...

    work page

  63. [63]

    A5, we con- sider the four corners of a rectangle in the (x, t) plane, see Fig

    Element Basis Function The rectangular element basis function is a piecewise polynomial with four terms, Ψi(x, t) =b 1,ix+b 2,it+b 3,ixt+b 4,i.(A5) To solve for the unknown coefficients in Eq. A5, we con- sider the four corners of a rectangle in the (x, t) plane, see Fig. 11. (0,0) (hx,0) (0,ht) (hx,ht) FIG. 11 In figure 11,h x andh t represent the step s...

    work page

  64. [64]

    We will explain the set up of the element stiffness matrices, display the matri- ces, and then discuss practical implementation into the simulation

    Element Stiffness Matrices Now that we have the basis functions, the element stiff- ness matrices can be determined. We will explain the set up of the element stiffness matrices, display the matri- ces, and then discuss practical implementation into the simulation. 17 In the two dimensional (1 + 1) case, the weak form of the equations of motion, equations...

    work page

  65. [65]

    Here,λis a parameter controlling the strength of the nonlinear potential

    Manufactured Solution Tests We consider the following coupled partial differential equations (PDEs) in (1 + 1) dimensions: ϕtt −ϕ xx +ϕ+V ϕ(ϕ, χ) = 0,(B1) χtt −χ xx +χ−V χ(ϕ, χ) = 0,(B2) where the highly nonlinear potential terms are defined as Vϕ(ϕ, χ) =−2λ ϕ A B1.5 Vχ(ϕ, χ) = 2λ χ C B1.5 with A=ϕ 2 −χ 2 + 1 B= ϕ2 −χ 2 −1 2 + 4ϕ2, C=ϕ 2 −χ 2 −1. Here,λis...

    work page