## Analytic Model for the Energy Spectrum of the Anharmonic Oscillator

Michel Caffarel 1, ∗

1 Laboratoire de Chimie et Physique Quantiques (UMR 5626), Universit´ de Toulouse, CNRS, UPS, France e

In a recent work we have proposed an original analytic expression for the partition function of the quartic oscillator. This partition function, which has a simple and compact form with no adjustable parameters , reproduces some key mathematical properties of the exact partition function and provides free energies accurate to a few percent over a wide range of temperatures and coupling constants. In this work, we present the derivation of the energy spectrum of this model. We also generalize our previous study limited to the quartic oscillator to the case of a general anharmonic oscillator. Numerical application for a potential of the form V ( x ) = ω 2 2 x 2 + gx 2 m show that the energy levels are obtained with a relative error of about a few percent, a precision which we consider to be quite satisfactory given the simplicity of the model, the absence of adjustable parameters, and the negligible computational cost.

The one-dimensional quantum anharmonic oscillator[1] plays an important role in quantum mechanics as a simple yet nontrivial model for describing nonlinear (anharmonic) effects. As such, it is commonly used in various scientific fields including molecular physics (rovibrational spectra), condensed matter physics (phonons), nuclear physics (collective vibrational motions of nuclei), and quantum field theory ( e.g. , ϕ 4 -theory, Higgs mechanism), to cite the main ones. No exact solution for a general anharmonic potential has been found so far, although some partial solutions have been developed for specific potentials, for example in the case of the so-called quasi-exactly-solvable problems; see Turbiner[2] and references therein. From a numerical point of view, the simplicity of the ordinary differential equation to solve allows very accurate solutions. To cite one numerical approach among many, let us mention the recently developed Lagrange Mesh method of Baye [3] as recently implemented by del Valle[4]. However, as for any physical model, having analytical solutions -even approximate ones- is important since it may lead a deeper insight into the nature of the problem, reveal underlying hidden structures that may not be evident from numerical solutions and, also, guide further analysis or generalizations. A great variety of approximate analytical approaches have been proposed, making it very difficult to provide an exhaustive account of the literature. In this work, we are more specifically interested in evaluating the energy levels of the anharmonic oscillator. Among the main approaches developed for this purpose, let us mention the semi-classical approaches, such as WKB[5, 6] or the phase-integral method based on the generalized Bohr-Sommerfeld quantization condition ([7] and references therein) and the methods based on the use of a variational approach either using parameterized excited wavefunctions[8] or path integrals[9, 10]. Other interesting approaches include the design of simple analytic expressions for the energy levels built from the weakand/or strong-coupling expansions[11, 12] and

∗ caffarel@irsamc.ups-tlse.fr

methods based on the tuning of boundary conditions[13] or eigenvector continuation [14].

The aim of this work is to present novel analytical expressions for the energy levels of the anharmonic oscillator. For that, we take advantage of a recently proposed partition function for the quartic oscillator,[15] which we extend here to an arbitrary potential. As shown in our previous study[15], this partition function provides an appealing simple and physically meaningful model for the exact solution. Indeed, it has a simple and compact form with no adjustable parameters , a desirable property for an analytical model. Furthermore, it reproduces some key mathematical properties of the exact partition function, thus supporting the idea that the important features of the non-trivial mathematical structure of the solution are, at least partially, accounted for.

These properties are the following: i) The harmonic and classical limits are exactly recovered, ii) the well-known divergence of the weak-coupling (RayleighSchr¨dinger) expansion of the energy is reproduced. o As for the exact solution, the energy corrections are found to be rational numbers and display a factorial-like growth in terms of the perturbational order, and iii) the functional form of the strong-coupling expansion is recovered.

From a quantitative point of view, the free energy is found to be accurate to a few percent over a wide range of temperatures and coupling constants. A similar precision is also obtained for the ground- and first-excited state energies.

To the best of our knowledge, no partition function in closed-form proposed so far is capable of simultaneously reproducing all these features (see, for example, references [9, 16-18] and a comparative study with other models in [15]).

Having at our disposal a simple and physically meaningful model for the partition function, it is of interest to derive the full energy spectrum of the model. This is the purpose of the present work. As we shall see, the energy levels obtained turn out to be accurate with a relative error of about a few percent, a precision which we consider to be quite satisfactory for such a simple model. Of course, a higher precision can be obtained without

difficulty by using some of the numerical/analytical approaches cited above. However, we emphasize that the aim of this work is not to achieve an ultimate precision in the spectrum but, instead, to propose a simple yet faithful analytical model for the spectrum of a general anharmonic oscillator that can be readily applied in various scientific contexts with essentially no computational cost.

The Hamiltonian considered here is as follows

$$H = - \frac { 1 } { 2 } \frac { d ^ { 2 } } { d x ^ { 2 } } + V ( x ) \quad \quad ( 1 ) \quad \text{i}$$

where V ( x ) is a rather general potential function bounded from below and verifying lim | x |→∞ V ( x ) = + ∞ (in other words, there are only discrete energies). In the numerical applications presented below, we will particularize the potential in the form

$$V ( x ) = \frac { \omega ^ { 2 } } { 2 } x ^ { 2 } + g x ^ { 2 m } \quad \quad ( 2 ) \quad \text{Fo} \\ \text{for} \\ \text{...} \colon$$

where ω 2 and g denote the harmonic force constant and coupling constant, respectively.

As a first step let us briefly recap the main steps leading to the model partition function.

i) The partition function is first expressed as a path integral in a standard way (see, e.g. , [19])

$$Z = \text{Tr} e ^ { - \beta H } = \lim _ { n \to \infty } Z _ { n } \quad \quad ( 3 ) \quad \text{iv} )$$

with

$$Z _ { n } & = \left ( \frac { 1 } { \sqrt { 2 \pi \tau } } \right ) ^ { n } \int _ { - \infty } ^ { \infty } d x _ { 1 } \dots \int _ { - \infty } ^ { \infty } d x _ { n } \\ & e ^ { - \frac { 1 } { 2 \tau } \sum _ { i = 1 } ^ { n } ( x _ { i + 1 } - x _ { i } ) ^ { 2 } } e ^ { - \tau \sum _ { i = 1 } ^ { n } V ( x _ { i } ) }. \quad \quad ( 4 ) \quad \text{with}$$

where τ = β n is the time-step and periodic conditions are used, x n +1 = x 1 .

ii) Second, the short-time anharmonic contribution, e -τV ( x ) , is approximated by a gaussian distribution centered on some position x ∗ (typically, the position of the lowest minimum of the potential; here, x ∗ = 0 in our applications) with an effective frequency ω g ( τ ), that is

