A description of the ab initio quantum chemistry package GAMESS is presented. Chemical systems containing atoms through radon can be treated with wave functions ranging from the simplest closed-shell case up to a general MCSCF case, permitting calculations at the necessary level of sophistication. Emphasis is given to novel features of the program. The parallelization strategy used in the RHF, ROHF, UHF, and GVB sections of the program is described, and detailed speecup results are given. Parallel calculations can be run on ordinary workstations as well as dedicated parallel machines. © John Wiley & Sons, Inc.

General atomic and molecular electronic structure system

"Chemists familiar with conventional quantum mechanics will applaud and benefit greatly from this particularly instructive, thorough and clearly written exposition of density functional theory: its basis, concepts, terms, implementation, and performance in diverse applications. Users of DFT for structure, energy, and molecular property computations, as well as reaction mechanism studies, are guided to the optimum choices of the most effective methods. Well done!" Paul von RaguE Schleyer "A conspicuous hole in the computational chemist's library is nicely filled by this book, which provides a wide-ranging and pragmatic view of the subject.[...It] should justifiably become the favorite text on the subject for practioneers who aim to use DFT to solve chemical problems." J. F. Stanton, J. Am. Chem. Soc. "The authors' aim is to guide the chemist through basic theoretical and related technical aspects of DFT at an easy-to-understand theoretical level. They succeed admirably." P. C. H. Mitchell, Appl. Organomet. Chem. "The authors have done an excellent service to the chemical community. [...] A Chemist's Guide to Density Functional Theory is exactly what the title suggests. It should be an invaluable source of insight and knowledge for many chemists using DFT approaches to solve chemical problems." M. Kaupp, Angew. Chem.

A Chemist's Guide to Density Functional Theory

We present a new hybrid meta exchange-correlation functional, called M05-2X, for thermochemistry, thermochemical kinetics, and noncovalent interactions. We also provide a full discussion of the new M05 functional, previously presented in a short communication. The M05 functional was parametrized including both metals and nonmetals, whereas M05-2X is a high-nonlocality functional with double the amount of nonlocal exchange (2X) that is parametrized only for nonmetals. In particular, M05 was parametrized against 35 data values, and M05-2X is parametrized against 34 data values. Both functionals, along with 28 other functionals, have been comparatively assessed against 234 data values:  the MGAE109/3 main-group atomization energy database, the IP13/3 ionization potential database, the EA13/3 electron affinity database, the HTBH38/4 database of barrier height for hydrogen-transfer reactions, five noncovalent databases, two databases involving metal−metal and metal−ligand bond energies, a dipole moment databas...

/pdf/design-of-density-functionals-by-combining-the-method-of-3dsa6vb61i.pdf

Design of Density Functionals by Combining the Method of Constraint Satisfaction with Parametrization for Thermochemistry, Thermochemical Kinetics, and Noncovalent Interactions.

A redundant internal coordinate system for optimizing molecular geometries is constructed from all bonds, all valence angles between bonded atoms, and all dihedral angles between bonded atoms. Redundancies are removed by using the generalized inverse of the G matrix; constraints can be added by using an appropriate projector. For minimizations, redundant internal coordinates provide substantial improvements in optimization efficiency over Cartesian and nonredundant internal coordinates, especially for flexible and polycyclic systems. Transition structure searches are also improved when redundant coordinates are used and when the initial steps are guided by the quadratic synchronous transit approach. © 1996 by John Wiley & Sons, Inc.

Using redundant internal coordinates to optimize equilibrium geometries and transition states

Advances in theory and algorithms for electronic structure calculations must be incorporated into program packages to enable them to become routinely used by the broader chemical community. This work reviews advances made over the past five years or so that constitute the major improvements contained in a new release of the Q-Chem quantum chemistry package, together with illustrative timings and applications. Specific developments discussed include fast methods for density functional theory calculations, linear scaling evaluation of energies, NMR chemical shifts and electric properties, fast auxiliary basis function methods for correlated energies and gradients, equation-of-motion coupled cluster methods for ground and excited states, geminal wavefunctions, embedding methods and techniques for exploring potential energy surfaces.

/pdf/advances-in-methods-and-algorithms-in-a-modern-quantum-29c5vso11m.pdf

Advances in methods and algorithms in a modern quantum chemistry program package

An algorithm for locating transition states designed for use in the ab initio program package GAUSSIAN 82 is presented. It is capable of locating transition states even if started in the wrong region of the energy surface, and, by incorporating the ideas on hessian mode following due to Cerjan and Miller, can locate transition states for alternative rearrangement/dissociation reactions from the same initial starting point. It can also be used to locate minima.

