Pith. sign in

REVIEW 1 major objections 2 minor 23 references

Translating the fast time in Van der Pol's equation produces an exponentially small phase shift that varies with the slow time and arises from nonlinearity.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · grok-4.3

2026-06-26 02:42 UTC pith:3VAYBCTM

load-bearing objection The key thing is an exponentially small nonlinearity-induced phase shift that breaks scale independence in multiple scales for Van der Pol, found via high-order expansion and optimal truncation. the 1 major comments →

arxiv 2606.27038 v1 pith:3VAYBCTM submitted 2026-06-25 math.DS nlin.PS

On the independence of the slow and fast scales in multiple-scale expansions, with application to Van der Pol's equation

Gregory Kozyreff , John R. King This is my paper
classification math.DS nlin.PS
keywords multiple scalesVan der Pol equationasymptotic expansionsexponentially small termsphase shiftslow-fast dynamicsnonlinear oscillatorsoptimal truncation
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper questions the standard assumption in multiple-scale perturbation methods that slow and fast time scales can be treated as independent variables. It applies the method to Van der Pol's equation, a model weakly nonlinear oscillator, and extends the expansion to arbitrarily high orders. Because the series diverges, optimal truncation is used to extract the effect of an initial constant shift in the fast coordinate. This procedure isolates a phase correction that is exponentially small in the small parameter yet depends on the slow coordinate. The correction is traced directly to the nonlinear term and is confirmed by direct numerical integration of the original equation.

Core claim

When the multiple-scale expansion of Van der Pol's equation is carried to high order and optimally truncated, an initial translation of the fast time coordinate generates a phase shift that is exponentially small in the perturbation parameter and varies with the slow time; the shift originates in the nonlinearity and therefore demonstrates that the slow and fast scales are not independent.

What carries the argument

Optimal truncation of the divergent asymptotic series obtained from the multiple-scale expansion, which isolates the exponentially small, slow-dependent phase correction induced by the cubic nonlinearity.

Load-bearing premise

That carrying the multiple-scale expansion to arbitrarily high order and then applying optimal truncation reliably detects and measures the exponentially small phase shift produced by the nonlinearity.

What would settle it

A direct numerical solution of Van der Pol's equation in which an imposed constant shift of the fast time fails to produce a phase difference whose magnitude scales exponentially with the small parameter and whose value changes with slow time.

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

If this is right

  • The formal independence of scales in multiple-scale analysis holds only up to algebraically small orders and is violated by exponentially small terms generated by nonlinearity.
  • The size of the induced phase shift increases when the nonlinearity is quadratic rather than cubic, suggesting a stronger coupling in other Hopf bifurcations.
  • Long-time predictions obtained from truncated multiple-scale expansions acquire an additional error that grows with the slow time through this phase correction.
  • The same optimal-truncation procedure can be used to quantify exponentially small corrections in other weakly nonlinear oscillators.
  • The result supplies an explicit mechanism by which nonlinearity couples scales that are formally separated.

Where Pith is reading between the lines

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

  • The same exponentially small coupling may appear in multiple-scale analyses of pattern-forming systems or slowly varying waves whenever nonlinearity is present.
  • Numerical integrators for stiff oscillatory problems could incorporate an analogous phase adjustment to reduce long-term drift.
  • The finding raises the possibility that scale separation assumptions in averaging methods for systems with more than two time scales also require exponentially small corrections.
  • Extensions to forced or damped oscillators would test whether the phase shift persists or is modified by additional terms.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

1 major / 2 minor

Summary. The paper examines the formal independence of slow and fast time scales in the method of multiple scales, using Van der Pol's equation as a model weakly nonlinear oscillator and Hopf bifurcation. It performs the expansion to arbitrarily high order, notes the resulting divergent series, and applies optimal truncation to identify an exponentially small phase shift in the fast coordinate that depends on the slow time and originates from the cubic nonlinearity; this shift is said to reconnect the scales. Numerical simulations are invoked to confirm the existence and scaling of the shift, with indications that the phenomenon extends to quadratic nonlinearities.

Significance. If the central claim holds, the work would demonstrate that exponentially small terms arising from nonlinearity can break the assumed independence of scales in multiple-scale expansions, offering a refined view of the method's accuracy for oscillators and suggesting broader implications for asymptotic treatments of Hopf bifurcations.

