Simple Modifications of the SCAN Meta-Generalized Gradient Approximation Functional

Revisiting Artifacts of Kohn-Sham Density Functionals for Biosimulation

We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn-Sham (KS) density functional theory functionals on simulation of small to medium-sized proteins in a vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO-LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the KS self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver [ Phys. Chem. Chem. Phys. 2024, 26, 6557] to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree-Fock and KS solutions; the charge delocalization not only can occur between charged fragments (such as in zwitterionic polypeptides) but also involves neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of some range-separated hybrid functionals but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping "frontier" orbitals. Caution should be exercised when unexpectedly small (or vanishing) HOMO-LUMO gaps and atypical SCF convergence patterns (e.g., oscillatory) are observed in KS DFT simulations in any context (bio or otherwise).

Solid-state performance of a meta-GGA screened hybrid density functional constructed from Pauli kinetic enhancement factor dependent semilocal exchange hole

The semilocal form of the exchange hole is highly useful in developing non-local range-separated hybrid density functionals for finite and extended systems. The way to construct the conventional exact exchange hole model is based on either the Taylor series expansion or the reverse engineering technique from the corresponding exchange energy functional. Although the latter technique is quite popular in context of generalized gradient approximation (GGA) functionals, the same for the meta-GGA functionals is not so much explored. Thus, in this study, we propose a reverse-engineered semilocal exchange hole of a meta-GGA functional, that depends only on the meta-GGA ingredient α (also known as the Pauli kinetic energy enhancement factor). The model is used subsequently to design the short-range-separated meta-GGA hybrid density functional. We show that the present method can be successfully applied for several challenging problems in the context of solids, especially for which the GGA based hybrid fails drastically. This assessment proves that the present functional is quite useful for materials sciences. Finally, we also use this method for several molecular test cases, where the results are also as comparative as its base semilocal functional.

Hirshfeld atom refinement based on projector augmented wave densities with periodic boundary conditions

Hirshfeld atom refinement (HAR) is an X-ray diffraction refinement method that, in numerous publications, has been shown to give H-atom bond lengths in close agreement with neutron diffraction derived values. Presented here is a first evaluation of an approach using densities derived from projector augmented wave (PAW) densities with three-dimensional periodic boundary conditions for HAR. The results show an improvement over refinements that neglect the crystal environment or treat it classically, while being on a par with non-periodic approximations for treating the solid-state environment quantum mechanically. A suite of functionals were evaluated for this purpose, showing that the SCAN and revSCAN functionals are most suited to these types of calculation.

New Light on an Old Story: Breaking Kasha's Rule in Phosphorescence Mechanism of Organic Boron Compounds and Molecule Design

In this work, the phosphorescence mechanism of (E)-3-(((4-nitrophenyl)imino)methyl)-2H-thiochroman-4-olate-BF2 compound (S-BF2) is investigated theoretically. The phosphorescence of S-BF2 has been reassigned to the second triplet state (T2) by the density matrix renormalization group (DMRG) method combined with the multi-configurational pair density functional theory (MCPDFT) to approach the limit of theoretical accuracy. The calculated radiative and non-radiative rate constants support the breakdown of Kasha's rule further. Our conclusion contradicts previous reports that phosphorescence comes from the first triplet state (T1). Based on the revised phosphorescence mechanism, we have purposefully designed some novel compounds in theory to enhance the phosphorescence efficiency from T2 by replacing substitute groups in S-BF2. Overall, both S-BF2 and newly designed high-efficiency molecules exhibit anti-Kasha T2 phosphorescence instead of the conventional T1 emission. This work provides a useful guidance for future design of high-efficiency green-emitting phosphors.

Benchmark test of a dispersion corrected revised Tao-Mo semilocal functional for thermochemistry, kinetics, and noncovalent interactions of molecules and solids

In the density functional theory, dispersion corrected semilocal approximations are often used to benchmark weekly interacting finite and extended systems. Here, the focus is on providing a broad overview of the performance of D3 dispersion corrected revised Tao-Mo (revTM) semilocal functionals [A. Patra et al., J. Chem. Phys. 153, 084 117 (2020)] for thermochemistry and kinetics of molecules, molecular crystals, ice polymorphs, metal-organic systems, atom/molecular adsorption on solids, water interacting with nano-materials, binding energies of layered materials, and properties of weekly and strongly bonded solids. We show that the most suitable "optimized power" function for the revTM functional needs a modification to make it suitable for properties related to the diverse nature of finite and extended systems. The present work is an extension of the previously proposed revTM+D3 method with the motivation to design and benchmark the dispersion corrected cost-effective method based on this semilocal approximation. We show that the revised revTM+D3 functional provides various general purpose molecular and solid properties with the closest to experimental findings than its predecessor. The present assessment and benchmarking can be practically useful for performing cost-effective method based simulations of various molecular and solid-state properties.

