γ It delay fps plane: xy xt yt
GPU based Monte-Carlo simulation of the 4 Path Integrals for free particles (n = 1024). it is the number of MC iterations per frame. The thermostat algorithm and linear congruential random number generator are used. To simplify visualization paths are clamped to (0,1) interval. Really it is just closed random walks on torus.

draft

Path Integral for free particles

    W = Z -1 e-βE,     β = 1/T,
    E = ∑ i (xi+1 + xi )2/n.

In one MC pass we concider only even or odd (independent) beads. Let for xi bead
    xo = (xi+1 + xi-1 )/2,     Δ = (xi+1 - xi-1 )/2,
    xi+1 = xo + Δ,     xi-1 = xo - Δ.

Then
    ε = (xi - xi+1 )2 + (xi - xi-1 )2 = (xi - xo + Δ)2 + (xi - xo - Δ)2 = 2(δ2 + Δ2),
where δ = xi - xo and
    w ~ e-2βδ2.
We use random perturbation in every bead
    δ = γ(∑ k6 Rk - 3)
to generate random path.
WebCL Demos     updated 30 Oct 2011