An official website of the United States government
Modeling L-edge spectra at X-ray wavelengths requires consideration of spin–orbit splitting of the 2p orbitals. We introduce a low-cost tool to compute core-level spectra that combines a spin–orbit mean-field description of the Breit–Pauli Hamiltonian with nonrelativistic excited states computed using the semi-empirical density-functional theory configuration-interaction singles (DFT/CIS) method, within the state-interaction approach. Our version of DFT/CIS was introduced recently for K-edge spectra and includes a semi-empirical correction to the core orbital energies, significantly reducing ad hoc shifts that are typically required when time-dependent (TD-)DFT is applied to core-level excitations. In combination with the core/valence separation approximation and spin–orbit couplings, the DFT/CIS method affords semiquantitative L-edge spectra at CIS cost. Spin–orbit coupling has a qualitative effect on the spectra, as demonstrated for a variety of 3d transition metal systems and main-group compounds. The use of different active orbital spaces helps to facilitate spectral assignments. We find that spin–orbit splitting has a negligible effect on M-edge spectra for 3d transition metal species.
more »
« less
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) d1 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 d1 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
- PAR ID:
- 10593490
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 156
- Issue:
- 18
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
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
-
We reformulate the thermally assisted-occupation density functional theory (TAO-DFT) into the Kohn–Sham single-determinant framework and construct two new post-self-consistent field (post-SCF) static correlation correction schemes, named rTAO and rTAO-1. In contrast to the original TAO-DFT with the density in an ensemble form, in which each orbital density is weighted with a fractional occupation number, the ground-state density is given by a single-determinant wavefunction, a regular Kohn–Sham (KS) density, and total ground state energy is expressed in the normal KS form with a static correlation energy formulated in terms of the KS orbitals. In post-SCF calculations with rTAO functionals, an efficient energy scanning to quantitatively determine θ is also proposed. The rTAOs provide a promising method to simulate systems with strong static correlation as original TAO, but simpler and more efficient. We show that both rTAO and rTAO-1 is capable of reproducing most results from TAO-DFT without the additional functional Eθ used in TAO-DFT. Furthermore, our numerical results support that, without the functional Eθ, both rTAO and rTAO-1 can capture correct static correlation profiles in various systems.more » « less
-
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