Accurate density functional made more versatile

We propose a one-electron self-interaction-free correlation energy functional compatible with the order-of-limit problem-free Tao-Mo (TM) semilocal functional (regTM) [J. Tao and Y. Mo, Phys. Rev. Lett. 117, 073001 (2016) and Patra et al., J. Chem. Phys. 153, 184112 (2020)] to be used for general purpose condensed matter physics and quantum chemistry. The assessment of the proposed functional for large classes of condensed matter and chemical systems shows its improvement in most cases compared to the TM functional, e.g., when applied to the relative energy difference of MnO₂ polymorphs. In this respect, the present exchange-correction functional, which incorporates the TM technique of the exchange hole model combined with the slowly varying density correction, can achieve broad applicability, being able to solve difficult solid-state problems.

Performance of the DLPNO-CCSD and recent DFT methods in the calculation of isotropic and dipolar contributions to 14N hyperfine coupling constants of nitroxide radicals

In the present study, the performance of a set of density functionals: BP86, PBE, OLYP, BEEF, PBEpow, TPSS, SCAN, PBEGXPBE, M06L, MN15L, B3LYP, PBE0, mPW1PW, B97, BHandHLYP, mPW1PW, B98, TPSS0, PBE1KCIS, SCAN0, M06, M06-2X, MN15, CAM-B3LYP, ωB97x, B2PLYP, and the B3LYP/N07D and PBE/N07D schemes in the calculation of the ¹⁴N anisotropic hyperfine coupling (HFC) constants of a set of 23 nitroxide radicals is evaluated. The results are compared with those obtained with the DLPNO-CCSD method and experimental HFC values. Harmonic contribution to the ¹⁴N HFC vibrational correction was calculated at the revPBE0/def2-TZVPP level and included in the evaluation. With the vibrational correction, the DLPNO-CCSD method yielded HFC values in good agreement with the experiment (mean absolute deviation (MAD) = 0.3 G for the dipole-dipole contribution and MAD = 0.8 G for the contact coupling contribution). The best DFT results are obtained using the M06 functional with MAD = 0.2 G for the dipole-dipole contribution and MAD = 0.7 G for the contact coupling contribution. In general, vibrational correction significantly improved most DFT functionals' performance but did not change its overall ranking.

Describing adsorption of benzene, thiophene, and xenon on coinage metals by using the Zaremba-Kohn theory-based model

Semilocal (SL) density functional approximations (DFAs) are widely applied but have limitations due to their inability to incorporate long-range van der Waals (vdW) interaction. Non-local functionals (vdW-DF, VV10, and rVV10) or empirical methods (DFT+D, DFT+vdW, and DFT+MBD) are used with SL-DFAs to account for such missing interaction. The physisorption of a molecule on the surface of the coinage metals (Cu, Ag, and Au) is a typical example of systems where vdW interaction is significant. However, it is difficult to find a general method that reasonably describes both adsorption energy and geometry of even the simple prototypes of cyclic and heterocyclic aromatic molecules such as benzene (C₆H₆) and thiophene (C₄H₄S), respectively, with reasonable accuracy. In this work, we present an alternative scheme based on Zaremba-Kohn theory, called DFT+vdW-dZK. We show that unlike other popular methods, DFT+vdW-dZK and particularly SCAN+vdW-dZK give an accurate description of the physisorption of a rare-gas atom (xenon) and two small albeit diverse prototype organic molecules on the (111) surfaces of the coinage metals.

A computational study of the properties of low- and high-index Pd, Cu and Zn surfaces

We report a detailed Density Functional Theory (DFT) based investigation of the structure and stability of bulk and surface structures for the Group 10-12 elements Pd, Cu and Zn, considering the effect of the choice of exchange-correlation density functional and computation parameters. For the initial bulk structures, the lattice parameter and cohesive energy are calculated, which are then augmented by calculation of surface energies and work functions for the lower-index surfaces. Of the 22 density functionals considered, we highlight the mBEEF density functional as providing the best overall agreement with experimental data. The optimal density functional choice is applied to the study of higher index surfaces for the three metals, and Wulff constructions performed for nanoparticles with a radius of 11 nm, commensurate with nanoparticle sizes commonly employed in catalytic chemistry. For Pd and Cu, the low-index (111) facet is dominant in the constructed nanoparticles, covering ∼50% of the surface, with (100) facets covering a further 10 to 25%; however, non-negligible coverage from higher index (332), (332) and (210) facets is also observed for Pd, and (322), (221) and (210) surfaces are observed for Cu. In contrast, only the (0001) and (10-10) facets are observed for Zn. Overall, our results highlight the need for careful validation of computational settings before performing extensive density functional theory investigations of surface properties and nanoparticle structures of metals.