major comments (1)
  1. [Abstract] Abstract: the central claim that optimal truncation isolates a unique nonlinearity-induced, slow-time-dependent phase shift heta( au) rests on the unverified assumption that the remainder after truncation is dominated by the Stokes phenomenon from the cubic term and is free of contamination by other exponentially small contributions or residual secular terms; no explicit remainder estimates, comparison to Borel summation, or hyperasymptotic analysis is supplied to support this.
minor comments (2)
  1. [Abstract] Abstract: 'one cannot failt to notice' contains a typographical error.
  2. [Abstract] Abstract: the statement that 'the calculation is carried out in sufficient detail to provide confidence in the generality' would benefit from explicit cross-references to the sections containing the high-order terms and truncation procedure.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thoughtful review of our manuscript. We address the major comment point by point below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim that optimal truncation isolates a unique nonlinearity-induced, slow-time-dependent phase shift heta( au) rests on the unverified assumption that the remainder after truncation is dominated by the Stokes phenomenon from the cubic term and is free of contamination by other exponentially small contributions or residual secular terms; no explicit remainder estimates, comparison to Borel summation, or hyperasymptotic analysis is supplied to support this.

    Authors: We acknowledge that the manuscript applies optimal truncation in the standard manner without supplying explicit remainder estimates, Borel summation, or hyperasymptotic analysis to rigorously isolate the Stokes contribution from the cubic term. The central claim is supported by the explicit high-order calculation, the identification of the phase shift, and direct numerical confirmation of its existence and scaling. In the revised manuscript we will insert a concise paragraph in the discussion section that states the working assumptions regarding the remainder and notes the absence of secular terms at the optimal truncation point. This addition will clarify the scope of the claim while preserving the formal and numerical results already presented. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper derives the exponentially small phase shift via direct multiple-scale expansion of Van der Pol's equation to arbitrary order followed by optimal truncation. This procedure operates on the original ODE without any fitted parameters renamed as predictions, self-definitional loops, or load-bearing self-citations. The result follows from the asymptotic construction applied to the differential equation itself and is therefore self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based solely on abstract; no explicit free parameters or invented entities described. The analysis relies on standard assumptions of asymptotic perturbation theory.

axioms (1)
  • domain assumption Multiple-scale expansions can be carried to arbitrarily large order while remaining amenable to optimal truncation.
    Invoked to reconnect the scales and extract the phase shift.

pith-pipeline@v0.9.1-grok · 5789 in / 1085 out tokens · 43696 ms · 2026-06-26T02:42:21.612738+00:00 · methodology

0 comments
read the original abstract

When implementing the method of multiple scales, one is traditionally instructed to treat the slow and fast time scales as if they were independent. Despite the intuitive motivation and the effectiveness of this perturbation method, one cannot failt to notice that these two scales relate to the same unique variable, so independence can only be formal. How sensible is it, then, to split a variable asymptotically into two (or more) independent ones? In this paper, we elucidate this issue with Van der Pol's equation, one of the simplest weakly nonlinear oscillators, as well as a simple example of a Hopf bifurcation. The discussion involves carrying the multiple-scale analysis up to arbitrarily large order and dealing with the divergent character of the resulting asymptotic series. Using the technique of optimal truncation, we re-connect the two scales. Specifically, we show that an initial translation of the fast coordinate leads to a non-trivial, exponentially small, phase shift that depends on the slow coordinate. This phase shift breaks the independence of the slow and fast scales and is found to result from the nonlinearity. Numerical simulations confirm its existence, as well as the predicted scaling. The calculation is carried out in sufficient detail to provide confidence in the generality of our result, both in its essence and in its form. In particular, we find strong indications that a Hopf bifurcation with a quadratic nonlinearity would lead to the same phenomenon, but with a larger magnitude.

Figures

Figures reproduced from arXiv: 2606.27038 by Gregory Kozyreff, John R. King.

