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.

Enhanced Reactivity for Aromatic Bromination via Halogen Bonding with Lactic Acid Derivatives

We report a new method for regioselective aromatic bromination using lactic acid derivatives as halogen bond acceptors with N-bromosuccinimide (NBS). Several structural analogues of lactic acid affect the efficiency of aromatic brominations, presumably via Lewis acid/base halogen-bonding interactions. Rate comparisons of aromatic brominations demonstrate the reactivity enhancement available via catalytic additives capable of halogen bonding. Computational results demonstrate that Lewis basic additives interact with NBS to increase the electropositive character of bromine prior to electrophilic transfer. An optimized procedure using catalytic mandelic acid under aqueous conditions at room temperature was developed to promote aromatic bromination on a variety of arene substrates with complete regioselectivity.

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.

New Photoelectron-Valence Electron Interactions Evident in the Photoelectron Spectrum of Gd2O. J Phys Chem A 2021;125:9892-9903. [PMID: 34730978 DOI: 10.1021/acs.jpca.1c07818]

Evidence of strong photoelectron-valence electron (PEVE) interactions has been observed in the anion photoelectron (PE) spectra of several lanthanide suboxide clusters, which are exceptionally complex from an electronic structure standpoint and are strongly correlated systems. The PE spectrum of Gd₂O^-, which should have relatively simple electronic structure because of its half-filled 4f subshell, exhibits numerous electronic transitions. The electron affinity determined from the spectrum is 0.26 eV. The intensities of transitions to excited states increase relative to the lower-energy states with lower photon energy, which is consistent with shakeup transitions driven by time-dependent electron-neutral interactions. A group of intense spectral features that lie between electron binding energies of 0.7 and 2.3 eV are assigned to transitions involving detachment of an electron from outer-valence σ_u and σ_g orbitals that have large Gd 6s contributions. The spectra show parallel transition manifolds in general, which is consistent with detachment from these orbitals. However, several distinct perpendicular transitions are observed adjacent to several of the vertical transitions. A possible explanation invoking interaction between the ejected electron and the high-spin neutral is proposed. Specifically, the angular momentum of electrons ejected from σ_u or σ_g orbitals, which is l = 1, can switch to l = 0, 2 with an associated change in the M_s of the remnant neutral, which is spin-orbit coupling between a free electron and the spin of a neutral.

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.

Versatile New Reagent for Nitrosation under Mild Conditions

Here we report a new chemical reagent for transnitrosation under mild experimental conditions. This new reagent is stable to air and moisture across a broad range of temperatures and is effective for transnitrosation in multiple solvents. Compared with traditional nitrosation methods, our reagent shows high functional group tolerance for substrates that are susceptible to oxidation or reversible transnitrosation. Several challenging nitroso compounds are accessed here for the first time, including ¹⁵N isotopologues. X-ray data confirm that two rotational isomers of the reagent are configurationally stable at room temperature, although only one isomer is effective for transnitrosation. Computational analysis describes the energetics of rotamer interconversion, including interesting geometry-dependent hybridization effects.

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.

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.

Dispersion-Controlled Regioselective Acid-Catalyzed Intramolecular Hydroindolation of cis-Methindolylstyrenes To Access Tetrahydrobenzo[ cd]indoles

Readily prepared cis-β-(α',α'-dimethyl)-4'-methindolylstyrenes undergo acid-catalyzed intramolecular hydroindolation to afford tetrahydrobenzo[ cd]indoles. Our experimental and computational investigations suggest that dispersive interactions between the indole and styrene preorganize substrates such that 6-membered ring formation is preferred, apparently via concerted protonation and C-C bond formation. When dispersion is attenuated (by a substituent or heteroatom), regioselectivity erodes and competing oligomerization predominates for cis substrates. Similarly, all trans-configured substrates that we evaluated failed to cyclize efficiently.

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.

Hydrogen-bonding interactions in peptide nucleic acid and deoxyribonucleic acid: a comparative study

Peptide nucleic acid (PNA) is a synthetic analogue of deoxyribonucleic acid (DNA) capable of tightly binding to itself and DNA with high specificity. Using hybrid density functional methods, hydrogen-bond (H-bond) strengths have been evaluated for isolated Watson-Crick base pairs, PNA base pairs, and charged as well as neutral DNA base pairs. Heterogeneous base pairs of PNA with charged and neutral DNA have also been investigated. The competing effects of short-range H-bonding and long-range Coulombic repulsions in charged DNA base pairs have been analyzed. Polarizable continuum models have been employed to evaluate solvation effects on the binding energies.

Accurate reaction paths using a Hessian based predictor-corrector integrator

Central to the theoretical description of a chemical reaction is the reaction pathway. The intrinsic reaction coordinate is defined as the steepest descent path in mass weighted Cartesian coordinates that connects the transition state to reactants and products. In this work, a new integrator for the steepest descent pathway is presented. This method is a Hessian based predictor-corrector algorithm that affords pathways comparable to our previous fourth order method at the cost of a second order approach. The proposed integrator is tested on an analytic surface, four moderately sized chemical reactions, and one larger organometallic system.

Structural and Electronic Behavior of Unprecedented Five-Coordinate Iron(III) and Gallium(III) Complexes with a New Phenol-Rich Electroactive Ligand

A new asymmetric pentadentate ligand was designed to impose low symmetry to trivalent ions. Five-coordinate Fe3+ and Ga3+ complexes were investigated by crystallographic, electrochemical, and electron paramagnetic resonance methods showing enhanced redox reversibility. Calculations were performed to account for the observed trends.

Structural, spectroscopic, and electrochemical behavior of trans-phenolato cobalt(iii) complexes of asymmetric NN′O ligands as archetypes for metallomesogens

In order to understand and predict structural, redox, magnetic, and optical properties of more complex and potentially mesogenic electroactive compounds such as [Co(III)(L(t-BuLC))2]ClO4 (1), five archetypical complexes of general formula [Co(III)(L(RA))2]ClO4, where R = H (2), tert-butyl (3), methoxy (4), nitro (5), and chloro (6), were obtained and studied by means of several spectrometric, spectroscopic, and electrochemical methods. The complexes 2, 4, and 6 were characterized by single-crystal X-ray diffraction, and show the metal center in an approximate D2h symmetry. Experimental results support the fact that the electron donating or withdrawing nature of the phenolate-appended substituents changes dramatically the redox and spectroscopic properties of these compounds. The 3d6 electronic configuration of the metal ion dominates the overall geometry adopted by these compounds with the phenolate rings occupying trans positions to one another. Formation of phenoxyl radicals has been observed for 1, 3, and 6, but irreversible ligand oxidation takes place upon bulk electrolysis. These data were compared to detailed B3LYP/6-31G (d)-level computational calculations and have been used to account for the results observed. A comparison between compound 1 and archetype 3, validates the approach of using archetypical models to study metal-containing soft materials.