The combination of multipartitioning of the Hamiltonian with canonical Van Vleck perturbation theory leads to a Hermitian variant of quasidegenerate N-electron valence perturbation theory

Consistent Second-Order Treatment of Spin-Orbit Coupling and Dynamic Correlation in Quasidegenerate N-Electron Valence Perturbation Theory

We present a formulation and implementation of second-order quasidegenerate N-electron valence perturbation theory (QDNEVPT2) that provides a balanced and accurate description of spin-orbit coupling and dynamic correlation effects in multiconfigurational electronic states. In our approach, the energies and wave functions of electronic states are computed by treating electron repulsion and spin-orbit coupling operators as equal perturbations to the nonrelativistic complete active-space wave functions, and their contributions are incorporated fully up to the second order. The spin-orbit effects are described using the Breit-Pauli (BP) or exact two-component Douglas-Kroll-Hess (DKH) Hamiltonians within spin-orbit mean-field approximation. The resulting second-order methods (BP2- and DKH2-QDNEVPT2) are capable of treating spin-orbit coupling effects in nearly degenerate electronic states by diagonalizing an effective Hamiltonian expanded in a compact non-relativistic basis. For a variety of atoms and small molecules across the entire periodic table, we demonstrate that DKH2-QDNEVPT2 is competitive in accuracy with variational two-component relativistic theories. BP2-QDNEVPT2 shows high accuracy for the second- and third-period elements, but its performance deteriorates for heavier atoms and molecules. We also consider the first-order spin-orbit QDNEVPT2 approximations (BP1- and DKH1-QDNEVPT2), among which DKH1-QDNEVPT2 is reliable but less accurate than DKH2-QDNEVPT2. Both DKH1- and DKH2-QDNEVPT2 hold promise as efficient and accurate electronic structure methods for treating electron correlation and spin-orbit coupling in a variety of applications.

Bistability between π-diradical open-shell and closed-shell states in indeno[1,2-a]fluorene

Indenofluorenes are non-benzenoid conjugated hydrocarbons that have received great interest owing to their unusual electronic structure and potential applications in nonlinear optics and photovoltaics. Here we report the generation of unsubstituted indeno[1,2-a]fluorene on various surfaces by the cleavage of two C-H bonds in 7,12-dihydroindeno[1,2-a]fluorene through voltage pulses applied by the tip of a combined scanning tunnelling microscope and atomic force microscope. On bilayer NaCl on Au(111), indeno[1,2-a]fluorene is in the neutral charge state, but it exhibits charge bistability between neutral and anionic states on the lower-workfunction surfaces of bilayer NaCl on Ag(111) and Cu(111). In the neutral state, indeno[1,2-a]fluorene exhibits one of two ground states: an open-shell π-diradical state, predicted to be a triplet by density functional and multireference many-body perturbation theory calculations, or a closed-shell state with a para-quinodimethane moiety in the as-indacene core. We observe switching between open- and closed-shell states of a single molecule by changing its adsorption site on NaCl.

Unravelling the 6sp ← 6s absorption spectra of Bi(III) complexes

We report a spectroscopic and computational study that investigates the absorption spectra of Bi(III) complexes, which often show an absorption band in the UV region (∼270-350 nm) due to 6sp ← 6s transitions. We investigated the spectra of three simple complexes, [BiCl5]2-, [BiCl6]3- and [Bi(DMSO)8]3+, which show absorption maxima at 334, 326 and 279 nm due to 3P1 ← 1S0 transitions. Theoretical calculations based on quasi-degenerate N-electron valence perturbation theory to second order (QD-NEVPT2) provide an accurate description of the absorption spectra when employing CAS(2,9) wave functions. We next investigated the absorption spectra of the [Bi(NOTA)] complex (H3NOTA = 1,4,7-triazacyclononane-1,4,7-triacetic acid), which forms ternary complexes [Bi(NOTA)X]- (X = Cl, Br or I) in the presence of excess halide in aqueous solutions. Halide binding has an important impact on the position of the 3P1 ← 1S0 transition, which shifts progressively to longer wavelengths from 282 nm ([Bi(NOTA)]) to 298 nm (X = Cl), 305 nm (X = Br) and 325 nm (X = I). Subsequent QD-NEVPT2 calculations indicate that this effect is related to the progressive stabilization of the spin-orbit free states associated with the 6s16p1 configuration on increasing the covalent character of the metal-ligand(s) bonds, rather than with significant differences in spin-orbit coupling (SOC). These studies provide valuable insight into the coordination chemistry of Bi(III), an ion with increasing interest in targeted alpha therapy due to the possible application of bismuth isotopes bismuth-212 (212Bi, t1/2 = 60.6 min) and bismuth-213 (213Bi, t1/2 = 45.6 min).

A trivalent 4f complex with two bis-silylamide ligands displaying slow magnetic relaxation

