Fractional occupation numbers and self-interaction correction-scaling methods with the Fermi-Löwdin orbital self-interaction correction approach.

We present a new assessment of the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) approach with an emphasis on its performance for predicting energies as a function of fractional occupation numbers (FONs) for various multielectron systems. Our approach is implemented in the massively parallelized NWChem quantum chemistry software package and has been benchmarked on the prediction of total energies, atomization energies, and ionization potentials of small molecules and relatively large aromatic systems. Within our study, we also derive an alternate expression for the FLO-SIC energy gradient expressed in terms of gradients of the Fermi-orbital eigenvalues and revisit how the FLO-SIC methodology can be seen as a constrained unitary transformation of the canonical Kohn-Sham orbitals. Finally, we conclude with calculations of energies as a function of FONs using various SIC-scaling methods to test the limits of the FLO-SIC formalism on a variety of multielectron systems. We find that these relatively simple scaling methods do improve the prediction of total energies of atomic systems as well as enhance the accuracy of energies as a function of FONs for other multielectron chemical species.

Time-dependent density functional methods for Raman spectra in open-shell systems.

We present an implementation of a time-dependent density functional theory (TD-DFT) linear response module in NWChem for unrestricted DFT calculations and apply it to the calculation of resonant Raman spectra in open-shell molecular systems using the short-time approximation. The new source code was validated and applied to simulate Raman spectra on several doublet organic radicals (e.g., benzyl, benzosemiquinone, TMPD, trans-stilbene anion and cation, and methyl viologen) and the metal complex copper phthalocyanine. We also introduce a divide-and-conquer approach for the evaluation of polarizabilities in relatively large systems (e.g., copper phthalocyanine). The implemented tool gives comparisons with experiment that are similar to what is commonly found for closed-shell systems, with good agreement for most features except for small frequency shifts, and occasionally large deviations for some modes that depend on the molecular system studied, experimental conditions not being accounted in the modeling such as solvation effects and extra solvent-based peaks, and approximations in the underlying theory. The approximations used in the quantum chemical modeling include (i) choice of exchange-correlation functional and basis set; (ii) harmonic approximation used in the frequency analysis to determine vibrational normal modes; and (iii) short-time approximation (omission of nuclear motion effects) used in calculating resonant Raman spectra.

Additional Insights between Fermi-Löwdin Orbital SIC and the Localization Equation Constraints in SIC-DFT.

This letter highlights additional mathematical relationships between the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) formalism and the localization equation constraints in SIC-DFT. We demonstrate this relationship analytically by highlighting symmetries in the mathematical expression for the gradient of EPZ-SIC, which has not been previously shown in the scientific literature. To complement our analytical derivation, we also present additional numerical tests that allow us to investigate a possible accelerated-convergence technique that could be used when solving the iterative FLO-SIC equations. Taken together, our results highlight the importance of satisfying the localization equation constraints for obtaining accurate DFT energies, which we demonstrate are nearly satisfied in the FLO-SIC formalism.

Halogen Bonding Interactions: Revised Benchmarks and a New Assessment of Exchange vs Dispersion.

We present a new analysis of exchange and dispersion effects for calculating halogen-bonding interactions in a wide variety of complex dimers (69 total) within the XB18 and XB51 benchmark sets. Contrary to previous work on these systems, we find that dispersion plays a more significant role than exact exchange in accurately calculating halogen-bonding interaction energies, which are further confirmed by extensive SAPT analyses. In particular, we find that even if the amount of exact exchange is nonempirically tuned to satisfy known DFT constraints, we still observe an overall improvement in predicting dissociation energies when dispersion corrections are applied, in stark contrast to previous studies ( Kozuch, S.; Martin, J. M. L. J. Chem. Theory Comput. 2013 , 9 , 1918 - 1931 ). In addition to these new analyses, we correct several (14) inconsistencies in the XB51 set, which is widely used in the scientific literature for developing and benchmarking various DFT methods. Together, these new analyses and revised benchmarks emphasize the importance of dispersion and provide corrected reference values that are essential for developing/parametrizing new DFT functionals, specifically for complex halogen-bonding interactions.

Correlating Li+-Solvation Structure and its Electrochemical Reaction Kinetics with Sulfur in Subnano Confinement.

Combining theoretical and experimental approaches, we investigate the solvation properties of Li+ ions in a series of ether solvents (dimethoxyethane, diglyme, triglyme, tetraglyme, and 15-crown-5) and their subsequent effects on the solid-state lithium-sulfur reactions in subnano confinement. The ab initio and classical molecular dynamics (MD) simulations predict Li+ ion solvation structures within ether solvents in excellent agreement with experimental evidence from electrospray ionization-mass spectroscopy. An excellent correlation is also established between the Li+-solvation binding energies from the ab initio MD simulations and the lithiation overpotentials obtained from galvanostatic intermittent titration techniques (GITT). These findings convincingly indicate that a stronger solvation binding energy imposes a higher lithiation overpotential of sulfur in subnano confinement. The mechanistic understanding achieved at the electronic and atomistic level of how Li+-solvation dictates its electrochemical reactions with sulfur in subnano confinement provides invaluable guidance in designing future electrolytes and electrodes for Li-sulfur chemistry.