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

arXiv:1201.1234 · Journal of Physics A Mathematical and Theoretical

WKB analysis of $\mathcal P$$\mathcal T$-symmetric Sturm–Liouville problems

Carl M. Bender, Hugh F Jones

Published Oct 23, 2012• 13 citations
arXivHTMLPDF

Abstract

Most studies of -symmetric quantum-mechanical Hamiltonians have considered the Schrödinger eigenvalue problem on an infinite domain. This paper examines the consequences of imposing the boundary conditions on a finite domain. As is the case with regular Hermitian Sturm–Liouville problems, the eigenvalues of the -symmetric Sturm–Liouville problem grow like n 2 for large n . However, the novelty is that a eigenvalue problem on a finite domain typically exhibits a sequence of critical points at which pairs of eigenvalues cease to be real and become complex conjugates of one another. For the potentials considered here this sequence of critical points is associated with a turning point on the imaginary axis in the complex plane. WKB analysis is used to calculate the asymptotic behaviours of the real eigenvalues and the locations of the critical points. The method turns out to be surprisingly accurate even at low energies. This article is part of a special issue of Journal of Physics A: Mathematical and Theoretical devoted to ‘Quantum physics with non-Hermitian operators’.

DOI
10.1088/1751-8113/45/44/444004

Related papers

Reflectionless potentials and symmetry

Zafar Ahmed, Carl M. Bender, Michael Berry

Large families of Hamiltonians that are non-Hermitian in the conventional sense have been found to have all eigenvalues real, a fact attributed to an unbroken symmetry. The corresponding quantum theories possess an unconventional scalar product. The eigenvalues are determined by differential equations with boundary conditions imposed in wedges in the complex plane. For a special class of such systems, it is possible to impose the -symmetric boundary conditions on the real axis, which lies on the edges of the wedges. The -symmetric spectrum can then be obtained by imposing the more transparent requirement that the potential be reflectionless.

arXiv:quant-ph/0508117•Journal of Physics A Mathematical and General•Sep 14, 2005•74 citations

Effects of Finite Boundaries on a One-Dimensional Harmonic Oscillator

Richard Vawter

In the study of the effects of finite boundaries on the magnetic properties of a solid, one encounters the problem of finding the energy eigenvalues of a one-dimensional linear harmonic oscillator located in a potential enclosure. The WKB method is employed to solve this problem for an oscillator whose center is fixed at some arbitrary position inside a potential box. Both numerical and analytical approximations of the WKB method are employed to find the energy eigenvalues over a wide range of the parameters of the problem. Numerical methods are also used on the exact series solutions of a bound oscillator to find the exact eigenvalues for the first few quantum states. A comparison of the two methods show that, in general, the WKB eigenvalues are accurate to much better than 1% except when the classical turning points are near the wall. Here the difference is of the order of 10% for the worst possible cases. The eigenvalues of a bound oscillator are shown to reduce to unbound-oscillator energy eigenvalues if the classical turning points are inside the potential enclosure and not near the walls. At the other extreme the eigenvalues are shown to become plane wave box eigenvalues when the separation of the classical turning points is large compared to the size of the enclosure. Also included is the relationship of the WKB solution to the Bohr-Sommerfeld quantization rule for the bound oscillator.

Physical Review•Oct 15, 1968•55 citations

Resonances from non-Hermitian quantum mechanical calculations

Nimrod Moiseyev

The Hermitian properties of the Hamiltonian are related not only to the operator itself but also to the functions on which it acts. Hermitian Hamiltonians operate on functions in the L 2 Hilbert space which correspond to boundary conditions which vanish at infinity. In this chapter, in order to move into the non-Hermitian domain, we will impose on the solutions to the time-independent Schrödinger equation (TISE) different boundary conditions which lead to solutions which can be associated with different types of the complex poles of the scattering matrix. These solutions will contain information which was not available within the scope of functions in L 2 . By imposing outgoing boundary conditions on the eigenfunctions of the timeindependent Hamiltonian complex eigenvalues, E res = ε − (i/2)Г, are obtained. These complex energies are associated with decaying resonance states which were discussed in the previous two chapters. The bound states (if they exist) appear as real eigenvalues since they result from exactly such outgoing boundary conditions which appear under the threshold energy. When incoming boundary conditions are imposed two kind of solution are obtained. One type of solution is the complex conjugates of the decay resonance solutions mentioned above. In scattering theory text books (see Taylor for example) the physical resonance solutions are associated with the poles of the scattering matrix which are embedded in the lower half of the complex energy plane. However, in nuclear physics the complex poles embedded in the upper half of the complex energy plane, so-called virtual states , are denoted as capture resonances .

