Natural ionization orbitals for interpreting electron detachment processes

Visualizing and characterizing excited states from time-dependent density functional theory

Time-dependent density functional theory (TD-DFT) is the most widely-used electronic structure method for excited states, due to a favorable combination of low cost and semi-quantitative accuracy in many contexts, even if there are well recognized limitations. This Perspective describes various ways in which excited states from TD-DFT calculations can be visualized and analyzed, both qualitatively and quantitatively. This includes not just orbitals and densities but also well-defined statistical measures of electron-hole separation and of Frenkel-type exciton delocalization. Emphasis is placed on mathematical connections between methods that have often been discussed separately. Particular attention is paid to charge-transfer diagnostics, which provide indicators of when TD-DFT may not be trustworthy due to its categorical failure to describe long-range electron transfer. Measures of exciton size and charge separation that are directly connected to the underlying transition density are recommended over more ad hoc metrics for quantifying charge-transfer character.

Comparison of Linear Response Theory, Projected Initial Maximum Overlap Method, and Molecular Dynamics-Based Vibronic Spectra: The Case of Methylene Blue

The simulation of optical spectra is essential to molecular characterization and, in many cases, critical for interpreting experimental spectra. The most common method for simulating vibronic absorption spectra relies on the geometry optimization and computation of normal modes for ground and excited electronic states. In this report, we show that the utilization of such a procedure within an adiabatic linear response (LR) theory framework may lead to state mixings and a breakdown of the Born-Oppenheimer approximation, resulting in a poor description of absorption spectra. In contrast, computing excited states via a self-consistent field method in conjunction with a maximum overlap model produces states that are not subject to such mixings. We show that this latter method produces vibronic spectra much more aligned with vertical gradient and molecular dynamics (MD) trajectory-based approaches. For the methylene blue chromophore, we compare vibronic absorption spectra computed with the following: an adiabatic Hessian approach with LR theory-optimized structures and normal modes, a vertical gradient procedure, the Hessian and normal modes of maximum overlap method-optimized structures, and excitation energy time-correlation functions generated from an MD trajectory. Because of mixing between the bright S₁ and dark S₂ surfaces near the S₁ minimum, computing the adiabatic Hessian with LR theory and time-dependent density functional theory with the B3LYP density functional predicts a large vibronic shoulder for the absorption spectrum that is not present for any of the other methods. Spectral densities are analyzed and we compare the behavior of the key normal mode that in LR theory strongly couples to the optical excitation while showing S₁/S₂ state mixings. Overall, our study provides a note of caution in computing vibronic spectra using the excited-state adiabatic Hessian of LR theory-optimized structures and also showcases three alternatives that are less sensitive to adiabatic state mixing effects.

Using projection operators with maximum overlap methods to simplify challenging self-consistent field optimization

Maximum overlap methods are effective tools for optimizing challenging ground- and excited-state wave functions using self-consistent field models such as Hartree-Fock and Kohn-Sham density functional theory. Nevertheless, such models have shown significant sensitivity to the user-defined initial guess of the target wave function. In this work, a projection operator framework is defined and used to provide a metric for non-aufbau orbital selection in maximum-overlap-methods. The resulting algorithms, termed the Projection-based Maximum Overlap Method (PMOM) and Projection-based Initial Maximum Overlap Method (PIMOM), are shown to perform exceptionally well when using simple user-defined target solutions based on occupied/virtual molecular orbital permutations. This work also presents a new metric that provides a simple and conceptually convenient measure of agreement between the desired target and the current or final SCF results during a calculation employing a maximum-overlap method.

Lanthanide Oxides: From Diatomics to High-Spin, Strongly Correlated Homo- and Heterometallic Clusters

