Paper detail

Local self-interaction correction method with a simple scaling factor

A recently proposed local self-interaction correction (LSIC) method [Zope \textit{et al.} J. Chem. Phys., 2019,{\bf 151}, 214108] when applied to the simplest local density approximation provides significant improvement over standard Perdew-Zunger SIC (PZSIC) for both equilibrium properties such as total or atomization energies as well as properties involving stretched bond such as barrier heights. The method uses an iso-orbital indicator to identify the single-electron regions. To demonstrate the LSIC method, Zope \textit{et al.} used the ratio $z_σ$ of von Weizsäcker $τ_σ^W$ and total kinetic energy densities $τ_σ$, ($z_σ= τ_σ^W/τ_σ$) as a scaling factor to scale the self-interaction correction. The present work further explores the LSIC method using a simpler scaling factor as a ratio of orbital and spin densities in place of the ratio of kinetic energy densities. We compute a wide array of both, equilibrium and non-equilibrium properties using the LSIC and orbital scaling methods using this simple scaling factor and compare them with previously reported results. Our study shows that the present results with simple scaling factor are comparable to those obtained by LSIC($z_σ$) for most properties but have slightly larger errors. We furthermore study the binding energies of small water clusters using both the scaling factors. Our results show that LSIC with $z_σ$ has limitation in predicting the binding energies of weakly bonded system due to the inability of $z_σ$ to distinguish weakly bonded region from slowly varying density region. LSIC when used with density ratio as a scaling factor, on the other hand, provides good description of water cluster binding energies, thus highlighting the appropriate choice of iso-orbital indicator.