Molecular dynamics simulation computational chemistry

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. 

The visualization of volumetric properties is more important in other scientific disciplines (e.g., computer tomography in medicine, or convection streams in geology). However, there are also some applications in chemistry (Figure 2-125d), among which only the distribution of water density in molecular dynamics simulations will be mentioned here. Computer visualization of this property is usually realized with two or three dimensional textures [203]. 

Field M J, P A Bash and M Karplus 1990. A Combined Quantum Mechanical and Molecular Mechanical Potential for Molecular Dynamics Simulations. Journal of Computational Chemistry 11 700-733. 

See van Gunsteren, W.F. Berendsen, H.J.C. Computer simulation of molecular dynamics-methodology, applications, and perspectives in chemistry Angewandre Chemie, International Edition in English, 29 992-1023, 1990, and Karplus, M. Petsko, G.A. Molecular dynamics simulations in biology Nature 347 631-639, 1990. 

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. 

Onufriev, A. Implicit solvent models in molecular dynamics simulations a brief overview. In Annual Reports in Computational Chemistry (eds R.A. Wheeler and D.C. Spellmeyer), Vol. 4, Elsevier, Amsterdam, 2008, pp. 125-37. 

Grassy, G., Yasri, A., Sans, P., Armbruster, A. M., Lahana, R. and Fauchere, J. L. (1997) Molecular dynamics simulations as a tool to investigate the three-dimensional diversity of peptide and pseudopeptide libraries. In Computer-assisted lead finding and optimization. Current tools for medicinal chemistry, Van de Waterbeemd, H., Testa, B., and Folkers, G. (eds.), Wiley-VCH, Basel, pp. 209-222. 

Hofer, T. S. Randolf, B. R. Rode, . M. In Challenges and Advances in Computational Chemistry and Physicsvol. 6 Ed. Canute, S. Molecular Dynamics Simulation Methods Including Quantum Effects. Springer Heidelberg, 2008. 

The credit load for die computational chemistry laboratory course requires that the average student should be able to complete almost all of the work required for the course within die time constraint of one four-hour laboratory period per week. This constraint limits the material covered in the course. Four principal computational methods have been identified as being of primary importance in the practice of chemistry and thus in the education of chemistry students (1) Monte Carlo Methods, (2) Molecular Mechanics Methods, (3) Molecular Dynamics Simulations, and (4) Quantum Chemical Calculations. Clearly, other important topics could be added when time permits. These four methods are developed as separate units, in each case beginning with die fundamental principles including simple programming and visualization, and building to the sophisticated application of the technique to a chemical problem. 

The first task of chemoinformatics is to transform chemical knowledge, such as molecular structures and chemical reactions, into computer-legible digital information. The digital representations of chemical information are the foundation for all chemoin-formatic manipulations in computer. There are many file formats for molecular information to be imported into and exported from computer. Some formats contain more information than others. Usually, intended applications will dictate which format is more suitable. For example, in a quantum chemistry calculation the molecular input file usually includes atomic symbols with three-dimensional (3D) atomic coordinates as the atomic positions, while a molecular dynamics simulation needs, in addition, atom types, bond status, and other relevant information for defining a force field. 

Crown ethers continue to be one of the most useful parts of supramolecular chemistry/91 From the beginning computations of metal ions complexes with synthetic ionophores/101 which have been aptly reviewed/111 emphasized the importance of including explicitly solvation in free energy calculations, also with ab initio calculations on calixarene complexes/121 Molecular dynamics simulations of 18-crown-6 ether complexes in aqueous solutions predict too low affinities, but at least correctly reproduce the sequence trend K+ > Rb+ > Cs+ > Na+. However, only the selection of K+ over Rb+ and Cs+ is ascribed to the cation size relative to that of the crown cavity, whereas K+ appears in these calculations to be selected over Na+ as consequence of the greater free energy penalty involved in displacing water molecules ftomNa/1131... 

