On the multi-reference nature of plutonium oxides: PuO22+, PuO2, PuO3 and PuO2(OH)2

Delving into the catalytic mechanism of molybdenum cofactors: a novel coupled cluster study

In this work, we use modern electronic structure methods to model the catalytic mechanism of different variants of the molybdenum cofactor (Moco). We investigate the dependence of various Moco model systems on structural relaxation and the importance of environmental effects for five critical points along the reaction coordinate with the DMSO and NO3- substrates. Furthermore, we scrutinize the performance of various coupled-cluster approaches for modeling the relative energies along the investigated reaction paths, focusing on several pair coupled cluster doubles (pCCD) flavors and conventional coupled cluster approximations. Moreover, we elucidate the Mo-O bond formation using orbital-based quantum information measures, which highlight the flow of σM-O bond formation and σN/S-O bond breaking. Our study shows that pCCD-based models are a viable alternative to conventional methods and offer us unique insights into the bonding situation along a reaction coordinate. Finally, this work highlights the importance of environmental effects or changes in the core and, consequently, in the model itself to elucidate the change in activity of different Moco variants.

Quantum Information-Assisted Complete Active Space Optimization (QICAS)

We propose an effective quantum information-assisted complete active space optimization scheme (QICAS). What sets QICAS apart from other correlation-based selection schemes is (i) the use of unique measures from quantum information that assess the correlation in electronic structures in an unambiguous and predictive manner and (ii) an orbital optimization step that minimizes the correlation discarded by the active space approximation. Equipped with these features, QICAS yields, for smaller correlated molecule, sets of optimized orbitals with respect to which the complete active space configuration interaction energy reaches the corresponding complete active space self-consistent field (CASSCF) energy within chemical accuracy. For more challenging systems such as the chromium dimer, QICAS offers an excellent starting point for CASSCF by greatly reducing the number of iterations required for numerical convergence. Accordingly, our study validates a profound empirical conjecture: the energetically optimal nonactive spaces are predominantly those that contain the least entanglement.

Geminal-Based Strategies for Modeling Large Building Blocks of Organic Electronic Materials

We elaborate on unconventional electronic structure methods based on geminals and their potential to advance the rapidly developing field of organic photovoltaics (OPVs). Specifically, we focus on the computational advantages of geminal-based methods over standard approaches and identify the critical aspects of OPV development. Examples are reliable and efficient computations of orbital energies, electronic spectra, and van der Waals interactions. Geminal-based models can also be combined with quantum embedding techniques and a quantum information analysis of orbital interactions to gain a fundamental understanding of the electronic structures and properties of realistic OPV building blocks. Furthermore, other organic components present in, for instance, dye-sensitized solar cells (DSSCs) represent another promising scope of application. Finally, we provide numerical examples predicting the properties of a small building block of OPV components and two carbazole-based dyes proposed as possible DSSC sensitizers.

A configuration interaction correction on top of pair coupled cluster doubles

Numerous numerical studies have shown that geminal-based methods are a promising direction to model strongly correlated systems with low computational costs. Several strategies have been introduced to capture the missing dynamical correlation effects, which typically exploit a posteriori corrections to account for correlation effects associated with broken-pair states or inter-geminal correlations. In this article, we scrutinize the accuracy of the pair coupled cluster doubles (pCCD) method extended by configuration interaction (CI) theory. Specifically, we benchmark various CI models, including, at most double excitations against selected CC corrections as well as conventional single-reference CC methods. A simple Davidson correction is also tested. The accuracy of the proposed pCCD-CI approaches is assessed for challenging small model systems such as the N₂ and F₂ dimers and various di- and triatomic actinide-containing compounds. In general, the proposed CI methods considerably improve spectroscopic constants compared to the conventional CCSD approach, provided a Davidson correction is included in the theoretical model. At the same time, their accuracy lies between those of the linearized frozen pCCD and frozen pCCD variants.

Resolving the π-assisted U-N σf-bond formation using quantum information theory

