Potential-energy surfaces of local excited states from subsystem- and selective Kohn-Sham-TDDFT

Calculating excited-state potential-energy surfaces for systems with a large number of close-lying excited states requires the identification of the relevant electronic transitions for several geometric structures. Time-dependent density functional theory (TDDFT) is very efficient in such calculations, but the assignment of local excited states of the active molecule can be difficult. We compare the results of the frozen-density embedding (FDE) method with those of standard Kohn-Sham density-functional theory (KS-DFT) and simpler QM/MM-type methods. The FDE results are found to be more accurate for the geometry dependence of excitation energies than classical models. We also discuss how selective iterative diagonalization schemes can be exploited to directly target specific excitations for different structures. Problems due to strongly interacting orbital transitions and possible solutions are discussed. Finally, we apply FDE and the selective KS-TDDFT to investigate the potential energy surface of a high-lying π → π^∗ excitation in a pyridine molecule approaching a silver cluster.     

An average-of-configuration method for using Kohn-Sham density functional theory in modeling ligand-field theory

The Amsterdam Density Functional (ADF) package has been used to constrain Kohn-Sham DFT in such a fashion that a transition from KS-DFT to ligand-field theory in the form of the parametrical d(q)() model is completely well-defined. A relationship is established between the strong-field approximation of the parametrical d(2) model for the tetrahedral complexes VCl(4)(-) and VBr(4)(-) and certain fixed-orbital ADF-computed energies. In this way values for all the parameters of the d(2)() model may be computed, thus allowing the ADF results to be expressed in terms of a KS-DFT energy matrix that can be diagonalized. This means that the KS-DFT deficiency with regard to computation of nondiagonal elements has been overcome and the KS-DFT eigenenergies have become available through the KS-DFT mimicking of the ligand-field plus repulsion model. By using mutually orthogonal strong-field energy matrices, the mimicking has been further elucidated. The computed values for the empirical parameters of VCl(4)(-) and VBr(4)(-) are in good agreement with experimental data. The spectrochemical and the nephelauxetic series have been computed by including the remaining halide complexes and the quantitatively special position of F(-)() among the halides corroborated for both series.     

A detailed test study of barrier heights for the HO2 + H2O + O3 reaction with various forms of multireference perturbation theory

We report an ab initio multireference perturbation theory investigation of the HO(2) + H(2)O + O(3) reaction, with particular emphasis on the barrier heights for two possible reaction mechanisms: oxygen abstraction and hydrogen abstraction, which are identified by two distinct saddle points. These saddle points and the corresponding pre-reactive complexes were optimized at the CASSCF(11,11) level while the single point energies were calculated with three different MRPT2 theories: MRMP, CASPT2, and SC-NEVPT2. Special attention has been drawn on the "intruder state" problem and the effect of its corrections on the relative energies. The results were then compared with single reference coupled-cluster methods and also with our recently obtained Kohn-Sham density functional theory (KS-DFT) calculations [L. P. Viegas and A. J. C. Varandas, Chem. Phys., (2011)]. It is found that the relative energies of the pre-reactive complexes have a very good agreement while the MRPT2 classical barrier heights are considerably higher than the KS-DFT ones, with the SC-NEVPT2 calculations having the highest energies between the MRPT2 methods. Possible explanations have been given to account for these differences.     

Extending ligand field molecular mechanics to modelling organometallic π-bonded systems: applications to ruthenium-arenes

The ligand field molecular mechanics method has been extended to treat η(6)-arene ligands coordinated to a ruthenium(II) centre by employing a dummy atom located at the centroid of the arene ring and distributing the forces on the dummy to the arene carbon atoms. Angular overlap model parameters based on orbital energies derived from Kohn-Sham density functional theory (KS-DFT) calculations show that, relative to the Ru-dummy vector, the arene behaves as a very strong π donor and weak σ donor. Based on KS-DFT geometries, partial atomic charges and potential energy scans for a series of homoleptic and half sandwich complexes spanning arene, am(m)ine, imine, pyridyl, hydride and chloride ligands, a new LFMM force field has been developed which accurately reproduces the KS-DFT data. This FF was validated against 47 half-sandwich complexes obtained from the Cambridge Structural Database which, after minor corrections to account for the systematic errors between our chosen functional (BP86) and the experimental structures, yields a 'structurally tuned' FF where 93% of the Ru-L contacts are reproduced to 0.05 ? or better and all bar two bond lengths are within 0.1 ? of experiment. Over half the systems have non-hydrogen-atom rmsds of less than 0.5 ?. Larger differences are usually due to rotation of the arene moiety which is shown by ligand field molecular dynamics (LFMD) simulations to be an inherently low-energy process. Comparisons between LFMD and Car-Parrinello MD for [Ru(p-cymene)(ethylenediamine)Cl](+)show that LFMD is equally accurate but much faster enabling modelling of dynamic properties which occur on a timescale beyond the scope of CPMD.     

Linear-scaling implementation of molecular response theory in self-consistent field electronic-structure theory

