Solving the nonlinear Klein-Gordon equation: semianalytical Galerkin method

Annibal D. de Figueiredo Neto; Caio C. Holanda Ribeiro; Luana L. Silva Ribeiro

REVIEW 2 major objections 4 minor 20 references

Solving the nonlinear Klein-Gordon equation: semianalytical Galerkin method

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

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

Pith's one-line read A Galerkin truncation of the nonlinear Klein-Gordon equation yields finite mechanical systems whose stationary points reproduce the field's exact stationary solutions, and whose dynamics tracks smooth time evolution with a computable residu

desk verdict Honest, correct paper: exact elliptic-function stationary solutions on an interval are the real contribution; the time-dependent Galerkin convergence is a numerical hint, not a proof, and the authors mostly say so. read the letter →

arxiv 2509.02925 v1 pith:YOAEAKEX submitted 2025-09-03 math-ph math.MP

Annibal D. de Figueiredo Neto , Caio C. Holanda Ribeiro , Luana L. Silva Ribeiro This is my paper · ORCID

classification math-phmath.MP MSC 35L7035C0535Q40

keywords nonlinearKlein-GordonequationGalerkinmethodJacobiellipticfunctionsstationarysolutionsmechanicalsystemsDirichletboundaryconditionsMexican-hatpotentialLagrangiantruncation

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

This paper sets out to show that a nonlinear field theory—a real scalar field in 1+1 dimensions with Dirichlet boundary conditions and a Mexican-hat potential—can be replaced, for many purposes, by a mechanical system of a few particles. The equivalence is built at the Lagrangian level, so every truncated system automatically preserves the conserved quantities of the field. The authors derive exact stationary solutions of the nonlinear Klein-Gordon equation in terms of Jacobi elliptic functions and show these correspond to the stationary points of the finite mechanical systems. They also define an explicit local error for the time-dependent approximation and give numerical evidence that the error shrinks as the number of particles grows. If correct, the method offers a controlled semianalytical route into a nonlinear field equation without relying on perturbation theory.

What carries the argument

The carrying mechanism is the Galerkin expansion φ(t,x)=Σ_n A_n(t) sin(nπx/ℓ), inserted directly into the field's Lagrangian. This turns the field theory into an infinite chain of particles with coordinates A_n, harmonic frequencies n²+λ, and quartic couplings governed by the tensor D_nmpq. Truncating at N modes yields the finite Lagrangian and equations of motion (45). Two ingredients make the analysis work: the Jacobi elliptic functions sn(u,k) and cn(u,k), whose second-order differential identities match the field equation exactly and produce the stationary solutions, and the potential U^(N), whose critical points are the stationary configurations of the truncated system. Because the appr

What would settle it

Run the truncated system with a shock-like or high-gradient initial condition and compute R^(N) for increasing N (say 10, 20, 40, 80); if the total error stops decreasing or grows with N, the claimed convergence fails. A cheaper probe is to extend the smooth initial-condition simulation of Fig. 7 beyond τ=10 and check whether A_1^(10) and A_1^(20) remain close, since the paper's evidence stops at τ=10.

Watch

Extended reading notes

Core claim

The central claim is that truncating the Galerkin expansion of the nonlinear Klein-Gordon equation at N modes produces a finite, Hamiltonian mechanical system whose stationary points approximate the analytic stationary solutions of the field and whose dynamics approximates the field's Cauchy problem. For V(φ)=β/4(φ²−φ₀²)², the stationary solutions are written explicitly: sn-type solutions when λ<0, with only finitely many nontrivial states, and cn-type solutions when λ>0, with infinitely many. The paper derives the algebraic conditions that fix the parameters of these elliptic functions from the Dirichlet boundary conditions, then defines the truncated field φ^(N) and the residual R^(N) (Eq.

Load-bearing premise

The Galerkin truncation converges to the true time-dependent field dynamics; the numerical evidence covers only smooth initial data over a finite time window and the authors state the approximation may fail once shock or other types of waves form.

Editorial extensions

If this is right

For smooth initial data with only a few Fourier modes excited, a handful of ordinary differential equations can replace the nonlinear PDE over moderate time scales; the simulations show N=10 and N=20 already track each other closely over 0<τ<10.
The exact Jacobi-elliptic stationary solutions give an analytic catalog of equilibrium configurations for any value of the potential parameter λ, including the inverted-potential case where infinitely many static solutions exist.
Because the approximation is made at the Lagrangian level, each truncated system is Hamiltonian and energy is conserved exactly, avoiding the conservation-law violations that can afflict direct discretizations of the field equation.
The explicit residual R^(N) provides a checkable error estimate, allowing one to estimate how many Galerkin modes are needed for a given initial condition and integration time.