Small lanthanide (Ln) oxide clusters present both experimental and theoretical challenges because of their partially filled, core-like 4f ⁿ orbitals, a feature that results in a plethora of close-lying and fundamentally similar electronic states. These clusters provide a bottom-up approach toward understanding the electronic structure of defective or doped bulk material but also can offer a challenge to the theorists to find a method robust enough to capture electronic structure patterns that emerge from within the 4f ⁿ (0 < n < 14) series. In this Feature Article, we explore the electronic structures of small lanthanide oxide clusters that deviate from bulk stoichiometry using anion photoelectron spectroscopy and supporting density functional theory calculations. We will describe the evolution of electronic structure with oxidation and how Ln_xO_y^- cluster reactivities can be correlated with specific Ln-local orbital occupancies. These strongly correlated systems offer additional insights into how interactions between electrons and electronically complex neutrals can lead to detachment transitions that lie outside of the sudden one-electron detachment approximation generally assumed in anion photoelectron spectroscopy. With a better understanding of how we can control nominally forbidden transitions to sample an array of spin states, we suggest that more in-depth studies on the magnetic states of these systems can be explored. Extending these studies to other Ln-based materials with hidden magnetic phases, along with sequentially ligated single molecule magnets, could advance current understanding of these systems.

Natural transition orbitals for complex two-component excited state calculations

While the natural transition orbital (NTO) method has allowed electronic excitations from time-dependent Hartree-Fock and density functional theory to be viewed in a traditional orbital picture, the extension to multicomponent molecular orbitals such as those used in relativistic two-component methods or generalized Hartree-Fock (GHF) or generalized Kohn-Sham (GKS) is less straightforward due to mixing of spin-components and the inherent inclusion of spin-flip transitions in time-dependent GHF/GKS. An extension of single-component NTOs to the two-component framework is presented, in addition to a brief discussion of the practical aspects of visualizing two-component complex orbitals. Unlike the single-component analog, the method explicitly describes the spin and frequently obtains solutions with several significant orbital pairs. The method is presented using calculations on a mercury atom and a CrO₂ Cl₂ complex.

Exceptionally Complex Electronic Structures of Lanthanide Oxides and Small Molecules

Lanthanide (Ln) oxide clusters and molecular systems provide a bottom-up look at the electronic structures of the bulk materials because of close parallels in the patterns of Ln 4f^N subshell occupancy between the molecular and bulk Ln₂O₃ size limits. At the same time, these clusters and molecules offer a challenge to the theory community to find appropriate and robust treatments for the 4f^N patterns across the Ln series. Anion photoelectron (PE) spectroscopy provides a powerful experimental tool for studying these systems, mapping the energies of the ground and low-lying excited states of the neutral relative to the initial anion state, providing spectroscopic patterns that reflect the Ln 4f^N occupancy. In this Account, we review our anion PE spectroscopic and computational studies on a range of small lanthanide molecules and cluster species. The PE spectra of LnO^- (Ln = Ce, Pr, Sm, Eu) diatomic molecules show spectroscopic signatures associated with detachment of an electron from what can be described as a diffuse Ln 6s-like orbital. While the spectra of all four diatomics share this common transition, the fine structure in the transition becomes more complex with increasing 4f occupancy. This effect reflects increased coupling between the electrons occupying the corelike 4f and diffuse 6s orbitals with increasing N. Understanding the PE spectra of these diatomics sets the stage for interpreting the spectra of polyatomic molecular and cluster species. In general, the results confirm that the partial 4f^N subshell occupancy is largely preserved between molecular and bulk oxides and borides. However, they also suggest that surfaces and edges of bulk materials may support a low-energy, diffuse Ln 6s band, in contrast to bulk interiors, in which the 6s band is destabilized relative to the 5d band. We also identify cases in which the molecular Ln centers have 4f^N+1 occupancy rather than bulklike 4f^N, which results in weaker Ln-O bonding. Specifically, Sm centers in mixed Ce-Sm oxides or in Sm_xO_y^- (y ≤ x) clusters have this higher 4f^N+1 occupancy. The PE spectra of these particular species exhibit a striking increase in the relative intensities of excited-state transitions with decreasing photon energy (resulting in lower photoelectron kinetic energy). This is opposite of what is expected on the basis of the threshold laws that govern photodetachment. We relate this phenomenon to strong electron-neutral interactions unique to these complex electronic structures. The time scale of the interaction, which shakes up the electronic configuration of the neutral, increases with decreasing electron momentum. From a computational standpoint, we point out that special care must be taken when considering Ln cluster and molecular systems toward the center of the Ln series (e.g., Sm, Eu), where treatment of electrons explicitly or using an effective core potential can yield conflicting results on competing subshell occupancies. However, despite the complex electronic structures associated with partially filled 4f^N subshells, we demonstrate that inexpensive and tractable calculations yield useful qualitative insight into the general electronic structural features.