Quantum Mechanical Out-Of-Time-Ordered-Correlators for the Anharmonic (Quartic) Oscillator

Out-of-time-ordered correlators (OTOCs) have been suggested as a means to study quantum chaotic behavior in various systems. In this work, I calculate OTOCs for the quantum mechanical anharmonic oscillator with quartic potential, which is classically integrable and has a Poisson-like energy-level distribution. For low temperature, OTOCs are periodic in time, similar to results for the harmonic oscillator and the particle in a box. For high temperature, OTOCs exhibit a rapid (but power-like) rise at early times, followed by saturation consistent with $2\langle x^2\rangle_T \langle p^2\rangle_T$ at late times. At high temperature, the spectral form factor decreases at early times, bounces back and then reaches a plateau with strong fluctuations.


I. INTRODUCTION
The framework to calculate OTOCs for quantum mechanical systems with general Hamiltonians has been set up in Ref. [1]. The simplest case is that of the harmonic oscillator, which can be treated analytically, and gives OTOCs that are purely oscillatory. More complicated examples that exhibit classical chaos such as the two-dimensional stadium billiard must be treated numerically and give OTOCs that are growing non-exponentially at early times followed by a saturation at late times [1].
Early-time exponential growth of OTOCs has been found in a system of non-linearly coupled oscillators [2], which is expected to exhibit quantum chaos. More recently, however, it has been found that OTOCs show exponential growth in systems that do not possess quantum chaos, such as the inverted harmonic oscillator [3,4].
In the present work, I study OTOCs and spectral form factors for the case of the onedimensional anharmonic quantum oscillator, which has not been done before.

II. OUT-OF-TIME-ORDERED CORRELATORS IN QUANTUM MECHANICS
This section largely follows the setup outlined in Ref. [1] to calculate both microcanonical and thermal OTOCs in quantum mechanics, which I will be using in the following. Given a Hamiltonian H = H(x,p), the OTOC is defined as where the subscript T denotes calculation of the expectation value in a heat bath of temperature T = 1 β . Specifically, introducing the microcanonical OTOC c n (t) as the OTOC for a fixed energy eigenstate |n as in [1] c n (t) = − n|[x(t), p(0)] 2 |n , this implies where H|n = E n |n . The microcanonical OTOC may be expressed as and if the Hamiltonian is of the form it was shown in Ref. [1] that this leads to As a consequence, knowledge of the energy eigenvalues E n and the matrix elements x nm is sufficient to calculate the quantum mechanical OTOC.

Spectral form factor in quantum mechanics
Another quantity that is related to OTOCs is the spectral form factor. Following [5], it that has been proposed as a handle on information loss and I will employ its normalized version defined as where Z(c) for complex c is the analytically continued partition function defined in Eq. (3).

III. DISCRETE "SOLUTIONS" FOR THE QUARTIC OSCILLATOR
Let me consider a Hamiltonian such as (5) with an anharmonic oscillator potential V (x) given by with N ≥ 1 not necessarily integer.
Widely known examples of this potential are the cases N = 1 (the harmonic oscillator), N = 2 (the quartic oscillator), N = 3 (the sextic oscillator) and N = ∞ (particle in a box). Only N = 1 and N = ∞ have known analytic solutions while for generic N analytic solutions for the energy spectrum of the anharmonic quantum oscillator remain unknown. In the following, I concentrate on the case N = 2 (the quartic oscillator) for which the Schrödinger equation becomes with ψ n (x) = n|ψ . Now consider the auxiliary problem which has a discrete solution spectrum that is spanned by the associated generalized Laguerre polynomials. Specifically, separating wave-functions into parity-even (+) and parity-odd (-) solutions, (10) is solved by where i = 0, 1, 2, . . . and the respective eigenvalues γ in (10) take the values γ + i = 1 + 3i, γ − i = 2 + 3i. Expanding wave-functions ψ n (x) of the original problem (9) in terms of the auxiliary functions thus leads to with K → ∞ and obvious modifications for x < 0. I will refer to (12) as the "Laguerretransform" of ψ(x), and to c i as the "Laguerre-coefficients". Plugging (12) back into (9), the Laguerre coefficients c Multiplying these equations with e − x 3 , respectively, changing variables to z = x 3 3 and integrating using [6] ∞ 0 where (a) i denotes the Pochhammer symbol and 3 F 2 is a generalized hypergeometric function evaluated at unit argument, leads to where i, j = 0, 1, 2, . . . K. Wavefunctions are suitably normalized if ∞ −∞ ψ 2 n (x) = 1, which leads to the following normalization conditions forc i : A full solution for the quartic oscillator requires finding the eigenvalues and eigenvectors for K → ∞, for which A, B become infinite matrices. Keeping K finite can either be viewed as an approximation to the full solution, or, to a solution at a discrete set of points.
Specifically, solving the eigenproblem (18) for fixed and finite K, leads to K parity-even and K parity-odd wave-functionsψ n (x) that each are exact solutions to (9) at 2K points x = ±x 1 , ±x 2 . . . , ±x K . In many respects, this is reminiscent of constructing continuum solutions out of Fourier series.
For pedagogical reasons, some of the above properties are discussed for analytically tractable case K=1 in the next section.
I close this section by pointing out that because of parity-symmetry, matrix elements x nm require n to be even and m to be odd (or vice-versa). Without loss of generality assuming n to be even, one finds using (15), [7, Eq. 16.4.3] Other matrix elements that will be needed in the following are For K=1, the eigenvalue problems (18) is analytically solvable. The matrices A, B become which leads to the energy eigenvalues and eigenvectorsc Note that x 03 is an order of magnitude smaller than the other matrix elements, such that neglecting contributions involving x 03 is a good approximation. As a consequence, to good approximation, such that This results suggests that c 0 (t) is the superposition of two harmonics plus a constant. I find that (28) is correct semi-quantitatively also when K → ∞, when replacing the numerical values E 0 , E 1 , E 2 , x 01 , x 21 accordingly from K = 1 to higher K.