A linear-scaling implementation of Hartree-Fock and Kohn-Sham self-consistent field theories for the calculation of frequency-dependent molecular response properties and excitation energies is presented, based on a nonredundant exponential parametrization of the one-electron density matrix in the atomic-orbital basis, avoiding the use of canonical orbitals. The response equations are solved iteratively, by an atomic-orbital subspace method equivalent to that of molecular-orbital theory. Important features of the subspace method are the use of paired trial vectors (to preserve the algebraic structure of the response equations), a nondiagonal preconditioner (for rapid convergence), and the generation of good initial guesses (for robust solution). As a result, the performance of the iterative method is the same as in canonical molecular-orbital theory, with five to ten iterations needed for convergence. As in traditional direct Hartree-Fock and Kohn-Sham theories, the calculations are dominated by the construction of the effective Fock/Kohn-Sham matrix, once in each iteration. Linear complexity is achieved by using sparse-matrix algebra, as illustrated in calculations of excitation energies and frequency-dependent polarizabilities of polyalanine peptides containing up to 1400 atoms.     

Tight-binding density functional theory: an approximate Kohn-Sham DFT scheme

The DFTB method is an approximate KS-DFT scheme with an LCAO representation of the KS orbitals, which can be derived within a variational treatment of an approximate KS energy functional. But it may also be related to cellular Wigner-Seitz methods and to the Harris functional. It is an approximate method, but it avoids any empirical parametrization by calculating the Hamiltonian and overlap matrices out of DFT-derived local orbitals (atomic orbitals, AO's). The method includes ab initio concepts in relating the Kohn-Sham orbitals of the atomic configuration to a minimal basis of the localized atomic valence orbitals of the atoms. Consistent with this approximation, the Hamiltonian matrix elements can strictly be restricted to a two-center representation. Taking advantage of the compensation of the so-called "double counting terms" and the nuclear repulsion energy in the DFT total energy expression, the energy may be approximated as a sum of the occupied KS single-particle energies and a repulsive energy, which can be obtained from DFT calculations in properly chosen reference systems. This relates the method to common standard "tight-binding" (TB) schemes, as they are well-known in solid-state physics. This approach defines the density-functional tight-binding (DFTB) method in its original (non-self-consistent) version.     

Four-component relativistic Kohn-Sham theory

A four-component relativistic implementation of Kohn-Sham theory for molecular systems is presented. The implementation is based on a nonredundant exponential parametrization of the Kohn-Sham energy, well suited to studies of molecular static and dynamic properties as well as of total electronic energies. Calculations are presented of the bond lengths and the harmonic and anharmonic vibrational frequencies of Au(2), Hg(2+)(2), HgAu(+), HgPt, and AuH. All calculations are based on the full four-component Dirac-Coulomb Hamiltonian, employing nonrelativistic local, gradient-corrected, and hybrid density functionals. The relevance of the Coulomb and Breit operators for the construction of relativistic functionals is discussed; it is argued that, at the relativistic level of density-functional theory and in the absence of a vector potential, the neglect of current functionals follows from the neglect of the Breit operator.     

Mimicking the two-dimensional spectrochemical series using density functional computations

With tetragonal distortions of tetrahedral d2 complexes as examples, nonadditive and additive ligand fields are compared computationally, using Kohn-Sham density functional theory (KS-DFT) and ligand-field theory to obtain 45 linear, parametrical equations. For each complex, a "data" reduction from three nonadditive-field parameters to two parameters of the additive field occurs. The complexes V and CrX4- (where X=F, Cl, Br, I) provide the two-dimensional spectrochemical series of the sigma and pi AOM parameters, which are known semiempirically for the halide ligands. The same parametrical results can be obtained from the Kohn-Sham orbital energies of the "average of configuration" computation.     

Massively parallel algorithm and implementation of RI‐MP2 energy calculation for peta‐scale many‐core supercomputers

A new parallel algorithm and its implementation for the RI‐MP2 energy calculation utilizing peta‐flop‐class many‐core supercomputers are presented. Some improvements from the previous algorithm (J. Chem. Theory Comput. 2013, 9, 5373) have been performed: (1) a dual‐level hierarchical parallelization scheme that enables the use of more than 10,000 Message Passing Interface (MPI) processes and (2) a new data communication scheme that reduces network communication overhead. A multi‐node and multi‐GPU implementation of the present algorithm is presented for calculations on a central processing unit (CPU)/graphics processing unit (GPU) hybrid supercomputer. Benchmark results of the new algorithm and its implementation using the K computer (CPU clustering system) and TSUBAME 2.5 (CPU/GPU hybrid system) demonstrate high efficiency. The peak performance of 3.1 PFLOPS is attained using 80,199 nodes of the K computer. The peak performance of the multi‐node and multi‐GPU implementation is 514 TFLOPS using 1349 nodes and 4047 GPUs of TSUBAME 2.5. © 2016 Wiley Periodicals, Inc.