Calculation of exchange couplings in the electronically excited state of molecular three-spin systems

Photogenerated molecular three-spin systems, composed of a chromophore and a covalently bound stable radical, are promising candidates for applications in the field of molecular spintronics. Through excitation with light, an excited doublet state and a quartet state are generated, whereby their energy difference depends on the exchange interaction JTR between the chromophore triplet state (T) and the stable radical (R). In order to establish design rules for new materials to be used in molecular spintronics devices, it is of great importance to gain knowledge on the magnitude of JTR as well as the factors influencing JTR on a molecular level. Here, we present a robust and reliable computational method to determine excited state exchange couplings in three-electron-three-centre systems based on a CASSCF/QD-NEVPT2 approach. The methodology is benchmarked and then applied to a series of molecules composed of a perylene chromophore covalently linked to various stable radicals. We calculate the phenomenological exchange interaction JTR between chromophore and radical, which can be compared directly to the experiment, but also illustrate how the individual exchange interactions Jij can be extracted using an effective Hamiltonian that corresponds to the Heisenberg–Dirac–Van-Vleck Hamiltonian. The latter procedure enables a more detailed analysis of the contributions to the exchange interaction JTR and yields additional insight that will be invaluable for future design optimisation

An excited state coupled-cluster study on indigo dyes

In the present study, the domain-based pair natural orbital implementation of the similarity-transformed equation of motion method is employed to reproduce the vibrationally resolved absorption spectra of indigo dyes. After an initial investigation of multireference, basis set and implicit solvent effects, our calculated 0–0 transition energies are compared to a benchmark set of experimental absorption band maxima. It is established that the agreement between our method and experimental results is well below the desired 0.1 eV threshold in virtually all cases and that the shift in excitation energies upon chemical substitution is also well reproduced. Finally, the entire spectra of some of the main components of the Tyrian purple dye mixture are reproduced and it is found that our computed spectra match the experimental ones without an empirical shift

A perturbative approach to multireference equation-of-motion coupled cluster

We introduce a variant of the multireference equation-of-motion coupled-cluster (MR-EOMCC) method where the amplitudes used for the similarity transformations are estimated from perturbation theory. Consequently, the new variant retains the many-body formalism, a reliance on at most two-body densities, and the state-universal character. As a non-iterative variant, computational costs are reduced, and no convergence difficulties with near-singular amplitudes can arise. Its performance was evaluated on several test sets covering transition metal atoms, small diatomics, and organic molecules against (near-)full CI quality reference data. We further highlight its efficacy on the weakly avoided crossing of LiF and place MR-EOMCC and the new variant into context with linear response theory. The accuracy of the variant was found to be at least on par with expectations for multireference perturbation theories, judging by the NEVPT2 method. The variant can be especially useful in multistate situations where the high accuracy of the iterative MR-EOMCC method is not required

Outer-Sphere Contributions to the Electronic Structure of Type Zero Copper Proteins

Bioinorganic canon states that active-site thiolate coordination promotes rapid electron transfer (ET) to and from type 1 copper proteins. In recent work, we have found that copper ET sites in proteins also can be constructed without thiolate ligation (called “type zero” sites). Here we report multifrequency electron paramagnetic resonance (EPR), magnetic circular dichroism (MCD), and nuclear magnetic resonance (NMR) spectroscopic data together with density functional theory (DFT) and spectroscopy-oriented configuration interaction (SORCI) calculations for type zero Pseudomonas aeruginosa azurin variants. Wild-type (type 1) and type zero copper centers experience virtually identical ligand fields. Moreover, O-donor covalency is enhanced in type zero centers relative that in the C112D (type 2) protein. At the same time, N-donor covalency is reduced in a similar fashion to type 1 centers. QM/MM and SORCI calculations show that the electronic structures of type zero and type 2 are intimately linked to the orientation and coordination mode of the carboxylate ligand, which in turn is influenced by outer-sphere hydrogen bonding

Molecular structure refinement by direct fitting of atomic coordinates to experimental ESR spectra

An attempt is made to bypass spectral analysis and fit internal coordinates of radicals directly to experimental liquid- and solid-state electron spin resonance (ESR) spectra. We take advantage of the recently introduced large-scale spin dynamics simulation algorithms and of the fact that the accuracy of quantum mechanical calculations of ESR parameters has improved to the point of quantitative correctness. Partial solutions are offered to the local minimum problem in spectral fitting and to the problem of spin interaction parameters (hyperfine couplings, chemical shifts, etc.) being very sensitive to distortions in molecular geometry.Comment: Submitted for publicatio

Ionization Energies and Redox Potentials of Hydrated Transition Metal Ions: Evaluation of Domain-Based Local Pair Natural Orbital Coupled Cluster Approaches

