Is the protactinium(V) mono-oxo bond weaker than what we thought?

The bond distance is the simplest and most obvious indicator of the nature of a given chemical bond. However, for rare chemistry, it may happen that it is not yet firmly established. In this communication, we will show that the formally-triple protactinium(V) mono-oxo bond is predicted to be longer than what was previously reported in the solid state and in solution, based on robust quantum mechanical calculations, supported by an extensive methodological study.

Uranium(iv) alkyl cations: synthesis, structures, comparison with thorium(iv) analogues, and the influence of arene-coordination on thermal stability and ethylene polymerization activity.

Reaction of [(XA)U(CHSiMe)] (1; XA = 4,5-bis(2,6-diisopropylanilido)-2,7-di--butyl-9,9-dimethylxanthene) with 1 equivalent of [PhC][B(CF)] in arene solvents afforded the arene-coordinated uranium alkyl cations, [(XA)U(CHSiMe)(η -arene)][B(CF)] {arene = benzene (2), toluene (3), bromobenzene (4) and fluorobenzene (5)}. Compounds 2, 3, and 5 were crystallographically characterized, and in all cases the arene is π-coordinated. Solution NMR studies of 2-5 suggest that the binding preferences of the [(XA)U(CHSiMe)] cation follow the order: toluene ≈ benzene > bromobenzene > fluorobenzene.

Frozen-Density Embedding for Including Environmental Effects in the Dirac-Kohn-Sham Theory: An Implementation Based on Density Fitting and Prototyping Techniques.

Frozen density embedding (FDE) represents an embedding scheme in which environmental effects are included from first-principles calculations by considering the surrounding system explicitly by means of its electron density. In the present paper, we extend the full four-component relativistic Dirac-Kohn-Sham (DKS) method, as implemented in the BERTHA code, to include environmental and confinement effects with the FDE scheme (DKS-in-DFT FDE). The implementation, based on the auxiliary density fitting techniques, has been enormously facilitated by BERTHA's python API (PyBERTHA), which facilitates the interoperability with other FDE implementations available through the PyADF framework.

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations.

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree-Fock virtual orbitals by natural orbitals (NOs). In this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller-Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al.

Simulating core electron binding energies of halogenated species adsorbed on ice surfaces and in solution relativistic quantum embedding calculations.

In this work, we investigate the effects of the environment on the X-ray photoelectron spectra of hydrogen chloride and chloride ions adsorbed on ice surfaces, as well as of chloride ions in water droplets. In our approach, we combine a density functional theory (DFT) description of the ice surface with that of halogen species using the recently developed relativistic core-valence separation equation of motion coupled cluster (CVS-EOM-IP-CCSD) the frozen density embedding formalism (FDE), to determine the K and L edges of chlorine. Our calculations, which incorporate temperature effects through snapshots from classical molecular dynamics simulations, are shown to reproduce experimental trends in the change of the core binding energies of Cl upon moving from a liquid (water droplets) to an interfacial (ice quasi-liquid layer) environment.

Environment Effects on X-Ray Absorption Spectra With Quantum Embedded Real-Time Time-Dependent Density Functional Theory Approaches.

In this work we implement the real-time time-dependent block-orthogonalized Manby-Miller embedding (rt-BOMME) approach alongside our previously developed real-time frozen density embedding time-dependent density functional theory (rt-TDDFT-in-DFT FDE) code, and investigate these methods' performance in reproducing X-ray absorption spectra (XAS) obtained with standard rt-TDDFT simulations, for model systems comprised of solvated fluoride and chloride ions ([X@ , X = F, Cl). We observe that for ground-state quantities such as core orbital energies, the BOMME approach shows significantly better agreement with supermolecular results than FDE for the strongly interacting fluoride system, while for chloride the two embedding approaches show more similar results. For the excited states, we see that while FDE (constrained not to have the environment densities relaxed in the ground state) is in good agreement with the reference calculations for the region around the K and L edges, and is capable of reproducing the splitting of the 1s (n + 1)p final states ( + 1 being the lowest virtual p orbital of the halides), it by and large fails to properly reproduce the 1s (n + 2)p states and misses the electronic states arising from excitation to orbitals with important contributions from the solvent.

Electronic spectra of ytterbium fluoride from relativistic electronic structure calculations.

