Density functional theory with alternative spin densities: application to magnetic systems with localized spins

Hybrid Density Functional Valence Bond Method with Multistate Treatment

Recently, a hybrid density functional valence bond (VB) method, λ-DFVB(U), has been proposed and shown to give accuracy that is comparable to that of CASPT2 in calculations of atomization energies, atomic excitation energies, and reaction barriers, while its computational cost is approximately the same as the valence bond self-consistent-field (VBSCF) method. However, the interaction between electronic states is not included in λ-DFVB(U) since the last step of λ-DFVB(U) is not a diagonalization of the Hamiltonian matrix on the electronic state basis. Therefore, λ-DFVB(U) gives the wrong topology of the potential energy surfaces (PESs) near the conical intersection region. In the present paper, we propose a novel hybrid density functional VB method with multistate treatment, named λ-DFVB(MS), in which an effective Hamiltonian matrix is constructed on the basis of the diabatic states obtained by the valence-bond-based compression approach for the diabatization scheme, and the interaction between electronic states can be included through the diagonalization of the effective Hamiltonian matrix. Test calculations show that λ-DFVB(MS) gives the correct topology of the PESs near the conical intersection region. We also show that the VBSCF wave function with selected VB structures can be applied as a reference in λ-DFVB(MS).

λ-DFVB(U): A hybrid density functional valence bond method based on unpaired electron density

In this paper, a hybrid density functional valence bond method based on unpaired electron density, called λ-DFVB(U), is presented, which is a combination of the valence bond self-consistent field (VBSCF) method and Kohn–Sham density functional theory. In λ-DFVB(U), the double-counting error of electron correlation is mitigated by a linear decomposition of the electron–electron interaction using a parameter λ, which is a function of an index based on the number of effectively unpaired electrons. In addition, λ-DFVB(U) is based on the approximation that correlation functionals in KS-DFT only cover dynamic correlation and exchange functionals mimic some amount of static correlation. Furthermore, effective spin densities constructed from unpaired density are used to address the symmetry dilemma problem in λ-DFVB(U). The method is applied to test calculations of atomization energies, atomic excitation energies, and reaction barriers. It is shown that the accuracy of λ-DFVB(U) is comparable to that of CASPT2, while its computational cost is approximately the same as VBSCF.

Spin-Flip Pair-Density Functional Theory: A Practical Approach To Treat Static and Dynamical Correlations in Large Molecules

We present a practical approach to treat static and dynamical correlation accurately in large multiconfigurational systems. The static correlation is taken into account by using the spin-flip approach, which is well-known for capturing static correlation accurately at low-computational expense. Unlike previous approaches to add dynamical correlation to spin-flip models which use perturbation theory or coupled-cluster theory, we explore the ability to use the on-top pair-density functional theory approaches recently developed by Gagliardi and co-workers (J. Comput. Theor. Chem., 2014, 10, 3669). External relaxations are performed in the spin-flip calculations through a restricted active space framework for which a truncation scheme for the orbitals used in the external excitation is presented. The performance of the approach is demonstrated by computing energy gaps between ground and excited states for diradicals, triradicals, and linear polyacene chains ranging from naphthalene to dodecacene. Accurate results are obtained using the new approach for these challenging open-shell molecular systems.

Short-range DFT energy correction to multiconfigurational wave functions for open-shell systems

Electronic structure methods emerging from the combination of multiconfigurational wave functions and density functional theory (DFT) aim to take advantage of the strengths of the two nearly antagonistic theories. One of the common strategies employed to merge wave function theory (WFT) with DFT relies on the range separation of the Coulomb operator in which DFT functionals take care of the short-distance part, while long-range inter-electronic interactions are evaluated by using the chosen wave function method (WFT-srDFT). In this work, we uncover the limitations of WFT-srDFT in the characterization of open-shell systems. We show that spin polarization effects have a major impact on the (short-range) DFT exchange energy and are of vital importance in order to provide a balanced description between closed and open-shell configurations. We introduce different strategies to account for spin polarization in the short range based on the definition of a spin polarized electron density and with the use of short-range exact exchange. We test the performance of these approaches in the dissociation of the hydrogen molecule, the calculation of energy gaps in spin-triplet atoms and molecular diradicals, and the characterization of low-lying states of the gallium dimer. Our results indicate that the use of short-range DFT correlation in combination with a (full-range) multiconfigurational wave function might be an excellent approach for the study of open-shell molecules and largely improves the performance of WFT and WFT-srDFT.