Cambridge University Press eBooks•Feb 17, 2011•0 citations

PT-symmetric potentials having continuous spectra

Zichao Wen, Carl M. Bender

One-dimensional PT-symmetric quantum-mechanical Hamiltonians having continuous spectra are studied. The Hamiltonians considered have the form $H=p^2+V(x)$, where $V(x)$ is odd in $x$, pure imaginary, and vanishes as $|x|\to\infty$. Five PT-symmetric potentials are studied: the Scarf-II potential $V_1(x)=iA_1\,{\rm sech}(x)\tanh(x)$, which decays exponentially for large $|x|$; the rational potentials $V_2(x)=iA_2\,x/(1+x^4)$ and $V_3(x)=iA_3\,x/(1+|x|^3)$, which decay algebraically for large $|x|$; the step-function potential $V_4(x)=iA_4\,{\rm sgn}(x)θ(2.5-|x|)$, which has compact support; the regulated Coulomb potential $V_5(x)=iA_5\,x/(1+x^2)$, which decays slowly as $|x|\to\infty$ and may be viewed as a long-range potential. The real parameters $A_n$ measure the strengths of these potentials. Numerical techniques for solving the time-independent Schrödinger eigenvalue problems associated with these potentials reveal that the spectra of the corresponding Hamiltonians exhibit universal properties. In general, the eigenvalues are partly real and partly complex. The real eigenvalues form the continuous part of the spectrum and the complex eigenvalues form the discrete part of the spectrum. The real eigenvalues range continuously in value from $0$ to $+\infty$. The complex eigenvalues occur in discrete complex-conjugate pairs and for $V_n(x)$ ($1\leq n\leq4$) the number of these pairs is finite and increases as the value of the strength parameter $A_n$ increases. However, for $V_5(x)$ there is an {\it infinite} sequence of discrete eigenvalues with a limit point at the origin. This sequence is complex, but it is similar to the Balmer series for the hydrogen atom because it has inverse-square convergence.

arXiv:2002.04398•Feb 4, 2020•3 citations

A numerical and analytical investigation into non-Hermitian Hamiltonians

Gert Jermia Cornelus Wessels

In this thesis we aim to show that the Schrödinger equation, which is a boundary eigenvalue problem, can have a discrete and real energy spectrum (eigenvalues) even when the Hamiltonian is non-Hermitian.After a brief introduction into non-Hermiticity, we will focus on solving the Schrödinger equation with a special class of non-Hermitian Hamiltonians, namely PTsymmetric Hamiltonians.PT -symmetric Hamiltonians have been discussed by various authors [1, 2, 3, 4, 5] with some of them focusing specifically on obtaining the real and discrete energy spectrum.Various methods for solving this problematic Schrödinger equation will be considered.After starting with perturbation theory, we will move on to numerical methods.Three different categories of methods will be discussed.First there is the shooting method based on a Runge-Kutta solver.Next, we investigate various implementations of the spectral method.Finally, we will look at the Riccati-Padé method, which is a numerical implemented analytical method.PT -symmetric potentials need to be solved along a contour in the complex plane.We will propose modifications to the numerical methods to handle this.After solving the widely documented PT -symmetric Hamiltonian H = p 2 -(ix) N with these methods, we give a discussion and comparison of the obtained results.Finally, we solve another PT -symmetric potential, illustrating the use of paths in the complex plane to obtain a real and discrete spectrum and their influence on the results.

•Mar 1, 2009•4 citations

PT-symmetry and Schrödinger operators. The double well case

Nawal Mecherout, Naima Boussekkine, Thierry Ramond +1 more

We study a class of PT-symmetric semiclassical Schrödinger operators, which are perturbations of a selfadjoint one. Here, we treat the case where the unperturbed operator has a double-well potential. In the simple well case, two of the authors have proved in \cite{BoMe14} that, when the potential is analytic, the eigenvalues stay real for a perturbation of size O(1). We show here, in the double-well case, that the eigenvalues stay real only for exponentially small perturbations, then bifurcate into the complex domain when the perturbation increases and we get precise asymptotic expansions. The proof uses complex WKB-analysis, leading to a fairly explicit quantization condition.

