pith. sign in

arxiv: 2603.02257 · v1 · submitted 2026-02-27 · 🪐 quant-ph

Rayleigh-Ritz Variational Method in The Complex Plane

M.W. AlMasri This is my paper

Pith reviewed 2026-05-15 17:54 UTC · model grok-4.3

classification 🪐 quant-ph
keywords Rayleigh-Ritz variational methodSegal-Bargmann spaceGaussian trial functionsnormalizability conditionanharmonic oscillatorvariational upper boundsquantum oscillatorsdisplaced Gaussians
0
0 comments X

The pith

Rayleigh-Ritz variational method in Segal-Bargmann space recovers exact harmonic oscillator energies using generalized Gaussians that satisfy |α| < 1/2.

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

The paper establishes that the Rayleigh-Ritz method can be applied rigorously in the complex Segal-Bargmann representation of quantum oscillators. It derives the condition |α| < 1/2 as the requirement for generalized Gaussian trial functions of the form e^{α z² + β z} to have finite norm through direct analysis of the integrals. When the trial family includes the true ground state, the variational procedure returns the exact energy for the harmonic oscillator. For the quartic anharmonic oscillator the same ansatz produces a cubic stationarity condition whose solutions generate perturbative energy series that incorporate wavefunction narrowing. Monomial trials supply simple rigorous upper bounds but cannot adjust their width, while displaced versions handle parity breaking in asymmetric potentials.

Core claim

We present a systematic study of the Rayleigh-Ritz variational method for quantum oscillators in the Segal-Bargmann space. We rigorously derive the normalizability condition |α| < 1/2 for generalized Gaussian trial functions ψ(z) = e^{α z² + β z} through convergence analysis of Gaussian integrals in the complex plane. Applications to the harmonic oscillator demonstrate exact recovery of the ground state in Segal-Bargmann space when the trial family contains the true solution. For the quartic anharmonic oscillator adaptive Gaussian ansätze in position space yield a cubic stationarity equation and perturbative energy expansions beyond first order, capturing anharmonic wavefunction narrowing.

What carries the argument

The normalizability condition |α| < 1/2 obtained from convergence analysis of the integrals ∫ |e^{α z² + β z}|^2 dμ(z) in the complex plane, which makes the generalized Gaussian trial functions admissible for the Rayleigh-Ritz procedure.

If this is right

  • Exact recovery of the ground-state energy for the harmonic oscillator when the trial family contains the exact solution.
  • Cubic stationarity equation for the quartic oscillator that produces energy expansions capturing wavefunction narrowing beyond first order.
  • Rigorous upper bounds E_n = n + 1/2 + (3λ/4)(2n² + 2n + 1) for excited states using monomial trials.
  • Displacement parameters become necessary to capture parity breaking and stabilization in asymmetric potentials such as x³ + x⁴.

Where Pith is reading between the lines

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

  • The same normalizability analysis could be applied to trial functions with higher powers of z in the exponent for more flexible approximations.
  • Numerical minimization of the variational energy with respect to α and β may extend the method beyond the perturbative regime for stronger anharmonicity.
  • The complex-plane formulation may simplify calculations for potentials that admit analytic continuation into the complex domain.

Load-bearing premise

The Gaussian integrals over the complex plane converge precisely when |α| < 1/2, which is required for the trial functions to be square-integrable and for the variational upper bounds to be valid.

What would settle it

A direct numerical or analytic evaluation of the norm integral diverging for any |α| > 1/2, or a variational energy for the harmonic oscillator that exceeds the known exact value 1/2 when the trial family includes the true Gaussian.

read the original abstract

