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

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.

Various contracted Gaussian basis sets for atoms up to Kr are presented which have been determined by optimizing atomic self‐consistent field ground state energies with respect to all basis set parameters, i.e., orbital exponents and contraction coefficients.

/pdf/fully-optimized-contracted-gaussian-basis-sets-for-atoms-li-2ce6znei10.pdf

Fully optimized contracted Gaussian basis sets for atoms Li to Kr

Efficient implementation of the gauge-independent atomic orbital method for NMR chemical shift calculations

The general expression for the exact forces on the nuclei (negative derivatives of the total energy with respect to the nuclear coordinates) is applied for Hartree-Fock wavefunctions. It is suggested that force constants should be calculated by differentiating the forces numerically. This method, called the force method, is numerically more accurate and requires less computation than the customary one of differentiating the energy numerically twice. It permits the quick determination of the equilibrium geometry by relaxing the nuclear coordinates until the forces vanish. The unreliability of the methods using the Hellmann-Feynman forces is re-emphasized. The question of which force constants can be best calculated ab initio is discussed.

Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules

Systematic ab initio gradient calculation of molecular geometries, force constants, and dipole moment derivatives is described. A new basis set, denoted 4-21, is presented for first-row atoms. It is nearly equivalent to the 4-31G set but requires less computational effort. Completely optimized Hartree-Fock geometries of 18 molecules are compared using several basis sets, with and without polarization functions. The question of the best representation of molecular force fields is discussed, and a set of standardized internal coordinates is suggested for future work. Quadratic and the most important cubic force constants and dipole moment derivatives of first-row hydrides are calculated using the 4-21 basis set, and the results are compared with those from other basis sets, including near-Hartree-Fock ones. Force-field calculations on larger molecules with the 4-21 basis are summarized. A general formulation of the rotational correction to dipole moment derivatives is given.

Systematic AB Initio Gradient Calculation of Molecular Geometries, Force Constants, and Dipole Moment Derivatives

The gradient geometry‐optimization procedure is reformulated in terms of redundant internal coordinates. By replacing the matrix inverse with the generalized inverse, the usual Newton–Raphson–type algorithms can be formulated in exactly the same way for redundant and nonredundant coordinates. Optimization in redundant coordinates is particularly useful for bridged polycyclic compounds and cage structures where it is difficult to define physically reasonable redundancy‐free internal coordinates. This procedure, already used for the geometry optimization of porphine, C20N4H14, is illustrated here at the ab initio self‐consistent‐field level for the four‐membered ring azetidine, for bicyclo[2.2.2]octane, and for the four‐ring system C16O2H22, the skeleton of taxol.

Peter Pulay

Papers

Efficient implementation of the gauge-independent atomic orbital method for NMR chemical shift calculations

Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules

Systematic AB Initio Gradient Calculation of Molecular Geometries, Force Constants, and Dipole Moment Derivatives

Geometry optimization in redundant internal coordinates

High accuracy benchmark calculations on the benzene dimer potential energy surface