How Does Bending the Uranyl Unit Influence Its Spectroscopy and Luminescence?

Bent uranyl complexes can be formed with chloride ligands and 1,10-phenanthroline (phen) ligands bound to the equatorial and axial planes of the uranyl(VI) moiety, as revealed by the crystal structures, IR and Raman spectroscopy, and quantum-chemical calculations. With the goal of probing the influence of chloride and phenanthroline coordination enforcing the bending on the absorption and emission spectra of this complex, spin-orbit time-dependent density functional theory calculations for the bare uranyl complexes as well as for the free UO2Cl2 subunit and the UO2Cl2(phen)2 complex were performed. The emission spectra have been fully simulated by ab initio methods and compared to experimental photoluminescence spectra, recorded for the first time for UO2Cl2(phen)2. Notably, the bending of uranyl in UO2Cl2 and UO2Cl2(phen)2 triggers excitations of the uranyl bending mode, yielding a denser luminescence spectrum.

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. The accuracy and numerical stability of this new implementation, also using different auxiliary fitting basis sets, has been demonstrated on the simple NH3-H2O system, in comparison with a reference nonrelativistic implementation. The computational performance has been evaluated on a series of gold clusters (Aun, with n = 2, 4, 8) embedded into an increasing number of water molecules (5, 10, 20, 40, and 80 water molecules). We found that the procedure scales approximately linearly both with the size of the frozen surrounding environment (consistent with the underpinnings of the FDE approach) and with the size of the active system (in line with the use of density fitting). Finally, we applied the code to a series of heavy (Rn) and super-heavy elements (Cn, Fl, Og) embedded in a C60 cage to explore the confinement effect induced by C60 on their electronic structure. We compare the results from our simulations, with respect to more-approximate models employed in the atomic physics literature. Our results indicate that the specific interactions described by FDE are able to improve upon the cruder approximations currently employed, and, thus, they provide a basis from which to generate more-realistic radial potentials for confined atoms.

