Quantum Mechanics-Based Computational Methods

Many other approaches for finding a correct structural model are possible. A short description of ab initio, density functional and semiempirical methods are included here. This information has been summarized from the paperback Chemistry with Computation An Introduction to Spartan The Spartan program is described in Section 4.5. Another description of computational chemistry including more mathematical treatments of quantum mechanical, molecular mechanical, and statistical mechanical methods is found in the Oxford Chemistry Primers volume Computational Chemistry  

---

Mainframes are large computers comprised of a cluster of tightly coupled machines or having multiple processors. These units will often be set up for specific applications as database servers, or for handling calculations such as those generated by quantum mechanics-based computational chemistry methods. 

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

Density functional theory (DFT),32 also a semi-empirical method, is capable of handling medium-sized systems of biological interest, and it is not limited to the second row of the periodic table. DFT has been used in the study of some small protein and peptide surfaces. Nevertheless, it is still limited by computer speed and memory. DFT offers a quantum mechanically based approach from a fundamentally different perspective, using electron density with an accuracy equivalent to post Hartree-Fock theory. The ideas have been around for many years,33 34 but only in the last ten years have numerous studies been published. DFT, compared to ab initio... 

Quantum mechanical methods can now be applied to systems with up to 1000 atoms 87 this capacity is not only from advances in computer technology but also from improvements in algorithms. Recent developments in reactive classical force fields promise to allow the study of significantly larger systems.88 Many approximations can also be made to yield a variety of methods, each of which can address a range of questions based on the inherent accuracy of the method chosen. We now discuss a range of quantum mechanical-based methods that one can use to answer specific questions regarding shock-induced detonation conditions. 

It is apparent that progress in our understanding of the properties of neutral heavy elements and their ions, including very highly ionized ones, as well as their role as constituents of molecules and solids, will depend on the development of theoretical methods and computational techniques, which are based on relativistic quantum mechanics. Fairly efficient methods of this kind have already been elaborated and many versions of relativistic codes for work with isolated atoms and ions are already available and in daily use by internationally known theoretical and experimental physicists and chemists [18, 54-57],... 

Semiempirical quantum mechanics. The computational effort in ab initio calculations increases as the fonrth power of the size of the basis set, and, therefore, its appfication to large molecnles is expensive in terms of time and computer resources. Consequently, semiempirical methods treating only the valence electrons, in which some integrals are ignored or replaced by empirically based parameters, have been developed. The various semiempirical parameterizations now in nse (MNDO, AM 1, PM3, etc.) have greatly increased the molecnlar size that is accessible to quantitative modeling methods and also the accnracy of the resnlts. 

Next the rate constant is calculated for each jump. In prior work [88] we used a harmonic approximation in the flexible degrees of freedom, with the quantum mechanics-based partition function [Eq. (20)]. Such a rate constant calculation is fast in comparison to tracking the IRC. Using a free energy method (Eq. (1.117) or sampling as in Ref. [Ill]) would yield more accurate rate constants at a cost of a larger, more time consuming calculation. Dynamical corrections [Eqs. (1.103) and (1.104)] would increase accuracy but add additional computational requirements. If the software and run time are available, these more advanced methods are recommended. 

The FDET formalism [45 8, 52, 56, 59] (cf also. Refs. [62-66]), provides basic equations for the variational treatment of a quantum-mechanical subsystem embedded in a given electronic density. In order to introduce the FDET-based computational methods for describing the molecular system A embedded in the environment B created by some other molecule(s), one introduces two types of electronic densities to represent the total system AB. The first one is the density of the subsystem A defined by embedded molecule(s), PA(r), which is typically represented using one the following auxiliary quantities (1) the occupied orbitals... 

While simulations reach into larger time spans, the inaccuracies of force fields become more apparent on the one hand properties based on free energies, which were never used for parametrization, are computed more accurately and discrepancies show up on the other hand longer simulations, particularly of proteins, show more subtle discrepancies that only appear after nanoseconds. Thus force fields are under constant revision as far as their parameters are concerned, and this process will continue. Unfortunately the form of the potentials is hardly considered and the refinement leads to an increasing number of distinct atom types with a proliferating number of parameters and a severe detoriation of transferability. The increased use of quantum mechanics to derive potentials will not really improve this situation ab initio quantum mechanics is not reliable enough on the level of kT, and on-the-fly use of quantum methods to derive forces, as in the Car-Parrinello method, is not likely to be applicable to very large systems in the foreseeable future. 

Ab initio methods, unlike either molecular mechanics or semi-empirical methods, use no experimental parameters in their computations. Instead, their computations are based solely on the laws of quantum mechanics—the first principles referred to in the name ah initio—and on the values of a small number of physical constants ... 

Atoms defined in this way can be treated as quantum-mechanically distinct systems, and their properties may be computed by integrating over these atomic basins. The resulting properties are well-defined and are based on physical observables. This approach also contrasts with traditional methods for population analysis in that it is independent of calculation method and basis set. 

The definition of "concepts" must be accompanied by explicit recipes for computing them is actual cases. There is no more space in theoretical chemistty for "driving forces", "effects, etc. not accompanied by specific rules for their quantification. The impact of a new "concept will be greater if the rules of quantifications are not restricted to ad hoc methods, but related to methods of general use in molecular quantum mechanics. A concept based exclusively on some specific features of a given method, e g. the extended Hiickel method, is less robust than a concurrent concept which may be quantified also using other levels of the theory. 

---