Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review A

Paths to optimization in the multistate Rayleigh-Ritz variational method: Applications to the double-well quantum anharmonic oscillator

R. F. Bishop, M. F. Flynn, M. C. Boscá, R. Guardiola

Published Dec 1, 1989• 21 citations
Publisher

Abstract

Both single-well and multiwell one-dimensional anharmonic oscillators have provided an extremely fruitful testing ground for various modern techniques of microscopic quantum many-body theory and quantum field theory. They have served both as (0+1)-dimensional field-theory models of their important counterparts in a higher number of dimensions, and as algebraically simple yet highly nontrivial problems on which various approximation or truncation schemes could be tested to very high orders. Such techniques have included perturbation theory and its many possible resummation schemes, and the coupled-cluster method. In the present paper the spotlight is focused on variational methods, and we concentrate on their application to the particularly demanding double-well quartic anharmonic oscillator. We examine in detail both single-state and multistate variational calculations based on either orthogonal (uncorrelated) or nonorthogonal (correlated) basis functions. Particular attention is paid to the question of the sensitivity and accuracy of the results with respect to any nonlinear parameters that characterize the basis states. Practical schemes are suggested for their choice and are shown to be capable of very high accuracy. Use of this simple model also allows detailed comparisons to be made with other techniques of a perturbative or perturbative-variational hybrid character.

DOI
10.1103/physreva.40.6154

Related papers

Variational and coupled-cluster calculations of the spectra of anharmonic oscillators

R. F. Bishop, M. F. Flynn

Using the general anharmonic oscillator as a case study, we examine the coupled-cluster method (CCM) in some detail. Emphasis is specially placed on the accuracy of the standard ground-state energy calculation and the excitation spectrum as derived via the Emrich ansatz. We are particularly interested in problems that can arise from large differences between the exact wave function and its model counterpart that forms the starting point of the CCM. To this end we begin with a variational Hartree approximation applied to three particular anharmonic-oscillator Hamiltonians that contain, respectively, pure quartic, equally weighted cubic and quartic, and double-well perturbations. These are chosen to provide increasingly stringent tests for the CCM in the above sense. Considerable attention is paid to the variational description of the double-well case, where the various possible solutions to the Hartree equations are considered and a further calculation on the energy-level splitting is performed. The CCM is then used to improve systematically upon our chosen starting point and is shown, particularly for the ground state, both to produce extremely accurate results and to be rather resilient to even gross deficiencies in the starting wave function. The main problem is the double-well system, where the CCM shows its lack of intrinsic inbuilt tunneling mechanisms via the absence of level splitting. Even here the CCM still produces very accurate average energies for very deep wells. While clearly needing some modification in such extreme cases, the CCM is shown to be quite adaptable and robust with regard to inaccurate starting wave functions.

Physical Review A•Sep 1, 1988•53 citations

Rayleigh-Ritz Variational Method in The Complex Plane

M. W. AlMasri

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 $|α| < \tfrac{1}{2}$ for generalized Gaussian trial functions $ψ(z) = e^{αz^2 + β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 + λx^4$), adaptive Gaussian ansätze in position space yield a cubic stationarity equation and perturbative energy expansions beyond first order, capturing anharmonic wavefunction narrowing. In contrast, monomial trial functions ($ψ_n(z) = z^n$) in the Segal--Bargmann space -- while providing rigorous upper bounds $E_n = n + \tfrac{1}{2} + \tfrac{3λ}{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μ}{4} - \tfrac{9λ^2}{4} + \cdots$).

quant-pharXiv:2603.02257•Physica Scripta•Mar 3, 2026•0 citations

Variational method for approximating energy levels

Yu Zhou, Jay D. Mancini, P. F. Meier +1 more

We present here a systematic scheme for improving the variational wave functions and corresponding energy levels for quantum systems. By expanding the wave function around a variational parameter value, a family of independent functions can be systematically generated. The eigenstates are then obtained by diagonalizing the Hamiltonian matrix within the basis and optimized with respect to the variational parameter. As a test, the ground state of the quartic anharmonic oscillator has been investigated, and it is found that this scheme converges more rapidly than the conventional Lanczos method and yields better approximations of the energy levels than other variational methods. The effectiveness of this scheme for larger systems remains to be seen.

Physical Review A•Apr 1, 1995•15 citations

Approximation of the Quantum Anharmonic Oscillator: Investigating the Accuracy of Effective Models for the Double Well Potential

Linda Nordmark, Jacob Pirinen

