Multireference Ab Initio Calculations Research Articles

Multireference Ab Initio Calculations on Excited Electronic States of Carbazole-Based Organic Compounds for Thermally Activated Delayed Fluorescence.

The emerging technology of organic light-emitting diodes takes advantage of the thermally activated delayed fluorescence (TADF) mechanism for improved efficiency. Carbazole-based organic molecules are suitable for TADF emission because of charge transfer excitations between the electron-donor carbazole and an electron-acceptor unit. Computational design of new TADF molecules with the desired properties is challenging because charge-transfer excitations cannot be predicted accurately by time-dependent density functional theory. Four groups of carbazole-based donor-acceptor molecules have been studied using multireference ab initio methods to understand the nature of excited electronic states. The state-averaged complete active space self-consistent field (SA-CASSCF) and the N-electron valence state perturbation theory (NEVPT2) were used to calculate energies and oscillator strengths for multiple excited electronic states. The number of active electrons and orbitals and the number of excited states included in state-averaged CASSCF were selected such that the accuracy of ab initio predictions could be improved systematically. The procedure introduced here for the calculation of multiple excited electronic states of TADF candidates can be used to accelerate the computational search for efficient TADF materials.

Oxygen-Donor Metalloligands Induce Slow Magnetization Relaxation in Zero Field for a Cobalt(II) Complex with {CoO4} Motif.

