skip to main content


Title: Molecular orbital projectors in non-empirical jmDFT recover exact conditions in transition-metal chemistry
Low-cost, non-empirical corrections to semi-local density functional theory are essential for accurately modeling transition-metal chemistry. Here, we demonstrate the judiciously modified density functional theory (jmDFT) approach with non-empirical U and J parameters obtained directly from frontier orbital energetics on a series of transition-metal complexes. We curate a set of nine representative Ti(III) and V(IV) d 1 transition-metal complexes and evaluate their flat-plane errors along the fractional spin and charge lines. We demonstrate that while jmDFT improves upon both DFT+U and semi-local DFT with the standard atomic orbital projectors (AOPs), it does so inefficiently. We rationalize these inefficiencies by quantifying hybridization in the relevant frontier orbitals. To overcome these limitations, we introduce a procedure for computing a molecular orbital projector (MOP) basis for use with jmDFT. We demonstrate this single set of d 1 MOPs to be suitable for nearly eliminating all energetic delocalization and static correlation errors. In all cases, MOP jmDFT outperforms AOP jmDFT, and it eliminates most flat-plane errors at non-empirical values. Unlike DFT+U or hybrid functionals, jmDFT nearly eliminates energetic delocalization and static correlation errors within a non-empirical framework.  more » « less
Award ID(s):
1846426
NSF-PAR ID:
10404794
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
156
Issue:
18
ISSN:
0021-9606
Page Range / eLocation ID:
184112
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Transition metal oxides such as BiVO 4 are promising photoelectrode materials for solar-to-fuel conversion applications. However, their performance is limited by the low carrier mobility (especially electron mobility) due to the formation of small polarons. Recent experimental studies have shown improved carrier mobility and conductivity by atomic doping; however the underlying mechanism is not understood. A fundamental atomistic-level understanding of the effects on small polaron transport is critical to future material design with high conductivity. We studied the small polaron hopping mobility in pristine and doped BiVO 4 by combining Landau–Zener theory and kinetic Monte Carlo (kMC) simulation fully from first-principles, and investigated the effect of dopant–polaron interactions on the mobility. We found that polarons are spontaneously formed at V in both pristine and Mo/W doped BiVO 4 , which can only be described correctly by density functional theory (DFT) with the Hubbard correction (DFT+U) or hybrid exchange-correlation functional but not local or semi-local functionals. We found that DFT+U and dielectric dependant hybrid (DDH) functionals give similar electron hopping barriers, which are also similar between the room temperature monoclinic phase and the tetragonal phase. The calculated electron mobility agrees well with experimental values, which is around 10 −4 cm 2 V −1 s −1 . We found that the electron polaron transport in BiVO 4 is neither fully adiabatic nor nonadiabatic, and the first and second nearest neighbor hoppings have significantly different electronic couplings between two hopping centers that lead to different adiabaticity and prefactors in the charge transfer rate, although they have similar hopping barriers. Without considering the detailed adiabaticity through Landau–Zener theory, one may get qualitatively wrong carrier mobility. We further computed polaron mobility in the presence of different dopants and showed that Cr substitution of V is an electron trap while Mo and W are “repulsive” centers, mainly due to the minimization of local lattice expansion by dopants and electron polarons. The dopants with “repulsive” interactions to polarons are promising for mobility improvement due to larger wavefunction overlap and delocalization of locally concentrated polarons. 
    more » « less
  2. Density functional theory (DFT) is widely used in transition-metal chemistry, yet essential properties such as spin-state energetics in transition-metal complexes (TMCs) are well known to be sensitive to the choice of the exchange-correlation functional. Increasing the amount of exchange in a functional typically shifts the preferred ground state in first-row TMCs from low-spin to high-spin by penalizing delocalization error, but the effect on properties of second-row complexes is less well known. We compare the exchange sensitivity of adiabatic spin-splitting energies in pairs of mononuclear 3d and 4d mid-row octahedral transition metal complexes. We analyze hundreds of complexes assembled from four metals in two oxidation states with ten small monodentate ligands that span a wide range of field strengths expected to favor a variety of ground states. We observe consistently lower but proportional sensitivity to exchange fraction among 4d TMCs with respect to their isovalent 3d TMC counterparts, leading to the largest difference in sensitivities for the strongest field ligands. The combined effect of reduced exchange sensitivities and the greater low-spin bias of most 4d TMCs means that while over one-third of 3d TMCs change ground states over a modest variation (ca. 0.0–0.3) in exchange fraction, almost no 4d TMCs do. Differences in delocalization, as judged through changes in the metal–ligand bond lengths of spin states, do not explain the distinct behavior of 4d TMCs. Instead, evaluation of potential energy curves in 3d and 4d TMCs reveals that higher exchange sensitivities in 3d TMCs are likely due to the opposing effect of exchange on the low-spin and high-spin states, whereas the effect on both spin states is more comparable in 4d TMCs. 
    more » « less
  3. Here we present the synthesis and characterization of a hybrid vanadium-organic coordination polymer with robust magnetic order, a Curie temperature T C of ∼110 K, a coercive field of ∼5 Oe at 5 K, and a maximum mass magnetization of about half that of the benchmark ferrimagnetic vanadium(tetracyanoethylene) ∼2 (V·(TCNE) ∼2 ). This material was prepared using a new tetracyano-substituted quinoidal organic small molecule 7 based on a tricyclic heterocycle 4-hexyl-4 H -pyrrolo[2,3- d :5,4- d ′]bis(thiazole) ( C6-PBTz ). Single crystal X-ray diffraction of the 2,6-diiodo derivative of the parent C6-PBTz , showed a disordered hexyl chain and a nearly linear arrangement of the substituents in positions 2 and 6 of the tricyclic core. Density functional theory (DFT) calculations indicate that C6-PBTz -based ligand 7 is a strong acceptor with an electron affinity larger than that of TCNE and several other ligands previously used in molecular magnets. This effect is due in part to the electron-deficient thiazole rings and extended delocalization of the frontier molecular orbitals. The ligand detailed in this study, a representative example of fused heterocycle aromatic cores with extended π conjugation, introduces new opportunities for structure–magnetic-property correlation studies where the chemistry of the tricyclic heterocycles can modulate the electronic properties and the substituent at the central N -position can vary the spatial characteristics of the magnetic polymer. 
    more » « less
  4. In order to shed light on metal-dependent mechanisms for O–O bond cleavage, and its microscopic reverse, we compare herein the electronic and geometric structures of O2-derived binuclear Co(III)– and Mn(III)–peroxo compounds. Binuclear metal peroxo complexes are proposed to form as intermediates during Mn-promoted photosynthetic H2O oxidation, as well as a Co-containing artificial leaf inspired by nature’s photosynthetic H2O oxidation catalyst. Crystallographic characterization of an extremely activated peroxo is made possible by working with substitution-inert, low-spin Co(III). Density functional theory (DFT) calculations show that the frontier orbitals of the Co(III)–peroxo compound differ noticeably from the analogous Mn(III)–peroxo compound. The highest occupied molecular orbital (HOMO) associated with the Co(III)–peroxo is more localized on the peroxo in an antibonding π*(O–O) orbital, whereas the HOMO of the structurally analogous Mn(III)–peroxo is delocalized over both the metal d-orbitals and peroxo π*(O–O) orbital. With low-spin d6 Co(III), filled t2g orbitals prevent π-back-donation from the doubly occupied antibonding π*(O–O) orbital onto the metal ion. This is not the case with high-spin d4 Mn(III), since these orbitals are half-filled. This weakens the peroxo O–O bond of the former relative to the latter. 
    more » « less
  5. null (Ed.)
    Accurate computational predictions of band gaps are of practical importance to the modeling and development of semiconductor technologies, such as (opto)electronic devices and photoelectrochemical cells. Among available electronic-structure methods, density-functional theory (DFT) with the Hubbard U correction (DFT+U) applied to band edge states is a computationally tractable approach to improve the accuracy of band gap predictions beyond that of DFT calculations based on (semi)local functionals. At variance with DFT approximations, which are not intended to describe optical band gaps and other excited-state properties, DFT+U can be interpreted as an approximate spectral-potential method when U is determined by imposing the piecewise linearity of the total energy with respect to electronic occupations in the Hubbard manifold (thus removing self-interaction errors in this subspace), thereby providing a (heuristic) justification for using DFT+U to predict band gaps. However, it is still frequent in the literature to determine the Hubbard U parameters semiempirically by tuning their values to reproduce experimental band gaps, which ultimately alters the description of other total-energy characteristics. Here, we present an extensive assessment of DFT+U band gaps computed using self-consistent ab initio U parameters obtained from density-functional perturbation theory to impose the aforementioned piecewise linearity of the total energy. The study is carried out on 20 compounds containing transition-metal or p-block (group III-IV) elements, including oxides, nitrides, sulfides, oxynitrides, and oxysulfides. By comparing DFT+U results obtained using nonorthogonalized and orthogonalized atomic orbitals as Hubbard projectors, we find that the predicted band gaps are extremely sensitive to the type of projector functions and that the orthogonalized projectors give the most accurate band gaps, in satisfactory agreement with experimental data. This work demonstrates that DFT+U may serve as a useful method for high-throughput workflows that require reliable band gap predictions at moderate computational cost. 
    more » « less