We model the potential energy profiles of the UO₂ (NCO)Cl₂^- → NUOCl₂^- + CO₂ reaction pathway [Y. Gong, V. Vallet, M. del Carmen Michelini, D. Rios and J. K. Gibson, J. Phys. Chem. A, 2014, 118, 325-330] using different pair coupled-cluster doubles (pCCD) methods. Specifically, we focus on pCCD and pCCD-tailored coupled cluster models in predicting relative energies for the various intermediates and transition states along the reaction coordinate. Furthermore, we augment our study on energetics with an orbital-pair correlation analysis of the complete reaction pathway that features two distinct paths. Our analysis of orbital correlations sheds new light on the formation and breaking of respective bonds between the uranium, oxygen, and nitrogen atoms along the reaction coordinates where the "yl" bond is broken and a nitrido compound formed. Specifically, the strengthening of the U-N σ_f-bond is assisted by a π-type interaction that is delocalized over the C-N-U backbone of the UO₂ (NCO)Cl₂^- complex.

Assessing the Accuracy of Tailored Coupled Cluster Methods Corrected by Electronic Wave Functions of Polynomial Cost

Tailored coupled cluster theory represents a computationally inexpensive way to describe static and dynamical electron correlation effects. In this work, we scrutinize the performance of various coupled cluster methods tailored by electronic wave functions of polynomial cost. Specifically, we focus on frozen-pair coupled cluster (fpCC) methods, which are tailored by pair-coupled cluster doubles (pCCD), and coupled cluster theory tailored by matrix product state wave functions optimized by the density matrix renormalization group (DMRG) algorithm. As test system, we selected a set of various small- and medium-sized molecules containing diatomics (N₂, F₂, C₂, CN⁺, CO, BN, BO⁺, and Cr₂) and molecules (ammonia, ethylene, cyclobutadiene, benzene, hydrogen chains, rings, and cuboids) for which the conventional single-reference coupled cluster singles and doubles (CCSD) method is not able to produce accurate results for spectroscopic constants, potential energy surfaces, and barrier heights. Most importantly, DMRG-tailored and pCCD-tailored approaches yield similar errors in spectroscopic constants and potential energy surfaces compared to accurate theoretical and/or experimental reference data. Although fpCC methods provide a reliable description for the dissociation pathway of molecules featuring single and quadruple bonds, they fail in the description of triple or hextuple bond-breaking processes or avoided crossing regions.

Considering Density Functional Approaches for Actinide Species: The An66 Molecule Set

The importance of spin-orbit effects on the predictions of energetic properties of actinide compounds has been considered for 18 different density functionals, comparing the spin-orbit and non-spin-orbit ("standard") forms of density functional theory (DFT). A set of enthalpies of formation for 66 small actinide (Th-Am) compounds-the An66 set, for which experimental data are available-have been investigated. The set includes actinide halides, oxides, and oxohalides in the general form AnO_mX_n, where n = 0-6, m = 0-3, and X = F, Cl, Br, or I. The impact of basis set choice was investigated, and to help account for the impact of relativity, the Stuttgart general and segmented contracted atomic natural orbital (ANO) basis sets paired with small core relativistic effective core potentials (RECP) as well as all-electron calculations utilizing the third-order Douglas-Kroll-Hess were considered.

Hydrolysis of Small Oxo/Hydroxo Molecules Containing High Oxidation State Actinides (Th, Pa, U, Np, Pu): A Computational Study

The energetics of hydrolysis reactions for high oxidation states of oxo/hydroxo monomeric actinide species (Th^IVO₂, Pa^IVO₂, U^IVO₂, Pa^VO₂(OH), U^VO₂(OH), U^VIO₃, Np^VIO₃, Np^VIIO₃(OH), and Pu^VIIO₃(OH)) were calculated at the CCSD(T) level. The first step is the formation of a Lewis acid/base adduct with H₂O (hydration), followed by a proton transfer to form a dihydroxide molecule (hydrolysis); this process is repeated until all oxo groups are hydrolyzed. The physisorption (hydration) for each H₂O addition was predicted to be exothermic, ca. -20 kcal/mol. The hydrolysis products are preferred energetically over the hydration products for the +IV and +V oxidation states. The compounds with An^VI are a turning point in terms of favoring hydration over hydrolysis. For An^VIIO₃(OH), hydration products are preferred, and only two waters can bind; the complete hydrolysis process is now endothermic, and the oxidation state for the An in An(OH)₇ is +VI with two OH groups each having one-half an electron. The natural bond order charges and the reaction energies provide insights into the nature of the hydrolysis/hydration processes. The actinide charges and bond ionicity generally decrease across the period. The ionic character decreases as the oxidation state and coordination number increase so that covalency increases moving to the right in the actinide period.