The best-performing single-molecule magnets (SMMs) have historically relied on pseudoaxial ligands delocalized across several coordinated atoms. This coordination environment has been found to elicit strong magnetic anisotropy, but lanthanide-based SMMs with low coordination numbers have remained synthetically elusive species. Here we report a cationic 4f complex bearing only two bis-silylamide ligands, Yb(III)[{N(SiMePh₂)₂}₂][Al{OC(CF₃)₃}₄], which exhibits slow relaxation of its magnetization. The combination of the bulky silylamide ligands and weakly coordinating [Al{OC(CF₃)₃}₄]^- anion provides a sterically hindered environment that suitably stabilizes the pseudotrigonal geometry necessary to elicit strong ground-state magnetic anisotropy. The resolution of the m_J states by luminescence spectroscopy is supported by ab initio calculations, which show a large ground-state splitting of approximately 1,850 cm^-1. These results provide a facile route to access a bis-silylamido Yb(III) complex, and further underline the desirability of axially coordinated ligands with well-localized charges for high-performing SMMs.

Elucidating the Mechanism of Excited-State Bond Homolysis in Nickel-Bipyridine Photoredox Catalysts

Ni 2,2'-bipyridine (bpy) complexes are commonly employed photoredox catalysts of bond-forming reactions in organic chemistry. However, the mechanisms by which they operate are still under investigation. One potential mode of catalysis is via entry into Ni(I)/Ni(III) cycles, which can be made possible by light-induced, excited-state Ni(II)-C bond homolysis. Here, we report experimental and computational analyses of a library of Ni(II)-bpy aryl halide complexes, Ni(^Rbpy)(^R'Ph)Cl (R = MeO, t-Bu, H, MeOOC; R' = CH₃, H, OMe, F, CF₃), to illuminate the mechanism of excited-state bond homolysis. At given excitation wavelengths, photochemical homolysis rate constants span 2 orders of magnitude across these structures and correlate linearly with Hammett parameters of both bpy and aryl ligands, reflecting structural control over key metal-to-ligand charge-transfer (MLCT) and ligand-to-metal charge-transfer (LMCT) excited-state potential energy surfaces (PESs). Temperature- and wavelength-dependent investigations reveal moderate excited-state barriers (ΔH^‡ ∼ 4 kcal mol^-1) and a minimum energy excitation threshold (∼55 kcal mol^-1, 525 nm), respectively. Correlations to electronic structure calculations further support a mechanism in which repulsive triplet excited-state PESs featuring a critical aryl-to-Ni LMCT lead to bond rupture. Structural control over excited-state PESs provides a rational approach to utilize photonic energy and leverage excited-state bond homolysis processes in synthetic chemistry.

First-Principles Study of Optical Absorption Energies, Ligand Field and Spin-Hamiltonian Parameters of Cr3+ Ions in Emeralds

Herein, we study the electronic structure, energies, and vibronic structure of optical d-d transitions of Cr³⁺ ions doped in beryl (Be₃Si₆Al₂O₁₈:Cr³⁺, emerald). A computational protocol is developed that combines periodic density functional theory (for modeling of the bulk crystalline lattice of emerald) and the multireference configuration interaction complete active space self-consistent field method supplemented with n-electron valence second-order perturbation theory (for the calculation of the energy levels, wave functions, and spin-Hamiltonian and ligand-field parameters of the trigonal Cr³⁺ centers in the [CrO₆]^9- clusters embedded in an extended point charge field). Ligand-field parameters were extracted from mapping the effective ligand-field Hamiltonian onto the full many-particle Hamiltonian from one side and from a direct fit to energies of computed d-d transitions on the other side. These have been analyzed using ab initio ligand-field theory. The quality of the theoretical predictions is critically assessed through a detailed comparison with the available experimental data.

Improvement of Ab Initio Ligand Field Theory by Means of Multistate Perturbation Theory

Over the last few years, ab initio ligand field theory (AILFT) has evolved into an important tool for the extraction of ligand field models from ab initio calculations. The inclusion of dynamic correlation on top of complete active space self-consistent field (CASSCF) reference functions, which is important for accurate results, was so far realized at the level of second-order N-electron valence state perturbation theory (NEVPT2). In this work, we introduce two alternative methods for the inclusion of dynamic correlation into AILFT calculations, the second-order dynamic correlation dressed complete active space method (DCD-CAS(2)) and the Hermitian quasi-degenerate NEVPT2 (HQD-NEVPT2). These methods belong to the class of multistate perturbation theory approaches, which allow for the mixing of CASSCF states under the effect of dynamic correlation (state-mixing). The two new versions of AILFT were tested for a diverse set of transition-metal complexes. It was found that the multistate methods have, compared to NEVPT2, an AILFT fit with smaller root mean square deviations (rmsds) between ab initio and AILFT energies. A comparison of AILFT excitation energies with the experiment shows that for some systems, the agreement gets better at the multistate level because of the smaller rmsds. However, for some systems, the agreement gets worse, which could be attributed to a cancellation of errors at the NEVPT2 level that is partly removed at the multistate level. An investigation of trends in the extracted ligand field parameters shows that at the multistate level, the ligand field splitting Δ gets larger, whereas the Racah parameters B and C get smaller and larger, respectively. An investigation of the reasons for the observed improvement for octahedral Cr^III halide complexes shows that the possibility of state-mixing relaxes constraints that are present at the NEVPT2 level and that keep Δ and B from following their individual preferences.