Hydrated transition metal ions are prototypical systems that can be used to model properties of transition metals in complex chemical environments. These seemingly simple systems present challenges for computational chemistry and are thus crucial in evaluations of quantum chemical methods for spin-state and redox energetics. In this work, we explore the applicability of the domain-based pair natural orbital implementation of coupled cluster (DLPNO-CC) theory to the calculation of ionization energies and redox potentials for hydrated ions of all first transition row (3d) metals in the 2+/3+ oxidation states, in connection with various solvation approaches. In terms of model definition, we investigate the construction of a minimally explicitly hydrated quantum cluster with a first and second hydration layer. We report on the convergence with respect to the coupled cluster expansion and the PNO space, as well as on the role of perturbative triple excitations. A recent implementation of the conductor-like polarizable continuum model (CPCM) for the DLPNO-CC approach is employed to determine self-consistent redox potentials at the coupled cluster level. Our results establish conditions for the convergence of DLPNO-CCSD(T) energetics and stress the absolute necessity to explicitly consider the second solvation sphere even when CPCM is used. The achievable accuracy for redox potentials of a practical DLPNO-based approach is, on average, 0.13 V. Furthermore, multilayer approaches that combine a higher-level DLPNO-CCSD(T) description of the first solvation sphere with a lower-level description of the second solvation layer are investigated. The present work establishes optimal and transferable methodological choices for employing DLPNO-based coupled cluster theory, the associated CPCM implementation, and cost-efficient multilayer derivatives of the approach for open-shell transition metal systems in complex environments

First-Principles Study of Optical Absorption Energies, Ligand Field and Spin-Hamiltonian Parameters of Cr³⁺ Ions in Emeralds

Herein, we study the electronic structure, energies, and vibronic structure of optical d-d transitions of Cr3+ ions doped in beryl (Be3Si6Al2O18:Cr3+, 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 Cr3+ centers in the [CrO6]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

Converged Structural and Spectroscopic Properties for Refined QM/MM Models of Azurin

Blue copper proteins continue to challenge experiment and theory with their electronic structure and spectroscopic properties that respond sensitively to the coordination environment of the copper ion. In this work, we report state-of-the art electronic structure studies for geometric and spectroscopic properties of the archetypal “Type I” copper protein azurin in its Cu(II) state. A hybrid quantum mechanics/molecular mechanics (QM/MM) approach is used, employing both density functional theory (DFT) and coupled cluster with singles, doubles, and perturbative triples (CCSD(T)) methods for the QM region, the latter method making use of the domain-based local pair natural orbital (DLPNO) approach. Models of increasing QM size are employed to investigate the convergence of critical geometric parameters. It is shown that convergence is slow and that a large QM region is critical for reproducing the short experimental Cu–SCys112 distance. The study of structural convergence is followed by investigation of spectroscopic parameters using both DFT and DLPNO-CC methods and comparing these to the experimental spectrum using simulations. The results allow us to examine for the first time the distribution of spin densities and hyperfine coupling constants at the coupled cluster level, leading us to revisit the experimental assignment of the 33S hyperfine splitting. The wavefunction-based approach to obtain spin-dependent properties of open-shell systems demonstrated here for the case of azurin is transferable and applicable to a large array of bioinorganic systems

Productive Alkyne Metathesis with “Canopy Catalysts” Mandates Pseudorotation

Molybdenum alkylidyne complexes of the “canopy catalyst” series define new standards in the field of alkyne metathesis. The tripodal ligand framework lowers the symmetry of the metallacyclobutadiene complex formed by [2 + 2] cycloaddition with the substrate and imposes constraints onto the productive [2 + 2] cycloreversion; pseudorotation corrects this handicap and makes catalytic turnover possible. A combined spectroscopic, crystallographic, and computational study provides insights into this unorthodox mechanism and uncovers the role that metallatetrahedrane complexes play in certain cases

Spectroscopic and Theoretical Study on Siloxy-Based Molybdenum and Tungsten Alkylidyne Catalysts for Alkyne Metathesis

A combined spectroscopic and theoretical study on triphenyl- and dimethyl-phenyl siloxy molybdenum and tungsten alkylidyne catalysts for alkyne metathesis is reported. Using NMR, X-ray, UV-VIS, Resonance Raman spectroscopy and DFT calculations, the influence of different ligand systems and metal centers on the geometric and electronic structure and thermochemistry of different intermediates is investigated, that is the starting alkylidyne and the derived metallacyclobutadiene and metallatetrahedrane. This includes a mechanistic and kinetic study on the formation and isomerization of metallacyclobutadienes and metallatetrahedranes. Upon changing from monodentate to tripodal siloxy ligands, higher steric strain is imposed, which modulates the relative energies of the different intermediates. Additionally, intramolecular dispersion interactions between bound substrate and ligand can be operative. Tungsten as the central metal leads to stronger M-C σ-bonds, which overstabilize the reaction intermediates and preclude effective turnover. Furthermore, kinetic modeling strongly suggests that metallatetrahedranes are off-cycle intermediates based on the high barriers for direct formation but low barriers for isomerization from metallacyclobutadienes. We infer from our findings that effective catalysis can only be achieved when factors that (over)stabilize intermediates, such as strong M-C bonds or large dispersion interactions, are prevented by appropriate catalyst design