V. NUMERICAL RESULTS
Except for the ground state where reasonable approximations can be found 1 , Eq. (18) is hard to solve analytically n the limit K → ∞. However, for generic K, (18) is amenable to efficient numerical solution using readily available eigenvalue packages. In practice, I use only results that do not exhibit such sensitivity to the naked eye. As an example, results for the eigenvalue spectrum of the quartic oscillator are shown in Fig. 1. For comparison, results from the WKB approximation [8,9] are shown in Fig. 1, which in my units becomes It is a well-known feature that the eigenvalue spectrum of the quartic oscillator is numerically close to the WKB spectrum for all but lowest-lying eigenvalues (see inset in Fig. 1). However, the same is not true for the eigenfunctions ψ n (x) which in the WKB approximation feature a singularity close to the classical turning point.
Using the numerically determined Laguerre coefficientsc i , I find that the matrix elements x nm for generic K are strongly peaked for nearest energy-levels only, e.g.
Therefore, I find for n ≥ 2 c n (t) |b n,n−2 (t)| 2 + |b n,n (t)| 2 + |b n+2,n (t)| 2 , and b nm (t) similarly dominated by nearest-neighbor energy levels. As a consequence, evaluation of c n (t) is dominated by matrix elements x ij with i, j n, which greatly reduces the computational cost whenever K becomes large. In practice, I determine and use all matrix elements up to x 64,65 when calculating OTOCs, resulting in a very good approximation for c n (t) with n ≤ 32. Examples for low-lying microcanonical OTOCs c n (t) are shown in Fig. 2. First note that c 0 (t) qualitatively agrees with the analytic result (28) in that it behaves as the superposition of two harmonics plus a constant. For higher n, c n (t) exhibit a rapid rise at early times, followed by a qualitatively similar behavior: a constant plus two harmonics. However, results for c n (t) shown in Fig. 2 Since Fig. 2 also indicates that the amplitude of the high-frequency harmonic contributing to c n (t) is considerably less than those of the low-frequency parts for n 1, the resulting behavior of c n (t) is that of a rapid rise at early times followed by a near-constant behavior at late times.
The resulting OTOC (3) for various temperatures is shown in Fig. 3. Thermal averaging of the microcanonical c n (t) seems to have the effect of further reducing the amplitude of the harmonic contributions such that at temperatures T ≥ 8, Fig. 3 suggests a rapid earlytime rise followed by a plateau at late times for O T (t). While the early-time rise is clearly not a simple exponential, this qualitative behavior of O T (t) (rapid rise, saturation) has been associated with quantum chaotic behavior in systems that exhibit chaos, cf. Ref. [2].
However, for the case of the quartic oscillator, the system is classically integrable, and hence not expected to exhibit quantum chaos.
For quantum chaotic system, it is expected that at late times [2] lim where O T again denotes the thermal expectation value of the operator O. Using the quantum virial theorem, it is straightforward to relate the expectation value of the momentumoperator squared to whereas the expectation value of the position operator is given by (23). The result from evaluating (33) compared to O T (t) is shown for higher temperatures in Fig. 3. This comparison indicates that the apparent saturation of O T (t) at high temperatures and late times is roughly consistent with (33).
Finally, to complete the picture I also show the spectral form factor (7) in Fig. 4. At low temperatures, g(β, t) seems to first decrease and then bounce back to its value at early times, while for high temperatures, g(β, t) first drops, then rises again followed by a plateau with strong fluctuations. The high temperature behavior of g(β, t) is not that dissimilar from what has been observed for a single random matrix, cf. Ref. [5]. is very much reminiscent of that from random matrices, cf. Ref. [5].

VI. SUMMARY AND CONCLUSIONS
In this work, I studied the quantum-mechanical out-of-time-ordered-correlators for quartic interaction potential. It was found that at low temperature, OTOCs are periodic, while at high temperature OTOCs exhibit rapid power-law growth followed by an apparent saturation consistent with the temperature-dependent value 2 x 2 T p 2 T . At high temperatures, the spectral form factor for this theory decreases at early times, followed by a bounce and a plateau with strong fluctuations.
In conclusion, many interesting observables for one-dimensional quantum mechanics with anharmonic oscillator potential can be readily calculated using a spectral (Laguerre) decomposition. This may have implications for our understanding of quantum chaos, in particular when paired with the search for theories with gravitational duals.
While the present study was performed for quartic oscillator potential, all of the steps presented here can be repeated for potentials V (x) ∝ x N with arbitrary (even non-integer) N, cf. Eq. (8). Since it may be interesting in the future to study properties such as the power-law rise in early-time OTOCs or the dip-time in the spectral form factor as a function of the potential index N , I made the numerical codes used to generate the plots in this study publicly available at [10].
such that using i = 0 in Eq. (18), the even energy-levels E 2n can be written as where c where using (18) and (A1) the ground-state coefficients and energy can be found from the recursion relatioñ Specifically, the starting point isc iteration q i = δ i,0 , such that (A2) implies The first non-trivial iteration gives where here and in the following p F q denotes a generalized hypergeometric function of argument unity.
The next iterationc iteration 2 The generalized hypergeometric function appearing in this sum is a non-terminating