pith. sign in

arxiv: 1907.10875 · v1 · pith:XYTT32PHnew · submitted 2019-07-25 · 🧮 math.PR · math-ph· math.MP

Stochastic perturbation of a cubic anharmonic oscillator

Enrico Bernardi , Alberto Lanconelli This is my paper

Pith reviewed 2026-05-24 16:32 UTC · model grok-4.3

classification 🧮 math.PR math-phmath.MP
keywords stochastic perturbationcubic anharmonic oscillatorJacobi elliptic functionsLamé equationwhite noiseperiodic solutionprobability boundperturbation expansion
0
0 comments X

The pith

The stochastically perturbed cubic anharmonic oscillator admits an explicit power-series solution in the noise strength whose coefficients solve perturbed Lamé equations via Jacobi elliptic functions, and the approximation stays close to a

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

The paper studies a cubic anharmonic oscillator driven by small additive white noise. It begins with initial conditions chosen so the noise-free equation has a periodic orbit. The solution is expanded formally as a power series in the diffusion coefficient. Each coefficient in the series satisfies a stochastic version of Lamé's equation, solved explicitly with Jacobi elliptic functions. A lower bound is derived on the probability that a finite truncation of the series remains near the deterministic periodic solution over a fixed interval, together with conditions under which the series converges.

Core claim

We perturb the Hamiltonian system of the cubic anharmonic oscillator with additive Gaussian white noise starting from initial conditions yielding a periodic solution. A formal expansion in powers of the diffusion coefficient is written, and the coefficients are shown to be unique strong solutions to stochastic perturbations of Lamé's equation. Explicit solutions are obtained in terms of Jacobi elliptic functions. A lower bound is proved for the probability that an approximation of the stochastic solution stays close to the deterministic solution, along with conditions for convergence of the expansion.

What carries the argument

Stochastic perturbations of Lamé's equation whose unique strong solutions are expressed via Jacobi elliptic functions, used to control the probability that the truncated series stays close to the deterministic periodic orbit.

If this is right

  • Under the stated convergence conditions the power series yields an approximate solution whose deviation from the deterministic orbit has probability bounded away from zero.
  • The explicit Jacobi elliptic expressions allow direct calculation of moments and other statistics for each term in the expansion.
  • The same expansion technique applies to any Hamiltonian system whose unperturbed periodic solutions are expressible by elliptic functions.

Where Pith is reading between the lines

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

  • The positive lower bound suggests that for sufficiently small noise the oscillator retains its periodic character over observable time scales, which could be checked in physical experiments with controlled random forcing.
  • The method supplies a template for constructing explicit perturbation series around other integrable nonlinear oscillators that admit elliptic-function solutions.
  • Numerical verification of the closeness probability can be performed by integrating the original SDE and comparing trajectories to the deterministic orbit, providing a direct test independent of the analytic bound.

Load-bearing premise

The initial conditions are chosen so that the noise-free equation possesses a periodic solution.

What would settle it

A Monte Carlo simulation of the stochastic differential equation that produces an empirical closeness probability strictly below the claimed lower bound for any truncation order and sufficiently small noise strength.

Figures

Figures reproduced from arXiv: 1907.10875 by Alberto Lanconelli, Enrico Bernardi.