arXiv:1502.06102•Feb 21, 2015•4 citations
WKB analysis of P\mathcal PPT\mathcal TT-symmetric Sturm–Liouville problems
Most studies of -symmetric quantum-mechanical Hamiltonians have considered the Schrödinger eigenvalue problem on an infinite domain. This paper examines the consequences of imposing the boundary conditions on a finite domain. As is the case with regular Hermitian Sturm–Liouville problems, the eigenvalues of the -symmetric Sturm–Liouville problem grow like n 2 for large n . However, the novelty is that a eigenvalue problem on a finite domain typically exhibits a sequence of critical points at which pairs of eigenvalues cease to be real and become complex conjugates of one another. For the potentials considered here this sequence of critical points is associated with a turning point on the imaginary axis in the complex plane. WKB analysis is used to calculate the asymptotic behaviours of the real eigenvalues and the locations of the critical points. The method turns out to be surprisingly accurate even at low energies. This article is part of a special issue of Journal of Physics A: Mathematical and Theoretical devoted to ‘Quantum physics with non-Hermitian operators’.
Large families of Hamiltonians that are non-Hermitian in the conventional sense have been found to have all eigenvalues real, a fact attributed to an unbroken symmetry. The corresponding quantum theories possess an unconventional scalar product. The eigenvalues are determined by differential equations with boundary conditions imposed in wedges in the complex plane. For a special class of such systems, it is possible to impose the -symmetric boundary conditions on the real axis, which lies on the edges of the wedges. The -symmetric spectrum can then be obtained by imposing the more transparent requirement that the potential be reflectionless.
In the study of the effects of finite boundaries on the magnetic properties of a solid, one encounters the problem of finding the energy eigenvalues of a one-dimensional linear harmonic oscillator located in a potential enclosure. The WKB method is employed to solve this problem for an oscillator whose center is fixed at some arbitrary position inside a potential box. Both numerical and analytical approximations of the WKB method are employed to find the energy eigenvalues over a wide range of the parameters of the problem. Numerical methods are also used on the exact series solutions of a bound oscillator to find the exact eigenvalues for the first few quantum states. A comparison of the two methods show that, in general, the WKB eigenvalues are accurate to much better than 1% except when the classical turning points are near the wall. Here the difference is of the order of 10% for the worst possible cases. The eigenvalues of a bound oscillator are shown to reduce to unbound-oscillator energy eigenvalues if the classical turning points are inside the potential enclosure and not near the walls. At the other extreme the eigenvalues are shown to become plane wave box eigenvalues when the separation of the classical turning points is large compared to the size of the enclosure. Also included is the relationship of the WKB solution to the Bohr-Sommerfeld quantization rule for the bound oscillator.
The Hermitian properties of the Hamiltonian are related not only to the operator itself but also to the functions on which it acts. Hermitian Hamiltonians operate on functions in the L 2 Hilbert space which correspond to boundary conditions which vanish at infinity. In this chapter, in order to move into the non-Hermitian domain, we will impose on the solutions to the time-independent Schrödinger equation (TISE) different boundary conditions which lead to solutions which can be associated with different types of the complex poles of the scattering matrix. These solutions will contain information which was not available within the scope of functions in L 2 . By imposing outgoing boundary conditions on the eigenfunctions of the timeindependent Hamiltonian complex eigenvalues, E res = ε − (i/2)Г, are obtained. These complex energies are associated with decaying resonance states which were discussed in the previous two chapters. The bound states (if they exist) appear as real eigenvalues since they result from exactly such outgoing boundary conditions which appear under the threshold energy. When incoming boundary conditions are imposed two kind of solution are obtained. One type of solution is the complex conjugates of the decay resonance solutions mentioned above. In scattering theory text books (see Taylor for example) the physical resonance solutions are associated with the poles of the scattering matrix which are embedded in the lower half of the complex energy plane. However, in nuclear physics the complex poles embedded in the upper half of the complex energy plane, so-called virtual states , are denoted as capture resonances .
One-dimensional PT-symmetric quantum-mechanical Hamiltonians having continuous spectra are studied. The Hamiltonians considered have the form H=p2+V(x)H=p^2+V(x)H=p2+V(x), where V(x)V(x)V(x) is odd in xxx, pure imaginary, and vanishes as ∣x∣→∞|x|\to\infty∣x∣→∞. Five PT-symmetric potentials are studied: the Scarf-II potential V1(x)=iA1 sech(x)tanh⁡(x)V_1(x)=iA_1\,{\rm sech}(x)\tanh(x)V1​(x)=iA1​sech(x)tanh(x), which decays exponentially for large ∣x∣|x|∣x∣; the rational potentials V2(x)=iA2 x/(1+x4)V_2(x)=iA_2\,x/(1+x^4)V2​(x)=iA2​x/(1+x4) and V3(x)=iA3 x/(1+∣x∣3)V_3(x)=iA_3\,x/(1+|x|^3)V3​(x)=iA3​x/(1+∣x∣3), which decay algebraically for large ∣x∣|x|∣x∣; the step-function potential V4(x)=iA4 sgn(x)θ(2.5−∣x∣)V_4(x)=iA_4\,{\rm sgn}(x)θ(2.5-|x|)V4​(x)=iA4​sgn(x)θ(2.5−∣x∣), which has compact support; the regulated Coulomb potential V5(x)=iA5 x/(1+x2)V_5(x)=iA_5\,x/(1+x^2)V5​(x)=iA5​x/(1+x2), which decays slowly as ∣x∣→∞|x|\to\infty∣x∣→∞ and may be viewed as a long-range potential. The real parameters AnA_nAn​ measure the strengths of these potentials. Numerical techniques for solving the time-independent Schrödinger eigenvalue problems associated with these potentials reveal that the spectra of the corresponding Hamiltonians exhibit universal properties. In general, the eigenvalues are partly real and partly complex. The real eigenvalues form the continuous part of the spectrum and the complex eigenvalues form the discrete part of the spectrum. The real eigenvalues range continuously in value from 000 to +∞+\infty+∞. The complex eigenvalues occur in discrete complex-conjugate pairs and for Vn(x)V_n(x)Vn​(x) (1≤n≤41\leq n\leq41≤n≤4) the number of these pairs is finite and increases as the value of the strength parameter AnA_nAn​ increases. However, for V5(x)V_5(x)V5​(x) there is an {\it infinite} sequence of discrete eigenvalues with a limit point at the origin. This sequence is complex, but it is similar to the Balmer series for the hydrogen atom because it has inverse-square convergence.
In this thesis we aim to show that the Schrödinger equation, which is a boundary eigenvalue problem, can have a discrete and real energy spectrum (eigenvalues) even when the Hamiltonian is non-Hermitian.After a brief introduction into non-Hermiticity, we will focus on solving the Schrödinger equation with a special class of non-Hermitian Hamiltonians, namely PTsymmetric Hamiltonians.PT -symmetric Hamiltonians have been discussed by various authors [1, 2, 3, 4, 5] with some of them focusing specifically on obtaining the real and discrete energy spectrum.Various methods for solving this problematic Schrödinger equation will be considered.After starting with perturbation theory, we will move on to numerical methods.Three different categories of methods will be discussed.First there is the shooting method based on a Runge-Kutta solver.Next, we investigate various implementations of the spectral method.Finally, we will look at the Riccati-Padé method, which is a numerical implemented analytical method.PT -symmetric potentials need to be solved along a contour in the complex plane.We will propose modifications to the numerical methods to handle this.After solving the widely documented PT -symmetric Hamiltonian H = p 2 -(ix) N with these methods, we give a discussion and comparison of the obtained results.Finally, we solve another PT -symmetric potential, illustrating the use of paths in the complex plane to obtain a real and discrete spectrum and their influence on the results.
We study a class of PT-symmetric semiclassical Schrödinger operators, which are perturbations of a selfadjoint one. Here, we treat the case where the unperturbed operator has a double-well potential. In the simple well case, two of the authors have proved in \cite{BoMe14} that, when the potential is analytic, the eigenvalues stay real for a perturbation of size O(1). We show here, in the double-well case, that the eigenvalues stay real only for exponentially small perturbations, then bifurcate into the complex domain when the perturbation increases and we get precise asymptotic expansions. The proof uses complex WKB-analysis, leading to a fairly explicit quantization condition.