The advent of high-speed computers, availability of sophisticated algorithms, and state-of-the-art computer graphics have made plausible the use of computationally intensive methods such as quantum mechanics, molecular mechanics, and molecular dynamics simulations to determine those physical and structural properties most commonly involved in molecular processes. The power of molecular modeling rests solidly on a variety of well-established scientific disciplines including computer science, theoretical chemistry, biochemistry, and biophysics. Molecular modeling has become an indispensable complementary tool for most experimental scientific research. 

Quantum chemical methods are valuable tools for studying atmospheric nucle-ation phenomena. Molecular geometries and binding energies computed using electronic structure methods can be used to determine potential parameters for classical molecular dynamic simulations, which in turn provide information on the dynamics and qualitative energetics of nucleation processes. Quantum chemistry calculations can also be used to obtain accurate and reliable information on the fundamental chemical and physical properties of molecular systems relevant to nucleation. Successful atmospheric applications include investigations on the hydration of sulfuric acid and the role of ammonia, sulfur trioxide and/or ions... 

Einally, we draw attention to several review articles in this area. In 1986 Lowdin 23 considered various aspects of the historical development of computational QC in view of the development of both conventional supercomputers and large-scale parallel computers. More recently, Weineri24 presented a discussion on the programming of parallel computers and their use in molecular dynamics simulations, free energy perturbation, and large scale ab initio calculations, as well as the use of very elaborate graphical display programs in chemistry research. We also note a review on the use of parallel processors in... 

The aim of this chapter is to provide the reader with an overview of the potential of modern computational chemistry in studying catalytic and electro-catalytic reactions. This will take us from state-of-the-art electronic structure calculations of metal-adsorbate interactions, through (ab initio) molecular dynamics simulations of solvent effects in electrode reactions, to lattice-gas-based Monte Carlo simulations of surface reactions taking place on catalyst surfaces. Rather than extensively discussing all the different types of studies that have been carried out, we focus on what we believe to be a few representative examples. We also point out the more general theory principles to be drawn from these studies, as well as refer to some of the relevant experimental literature that supports these conclusions. Examples are primarily taken from our own work other recent review papers, mainly focused on gas-phase catalysis, can be found in [1-3]. 

The next section gives a brief overview of the main computational techniques currently applied to catalytic problems. These techniques include ab initio electronic structure calculations, (ab initio) molecular dynamics, and Monte Carlo methods. The next three sections are devoted to particular applications of these techniques to catalytic and electrocatalytic issues. We focus on the interaction of CO and hydrogen with metal and alloy surfaces, both from quantum-chemical and statistical-mechanical points of view, as these processes play an important role in fuel-cell catalysis. We also demonstrate the role of the solvent in electrocatalytic bondbreaking reactions, using molecular dynamics simulations as well as extensive electronic structure and ab initio molecular dynamics calculations. Monte Carlo simulations illustrate the importance of lateral interactions, mixing, and surface diffusion in obtaining a correct kinetic description of catalytic processes. Finally, we summarize the main conclusions and give an outlook of the role of computational chemistry in catalysis and electrocatalysis. 

V. van Speybroeck and R. J. Meier, A recent development in computational chemistry Chemical reactions from first principles molecular dynamics simulations, Chem. Soc. Rev., 32 (2003) 151-157. 

So far, our discussion has focussed on stationary quantum chemical methods, which yield results for fixed atomic nuclei, i.e. for frozen molecular structures like minimum structures on the Born-Oppenheimer potential energy surface. Processes in supramolecular assemblies usually feature prominent dynamical effects, which can only be captured through explicit molecular dynamics or Monte Carlo simulations [95-98]. Molecular dynamics simulations proved to be a useful tool for studying the detailed microscopic dynamic behavior of many-particle systems as present in physics, chemistry and biology. The aim of molecular dynamics is to study a system by recreating it on the computer as close to nature as possible, i.e. by simulating the dynamics of a system in all microscopic detail over a physical length of time relevant to properties of interest. 

---