$$\frac { e ^ { - \tau V ( x ) } } { \int _ { - \infty } ^ { \infty } d x e ^ { - \tau V ( x ) } } \sim \frac { e ^ { - \tau \frac { 1 } { 2 } \omega _ { g } ^ { 2 } ( \tau ) ( x - x ^ { * } ) ^ { 2 } } } { \int _ { - \infty } ^ { \infty } d x e ^ { - \tau \frac { 1 } { 2 } \omega _ { g } ^ { 2 } ( \tau ) ( x - x ^ { * } ) ^ { 2 } } }. \quad ( 5 ) \quad \text{ner}$$

To set the frequency ω g ( τ ) we have proposed to impose to the two distributions to have the same variance. In the case of the quartic oscillator treated previously,[15] it leads to ω g = ω √ [ B 4 g τω 4 ] where B is some function. In the more general case V ( x ) = ω 2 2 x 2 + gx 2 m , we get

$$\omega _ { g } ( \tau ) = \omega \sqrt { B ^ { ( m ) } \left [ \frac { 2 ^ { m } g } { \tau ^ { m - 1 } \omega ^ { 2 m } } \right ] } \quad \quad ( 6 ) \quad \text{ of } \colon$$

where the parameter-free function B ( m ) ( x ) is given by

$$B ^ { ( m ) } ( x ) = \frac { 1 } { 2 } \frac { \int _ { - \infty } ^ { \infty } d y \, e ^ { - y ^ { 2 } - x y ^ { 2 m } } } { \int _ { - \infty } ^ { \infty } d y \, y ^ { 2 } e ^ { - y ^ { 2 } - x y ^ { 2 m } } }. \quad \quad ( 7 )$$

For a general potential V the formula writes

$$\omega _ { g } ( \tau ) = \sqrt { \frac { 1 } { \tau } \frac { \int _ { - \infty } ^ { \infty } d x e ^ { - \tau V ( x ) } } { \int _ { - \infty } ^ { \infty } d x ( x - x ^ { * } ) ^ { 2 } e ^ { - \tau V ( x ) } } } \quad ( 8 )$$

iii) The gaussian approximation being made, the infiniten limit of Z n is no longer defined, lim n →∞ Z n = + ∞ . To circumvent this problem, we have proposed to introduce a Principle of Minimal Sensitivity (PMS) for the path integral. More precisely, we impose to the path integral to minimally depend on the effective frequency used in the gaussian approximation, that is, ∂Z n ∂ω g ( τ ) = 0. For a given temperature, the PMS condition holds only for a unique value of n denoted as n c ( β ) [and, thus, a unique time-step, τ c ( β ) = β n c ( β ) ]. We have then proposed to define the model partition function as the value of Z n at this 'optimal' value of n , Z ≡ Z n c ( β ) . In particular, the ill-defined n → ∞ -limit is avoided. The nonlinear implicit equation for n c ( β ) resulting from the PMS condition writes

$$n _ { c } ( \beta ) = \frac { \beta \omega _ { g } [ \tau _ { c } ( \beta ) ] } { 2 } \coth \frac { \beta \omega _ { g } [ \tau _ { c } ( \beta ) ] } { 2 }. \quad \ ( 9 )$$

iv) Finally, the analytical expression of the model partition function is given by

$$Z = \frac { C ( \beta ) ^ { n _ { c } ( \beta ) } } { e ^ { \frac { \beta \omega _ { g } [ \tau _ { c } ( \beta ) ] } { 2 } } - e ^ { - \frac { \beta \omega _ { g } [ \tau _ { c } ( \beta ) ] } { 2 } } } \quad \quad ( 1 0 )$$

with and

$$I _ { g } ( \beta ) = \int _ { - \infty } ^ { \infty } d x e ^ { - \tau _ { c } ( \beta ) V ( x ) }, \text{ \quad \ \ } ( 1 2 )$$

where the time-step τ c ( β ) = β n c ( β ) is obtained from Eq.(9).

At this point, an important remark is in order. Introducing a gaussian approximation of a nongaussian quantity is a standard practice in physics (harmonic phonons, Gaussian approximation for the Ginzburg-Landau action, etc.) In short, it is done by restricting to the second-order a Taylor expansion of some Hamiltonian or action. Here, our gaussian approximation is very different in nature. Instead of approximating the potential V ( x ) by a quadratic potential, we approximate the quantity e -τV ( x ) by a

$$C ( \beta ) = \sqrt { \frac { \omega _ { g } [ \tau _ { c } ( \beta ) ] } { \pi \coth \frac { \beta \omega _ { g } [ \tau _ { c } ( \beta ) ] } { 2 } } } I _ { g } ( \beta ) \quad \quad ( 1 1 )$$

gaussian distribution with an effective frequency which depends explicitly on the time-step . It is this dependence on the time-step that makes Z n to diverge in the largen limit, but which, after application of the PMS condition, allows to fix n at a finite value n c ( β ) and, then, leads to an accurate model for the exact partition function. Approximating V ( x ) using a quadratic potential independent on the time-step would merely lead to the partition function of a simple harmonic oscillator and, thus, to a very poor model for the exact partition function.

A. Ground-state energy. The ground-state energy is obtained from the largeβ behavior of Z

$$E _ { 0 } ( g ) = \lim _ { \beta \to \infty } - \frac { 1 } { \beta } \ln Z. \quad \quad ( 1 3 ) \quad \ t h e _ { \text{rat} }$$

At large β , the (unique) solution n c ( β ) of Eqs.(6) and (9) is proportional to β and given by

$$n _ { c } ( \beta ) = \frac { \beta \bar { \omega } _ { g } } { 2 }, \quad \quad \quad ( 1 4 ) \quad \quad \text{by}$$

where ¯ ω g = ω g [ τ c (+ ∞ )] is the solution of the implicit equation given by

$$\bar { \omega } _ { g } = \omega \sqrt { B ^ { ( m ) } \left ( \frac { 2 g \bar { \omega } _ { g } ^ { m - 1 } } { \omega ^ { 2 m } } \right ) }. \quad \quad ( 1 5 ) \quad \text{Let}$$

Using the zero-temperature limit of the free energy and the expression for the partition function we get

$$E _ { 0 } ( g ) = \frac { \bar { \omega } _ { g } } { 2 } \left [ 1 - \ln \left ( \sqrt { \frac { \bar { \omega } _ { g } } { \pi } } \bar { I } _ { g } \right ) \right ] \quad \text{(1)} \quad \text{ders} \\ \text{the}$$

where

