It is well known that noisy forces can be used in different way for obtaining, following a first order stochastic differential equation, the Canonical distribution. For instance it is possible to use the Langevin dynamic defined by:

where is a random noise with variance and zero mean. It is easy to show, using the Fokker-Plank equation associated to this equation, that in order to obtain the usual Boltzmann distribution the matrix has to be chosen as . The problem to obtain the desired canonical distribution may be therefore solved in this way. In QMC one can calculate the covariance matrix and then invert this matrix to obtain and continue the dynamics. This method is unfortunately very unstable, because the matrix can be ill-defined because of statistical fluctuations. Moreover it is not possible to estimate the error on the temperature simulated.

Here we present a new method that uses these QMC forces to perform a Molecular Dynamics at finite temperature. In the past the major problem of using QMC to perform

In our simulation there exists a correlated noise associated to the forces. We rely on the central limit theorem implying the noise in all component of the forces evaluated by QMC is Gaussian with a given covariance matrix. We used the Jackknife re-sampling method (see Appendix A) to estimate the covariance matrix.

The idea of the GLQ is to use this noise to produce a given finite temperature using a Generalized Langevin Equation. The use of the Generalized Langevin Equation (GLE) as thermostat is not new. In the past some authors have used this approach to simulate different systems. This method was applied for the fist time by Schneider and Stoll (87), to study distortive-phase transitions. Later the GLE was used to simulate different systems and also to stabilize the usual Molecular Dynamic method (88).