# Cellular automata rules for quantum mechanics

My limited understanding of quantum theory is that a quantum system is completely described by its wave function, which deterministically evolves according to Schrödinger's equation until wave function collapse.

Schrödinger's equation is essentially just a differential equation, which means (as far as I can see) that the state at every point in space changes depending only on its surroundings (i.e local rules). Correct me if I'm wrong but it's at the point of wave function collapse that the non-locality comes in, which is what cellular autonoma (CA) models can't deal with very well.

I'm not hypothesizing or suggesting that the universe IS a cellular autonoma (as Wolfram has suggested), but what would a CA model of Schrödinger's equation be like? (Ignoring wave function collapse for the minute)

Providing my understanding described above isn't wrong, I want to know, for the CA model:

1. What values describe the state at a given point? (i.e. what information would be in each CA cell?)
2. What rules are followed every discrete time tick to progress these values?

I ask because I've gained a good intuition (an "understanding" that I'm happy with at least!) of heat flow, wave propagation and the rules of fluid dynamics through interpreting their differential equations as local cellular automata rules, and writing little toy simulation programs to see how they behave and look. I'm a computer science/mathematics undergraduate and I've always said that I don't understand something until I can write a program that simulates it, so being a little too ambitious for my own good I want to write a simulation for a QM system! My background is I've read a few Penrose/Feynman books, and I've got the complete Feynman lectures on physics (not read them yet...), and physics has always been a side hobby, but I'm not too hot yet on the mathematical language that is used for QM.

-
What do you mean by a cellular autonoma model? Aren't you describing essentially numerical solution to a differential equation using some sort of finite difference method? –  user758556 Aug 10 '12 at 1:20
Probably he's referring to Gerard 't Hooft series of papers on cellular automata and quantum mechanics. –  linello Aug 16 '12 at 10:19
Quantum Hexodynamics: science20.com/hammock_physicist/quantum_hexodynamics –  Johannes Mar 12 at 14:47

There is nothing stopping you from interpreting the Schrodinger equation as rules for a cellular automata, in fact, the Schrodinger equation has the same form as the diffusion equation, but evolving in imaginary time. Let's write down some rules.

$$i \hbar \dot \psi = H \psi$$ assuming we have a time independent potential $$i \hbar \dot \psi = -\frac{\hbar^2}{2m} \psi'' + V \psi$$ Taking first order differences for the derivatives, we create

$$i \hbar \frac{1}{dt} \left( \psi^{t+1}_x - \psi^t_x \right) = -\frac{\hbar^2}{2m} \frac{1}{dx^2} \left( \psi^t_{x-1} - 2 \psi^t_{x} + \psi^t_{x+1} \right) + V \psi^t_x$$ rearranging we obtain $$\psi^{t+1}_x = (1 - i\bar V) \psi_x + i \alpha \left( \psi_{x-1} - 2 \psi_x + \psi_{x+1} \right)$$ where I have dropped the time index on the right, $\alpha = h dt / 2m dx^2$ and $\bar V = V dt / \hbar$, or perhaps more symbolically

$$\psi^{t+1} = \psi + i \alpha \ \text{lap}(\psi) - i \bar V \psi$$ where $\text{lap}$ stands for the lattice Laplacian operator (so that we can easily generalize to higher dimensions).

Compare this to the standard diffusion equation, for which we have $$\phi^{t+1} = \phi - D \ \text{lap}(\phi)$$

The Schrodinger equation appears very similar, although the field is complex and its change at each location is a $i$ rotation of its current state Laplacian and potential.

I will end with a cautionary note. I tried to implement this simple cellular automata, with hopes of adding a movie of the solution of a harmonic oscillator or something similar, but as presented, the automata is horribly ill conditioned and quickly blows up with a checkerboard mode. This is due to the fact that we have tried to describe the first order differential equation with a simple forward euler evolution, which doesn't have good convergence properties. By using RK4 (Runge Kutta 4th order), I was able to get the automata to behave nicely, but it wasn't all that inspiring of a movie, so I'll leave it as an exercise to the reader to implement this for themselves.

I will also note that this isn't really all that novel. Solving quantum systems by discretizing the Schrodinger equation is fairly common practice, as it is for most numerical solutions to differential equations generally, but people usually don't think or talk about these numerical solutions as cellular automata, although you could think of them as such if it aids in understanding.

-