We report an investigation of the low-lying excited states of the YbF molecule-a candidate molecule for experimental measurements of the electron electric dipole moment-with 2-component based multi-reference configuration interaction (MRCI), equation of motion coupled cluster (EOM-CCSD) and the extrapolated intermediate Hamiltonian Fock-space coupled cluster (XIHFS-CCSD). Specifically, we address the question of the nature of these low-lying states in terms of configurations containing filled or partially-filled Yb 4f shells. We show that while it does not appear possible to carry out calculations with both kinds of configurations contained in the same active space, reliable information can be extracted from different sectors of Fock space-that is, by performing electron attachment and detachment IHFS-CCSD and EOM-CCSD calculation on the closed-shell YbF and YbF species, respectively.

The DIRAC code for relativistic molecular calculations.

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module.

On the multi-reference nature of plutonium oxides: PuO, PuO, PuO and PuO(OH).

Actinide-containing complexes present formidable challenges for electronic structure methods due to the large number of degenerate or quasi-degenerate electronic states arising from partially occupied 5f and 6d shells. Conventional multi-reference methods can treat active spaces that are often at the upper limit of what is required for a proper treatment of species with complex electronic structures, leaving no room for verifying their suitability. In this work we address the issue of properly defining the active spaces in such calculations, and introduce a protocol to determine optimal active spaces based on the use of the Density Matrix Renormalization Group algorithm and concepts of quantum information theory.

Structural, dynamical, and transport properties of the hydrated halides: How do At(-) bulk properties compare with those of the other halides, from F(-) to I(-)?

The properties of halides from the lightest, fluoride (F(-)), to the heaviest, astatide (At(-)), have been studied in water using a polarizable force-field approach based on molecular dynamics (MD) simulations at the 10 ns scale. The selected force-field explicitly treats the cooperativity within the halide-water hydrogen bond networks. The force-field parameters have been adjusted to ab initio data on anion/water clusters computed at the relativistic Möller-Plesset second-order perturbation theory level of theory.

Effective bond orders from two-step spin-orbit coupling approaches: the I2, At2, IO(+), and AtO(+) case studies.

The nature of chemical bonds in heavy main-group diatomics is discussed from the viewpoint of effective bond orders, which are computed from spin-orbit wave functions resulting from spin-orbit configuration interaction calculations. The reliability of the relativistic correlated wave functions obtained in such two-step spin-orbit coupling frameworks is assessed by benchmark studies of the spectroscopic constants with respect to either experimental data, or state-of-the-art fully relativistic correlated calculations. The I2, At2, IO(+), and AtO(+) species are considered, and differences and similarities between the astatine and iodine elements are highlighted.

Communication: Relativistic Fock-space coupled cluster study of small building blocks of larger uranium complexes.

We present a study of the electronic structure of the [UO2](+), [UO2](2 +), [UO2](3 +), NUO, [NUO](+), [NUO](2 +), [NUN](-), NUN, and [NUN](+) molecules with the intermediate Hamiltonian Fock-space coupled cluster method. The accuracy of mean-field approaches based on the eXact-2-Component Hamiltonian to incorporate spin-orbit coupling and Gaunt interactions are compared to results obtained with the Dirac-Coulomb Hamiltonian. Furthermore, we assess the reliability of calculations employing approximate density functionals in describing electronic spectra and quantities useful in rationalizing Uranium (VI) species reactivity (hardness, electronegativity, and electrophilicity).

Solvatochromic shifts from coupled-cluster theory embedded in density functional theory.