The report investigates how the quantum mechanical model of a particle in one dimension underthe influence of the double well potential V (x) = λ(x^2 − x_0^2)^2 can be approximated using aneffective model, where the system is approximated with a system of truncated dimension. Theinfinite dimensional space of energy eigenfunctions to the double well potential was estimatedby two finite dimensional spaces, ones which were spanned by the ground- and first excitedeigenstates of the harmonic oscillators corresponding to the two minima of the potential. Theenergy eigenvalues given by these models were then studied. Additionally, the effects on thepartition function, the internal energy and the heat capacity of the system were checked tofurther study the accuracy of the models. To determine the accuracy of the approximation ahigh order, highly resolved finite element method was used for comparison as stand in for anexact, analytical solution. All investigated quantities were found to be better approximated forlarger λ and lower temperatures, with the effective model of larger dimension achieving a betterapproximation. The result with respect to λ likely results from that, for low energy, the doublewell potential resembles two harmonic oscillator potentials more the higher the energy barrier,i.e. the higher the λ. The effective models being better approximations for lower temperaturescan be explained by the fact that at low temperature, the lowest energies have the strongestimpact on the behaviour of the system.

•Jan 1, 2026•0 citations

Matrix mechanics as a practical tool in quantum theory: The anharmonic oscillator

Ching‐Teh Li, Abraham Klein, Franz R. Krejs