Simulating core electron binding energies of halogenated species adsorbed on ice surfaces and in solution via 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) via the frozen density embedding formalism (FDE), to determine the K and L1,2,3 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. Our simulations yield water valence band binding energies in good agreement with experiment, which vary little between the droplets and the ice surface. For halide core binding energies there is an overall trend for overestimating experimental values, though good agreement between theory and experiment is found for Cl- in water droplets and on ice. For HCl on the other hand there are significant discrepancies between experimental and calculated core binding energies when we consider structural models that maintain the H-Cl bond more or less intact. An analysis of models that allow for pre-dissociated and dissociated structures suggests that experimentally observed chemical shifts in binding energies between Cl- and HCl would require that H+ (in the form of H3O+) and Cl- are separated by roughly 4-6 Å.

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. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties-electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F-Ts), and parity-violating energy differences for H2Z2 (Z = O-Se)-we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

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@ ( H 2 O ) 8 - , 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 L1 edges, and is capable of reproducing the splitting of the 1s1 (n + 1)p1 final states (n + 1 being the lowest virtual p orbital of the halides), it by and large fails to properly reproduce the 1s1 (n + 2)p1 states and misses the electronic states arising from excitation to orbitals with important contributions from the solvent. The BOMME results, on the other hand, provide a faithful qualitative representation of the spectra in all energy regions considered, though its intrinsic approximation of employing a lower-accuracy exchange-correlation functional for the environment induces non-negligible shifts in peak positions for the excitations from the halide to the environment. Our results thus confirm that QM/QM embedding approaches are viable alternatives to standard real-time simulations of X-ray absorption spectra of species in complex or confined environments.

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. From these calculations we predict Ω = 1/2, 3/2 states, arising from the 4f13σ26s, 4f145d1/6p1, and 4f135d1σ16s configurations to be able to interact as they appear in the same energy range around the ground-state equilibrium geometry. As these states are generated from different sectors of Fock space, they are almost orthogonal and provide complementary descriptions of parts of the excited state manifold. To obtain a comprehensive picture, we introduce a simple adiabatization model to extract energies of interacting Ω = 1/2, 3/2 states that can be compared to experimental observations.

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. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.

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

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. We apply the protocol to elucidate the electronic structure and bonding mechanism of volatile plutonium oxides (PuO3 and PuO2(OH)2), species associated with nuclear safety issues for which little is known about the electronic structure and energetics. We show how, within a scalar relativistic framework, orbital-pair correlations can be used to guide the definition of optimal active spaces which provide an accurate description of static/non-dynamic electron correlation, as well as to analyse the chemical bonding beyond a simple orbital model. From this bonding analysis we are able to show that the addition of oxo- or hydroxo-groups to the plutonium dioxide species considerably changes the π-bonding mechanism with respect to the bare triatomics, resulting in bent structures with a considerable multi-reference character.

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. The anion static polarizabilities of the two heaviest halides, I(-) and At(-), were computed in the gas phase using large and diffuse atomic basis sets, and taking into account both electron correlation and spin-orbit coupling within a four-component framework. Our MD simulation results show the solvation properties of I(-) and At(-) in aqueous phase to be very close. For instance, their first hydration shells are structured and encompass 9.2 and 9.1 water molecules at about 3.70 ± 0.05 Å, respectively. These values have to be compared to the F(-), Cl(-), and Br(-) ones, i.e., 6.3, 8.4, and 9.0 water molecules at 2.74, 3.38, and 3.55 Å, respectively. Moreover our computations predict the solvation free energy of At(-) in liquid water at ambient conditions to be 68 kcal mol(-1), a value also close the I(-) one, about 70 kcal mol(-1). In all, our simulation results for I(-) are in excellent agreement with the latest neutron- and X-ray diffraction studies. Those for the At(-) ion are predictive, as no theoretical or experimental data are available to date.

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. In particular, we demonstrate that spin-orbit coupling weakens the covalent character of the bond in At2 even more than electron correlation, making the consideration of spin-orbit coupling compulsory for discussing chemical bonding in heavy (6p) main group element systems.

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. P. Gomes, and L. Visscher, J. Chem. Phys. 136, 044104 (2012)], in this work we report a first implementation of an embedded coupled-cluster in density-functional theory (CC-in-DFT) scheme for electronic excitations, where only the response of the active subsystem is taken into account. The formalism is applied to the calculation of coupled-cluster excitation energies of water and uracil in aqueous solution. We find that the CC-in-DFT results are in good agreement with reference calculations and experimental results. The accuracy of calculations is mainly sensitive to factors influencing the correlation treatment (basis set quality, truncation of the cluster operator) and to the embedding treatment of the ground-state (choice of density functionals). This allows for efficient approximations at the excited state calculation step without compromising the accuracy. This approximate scheme makes it possible to use a first principles approach to investigate environment effects with specific interactions at coupled-cluster level of theory at a cost comparable to that of calculations of the individual subsystems in vacuum.

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. Chem. Phys., 2008, 10, 5353] as a way to construct models of increasing sophistication for describing the low-lying electronic absorption spectra of UO2(2+) in the Cs2UO2Cl4 crystal. We find that a relatively simple embedding model, in which all but the UO2(2+) unit are represented by an embedding potential, can already describe the main spectral features and the main environment effects can be attributed to the four chloride ions situated at the UO2(2+) equatorial plane. Contributions from species further away, albeit small, are found to be important for reaching a close agreement with experimentally observed quantities such as the excited states' relative positions. These findings suggest that such an embedding approach is a viable alternative to supermolecular calculations employing larger models of actinyl species in condensed phase. Nevertheless, we observe a slight red shift of the excitation energies calculated with our models compared to experimental results, and attribute this discrepancy to inaccuracies in the underlying structural parameters.

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. This data is used to test the suitability of the FDE approach to model the influence of different matrices on the vertical electronic transitions of this molecule. The most suitable potential was chosen to perform relativistic wave function theory in density functional theory calculations to study the vertical electronic spectra of the CUO and CUONg(4) with the IHFSCC-SD method.

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. Phys. 131, 084104 (2009)] as special cases. This formulation, which is based on the time-averaged quasi-energy formalism, makes use of the variation Lagrangian techniques to allow the use of non-variational (in particular: coupled cluster) wave-function-based methods. We show how, in the time-independent limit, we naturally obtain expressions for the ground-state DFT-in-DFT and WFT-in-DFT embedding via a local potential. We furthermore provide working equations for the special case in which coupled cluster theory is used to obtain the density and excitation energies of the active subsystem. A sample application is given to demonstrate the method.

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. For this purpose we have developed PYADF, a scripting framework for quantum chemistry. PYADF handles all steps necessary in a typical workflow in quantum chemistry and is easily extensible due to its object-oriented implementation in the Python programming language. We give an overview of the capabilities of PYADF and illustrate its usefulness in quantum-chemical multiscale simulations with a number of examples taken from recent applications.

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. A 2-component relativistic computational procedure is proposed which bridges the relativistic and non-relativistic approaches to the calculation of parity violation in chiral molecules and allows us to explore the single-center theorem in a variational setting.

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. Time-dependent density functional theory with the statistical averaging of model orbital potential functional is off the mark for this system. Another choice of functional may improve performance with respect to vertical excitation energies and spectroscopic constants, but all functionals are likely to face instability problems away from the equilibrium region. The Fock-space coupled cluster method was shown to perform clearly best in regions not too far from equilibrium but is plagued by convergence problems toward the dissociation limit due to intruder states. CASPT2 shows good performance at significantly lower computational cost, but is quite sensitive to symmetry breaking. We furthermore observe spikes in the CASPT2 potential curves away from equilibrium, signaling intruder state problems that we were unable to curb through the use of level shifts. Multireference configuration interaction is, in principle, a viable option, but its computational cost in the present case prohibits use other than for benchmarking purposes.

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. The vibrational transition frequency differences between the enantiomeric pairs due to PV have been calculated with two- and four-component relativistic Hamiltonians using Hartree-Fock (HF) and density functional theory (DFT). For three complexes, including one synthesized in enantioenriched form, our HF calculations predict frequency differences above the present resolution limit of 1 Hz. These results confirm the order of magnitude for the calculated HF PV vibrational frequency differences reported earlier for this class of compounds [P. Schwerdtfeger and R. Bast, J. Am. Chem. Soc., 2004, 126, 1652]. However, at the DFT level the PV vibrational frequency differences are in some cases reduced by an order of magnitude, but are still within the sensitivity of 0.01 Hz, which is the anticipated sensitivity in a new proposed experiment. We therefore believe that the present study represents an important step towards the experimental observation of PV in molecular systems, and emphasizes the extreme sensitivity of the PV vibrational frequency difference to the chemical environment around the rhenium center.