Reading between the lines

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

Extension: for other boundary conditions, such as Neumann, or for other nonlinear potentials, the same Lagrangian truncation should apply, but convergence is likely to be slower if the basis functions do not respect the boundary conditions; the paper notes the Neumann case explicitly.
Extension: the correspondence between stationary solutions of the field and critical points of the truncated potential is special to this potential and basis; a natural test is whether spurious critical points appear for generic nonlinear potentials.
Extension: the reduced systems, particularly the N=3 case with stable and unstable fixed points, could be used to study energy transfer or prethermalization in the field, a direction the paper does not pursue.

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.

Desk Editor's Note

read the letter

Bottom line: this is a careful, honest paper that does one genuinely new thing—exact stationary solutions of the nonlinear KG equation on a finite interval with Dirichlet boundary conditions, written in terms of Jacobi elliptic functions, with a count of how many exist for each λ—and then shows that the Galerkin-truncated mechanical systems have those exact solutions as stationary points. The elliptic-function work is real; I checked the algebra for sn and cn, and the matching conditions (Eqs. 28–32, 36–41) are correct. They do not pretend the traveling-wave ansatz is new; they cite Ates and Inc. The formula for D_nmpq in Eq. (14) is quoted without derivation, which is a small gap, but it is checkable and the tables in III and IV line up with it.

The soft spot is the time-dependent claim. The paper honestly calls it "hints toward convergence," and the residual R^(N) in Eq. (50) is a residual, not an error bound. A small residual does not imply a small solution error without an a priori or a posteriori argument, and the authors do not supply one. The numerics cover one initial condition (first four modes excited), one value λ = −10, and 0 < τ < 10. That is enough to support a conjecture, not a theorem, and the stress-test's worry about λ > 0 is legitimate: with the potential unbounded below, the linearized error equation can be unstable, so residual control would need to be even more careful. The authors themselves note that the approximation may fail for shocks. So the time-dependent part is a limitation they acknowledge, not an overclaim.

What else? The Lagrangian-level Galerkin truncation conserves energy by Noether, which is a genuine advantage over some finite-difference schemes (they cite the Gross-Pitaevskii problem). The qualitative analysis of the N = 3 system, with the potential plots, is nice pedagogy but not deep.

I would send this to a referee. It is clean, honest, and the stationary part is new and correct. The main revision request would be: (1) add a derivation of Eq. (14), (2) run a wider convergence study—several initial conditions, longer times, both signs of λ, and maybe an independent numerical solver as benchmark—and (3) either state clearly that the time-dependent convergence is empirical or give a real error estimate. None of this is fatal; the paper already frames itself as semianalytical and does not claim a proof.

Reader: anyone working on spectral methods for nonlinear field equations, or with an interest in mechanical analogies for scalar fields. I would bring it to a reading group as a good example of how to present a Galerkin method with exact checks.

Referee Report

2 major / 4 minor

Summary. The manuscript develops a Galerkin reduction of the 1+1 nonlinear Klein-Gordon equation with a Mexican-hat potential and Dirichlet boundary conditions. The field is expanded in sine modes, yielding an infinite mechanical system with quartic mode couplings. The authors derive exact stationary solutions in terms of Jacobi elliptic functions for both signs of the nonlinearity, after a one-parameter rescaling λ. They then truncate the mechanical system to N degrees of freedom and show numerically that its stationary points reproduce the exact stationary solutions. For the time-dependent Cauchy problem, they define a local residual R^(N) and present simulations for λ=-10 and one initial condition that excite only the first four modes, arguing that the method converges as N increases.

Significance. If established, the Lagrangian Galerkin framework is attractive: it preserves conservation laws by construction, gives exact stationary solutions in closed form, and reduces the field theory to small mechanical systems that can be studied with ODE tools. The stationary-solution part is largely sound and is a useful contribution. The explicitly stated limitation that the approximation may fail for shocks or other nonsmooth waves, together with the absence of an error bound, means the time-dependent claim is only a numerical hint. The incorrect formula for the mode-coupling tensor in Eq. (14) must be fixed before the paper can be relied upon.

major comments (2)

