We study the effect of self-interaction errors on the barrier heights of chemical reactions. For this purpose, we use the well-known Perdew–Zunger self-interaction-correction (PZSIC) [J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981)] as well as two variations of the recently developed, locally scaled self-interaction correction (LSIC) [Zope et al., J. Chem. Phys. 151, 214108 (2019)] to study the barrier heights of the BH76 benchmark dataset. Our results show that both PZSIC and especially the LSIC methods improve the barrier heights relative to the local density approximation (LDA). The version of LSIC that uses the iso-orbital indicator z as a scaling factor gives a more consistent improvement than an alternative version that uses an orbital-dependent factor w based on the ratio of orbital densities to the total electron density. We show that LDA energies evaluated using the self-consistent and self-interaction-free PZSIC densities can be used to assess density-driven errors. The LDA reaction barrier errors for the BH76 set are found to contain significant density-driven errors for all types of reactions contained in the set, but the corrections due to adding SIC to the functional are much larger than those stemming from the density for the hydrogen transfer reactions and of roughly equal size for the non-hydrogen transfer reactions.

Density functional theory (DFT) suffers from selfinteraction errors (SIEs) that generally result in the underestimation of chemical reaction barrier heights. This is commonly attributed to the tendency of density functional approximations to overstabilize delocalized densities that typically occur in the stretched bonds of transition state structures. The Perdew−Zunger self-interaction correction (PZSIC) and locally scaled selfinteraction correction (LSIC) improve the prediction of barrier heights of chemical reactions, with LSIC giving better accuracy than PZSIC on average. These methods employ an orbital-byorbital correction scheme to remove the one-electron SIE. In the context of barrier heights, this allows an analysis of how the selfinteraction correction (SIC) for each orbital contributes to the calculated barriers using Fermi−Loẅdin orbitals (FLOs). We hypothesize that the SIC contribution to the reaction barrier comes mainly from a limited number of orbitals that are directly involved in bond-breaking and bond-making in the reaction transition state. We call these participant orbitals (POs), in contrast to spectator orbitals (SOs) which are not directly involved in changes to the bonding. Our hypothesis is that ΔE Total SIC ≈ ΔE PO SIC , where ΔE Total SIC is the difference in the SIC corrections for the reactants or products and the transition state. We test this hypothesis for the reaction barriers of the BH76 benchmark set of reactions. We find that the stretched-bond orbitals indeed make the largest individual SIC contributions to the barriers. These contributions increase the barrier heights relative to LSDA, which underpredicts the barrier. However, the full stretched-bond hypothesis does not hold in all cases for either PZSIC or LSIC. There are many cases where the total SIC contribution from the SOs is significant and cannot be ignored. The size of the SIC contribution to the barrier height is a key indicator. A large SIC correction is correlated to a large LSDA error in the barrier, showing that PZSIC properly gives larger corrections when corrections are needed most. A comparison of the performance of PZSIC and LSIC shows that the two methods have similar accuracy for reactions with large LSDA errors, but LSIC is clearly better for reactions with small errors. We trace this to an improved description of reaction energies in LSIC.

We examine the role of self-interaction error (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular, we analyze the effect of scaling down the selfinteraction correction (SIC) for three nonempirical density functional approximations (DFAs) namely, the local spin density approximation, the Perdew−Burke−Ernzerhof generalized gradient approximation, and the recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods:

Recently, pressure has been used to induce structural and magnetic phase transitions in many layered quantum materials whose layers are linked by van der Waals forces. Materials with such weakly held layers allow for relatively easy manipulation of the superexchange mechanism, which can give rise to novel magnetic behavior. Using hydrostatic pressure as a disorderless means to manipulate the interlayer coupling, we applied pressure on two quasi-2D sister compounds, namely, Cr 2 Si 2 Te 6 (CST) and Mn 3 Si 2 Te 6 (MST), up to ∼1 GPa. Magnetic property measurements with the application of pressure revealed that the ferromagnetic transition temperature decreases in CST, while the opposite occurs for the ferrimagnetic MST. In MST, magnetization decreases with the increase in pressure, and such a trend is not clearly observed within the pressure range studied for CST. The overall pressure effect on magnetic characteristics such as exchange couplings and magnetic anisotropy energies is also examined theoretically using density functional theory. Exchange coupling in MST is strongly frustrated, and the first nearest neighbor interaction is the most dominant of the components with the strongest pressure dependence. In CST, the exchange coupling parameters exhibit very little dependence on pressure. This combined experimental and theoretical work has the potential to expand to other relevant quantum materials.