$$\bar { I } _ { g } = \int _ { - \infty } ^ { \infty } d x e ^ { - \frac { 2 } { \omega _ { g } } V ( x ) }. \quad \quad ( 1 7 ) \quad \begin{smallmatrix} \Delta _ { \omega _ { g } } ^ { \cdot } \\ \omega _ { g } | \end{smallmatrix}$$

B. Excited-state energies. The exact partition function of the quartic oscillator decomposes as a discrete sum of exponentials

$$Z = \sum _ { n = 0 } ^ { \infty } e ^ { - \beta E _ { n } } \quad \quad \quad ( 1 8 ) \quad \text{th}.$$

where E n are the excited-state energies. As just seen, the ground-state energy E 0 is obtained by extracting the leading exponential component of the PF at large β . To get excited-state energies subleading components are to be evaluated.

In the following we will show that the model partition function actually does not write as a sum of simple exponentials, as it should be for the exact PF, but, instead, as a sum of exponentials multiplied by a polynomial term as follows

$$Z = \sum _ { n = 0 } ^ { \infty } P _ { n } ( \beta ) e ^ { - \beta ( E _ { 0 } + n \bar { \omega } _ { g } ) }$$

where P n ( β ) is a polynomial of degree n in β with P n (0) = 1. The partition function being no longer expressed as a sum of simple exponentials, the definition of what is meant by excited-state energies becomes problematic. Comparisons with the 'exact' numerical spectrum show that defining the excited energies as the exponents of the exponential contributions, that is,

$$E _ { n } = E _ { 0 } + n \bar { \omega } _ { g } \text{ \quad \ \ } ( 2 0 )$$

gives a very poor approximation of the exact spectrum. Note that this is not surprising since, in such a case, the energy differences between successive states would remain constant, a property which is clearly wrong for the exact spectrum. In sharp contrast, we have found that modifying the partition function by exponentiating the linear contribution of the polynomials and incorporating it into the exponential part leads to remarkably good energy levels. Precisely, we propose to replace the polynomial

$$P _ { n } ( \beta ) = 1 + P _ { n 1 } \beta + \dots + P _ { n n } \beta ^ { n } \quad \ \ ( 2 1 )$$

by e βP n 1 and to define the excited-state energies of our model as

$$E _ { n } = E _ { 0 } + n \bar { \omega } _ { g } - P _ { n 1 }. \text{ \quad \ \ } ( 2 2 )$$

Let us emphasize that this replacement should not be considered as a quantitative approximation, but rather, as a 'minimal' modification of the model to impose to the partition function to have the exact expansion, Eq.(18). Unfortunately, we have not been able to understand why this simple additional prescription to our model is so effective in giving accurate energy levels (see, the figures to follow). However, it should be considered as a salient result of this work.

To derive the functional form of the partition function, Eq.(19), we first need to expand the effective frequency ω g [ τ c ( β )] as a power series of exponential-like contributions.

̸

Let us first consider the case ω = 0. Using Eqs.(6) and (9), and introducing the variable y defined as

$$y = e ^ { - \beta \omega _ { g } }$$

the implicit equation obeyed by ω g is rewritten as

$$\frac { 1, } { g } \sum _ { \substack { g \\ \omega } } \omega _ { g } = \omega \sqrt { B \left [ \frac { 2 g } { \omega ^ { 2 m } } \omega _ { g } ^ { m - 1 } \left ( 1 + 2 \sum _ { n = 1 } ^ { \infty } y ^ { n } \right ) ^ { m - 1 } \right ] } \quad ( 2 4 )$$

where the coth function has been expanded in power series of y . Note that the superscript ( m ) has been removed from B ( m ) ( x ) to simplify the notation. This will also be the case in the following for most quantities depending on m when no confusion is possible. The function B being infinitely differentiable, ω g can be expanded in powers of y

$$\omega _ { g } = \bar { \omega } _ { g } + \sum _ { n = 1 } ^ { \infty } \omega _ { n } y ^ { n }. \text{ \quad \ \ } ( 2 5 )$$

Remark that the equation obeyed by ω g has no explicit dependence on β (the dependence on the temperature is only through the variable y ). Accordingly, the coefficients ω n are independent of β . Their evaluation is done i) by introducing the expansion of ω g in Eq.(24), ii) by Taylor-expanding the function B at x 0 = 2 g ω 2 m ¯ ω m -1 g , and, finally, by identifying the contributions corresponding to a given power of y . After some algebra, an explicit expression for the ω n 's can be derived. For example, the first coefficient, ω 1 's is given by

$$\omega _ { 1 } = \bar { \omega } _ { g } \frac { \frac { B _ { 1 } } { B _ { 0 } } x _ { 0 } ( m - 1 ) } { 1 - \frac { \frac { B _ { 1 } } { B _ { 0 } } x _ { 0 } } { 2 } ( m - 1 ) } \quad \quad ( 2 6 ) \quad \text{To}$$

where B 0 = B ( m ) ( x 0 ) and B 1 = dB ( m ) dx ( x 0 ). For n ≥ 2, ω n can be expressed as a function of the preceding coefficients, ω p with 1 ≤ p ≤ n -1. The explicit formula is given in Appendix A.

When ω = 0 the coefficients ω k are more easily derived. From the general expression of the effective frequency, Eq.(8), we get ω g ( τ ) = c 0 ( g τ ) 1 4 with c 0 = √ ∫ dxe -x 4 ∫ dxx e 2 -x 4 . Using the PMS condition, Eq.(9), we obtain

$$\omega _ { k } = \bar { \omega } _ { g } \sum _ { l = 1 } ^ { k } 2 ^ { l } \binom { \frac { 1 } { 3 } } { l } \binom { k - 1 } { l - 1 } \ k \geq 1 \quad \ ( 2 7 )$$

with

$$\bar { \omega } _ { g } = c _ { 0 } ^ { \frac { 4 } { 3 } } \left ( \frac { g } { 2 } \right ) ^ { \frac { 1 } { 3 } }. \quad \quad \quad ( 2 8 )$$

In this ω = 0-case, the ground-state energy can be explicitly written as

$$E _ { 0 } = \left [ \frac { \Gamma \left ( \frac { 5 } { 4 } \right ) } { \Gamma \left ( \frac { 3 } { 4 } \right ) } \right ] ^ { \frac { 2 } { 3 } } \left ( 1 - \log \sqrt { \frac { \left [ 2 \Gamma \left ( \frac { 5 } { 4 } \right ) \right ] ^ { 3 } } { \pi \Gamma \left ( \frac { 3 } { 4 } \right ) } } \right ) g ^ { \frac { 1 } { 3 } }. \quad ( 2 9 )$$

