Diagrammatic Multiplet-Sum Method (MSM) Density-Functional Theory (DFT): II. Completion of the Two-Orbital Two-Electron Model (TOTEM) with an Application to the Avoided Crossing in Lithium Hydride (LiH)

TL;DR

This paper extends the diagrammatic multiplet-sum method within density-functional theory to fully treat the two-orbital two-electron model, successfully describing the avoided crossing in lithium hydride and demonstrating consistent results across different functionals.

Contribution

It completes the development of the diagrammatic MSM-DFT for the full TOTEM, enabling accurate modeling of avoided crossings without parameter dependence.

Findings

LiH dissociates into neutral atoms with the method

All tested functionals produce similar results

Coupling element magnitude matches high-quality estimates

Abstract

The Ziegler-Rauk-Baerends multiplet sum method (MSM) assumes that density-functional theory (DFT) provides a good description of states dominated by a single determinant. It then uses symmetry to add static correlation to DFT. In our previous article (Article I) [J. Chem. Phys. 159, 244306 (2023)], we introduced diagrammatic MSM-DFT as a tool to aid in extending MSM-DFT to include the nondynamic correlation needed for making and breaking bonds even in the absence of symmetry. An attractive feature of this approach is that no functional-dependent parameters need to be introduced, though choices are needed in making correspondances between wave function theory (WFT) and MSM-DFT diagrams. The preliminary examples in Article I used the two-orbital two-electron model (TOTEM) less completely than could have been the case. Diagrammatic MSM-DFT is extended here to treat the full TOTEM and it is shown that the unsymmetric lithium hydride (LiH) molecule dissociates into neutral atoms when diagrammatic MSM-DFT techniques are used to introduce a proper description of the avoided crossing between ionic bonding and covalent bonding states.The method is tested for Hartree-Fock and for three functionals (LDA, PW91, and B3LYP). All the functionals yield similar results as should be expected for a properly-formulated parameter-free theory. Agreement with available estimates show that the magnitude of the coupling element introduced here is excellent. However more work will be needed to obtain quantitative agreement between our diagrammatic MSM-DFT ground-state potential energy curve and that found from high-quality ab initio calculations