We present a systematic study of the Rayleigh--Ritz variational method for quantum oscillators in the Segal--Bargmann space. We rigorously derive the normalizability condition $|\alpha| < \tfrac{1}{2}$ for generalized Gaussian trial functions $\psi(z) = e^{\alpha z^2 + \beta z}$ through convergence analysis of Gaussian integrals in the complex plane. Applications to the harmonic oscillator demonstrate exact recovery of the ground state in Segal--Bargmann space when the trial family contains the true solution. For the quartic anharmonic oscillator ($\hat{H} = -\tfrac{1}{2}\partial_x^2 + \tfrac{1}{2}x^2 + \lambda x^4$), adaptive Gaussian ans\"atze in position space yield a cubic stationarity equation and perturbative energy expansions beyond first order, capturing anharmonic wavefunction narrowing. In contrast, monomial trial functions ($\psi_n(z) = z^n$) in the Segal--Bargmann space -- while providing rigorous upper bounds $E_n = n + \tfrac{1}{2} + \tfrac{3\lambda}{4}(2n^2 + 2n + 1)$ for excited states -- lack width adaptability and are limited to first-order accuracy for ground-state calculations. We further analyze displaced Gaussians and displaced monomials for asymmetric potentials (e.g., $x^3 + x^4$), showing that displacement parameters are essential to capture parity breaking and stabilization effects ($E_0 \approx \tfrac{1}{2} + \tfrac{3\mu}{4} - \tfrac{9\lambda^2}{4} + \cdots$).

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

0 major / 3 minor

Summary. The manuscript applies the Rayleigh-Ritz variational method to quantum oscillators in the Segal-Bargmann space. It rigorously derives the normalizability condition |α| < 1/2 for generalized Gaussian trial functions ψ(z) = exp(α z² + β z) via convergence analysis of the associated Gaussian integrals. For the harmonic oscillator it recovers the exact ground state when the trial family contains the true solution. For the quartic anharmonic oscillator it obtains a cubic stationarity equation from adaptive Gaussians together with perturbative energy expansions, while monomial trials furnish rigorous upper bounds but lack width adaptability; displaced trials are examined for asymmetric potentials such as x³ + x⁴.

Significance. If the derivations hold, the work supplies a rigorous foundation for variational calculations in the complex Segal-Bargmann space, with exact recovery for the harmonic oscillator and perturbative results beyond first order for anharmonic cases. The explicit normalizability bound obtained from the quadratic-form analysis of the exponent and the cubic stationarity equation are concrete strengths that could facilitate both analytical and numerical studies of oscillator problems.

minor comments (3)
  1. [Abstract] Abstract: the perturbative expansion E_0 ≈ ½ + 3μ/4 − 9λ²/4 + ⋯ is stated without specifying the orders in λ and μ or indicating whether additional terms are retained; a brief clarification of the truncation would improve readability.
  2. [Monomial trial functions] The monomial upper bound E_n = n + ½ + (3λ/4)(2n² + 2n + 1) is presented as rigorous; it would be useful to state explicitly whether this expression is exact to first order in λ or includes higher-order contributions from the variational procedure.
  3. [Introduction / Segal-Bargmann space] Notation: the measure d²z/π in the Segal-Bargmann inner product is used without a reminder of its normalization; a single sentence recalling the standard convention would aid readers unfamiliar with the space.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive evaluation of our manuscript and the recommendation for minor revision. The recognition of the rigorous normalizability bound, exact recovery for the harmonic oscillator, and the cubic stationarity equation for the anharmonic case is appreciated. No specific major comments were provided in the report, so we will incorporate any minor clarifications on the derivations in the revised version.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The central derivation of the normalizability condition |α| < 1/2 proceeds via direct quadratic-form analysis of the exponent 2 Re(α z²) - |z|² in the Segal-Bargmann integral, producing the matrix trace -2 < 0 and determinant 1 - 4|α|² > 0. This is a standard, parameter-free mathematical step that does not reduce to any fitted input, self-citation, or ansatz imported from prior work. The subsequent variational applications to oscillators use this condition as an independent prerequisite rather than re-deriving it from the results. No load-bearing step collapses by construction to the paper's own inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claims rest on the standard Rayleigh-Ritz variational principle together with convergence of complex Gaussian integrals; no new physical entities are introduced and the variational parameters are optimized rather than fitted to external data.

free parameters (1)
  • alpha and beta in Gaussian trial function
    Variational parameters that are optimized to stationarize the energy expectation value; they are not pre-chosen constants or data-fitted coefficients.