Most 3d metal-based single-molecule magnets (SMMs) use N-ligands or ligands with even softer donors to impart a particular coordination geometry and increase the zero-field splitting parameter |D|, while complexes with hard O-donor ligands showing slow magnetization relaxation are rare. Here, we report that a diamagnetic NiII complex of a tetradentate ligand featuring two N-heterocyclic carbene and two alkoxide-O donors, [LO,ONi], can serve as a {O,O'}-chelating metalloligand to give a trinuclear complex [(LO,ONi)Co(LO,ONi)](OTf)2 (2) with an elongated tetrahedral {CoIIO4} core, D = -74.3 cm-1, and a spin reversal barrier Ueff = 86.9 cm-1 in the absence of an external dc field. The influence of diamagnetic NiII on the electronic structure of the {CoO4} unit in comparison to [Co(OPh)4]2- (A) has been probed with multireference ab initio calculations. These reveal a contrapolarizing effect of the NiII, which forms stronger metal-alkoxide bonds than the central CoII, inducing a change in ligand field splitting and a 5-fold increase in the magnetic anisotropy in 2 compared to A, with an easy magnetization axis along the Ni-Co-Ni vector. This demonstrates a strategy to enhance the SMM properties of 3d metal complexes with hard O-donors by modulating the ligand field character via the coordination of diamagnetic ions and the benefit of robust metalloligands in that regard.

Doubly charged dimers and trimers of heavy noble gases.

Many doubly charged heteronuclear dimers are metastable or even thermodynamically stable with respect to charge separation. Homonuclear dicationic dimers, however, are more difficult to form. He22+ was the first noble gas dimer predicted to be metastable and, decades later, observed. Ne22+ is the only other dicationic noble gas dimer that has been detected so far. Here, we present a novel approach to form fragile dicationic species, by post-ionization of singly charged ions that are embedded in helium nanodroplets (HNDs). Bare ions are then extracted by colliding the HNDs with helium gas. We detect homonuclear doubly charged dimers and trimers of krypton and xenon, but not argon. Our multi-reference ab initio calculations confirm the stability of Kr22+, Kr32+, Xe22+, Xe32+, and Ar22+, but put the stability of Ar32+ towards dissociation to Ar+ + Ar2+ into question.

Electronic Structure and Excited-State Dynamics of the NIR-II Emissive Molybdenum(III) Analogue to the Molecular Ruby.

Photoactive chromium(III) complexes saw a conceptual breakthrough with the discovery of the prototypical molecular ruby mer-[Cr(ddpd)2]3+ (ddpd = N,N'-dimethyl-N,N'-dipyridin-2-ylpyridine-2,6-diamine), which shows intense long-lived near-infrared (NIR) phosphorescence from metal-centered spin-flip states. In contrast to the numerous studies on chromium(III) photophysics, only 10 luminescent molybdenum(III) complexes have been reported so far. Here, we present the synthesis and characterization of mer-MoX3(ddpd) (1, X = Cl; 2, X = Br) and cisfac-[Mo(ddpd)2]3+ (cisfac-[3]3+), an isomeric heavy homologue of the prototypical molecular ruby. For cisfac-[3]3+, we found strong zero-field splitting using magnetic susceptibility measurements and electron paramagnetic resonance spectroscopy. Electronic spectra covering the spin-forbidden transitions show that the spin-flip states in mer-1, mer-2, and cisfac-[3]3+ are much lower in energy than those in comparable chromium(III) compounds. While all three complexes show weak spin-flip phosphorescence in NIR-II, the emission of cisfac-[3]3+ peaking at 1550 nm is particularly low in energy. Femtosecond transient absorption spectroscopy reveals a short excited-state lifetime of 1.4 ns, 6 orders of magnitude shorter than that of mer-[Cr(ddpd)2]3+. Using density functional theory and ab initio multireference calculations, we break down the reasons for this disparity and derive principles for the design of future stable photoactive molybdenum(III) complexes.

Long-range ab initio potential energy surface for the ground electronic state of the ozone molecule with the accurate dissociation asymptote

In this work an excellent agreement was achieved between the empirical (9230.08 ± 4.85 cm−1) and theoretical (9230.70 cm−1) values of the first dissociation limit (De) of O3. Ab initio multireference calculations were performed for the first time using the extended complete active space [CAS(24,15)] reference in the icMRCI(Q) method combined with the extrapolation to the complete basis set limit [CBS(56Z)] and with the inclusion of the spin–orbit coupling between 9 multiplet electronic states approaching De. Finally, the full-dimensional potential energy surface was developed in the scattering coordinates for the long-range region toward De at this level of theory.

Structural and Magnetization Dynamics of Borohydride-Bridged Rare-Earth Metallocenium Cations.

The structure and magnetic properties of the bimetallic borohydride-bridged dysprosocenium compound [{(η5-Cpttt)(η5-CpMe4t)Dy}2(μ:κ2:κ2-BH4)]+[B(C6F5)4]- ([3Dy][B(C6F5)4]) are reported along with the solution-phase dynamics of the isostructural yttrium and lutetium analogues (Cpttt is 1,2,4-tri(tert-butyl)cyclopentadienyl, CpMe4t is tetramethyl(tert-butyl)cyclopentadienyl). The synthesis of [3M][B(C6F5)4] was accomplished in the 2:1 stoichiometric reactions of [(η5-Cpttt)(η5-CpMe4t)Dy(BH4)] (2M) with [CPh3][B(C6F5)4], with the metallocenes 2M obtained from reactions of the half-sandwich complexes [(η5-Cpttt)M(BH4)2(THF)] (1M) (M = Y, Dy, Lu) with NaCpMe4t. Crystallographic studies show significant lengthening of the M···B distance on moving through the series 1M, 2M, and 3M, with essentially linear {M···B···M} bridges in 3M. Multinuclear NMR spectroscopy indicates restricted rotation of the Cpttt ligands in 3Y and 3Lu in solution. The single-molecule magnet (SMM) properties of [3M][B(C6F5)4] are characterized by Raman and Orbach processes, with an effective barrier of 533(18) cm-1 and relaxation via the second-excited Kramers doublet. Although quantum tunneling of the magnetization (QTM) was not observed for [3M][B(C6F5)4], it was, surprisingly, found in its magnetically dilute version, which has a very similar barrier of Ueff = 499(21) cm-1. Consistent with this observation, slightly wider openings of the magnetic hysteresis loop at 2 K are found for [3M][B(C6F5)4] but not for the diluted analogue. The dynamic magnetic properties of the dysprosium SMMs and the role of exchange interactions in 3Dy are interpreted with the aid of multireference ab initio calculations.

The Role of the Redox Non-Innocent Hydroxyl Ligand in the Activation of O2 Performed by [Ni(H)(OH)

The cationic complex [Ni(H)(OH)]+ was previously found to activate dioxygen and methane in gas phase under single collision conditions. These remarkable reactivities were thought to originate from a non-classical electronic structure, where the Ni-center adopts a Ni(II), instead of the classically expected Ni(III) oxidation state by formally accepting an electron from the hydroxo ligand, which formally becomes a hydroxyl radical in the process. Such radicaloid oxygen moieties are envisioned to easily react with otherwise inert substrates, mimicking familiar reactivities of free radicals. In this study, the reductive activation of dioxygen by [Ni(H)(OH)]+ to afford the hydroperoxo species was investigated using coupled cluster, multireference ab initio and density functional theory calculations. Orbital and wave function analyses indicate that O2 binding tranforms the aforementioned non-classical electronic structure to a classical Ni(III)-hydroxyl system, before O2 reduction takes place. Remarkably, we found no evidence for a direct involvement of the radicaloid hydroxyl in the reaction with O2 , as is often assumed. The function of the redox non-innocent character of the activator complex is to protect the reactive electronic structure until the complex engages O2 , upon which a dramatic electronic reorganization releases internal energy and drives the chemical reaction to completion.

Harnessing Colossal Magnetic Anisotropy in Sandwiched 3d2-Metallocenes.

Single-molecule magnets are gaining attention in recent years with the growing focus on achieving higher barriers of magnetization reversal. Metallocenes, owing to their unique sandwiched structure, assure themselves as plausible molecular systems for the development of novel single-molecule magnets (SMMs). Here in this work, we have explicitly investigated metallocenes of first-row transition elements, along with their one-electron-oxidized (cationic) and -reduced (anionic) analogues, for their magnetic anisotropies by adopting multireference ab initio calculations. Herein, we report a high magnetic anisotropy for 3d2 systems among all 3d-metallocenes.

Characterization of a strong covalent Th3+\u2013Th3+ bond inside an Ih(7)-C80 fullerene cage

The nature of the actinide-actinide bonds is of fundamental importance to understand the electronic structure of the 5f elements. It has attracted considerable theoretical attention, but little is known experimentally as the synthesis of these chemical bonds remains extremely challenging. Herein, we report a strong covalent Th-Th bond formed between two rarely accessible Th3+ ions, stabilized inside a fullerene cage nanocontainer as Th2@Ih(7)-C80. This compound is synthesized using the arc-discharge method and fully characterized using several techniques. The single-crystal X-Ray diffraction analysis determines that the two Th atoms are separated by 3.816 Å. Both experimental and quantum-chemical results show that the two Th atoms have formal charges of +3 and confirm the presence of a strong covalent Th-Th bond inside Ih(7)-C80. Moreover, density functional theory and ab initio multireference calculations suggest that the overlap between the 7s/6d hybrid thorium orbitals is so large that the bond still exists at Th-Th separations larger than 6 Å. This work demonstrates the authenticity of covalent actinide metal-metal bonds in a stable compound and deepens our fundamental understanding of f element metal bonds.

Multi-reference ab initio calculations of Hg spectral data and analysis of magic and zero-magic wavelengths.

We have identified magic wavelengths for 1S0 ↔ 3P1,2 (mJ = 0) transitions and zero-magic wavelengths for the 3P1,2 (mJ = 0) states of 200Hg atoms, analysed the robustness of the magic conditions with respect to wavelength and polarization imperfections. We show that the most experimentally feasible magic wavelength for the 1S0 ↔ 3P2 transition is 351.8 nm of π polarized light. Relevant transition wavelengths and transition strengths are calculated using the state-of-the-art Complete Active Space Self-Consistent-Field (CASSCF) method with a perturbative inclusion of spin-orbit coupling. The transition wavelengths are a posteriori corrected for the dynamical energy using the second-order perturbation theory.

Multireference Description of Nickel-Aryl Homolytic Bond Dissociation Processes in Photoredox Catalysis.

Multireference electronic structure calculations consistent with known experimental data have elucidated a novel mechanism for photo-triggered Ni(II)-C homolytic bond dissociation in Ni 2,2'-bipyridine (bpy) photoredox catalysts. Previously, a thermally assisted dissociation from the lowest energy triplet ligand field excited state was proposed and supported by density functional theory (DFT) calculations that reveal a barrier of ∼30 kcal mol-1. In contrast, multireference ab initio calculations suggest that this process is disfavored, with barrier heights of ∼70 kcal mol-1, and highlight important ligand noninnocent and multiconfigurational contributions to excited state relaxation and bond dissociation processes that are not captured with DFT. In the multireference description, photo-triggered Ni(II)-C homolytic bond dissociation occurs via initial population of a singlet Ni(II)-to-bpy metal-to-ligand charge transfer (1MLCT) excited state, followed by intersystem crossing and aryl-to-Ni(III) charge transfer, overall a formal two-electron transfer process driven by a single photon. This results in repulsive triplet excited states from which spontaneous homolytic bond dissociation can occur, effectively competing with relaxation to the lowest energy nondissociative triplet Ni(II) ligand field excited state. These findings guide important electronic structure considerations for the experimental and computational elucidation of the mechanisms of ground and excited state cross-coupling catalysis mediated by Ni heteroaromatic complexes.

Probing the Magnetic Anisotropy of Co(II) Complexes Featuring Redox-Active Ligands.

Coordination complexes that possess large magnetic anisotropy (otherwise known as zero-field splitting, ZFS) have possible applications in the field of magnetic materials, including single molecule magnets (SMMs). Previous studies have explored the role of coordination number and geometry in controlling the magnetic anisotropy and SMM behavior of high-spin (S = 3/2) Co(II) complexes. Building upon these efforts, the present work examines the impact of ligand oxidation state and structural distortions on the spin states and ZFS parameters of pentacoordinate Co(II) complexes. The five complexes included in this study (1-5) have the general formula, [Co(TpPh2)(LX,Y)]n+ (X = O, S; Y = N, O; n = 0 or 1), where TpPh2 is the scorpionate ligand hydrotris(3,5-diphenyl-pyrazolyl)borate(1-) and LX,Y are bidentate dioxolene-type ligands that can access multiple oxidation states. The specific LX,Y ligands used herein are 4,6-di-tert-butyl substituted o-aminophenolate and o-aminothiophenolate (1 and 2, respectively), o-iminosemiquinonate and o-semiquinonate radicals (3 and 4, respectively), and o-iminobenzoquinone (5). Each complex exhibits a distorted trigonal bipyramidal geometry, as revealed by single-crystal X-ray diffraction. Direct current (dc) magnetic susceptibility experiments confirmed that the complexes with closed-shell ligands (1, 2, and 5) possess S = 3/2 ground states with negative D-values (easy-axis anisotropy) of -41, -78, and -30 cm-1, respectively. For 3 and 4, antiferromagnetic coupling between the Co(II) center and o-(imino)semiquinonate radical ligand results in S = 1 ground states that likewise exhibit very large and negative anisotropy (-100 > D > -140 cm-1). Notably, ZFS was measured directly for each complex using far-infrared magnetic spectroscopy (FIRMS). In combination with high-frequency and -field electron paramagnetic resonance (HFEPR) studies, these techniques provided precise spin-Hamiltonian parameters for complexes 1, 2, and 5. Multireference ab initio calculations, using the CASSCF/NEVPT2 approach, indicate that the strongly negative anisotropies of these Co(II) complexes arise primarily from distortions in the equatorial plane due to constrictions imposed by the TpPh2 ligand. This effect is further amplified by cobalt(II)-radical exchange interactions in 3 and 4.

Phenanthroline chromophore as efficient antenna for Tb3+ green luminescence: A theoretical study

The photophysical processes, from excitation in the UV range, energy migration pathways to luminescence of Tb(phen)2(NO3)3 complex in the Vis region have been studied by theoretical approaches and spectroscopic measurements. The UV–Vis absorption, diffuse reflectance, excitation and luminescence spectra indicated that the ligand 1,10-phenanthroline (phen) was a suitable antenna chromophore for Tb3+ luminescence at excitation of wavelength range 350–352 nm, corresponding to the S1 state of phen. The energy level diagram and character of the singlet and triplet excited states of the ligand-chromophore and Tb(phen)2(NO3)3 in gas phase and polar solution were obtained by the DFT/TDDFT/ωB97xD/B1 method and multireference ab-initio calculations with perturbative spin-orbit coupling corrections. It was also established that the vibro-electronic contribution to the excitation energy is important for correct reproduction and interpretation of the experimental absorption and emission spectra of Tb(phen)2(NO3)3. The rate constants for competitive singlet-triplet intersystem crossing relaxation transitions, as well as radiative and non-radiative deactivation of the ligand-centered excited states were qualitatively estimated and the fastest electronic relaxation mechanism was revealed. The excitation channel S0→S1→T2→T1→5D4 was predicted as the most probable one for effective phen→Tb3+ energy transfer. The sensitization results in a high luminescent quantum yield of 13%, which suggests energy transfer rates, sufficiently higher than the rates of radiative and non-radiative depopulation of the T1 state.

Photochemical activation of carbon dioxide in Mg+(CO2)(H2O)0,1

We combine multi-reference ab initio calculations with UV–VIS action spectroscopy to study photochemical activation of CO2 on a singly charged magnesium ion, [MgCO2(H2O)0,1]+, as a model system for the metal/ligand interactions relevant in CO2 photochemistry. For the non-hydrated species, two separated Mg+ 3s–3p bands are observed within 5.0 eV. The low-energy band splits upon hydration with one water molecule. [Mg(CO2)]+ decomposes highly state-selectively, predominantly via multiphoton processes. Within the low-energy band, CO2 is exclusively lost within the excited state manifold. For the high-energy band, an additional pathway becomes accessible: the CO2 ligand is activated via a charge transfer, with photochemistry taking place on the CO2– moiety eventually leading to a loss of CO after absorption of a second photon. Upon hydration, already excitation into the first and second excited state leads to CO2 activation in the excited state minimum; however, CO2 predominantly evaporates upon fluorescence or absorption of another photon.

Axial Ligation in Ytterbium(III) DOTAM Complexes Rationalized with Multireference and Ligand-Field ab Initio Calculations.

The nature of the axial ligand coordinated to the Yb3+ ion in [Yb(DOTAM)]3+ has profound consequences in the magnetic anisotropy and optical properties of the complex, as evidenced by 1H NMR and UV-vis spectroscopies. The pseudocontact shifts of 1H nuclei and the 2F5/2 ← 2F7/2 absorption band were found to be very sensitive to the nature of the axial ligand (MeOH, H2O, MeOH, or F-). The energy levels of the 2F5/2 and 2F7/2 manifolds in [Yb(DOTAM)(X)]3+ (X = MeOH, H2O, or dimethyl sulfoxide (DMSO)) and [Yb(DOTAM)F]2+ complexes were assigned from the analysis of the optical spectra and ab initio calculations based on CASSCF wave functions that considered dynamic correlation through perturbation theory (NEVPT2) and spin-orbit coupling effects. The magnetic anisotropies obtained with ab initio calculations are in good agreement with the experimental values derived from 1H NMR spectral data, though for the [Yb(DOTAM)(H2O)]3+ and [Yb(DOTAM)F]2+ complexes, the explicit inclusion of a few second-sphere water molecules is required to improve the calculated data. Crystal-field calculations show that the observed pseudocontact shifts do not correlate well with the crystal-field parameter B20, as predicted by Bleaney's theory. The change in the sign of the magnetic anisotropy from prolate (X = MeOH, H2O, or DMSO) to oblate in [Yb(DOTAM)F]2+ is related to the relative energies of the 4fz3 orbital and the 4fx3/4fy3 pair, which are affected by the coordination ability of the axial ligand.

Ab Initio Analysis of Metal-Ligand Bonding in An(COT)2 with An=Th, U in Their Ground- and Core-Excited States.

Relativistic multireference ab initio wave function calculations with the restricted active space second-order perturbation theory (RASPT2) were performed on thorocene and uranocene to determine the actinide N4,5 -edge and carbon K-edge X-ray absorption near-edge structure (XANES) intensities and the metal-ligand orbital mixing in the ground state and core-excited states. Calculated spectral intensities show very good agreement with the experiments and therefore allow detailed and unambiguous assignment of the observed spectral features. φ-type covalent bonding or antibonding interactions are observed for thorocene in the core-excited states, though not in the ground state. This is because the molecular orbital of φ symmetry, which is the in-phase combination of the ligand Lφ and the Th 5fφ orbitals, can be populated with electrons in core-excited states, whereas it is essentially unoccupied in the ground state. For uranocene, the XANES spectra do not reveal much information beyond multiplet broadening, despite the presence of distinct peaks in the spectra. Every core-excited peak is best characterized by its own set of bond orbitals, as the excited state covalency is clearly different from the ground state covalency.

Accurate Potential Energy Surface for Quartet State HN2 and Interplay of N(4S) + NH(X̃3Σ–) versus H + N2(A3Σu+) Reactions

A global potential energy surface for the lowest quartet state of HN2 is reported for the first time from accurate multireference ab initio calculations extrapolated to the complete basis set limit using the double many-body expansion method. All its stationary points are characterized, with the lowest quartet of HN2 predicted to have a bent global minimum 36 kcal mol-1 below the N(4S) + NH(X̃3Σ-) asymptote, from which it is barrierlessly achievable. The entire set of calculated ab initio points has been fitted for energies up to 1000 kcal mol-1 above the global minimum with an RMSD of 0.89 kcal mol-1, a gap comprising all identified stationary points. Special care is taken in modeling the involved long-range forces and cusps caused by crossing seams. The novel PES prompts for the calculation of rate constants for several unexplored reactions that are relevant for combustion, plasma, and atmospheric chemistry.

Control of magnetic anisotropy by macrocyclic ligand distortion in a family of DyIII and ErIII single molecule magnets.

A family of hexaazamacrocyclic lanthanide complexes, [Ln(Ln)(NCS)3] (LnIII = Dy, Er; n = 1-3) has been synthesized and characterized by single-crystal X-ray diffraction, magnetic measurements and ab initio calculations. Macrocyclic ligands (Ln) differ in the lateral spacers, which are aliphatic chains with two and three carbons (for Ln, n = 1 and 2, respectively), and an aromatic ring for Ln = 3. Modification of the macrocycle spacer tunes planarity and rigidity of the equatorial coordination for both oblate (Dy) and prolate (Er) lanthanide ions. Ac-susceptibility studies showed that four of the six complexes are field induced single molecule magnets (SMMs). Trends in magnetic relaxation properties are rationalized with the aid of ab initio multireference calculations, highlighting the combined influence of macrocycle planarity, lanthanide electronic density distribution and intermolecular interactions for the achievement of slow demagnetization.

Masked Radicals: Iron Complexes of Trityl, Benzophenone, and Phenylacetylene.

We report the first Fe─CPh3 complex, and show that the long Fe─C bond can be disrupted by neutral π-acceptor ligands (benzophenone and phenylacetylene) to release the triphenylmethyl radical. The products are formally iron(I) complexes, but X-ray absorption spectroscopy coupled with density functional and multireference ab initio calculations indicates that the best description of all the complexes is iron(II). In the formally iron(I) complexes, this does not imply that the π-acceptor ligand has radical character, because the iron(II) description arises from doubly-occupied frontier molecular orbitals that are shared equitably by the iron and the π-acceptor ligand, and the unpaired electrons lie on the metal. Despite the lack of substantial radical character on the ligands, alkyne and ketone fragments can couple to form a high-spin iron(III) complex with a cyclized metalladihydrofuran core.

M−O Bonding Beyond the Oxo Wall: Spectroscopy and Reactivity of Cobalt(III)‐Oxyl and Cobalt(III)‐Oxo Complexes

Terminal oxo complexes of late transition metals are frequently proposed reactive intermediates. However, they are scarcely known beyond Group 8. Using mass spectrometry, we prepared and characterized two such complexes: [(N4Py)CoIII(O)]+ (1) and [(N4Py)CoIV(O)]2+ (2). Infrared photodissociation spectroscopy revealed that the Co−O bond in 1 is rather strong, in accordance with its lack of chemical reactivity. On the contrary, 2 has a very weak Co−O bond characterized by a stretching frequency of ≤659 cm−1. Accordingly, 2 can abstract hydrogen atoms from non‐activated secondary alkanes. Previously, this reactivity has only been observed in the gas phase for small, coordinatively unsaturated metal complexes. Multireference ab‐initio calculations suggest that 2, formally a cobalt(IV)‐oxo complex, is best described as cobalt(III)‐oxyl. Our results provide important data on changes to metal‐oxo bonding behind the oxo wall and show that cobalt‐oxo complexes are promising targets for developing highly active C−H oxidation catalysts.