With recent advances in electronic structure methods, first-principles calculations of electronic response properties, such as linear and nonlinear polarizabilities, have become possible for molecules with more than 100 atoms. Basis set incompleteness is typically the main source of error in such calculations since traditional diffuse augmented basis sets are too costly to use or suffer from near linear dependence. To address this problem, we construct the first comprehensive set of property-optimized augmented basis sets for elements H-Rn except lanthanides. The new basis sets build on the Karlsruhe segmented contracted basis sets of split-valence to quadruple-zeta valence quality and add a small number of moderately diffuse basis functions. The exponents are determined variationally by maximization of atomic Hartree-Fock polarizabilities using analytical derivative methods. The performance of the resulting basis sets is assessed using a set of 313 molecular static Hartree-Fock polarizabilities. The mean absolute basis set errors are 3.6%, 1.1%, and 0.3% for property-optimized basis sets of split-valence, triple-zeta, and quadruple-zeta valence quality, respectively. Density functional and second-order Møller-Plesset polarizabilities show similar basis set convergence. We demonstrate the efficiency of our basis sets by computing static polarizabilities of icosahedral fullerenes up to C(720) using hybrid density functional theory.

Turbomole is a highly optimized software package for large-scale quantum chemical simulations of molecules, clusters, and periodic solids. It uses Gaussian basis sets and specializes on predictive electronic structure methods with excellent cost to performance characteristics, such as (time-dependent) density functional theory (TDDFT), second-order Møller-Plesset theory, and explicitly correlated coupled cluster (CC) methods. These methods are combined with ultra-efficient and numerically stable algorithms such as integral-direct and Laplace transform methods, resolution-of-the-identity, pair natural orbitals, and fast multipole and low-order scaling techniques. Apart from energies and structures, a variety of optical, electric, and magnetic properties are accessible from analytical energy derivatives for electronic ground and excited states. Recent additions include post-Kohn-Sham calculations within the random-phase approximation, periodic calculations, spin-orbit couplings, explicitly correlated CC singles doubles and perturbative triples methods, CC singles doubles excitation energies, and nonadiabatic molecular dynamics simulations using TDDFT. A dedicated graphical user interface and a user support network are also available.

This work presents theory, implementation, and validation of excited state properties obtained from time-dependent density functional theory (TDDFT). Based on a fully variational expression for the excited state energy, a compact derivation of first order properties is given. We report an implementation of analytic excited state gradients and charge moments for local, gradient corrected, and hybrid functionals, as well as for the configuration interaction singles (CIS) and time-dependent Hartree–Fock (TDHF) methods. By exploiting analogies to ground state energy and gradient calculations, efficient techniques can be transferred to excited state methods. Benchmark results demonstrate that, for low-lying excited states, geometry optimizations are not substantially more expensive than for the ground state, independent of the molecular size. We assess the quality of calculated adiabatic excitation energies, structures, dipole moments, and vibrational frequencies by comparison with accurate experimental data for a variety of excited states and molecules. Similar trends are observed for adiabatic excitation energies as for vertical ones. TDDFT is more robust than CIS and TDHF, in particular, for geometries differing significantly from the ground state minimum. The TDDFT excited state structures, dipole moments, and vibrational frequencies are of a remarkably high quality, which is comparable to that obtained in ground state density functional calculations. Thus, yielding considerably more accurate results at similar computational cost, TDDFT rivals CIS as a standard method for calculating excited state properties in larger molecules.

No abstract

