Quantum Monte Carlo algorithm
The Algorithm for performing a variational Monte Carlo calculations runs thus as this
- Initialisation: Fix the number of Monte Carlo steps. Choose an initial R and variational parameters α and calculate |ψαT(R)|2.
- Initialise the energy and the variance and start the Monte Carlo calculation.
- Calculate a trial position Rp=R+r∗step where r is a random variable r∈[0,1].
- Metropolis algorithm to accept or reject this move w=P(Rp)/P(R).
- If the step is accepted, then we set R=Rp.
- Update averages
- Finish and compute final averages.
Observe that the jumping in space is governed by the variable step. This is Called brute-force sampling.
Need importance sampling to get more relevant sampling, see lectures below.