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.
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 Demosupdated 30 Oct 2011