An algorithm for the location of transition states

The performance of four commonly used density functionals (VWN, BLYP, BP91, and Becke's original three-parameter approximation to the adiabatic connection formula, referred to herein as the adiabatic connection method or ACM) was studied with a series of six Gaussian-type atomic basis sets [DZP, 6–31G**, DZVP, TZVP, TZ2P, and uncontracted aug-cc-pVTZ (UCC)]. The geometries and dipole moments of over 100 first-row and second-row molecules and reaction energies of over 300 chemical reactions involving such molecules were computed using each of the four density functionals in combination with each of the six basis sets. The results were compared to experimentally determined values. Based on overall mean absolute theory versus experiment errors, it was found that ACM is the best choice for predictions of both energies of reaction [overall mean absolute theory versus experiment error (MATvEE) of 4.7 kcal/mol with our most complete (UCC) basis set] and molecular geometries (overall MATvEE of 0.92 pm for bond distances and 0.88° for bond angles with the UCC basis set). For routine calculations with moderate basis sets (those of double-ζ type: DZP, 6–31G**, and DZVP) the DZVP basis set was, on average, the best choice. There were, however, examples of reactions where significantly larger basis sets were required to achieve reasonable accuracy (errors ≤ 5 kcal/mol). For dipole moments, ACM, BP91, and BLYP performed comparably (overall MATvEE of 0.071, 0.067, and 0.059 debye, respectively, with the UCC basis set). Basis sets that include additional polarization functions and diffuse functions were found to be important for accurate density functional theory predictions of dipole moments. © 1997 by John Wiley & Sons, Inc.

Molecular energies and properties from density functional theory: Exploring basis set dependence of Kohn—Sham equation using several density functionals

We have investigated the performance of the new optimized exchange functional (OPTX) developed by Handy and Cohen [Mol. Phys. 99, 403 (2001)] for predicting geometries, heats of reaction, and barrier heights for twelve organic reactions (six closed-shell and six radical). OPTX has been used in conjunction with, among others, the well-known Lee–Yang–Parr (LYP) correlational functional to form two new functionals, OLYP and O3LYP. These are similar to the well-established BLYP and B3LYP functionals, respectively, with OPTX replacing the standard Becke exchange functional, B88. Our results strongly support claims made by their developers that OLYP is superior to BLYP, and essentially renders it obsolete. The computed OLYP heats of reaction, barrier heights, and even molecular geometries (with larger basis sets), are comparable with, if not better than, the corresponding B3LYP values. The O3LYP functional is overall better than B3LYP, albeit not by much. Both OLYP and O3LYP are among the best functionals currently available; the performance of OLYP in particular is noteworthy given that this functional includes no exact exchange.

Assessment of the Handy–Cohen optimized exchange density functional for organic reactions

A numerical algorithm for locating both minima and transition states designed for use in the ab initio program package GAUSSIAN 82 is presented It is based on the RFO method of Simons and coworkers and is effectively the numerical version of an analytical algorithm (OPT = EF) previously published in this journal The algorithm is designed to make maximum use of external second derivative information obtained from prior optimizations at lower levels of theory It can be used with any wave function for which an energy can be calculated and is about two to three times faster than the default DFP algorithm (OPT = FP) supplied with GAUSSIAN 82

An algorithm for geometry optimization without analytical gradients

A comparison is made between geometry optimization in Cartesian coordinates, using an appropriate initial Hessian, and natural internal coordinates. Results on 33 different molecules covering a wide range of symmetries and structural types demonstrate that both coordinate systems are of comparable efficiency. There is a marked tendency for natural internals to converge to global minima whereas Cartesian optimizations converge to the local minimum closest to the starting geometry. Because they can now be generated automatically from input Cartesians, natural internals are to be preferred over Z-matrix coordinates. General optimization strategies using internal coordinates and/or Cartesians are discussed for both unconstrained and constrained optimization. © John Wiley & Sons, Inc.

Jon Baker

Papers

An algorithm for the location of transition states

Molecular energies and properties from density functional theory: Exploring basis set dependence of Kohn—Sham equation using several density functionals

Assessment of the Handy–Cohen optimized exchange density functional for organic reactions

An algorithm for geometry optimization without analytical gradients

Techniques for geometry optimization: a comparison of Cartesian and natural internal coordinates