A Valence-Bond-Based Multiconfigurational Density Functional Theory: The λ-DFVB Method Revisited

A recently developed valence-bond-based multireference density functional theory, named λ-DFVB, is revisited in this paper. λ-DFVB remedies the double-counting error of electron correlation by decomposing the electron-electron interactions into the wave function term and density functional term with a variable parameter λ. The λ value is defined as a function of the free valence index in our previous scheme, denoted as λ-DFVB(K) in this paper. Here we revisit the λ-DFVB method and present a new scheme based on natural orbital occupation numbers (NOONs) for parameter λ, named λ-DFVB(IS), to simplify the process of λ-DFVB calculation. In λ-DFVB(IS), the parameter λ is defined as a function of NOONs, which are straightforwardly determined from the many-electron wave function of the molecule. Furthermore, λ-DFVB(IS) does not involve further self-consistent field calculation after performing the valence bond self-consistent field (VBSCF) calculation, and thus, the computational effort in λ-DFVB(IS) is approximately the same as the VBSCF method, greatly reduced from λ-DFVB(K). The performance of λ-DFVB(IS) was investigated on a broader range of molecular properties, including equilibrium bond lengths and dissociation energies, atomization energies, atomic excitation energies, and chemical reaction barriers. The computational results show that λ-DFVB(IS) is more robust without losing accuracy and comparable in accuracy to high-level multireference wave function methods, such as CASPT2.

λ-Density Functional Valence Bond: A Valence Bond-Based Multiconfigurational Density Functional Theory With a Single Variable Hybrid Parameter

A new valence bond (VB)-based multireference density functional theory (MRDFT) method, named λ-DFVB, is presented in this paper. The method follows the idea of the hybrid multireference density functional method theory proposed by Sharkas et al. (2012). λ-DFVB combines the valence bond self-consistent field (VBSCF) method with Kohn–Sham density functional theory (KS-DFT) by decomposing the electron–electron interactions with a hybrid parameter λ. Different from the Toulouse's scheme, the hybrid parameter λ in λ-DFVB is variable, defined as a function of a multireference character of a molecular system. Furthermore, the E_C correlation energy of a leading determinant is introduced to ensure size consistency at the dissociation limit. Satisfactory results of test calculations, including potential energy surfaces, bond dissociation energies, reaction barriers, and singlet–triplet energy gaps, show the potential capability of λ-DFVB for molecular systems with strong correlation.

Combining Wave Function Methods with Density Functional Theory for Excited States

We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.

Hamiltonian Matrix Correction Based Density Functional Valence Bond Method

In this work, a valence bond type multireference density functional theory (MRDFT) method, called the Hamiltonian matrix correction based density functional valence bond method (hc-DFVB), is presented. In hc-DFVB, the static electronic correlation is considered by the valence bond self-consistent field (VBSCF) strategy, while the dynamic correlation energy is taken into account by Kohn-Sham density functional theory (KS-DFT). Different from our previous version of DFVB (J. Chem. Theory Comput. 2012, 8, 1608), hc-DFVB corrects the dynamic correlation energy with a Hamiltonian correction matrix, improving the functional adaptability and computational accuracy. The method was tested for various physical and chemical properties, including spectroscopic constants, bond dissociation energies, reaction barriers, and singlet-triplet gaps. The accuracy of hc-DFVB matches that of KS-DFT and high level molecular orbital (MO) methods quite well. Furthermore, hc-DFVB keeps the advantages of VB methods, which are able to provide clear interpretations and chemical insights with compact wave functions.

DFVB: A Density-Functional-Based Valence Bond Method

A new ab initio valence bond method with density-functional-based correlation correction, so-called DFVB, is presented. In the DFVB method, the dynamic correlation energy is taken into account by use of density correlation functional(s), while the static correlation energy is covered by the VBSCF wave function. Owing to incorporation of DFT methods, DFVB provides an economic route to improving the accuracy of ab initio VB theory. Various tests of the method are presented, including the spectroscopic parameters of a series of diatomic molecules, the dipole moments of the NF molecule for different electronic states, and the singlet-triplet gaps of the diradical species, chemical reactions barriers, and total charge-shift resonance energies. These tests show that DFVB is capable of providing high accuracy with relatively low computational cost by comparison to the currently existing post-VBSCF methods.