axioms (1)
  • domain assumption Gaussian integrals over the complex plane converge when |α| < 1/2
    Invoked to guarantee normalizability of the trial functions and validity of the variational procedure.

pith-pipeline@v0.9.0 · 5597 in / 1336 out tokens · 64040 ms · 2026-05-15T17:54:41.181471+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

42 extracted references · 42 canonical work pages

  1. [1]

    Position Representation Solution The Hamiltonian for a one-dimensional quantum har- monic oscillator is given by ˆH=− ℏ2 2m d2 dx2 + 1 2 mω2x2.(14) The exact ground-state energy is known to be E0 = 1 2 ℏω,(15) with the corresponding normalized wavefunction ψ0(x) = mω πℏ 1/4 exp − mωx2 2ℏ .(16) We now apply the Rayleigh–Ritz variational method to approxima...

    work page

  2. [2]

    The harmonic oscillator Hamiltonian is ˆH=ℏω(ˆa †ˆa+1 2)

    Complex Plane Solution Consider the normalized coherent state as a trial func- tion in the Segal–Bargmann space: ψα(z) = exp αz− 1 2 |α|2 , α∈C,(27) which corresponds to the standard coherent state|α⟩= e−|α|2/2eαˆa† |0⟩. The harmonic oscillator Hamiltonian is ˆH=ℏω(ˆa †ˆa+1 2). Using⟨ˆa†ˆa⟩=|α|2 for coherent states, the energy expectation value is E(α) =ℏ...

    work page

  3. [3]

    Position Space Solution We apply the Rayleigh–Ritz variational method to the quartic anharmonic oscillator, whose Hamiltonian (in unitsℏ=m=ω= 1) is ˆH=− 1 2 d2 dx2 + 1 2 x2 +λx 4, λ >0.(32) In position space, we adopt the normalized Gaussian trial wavefunction ψ(x;α) = α π 1/4 e−αx2/2, α >0,(33) which respects the even parity and asymptotic decay of the t...

    work page

  4. [4]

    Using⟨x⟩= √ 2αand Wick’s theorem for Gaussian 6 states, the fourth moment is⟨x 4⟩= 4α 4 + 6α2 + 3

    Complex Plane Solution In the Segal–Bargmann space, the Hamiltonian for the quartic anharmonic oscillator (in unitsℏ=m=ω= 1) is ˆH=z d dz + 1 2 +λ z+∂ z√ 2 4 .(38) We consider three trial functions: (i) The coherent stateψ α(z) =e αz−|α|2/2 withα∈ R. Using⟨x⟩= √ 2αand Wick’s theorem for Gaussian 6 states, the fourth moment is⟨x 4⟩= 4α 4 + 6α2 + 3

    work page

  5. [5]

    Forλ >0, the only real solution isα= 0, yieldingE= 1 2 + 3 4 λ, identical to first-order perturbation theory

    The energy functional is E(α) =α 2 + 1 2 +λ 4α4 + 6α2 + 3 4 .(39) Minimization givesdE/dα= 2α+λ(16α 3+12α) = 2α[1+ λ(8α2+6)] = 0. Forλ >0, the only real solution isα= 0, yieldingE= 1 2 + 3 4 λ, identical to first-order perturbation theory. This reflects the symmetry of the potential: the ground state remains centered atx= 0. (ii) The monomialψ(z) =z(first...

    work page

  6. [6]

    Nevertheless, it provides a valid variational upper bound

    This ansatz introduces phase-space anisotropy (⟨x2⟩ ̸=⟨p 2⟩), which is not present in the true ground state of an isotropic potential. Nevertheless, it provides a valid variational upper bound. The norm is ⟨ψ|ψ⟩= (1−4α 2)−1/2, and expectation values respecting z→ −zsymmetry (only even powers ofα) are: ⟨z∂z⟩= 4α2 1−4α 2 ,(40) ⟨x2⟩= 1 2 1 + 4α2 1−4α 2 ,(41)...

    work page

  7. [7]

    Normalizability of Generalized Gaussian T rial F unctions We derive the condition|α|< 1 2 for the general- ized Gaussianψ(z) =e αz2+βz to belong to the Segal– Bargmann space. The squared norm is ∥ψ∥2 = 1 π Z C e2 Re(αz2+βz)−|z| 2 d2z.(A1) Writingz=x+iywithx, y∈Rand decomposingα= a+ib,β=c+idfor reala, b, c, d, we compute: Re(αz2) =a(x 2 −y 2)−2bxy,(A2) Re(...

    work page

  8. [8]

    The linear termv ⊤xshifts the Gaussian center but does not affect convergence, soβ∈C remains unrestricted. 9

    work page

  9. [9]

    Phase-space anisotropy of generalized Gaussian trial functions To rigorously establish that the generalized Gaus- sian trial functionψ(z) =e αz2+βz induces phase-space anisotropy forα̸= 0, we compute expectation values of position and momentum squared in the Segal–Bargmann space. Settingβ= 0 without loss of generality (as dis- placement does not affect an...

    work page

  10. [10]

    Forx= (z+∂ z)/ √ 2: ⟨x2⟩= 1 2 ⟨z∂z⟩+⟨z∂ z + 1⟩ = 1 2 2⟨z∂z⟩+ 1 = 1 2 1 + 4α2 1−4α 2 ,(A23) using⟨z 2⟩=⟨∂ 2 z ⟩= 0 and [∂ z, z] = 1

    Expectation V alue of⟨x 4⟩for Generalized Gaussian (Squeezed-State) T rial F unctions Forψ(z) =e αz2 with realαand|α|< 1 2, the norm is ⟨ψ|ψ⟩= (1−4α 2)−1/2.(A21) For the number operator expectation: ⟨z∂z⟩= 2α ⟨ψ|ψ⟩ 1 π Z C |z|2e2α(x2−y2)−|z|2 d2z = 4α2 1−4α 2 ,(A22) which contains only even powers ofα, respectingz→ −z symmetry. Forx= (z+∂ z)/ √ 2: ⟨x2⟩= 1...

    work page

  11. [11]

    Monomial Expectation V alues for Quartic Potential Forϕ n(z) =z n/ √ n! andx= (z+∂ z)/ √ 2, we expand x4 = 1 4(z+∂ z)4.(A25) Using [∂z, z] = 1 and normal ordering :· · ·:: z3∂z =:z 3∂z : +3z2,(A26) z2∂2 z =:z 2∂2 z : +4z∂z + 2.(A27) Only c-number terms contribute to expectation values: ⟨z3∂z⟩= 3n,(A28) ⟨z2∂2 z ⟩= 4n+ 2.(A29) Summing all contributions and ...

    work page

  12. [12]

    (A31) becomesu 3 +v 3 = 6λwithu 3v3 = 1

    Solution of the Cubic Stationarity Equation The stationarity condition yields α3 −α−6λ= 0.(A31) Using Cardano’s substitutionα=u+vwithuv= 1 3, Eq. (A31) becomesu 3 +v 3 = 6λwithu 3v3 = 1

    work page

  13. [13]

    The discriminant is ∆ = 36λ2 − 4 27 = 4 9λ2 − 1 27 ,(A32) so u3, v3 = 3λ± q 9λ2 − 1 27 .(A33) Taking real cube roots gives αopt = 3 r 3λ+ q 9λ2 − 1 27 + 3 r 3λ− q 9λ2 − 1 27 .(A34)

    Thus u3 andv 3 are roots oft 2 −6λt+ 1 27 = 0. The discriminant is ∆ = 36λ2 − 4 27 = 4 9λ2 − 1 27 ,(A32) so u3, v3 = 3λ± q 9λ2 − 1 27 .(A33) Taking real cube roots gives αopt = 3 r 3λ+ q 9λ2 − 1 27 + 3 r 3λ− q 9λ2 − 1 27 .(A34)

    work page

  14. [14]

    Solving pertur- batively withδ=a 1λ+a 2λ2 +· · ·: •O(λ):a 1 = 3 2, •O(λ 2):a 2 =− 81 16

    Perturbative Energy Expansion for Quartic Oscillator Substitutingα= 1 +δwithδ=O(λ) intoE(α) = α 4 + 1 4α + 3λ 4α2 : E= 1 2 + δ2 2 + 3λ 4 − 3λδ 2 +O(λ 3, δ3).(A35) Minimizing givesδ− 3λ 2 +O(λδ, δ 2) = 0. Solving pertur- batively withδ=a 1λ+a 2λ2 +· · ·: •O(λ):a 1 = 3 2, •O(λ 2):a 2 =− 81 16. Thusδ= 3 2 λ− 81 16 λ2 +O(λ 3). Substituting back: E= 1 2 + 3λ 4...

    work page

  15. [15]

    Energy F unctional for Displaced Monomials For the displaced coherent stateψ 0(z;γ) =e γz−|γ| 2/2 with realγ, Wick’s theorem gives: ⟨x⟩= √ 2γ,(A37) ⟨x2⟩= 2γ 2 + 1 2 ,(A38) ⟨x3⟩= 2 √ 2γ 3 + 3 √ 2γ,(A39) ⟨x4⟩= 4γ 4 + 6γ2 + 3 4 .(A40) The harmonic energy is⟨ ˆH0⟩=γ 2 + 1

    work page

  16. [16]

    Combining all terms: E0(γ) =γ 2 + 1 2 +λ 2 √ 2γ 3 + 3 √ 2γ +µ 4γ4 + 6γ2 + 3 4 .(A41) Forλ≪1, perturbation theory withγ=b 1λ+O(λ 2) yieldsb 1 =− 3 2 and E0,min = 1 2 + 3µ 4 − 9λ2 4 +O(λ 3),(A42) confirming displacement-induced stabilization

    work page

  17. [17]

    In the Segal–Bargmann space, parity acts asz7→ −z,∂ z 7→ −∂z

    Symmetry Constraints in Holomorphic V ariational Calculations For Hamiltonians with even parity ([ ˆH, ˆP] = 0 where ˆP x ˆP −1 =−x), the ground state has definite parity. In the Segal–Bargmann space, parity acts asz7→ −z,∂ z 7→ −∂z. Trial functions must respect this symmetry to avoid unphysical contributions. Considerψ(z;α) =e αz2 with realα: •Underz→ −z...

    work page

  18. [18]

    A. B. Migdal, and V. P. Krainov,Approximation Methods in Quantum Mechanics, W. A. Benjamin, (1969)

    work page 1969

  19. [19]

    A. S. Davydov,Quantum Mechanics, 2nd edition, Perg- amon Press (1976)

    work page 1976

  20. [20]

    A. L. Kuzemsky,Statistical Mechanics And The Physics Of Many-particle Model Systems, World Scientific, Pub- lishing (2017)

    work page 2017

  21. [21]

    F. M. Fern´ andez,On the Rayleigh-Ritz variational method,arXiv:2206.05122(2022)

    work page arXiv 2022

  22. [22]

    F. M. Fern´ andez,On the Raleigh–Ritz variational method. Non-orthogonal basis set, Journal of Mathemat- ical Chemistry62, 2083–2089 (2024)

    work page 2083

  23. [23]

    Borowski, K

    P. Borowski, K. Kowalska, and P. Pikuta,The Rayleigh–Ritz Variation Method: An Illustrative Applica- tion to Anharmonicity of the Hydrogen Chloride Molecule Vibrations, Journal of Chemical Education102(2025) 1529

    work page 2025

  24. [24]

    Echenique, J

    P. Echenique, J. L. Alonso,A mathematical and compu- tational review of Hartree-Fock SCF methods in Quantum Chemistry, Molecular Physics105(2007) 3057-3098

    work page 2007

  25. [25]

    W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Ra- jagopal,Quantum Monte Carlo simulations of solids, Rev. Mod. Phys.73, 33 (2001)

    work page 2001

  26. [26]

    Lipkin, N

    N. Lipkin, N. Moiseyev, and E. Br¨ andas,Resonances by the exterior-scaling method within the framework of the finite-basis-set approximation, Phys. Rev. A40, 549 (1989)

    work page 1989

  27. [27]

    Moiseyev,Quantum theory of resonances: calculating energies, widths and cross-sections by complex scaling, Physics Reports302, (1998) 212

    N. Moiseyev,Quantum theory of resonances: calculating energies, widths and cross-sections by complex scaling, Physics Reports302, (1998) 212

    work page 1998

  28. [28]

    Kuro´ s, P

    A. Kuro´ s, P. Ko´ scik, A. Okopi´ nska,Determination of Resonances by the Optimized Spectral Approach, J. Phys. A: Math. Theor.46(2013) 085303

    work page 2013

  29. [29]

    Imakura, K

    A. Imakura, K. Morikuni, and A. Takayasu,Verified eigenvalue and eigenvector computations using complex moments and the Rayleigh–Ritz procedure for general- ized Hermitian eigenvalue problems, Journal of Compu- tational and Applied Mathematics424(2023) 114994

    work page 2023

  30. [30]

    F. M. Fern´ andez, J. Garcia,Highly accurate calculation of the resonances in the Stark effect in hydrogen, Applied Mathematics and Computation317(2017) 550–565

    work page 2017

  31. [31]

    Kuznetsova A., et al ,Spectroscopy of multielectron atomic systems in a DC electric field, Advances in Quan- tum Chemistry78, 287–306 (2019)

    work page 2019

  32. [32]

    Fleig,Variational calculation of the hyperfine Stark ef- fect in atomic systems, arXiv:2409.08787 [physics.atom- ph], 2024

    T. Fleig,Variational calculation of the hyperfine Stark ef- fect in atomic systems, arXiv:2409.08787 [physics.atom- ph], 2024

    work page arXiv 2024

  33. [33]

    Singh, P

    A. Singh, P. Siwach, P. Arumugam,Quantum simu- lations of nuclear resonances with variational methods arXiv:2504.11685 [quant-ph], 2025

    work page arXiv 2025

  34. [34]

    I. E. Segal,Mathematical problems of relativistic physics, in Kac, M. (ed.), Proceedings of the Summer Seminar, Boulder, Colorado, 1960, Vol. II, Lectures in Applied Mathematics, American Mathematical Society (1963)

    work page 1960

  35. [35]

    Bargmann,On a Hilbert space of analytic functions and an associated integral transform part I, Commun

    V. Bargmann,On a Hilbert space of analytic functions and an associated integral transform part I, Commun. Pure. Appl. Math.14(3): 187 (1961)

    work page 1961

  36. [36]

    Voros,Wentzel-Kramers-Brillouin method in the Bargmann representation, Phys

    A. Voros,Wentzel-Kramers-Brillouin method in the Bargmann representation, Phys. Rev. A40, 6814 (1989)

    work page 1989

  37. [37]

    Vourdas, and R

    A. Vourdas, and R. F. Bishop,Thermal coherent states in the Bargmann representation, Phys. Rev. A50, 3331 (1994)

    work page 1994

  38. [38]

    M. W. AlMasri, and M. R. B. Wahiddin,Bargmann representation of quantum absorption refrigerators, Rep. Math. Phys.89(2), Pages 185-198 (2022)

    work page 2022

  39. [39]

    Perelomov,Generalized Coherent States and Their Applications, Springer Berlin, Heidelberg (1986)

    A. Perelomov,Generalized Coherent States and Their Applications, Springer Berlin, Heidelberg (1986)

    work page 1986

  40. [40]

    G. B. Folland,Harmonic Analysis in Phase Space, Princeton University Press (1989)

    work page 1989

  41. [41]

    Zinn-Justin,Path Integrals in Quantum Mechanics, Oxford Graduate Texts, Oxford University Press (2010)

    J. Zinn-Justin,Path Integrals in Quantum Mechanics, Oxford Graduate Texts, Oxford University Press (2010)

    work page 2010

  42. [42]

    P. K. Patnaik,Rayleigh-Schr¨ odinger perturbation theory for the anharmonic oscillator, Phys. Rev. D35, 1234 ( 1987)

    work page 1987