Mixing of Quantum Particles Influenced by Nodal Points in the Causal Interpretation

According to the causal interpretation of quantum theory, particles move in a deterministic orbit under the influence of the external potential plus a quantum potential generated by the wavefunction. These orbits can be examined with the tools of classical chaos theory. In this Demonstration a fractal basin boundary is found in an unnormalized, isotopic, two-dimensional harmonic oscillator with a constant phase shift, where the trajectories show chaos. Five moving nodal points (shown in gray) act like attractors protected by limit cycles or as a repeller, in which the quantum trajectories avoid approaching nodal points too closely. Further nodal points form saddle points with respect to the velocity, where the trajectories of the quantum particles become unstable and bifurcation of the trajectories occurs [1]. The attractors generate limited zones of configuration space to which the trajectories (originating from typical initial conditions) tend as time increases. The set of the initial conditions that converge to an attractor is called its basin; the boundary with fractal structures is called the fractal basin boundary or basin of attraction [2]. For each such nodal point, its basin of attraction depends on the set of initial conditions leading to longtime behavior of the trajectory that approaches that attractor. These boundary regions of the attracting nodal points have fractal structures. The graphics show the zones of attraction in relation to the starting points of the trajectories at given times and the trajectories in the two-dimensional configuration space. Each color defines an approximate attractor zone, to which the trajectories converge as time increases.


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


The wavefunction obeys the Schrödinger equation:
with , , and so on. The unnormalized wavefunction is:
with the eigenfunctions , where , are Hermite polynomials. The parameters , , are arbitrary constants (here: ), and , are constant phase shifts (), and the quantum numbers with . The velocity field is calculated from the gradient of the phase from the total wavefunction in the eikonal form .
Due to limitations of CPU power, AccuracyGoal, PrecisionGoal,and MaxSteps are decreased and the estimation of the nodal points are deactivated. So, be careful: the orbits have only qualitative character because of the inaccuracy of the numerically integrated functions.
[1] C. Efthymiopoulos, C. Kalapotharakos, and G. Contopoulos, "Origin of Chaos near Critical Points of the Quantum Flow," Physical Review E, 79, 2009 pp. 036203–1–036203–18. doi:10.1103/PhysRevE.79.036203, arXiv:0903.2655 [quant-ph].
[2] H-O. Peitgen and P. H. Richter, The Beauty of Fractals, Berlin: Springer-Verlag, 1986.
    • 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+