Figure 1
Figure 1. Figure 1: Typical two-scale behaviours resulting from differential equations: (a) a weakly amplified linear oscillator, governed by Y ′′(t) − ϵY ′ (t) + Y (t) = 0 and (b) Van der Pol’s equation, Y ′′(t) − ϵ[1 − Y (t) 2 ]Y ′ (t) + Y (t) = 0 with ϵ = 0.05. The dashed lines are the slowly evolving amplitudes of oscillation, as obtained by a routine multiple-scale analysis [1]. 1. Introduction Many autonomous differenti… view at source ↗
Figure 2
Figure 2. Figure 2: Families of late terms and their Stokes lines in the “slow” complex plane. remains relevant to our consideration is ∆k = 2, corresponding to k = 1, 3. We thus have, for large order j, and taking only the singularity τ0 into account: yj ∼ X k=1,3 Aj,ke ik(θ+ϕ) ∼ X k=1,3  i 2 (τ0 − τ ) j+α Γ(j + α) h fk(τ ) + O  j −1 i e ik(θ+ϕ) , (3.6) Using the above expression, we deduce that Cj,k ∼  i 2 (τ0 − τ ) j… view at source ↗
Figure 3
Figure 3. Figure 3: Numerically computed phase shift of the oscillations after the establishment of the nonlinear oscillations, as a function of the initial phase, ϕ, see (4.1), for three values of ϵ. If the evolution of the complex amplitude of oscillations was independent of ϕ, this phase shift should be constant. However it oscillates with ϕ in agreement with (3.32). We define on these graphs ∆(ϵ) as the maximum variation … view at source ↗
Figure 4
Figure 4. Figure 4: Comparison between the analytical prediction (3.33) (dashed line) and values drawn from [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Convergence and accelerated convergence of the series Λ (0) j towards the value Λ. we find it much more convenient to successively compute Yj = 2 X j+1 k=1  Bj,kz k + Bj,−kz −k  ,  Y 2  j = X j j ′=0 Yj ′Yj−j ′ ,  Y 3  j = X j j ′=0  Y 2  j ′ Yj−j ′ , (A 6) and extract Dj,k as the coefficient of z k in  Y 3  j . Such manipulations can easily be done with a computer algebra system such as Mathemat… view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

23 extracted references · 13 canonical work pages · 1 internal anchor

  1. [1]

    (to appear in 2026)Practical asymptotics

    Kozyreff G. (to appear in 2026)Practical asymptotics. Cambridge University Press

    2026

  2. [2]

    1999Advanced Mathematical Methods for Scientists and Engineers

    Bender CM, Orszag SA. 1999Advanced Mathematical Methods for Scientists and Engineers. Springer-Verlag

  3. [3]

    1991Perturbation methods

    Hinch EJ. 1991Perturbation methods. Cambridge University Press

  4. [4]

    2023 Speed of wave packets and the nonlinear Schrödinger equation.Phys

    Kozyreff G. 2023 Speed of wave packets and the nonlinear Schrödinger equation.Phys. Rev. E

    2023

  5. [5]

    (10.1103/physreve.107.014219)

    work page doi:10.1103/physreve.107.014219

  6. [6]

    1997 On asymmetric gravity-capillary solitary waves.J

    Yang TS, Akylas TR. 1997 On asymmetric gravity-capillary solitary waves.J. Fluid Mech.330, 215

    1997

  7. [7]

    1999 Effects of exponentially small terms in the perturbation approach to localized buckling.Proc

    Wadee MK, Bassom AP . 1999 Effects of exponentially small terms in the perturbation approach to localized buckling.Proc. R. Soc. A455, 2351. (10.1098/rspa.1999.0407) 13rspa.royalsocietypublishing.org Proc R Soc A 0000000

    work page doi:10.1098/rspa.1999.0407 1999

  8. [8]

    2017 New Aspects of Exponential Asymptotics in Multiple-Scale Nonlinear Wave Problems.Studies in Applied Mathematics139, 223–247

    Nixon SD, Akylas TR, Yang J. 2017 New Aspects of Exponential Asymptotics in Multiple-Scale Nonlinear Wave Problems.Studies in Applied Mathematics139, 223–247. (10.1111/sapm.12179)

    work page doi:10.1111/sapm.12179 2017

  9. [9]

    2003 Beyond-all-order effects in multiple-scales asymptotics: travelling-wave solutions to the Kuramoto-Sivashinsky equation.J

    Adams KL, King JR, Tew RH. 2003 Beyond-all-order effects in multiple-scales asymptotics: travelling-wave solutions to the Kuramoto-Sivashinsky equation.J. Eng. Math.45, 197

    2003

  10. [10]

    2006 Asymptotics of Large Bound States of Localized Structures

    Kozyreff G, Chapman SJ. 2006 Asymptotics of Large Bound States of Localized Structures. Phys. Rev. Lett.97, 044502

    2006

  11. [11]

    2009 Exponential asymptotics of localised patterns and snaking bifurcation diagrams.Physica D238, 319

    Chapman SJ, Kozyreff G. 2009 Exponential asymptotics of localised patterns and snaking bifurcation diagrams.Physica D238, 319. (10.1016/j.physd.2008.10.005)

    work page doi:10.1016/j.physd.2008.10.005 2009

  12. [12]

    2009 Localized Beating between Dynamically Generated Frequencies.Phys

    Kozyreff G, Tlidi M, Mussot A, Louvergneaux E, Taki M, Vladimirov AG. 2009 Localized Beating between Dynamically Generated Frequencies.Phys. Rev. Lett.102, 043905. (10.1103/PhysRevLett.102.043905)

    work page doi:10.1103/physrevlett.102.043905 2009

  13. [13]

    2011 Exponential asymptotics of homoclinic snaking.Nonlinearity24, 3323

    Dean AD, Matthews PC, Cox SM, King JR. 2011 Exponential asymptotics of homoclinic snaking.Nonlinearity24, 3323. (10.1088/0951-7715/24/12/003)

    work page doi:10.1088/0951-7715/24/12/003 2011

  14. [14]

    2015 Orientation-Dependent Pinning and Homoclinic Snaking on a Planar Lattice.SIAM Journal on Applied Dynamical Systems14, 481–521

    Dean AD, Matthews PC, Cox SM, King JR. 2015 Orientation-Dependent Pinning and Homoclinic Snaking on a Planar Lattice.SIAM Journal on Applied Dynamical Systems14, 481–521. (10.1137/140966897)

    work page doi:10.1137/140966897 2015

  15. [15]

    2019 Beyond all order asymptotics for homoclinic snaking in a Schnakenberg system.Nonlinearity32, 2667–2693

    de Witt H. 2019 Beyond all order asymptotics for homoclinic snaking in a Schnakenberg system.Nonlinearity32, 2667–2693. (10.1088/1361-6544/ab0b1d)

    work page doi:10.1088/1361-6544/ab0b1d 2019

  16. [16]

    https://doi.org/10.48550/arXiv

    Villar-Sepúlveda E. 2025 Beyond-all-order asymptotics for homoclinic snaking of localised patterns in reaction-transport systems. (10.48550/ARXIV .2501.02698)

    work page internal anchor Pith review doi:10.48550/arxiv 2025

  17. [17]

    1979 Exact coefficients of the limit cycle in van der Pol’s equation

    Deprit A, Schmidt DS. 1979 Exact coefficients of the limit cycle in van der Pol’s equation. Journal of Research of the National Bureau of standards84, 293

    1979

  18. [18]

    1982 Power Series Expansions for the Frequency and Period of the Limit Cycle of the Van Der Pol Equation.SIAM Journal on Applied Mathematics42, 678–693

    Andersen CM, Geer JF. 1982 Power Series Expansions for the Frequency and Period of the Limit Cycle of the Van Der Pol Equation.SIAM Journal on Applied Mathematics42, 678–693. (10.1137/0142047)

    work page doi:10.1137/0142047 1982

  19. [19]

    1998 The Periodic Solution of Van Der Pol’s Equation.SIAM Journal on Applied Mathematics59, 156–171

    Buonomo A. 1998 The Periodic Solution of Van Der Pol’s Equation.SIAM Journal on Applied Mathematics59, 156–171. (10.1137/s0036139997319797)

    work page doi:10.1137/s0036139997319797 1998

  20. [20]

    2018 High order analysis of the limit cycle of the van der Pol oscillator.Journal of Mathematical Physics59

    Amore P , Boyd JP , Fernandez FM. 2018 High order analysis of the limit cycle of the van der Pol oscillator.Journal of Mathematical Physics59. (10.1063/1.5016961)

    work page doi:10.1063/1.5016961 2018

  21. [21]

    1973Asymptotic expansions: their derivation and interpretation

    Dingle R. 1973Asymptotic expansions: their derivation and interpretation. Academic Press Inc

  22. [22]

    1995 Stokes Phenomenon and Matched Asymptotic Expansions.SIAM Journal on Applied Mathematics55, 1469–1483

    Olde Daalhuis AB, Chapman SJ, King JR, Ockendon JR, Tew RH. 1995 Stokes Phenomenon and Matched Asymptotic Expansions.SIAM Journal on Applied Mathematics55, 1469–1483. (10.1137/s0036139994261769)

    work page doi:10.1137/s0036139994261769 1995

  23. [23]

    1989 Uniform asymptotic smoothing of Stokes’s discontinuities.Proc

    Berry MV . 1989 Uniform asymptotic smoothing of Stokes’s discontinuities.Proc. R. Soc. A422, 7–21

    1989