We present Gaussian basis sets of quadruple zeta valence quality with a segmented contraction scheme for atoms H to Kr. This extends earlier work on segmented contracted split valence (SV) and triple zeta valence (TZV) basis sets. Contraction coefficients and orbital exponents are fully optimized in atomic Hartree–Fock (HF) calculations. As opposed to other quadruple zeta basis sets, the basis set errors in atomic ground-state HF energies are less than 1 mEh and increase smoothly across the Periodic Table, while the number of primitives is comparably small. Polarization functions are taken partly from previous work, partly optimized in atomic MP2 calculations, and for a few cases determined at the HF level for excited atomic states nearly degenerate with the ground state. This leads to basis sets denoted QZVP for HF and density functional theory (DFT) calculations, and for some atoms to a larger basis recommended for correlated treatments, QZVPP. We assess the performance of the basis sets in molecular HF, DFT, and MP2 calculations for a sample of diatomic and small polyatomic molecules by a comparison of energies, bond lengths, and dipole moments with results obtained numerically or using very large basis sets. It is shown that basis sets of quadruple zeta quality are necessary to achieve an accuracy of 1 kcal/mol per bond in HF and DFT atomization energies. For compounds containing third row as well as alkaline and earth alkaline metals it is demonstrated that the inclusion of high-lying core orbitals in the active space can be necessary for accurate correlated treatments. The QZVPP basis sets provide sufficient flexibility to polarize the core in those cases. All test calculations indicate that the new basis sets lead to consistent accuracies in HF, DFT, or correlated treatments even in critical cases where other basis sets may show deficiencies.

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

A combined experimental and theoretical study of small gold cluster anions is performed. The experimental effort consists of ion mobility measurements that lead to the assignment of the collision cross sections for the different cluster sizes at room temperature. The theoretical study is based on ab initio molecular dynamics calculations with the goal to find energetically favorable candidate structures. By comparison of the theoretical results with the measured collision cross sections as well as vertical detachment energies ͑VDEs͒ from the literature, we assign structures for the small Au n Ϫ ions (nϽ13) and locate the transition from planar to three-dimensional structures. While a unique assignment based on the observed VDEs alone is generally not possible, the collision cross sections provide a direct and rather sensitive measure of the cluster structure. In contrast to what was expected from other metal clusters and previous theoretical studies, the structural transition occurs at an unusually large cluster size of twelve atoms.

We investigate the performance of contemporary semilocal and hybrid density functionals for bond energetics, structures, dipole moments, and harmonic frequencies of 3d transition-metal (TM) compounds by comparison with gas-phase experiments. Special attention is given to the nonempirical metageneralized gradient approximation (meta-GGA) of Tao, Perdew, Staroverov, and Scuseria (TPSS) [Phys. Rev. Lett. 91, 146401 (2003)], which has been implemented in TURBOMOLE for the present work. Trends and error patterns for classes of homologous compounds are analyzed, including dimers, monohydrides, mononitrides, monoxides, monofluorides, polyatomic oxides and halogenides, carbonyls, and complexes with organic pi ligands such as benzene and cyclopentadienyl. Weakly bound systems such as Ca(2), Mn(2), and Zn(2) are discussed. We propose a reference set of reaction energies for benchmark purposes. Our all-electron results with quadruple zeta valence basis sets validate semilocal density-functional theory as the workhorse of computational TM chemistry. Typical errors in bond energies are substantially larger than in (organic) main group chemistry, however. The Becke-Perdew'86 [Phys. Rev. A 38, 3098 (1988); Phys. Rev. B 33, 8822 (1986)] GGA and the TPSS meta-GGA have the best price/performance ratio, while the TPSS hybrid functional achieves a slightly lower mean absolute error in bond energies. The popular Becke three-parameter hybrid B3LYP underbinds significantly and tends to overestimate bond distances; we give a possible explanation for this. We further show that hybrid mixing does not reduce the width of the error distribution on our reference set. The error of a functional for the s-d transfer energy of a TM atom does not predict its error for TM bond energies and bond lengths. For semilocal functionals, self-interaction error in one- and three-electron bonds appears to be a major source of error in TM reaction energies. Nevertheless, TPSS predicts the correct ground-state symmetry in the vast majority of cases and rarely fails qualitatively. This further confirms TPSS as a general purpose functional that works throughout the periodic table. We also give workstation timing comparisons for the 645-atom protein crambin.