Plutonium complexes in water: new approach to ab initio modeling

Geometry optimization and the electronic structure calculations of Pu Z+ complexes (Z = 3–6) in water solution have been performed, within the framework of the DMol3 and Relativistic Discrete-Variational (RDV) methods. For the simulation of Pu Z+ molecular environment in aqueous solution we used 22 and 32 water molecules randomly distributed around cation. To model the effect of bulk solvent environment we used COSMO (Conductor-like Screening Model) potential for water (ε = 78.54). The obtained results showed that this approach allows the modeling of water dissociation and the formation of hydrolysis products. Our previously suggested scheme for the calculation of interaction energies between selected fragments of multi-molecular systems provides the quantitative estimation of the interaction strengths between plutonium in various oxidation states and each ligand in the first and second coordination shells in water solution.

Molecular oxides of high-valent actinides

The past decade has been very productive in the field of actinide (An) oxides containing high-valent An. Novel gas-phase experimental and an impressive number of theoretical studies have been performed, mostly on pure oxides or oxides extended with other ligands. The review covers the structural properties of molecular An oxides with high (An≥V) oxidation states. The presented compounds include the actinide dioxide cations [AnO2]+ and [AnO2]2+, neutral and ionic AnOx (x = 3–6), oxides with more than one An atom like neutral dimers, trimers and dimers from cation–cation interactions, as well as large U-oxide clusters observed very recently in the gaseous phase.

Remarkably enhanced direct dissolution of plutonium oxide in task-specific ionic liquid: insights from electrochemical and theoretical investigations

The present work envisages an approach for direct dissolution of PuO2 in a task-specific ionic liquid (TSIL). An attractive possibility to electrodeposit plutonium from the mixture of TSIL and PuO2 has been explored further. The carboxyl functional group attached to the TSIL plays a key role in facilitating the dissolution of plutonium ions.

Elucidating cation–cation interactions in neptunyl dications using multi-reference ab initio theory

Understanding the binding mechanism in neptunyl clusters formed due to cation–cation interactions is of crucial importance in nuclear waste reprocessing and related areas of research.

Electron correlation effects of the ThO and ThS molecules in the spinor basis. A relativistic coupled cluster study of ground and excited states properties

We present a comprehensive relativistic coupled cluster study of the electronic structures of the ThO and ThS molecules in the spinor basis. Specifically, we use the single-reference coupled cluster and the multi-reference Fock Space Coupled Cluster (FSCC) methods to model their ground and electronically-excited states. Two variants of the FSCC method have been investigated: (a) one where the electronic spectrum is obtained from sector (1,1) of the Fock space, and (b) another where the excited states come from the doubly attached electronic states to the doubly charged systems (ThO2+ and ThS2+), that is, from sector (0,2) of the Fock space. Our study provides a reliable set of spectroscopic parameters such as bond lengths, excitation energies, and vibrational frequencies, as well as a detailed analysis of the electron correlation effects in the ThO and ThS molecules. Finally, we examine the first ionization potential and electron affinity of the above mentioned molecules.

Quantum Chemical Study of the Water Exchange Mechanism of the Neptunyl(VI) and -(V) Aqua Ions

The water exchange reaction of the neptunyl(VI) and -(V) aqua ions was investigated with quantum chemical methods. Associative ( a) and dissociative ( d) exchange mechanisms were studied. The geometries and vibrational frequencies were computed with density functional theory (DFT) and multiconfiguration self-consistent field (MCSCF) methods. The Gibbs activation energies (Δ G^‡) for a and d pathways were calculated with DFT and wave function theory (WFT), extended general multiconfiguration quasi-degenerate second-order perturbation theory (XGMC-QDPT2) including spin-orbit (SO) coupling at the SO configuration interaction (CI) level. The DFT-WFT agreement for Δ G^‡ was poor for the investigated functionals except for LC-BOP-LRD. Due to ligand-field effects, Δ G^‡ for the associatively activated exchange reaction of NpO₂(OH₂)₅²⁺ with a f_δ¹ electron configuration is higher than for the actinyl(VI) aqua ions of U, Pu, and Am exhibiting f⁰, f_δ², and f_δ²f_ϕ¹ electron configurations.