Molecular systems mechanics

Like the geometry of Euclid and the mechanics of Newton, quantum mechanics is an axiomatic subject. By making several assertions, or postulates, about the mathematical properties of and physical interpretation associated with solutions to the Scluodinger equation, the subject of quantum mechanics can be applied to understand behaviour in atomic and molecular systems. The fust of these postulates is ... 

A quantum mechanical treatment of molecular systems usually starts with the Bom-Oppenlieimer approximation, i.e., the separation of the electronic and nuclear degrees of freedom. This is a very good approximation for well separated electronic states. The expectation value of the total energy in this case is a fiinction of the nuclear coordinates and the parameters in the electronic wavefunction, e.g., orbital coefficients. The wavefiinction parameters are most often detennined by tire variation theorem the electronic energy is made stationary (in the most important ground-state case it is minimized) with respect to them. The... 

Clusters are intennediates bridging the properties of the atoms and the bulk. They can be viewed as novel molecules, but different from ordinary molecules, in that they can have various compositions and multiple shapes. Bare clusters are usually quite reactive and unstable against aggregation and have to be studied in vacuum or inert matrices. Interest in clusters comes from a wide range of fields. Clusters are used as models to investigate surface and bulk properties [2]. Since most catalysts are dispersed metal particles [3], isolated clusters provide ideal systems to understand catalytic mechanisms. The versatility of their shapes and compositions make clusters novel molecular systems to extend our concept of chemical bonding, stmcture and dynamics. Stable clusters or passivated clusters can be used as building blocks for new materials or new electronic devices [4] and this aspect has now led to a whole new direction of research into nanoparticles and quantum dots (see chapter C2.17). As the size of electronic devices approaches ever smaller dimensions [5], the new chemical and physical properties of clusters will be relevant to the future of the electronics industry. 

Traditionally, for molecular systems, one proceeds by considering the electronic Hamiltonian which consists of the quantum mechanical operators for the kinetic energy of the electrons, their mutual Coulombic repulsions, and... 

Singh, U.C., Kollman, P.A. A combined ab initio quantum mechanical and molecular mechanical method for carrying out simulations on complex molecular systems Applications to the CH3CI 4- Cl exchange reaction and gas phase protonation of polyethers. J. Comput. Chem. 7 (1986) 718-730. 

The principal idea behind the CSP approach is to use input from Classical Molecular Dynamics simulations, carried out for the process of interest as a first preliminary step, in order to simplify a quantum mechanical calculation, implemented in a subsequent, second step. This takes advantage of the fact that classical dynamics offers a reasonable description of many properties of molecular systems, in particular of average quantities. More specifically, the method uses classical MD simulations in order to determine effective... 

In molecular dynamics applications there is a growing interest in mixed quantum-classical models various kinds of which have been proposed in the current literature. We will concentrate on two of these models the adiabatic or time-dependent Born-Oppenheimer (BO) model, [8, 13], and the so-called QCMD model. Both models describe most atoms of the molecular system by the means of classical mechanics but an important, small portion of the system by the means of a wavefunction. In the BO model this wavefunction is adiabatically coupled to the classical motion while the QCMD model consists of a singularly perturbed Schrddinger equation nonlinearly coupled to classical Newtonian equations, 2.2. 

A Fortran90 library for the simulation of molecular systems using molecular mechanics (MM) and hybrid quantum mechanics/molecular mechanics (QM)/ MM) potential energy functions. http //www.ibs.fr/ext/labos/LDM/projet6/... 

IlyperChem assumes that it is easiest for yon to just use subset selection to select that portion oT the molecular system that is to be treated quantum mechanically. You can then extend the initial selection to form a convenient and universally acceptable boundary. Thus, you make a simple selection of atoms for the first pass at selecting the quantum mechanical portion. The selected atoms are quantum atomsand the iinselecied atoms are classical atoms. 

An N-atom molecular system may he described by dX Cartesian coordinates. Six independent coordinates (five for linear molecules, three fora single atom) describe translation and rotation of the system as a whole. The remaining coordinates describe the nioleciiUir configuration and the internal structure. Whether you use molecular mechanics, quantum mechanics, or a specific computational method (AMBER, CXDO. etc.), yon can ask for the energy of the system at a specified configuration. This is called a single poin t calculation. ... 

There is a lot of confusion over the meaning of the terms theoretical chemistry, computational chemistry and molecular modelling. Indeed, many practitioners use all three labels to describe aspects of their research, as the occasion demands "Theoretical chemistry is often considered synonymous with quantum mechanics, whereas computational chemistry encompasses not only quantum mechanics but also molecular mechaiucs, minimisation, simulations, conformational analysis and other computer-based methods for understanding and predicting the behaviour of molecular systems. Molecular modellers use all of these methods and so we shall not concern ourselves with semantics but rather shall consider any theoretical or computational tecluiique that provides insight into the behaviour of molecular systems to be an example of molecular modelling. If a distinction has to be... 