Now, to proceed, let us introduce the following new variable y 0

$$y _ { 0 } = e ^ { - \beta \bar { \omega } _ { 9 } }. \quad \quad \ \ ( 3 0 ) \quad \ \ a$$

The next step consists in expressing the partition function as a power series in y 0 , thus leading to Eq.(19). For that, we first need to derive the expansion of the effective frequency in terms of y 0 . We have

$$\omega _ { g } = \bar { \omega } _ { g } + \sum _ { n = 1 } ^ { \infty } \omega _ { n } y _ { 0 } ^ { n } e ^ { - \beta n [ \omega _ { g } - \bar { \omega } _ { g } ] }. \quad \quad ( 3 1 )$$

By expanding the exponential in power series and by simple inspection, the form of the solution is

$$\omega _ { g } = \bar { \omega } _ { g } + \sum _ { n = 1 } ^ { \infty } Q _ { n } ( \beta ) y _ { 0 } ^ { n } \text{ \quad \ \ } ( 3 2 ) \text{ \quad,}$$

where Q n ( β ) are polynomials of degree n -1 in β . The polynomials can be evaluated by differentiation

$$Q _ { n } ( \beta ) = \frac { 1 } { n! } \frac { \partial ^ { n } \omega _ { g } } { \partial y _ { 0 } ^ { n } } ( y _ { 0 } = 0 ). \text{ \quad \ \ } ( 3 3 )$$

Using Eq.(31) and the Leibniz formula for derivatives, the n -th derivative of ω g writes

$$\text{le, the} \\ \quad & \frac { 1 } { n! } \frac { \partial ^ { n } \omega _ { g } } { \partial y _ { 0 } ^ { n } } ( y _ { 0 } = 0 ) = \sum _ { k = 1 } ^ { n } \omega _ { k } \frac { 1 } { ( n - k )! } \frac { \partial ^ { n - k } } { \partial y _ { 0 } ^ { n - k } } \left [ e ^ { - \beta k ( \omega _ { g } - \bar { \omega } _ { g } ) } \right ] ( y _ { 0 } = 0 ). \\ \quad & \text{$T$} \text{ procedure two make one of the $E$} \, \text{$di$} \, \text{Runno formula}.$$

