Bohmian Quantum Trajectories for Coherent States of the Pöschl-Teller Potential

The trigonometric Pöschl–Teller potential (proposed for the first time in 1933 [1]) approximates the diatomic molecular potential energy. In the causal interpretation, superposed states, constant phase shifts, and non-factorizability are necessary for quantum motion to occur. Under certain circumstances, the motion might be chaotic (see Related Links). This Demonstration studies a two-dimensional version of the trigonometric Pöschl–Teller potential, which is in a special superposition state. If the perturbation term becomes zero (), the wave density (given by the square of the Schrödinger wavefunction) is not time dependent, but the corresponding autonomous differential equation system for the velocity field leads to a periodic motion in the configuration space. With the perturbation term, the Bohmian trajectory forms periodic, quasi-periodic, weak ergodic, or chaotic curves while interacting with the nodal points. Further investigations are necessary to capture the full dynamics of this system. The graphic shows the squared wavefunction, the trajectory (white), the nodal points (blue), the contour lines of the potential (gray), and the velocity field (red).


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


Associated Legendre polynomials arise as the solution of the Schrödinger equation:
with , , , , and so on. A degenerate, unnormalized wavefunction with time period for for the two-dimensional case can be expressed by:
where , are eigenfunctions, and are permuted eigenenergies of the corresponding stationary one-dimensional Schrödinger equation with . In the wavefunction , the perturbation term is given by . The eigenfunctions are defined by
where , are associated Legendre polynomials. The parameter is an arbitrary constant, is a constant phase shift (here ), and are the quantum numbers with and . The wavefunction is taken from [2].
For this Demonstration the three cases are distinguished by:
The velocity field is calculated from the gradient of the phase from the total wavefunction in the eikonal form (often called polar form) . For the three cases, the phase functions are:
For , the velocity field becomes time independent (autonomous).
In the program, if PlotPoints, AccuracyGoal, PrecisionGoal, and MaxSteps are increased, the results will be more accurate.
[1] G. Pöschl, E. Teller, "Bemerkungen zur Quantenmechanik des anharmonischen Oszillators," Zeitschrift für Physik, 83 (3–4), 1933 pp. 143–151, doi:10.1007/BF01331132.
[2] M. Trott, The Mathematica GuideBook for Symbolics, New York: Springer, 2006.
[3] "Bohmian-Mechanics.net." (Sep 26, 2014) www.bohmian-mechanics.net/index.html.
[4] S. Goldstein. "Bohmian Mechanics." The Stanford Encyclopedia of Philosophy. (Sep 26, 2014) plato.stanford.edu/entries/qm-bohm.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2017 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+