[§II, Eq. (14)] The formula for D_{nmpq} is incorrect as written. For n=m=p=q=1, direct integration with φ_n = sqrt(2/ℓ) sin(nπx/ℓ) gives D_1111 = ℓ ∫ φ_1^4 dx = 3/2, while Eq. (14) yields -1/ℓ². The expression is also dimensionally inconsistent. The correct formula involves Kronecker deltas of |n-m| with |p-q| and with p+q, with special handling of the zero mode; the printed N_{n-m} and N_{n+m} prefactors are not the correct zero-mode factors. Since D enters the Lagrangian (16), the equations of motion (18), and the truncated systems (45), this is a central error and must be corrected and accompanied by a derivation.
[§IV.B, Eq. (50)] The residual R^(N) is not shown to control the error between the Galerkin solution ϕ^(N) and the exact solution of Eq. (49). A small residual does not imply a small solution error without a stability or continuity estimate for the nonlinear evolution operator. For λ>0 the potential (4) is unbounded below, so the linearized error equation can be unstable; the provided numerical evidence is restricted to λ=-10, a single initial condition with only the first four modes excited, and the interval 0<τ<10. The authors themselves state in §IV that the approximation "might not hold when shock or other types of waves exist." Given these limitations, the claim in §IV.B that the method "provides a good approximation" is stronger than the evidence supports. Either an a posteriori error bound, or a systematic convergence study varying λ, initial conditions, and time scales, is needed to justify the ti

minor comments (4)

[§II, Eq. (14)] No derivation of D_{nmpq} is given; it should at least be sketched or referenced, especially since the printed formula is wrong.
[§II, Sec. IV.B] The definition of R^(N) in Eq. (50) uses a norm without specifying which norm is meant; later Rbar^(N) is integrated, so the notation should be clarified.
[Figs. 2 and 4] The axes are not labeled; the curves representing multiples of 2K/π and the λ-dependent functions should be identified directly or in the captions.
[General] There are minor typos, e.g., "is β <0" on p. 2. The paper would also benefit from a statement on reproducibility; no code or numerical parameter details are provided for the simulations in Figs. 7–9.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: stationary solutions follow from a direct elliptic-function ansatz, and the Galerkin finite systems are independent approximations whose agreement is verified numerically.

full rationale

The paper's central derivation is self-contained. The exact stationary solutions of Eqs. (26)/(34) are obtained by substituting a Jacobi elliptic function ansatz into the stationary ODE (21) and imposing Dirichlet boundary conditions; no fitted parameter or target solution is fed back into the construction (the constants k1,k2,k3,k4 and q1,...,q4 are fixed by the differential identities (27)/(35) and the boundary conditions (30)/(39)). The finite mechanical systems of Sec. IV are obtained by truncating the same Lagrangian (16) at N modes, so their stationary points are, by construction, Galerkin approximations to the exact solutions; Tables I-VI and Fig. 6 compare the two independently computed quantities and show convergence, which is verification rather than circularity. The time-dependent claim is supported only by numerical residual/evolution studies (Eq. (50), Figs. 7-9), and the authors explicitly concede in Sec. IV that the approximation 'might not hold when shock or other types of waves exist'; that is an unproved-convergence limitation, not a circular derivation. The only self-citation ([7], by co-author Ribeiro) appears in the introduction as a general reference to semianalytical methods and is not load-bearing for any result. Therefore no step in the derivation reduces to its own input.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central claim rests on standard spectral theory and elliptic function identities, plus the untested assumption that truncation converges for time-dependent solutions. No free parameters are fitted; the elliptic modulus is determined by solving a transcendental equation. No new entities are introduced.

assumptions (4)

standard math The sine functions sin(nπx/l) form a complete basis for functions with Dirichlet BC on [0,l], and interactions are computed via the D tensor.
Invoked in Section II A for the Galerkin expansion; standard Sturm-Liouville completeness.
standard math Jacobi elliptic function identities (d² sn/du² and d² cn/du²) and zero locations are correct.
Used in Section III to verify the ansatz; cited to standard references.
domain assumption The classical field equation is well defined even when the potential is unbounded below (β<0), with runaway solutions ignored.
Section II, paragraph after Eq. (4), 'runaway solutions are of no concern'; this is a physical assumption that the field stays in a bounded region.
ad hoc to paper The Galerkin truncation converges for smooth solutions of the time-dependent problem.
Assumed in Section IV B and supported only by numerical evidence; the authors explicitly state it may fail for shocks.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Solving the nonlinear Klein-Gordon equation: semianalytical Galerkin method." pith.science (2026). https://pith.science/paper/YOAEAKEX