Accurate and Numerically Efficient r2SCAN Meta-Generalized Gradient Approximation

The recently proposed rSCAN functional [ J. Chem. Phys. 2019 150, 161101] is a regularized form of the SCAN functional [ Phys. Rev. Lett. 2015 115, 036402] that improves SCAN's numerical performance at the expense of breaking constraints known from the exact exchange-correlation functional. We construct a new meta-generalized gradient approximation by restoring exact constraint adherence to rSCAN. The resulting functional maintains rSCAN's numerical performance while restoring the transferable accuracy of SCAN.

Assessing the effect of regularization on the molecular properties predicted by SCAN and self-interaction corrected SCAN meta-GGA

Recent regularization of the SCAN meta-GGA functional (rSCAN) has simplified the numerical complexities of the SCAN functional, alleviating SCAN's stringent demand on the numerical integration grids to some extent. The regularization of rSCAN, however, results in the breaking of some constraints such as the uniform electron gas limit, the slowly varying density limit, and coordinate scaling of the iso-orbital indicator. Here, we assess the effects of regularization on the electronic, structural, vibrational, and magnetic properties of molecules by comparing the SCAN and rSCAN predictions. The properties studied include atomic energies, atomization energies, ionization potentials, electron affinities, barrier heights, infrared intensities, dissociation and reaction energies, spin moments of molecular magnets, and isomer ordering of water clusters. Our results show that rSCAN requires less dense numerical grids and gives very similar results to those of SCAN for all properties examined with the exception of atomization energies, which are worsened in rSCAN. We also examine the performance of self-interaction-corrected (SIC) rSCAN with respect to SIC-SCAN using the Perdew-Zunger (PZ) SIC method. The PZSIC method uses orbital densities to compute one-electron self-interaction errors and places an even more stringent demand on numerical grids. Our results show that SIC-rSCAN gives marginally better performance than SIC-SCAN for almost all properties studied in this work with numerical grids that are on average half or less as dense as that needed for SIC-SCAN.

First-Principles Calculation of Transition Metal Hyperfine Coupling Constants with the Strongly Constrained and Appropriately Normed (SCAN) Density Functional and its Hybrid Variants

Density functional theory (DFT) is used extensively for the first-principles calculation of hyperfine coupling constants in both main-group and transition metal systems. As with many other properties, the performance of DFT for hyperfine coupling constants is of variable quality, particularly for transition metal complexes, because it strongly depends on the nature of the chemical system and the type of approximation to the exchange-correlation functional. Recently, a meta-generalized-gradient approximation (mGGA) functional was proposed that obeys all known exact constraints for such a method, known as the Strongly Constrained and Appropriately Normed (SCAN) functional. In view of its theoretically superior formulation a benchmark set of complexes is used to assess the performance of SCAN for the challenging case of transition metal hyperfine coupling constants. In addition, two global hybrid versions of the functional, SCANh and SCAN0, are described and tested. The values computed with the new functionals are compared with experiment and with those of other DFT approximations. Although the original SCAN and the SCAN-based hybrids may offer improved hyperfine coupling constants for specific systems, no uniform improvement is observed. On the contrary, there are specific cases where the new functionals fail badly due to a flawed description of the underlying electronic structure. Therefore, despite these methodological advances, systematically accurate and system-independent prediction of transition metal hyperfine coupling constants with DFT remains an unmet challenge.

Origin of structural stability of ScH3 molecular nanowires and their chemical-bonding behavior: Correlation effects of the Sc 3d electrons

A new stable transition-metal trihydride (ScH₃) molecular nanowire was recently reported by Li et al. [J. Am. Chem. Soc. 139, 6290-6293 (2017)]. Of the two typical structures (T-ScH₃ and O-ScH₃), T-ScH₃ is more stable than O-ScH₃. However, the reason why O-ScH₃ is less stable than T-ScH₃ was not known. Using Perdew-Burke-Ernzerhof (PBE), PBE+U, SCAN, and HSE06, as well as crystal orbital Hamilton populations (COHPs), we investigate the orbital-projected band structures and chemical bonding of T-ScH₃ and O-ScH₃. It is found that the energies calculated by PBE, SCAN, and HSE06 indeed reveal that T-ScH₃ is more stable than O-ScH₃, and there is no occupied antibonding state at the Fermi level of the COHP curves of T-ScH₃, supporting the stable Sc-H bonding of T-ScH₃. To the contrary, the Sc-H bonding of O-ScH₃ is unstable because there exist occupied antibonding states at the Fermi level of the COHP curves of O-ScH₃. We found that the results of PBE+U are consistent with those of PBE, SCAN, and HSE06 in the case of U < U_c. However, when U > U_c, the results of PBE+U are opposite to those of PBE, SCAN, and HSE06.