Chaotic Quantum Motion of Two Particles in a 3D Harmonic Oscillator Potential

A system with three degrees of freedom, consisting of a superposition of three coherent stationary eigenfunctions with commensurate energy eigenvalues and a constant relative phase, can exhibit chaotic motion in the de Broglie–Bohm formulation of quantum mechanics (see Quantum Motion of Two Particles in a 3D Trigonometric Pöschl–Teller Potential). We consider here an analog using a three-dimensional harmonic-oscillator potential. In this case, the velocities of the particles are autonomous, with a complex, chaotic trajectory structure. Two particles are placed randomly, separated by an initial distance , on the boundary of the harmonic potential.
The dynamic behavior for such a system is quite complex. Some of the curves are closed and periodic, while others are quasi-periodic. In the region of nodal points of the wavefunction, the trajectories apparently become accelerated and chaotic. The parameters have to be chosen carefully, because of the singularities in the velocities and the resulting large oscillations, which can lead to very unstable trajectories. The motion originates from the relative phase of the total wavefunction, which has no analog in classical particle mechanics. Further investigation to capture the full dynamics of the system is necessary.
The graphics show three-dimensional contour plots of the squared wavefunction (if enabled) and two initially neighboring trajectories. Black points mark the initial positions of the two quantum particles and green points the actual positions. Blue points indicate the nodal point structure.


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


Associated Hermite polynomials arise as the solution of the Schrödinger equation:
with , , and so on. A degenerate, unnormalized, complex-valued wavefunction for the three-dimensional case can be given by:
where , , are eigenfunctions, and are permuted eigenenergies of the corresponding stationary one-dimensional Schrödinger equation with . The eigenfunctions are defined by
where , , are Hermite polynomials. The parameter is a constant phase shift. The eigenvalues' numbers depend on the three quantum numbers .
In this Demonstration, the wavefunction is defined by:
In this case, the square of the Schrödinger wavefunction , where is its complex conjugate, is not time dependent:
The velocity field is calculated from the gradient of the phase from the total wavefunction in the eikonal form (often called polar form) . The time-dependent phase function from the total wavefunction is:
The corresponding velocity field becomes time independent (autonomous) because of the time-independent gradient of the phase function.
In the program, if PlotPoints, AccuracyGoal, PrecisionGoal, MaxSteps, and MaxIterations are enabled, increasing them will give more accurate results.
[1] "Bohmian-Mechanics.net." (Jul 30, 2015) www.bohmian-mechanics.net/index.html.
[2] S. Goldstein. "Bohmian Mechanics." The Stanford Encyclopedia of Philosophy. (Jul 30, 2015) 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 © 2018 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+