@misc{pith2026250902925,
  author       = {Pith},
  title        = {Pith review of: Solving the nonlinear Klein-Gordon equation: semianalytical Galerkin method},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YOAEAKEX}},
  note         = {Machine review of arXiv:2509.02925}
}

abstract

In this work, approximate solutions to the nonlinear Klein-Gordon equation are constructed by means of the Galerkin method. Specifically, it is shown how the dynamics of a real scalar field in $1+1$ dimensions subjected to Dirichlet boundary conditions and Mexican-hat-like potentials can be approximated by mechanical systems with a few particles. Because the approximation is performed at a Lagrangian level, one of the advantages of this method is the control over conservation laws present in the field theory, which are captured by the finite mechanical systems. Among the results, exact stationary solutions for the nonlinear KG equation are found in terms of Jacobi elliptic functions, which are shown to correspond to stationary configurations of the mechanical systems. Furthermore, numerical simulations are provided, giving hints towards the convergence of the method.

Figures

Figures reproduced from arXiv: 2509.02925 by the authors.

**Figure 2.** FIG. 2. Graphical determination of the values of [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 4.** FIG. 4. Representation of the solutions for [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗

Figures from the paper (6 more)

**Figure 5.** Figure 5: FIG. 5. The first three stationary solutions, according to the [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]

**Figure 6.** Figure 6: FIG. 6. Approximations to the function [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]

**Figure 7.** Figure 7: depicts the value of ϕ (N) for N = 4, 5, and 10, at τ = 1. The parameters are such that ϕ (N) (0, x) = ϕ(0, x), and A1 = A2 = A4 = −A3 = 1, An = 0 for n > 4, with A˙ (N) n (0) = 0 for all n. This means that only the first four particles of the finite mechanical systems…

**Figure 10.** Figure 10: FIG. 10. Potential energy [PITH_FULL_IMAGE:figures/full_fig_p009_10.png]

**Figure 9.** Figure 9: FIG. 9. The local and total errors of several simulations. The [PITH_FULL_IMAGE:figures/full_fig_p009_9.png]

**Figure 11.** Figure 11: FIG. 11. Potential energy [PITH_FULL_IMAGE:figures/full_fig_p010_11.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references · 20 canonical work pages

[1]

Furthermore, if β ̸= 0, these particles interact with coupling constant modulated by the tensor Dnmpq. This is, in fact, the core idea behind the standard model of particle physics, in which fields become platforms to describe particle physics after suitable bases of eigenfunctions are fixed. It is in- structive to work with normalized variables by defini...

work page
[2]

P. W. Higgs, Broken Symmetries and the Masses of Gauge Bosons, Phys. Rev. Lett. 13, 508 (1964)

work page 1964
[3]

M. K. Gaillard, P. D. Grannis, and F. J. Sciulli, The standard model of particle physics, Rev. Mod. Phys. 71, S96 (1999)

work page 1999
[4]

sufficiently

(37) Also, by imposing that q2 must be real, we find that 1 2 < q2 4 < 1. (38) The rightmost inequality is necessary to guarantee that cn(u, q4) is real for all u [14]. The constants qi can be fixed by specifying the value of the parameter λ. For Dirichlet BC we find that cn(q3, q4) = cn (πq2 + q3, q4) = 0, (39) implying q3 = (2n′ + 1)K(q4), (40) q2 = |λ|...

work page 2023
[5]

A. M. Sirunyan, A. Tumasyan, W. Adam, F. Ambrogi, E. Asilar, T. Bergauer, J. Brandstetter, M. Dragicevic, J. Er¨ o,et al. (CMS Collaboration), Observation of Higgs Boson Decay to Bottom Quarks, Phys. Rev. Lett. 121, 121801 (2018)

work page 2018
[6]

X. Fan, T. G. Myers, B. A. D. Sukra, and G. Gabrielse, Measurement of the Electron Magnetic Moment, Phys. Rev. Lett. 130, 071801 (2023)

work page 2023
[7]

M. E. Peskin and D. V. Schroeder, An introduction to quantum field theory (Westview, Boulder, CO, 1995)

work page 1995
[8]

Kukuljan, S

I. Kukuljan, S. Sotiriadis, and G. Takacs, Correlation Functions of the Quantum Sine-Gordon Model in and out of Equilibrium, Phys. Rev. Lett. 121, 110402 (2018)

work page 2018

Show all 20 references

[9]

C. C. H. Ribeiro and U. R. Fischer, Nonlocal field the- ory of quasiparticle scattering in dipolar Bose-Einstein condensates, SciPost Phys. Core 6, 003 (2023)

2023
[10]

C. A. J. Fletcher, Computational Galerkin Methods , Springer Series in Computational Physics (Springer- Verlag, New York, 1984)

1984
[11]

B. G. Galerkin, Rods and plates : series in some ques- tions of elastic equilibrium of rods and plates (National Technical Information Service, Springfield, Va, 1968)

1968
[12]

E. N. Lorenz, Deterministic Nonperiodic Flow, Journal of Atmospheric Sciences 20, 130 (1963)

1963
[13]

M. A. Al-Gwaiz, Sturm-Liouville Theory and Its Applica- tions, 1st ed., Springer Undergraduate Mathematics Se- ries, Vol. 1 (Springer London, 2008)

2008
[14]

Z. X. Wang and D. R. Guo, Special Functions, 2nd ed. (World Scientific, Singapore, 2000)

2000
[15]

D. F. Lawden, Elliptic Functions and Applications , Ap- plied Mathematical Sciences, Vol. 80 (Springer-Verlag, New York, 1989)

1989
[16]

I. S. Gradshteyn, I. M. Ryzhik, D. Zwillinger, and V. Moll, Table of integrals, series, and products; 8th ed. (Academic Press, Amsterdam, 2015)

2015
[17]

Ates and M

E. Ates and M. Inc, Travelling wave solutions of general- ized Klein–Gordon equations using Jacobi elliptic func- tions, Nonlinear Dynamics 88, 2281 (2017)

2017
[18]

Paiva, Interaction of Dirac δ-waves in the nonlinear Klein-Gordon equation, Journal of Differential Equations 270, 1196 (2021)

A. Paiva, Interaction of Dirac δ-waves in the nonlinear Klein-Gordon equation, Journal of Differential Equations 270, 1196 (2021)

2021
[19]

Saltzman, Finite Amplitude Free Convection as an Ini- tial Value Problem—I, Journal of Atmospheric Sciences 19, 329 (1962)

B. Saltzman, Finite Amplitude Free Convection as an Ini- tial Value Problem—I, Journal of Atmospheric Sciences 19, 329 (1962)

1962
[20]

W. Bao, D. Jaksch, and P. A. Markowich, Numerical so- lution of the Gross–Pitaevskii equation for Bose–Einstein condensation, Journal of Computational Physics 187, 318 (2003)

2003

Pith tools

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

[1] [1]

Furthermore, if β ̸= 0, these particles interact with coupling constant modulated by the tensor Dnmpq. This is, in fact, the core idea behind the standard model of particle physics, in which fields become platforms to describe particle physics after suitable bases of eigenfunctions are fixed. It is in- structive to work with normalized variables by defini...

work page

[2] [2]

P. W. Higgs, Broken Symmetries and the Masses of Gauge Bosons, Phys. Rev. Lett. 13, 508 (1964)

work page 1964

[3] [3]

M. K. Gaillard, P. D. Grannis, and F. J. Sciulli, The standard model of particle physics, Rev. Mod. Phys. 71, S96 (1999)

work page 1999

[4] [4]

sufficiently

(37) Also, by imposing that q2 must be real, we find that 1 2 < q2 4 < 1. (38) The rightmost inequality is necessary to guarantee that cn(u, q4) is real for all u [14]. The constants qi can be fixed by specifying the value of the parameter λ. For Dirichlet BC we find that cn(q3, q4) = cn (πq2 + q3, q4) = 0, (39) implying q3 = (2n′ + 1)K(q4), (40) q2 = |λ|...

work page 2023

[5] [5]

A. M. Sirunyan, A. Tumasyan, W. Adam, F. Ambrogi, E. Asilar, T. Bergauer, J. Brandstetter, M. Dragicevic, J. Er¨ o,et al. (CMS Collaboration), Observation of Higgs Boson Decay to Bottom Quarks, Phys. Rev. Lett. 121, 121801 (2018)

work page 2018

[6] [6]

X. Fan, T. G. Myers, B. A. D. Sukra, and G. Gabrielse, Measurement of the Electron Magnetic Moment, Phys. Rev. Lett. 130, 071801 (2023)

work page 2023

[7] [7]

M. E. Peskin and D. V. Schroeder, An introduction to quantum field theory (Westview, Boulder, CO, 1995)

work page 1995

[8] [8]

Kukuljan, S

I. Kukuljan, S. Sotiriadis, and G. Takacs, Correlation Functions of the Quantum Sine-Gordon Model in and out of Equilibrium, Phys. Rev. Lett. 121, 110402 (2018)

work page 2018