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.

Adiabatic time-dependent density functional methods for excited state properties

Efficient, approximate and parallel Hartree–Fock and hybrid DFT calculations. A ‘chain-of-spheres’ algorithm for the Hartree–Fock exchange

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...

/pdf/property-optimized-gaussian-basis-sets-for-molecular-80kgu7ex4o.pdf

Property-optimized gaussian basis sets for molecular response calculations.

Prediction of molecular properties and molecular spectroscopy with density functional theory: From fundamental theory to exchange-coupling

We propose and discuss an efficient scheme for the in silico sampling for parts of the molecular chemical space by semiempirical tight-binding methods combined with a meta-dynamics driven search algorithm. The focus of this work is set on the generation of proper thermodynamic ensembles at a quantum chemical level for conformers, but similar procedures for protonation states, tautomerism and non-covalent complex geometries are also discussed. The conformational ensembles consisting of all significantly populated minimum energy structures normally form the basis of further, mostly DFT computational work, such as the calculation of spectra or macroscopic properties. By using basic quantum chemical methods, electronic effects or possible bond breaking/formation are accounted for and a very reasonable initial energetic ranking of the candidate structures is obtained. Due to the huge computational speedup gained by the fast low-cost quantum chemical methods, overall short computation times even for systems with hundreds of atoms (typically drug-sized molecules) are achieved. Furthermore, specialized applications, such as sampling with implicit solvation models or constrained conformational sampling for transition-states, metal-, surface-, or noncovalently bound complexes are discussed, opening many possible applications in modern computational chemistry and drug discovery. The procedures have been implemented in a freely available computer code called CREST, that makes use of the fast and reliable GFNn-xTB methods.

H Öhm

Papers

Auxiliary basis sets to approximate Coulomb potentials (Chem. Phys. Letters 240 (1995) 283-290)