Quantum Monte Carlo simulation

Inserting equation (6-14) into equation (6-12) retrieves the p4/3 dependence of the exchange energy indicated in equation (3-5). This exchange functional is frequently called Slater exchange and is abbreviated by S. No such explicit expression is known for the correlation part, ec. However, highly accurate numerical quantum Monte-Carlo simulations of the homogeneous electron gas are available from the work of Ceperly and Alder, 1980. 

Marchi, M. Sprik, M. Klein, M. L., Calculation of the free energy of electron solvation in liquid ammonia using a path integral quantum Monte Carlo simulation, J. Phys. Chem. 1988, 92, 3625-3629... 

A different approach simulates the thermodynamic parameters of a finite spin system by using Monte Carlo statistics. Both classical spin and quantum spin systems of very large dimension can be simulated, and Monte Carlo many-body simulations are especially suited to fit a spin ensemble with defined interaction energies to match experimental data. In the case of classical spins, the simulations involve solving the equations of motion governing the orientations of the individual unit vectors, coupled to a heat reservoir, that take the form of coupled deterministic nonlinear differential equations.23 Quantum Monte Carlo involves the direct representation of many-body effects in a wavefunction. Note that quantum Monte Carlo simulations are inherently limited in that spin-frustrated systems can only be described at high temperatures.24... 

The accuracy of diffusion quantum Monte Carlo simulations in the determination of molecular equilibrium structures. 

Fig. 9 Temperature dependence of the directional order parameter obtained by the quantum Monte-Carlo simulation [4]...

In two dimensions, the ground-state phase diagram has been obtained by quantum Monte Carlo simulation [118,119]. There is a little difference in the vicinity of the Heisenberg point (/ = J ) from that in one dimension, but they are qualitatively quite similar. 

In this chapter we will review recent developments in the simulation of lattice (and continuum) models by classical and quantum Monte Carlo simulations. Unbiased numerical methods are required to obtain reliable results for classical and quantum lattice model when interactions or fluctuations are strong, especially in the vicinity of phase transitions, in frustrated models and in systems where quantum effects are important. For classical systems, molecular dynamics or the Monte Carlo method are the methods of choice since they can treat large systems. 

To perform a quantum Monte Carlo simulation on the world line representation, update moves that are ergodic and fulfill detailed balance are required. The simplest types of moves are again local updates. Since magnetization conservation prohibits the breaking of world lines, the local updates need to move world lines instead of just changing local states as in a classical model. 

In this section we will present generalizations of extended ensemble simulations to world line quantum Monte Carlo simulations, in particular ... 

M. Troyer and U.-J. Wiese (2005) Computational complexity and fundamental limitations to fermionic quantum Monte Carlo simulations. Phys. Rev. Lett. 94, p. 170201... 

A. Sandvik and J. Kurkijarvi (1991) Quantum Monte Carlo simulation method for spin systems. Phys. Rev. B 43, p. 5950... 

K. Delaney, C. Pierleoni and D.M. Ceperley (2006) Quantum Monte Carlo Simulation of the High-Pressure Molecular-Atomic Transition in Fluid Hydrogen. cond-mat/0603750, submitted to Phys. Rev. Letts.,... 

Solvent effects were found to have minimal influence on the excitation energies of phenol in aqueous solution using a quantum Monte Carlo simulation , which is in line with experimental observations on its absorption spectra " . Reaction field calculations of the excitation energy also showed a small shift in a solution continuum, in qualitative agreement with fluorescent studies of clusters of phenol with increasing number of water molecules " . The largest fluorescent shift of 2100 cm was observed in cyclohexane. 

Quantum Monte Carlo. - Mella et al. analysed the electron densities of both LiH and [LiH,e+] sampled during Quantum Monte Carlo simulations and saw an almost complete charge transfer from the Li atom basin, which bears a net positive charge, to the H atomic basin. 

Foulkes WMC et al (2001) Quantum Monte Carlo simulations of solids (review article). Rev Mod Phys 73 33... 

Measurements of the water proton mean kinetic energy and of the proton momentum distribution provide a richness of information about the potential surface that the proton experiences [3], including the effects of hydrogen bonding, thus complementing microscopic structural studies [7-9] and allowing a direct comparison with quantum Monte Carlo simulations [10,11]. 

Sandvik performed quantum Monte Carlo simulations of the Heisenberg Hamiltonian on the critical infinite percolation cluster (p = pp) using the stochastic series expansion method with operator loop update. 

To provide quantitative predictions about how to detect the superfluid-insulator transition in these experiments, Kashurnikov, Prokofev and Svistunov performed quantum Monte Carlo simulations of the singleparticle density matrix p,y = (I I /). They used the Bose-Hubbard model with harmonic confining potential and carried out world-line Monte Carlo simulations with the continuous-time Worm algorithm. The diagonal elements of the density matrix provide the real-space particle density, and... 

Despite this progress, the utility of quantum Monte Carlo simulations for studying quantum phase transitions in fermionic systems is still rather limited. Many of the most interesting problems, such as the ferromagnetic and antiferromagnetic quantum phase transitions in transition-metal compotmds and heavy-fermion materials, are still too complex to be attacked directly by microscopic quantum Monte Carlo methods. 

Katsufuji et al. used their results to estimate the critical value of the bandwidth-controlled Mott-Hubbard transition, (U/W)c 0.97, in the Ai-xCaxTiOs system. This value is below that estimated from quantum Monte Carlo simulations of the half-filled single band Hubbard model [72], (U/W)c 1.5 however a proper comparison between theory and experiment requires more careful consideration of both orbital degeneracy [149], which increases the critical value (U/W)c, and the deviation of the system from halffilling [150]. It is also worth noting for comparison that Ahn et al. estimate a value (U/W)c 2.7 from their measurements of the (Ca,Sr)Ru03 system [73]. 

As far as the correlation part is concerned, no such explicit expression for this term is known. However, there are highly accurate numerical quantum Monte-Carlo simulations of the uniform electron gas [21] from which several authors have derived analytical expressions by means of sophisticated interpolation schemes. One of the most widely used representations for this term is the one developed by Vosko, Wilk, and Nusair (VWN) [22]. 

G. B. Bachelet, D. M. Ceperley, and M. G. B. Chiocchetti, Phys. Rev. Lett., 62,2088 (1989). Novel Pseudo-Hamiltonian for Quantum Monte Carlo Simulations. 

M. Lao and P. A. Christiansen, /. Chem. Phys., 96, 2161 (1992). Relativistic Effective Potential Quantum Monte Carlo Simulations for Ne. 

G. Toth and G. Naray-Szabo, /. Chem. Phys., 100, 3742 (1994). Novel Semiempirical Method for Quantum Monte Carlo Simulation Application to Amorphous Silicon. 

M. Lewerenz and R. O. Watts, Mol. Rhys., 81,1075 (1994). Quantum Monte Carlo Simulation of Molecular Vibrations. Application to Formaldehyde. 

The layout is as follows the review begins with a brief recap of basis sets and programming strategy in the next two sections. Atom pairs are the physical entity used for integral evaluation, both in the Poisson equation technique and the Coulomb resolution. A case study of molecular interaction by Quantum Monte Carlo simulation, using ETOs is followed by accurate NMR chemical shift evaluation. 

---