Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set

A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.

Semiempirical GGA-type density functional constructed with a long-range dispersion correction.

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

We present an analysis of the performances of a parameter free density functional model (PBE0) obtained combining the so called PBE generalized gradient functional with a predefined amount of exact exchange. The results obtained for structural, thermodynamic, kinetic and spectroscopic (magnetic, infrared and electronic) properties are satisfactory and not far from those delivered by the most reliable functionals including heavy parameterization. The way in which the functional is derived and the lack of empirical parameters fitted to specific properties make the PBE0 model a widely applicable method for both quantum chemistry and condensed matter physics.

Toward reliable density functional methods without adjustable parameters: The PBE0 model

6.2.2. Definition of Effective Properties 3064 6.3. Response Properties to Magnetic Fields 3066 6.3.1. Nuclear Shielding 3066 6.3.2. Indirect Spin−Spin Coupling 3067 6.3.3. EPR Parameters 3068 6.4. Properties of Chiral Systems 3069 6.4.1. Electronic Circular Dichroism (ECD) 3069 6.4.2. Optical Rotation (OR) 3069 6.4.3. VCD and VROA 3070 7. Continuum and Discrete Models 3071 7.1. Continuum Methods within MD and MC Simulations 3072

/pdf/quantum-mechanical-continuum-solvation-models-3edv963g14.pdf

Quantum mechanical continuum solvation models.

Convergence acceleration of iterative sequences. the case of scf iteration

Transferable Scaling Factors for Density Functional Derived Vibrational Force Fields

Fully optimized geometries, complete in- and out-of-plane force fields, and dipole moment derivatives have been
calculated for the title compounds at the ab initio HartreeFock level using the 4-21 Gaussian basis set. The theoretical information is combined with experimental data by fitting the calculated force constants through a few parameters to the observed frequencies to obtain the final, scaled quantum mechanical (SQM) force fields. Recommendations for a standard procedure of this type are given. The SQM force fields give excellent reproduction of the fundamental frequencies and are considered as approaching the best accuracy which can be achieved in a harmonic treatment. The infrared intensities obtained at this level of theory are only qualitative estimates, but they are still useful for making assignments more reliable.

Combination of theoretical ab initio and experimental information to obtain reliable harmonic force constants. Scaled quantum mechanical (QM) force fields for glyoxal, acrolein, butadiene, formaldehyde, and ethylene

Can (semi) local density functional theory account for the london dispersion forces

Two suggestions are made to increase the efficiency and accuracy of ab initio optimization of molecular geometries. To improve the convergence of the optimization, a set of internal coordinates, the natural valence coordinates, is suggested. These coordinates originate from vibrational spectroscopy and reduce both harmonic and anharmonic coupling terms in the potential function as much as possible in a purely geometrical definition. The natural valence coordinates are local, eliminate most redundancies, and conform to local pseudosymmetry. Special attention has been paid to ring systems. A computer program has been included in our program system TX90 to generate the natural internal coordinates automatically. The usefulness of these coordinates is demonstrated by numerous examples of ab initio geometry optimization. Starting with a geometry preoptimized by molecular mechanics and using a simple diagonal estimate of the Hessian in conjunction with the GDIIS optimization technique, we usually achieved convergence in 8-15 steps, even for large molecules. It is demonstrated that, due to the reduction in anharmonic couplings, natural coordinates are superior to Cartesian or other simple internal coordinates, even when an accurate initial Hessian is available. Constrained optimization and the location of transition states are also discussed. The gradient optimization method has been generalized to handle redundancies; this is necessary in some complex polycyclic molecules and is illustrated on, among others, the porphine molecule. To increase the accuracy of relatively low-level calculations, empirical corrections to ab initio SCF geometries are suggested in the form of “offset forces” acting along bonds. We recommend offset forces for the most important bonds, to be used with the 4-21G(*) and the 6-31G* basis sets. Based on 130 comparisons, the mean absolute error between theoretical and experimental bond lengths is reduced this way from 0.014 to 0.005 A.

Peter Pulay

Papers

Convergence acceleration of iterative sequences. the case of scf iteration

Transferable Scaling Factors for Density Functional Derived Vibrational Force Fields

Combination of theoretical ab initio and experimental information to obtain reliable harmonic force constants. Scaled quantum mechanical (QM) force fields for glyoxal, acrolein, butadiene, formaldehyde, and ethylene

Can (semi) local density functional theory account for the london dispersion forces

The calculation of ab initio molecular geometries: efficient optimization by natural internal coordinates and empirical correction by offset forces