Quantum Monte-Carlo techniques

The approach to the evaluation of vibrational spectra described above is based on classical simulations for which quantum corrections are possible. The incorporation of quantum effects directly in simulations of large molecular systems is one of the most challenging areas in theoretical chemistry today. The development of quantum simulation methods is particularly important in the area of molecular spectroscopy for which quantum effects can be important and where the goal is to use simulations to help understand the structural and dynamical origins of changes in spectral lineshapes with environmental variables such as the temperature. The direct evaluation of quantum time- correlation functions for anharmonic systems is extremely difficult. Our initial approach to the evaluation of finite temperature anharmonic effects on vibrational lineshapes is derived from the fact that the moments of the vibrational lineshape spectrum can be expressed as functions of expectation values of positional and momentum operators. These expectation values can be evaluated using extremely efficient quantum Monte-Carlo techniques. The main points are summarized below. 

Quantum Monte Carlo techniques have considerable potential for application to problems involving open d or f shells where the treatment of electron correlation has proven particularly difficult. However if is to be a viable alternative one must be able to limit the simulations to small numbers of electrons and in addition relativeity must be included. Relativistic effective potentials offer one avenue (at the present time the only avenue) for achieving these conditions. However, as we have indicated, REPs do introduce carpi icat ions. 

In the investigation of Ortiz et al. [104], a stochastic method is presented which can handle complex Hermitian Hamiltonians where time-reversal invariance is broken explicitly. These workers fix the phase of the wave function and show that the equation for the modulus can be solved using quantum Monte Carlo techniques. Then, any choice for its phase affords a variational upper bound for the ground-state energy of the system. These authors apply this fixed phase method to the 2D electron fluid in an applied magnetic field with generalized periodic boundary conditions. 

Finally, can we dare to ask what is the future of first-principle MD It would be hard to be highly predictive. However we would like to quote the following directions of research QM/MM methods to treat quantum systems in an environment [92-94,225,226,269-272], Gaussian basis sets [23,30,38, 63,110,172] or Gaussian augmented plane waves methods [168] in search for order N methods [273,274] etc. Also, in order to go beyond Density Functional Theory, Quantum-Monte Carlo techniques are very attractive [119]. Some of these topics are already well-advanced and are discussed here in this book. 

There are several facts that speak in favor of the correctness of Eq. (124). Nooijen bases his reasoning [92] on the fact that the number of two-body coefficients is identical to the number of components of the Nakatsuji two-particle density equation [151], which is, in turn, equivalent to the time-independent Schrodinger equation for Hamiltonians containing up to two-body terms. The problem that remains open is the solubility of a rather complicated exponential variant of the Nakatsuji density equation, which forms an essential part of Nooijen s reasoning (cf. Eq. (11) in Ref. [92]). Van Voorhis and Head-Gordon base their reasoning [97] on the fact (exploited in Quantum Monte Carlo techniques) that one can always obtain the exact wave function by considering the expression... 

Table 5.2. Experimental vibrational redshifts for DF and HF with sequential addition of argon solvent atoms. Also shown are redshifts calculated using diffusion quantum Monte Carlo techniques from Ref. 66 and bound state variational calculations by Ernesti and Hutson from Refs. 9,11. The two columns reflect the values calculated within the approximation of pairwise additivity, and including the corrective three-body terms as described more fully in the text. 

Another branch of computational quantum mechanics, quantum Monte Carlo, is described in Chapter 3 by Professor James B. Anderson. Quantum Monte Carlo techniques, such as variational, diffusion, and Green s function, are explained, along with applications to atoms, molecules, clusters, liquids, and solids. Quantum Monte Carlo is not as widely used as other approaches to solving the Schrodinger equation for the electronic structure of a system, and the programs for running these calculations are not as user friendly as those based on the Hartree-Fock approach. This chapter sheds much needed light on the topic. 

PW91 (Perdew, Wang 1991) a gradient corrected DFT method QCI (quadratic conhguration interaction) a correlated ah initio method QMC (quantum Monte Carlo) an explicitly correlated ah initio method QM/MM a technique in which orbital-based calculations and molecular mechanics calculations are combined into one calculation QSAR (quantitative structure-activity relationship) a technique for computing chemical properties, particularly as applied to biological activity QSPR (quantitative structure-property relationship) a technique for computing chemical properties... 

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). 

The exchange-correlation functional for the uniform electron gas is known to high precision for all values of the electron density, n. For some regimes, these results were determined from careful quantum Monte Carlo calculations, a computationally intensive technique that can converge to the exact solution of the Schrodinger equation. Practical LDA functionals use a continuous function that accurately fits the known values of gas(/i). Several different... 

---