Quantum-mechanical computer methods

Quantum-mechanical studies on the tautomerism of heterocyclic compounds involve, in general, two aspects. The first deals with the prediction of physicochemical properties of defined tautomeric forms (e.g., ultraviolet spectra, dipole moments, ionization potentials, etc.). This seems to be easy to handle. Using any semiempirical or nonempirical quantum-mechanical computational method, depending on approximations involved in the method, we are able to calculate properties that, more or less, agree with experimental values. Calculations of this type do not contribute to a direct estimation of the relative stability of the tautomers, however they are particularly important for cases in which a tautomeric form of a compound is so rare that it is not possible to measure it directly. 

High level quantum mechanical computational methods (loosely called ab initio methods) can be used to predict molecular properties fairly accurately. These calculations typically produce much more information than can be incorporated into a classical, effective potential model. 

With the powerful quantum mechanical computational methods that are presently available, it is possible to obtain accurate at initio electronic wavefunctions for modest to large size systems. The availability of powerful computers, especially workstations, and vector and massively parallel supercomputers, makes it possible to take extensive advantage of these computational methods. The recent advances in both computer architecture and numerical methods have been truly breathtaking and we are now able to obtain reasonably accurate solutions for isolated molecules as well as for cluster models of condensed matter. In this review, our concern is for the use of cluster models to describe chemistry at solid surfaces in particular, the chemical bond formed between adsorbates and solid substrates. 

The preferable theoretical tools for the description of dynamical processes in systems of a few atoms are certainly quantum mechanical calculations. There is a large arsenal of powerful, well established methods for quantum mechanical computations of processes such as photoexcitation, photodissociation, inelastic scattering and reactive collisions for systems having, in the present state-of-the-art, up to three or four atoms, typically. " Both time-dependent and time-independent numerically exact algorithms are available for many of the processes, so in cases where potential surfaces of good accuracy are available, excellent quantitative agreement with experiment is generally obtained. In addition to the full quantum-mechanical methods, sophisticated semiclassical approximations have been developed that for many cases are essentially of near-quantitative accuracy and certainly at a level sufficient for the interpretation of most experiments.These methods also are com-... 

Despite advent of theoretical methods and techniques and faster computers, no single theoretical method seems to be capable of reliable computational studies of reactivities of biocatalysts. Ab initio quantum mechanical (QM) methods may be accurate but are still too expensive to apply to large systems like biocatalysts. Semi-empirical quantum methods are not as accurate but are faster, but may not be fast enough for long time simulation of large molecular systems. Molecular mechanics (MM) force field methods are not usually capable of dealing with bond-breaking and formation... 

The currently available quantum chemical computational methods and computer programs have not been utilized to their potential in elucidating the electronic origin of zeolite properties. As more and more physico-chemical methods are used successfully for the description and characterization of zeolites, (e.g. (42-45)), more questions will also arise where computational quantum chemistry may have a useful contribution towards the answer, e.g. in connection with combined approaches where zeolites and metal-metal bonded systems (e.g. (46,47)) are used in combination. The spectacular recent and projected future improvements in computer technology are bound to enlarge the scope of quantum chemical studies on zeolites. Detailed studies on optimum intercavity locations for a variety of molecules, and calculations on conformation analysis and reaction mechanism in zeolite cavities are among the promises what an extrapolation of current developments in computational quantum chemistry and computer technology holds out for zeolite chemistry. 

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. 

The size of the atoms and the rigidity of the bonds, bond angles, torsions, etc. are determined empirically, that is, they are chosen to reproduce experimental data. Electrons are not part of the MM description, and as a result, several key chemical phenomena cannot be reproduced by this method. Nevertheless, MM methods are orders of magnitude cheaper from a computational point of view than quantum mechanical (QM) methods, and because of this, they have found a preferential position in a number of areas of computational chemistry, like conformational analysis of organic compounds or molecular dynamics. 

A2 is precisely the 2-RDMC, and from Eq. (15) we note that expectation values for the composite A + B system can be computed using either D2 alone, or Di = Ai together with A2. Erom the standpoint of exact quantum mechanics, either method yields exactly the same expectation value and, in particular, both methods respect the extensivity of the electronic energy. If D2 is calculated by means of approximate quantum mechanics, however, one cannot generally expect that extensivity will be preserved, since exchange terms mingle the coordinates on different subsystems, and exact cancellation cannot be anticipated unless built in from the start. Methods that respect this separability by construction are said to be size-consistent [40-42]. 

Thdry, V., Rinaldi, D., Rivail, J.-L., Maigret, B., and Ferenczy, G. G. 1994. Quantum Mechanical Computations on Very Large Molecular Systems The Local Self-consistent Field Method , J. Comput. Chem., 15, 269. 

The numerical evaluation of the energies of orbitals and states is fundamentally a matter of making quantum mechanical computations. As indicated in Chapter 1, quantum mechanics per se is not the subject of this book, and indeed we have tried in general to avoid any detailed treatment of methods for solving the wave equation, emphasis being placed on the properties that the wave functions must have purely for reasons of symmetry and irrespective of their explicit analytical form. However, this discussion of the symmetry aspects of ligand field theory would be artificial and unsatisfying without some... 

V. Thery, D. Rinaldi, J. L. Rivail, B. Maigret and G. G. Ferenczy, Quantum-mechanical computations on very large molecular-systems - the local self-consistent-field method, J. Comput. Chem., 15 (1994) 269-282. 

Methods relying on parametrised potential functions for the description of energy hypersurfaces are commonly referred to as molecular mechanics (MM) or classical mechanics and these methods have a long tradition in computational chemistry. Entire data sets of balanced potential functions and their respective parameters are referred to as force fields [23,24,25], The key advantage of MM methods is the low computational demand compared to quantum mechanical computations. 

The analysis of potential energy surfaces may be of importance for both molecular-mechanical and quantum-mechanical computations. However, due to the fact that thousands of structures instead of only one need to be optimized, the methods briefly described here are only routinely used with force field calculations 35-1. 

The thermodynamic data as well as the detonation parameters can nowadays be very reliably obtained by using quantum-mechanical computer calculations. On the one hand it is important to check experimental results, and on the other hand and even more importantly - it is important to predict the properties of potential new energetic materials without any prior experimental parameters, for example during the planning of synthetic work. Moreover, such computational methods are ideal for the estimation of the detonation parameters of newly synthesized compounds, which have not been obtained in the 50 100 g quantities which are necessary for the experimental determination of such detonation parameters (e.g. detonation velocity). 

An obvious first concern is what would be a reasonable quantum mechanical approach (method and basis set) for computing NMR chemical shifts. We will consider this in the next section. 

Of the dynamical techniques available the most rigorous and informative are the quantum mechanical dynamics methods. These methods are, however, the most sophisticated and computationally intensive to employ. Two of the most widely used quantum dynamics techniques are quantum scattering (QS) [35] and wavepacket (WP) [125] analysis. 

---