Modeling the Photoelectron Spectra of MoNbO2– Accounting for Spin Contamination in Density Functional Theory

Mapping spin contamination-free potential energy surfaces using restricted open-shell methods with Grassmannians

The Lagrange-based Grassmann interpolation (G-Int) method has been extended for open-shell systems using restricted open-shell (RO) methods. The performance of this method was assessed in constructing potential energy surfaces (PESs) for vanadium(II) oxide, benzyl radical, and methanesulfenyl chloride radical cation. The density matrices generated by G-Int when used as initial guesses for self-consistent field (SCF) calculations, exhibit superior performance compared to other traditional SCF initial guess schemes, such as SADMO, GWH, and CORE. Additionally, the energy obtained from the G-Int scheme satisfies the variational principle and outperforms the direct energy-based Lagrange interpolation approach. In the case of methanesulfenyl chloride radical cation, a unique example with a flat PES at the end region along the H-C-S-Cl dihedral angle, the use of an equally-spaced grid sampling leads to significant oscillations near the end of the interval due to the effects of Runge's phenomenon. Introducing an unequally-spaced grid sampling based on a scaled Gauss-Chebyshev quadrature effectively mitigated the Runge's phenomenon, making it suitable for combining with G-Int in constructing PESs for general applications. Thus, G-Int provides an efficient and robust strategy for building spin contamination-free PESs with consistent accuracy.

Good Vibrations: Calculating Excited-State Frequencies Using Ground-State Self-Consistent Field Models

The use of Δ-self-consistent field (SCF) approaches for studying excited electronic states has received a renewed interest in recent years. In this work, the use of this scheme for calculating excited-state vibrational frequencies is examined. Results from Δ-SCF calculations for a set of representative molecules are compared with those obtained using configuration interaction with single substitutions (CIS) and time-dependent density functional theory (TD-DFT) methods. The use of an approximate spin purification model is also considered for cases where the excited-state SCF solution is spin-contaminated. The results of this work demonstrate that an SCF-based description of an excited-state potential energy surface can be an accurate and cost-effective alternative to CIS and TD-DFT methods.

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.

A Study of NbMo and NbMo- by Anion Photoelectron Spectroscopy

Photoelectron spectra of the niobium-molybdenum diatomic anion, obtained at 488 and 514 nm, display vibrationally resolved transitions from the ground state and one excited electronic state of the anion to the ground state and one excited electronic state of the neutral molecule. The electron affinity of NbMo is measured to be 1.130 ± 0.005 eV. Its ²Δ_3/2 spin-orbit component is observed to lie 870 ± 20 cm^-1 above its previously identified ²Δ_5/2 ground state. For ⁹³Nb⁹⁸Mo, vibrational energies measured for levels up to v = 4 for the ²Δ_5/2 and ²Δ_3/2 states give harmonic frequency and anharmonicity constant values of ω_e = 492 ± 12 cm^-1 and ω_ex_e = 8.0 ± 3.2 cm^-1, the former value corresponding to a force constant of 6.80 ± 0.35 mdyn/Å. These two vibrational parameters suggest a bond dissociation energy that is too low by at least a factor of 3, indicating that the ground state potential energy curve of NbMo deviates markedly from a Morse potential at higher energies. An excited electronic state of NbMo, assigned as a ²Σ⁺ state, is observed at 2900 ± 25 cm^-1 (T₀). Vibrational energies up to v = 8 in this excited state give values of ω_e = 544 ± 8 cm^-1 and ω_ex_e = 1.9 ± 1.2 cm^-1 for ⁹³Nb⁹⁸Mo. The former value corresponds to a high vibrational force constant of 8.30 ± 0.25 mdyn/Å. Both doublet states of the neutral molecule are accessed from the anion ground state, which is assigned as ¹Σ⁺. For the ⁹³Nb⁹⁸Mo^- anion, the fundamental vibrational frequency (ΔG_1/2) is 484 ± 15 cm^-1. Electron affinity data indicate that the bond dissociation energy of NbMo^- is 0.213 ± 0.005 eV greater than that of neutral NbMo, whose previously reported value then gives D₀ = 4.85 ± 0.27 eV for the anion. An excited state of the anion lying 3050 ± 25 cm^-1 (T₀) above its ground state is assigned as ³Δ, and the energies of its spin-orbit components above the ³Δ₃ lowest energy level are measured to be 450 ± 20 cm^-1 (³Δ₂) and 1100 ± 20 cm^-1 (³Δ₁). Their uneven spacing suggests that the energy of the ³Δ₂ level is lowered by interaction with a higher energy Ω = 2 anion state. The vibrational frequency (ΔG_1/2) for the ³Δ₁ and ³Δ₂ states is measured to be 433 ± 20 cm^-1. Bond length differences among the observed states are estimated from Franck-Condon fits to vibrational band intensity profiles. When combined with the previously reported NbMo bond length, these provide bond length estimates for the ground state of the anion (1.940 ± 0.025 Å) and for the observed excited states. These species offer extreme examples of multiple metal-metal bonding, with formal bond orders of 5¹/₂ for the ²Δ ground and ²Σ⁺ excited doublet states of NbMo, 6 for the singlet ground state of the anion, and 5 for its low-lying triplet state. The relationships among their bonding properties and those of related multiply bonded transition metal dimers are discussed.

Assessing the Calculation of Exchange Coupling Constants and Spin Crossover Gaps Using the Approximate Projection Model To Improve Density Functional Calculations

This work evaluates the quality of exchange coupling constant and spin crossover gap calculations using density functional theory corrected by the approximate projection model. Results show that improvements using the approximate projection model range from modest to significant. This study demonstrates that, at least for the class of systems examined here, spin projection generally improves the quality of density functional theory calculations of J-coupling constants and spin crossover gaps. Furthermore, it is shown that spin projection can be important for both geometry optimization and energy evaluations. The approximate projection model provides an affordable and practical approach for effectively correcting spin-contamination errors in such calculations.

On the linear geometry of lanthanide hydroxide (Ln-OH, Ln = La–Lu)

Density functional theory predicts that lanthanide hydroxides are linear, with the lanthanide-hydroxide bond being characterized as a covalent triple bond.