In our treatment of molecular systems we first show how to determine the energy for a given iva efunction, and then demonstrate how to calculate the wavefunction for a specific nuclear geometry. In the most popular kind of quantum mechanical calculations performed on molecules each molecular spin orbital is expressed as a linear combination of atomic orhilals (the LCAO approach ). Thus each molecular orbital can be written as a summation of the following form ... 

Singh U C and P A Kollman 1986. A Combined Ab Initio Quantum Mechanical and Molecule Mechanical Method for Carrying out Simulations on Complex Molecular Systems Applicatior to the CHsQ + Cr Exchange Reaction and Gas Phase Protonation of Polyethers. Journal Computational Chemistry 7 718-730. 

We shall examine the simplest possible molecular orbital problem, calculation of the bond energy and bond length of the hydrogen molecule ion Hj. Although of no practical significance, is of theoretical importance because the complete quantum mechanical calculation of its bond energy can be canied out by both exact and approximate methods. This pemiits comparison of the exact quantum mechanical solution with the solution obtained by various approximate techniques so that a judgment can be made as to the efficacy of the approximate methods. Exact quantum mechanical calculations cannot be carried out on more complicated molecular systems, hence the importance of the one exact molecular solution we do have. We wish to have a three-way comparison i) exact theoretical, ii) experimental, and iii) approximate theoretical. 

Molecular dynamics is a simulation of the time-dependent behavior of a molecular system, such as vibrational motion or Brownian motion. It requires a way to compute the energy of the system, most often using a molecular mechanics calculation. This energy expression is used to compute the forces on the atoms for any given geometry. The steps in a molecular dynamics simulation of an equilibrium system are as follows ... 

It is possible to use computational techniques to gain insight into the vibrational motion of molecules. There are a number of computational methods available that have varying degrees of accuracy. These methods can be powerful tools if the user is aware of their strengths and weaknesses. The user is advised to use ah initio or DFT calculations with an appropriate scale factor if at all possible. Anharmonic corrections should be considered only if very-high-accuracy results are necessary. Semiempirical and molecular mechanics methods should be tried cautiously when the molecular system prevents using the other methods mentioned. 

For reactions between atoms, the computation needs to model only the translational energy of impact. For molecular reactions, there are internal energies to be included in the calculation. These internal energies are vibrational and rotational motions, which have quantized energy levels. Even with these corrections included, rate constant calculations tend to lose accuracy as the complexity of the molecular system and reaction mechanism increases. 

HyperChem uses two types of methods in calculations molecular mechanics and quantum mechanics. The quantum mechanics methods implemented in HyperChem include semi-empirical quantum mechanics method and ab initio quantum mechanics method. The molecular mechanics and semi-empirical quantum mechanics methods have several advantages over ab initio methods. Most importantly, these methods are fast. While this may not be important for small molecules, it is certainly important for biomolecules. Another advantage is that for specific and well-parameterized molecular systems, these methods can calculate values that are closer to experiment than lower level ab initio techniques. 

Semi-empirical quantum mechanics methods have evolved over the last three decades. Using today s microcomputers, they can produce meaningful, often quantitative, results for large molecular systems. The roots of the methods lie in the theory of % electrons, now largely superseded by all-valence electron theories. 

You can perform quantum mechanical calculations on a part of a molecular system, such as a solute, while using molecular mechanics for the rest of the system, such as the solvent surrounding the solute. This boundary technique is available in HyperChem for all quantum mechanical methods. It is somewhat less complete with ab initio calculations than with semi-empirical calculations, however. With ab initio calculations the boundary must occur between molecules rather than inside a molecule. 

The algorithms of the mixed classical-quantum model used in HyperChem are different for semi-empirical and ab mi/io methods. The semi-empirical methods in HyperChem treat boundary atoms (atoms that are used to terminate a subset quantum mechanical region inside a single molecule) as specially parameterized pseudofluorine atoms. However, HyperChem will not carry on mixed model calculations, using ab initio quantum mechanical methods, if there are any boundary atoms in the molecular system. Thus, if you would like to compute a wavefunction for only a portion of a molecular system using ab initio methods, you must select single or multiple isolated molecules as your selected quantum mechanical region, without any boundary atoms. 

A Hamiltonian is the quantum mechanical description of an energy contribution. The exact Hamiltonian for a molecular system is ... 

---