Combining two-body density correlation functionals with multiconfigurational wave functions using natural orbitals and occupation numbers

We propose a procedure that combines multiconfigurational (MC) wave functions with two-body density correlation functionals by transforming the latter into functionals of the MC natural orbitals and occupation numbers. The method is tested with the spectroscopic constants of a set of 11 diatomics, the diradical-involved automerization barrier of cyclobutadiene, the energy difference between triplet and open-shell singlet states in He and the methylene molecule, and the magnetic coupling constants of several systems, such as NiO, KNiF(3), K(2)NiF(4), La(2)CuO(4), alpha-4-dehydrotoluene, 1,1('),5,5(')-tetramethyl-6,6(')-dioxo-3,3(')-biverdazyl, [Cu(2)Cl(6)](-2), copper(II) acetate monohidrate and H-He-H. The procedure is applied to the Colle-Salvetti [Theor. Chim. Acta 37, 329 (1975); 53, 55 (1979)], functional and to a size-consistent functional depending on the on-top pair density (F1-5-N(eff)). On average, the best results are provided by the transformed F1-5-N(eff) [J. Chem. Phys. 114, 2022 (2001)] functional.

Merging multiconfigurational wavefunctions and correlation functionals to predict magnetic coupling constants

We study the performance of different approaches that combine multiconfigurational wavefunctions with correlation functionals for the calculation of magnetic coupling constants of several materials and molecules. The systems under study include four antiferromagnetic materials: NiO, KNiF(3), K(2)NiF(4) and La(2)CuO(4); two biradicals: alpha-4-Dehydrotoluene and 1,1',5,5'-Tetramethyl-6,6'-dioxo-3,3'-biverdazyl; two molecular complexes: [Cu(2)Cl(6)](-2) and Copper(II) acetate monohidrate; and the prototypical H-He-H system. On average, the best results are obtained with a recently proposed method [Phys. Rev. A 75, 012503 (2007)] that estimates the correlation energy of density functionals from a pair of alternative spin densities built from the natural orbitals and occupation numbers of the multiconfigurational wavefunction.

Stable Tin (n = 2−15) Clusters and Their Geometries: DFT Calculations

We present a detailed structural analysis for small Tin (n = 2-15) clusters based on ab initio quantum mechanical calculations of their binding energies, frontier orbital gaps, and second energy derivatives. Local density approximation calculations revealed that while the smaller clusters (n < or = 8) prefer hexagonal atomic arrays with bulklike crystal symmetry, the bigger clusters prefer pentagonal atomic arrays. From the stability criteria of the magic number clusters we could identify three magic number clusters Ti7, Ti13, and Ti15. While the most stable configuration of Ti7 is a decahedral bipyramid induced by tetrahedral atomic arrays, the most stable configuration of Ti13 is an icosahedron. The other stable cluster Ti15 takes a closed-shell icosahedron-like configuration with both pentagonal and hexagonal symmetries. The stability of the Tin clusters strongly depends on their geometries and charge states. The HOMO-LUMO gap of the Tin clusters approaches its bulk value for n > 8. While there is not much difference between the HOMO and LUMO isosurface charge distributions for the Ti7 and Ti13 clusters in their most stable configurations, they are very different in the case of Ti15. Such a distinct charge distribution in Ti15 indicates its singular chemical selectivity over the other two magic number clusters.

Performance of the τ-dependent functionals in predicting the magnetic coupling of ionic antiferromagnetic insulators

The performance of some kinetic energy density (tau) dependent functionals in predicting the effective Heisenberg exchange has been explored using the KNiF3 and K2NiF4 insulators as case examples. Our results show that this new generation of functionals represents an important improvement with respect to the current local and gradient corrected functionals yielding a semi-quantitative description of the antiferromagnetic coupling without the need of hybrid approaches thus avoiding the calculation of exact, Hartree-Fock exchange. This feature opens a wide field of application especially in solid state.