Clusters quantum Monte Carlo

Grossman J C and Mitas L 1995 Quantum Monte Carlo determination of electronic and structural properties of Si clusters (n 20) Phys. Rev. Lett. 74 1323... 

James B. Anderson, Quantum Monte Carlo Atoms, Molecules, Clusters, Liquids, and Solids. 

Static charge-density susceptibilities have been computed ab initio by Li et al (38). The frequency-dependent susceptibility x(r, r cd) can be calculated within density functional theory, using methods developed by Ando (39 Zang-will and Soven (40 Gross and Kohn (4I and van Gisbergen, Snijders, and Baerends (42). In ab initio work, x(r, r co) can be determined by use of time-dependent perturbation techniques, pseudo-state methods (43-49), quantum Monte Carlo calculations (50-52), or by explicit construction of the linear response function in coupled cluster theory (53). Then the imaginary-frequency susceptibility can be obtained by analytic continuation from the susceptibility at real frequencies, or by a direct replacement co ico, where possible (for example, in pseudo-state expressions). 

Two coupling modes are considered for the Pdj CO cluster the first mode (denoted as h) represents vibration of the rigid CO molecule with respect to the transition metal surface. The second mode is either the Pd-Pd vibration wi in the plane of Pd surface atoms (r) or out-of-plane stretch of the surface/sub-surface Pd-Pd bond (z). The total energy surfaces (h,r) and (h,z) are calculated for discrete points and then fitted to a fourth order polynomial. Variational and Quantum Monte Carlo (QMC) methods were subsequently applied to calculate the ground and first excited vibrational states of each two-dimensional potential surfaces. The results of the vibrational frequences (o using both the variational and QMC approach are displayed in Table II. 

Two methods are in common use for simulating molecular liquids the Monte Carlo method (MC) and molecular dynamics calculations (MD). Both depend on the availability of reasonably accurate potential energy surfaces and both are based on statistical classical mechanics, taking no account of quantum effects. In the past 10-15 years quantum Monte Carlo methods (QMC) have been developed that allow intramolecular degrees of freedom to be studied, but because of the computational complexity of this approach results have only been reported for water clusters. 

We examine the orbital compass model by utilizing the quantum Monte-Carlo method is a finite-size cluster [4], The simulations have been performed on a square lattice of Lx L sites with periodic-boundary conditions. 

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. 

D. M. Benoit, A. X. Chavagnac, and D. C. Clary, Speed improvement of diffusion quantum Monte Carlo calculations on weakly bound clusters, Chem. Phys. Lett. 283, 269-276 (1998). 

Many theoretical determinations have been proposed for small clusters. In Monte Carlo calculations, we showed the crucial role of the three-body interactions to describe clusters with N=2 and 3 (59hj). For ab initio quantum chemistry calculations, it is seen that small basis sets systematically overestimate the clustering energies (9, 52). In this case, the Zero-Point Energy correction is large (9), and seems less important with larger basis sets (63b). The role of the correlation contribution is not clear and seems to depend on the clusters considered (52, 63b, 69). The other determinations, from semi-empirical quantum chemical methods (8, 55) or analytical models (22, 50, 67a) are of variable accuracy. 

Figure 3.1 Size dependence of cohesive energies per atom (CE/n) of mercury clusters Hgn from calculations using a large-core EC-PP and CPP for Hg. Valence correlation is accounted for either within die hybrid model approach (HM) by a pair-potential adjusted for Hg2 or by pure-diffusion quantum Monte Carlo (PDMC) calculations (Wang etal 2000).

In this short review we have pointed out only very few of the basic issues involving the simulation of chemical systems with Quantum Monte Carlo. What has been achieved in the last few years is remarkable very precise calculations of small molecules, the most accurate calculations of the electron gas, silicon and carbon clusters, solids, and simulations of hydrogen at temperatures when bonds are forming. New methods have been developed as well high-accuracy trial wavefunctions for atoms, molecules, and solids, treatment of atomic cores, and the generalization of path-integral Monte Carlo to treat many-electron systems at positive temperatures. 

In this review, almost all of the simulations we have described use only classical mechanics to describe the nuclear motion of the reaction system. However, a more accurate analysis of many reactions, including some of the ones that have already been simulated via purely classical mechanics, will ultimately require some infusion of quantum mechanical methods. This infusion has already taken place in several different types of reaction dynamics electron transfer in solution, > i> 2 HI photodissociation in rare gas clusters and solids,i i 22 >2 ° I2 photodissociation in Ar fluid,and the dynamics of electron solvation.22-24 Since calculation of the quantum dynamics of a full solvent is at present too time-consuming, all of these calculations involve a quantum solute in a classical solvent. (For a system where the solvent is treated quantum mechanically, see the quantum Monte Carlo treatment of an electron transfer reaction in water by Bader et al. O) As more complex reaaions are investigated, the techniques used in these studies will need to be extended to take into account effects involving electron dynamics such as curve crossing, the interaction of multiple electronic surfaces and other breakdowns of the Born-Oppenheimer approximation, the effect of solvent and solute polarization, and ultimately the actual detailed dynamics of the time evolution of the electronic degrees of freedom. 

Anderson JB (1999) Quantum Monte Carlo atoms, molecules, clusters, liquids and solids. In lipkowitz KB, Boyd DB (eds) Reviews in computational chemistry. Wiley, New York... 

---