Abstract:When the spins of molecular orbitals are allowed to be aligned with different directions in space rather than being aligned collinearly, the resulting noncollinear spin orbitals add extra flexibility to variational optimization of the orbitals, and solutions obtained with collinear spin orbitals may be unstable with respect to becoming noncollinear in the expanded variational space. The goal of the present work is to explore whether and in what way the molecular orbitals of the Kohn-Sham density functional theory become noncollinear when fully optimized for multi-reference molecules, transition states, and reaction paths. (We note that a noncollinear determinant has intermediate flexibility between a collinear determinant and a linear combination of many collinear determinants with completely independent coefficients. However, the Kohn-Sham method is defined to involve the variational optimization of a single determinant, and a noncollinear determinant represents the limit of complete optimization in the Kohn-Sham scheme.) We compare the results obtained with the noncollinear Kohn-Sham (NKS) scheme to those obtained with the widely used unrestricted Kohn-Sham (UKS) scheme for two types of multi-reference systems. For the dissociation of the MnO and NiO transition metal oxides, we find UKS fails to dissociate to the ground states of neutral atoms, while NKS dissociates to the correct limit and predicts potential energy curves that vary smoothly at intermediate bond lengths. This is due to the instability of UKS solutions at large bond distances. For barrier heights of O3, BeH2, and H4, NKS is shown to stabilize the multi-reference transition states by expanding the variational space. Although the errors vary because they are closely coupled with the capability of the employed exchange-correlation functionals in treating the multi-configurational states, these findings demonstrate that results with collinear spin orbitals should be further scrutinized, and future development of exchange-correlation functionals for multi-reference systems should incorporate the flexibilities of NKS.

Approach for noncollinear GGA kernels in closed-shell systems

Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation.

On the challenge to improve the density response with unusual gradient approximations

Noncollinear and Spin-Flip TDDFT in Multicollinear Approach

Two-dimensional limit of exchange-correlation energy functional approximations in density functional theory

A discontinuous functional for linear response time-dependent density functional theory: the exact-exchange kernel and approximate forms

Kernels for non interacting fermions via a Green’s function approach with applications to step potentials

Effect of orbital-overlap dependence in density functionals

Nonlocal Kinetic Energy Density Functionals for Isolated Systems Obtained Via Local Density Approximation Kernels

Ultranonlocality and accurate band gaps from a meta-generalized gradient approximation

Generalized Gradient Approximation Made Simple.

On the Efficient Evaluation of the Exchange Correlation Potential on Graphics Processing Unit Clusters

Gaussian expansion of Yukawa non-local kinetic energy functionals: application to metal clusters

Noncollinear Spin States for Density Functional Calculations of Open-Shell and Multi-Configurational Systems: Dissociation of MnO and NiO and Barrier Heights of O3, BeH2, and H4

Real-Time TDDFT Using Noncollinear Functionals

Non-equilibrium Green's function approach to inhomogeneous quantum many-body systems using the Generalized Kadanoff Baym Ansatz

Beyond the RPA and GW methods with adiabatic xc-kernels for accurate ground state and quasiparticle energies

Accurate Interatomic Force Fields via Machine Learning with Covariant Kernels

Nonlocal kinetic energy functionals by functional integration

A General and Transferable Local Hybrid Functional for Electronic Structure Theory and Many-Fermion Approaches

A simple generalized gradient approximation for the noninteracting kinetic energy density functional