Computer simulations Quantum Monte Carlo

Beyond the clusters, to microscopically model a reaction in solution, we need to include a very big number of solvent molecules in the system to represent the bulk. The problem stems from the fact that it is computationally impossible, with our current capabilities, to locate the transition state structure of the reaction on the complete quantum mechanical potential energy hypersurface, if all the degrees of freedom are explicitly included. Moreover, the effect of thermal statistical averaging should be incorporated. Then, classical mechanical computer simulation techniques (Monte Carlo or Molecular Dynamics) appear to be the most suitable procedures to attack the above problems. In short, and applied to the computer simulation of chemical reactions in solution, the Monte Carlo [18-21] technique is a numerical method in the frame of the classical Statistical Mechanics, which allows to generate a set of system configurations... 

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. 

Though the above argument can leave no doubt that in the jellium model there will be a localized assembly of electrons, i.e. a Wigner crystal, in the extremely low density limit, the actual analytic calculation of when the electron liquid, at absolute zero of temperature, freezes as the density is lowered has proved very delicate [20]. Eventually, this matter was settled using quantum Monte Carlo computer simulation by Ceperley and Alder [38], They found in this way that the crystallization first occurred at rs = 100. Herman and March [39] subsequently pointed out that, for the Wigner crystal phase, the theoretical expression [40,41]... 

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

The United States is currently a leader in most areas of theoretical/ computational chemistry. In basic theory, Europe has many talented young investigators. Within the next 10 years, given these demographics, the U.S. leadership will be challenged by Europe in electronic structure and basic theory development. This trend does not characterize the entire field of theoretical chemistry. For example, Monte Carlo and molecular dynamics simulation methods were invented in the United States, and to this day, the United States maintains a strong position, especially in quantum Monte Carlo calculations. 

Increases in computer power and improvements in algorithms have greatly extended the range of applicability of classical molecular simulation methods. In addition, the recent development of Internal Coordinate Quantum Monte Carlo (ICQMC) has allowed the direct comparison of classical simulations and quantum mechanical results for some systems. In particular, it has provided new insights into the zero point energy problem in many body systems. Classical studies of non-linear dynamics and chaos will be compared to ICQMC results for several systems of interest to nanotechnology applications. The ramifications of these studies for nanotechnology applications will be discussed. 

Because many details of the dynamics and structure of chemical systems cannot be directly observed, molecular simulation methods such as molecular dynamics (MD) [1-31, molecular mechanics (MM) [4], and classical and quantum Monte Carlo [5,6] are extremely valuable tools for making sense of experimental results. In the context of nanotechnology, molecular simulation is crucial for studying the feasibility of proposed directions of research and development [7], With the rapid improvement in computing power and algorithms, the capabilities and range of applicability of molecular simulation have dramatically increased over the past decade. 

Some of the most widely used computational approaches will be briefly described below, namely some quantum chemical methods, classical simulations by Monte Carlo and Molecular Dynamics techniques and a few mesoscale methods. 

Within the past decade much progress has also been made in experimental realizations of quantum computing hardware. Many architectures have been proposed based on a variety of physical hardware. On a small scale, quantum information has been stored and manipulated in superconducting quantum bits (qubits) [4,5], trapped ions [6,7], electron spins [8-11], nuclear spins in the liquid or solid state [12], and other systems. On the theoretical side, new quantum algorithms have recently been found, exhibiting significant pol momial speedups on quantum computers for solution of sparse linear equations or differential equations [13,14], quantum Monte Carlo problems [15], and classical simulated annealing problems [16]. 

This is a class of algorithms which makes feasible on contemporary computers an exact Monte Carlo solution of the Schrodinger equation. It is exact in the sense that as the number of steps of the random walk becomes large the computed energy tends toward the ground state energy of a finite system of bosons. It shares with all Monte Carlo calculations the problem of statistical errors and (sometimes) bias. In the simulations of extensive systems, in addition, there is the approximation of a uniform fluid by a finite portion with (say) periodic boundary conditions. The latter approximation appears to be less serious in quantum calculations than in corresponding classical ones. 

Progress in the theoretical description of reaction rates in solution of course correlates strongly with that in other theoretical disciplines, in particular those which have profited most from the enonnous advances in computing power such as quantum chemistry and equilibrium as well as non-equilibrium statistical mechanics of liquid solutions where Monte Carlo and molecular dynamics simulations in many cases have taken on the traditional role of experunents, as they allow the detailed investigation of the influence of intra- and intemiolecular potential parameters on the microscopic dynamics not accessible to measurements in the laboratory. No attempt, however, will be made here to address these areas in more than a cursory way, and the interested reader is referred to the corresponding chapters of the encyclopedia. 

Recently, molecular dynamics and Monte Carlo calculations with quantum mechanical energy computation methods have begun to appear in the literature. These are probably some of the most computationally intensive simulations being done in the world at this time. 

The computation of quantum many-body effects requires additional effort compared to classical cases. This holds in particular if strong collective phenomena such as phase transitions are considered. The path integral approach to critical phenomena allows the computation of collective phenomena at constant temperature — a condition which is preferred experimentally. Due to the link of path integrals to the partition function in statistical physics, methods from the latter — such as Monte Carlo simulation techniques — can be used for efficient computation of quantum effects. 

MD simulations in expHcit solvents are stiU beyond the scope of the current computational power for screening of a large number of molecules. However, mining powerful quantum chemical parameters to predict log P via this approach remains a challenging task. QikProp [42] is based on a study [3] which used Monte Carlo simulations to calculate 11 parameters, including solute-solvent energies, solute dipole moment, number of solute-solvent interactions at different cutoff values, number of H-bond donors and acceptors (HBDN and HBAQ and some of their variations. These parameters made it possible to estimate a number of free energies of solvation of chemicals in hexadecane, octanol, water as well as octanol-water distribution coefficients. The equation calculated for the octanol-water coefficient is ... 

All the macroscopic properties of polymers depend on a number of different factors prominent among them are the chemical structures as well as the arrangement of the macromolecules in a dense packing [1-6]. The relationships between the microscopic details and the macroscopic properties are the topics of interest here. In principle, computer simulation is a universal tool for deriving the macroscopic properties of materials from the microscopic input [7-14]. Starting from the chemical structure, quantum mechanical methods and spectroscopic information yield effective potentials that are used in Monte Carlo (MC) and molecular dynamics (MD) simulations in order to study the structure and dynamics of these materials on the relevant length scales and time scales, and to characterize the resulting thermal and mechanical proper-... 

Equation (4-5) can be directly utilized in statistical mechanical Monte Carlo and molecular dynamics simulations by choosing an appropriate QM model, balancing computational efficiency and accuracy, and MM force fields for biomacromolecules and the solvent water. Our group has extensively explored various QM/MM methods using different quantum models, ranging from semiempirical methods to ab initio molecular orbital and valence bond theories to density functional theory, applied to a wide range of applications in chemistry and biology. Some of these studies have been discussed before and they are not emphasized in this article. We focus on developments that have not been often discussed. 

---