Figure 1
Figure 1. Figure 1: Graph of Hamiltonian with c = −1, a = 1 Then, slightly changing our notations, we start directly with the three real roots of the polynomial and we denote them by c, −a − c, a with c < 0 < a; imposing x 3 /3 − Bx + H(y, η) = (x − a)(x + a + c)(x − c)/3 we get B = (a 2 + c 2 + ac)/3 and H(y, η) = ac(a + c)/3. Therefore, the Hamiltonian system we are going to analyze is ( x˙ 0(t) = ξ0(t), x0(0) = y ˙ξ0(t) = … view at source ↗
Figure 2
Figure 2. Figure 2: Graph of (2.7) with c = −1, a = 1 Here the oval part in [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Graph of H(y, η) = 0 2.2 The equation for x1: a stochastic Lam´e’s equation We now want to solve x¨1(t) = 2x0(t)x1(t) + B˙(t), x1(0) = 0 ˙x1(0) = 0 (2.9) which is equivalent to the system of stochastic differential equations ( dx1(t) = ξ1(t)dt, x1(0) = 0 dξ1(t) = 2x0(t)x1(t)dt + dB(t), ξ1(0) = 0 7 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Graph of u1 with q = 2/ √ 5 with C(t, q) : = −dn(t, q)  −1 + q 2 +  2 + q 2 (−5 + (3 − 2q 2 )q 2 )  cn2 (t, q) +2q 2  1 + (−1 + q 2 )q 2  cn4 (t, q) +  (2 − q 2 )(−1 + q 2 )x + 2(q 4 − q 2 + 1)E(t, q)  ×sn(t, q)cn(t, q)dn(t, q)  and D(t, q) := (−1 + q 2 ) 2 sn(t, q)cn(t, q)dn2 (t, q) we get that u2(t) = C(t, q) + D(t, q)u1(t, q) (2.15) is a second independent solution of (2.11). Here C(t, q) := α0(… view at source ↗
Figure 5
Figure 5. Figure 5: Graph of u2 with q = 2/ √ 5 with E(q) denotes the complete elliptic integral of the second kind. This coefficient behaves like this when 0 < q < 1: 0.2 0.4 0.6 0.8 1.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Graph of µ(q) with q ∈ (0, 1) Going back to equation (2.10) we have found two linearly independent solutions: w1(t) = u1 r a − c 6 t + iy ! and w2(t) = u2 r a − c 6 t + iy ! . It is easy to verify that the Wronskian determinant of (u1, u2) is −(1 − q 2 ) 2 6= 0 and that of (w1, w2) is − qa−c 6 (1 − q 2 ) 2 Without loss of generality multiplying w1 and w2 by suitable constants we may assume that their Wrons… view at source ↗
Figure 7
Figure 7. Figure 7: Graph of three paths of process (2.17) 12 [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
read the original abstract

We perturb with an additive Gaussian white noise the Hamiltonian system associated to a cubic anharmonic oscillator. The stochastic system is assumed to start from initial conditions that guarantee the existence of a periodic solution for the unperturbed equation. We write a formal expansion in powers of the diffusion parameter for the candidate solution and analyze the probabilistic properties of the sequence of the coefficients. It turns out that such coefficients are the unique strong solutions of stochastic perturbations of the famous Lam\'e's equation. We obtain explicit solutions in terms of Jacobi elliptic functions and prove a lower bound for the probability that an approximated version of the solution of the stochastic system stay close to the solution of the deterministic problem. Conditions for the convergence of the expansion are also provided.

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 / 1 minor

Summary. The paper perturbs the Hamiltonian system of a cubic anharmonic oscillator by additive Gaussian white noise. Starting from initial conditions yielding a periodic orbit in the deterministic case, it constructs a formal power-series expansion in the diffusion coefficient. The coefficients are identified as unique strong solutions to stochastic perturbations of Lamé's equation; explicit solutions in Jacobi elliptic functions are claimed for these coefficients. A lower bound is proved on the probability that a truncated version of the expansion remains close to the deterministic solution, together with conditions for convergence of the series.

Significance. If the pathwise closed-form expressions hold, the work supplies an explicit, verifiable probabilistic estimate of closeness between stochastic and deterministic trajectories for a nonlinear oscillator, which is rare. The reduction to perturbed Lamé equations and the use of elliptic functions give a concrete handle on the expansion that could support further quantitative results in stochastic perturbation theory.

major comments (1)
  1. [derivation of explicit solutions for stochastic coefficients] The section deriving the explicit solutions (following the identification of the coefficients as strong solutions of the stochastic Lamé perturbations): the deterministic Lamé equation is solved by Jacobi elliptic functions, but the SDE version includes stochastic integrals. The manuscript must show explicitly how these integrals are absorbed or cancel so that the closed-form expression holds almost surely (not merely formally or in expectation). This step is load-bearing for the subsequent probability lower bound, which is stated for the approximation constructed from these explicit coefficients.
minor comments (1)
  1. [Introduction] Clarify the precise statement of the initial conditions in the introduction; the current phrasing leaves open whether they are deterministic or random.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and for identifying the need for greater explicitness in the derivation of the closed-form solutions. We address the single major comment below.

read point-by-point responses

  1. Referee: The section deriving the explicit solutions (following the identification of the coefficients as strong solutions of the stochastic Lamé perturbations): the deterministic Lamé equation is solved by Jacobi elliptic functions, but the SDE version includes stochastic integrals. The manuscript must show explicitly how these integrals are absorbed or cancel so that the closed-form expression holds almost surely (not merely formally or in expectation). This step is load-bearing for the subsequent probability lower bound, which is stated for the approximation constructed from these explicit coefficients.

    Authors: We agree that the manuscript would benefit from a more detailed, step-by-step verification that the stochastic integrals cancel almost surely, so that the Jacobi-elliptic closed forms hold pathwise. In the revised version we will insert an expanded derivation (new subsection after the identification of the coefficients as strong solutions) that substitutes the candidate elliptic expressions into the SDE, computes the resulting stochastic integrals explicitly via the variation-of-parameters formula adapted to the Lamé operator, and verifies their cancellation using the specific form of the additive noise and the Wronskian identity for the elliptic fundamental solutions. This will be done before invoking the probability bound, thereby removing any ambiguity about whether the expressions are merely formal. revision: yes

Circularity Check

0 steps flagged

No circularity: power-series ansatz yields independent linear SDEs whose solutions are analyzed separately from the target probability bound.

full rationale

The derivation begins from the deterministic periodic orbit of the cubic oscillator, introduces a formal power-series expansion in the diffusion parameter, and obtains a sequence of linear stochastic equations for the coefficients (stochastic perturbations of Lamé's equation). The paper then claims explicit Jacobi-elliptic representations for those coefficients and uses the resulting approximation to prove a lower bound on the probability of closeness to the deterministic solution. No equation or step equates the target probability bound to a fitted quantity defined from the same data, nor does any load-bearing premise reduce to a self-citation or ansatz smuggled from prior work by the same authors. The stochastic analysis and convergence conditions are presented as independent consequences of the linear SDE structure, making the chain self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review; the only explicit premise is the existence of the deterministic periodic orbit under chosen initial data. No free parameters or invented entities are mentioned.

axioms (1)
  • domain assumption Existence of a periodic solution for the unperturbed cubic anharmonic oscillator under the chosen initial conditions.
    Stated directly in the abstract as the starting point for the stochastic perturbation.

pith-pipeline@v0.9.0 · 5641 in / 1358 out tokens · 23369 ms · 2026-05-24T16:32:18.221281+00:00 · methodology

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

22 extracted references · 22 canonical work pages

  1. [1]

    Albeverio, A

    S. Albeverio, A. Hilbert and E. Zehnder, Hamiltonian systems with a stochastic force: nonlinear versus linear, and a Girsanov formula, Stochastics Stochastics Rep. 39 2-3 (1992) 159-188

    work page 1992

  2. [2]

    Albeverio, A

    S. Albeverio, A. Hilbert and V. Kolokoltsov, Estimates uniform in time for the transition probability of diffusions with small drift and for stochastically perturbed Newton equations, J. Theoret. Probab. 12 2 (1999) 293-300

    work page 1999

  3. [3]

    Arscott, Periodic Differential Equations, The Macmillan Company, New York, 1964

    F.M. Arscott, Periodic Differential Equations, The Macmillan Company, New York, 1964

    work page 1964

  4. [4]

    Arscott and I.M.Khabaza, Table of Lam´ e’s Polynomials, Pergamon Press, Oxford, London, New York, Paris, 1962

    F.M. Arscott and I.M.Khabaza, Table of Lam´ e’s Polynomials, Pergamon Press, Oxford, London, New York, Paris, 1962

    work page 1962

  5. [5]

    Bernardi and T

    E. Bernardi and T. Nishitani, On the Cauchy problem for non-effectively hyperbolic operators, the Gevrey 5 well-posedness, J. d’Analyse Math. 105 (2008) 197–240

    work page 2008

  6. [6]

    Bernardi and T

    E. Bernardi and T. Nishitani, On the Cauchy problem for non-effectively hyperbolic operators, the Gevrey 4 well-posedness, Kyoto J. Math. 51 (2011) 767–810

    work page 2011

  7. [7]

    Delabaere and D.T

    E. Delabaere and D.T. Trinh, Spectral analysis of the complex cubic oscillator, J. Phys. A: Math.gen. 33 2 (2000) 8771-8796

    work page 2000

  8. [8]

    Ferreira and J

    E.M. Ferreira and J. Sesma, Global solution of the cubic oscillator, J. of Phys. A: Math. 47 (2014) 415306

    work page 2014

  9. [9]

    Gardiner, Handbook of Stochastic Methods , Third Edition, Springer-Verlag, New York, 2004

    C.W. Gardiner, Handbook of Stochastic Methods , Third Edition, Springer-Verlag, New York, 2004

    work page 2004

  10. [10]

    Gradshteyn and I.M

    I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products, Elsevier, 7th Edition, 2007

    work page 2007

  11. [11]

    H¨ ormander, The Cauchy problem for differential equations with double charac- teristics, Journal D’Analyse Math´ ematique, 32, (1977), 118-196

    L. H¨ ormander, The Cauchy problem for differential equations with double charac- teristics, Journal D’Analyse Math´ ematique, 32, (1977), 118-196

    work page 1977

  12. [12]

    Ikeda and S.Watanabe, Stochastic Differential Equations and Diffusion Pro- cesses, North Holland, Amsterdam, New York, Oxford, Kodansha, 1981

    N. Ikeda and S.Watanabe, Stochastic Differential Equations and Diffusion Pro- cesses, North Holland, Amsterdam, New York, Oxford, Kodansha, 1981

    work page 1981

  13. [13]

    Karatzas and S

    I. Karatzas and S. E. Shreve, Brownian motion and stochastic calculus , Springer- Verlag, New York, 1991 21

    work page 1991

  14. [14]

    Markus and A

    L. Markus and A. Weerasinghe, Stochastic oscillators, J. Differential Equations 71 2 (1988) 288-314

    work page 1988

  15. [15]

    Weerasinghe, Stochastic nonlinear oscillators, Adv

    L Markus and A. Weerasinghe, Stochastic nonlinear oscillators, Adv. in Appl. Probab. 25 3 (1993) 649-666

    work page 1993

  16. [16]

    T.Nishitani, A simple proof of the existence of tangent bicharacteristics for nonef- fectively hyperbolic operators Kyoto J. Math.. 55 (2015) 2 281–297

    work page 2015

  17. [17]

    Frank W.J.Olver and Daniel W.Lozier, NIST Handbook of Mathematical Functions, Cambridge University Press, 2010

    work page 2010

  18. [18]

    Revuz, M

    D. Revuz, M. Yor, Continuous Martingales and Brownian Motion , Third Edition, Springer-Verlag, Berlin, 1999

    work page 1999

  19. [19]

    Volker, Four Remarks on Eigenvalues of Lam´ e’s Equation, Analysis and Appli- cations 2 (2004) 161-175

    H. Volker, Four Remarks on Eigenvalues of Lam´ e’s Equation, Analysis and Appli- cations 2 (2004) 161-175

    work page 2004

  20. [20]

    Yakubovich and V.M

    V.A. Yakubovich and V.M. Starzhinskii, Linear Differential Equations With Peri- odic Coefficients vol.1, John Wiley & Sons, New York, 1975

    work page 1975

  21. [21]

    http://dlmf.nist.gov/22

    work page

  22. [22]

    http://dlmf.nist.gov/29 22

    work page