Schrödinger's equation in a hexagon (3D rendering)
Channel: Nils Berglund
Date: 13/06/2022

This simulation is a 3D rendering of a solution to Schrödinger's equation in a regular hexagon, analogous to the 2D simulations https://youtu.be/8WTIjRWjG1o , https://youtu.be/OJGTXVK3lpk and https://youtu.be/9Hnvd0ZAlPo . It shows the same evolution with two different representations:
Real part: 0:00
Phase: 1:02
In the first part, both the z-coordinate and the color hue represent the real part of the wave function. In the second part, the z-coordinate and the luminosity represent the modulus squared of the wave function, which gives the probability density of finding the particle at any given point. The color hue shows the phase of the wave function. In other words, if psi(t,x) denotes the complex wave function, the z-coordinate and luminosity depend in |psi(t,x)|^2, while the hue depends on the angle phi such that psi(t,x) = |psi(t,x)|*exp(i*phi(t,x)).  The initial state is a Gaussian wave packet directed at one corner of the hexagon. It is interesting to compare the quantum evolution near the corners to the classical one, shown for instance in the videos https://youtu.be/q1Rmjhp0UcQ and https://youtu.be/GkCYmRX8hCI

Render time: 50 minutes 53 seconds
Color scheme: Part 1 -  Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firing
https://github.com/BIDS/colormap
Part 2 - Turbo, by Anton Mikhailov
https://gist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a

Music: "Shadowing" by Corbyn Kites@corbynkites8006

See also https://images.math.cnrs.fr/Des-ondes-dans-mon-billard-partie-I.html for more explanations (in French) on  a few previous simulations of wave equations. 

The simulation solves the Schrödinger equation by discretization. The algorithm is adapted from the paper https://hplgit.github.io/fdm-book/doc/pub/wave/pdf/wave-4print.pdf 
C code: https://github.com/nilsberglund-orleans/YouTube-simulations
https://www.idpoisson.fr/berglund/software.html
Many thanks to my colleague Marco Mancini for helping me to accelerate my code!