For an anharmonic oscillator, it follows from the number of zeros of the exact wave functions as functions of the quantum number $n$ that the matrix elements $〈n|x|{n}^{\ensuremath{'}}〉$ and $〈n|p|{n}^{\ensuremath{'}}〉$ should be rapidly decreasing functions of $|n\ensuremath{-}{n}^{\ensuremath{'}}|$. Matrix elements of polynomials in $x$ and $p$ should therefore be well approximated by a finite number of terms in their sum-rule decomposition. From the matrix elements of the equations of motion for $x$ and for $p$ and of the commutator $|x,p|$, one thereby obtains closed sets of nonlinear algebraic equations to characterize subspaces of the Hilbert space of exact eigenfunctions. The approximations are also derived from a novel variational principle, and numerous variant approximation schemes are suggested. Essentially exact numerical results are obtained and compared with previous work. The broad applicability of the techniques is emphasized.

Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Oct 15, 1975•13 citations

Anharmonic oscillator as a test of the coupled-cluster method

Ulrich B. Kaulfuss, M. Altenbokum

Recent applications of the coupled-cluster method to the (${\ensuremath{\varphi}}^{4}$${)}_{2}$ quantum field theory showed the necessity of a detailed analysis of its convergence behavior. The anharmonic oscillator has always served as a good test case for different approximation schemes. We treat the model both in the maximum-overlap condition and in the Hartree approximation. The ground-state energy is reproduced very well for all values of the coupling strength and already for a low-order truncation scheme. The expansion in correlation amplitudes can be carried out in extremely high order and shows a divergent tendency of the amplitudes. Introducing temperature dependence allows us to select the stable ground state out of a variety of solutions of the hierarchy of equations.

Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Jun 15, 1986•35 citations
Both single-well and multiwell one-dimensional anharmonic oscillators have provided an extremely fruitful testing ground for various modern techniques of microscopic quantum many-body theory and quantum field theory. They have served both as (0+1)-dimensional field-theory models of their important counterparts in a higher number of dimensions, and as algebraically simple yet highly nontrivial problems on which various approximation or truncation schemes could be tested to very high orders. Such techniques have included perturbation theory and its many possible resummation schemes, and the coupled-cluster method. In the present paper the spotlight is focused on variational methods, and we concentrate on their application to the particularly demanding double-well quartic anharmonic oscillator. We examine in detail both single-state and multistate variational calculations based on either orthogonal (uncorrelated) or nonorthogonal (correlated) basis functions. Particular attention is paid to the question of the sensitivity and accuracy of the results with respect to any nonlinear parameters that characterize the basis states. Practical schemes are suggested for their choice and are shown to be capable of very high accuracy. Use of this simple model also allows detailed comparisons to be made with other techniques of a perturbative or perturbative-variational hybrid character.
Using the general anharmonic oscillator as a case study, we examine the coupled-cluster method (CCM) in some detail. Emphasis is specially placed on the accuracy of the standard ground-state energy calculation and the excitation spectrum as derived via the Emrich ansatz. We are particularly interested in problems that can arise from large differences between the exact wave function and its model counterpart that forms the starting point of the CCM. To this end we begin with a variational Hartree approximation applied to three particular anharmonic-oscillator Hamiltonians that contain, respectively, pure quartic, equally weighted cubic and quartic, and double-well perturbations. These are chosen to provide increasingly stringent tests for the CCM in the above sense. Considerable attention is paid to the variational description of the double-well case, where the various possible solutions to the Hartree equations are considered and a further calculation on the energy-level splitting is performed. The CCM is then used to improve systematically upon our chosen starting point and is shown, particularly for the ground state, both to produce extremely accurate results and to be rather resilient to even gross deficiencies in the starting wave function. The main problem is the double-well system, where the CCM shows its lack of intrinsic inbuilt tunneling mechanisms via the absence of level splitting. Even here the CCM still produces very accurate average energies for very deep wells. While clearly needing some modification in such extreme cases, the CCM is shown to be quite adaptable and robust with regard to inaccurate starting wave functions.
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 ∣α∣<12|α| < \tfrac{1}{2}∣α∣<21​ for generalized Gaussian trial functions ψ(z)=eαz2+βzψ(z) = e^{αz^2 + βz}ψ(z)=eαz2+β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 (H^=−12∂x2+12x2+λx4\hat{H} = -\tfrac{1}{2}\partial_x^2 + \tfrac{1}{2}x^2 + λx^4H^=−21​∂x2​+21​x2+λx4), adaptive Gaussian ansätze in position space yield a cubic stationarity equation and perturbative energy expansions beyond first order, capturing anharmonic wavefunction narrowing. In contrast, monomial trial functions (ψn(z)=znψ_n(z) = z^nψn​(z)=zn) in the Segal--Bargmann space -- while providing rigorous upper bounds En=n+12+3λ4(2n2+2n+1)E_n = n + \tfrac{1}{2} + \tfrac{3λ}{4}(2n^2 + 2n + 1)En​=n+21​+43λ​(2n2+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., x3+x4x^3 + x^4x3+x4), showing that displacement parameters are essential to capture parity breaking and stabilization effects (E0≈12+3μ4−9λ24+⋯E_0 \approx \tfrac{1}{2} + \tfrac{3μ}{4} - \tfrac{9λ^2}{4} + \cdotsE0​≈21​+43μ​−49λ2​+⋯).
We present here a systematic scheme for improving the variational wave functions and corresponding energy levels for quantum systems. By expanding the wave function around a variational parameter value, a family of independent functions can be systematically generated. The eigenstates are then obtained by diagonalizing the Hamiltonian matrix within the basis and optimized with respect to the variational parameter. As a test, the ground state of the quartic anharmonic oscillator has been investigated, and it is found that this scheme converges more rapidly than the conventional Lanczos method and yields better approximations of the energy levels than other variational methods. The effectiveness of this scheme for larger systems remains to be seen.
The report investigates how the quantum mechanical model of a particle in one dimension underthe influence of the double well potential V (x) = λ(x^2 − x_0^2)^2 can be approximated using aneffective model, where the system is approximated with a system of truncated dimension. Theinfinite dimensional space of energy eigenfunctions to the double well potential was estimatedby two finite dimensional spaces, ones which were spanned by the ground- and first excitedeigenstates of the harmonic oscillators corresponding to the two minima of the potential. Theenergy eigenvalues given by these models were then studied. Additionally, the effects on thepartition function, the internal energy and the heat capacity of the system were checked tofurther study the accuracy of the models. To determine the accuracy of the approximation ahigh order, highly resolved finite element method was used for comparison as stand in for anexact, analytical solution. All investigated quantities were found to be better approximated forlarger λ and lower temperatures, with the effective model of larger dimension achieving a betterapproximation. The result with respect to λ likely results from that, for low energy, the doublewell potential resembles two harmonic oscillator potentials more the higher the energy barrier,i.e. the higher the λ. The effective models being better approximations for lower temperaturescan be explained by the fact that at low temperature, the lowest energies have the strongestimpact on the behaviour of the system.
For an anharmonic oscillator, it follows from the number of zeros of the exact wave functions as functions of the quantum number nnn that the matrix elements $〈n|x|{n}^{\ensuremath{'}}〉$ and $〈n|p|{n}^{\ensuremath{'}}〉$ should be rapidly decreasing functions of $|n\ensuremath{-}{n}^{\ensuremath{'}}|$. Matrix elements of polynomials in xxx and ppp should therefore be well approximated by a finite number of terms in their sum-rule decomposition. From the matrix elements of the equations of motion for xxx and for ppp and of the commutator ∣x,p∣|x,p|∣x,p∣, one thereby obtains closed sets of nonlinear algebraic equations to characterize subspaces of the Hilbert space of exact eigenfunctions. The approximations are also derived from a novel variational principle, and numerous variant approximation schemes are suggested. Essentially exact numerical results are obtained and compared with previous work. The broad applicability of the techniques is emphasized.
Recent applications of the coupled-cluster method to the (${\ensuremath{\varphi}}^{4}$)2{)}_{2})2​ quantum field theory showed the necessity of a detailed analysis of its convergence behavior. The anharmonic oscillator has always served as a good test case for different approximation schemes. We treat the model both in the maximum-overlap condition and in the Hartree approximation. The ground-state energy is reproduced very well for all values of the coupling strength and already for a low-order truncation scheme. The expansion in correlation amplitudes can be carried out in extremely high order and shows a divergent tendency of the amplitudes. Introducing temperature dependence allows us to select the stable ground state out of a variety of solutions of the hierarchy of equations.