To proceed we make use of the Fa` di Bruno formula a

$$\underset { \substack { r \ n \geq \\ \coding } } { \underset { \substack { \longrightarrow \\ \text{remula} } } { \sum } } \frac { \partial ^ { n } } { \partial x ^ { n } } e ^ { f ( x ) } = e ^ { f ( x ) } \sum _ { m _ { 1 }, m _ { 2 } \dots, m _ { n } } ^ { \prime } \frac { n! } { m _ { 1 }! m _ { 2 }! \dots m _ { n }! } \prod _ { j = 1 } ^ { n } \left [ \frac { f ^ { ( j ) } ( x ) } { j! } \right ] ^ { m _ { j } } \\ \text{where the prime on the sum indicates the sum with the } \text{$\quad$}$$

where the prime on the sum indicates summing with the constraint

$$1 m _ { 1 } + 2 m _ { 2 } + \, \dots \, n m _ { n } = n. \quad \quad ( 3 6 )$$

We then have

$$\frac { 1 } { n! } \frac { \partial ^ { n } \omega _ { g } ( \beta ) } { \partial y _ { 0 } ^ { n } } ( y _ { 0 } = 0 ) = \sum _ { k = 1 } ^ { n } \omega _ { k } \frac { 1 } { ( n - k )! }$$

$$\times \sum _ { m _ { 1 }, m _ { 2 } \dots, m _ { n - k } } ^ { \prime } \frac { ( n - k )! } { m _ { 1 }! m _ { 2 }! \dots m _ { n - k }! } \prod _ { j = 1 } ^ { n - k } \left [ \frac { ( - k \beta ) ( \omega _ { g } ) ^ { ( j ) } ( x ) } { j! } \right ] ^ { m _ { j } } ( y _ { 0 } = 0 )$$

and, finally

$$Q _ { n } ( \beta ) = \omega _ { n } + \sum ^ { n - 1 } \omega _ { k }$$

$$k = 1$$

$$\text{new} \quad \times \sum _ { m _ { 1 }, m _ { 2 } \dots, m _ { n - k } } ^ { \prime } \frac { \left [ - k \beta Q _ { 1 } ( \beta ) \right ] ^ { m _ { 1 } } } { m _ { 1 }! } \dots \frac { \left [ - k \beta Q _ { n - k } ( \beta ) \right ] ^ { m _ { n - k } } } { m _ { n - k }! },$$

a relation which allows to evaluate iteratively the polynomials Q n ( β ). From this expression, we see that Q n ( β ) are indeed polynomials of degree n -1 as stated above. Let us give the first five polynomials

$$Q _ { 1 } ( \beta ) = \omega _ { 1 }$$

$$Q _ { 2 } ( \beta ) = \omega _ { 2 } - \beta \omega _ { 1 } ^ { 2 }$$

$$Q _ { 3 } ( \beta ) = \omega _ { 3 } - 3 \beta \omega _ { 1 } \omega _ { 2 } + \frac { 3 } { 2 } \beta ^ { 2 } \omega _ { 1 } ^ { 3 }$$

$$Q _ { 4 } ( \beta ) = \omega _ { 4 } - 2 \beta ( 2 \omega _ { 1 } \omega _ { 3 } + \omega _ { 2 } ^ { 2 } ) + 8 \beta ^ { 2 } \omega _ { 2 } \omega _ { 1 } ^ { 2 } - \frac { 8 } { 3 } \omega _ { 1 } ^ { 4 } \beta ^ { 3 }$$

$$Q _ { 5 } ( \beta ) = \omega _ { 5 } - 5 \beta ( \omega _ { 2 } \omega _ { 3 } + \omega _ { 1 } \omega _ { 4 } ) + \frac { 2 5 } { 2 } \beta ^ { 2 } ( \omega _ { 3 } \omega _ { 1 } ^ { 2 } + \omega _ { 1 } \omega _ { 2 } ^ { 2 } ) \quad \text{Let} \quad \text{$\alpha$} \quad \text{$\alpha$}$$

$$- \frac { 1 2 5 \beta ^ { 3 } } { 6 } \omega _ { 1 } ^ { 3 } \omega _ { 2 } + \frac { 1 2 5 } { 2 4 } \omega _ { 1 } ^ { 5 } \beta ^ { 4 }$$

We are now ready to derive the expansion of the partition function in the variable y 0 . It is convenient to introduce the following three quantities

$$\Delta \omega _ { g } \equiv \omega _ { g } - \bar { \omega } _ { g }, \quad \quad \ \ ( 3 7 ) \quad \stackrel { \cdots } { \text{th} }.$$

$$\Delta T \equiv \ln C ( \beta ) - \ln \sqrt { \frac { \bar { \omega } _ { g } } { \pi } } \bar { I } _ { g }, \quad \quad ( 3 8 ) \quad \text{ we }.$$

and

$$\Delta R \equiv \frac { 1 } { 2 } ( \coth \frac { \beta \omega _ { g } } { 2 } - 1 ). \quad \quad ( 3 9 ) \quad \text{$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$\,$$$

The partition function can be written as

$$Z = e ^ { - \beta E _ { 0 } } P.$$

Simple algebra shows that P can be expressed as

$$P = e ^ { \beta S } [ 1 + \Delta R ] \quad \quad \ \ ( 4 1 ) \quad \text{for}$$

where

$$S = c _ { 1 } \Delta \omega _ { g } + c _ { 2 } \Delta T + c _ { 3 } \Delta R$$

$$& + c _ { 4 } \Delta \omega _ { g } \Delta T + c _ { 5 } \Delta \omega _ { g } \Delta R + c _ { 6 } \Delta T \Delta R + c _ { 7 } \Delta \omega _ { g } \Delta T \Delta R \\ & \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots.$$

with the following coefficients

$$c _ { 1 } = \frac { 1 } { 2 } \left ( \ln \sqrt { \frac { \bar { \omega } _ { g } } { \pi } } \bar { I } _ { g } - 1 \right )$$

$$c _ { 2 } = \frac { \bar { \omega } _ { g } } { 2 }$$

$$c _ { 3 } = \bar { \omega } _ { g } \ln \sqrt { \frac { \bar { \omega } _ { g } } { \pi } } \bar { I } _ { g }$$

$$c _ { 4 } = \frac { 1 } { 2 }$$

$$c _ { 5 } = \ln \sqrt { \frac { \bar { \omega } _ { g } } { \pi } } \bar { I } _ { g }$$

$$c _ { 6 } = \bar { \omega } _ { g }$$

c 7

$$\dot { \iota } + c _ { 6 } \Delta T \Delta R + c \, \\ \text{ents} \, \\ \nolimits \sqrt { \frac { \varpi _ { g } } { \pi } } \dot { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \t \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \t    \t                                                                                                                                                                                                         }                                                                                                                                                                                                        -1                                                                                                                                                                                                        +1                                                                                                                                                                                                       -2                                                                                                                                                                                                       
                                                                                                                                                                                                       .                                                                                                                                                                                                        |                                                                                                                                                                                                        \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde } { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \that { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \t                                                                                                                                                                                                      -1
                                                                                                                                                                                                       -0                                                                                                                                                                                                        "                                                                                                                                                                                                        )                                                                                                                                                                                                        <1$$

Let us introduce the following form for the expansion of a quantity X

$$X = \sum _ { n = 1 } ^ { \infty } X _ { n } ( \beta ) y _ { 0 } ^ { n }$$

where X n is a polynomial of order n -1 in β . Note that this form is stable by multiplication by a scalar, addition, multiplication, and exponentiation. We have already seen that ∆ ω g admits such a representation. From the relation

$$y = y _ { 0 } e ^ { - \beta \Delta \omega _ { 9 } }$$

we see that it is also the case for y . The quantities ∆ T and ∆ R can be written as a function of the variable y only. By expanding these quantities in terms of y we thus find that they also both admit this representation. From (42) it is also true for S . Finally, the partition function, Eqs.(40,41), can be written as

$$Z = e ^ { - \beta E _ { 0 } } \left [ 1 + \sum _ { n = 1 } ^ { \infty } P _ { n } ( \beta ) y _ { 0 } ^ { n } \right ] \quad \quad ( 4 5 )$$

where P n ( β ) is a polynomial of order n in β . This is the form given in Eq.(19).

To get the excited-state energies the quantities P n 1 need to be evaluated, see Eq.(22). Using the formula

$$P _ { n 1 } = \frac { 1 } { n! } \frac { \partial ^ { n + 1 } P ( \beta ) } { \partial \beta \partial y _ { 0 } ^ { n } } ( \beta = 0, y _ { 0 } = 0 )$$

we get

$$P _ { n 1 } = S _ { n 0 } + \Delta R _ { n 1 } + \sum _ { k + l = n \ k \geq 1 } S _ { k 0 } \Delta R _ { l 0 }. \quad ( 4 6 )$$

where S nk and ∆ nk are the coefficients of the n -th order polynomial in the representation, Eq.(43), for S and ∆ R , respectively. Using Eq.(42) and Q n 0 = ω n we have

$$S _ { n 0 } = c _ { 1 } \omega _ { n } + c _ { 2 } \Delta T _ { n 0 } + c _ { 3 } \Delta R _ { n 0 }$$

$$+ c _ { 4 } \sum _ { k _ { 1 } + k _ { 2 } = n } \omega _ { k _ { 1 } } \Delta T _ { k _ { 2 } 0 } + c _ { 5 } \sum _ { k _ { 1 } + k _ { 2 } = n } \omega _ { k 1 } \Delta R _ { k _ { 2 } 0 }$$

$$+ c _ { 6 } \sum _ { k _ { 1 } + k _ { 2 } = n } \Delta T _ { k _ { 1 } 0 } \Delta R _ { k _ { 2 } 0 }$$

$$+ c _ { 7 } \sum _ { k _ { 1 } + k _ { 2 } = n } \left ( \sum _ { k _ { 3 } + k _ { 4 } = k _ { 1 } } \omega _ { k _ { 3 } } \Delta T _ { k _ { 4 } 0 } \right ) \Delta R _ { k _ { 2 } 0 } \quad ( 4 7 )$$

where ∆ T nk are the coefficients of the polynomial for ∆ . T In each sum above, k i ≥ 1. Let us evaluate the

quantities ∆ R n 0 , ∆ R n 1 , and ∆ T n 0 .

Starting from

$$\Delta R = \sum _ { n = 1 } ^ { \infty } y _ { 0 } ^ { n } e ^ { - \beta n \Delta \omega _ { g } }$$

$$\text{where} \, \bar { \tau } _ { c } = \tau _ { c } ( \infty ) = \frac { 2 } { \bar { \omega } _ { g } } \text{ and}$$

$$\bar { I } _ { n } \, \equiv \int _ { - \infty } ^ { \infty } d x x ^ { n } e ^ { - \bar { \tau } _ { c } V ( x ) }.$$

After some algebra we have

$$\Delta R = \sum _ { n = 1 } ^ { \infty } y _ { 0 } ^ { n } \sum _ { k = 1 } ^ { \infty } \frac { ( - \beta n ) ^ { k } } { k! } \left ( \sum _ { l = 1 } ^ { \infty } Q _ { l } ( \beta ) y _ { 0 } ^ { l } \right ) ^ { k } = \sum _ { n = 1 } ^ { \infty } y _ { 0 } ^ { n } \sum _ { k = 1 } \alpha _ { k } y _ { 0 } ^ { k } \text{ with } ^ { \lambda }$$

$$\Delta \tau _ { c } = \frac { 2 } { \bar { \omega } _ { g } } \sum _ { n = 1 } ^ { \infty } ( - 1 ) ^ { n } \left [ \sum _ { k = 1 } ^ { \infty } X _ { k } ( \beta ) y _ { 0 } ^ { k } \right ] ^ { n }$$

with

$$\alpha _ { k } = \sum _ { l = 1 } ^ { k } \frac { ( - \beta n ) ^ { l } } { l! } \sum _ { k _ { 1 } + \dots + k _ { l } = k } Q _ { k _ { 1 } } ( \beta ) \dots Q _ { k _ { l } } ( \beta )$$

$$X _ { k } ( \beta ) = 2 \Delta R _ { k } ( \beta ) + \frac { Q _ { k } ( \beta ) } { \bar { \omega } _ { g } } + \frac { 2 } { \bar { \omega } _ { g } } \sum _ { l + m = k \, l \geq 1 } Q _ { l } ( \beta ) Q _ { m } ( \beta ).$$

We then have which gives

$$\Delta R _ { n 0 } = 1, \ \Delta R _ { 1 1 } = 0, \ \text{and} \ \Delta R _ { n 1 } = - \sum _ { k = 1 } ^ { n - 1 } k \, \omega _ { n - k } \ \text{for} \ n \geq 2 \ \text{with} \\ \intertext { f o t n e n o w e c o l u l t a } \Delta T _ { n } \ \Delta T _ { n } \ \text{and} \ \text{and} \ \text{where} \ \text{$\mathcal{E}$} \ \text{$\mathcal{E}$} } Y _ { n } = \frac { \cdot } { \bar { u } }$$

Let us now calculate ∆ T n 0 . ∆ , as defined by Eq.(38), T is decomposed as

Then

Thus

Finally, with

$$\Delta B _ { n } = \sum _ { k = 1 } ^ { n } \frac { ( - 1 ) ^ { k } } { k \bar { I } _ { g } ^ { k } } \sum _ { l _ { 1 } + \dots + l _ { k } = n } Z _ { l _ { 1 } } \dots Z _ { l _ { k } }.$$

The quantities ∆ B n 0 are then evaluated by taking β = 0 in the preceding expressions.

$$\Delta \tau _ { c } = \sum _ { n = 1 } ^ { \infty } Y _ { n } ( \beta ) y _ { 0 } ^ { n }$$

$$Y _ { n } = \frac { 2 } { \bar { \omega } _ { g } } \sum _ { k = 1 } ^ { n } ( - 1 ) ^ { k } \sum _ { l _ { 1 } + \dots + l _ { k } = n } X _ { l _ { 1 } } \dots X _ { l _ { k } }$$

$$I _ { g } = \bar { I } _ { g } + \sum _ { n = 1 } ^ { \infty } \frac { ( - 1 ) ^ { n } \bar { I } _ { n } } { n! } \left ( \sum _ { k = 1 } ^ { \infty } Y _ { k } ( \beta ) y _ { 0 } ^ { k } \right ) ^ { n }$$

$$I _ { g } = \bar { I } _ { g } + \sum _ { n = 1 } ^ { \infty } Z _ { n } ( \beta ) y _ { 0 } ^ { n }$$

$$Z _ { n } = \sum _ { k = 1 } ^ { n } ( - 1 ) ^ { k } \frac { \bar { I } _ { k } } { k! } \sum _ { l _ { 1 } + \dots + l _ { k } = n } Y _ { l _ { 1 } } \dots Y _ { l _ { k } }$$

$$\ln I _ { g } = \ln \bar { I } _ { g } + \ln \left ( 1 + \sum _ { n = 1 } ^ { \infty } \frac { Z _ { n } } { \bar { I } _ { g } } y _ { 0 } ^ { n } \right ).$$

$$\ln I _ { g } = \ln \bar { I } _ { g } + \sum _ { n = 1 } ^ { \infty } \Delta B _ { n } ( \beta ) y _ { 0 } ^ { n }$$

$$\Delta T = \Delta A + \Delta B$$

with and

$$\Delta B = \ln I _ { g } ( \beta ) - \ln \bar { I } _ { g }$$

Let us begin with ∆ A which can be written as

$$\Delta A = \frac { 1 } { 2 } \ln \left [ \frac { 1 + \sum _ { n = 1 } ^ { \infty } \frac { Q _ { n } ( \beta ) } { \bar { \omega } _ { g } } y _ { 0 } ^ { n } } { 1 + \sum _ { n = 1 } ^ { \infty } 2 \Delta R _ { n } ( \beta ) y _ { 0 } ^ { n } } \right ].$$

Expanding the two logarithmic terms, we get after some algebra

$$\Delta A _ { n 0 } = \frac { 1 } { 2 } \sum _ { k = 1 } ^ { n } \frac { ( - 1 ) ^ { k - 1 } } { k } \sum _ { l _ { 1 } + \dots + l _ { k } = n } \left [ \left ( \frac { \omega _ { l _ { 1 } } } { \bar { \omega } _ { g } } \right ) \dots \left ( \frac { \omega _ { l _ { k } } } { \bar { \omega } _ { g } } \right ) - 2 ^ { k } \right ] ^ { \cdot \cdot }$$

We are now left wit the calculation of ∆ B n 0 . We have

$$I _ { g } ( \beta ) = \sum _ { n = 0 } ^ { \infty } \frac { ( - 1 ) ^ { n } \bar { I } _ { n } } { n! } \Delta \tau _ { c } ^ { n }$$

with

$$\Delta A = \ln \sqrt { \frac { \omega _ { g } ( \beta ) } { \pi \coth \frac { \beta \omega _ { g } [ \tau _ { c } ( \beta ) ] } { 2 } } } - \ln \sqrt { \frac { \bar { \omega } _ { g } } { \pi } }$$

$$\Delta \tau _ { c } \equiv \tau _ { c } - \bar { \tau } _ { c }$$

Figure 1 shows the ground-state energy and the first eight excited state energies of the quartic oscillator with ω = 1 versus the coupling constant g . The 'exact' numerical energies, obtained by diagonalization of the Hamiltonian in a sufficiently large Gaussian basis set, are shown as solid lines. At the scale of the figure, the energy levels and their overall behavior as functions of g and n are well reproduced and appear satisfactory. Quantitative results are presented in Table I, which reports the relative errors ϵ in the computed energies. We present results for three different potentials: the first two are the quartic ( m = 2, ω = 1) and sextic ( m = 3, ω = 1) anharmonic potentials, respectively. The third is the 'pure' quartic potential, defined by the absence of a quadratic term ( ω = 0). Results are given for five values of g ( g = 0 1 1 10 40, and 400), spanning the weak- to . , , , strong-coupling regimes. A first remark is that relative errors are all negative. In other words, the computed energies are always smaller than the exact ones. Unfortunately, we were not able to understand the origin of this interesting observed property. A second remark is that relative errors are all of the order of a few percent (up to 10% in the worst case of the sextic potential at large g ). We consider this level of accuracy as quite satisfactory in view of the simplicity of the model. Let us insist on the fact that the model partition function has a particularly simple and compact form, Eq.(10) with Eqs.(11) and (12); and, most importantly, no adjustable parameters. In addition, the computational cost for calculating the energies is negligible. Another remark is that for each potential the relative error is maximal for the two lowest energies. Quite remarkably, the errors on the higher energies ( n ≥ 2) are almost constant and, also, nearly independent on the value of g . This is a quite interesting feature of the model. The comparison between the three different potentials is instructive. As expected, the accuracy reached for the sextic oscillator is inferior to that obtained for the quartic oscillator, a consequence of the greater anharmonicity of the former potential. This is also the case at small values of g for the pure quartic oscillator in which the quadratic term has been removed. However, for large g 's the errors become nearly identical. A striking feature of the energies of the pure quartic potential is that the relative errors are independent on g (actually, all digits of the errors are identical for the different g 's, a result not shown here). This result is explained as follows. By a simple rescaling of the Schr¨dinger equation, the exact energies o E n ( g ) of the pure quartic oscillator can be shown to scale as g 1 3 . Quite satisfactorily, it is also the case for our model. This can be proved by noting that ¯ ω g , Eq.(28), and the ω k 's, Eq.(27) scale also as g 1 3 and, then, by invoking the series of equations, (22),(46),(47) and (48). The exact and model energies having the same scaling in g , the relative errors are thus independent on g .

Figure 1. Energy levels of the quartic oscillator as a function of g for n = 0 to n = 8. Exact results given by the solid lines. Harmonic frequency, ω = 1

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01146v2-figure-017dd1e4-005f-4a8b-8674-8f9ddfc89583.jpg)

## ACKNOWLEDGMENTS

I would like to thank the Centre National de la Recherche Scientifique (CNRS) for its continued support. I also acknowledge funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No. 863481).

Table I. Relative errors ϵ (in %) on the energy levels E n for different values of g and potentials V [ a,n ] ( x ) = ax 2 + gx n ; ϵ = ( E m odel -E ex ) /E e x .

$$\text{Table I. Relative errors} \, e (in %) on the energy levels E _ { n } \, \text{for} \\ \text{different values of} \, g \, \text{ and potentials} \, V _ { [ a, n ] } ( x ) = a x ^ { 2 } + g x ^ { \prime } ; \\ \epsilon = ( E _ { \text{model} } - E _ { e x } ) / E _ { e x }. \\ \text{$E_{0}$} \, = E _ { 1 } \, = E _ { 2 } \, = E _ { 3 } \, = E _ { 4 } \, = E _ { 5 } \, = E _ { 6 } \, = E _ { 7 } \, = E _ { 8 } \\ g = 0. 1 \\ V _ { [ \frac { 3 } { 2 }, 4 ] } - 0. 5 \, - 0. 4 \, - 0. 5 \, - 0. 6 \, - 0. 7 \, - 0. 8 \, - 0. 8 \, - 0. 9 \, - 0. 9 \\ V _ { [ \frac { 3 } { 2 }, 6 ] } - 2. 5 \, - 2. 6 \, - 2. 7 \, - 2. 9 \, - 3. 1 \, - 3. 2 \, - 3. 3 \, - 3. 4 \, - 3. 4 \\ V _ { [ 0, 4 ] } - 4. 2 \, - 2. 0 \, - 1. 6 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \quad \text{$\theta$} \\ V _ { [ \frac { 3 } { 2 }, 4 ] } - 3. 7 \, - 1. 9 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \, - 1. 6 \quad \text{$\theta$} \\ V _ { [ \frac { 3 } { 2 }, 4 ] } - 4. 0 \, - 2. 0 \, - 1. 6 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \, - 1. 7 \quad \text{where} \\ \text{Appendix A: Formula for evaluating} \, \omega ^ { \prime } s$$

## Appendix A: Formula for evaluating ω n 's

In this appendix, we give the explicit formulas for the coefficients of the expansion of ω g in powers of y written as

$$\omega _ { g } = \bar { \omega } _ { g } + \sum _ { n = 1 } ^ { \infty } \omega _ { n } y ^ { n }. \quad \quad ( A 1 )$$

Note that to simplify the notation the dependence on m of the various quantities considered here will be omitted.

Let us define now the following quantities

$$B _ { n } \equiv \frac { d ^ { n } B ^ { ( m ) } } { d x ^ { n } } ( x _ { 0 } ) \quad \text{(A2)}$$

$$x _ { 0 } = \frac { 2 g } { \omega ^ { 2 m } } \bar { \omega } _ { g } ^ { m - 1 } \quad \quad \text{(A3)}$$

$$\alpha _ { n } = \frac { 1 } { n! } \frac { B _ { n } } { B _ { 0 } } x _ { 0 } ^ { n } \quad \quad \text{(A4)} \quad m$$

and

$$\beta _ { n } = \frac { \frac { 1 } { 2 } \left ( \frac { 1 } { 2 } - 1 \right ) \dots \left ( \frac { 1 } { 2 } - ( n - 1 ) \right ) } { n! } \quad \quad ( A 5 )$$

¯ ω g is given by

$$\bar { \omega } _ { g } = \omega \sqrt { B _ { 0 } } \quad \text{(A6)}$$

The first coefficient ω 1 is given by

$$\omega _ { 1 } = A \alpha _ { 1 } ( m - 1 ) \text{ \quad \ \ } ( A 7 )$$

with

$$A = \frac { \bar { \omega } _ { g } } { 1 - \frac { \alpha _ { 1 } ( m - 1 ) } { 2 } } \, \quad \quad ( A 8 )$$

For n ≥ 2 the ω n 's are function of the previous coefficients ω p with 1 ≤ p ≤ n -1 as follows

$$\omega _ { n } = A \left [ \alpha _ { 1 } ( m - 1 ) S _ { n } + \frac { 1 } { 2 } \alpha _ { 1 } T _ { n } + U _ { n } + V _ { n } \right ] \quad ( \mathbf A 9 )$$

where

$$S _ { n } = 1 + \sum _ { k + l = n \, k \geq 1 } \frac { \omega _ { k } } { \bar { \omega } _ { g } } \quad \quad ( A 1 0 )$$

$$T _ { n } = \sum _ { k = 2 } ^ { m - 1 } \binom { m - 1 } { k } \sum _ { l _ { 1 } + \dots + l _ { k } = n } a _ { l _ { i } \geq 1 } \dots a _ { l _ { k } } \quad ( \mathbf A 1 1 )$$

with

$$a _ { n } = \frac { \omega _ { n } } { \bar { \omega } _ { g } } + 2 S _ { n } \ n \geq 1 \quad \quad ( \text{A12} )$$

$$U _ { n } = \frac { 1 } { 2 } \sum _ { k = 2 } ^ { n } \alpha _ { k } \sum _ { l _ { 1 } + \dots + l _ { k } = n } b _ { l _ { i } \geq 1 } \dots b _ { l _ { k } } \quad \ ( A 1 3 )$$

$$b _ { n } = ( m - 1 ) a _ { n } + T _ { n } \ n \geq 1 \quad \ \ ( \text{A14} )$$

$$V _ { n } = \sum _ { k = 2 } ^ { n } \beta _ { k } \sum _ { l _ { 1 } + \dots + l _ { k } = n } c _ { l _ { i } } \dots c _ { l _ { k } } \quad \text{(A15)}$$

with and

with

$$c _ { n } = \alpha _ { 1 } b _ { n } + 2 U _ { n } \ n \geq 1 \quad \quad ( A 1 6 )$$

When no index fulfills the constraint in the sum, the corresponding quantity is equal to zero. Here, it means that S 1 = 1, T 1 = 0, U 1 = 0, V 1 = 0, and T n = 0 for m = 2.

- [1] A. Turbiner and J.C. Del Valle Rosale. Quantum Anharmonic Oscillator . World Scientific, 2023.
- [2] Alexander V. Turbiner. One-dimensional quasi-exactly solvable schr¨dinger equations. o Physics Reports , 642:171, 2016.
- [3] Daniel Baye. The lagrange-mesh method. Physics Reports , 565:1-107, 2015. The Lagrange-mesh method.
- [4] J. C. del Valle. Solving the one-dimensional timeindependent schr¨dinger o equation with high accuracy: The lagrangemesh mathematica ® package. International Journal of Modern Physics C , 35(01):2450011, 2024.
- [5] Carl M. Bender, Kaare Olaussen, and Paul S. Wang. Numerological analysis of the wkb approximation in large order. Phys. Rev. D , 16:1740-1748, Sep 1977.
- [6] R. N. Kesarwani and Y. P. Varshni. Eigenvalues of an anharmonic oscillator. Journal of Mathematical Physics , 22(9):1983-1989, 09 1981.
- [7] M. Lakshmanan, F. Karlsson, and P. O. Fr¨man. Phaseo integral calculation of the energy levels of a quantal anharmonic oscillator. Phys. Rev. D , 24:2586-2598, Nov 1981.
- [8] C.C. Yao and K.S. Cheng. A simple variational method for excited states of coupled anharmonic oscillators. Nuov Cim B , 111:645-652, 1996.
- [9] R. P. Feynman and H. Kleinert. Effective classical partition functions. Phys. Rev. A , 34:5080-5084, Dec 1986.
- [10] H. Kleinert. Improving the variational approach to path integrals. Physics Letters B , 280(3):251-260, 1992.
- [11] Ranjan Bhattacharya, Dhiranjan Roy, and Siddhartha Bhowmick. Simple systematics in the ground state energies of quantum anharmonic oscillators. Physics Letters A , 244(1):9-12, 1998.
- [12] Ananda Dasgupta, Dhiranjan Roy, and Ranjan Bhattacharya. Simple systematics in the energy eigenvalues of quantum anharmonic oscillators. Journal of Physics A: Mathematical and Theoretical , 40(4):773, jan 2007.
- [13] David Leonard and Paul Mansfield. Solving the anharmonic oscillator: tuning the boundary condition. Journal of Physics A: Mathematical and Theoretical , 40(33):10291, aug 2007.
- [14] M. Companys Franzke, A. Tichai, K. Hebeler, and A. Schwenk. Excited states from eigenvector continuation: The anharmonic oscillator. Physics Letters B , 830:137101, 2022.
- [15] M. Caffarel. Path integral for the quartic oscillator: An accurate analytic expression for the partition function, arxiv quant-ph arXiv:2312.09859v3., 2023.
- [16] R.P. Feynman and A.R. Hibbs. Quantum Mechanics and Paths Integrals . McGraw-Hill, New York, 1965.
- [17] Riccardo Giachetti and Valerio Tognetti. Variational approach to quantum statistical mechanics of nonlinear systems with application to sine-gordon chains. Phys. Rev. Lett. , 55:912-915, Aug 1985.
- [18] H. B¨ uttner and N. Flytzanis. Effective free energies. Phys. Rev. A , 36:3443-3445, Oct 1987.
- [19] L.S. Schulman. Techniques And Applications of Path Integration . Dover Publications Inc., 2005.
Michel Caffarel
Analytic Model for the Energy Spectrum of the Anharmonic Oscillator
In a recent work we have proposed an original analytic expression for the
partition function of the quartic oscillator. This partition function, which
has a simple and compact form with {\it no adjustable parameters}, reproduces
some key mathematical properties of the exact partition function and provides
free energies accurate to a few percent over a wide range of temperatures and
coupling constants. In this work, we present the derivation of the energy
spectrum of this model. We also generalize our previous study limited to the
quartic oscillator to the case of a general anharmonic oscillator. Numerical
application for a potential of the form $V(x)=\frac{\omega^2}{2} x^2 + g
x^{2m}$ show that the energy levels are obtained with a relative error of about
a few percent, a precision which we consider to be quite satisfactory given the
simplicity of the model, the absence of adjustable parameters, and the
negligible computational cost.