Building on the framework recently reported for determining general response properties for frozen-density embedding [S. Höfener, A. S.

Towards systematically improvable models for actinides in condensed phase: the electronic spectrum of uranyl in Cs2UO2Cl4 as a test case.

In this work we explore the use of frozen density embedding [Gomes et al., Phys. Chem.

The electronic spectrum of CUONg4 (Ng = Ne, Ar, Kr, Xe): new insights in the interaction of the CUO molecule with noble gas matrices.

The electronic spectrum of the CUO molecule was investigated with the IHFSCC-SD (intermediate Hamiltonian Fock-space coupled cluster with singles and doubles) method and with TD-DFT (time-dependent density functional theory) employing the PBE and PBE0 exchange-correlation functionals. The importance of both spin-orbit coupling and correlation effects on the low-lying excited-states of this molecule are analyzed and discussed. Noble gas matrix effects on the energy ordering and vibrational frequencies of the lowest electronic states of the CUO molecule were investigated with density functional theory (DFT) and TD-DFT in a supermolecular as well as a frozen density embedding (FDE) subsystem approach.

Molecular properties via a subsystem density functional theory formulation: a common framework for electronic embedding.

In this article, we present a consistent derivation of a density functional theory (DFT) based embedding method which encompasses wave-function theory-in-DFT (WFT-in-DFT) and the DFT-based subsystem formulation of response theory (DFT-in-DFT) by Neugebauer [J. Neugebauer, J. Chem.

PyADF--a scripting framework for multiscale quantum chemistry.

Applications of quantum chemistry have evolved from single or a few calculations to more complicated workflows, in which a series of interrelated computational tasks is performed. In particular multiscale simulations, which combine different levels of accuracy, typically require a large number of individual calculations that depend on each other. Consequently, there is a need to automate such workflows.

Electronic spectroscopy of UO2(2+), NUO(+) and NUN: an evaluation of time-dependent density functional theory for actinides.

The performance of the time-dependent density functional theory (TDDFT) approach has been evaluated for the electronic spectrum of the UO(2)(2+), NUO(+) and NUN molecules. Different exchange-correlation functionals (LDA, PBE, BLYP, B3LYP, PBE0, M06, M06-L, M06-2X, CAM-B3LYP) and the SAOP model potential have been investigated, as has the relative importance of the adiabatic local density approximation (ALDA) to the exchange-correlation kernel. The vertical excitation energies have been compared with reference data obtained using accurate wave-function theory (WFT) methods.

Analysis of parity violation in chiral molecules.

In order to guide the experimental search for parity violation in molecular systems, in part motivated by the possible link to biomolecular homochirality, we present a detailed analysis in a relativistic framework of the mechanism behind the tiny energy difference between enantiomers induced by the weak force. A decomposition of the molecular expectation value into atomic contributions reveals that the effect can be thought of as arising from a specific mixing of valence s(1/2) and p(1/2) orbitals on a single center induced by a chiral molecular field. The intra-atomic nature of the effect is further illustrated by visualization of the electron chirality density and suggests that a simple model for parity violation in molecules may be constructed by combining pre-calculated atomic quantities with simple bonding models.

The electronic structure of the triiodide ion from relativistic correlated calculations: a comparison of different methodologies.

The triiodide ion I(3)(-) exhibits a complex photodissociation behavior, the dynamics of which are not yet fully understood. As a first step toward determining the full potential energy surfaces of this species for subsequent simulations of its dissociation processes, we investigate the performance of different electronic structure methods [time-dependent density functional theory, complete active space perturbation theory to second order (CASPT2), Fock-space coupled cluster and multireference configuration interaction] in describing the ground and excited states of the triiodide ion along the symmetrical dissociation path. All methods apart from CASPT2 include scalar relativity and spin-orbit coupling in the orbital optimization, providing useful benchmark data for the more common two-step approaches in which spin-orbit coupling is introduced in the configuration interaction.

Chiral oxorhenium(V) complexes as candidates for the experimental observation of molecular parity violation: a structural, synthetic and theoretical study.

We report the synthesis and resolution of a series of new chiral "3 + 1" oxorhenium(V) complexes, designed for high-resolution laser spectroscopy experiments probing molecular parity-violation (PV) effects in the Re=O stretching mode frequency. These complexes display a particularly simple chemical structure, with the rhenium atom as the stereogenic center, and show large PV energy differences according to our calculations. They were obtained in the racemic and enantioenriched forms, in the latter case by using either semi-preparative chiral HPLC resolution or enantioselective synthesis.

Benchmarking electronic structure calculations on the bare UO2(2+) ion: how different are single and multireference electron correlation methods?

In a recent investigation by some of us on the spectrum of the uranyl (UO22+) ion [Ral, F.; Vallet, V.; Marian, C.

Calculation of local excitations in large systems by embedding wave-function theory in density-functional theory.

We present a simple and efficient embedding scheme for the wave-function based calculation of the energies of local excitations in large systems. By introducing an embedding potential obtained from density-functional theory (DFT) it is possible to describe the effect of an environment on local excitations of an embedded system in wave-function theory (WFT) calculations of the excitation energies. We outline the implementation of such a WFT-in-DFT embedding procedure employing the ADF, Dalton and DIRAC codes, where the embedded subsystem is treated with coupled cluster methods.

Analysis of the segmented contraction of basis functions using density matrix theory.

A particular formulation based on density matrix (DM) theory at the Hartree-Fock level of theory and the description of the atomic orbitals as integral transforms is introduced. This formulation leads to a continuous representation of the density matrices as functions of a generator coordinate and to the possibility of plotting either the continuous or discrete density matrices as functions of the exponents of primitive Gaussian basis functions. The analysis of these diagrams provides useful information allowing: (a) the determination of the most important primitives for a given orbital, (b) the core-valence separation, and (c) support for the development of contracted basis sets by the segmented method.