Relativistic Coupled Cluster Calculations with Variational Quantum Electrodynamics Resolve the Discrepancy between Experiment and Theory Concerning the Electron Affinity and Ionization Potential of Gold

14-Membered Macrocyclic β-Diiminato Gold(II) - A New Member for the Gold(II) Complex Family?

The chemistry of molecular gold compounds is dominated by the oxidation states +I and +III. For the intermediate oxidation state +II with 5d9 electron configuration, dimerization or disproportionation of the gold(II) radicals is favored, so that only a few mononuclear gold(II) complexes have been isolated to date. The present study addresses the one-electron reduction of the macrocyclic gold(III) complex [AuIIIL]+ of the innocent β-diiminato ligand L2- with a 14-membered macrocycle (L2-=5,7,12,14-tetramethyl-1,4,8,11-tetraazacyclotetradeca-5,7,12,14-tetraenato). Electrochemistry, spectroelectrochemistry and chemical reduction of [AuIIIL]+ monitored by UV/Vis, NMR and EPR spectroscopy together with density functional theory calculations reveal disproportionation of the initially generated but elusive gold(II) complex AuIIL and provide guidelines for prospective stable mononuclear tetraazamacrocyclic gold(II) complexes.

Bound-State Relativistic Quantum Electrodynamics: A Perspective for Precision Physics with Atoms and Molecules

Precision physics aims to use atoms and molecules to test and develop the fundamental theory of matter, possibly beyond the Standard Model. Most of the atomic and molecular phenomena are described by the quantum electrodynamics (QED) sector of the Standard Model. Do we have the computational tools, algorithms, and practical equations for the most possible complete computation of atoms and molecules within the QED sector? What is the fundamental equation to start with? Is it still Schrödinger's wave equation for molecular matter, or is there anything beyond that? This paper provides a concise overview of the relativistic QED framework and recent numerical developments targeting precision physics and spectroscopy applications with common features of the robust and successful relativistic quantum chemistry methodology.

Unlocking the catalytic potential of gold(II) complexes: a comprehensive reassessment

Gold(II) complexes, unlike their gold(I) and gold(III) counterparts, have been sparsely employed in the field of catalysis. This is primarily due to the challenges associated with isolating and characterising these open-shell species. However, these complexes offer a wide range of possibilities. On one hand, this intermediate oxidation state has proven to be more easily accessible through reduction and oxidation processes compared to the gold(I)/gold(III) redox couple, thereby facilitating potential homo-coupling and cross-coupling reactions. On the other hand, gold(II) exhibits Lewis acid behaviour, bridging the characteristics of the soft acid gold(I) and the hard acid gold(III). In this review, we focus on mono- and dinuclear gold(II) complexes, whether they are isolated and well-studied or proposed as intermediates in cross-coupling reactions induced by the action of oxidants or light. We delve into the unique reactivity and potential applications of these gold(II) species, shedding light on their role in this field. This comprehensive exploration aims to underscore the latent promise of gold(II) complexes in catalysis, offering insights into their structural and mechanistic aspects while highlighting their relevance in contemporary chemical transformations.

Effective quantum electrodynamics: One-dimensional model of the relativistic hydrogen-like atom

We consider a one-dimensional effective quantum electrodynamics (QED) model of the relativistic hydrogen-like atom using delta-potential interactions. We discuss the general exact theory and the Hartree-Fock approximation. The present one-dimensional effective QED model shares the essential physical feature of the three-dimensional theory: the nuclear charge polarizes the vacuum state (creation of electron-positron pairs), which results in a QED Lamb-type shift of the bound-state energy. Yet, this 1D effective QED model eliminates some of the most serious technical difficulties of the three-dimensional theory coming from renormalization. We show how to calculate the vacuum-polarization density at zeroth order in the two-particle interaction and the QED Lamb-type shift of the bound-state energy at first order in the two-particle interaction. The present work may be considered a step toward the development of a quantum-chemistry effective QED theory of atoms and molecules.

Rigorous Negative Ion Binding Energies in Low-Energy Electron Elastic Collisions with Heavy Multi-Electron Atoms and Fullerene Molecules: Validation of Electron Affinities

Dramatically sharp resonances manifesting stable negative ion formation characterize Regge pole-calculated low-energy electron elastic total cross sections (TCSs) of heavy multi-electron systems. The novelty of the Regge pole analysis is in the extraction of rigorous and unambiguous negative ion binding energies (BEs), corresponding to the measured electron affinities (EAs) of the investigated multi-electron systems. The measured EAs have engendered the crucial question: is the EA of multi-electron atoms and fullerene molecules identified with the BE of the attached electron in the ground, metastable or excited state of the formed negative ion during a collision? Inconsistencies in the meaning of the measured EAs are elucidated and new EA values for Bk, Cf, Fm, and Lr are presented.

Diffuse Basis Functions for Relativistic s and d Block Gaussian Basis Sets

Diffuse s, p, and d functions have been optimized for use with previously reported relativistic basis sets for the s and d blocks of the periodic table. The functions were optimized on the 4:1 weighted average of the s² and p² configurations of the anion, with the d shell in the dⁿ⁺¹ configuration for the d blocks. Exponents were extrapolated for groups 2 and 12, which have unstable or weakly bound anions. The diffuse basis sets have been tested by application to calculations of electron affinities of the group 11 elements (Cu, Ag, and Au), double electron affinities of the group 11 monocations, and potential energy curves of Mg₂ and Ca₂ van der Waals dimers, as well as some response properties of the group 1 anions (Rb^-, Cs^-, and Fr^-), the group 2 elements (Sr, Ba, and Ra), and RbLi, CsLi, and FrLi molecules.

First-Principles Study of Stability and N2 Activation on the Octahedron RuRh Clusters

The geometric and electronic structures of different octahedron RuRh clusters are studied using density functional theory calculations. The binding energy, electronic structure, and energy gap of the clusters have been obtained to determine the possible stable structures. The results show that the Ru4Rh2 cluster is the most stable structure which has D4h symmetry with the largest ionization potential, smallest affinity energy and larger energy gap. Furthermore, the information on adsorption and dissociation of multiple nitrogen molecules and the density of state for the octahedral Ru4Rh2 cluster is analyzed. The dissociation barrier of three nitrogen molecules further decreases to 1.18 eV with an increase in the number of N2 molecules. The co-adsorption of multiple N2 molecules facilitates the dissociation of N2 on the Ru4Rh2 cluster. The strong interaction between the antibonding orbital of N2 and the d orbital of the Ru4Rh2 cluster is illustrated by calculating and analyzing the results of PDOS, which stretches the N−N bond length and reduces the activation energy to dissociation. The antibonding orbital of the nitrogen molecule shows distinct and unique catalytic activity for the dissociation of the adsorbed nitrogen molecule on the octahedral Ru4Rh2 cluster.

Electronic Structure of Lr+ (Z = 103) from Ab Initio Calculations

The four-component relativistic Dirac–Coulomb Hamiltonian and the multireference configuration interaction (MRCI) model were used to provide the reliable energy levels and spectroscopic properties of the Lr+ ion and the Lu+ homolog. The energy spectrum of Lr+ is very similar to that of the Lu+ homolog, with the multiplet manifold of the 7s2, 6d17s1 and 7s17p1 configurations as the ground and low-lying excited states. The results are discussed in light of earlier findings utilizing different theoretical models. Overall, the MRCI model can reliably predict the energy levels and properties and bring new insight into experiments with superheavy ions.

Accurate ab initio calculations of RaF electronic structure appeal to more laser-spectroscopical measurements

Recently, a breakthrough has been achieved in laser-spectroscopic studies of short-lived radioactive compounds with the first measurements of the radium monofluoride molecule (RaF) UV/vis spectra. We report results from high-accuracy ab initio calculations of the RaF electronic structure for ground and low-lying excited electronic states. Two different methods agree excellently with experimental excitation energies from the electronic ground state to the ²Π_1/2 and ²Π_3/2 states, but lead consistently and unambiguously to deviations from experimental-based adiabatic transition energy estimates for the ²Σ_1/2 excited electronic state, and show that more measurements are needed to clarify spectroscopic assignment of the ²Δ state.

The role of QED effects in transition energies of heavy-atom alkaline earth monofluoride molecules: A theoretical study of Ba+, BaF, RaF, and E120F

Heavy-atom alkaline earth monofluoride molecules are considered as prospective systems to study spatial parity or spatial parity and time-reversal symmetry violating effects such as the nuclear anapole moment or the electron electric dipole moment. A comprehensive and highly accurate theoretical study of the electronic structure properties and transition energies in such systems can simplify the preparation and interpretation of the experiments. However, almost no attempts to calculate quantum electrodynamics (QED) effects' contribution into characteristics of these neutral heavy-atom molecules have been performed. Recently, we have formulated and implemented such an approach to calculate QED contributions to transition energies of molecules [L. V. Skripnikov, J. Chem. Phys. 154, 201101 (2021)]. In this paper, we perform a benchmark theoretical study of the transition energies in the Ba⁺ cation and BaF molecule. The deviation of the calculated values from the experimental ones is of the order 10 cm^-1 and is more than an order of magnitude better than the "chemical accuracy," 350 cm^-1. The achievement of such an agreement has been provided, in particular, by the inclusion of the QED effects. The latter appeared to be not less important than the high-order correlation effects beyond the coupled cluster with single, double, and perturbative triple cluster amplitude level. We compare the role of QED effects for transition energies with heavier molecules-RaF and E120F, where E120 is the superheavy Z = 120 homolog of Ra.

Coupled Cluster Studies of Platinum-Actinide Interactions. Thermochemistry of PtAnOn+ (n = 0-2 and An = U, Np, Pu)

Accurate Pt-An bond dissociation enthalpies (BDEs) for PtAnOⁿ⁺ (An = U, Np, Pu and n = 0-2) and the corresponding enthalpies for the Pt + OAnOⁿ⁺ substitution reactions have been studied for the first time using an accurate composite coupled cluster approach. Analogous O-AnOⁿ⁺ bond dissociation enthalpies are also presented. To make the study possible, new correlation consistent basis sets optimized using the all-electron third-order Douglas-Kroll-Hess (DKH3) scalar relativistic Hamiltonian are developed and reported for Pt and Au, with accompanying benchmark calculations of their atomic ionization potentials to demonstrate the effectiveness of the new basis sets. For the charged PtAnOⁿ⁺ species (n = 1, 2), a low-spin state (LSS) for which the Pt-An σ bond is doubly occupied is studied together with a high-spin state (HSS) obtained by unpairing the σ bond orbital and placing one electron into the An 5f shell. The relative energies of the two spin states have been compared and qualitatively assessed via natural population and natural bond analyses. The enthalpies for the Pt substitution reactions, i.e., Pt + OAnOⁿ⁺ → PtAnOⁿ⁺ + O, are calculated to range from about 14-62 kcal/mol, and the Pt-AnOⁿ⁺ bond dissociation enthalpies range from about 78-149 kcal/mol for the ground electronic states. For the PtAnO⁺ species, the LSSs were all predicted to be the ground state, whereas the PtAnO²⁺ molecules all favored the HSSs. The prediction for PtUO²⁺ is consistent with previous theoretical findings. The natural bond orbital analyses indicate a triple bond between An and O, with a double to quadruple bond between the An and Pt.

Approaching meV level for transition energies in the radium monofluoride molecule RaF and radium cation Ra+ by including quantum-electrodynamics effects

Highly accurate theoretical predictions of transition energies in the radium monofluoride molecule, ²²⁶RaF, and radium cation, ²²⁶Ra⁺, are reported. The considered transition X²Σ_1/2 → A²Π_1/2 in RaF is one of the main features of this molecule and can be used to laser-cool RaF for a subsequent measurement of the electron electric dipole moment. For molecular and atomic predictions, we go beyond the Dirac-Coulomb Hamiltonian and treat high-order electron correlation effects within the coupled cluster theory with the inclusion of quadruple and ever higher amplitudes. The effects of quantum electrodynamics (QED) are included non-perturbatively using the model QED operator that is now implemented for molecules. It is shown that the inclusion of the QED effects in molecular and atomic calculations is a key ingredient in resolving the discrepancy between the theoretical values obtained within the Dirac-Coulomb-Breit Hamiltonian and the experiment. The remaining deviation from the experimental values is within a few meV. This is more than an order of magnitude better than the "chemical accuracy," 1 kcal/mol = 43 meV, that is usually considered as a guiding thread in theoretical molecular physics.

Ab initio Calculations of Energy Levels in Be-Like Xenon: Strong Interference between Electron-Correlation and QED Effects

The strong mixing of close levels with two valence electrons in Be-like xenon greatly complicates ab initio QED calculations beyond the first-order approximation. Because of a strong interplay between the electron-electron correlation and QED effects, the standard single-level perturbative QED approach may fail, even if it takes into account the second-order screened QED diagrams. In the present Letter, the corresponding obstacles are overcome by working out the QED perturbation theory for quasidegenerate states. The contributions of all the Feynman diagrams up to the second order are taken into account. The many-electron QED effects are rigorously evaluated in the framework of the extended Furry picture to all orders in the nuclear-strength parameter αZ. The higher-order electron-correlation effects are considered within the Breit approximation. The nuclear recoil effect is accounted for as well. The developed approach is applied to high-precision QED calculations of the ground and singly excited energy levels in Be-like xenon. The most accurate theoretical predictions for the binding and excitation energies are obtained. These results deviate from the most precise experimental value by 3σ but perfectly agree with a more recent measurement.

Compound-tunable embedding potential: which oxidation state of uranium and thorium as point defects in xenotime is favorable?

Modern strategies for the safe handling of high level waste (HLW) and its long-term disposal in deep geological formations include the immobilization of radionuclides in the form of mineral-like matrices. The most promising matrices for the immobilization of actinides are ceramic forms of waste based on phosphate minerals such as monazite, xenotime, and cheralite. However, the mechanism of substitution of lanthanides and Y by actinides in phosphate minerals is not entirely clear. We formulated a theoretical model, compound-tunable embedding potential (CTEP), that allows one to predict properties of such crystals with point defects. The reliability of the model is validated by a good agreement of calculated geometry parameters with available experimental data. The substitution of Y in the xenotime crystal by Th and U is studied by relativistic DFT in the framework of the CTEP method, based on constructing the embedding potential as the linear combination of short-range "electron-free" spherical "tunable" pseudopotentials. It is shown on the basis of the proposed model that oxidation state +3 is energetically more profitable than +4 not only for thorium but also for uranium as solitary point defects. This atypical oxidation state of U in the mineral is discussed.

The electron affinity of astatine

One of the most important properties influencing the chemical behavior of an element is the electron affinity (EA). Among the remaining elements with unknown EA is astatine, where one of its isotopes, 211At, is remarkably well suited for targeted radionuclide therapy of cancer. With the At- anion being involved in many aspects of current astatine labeling protocols, the knowledge of the electron affinity of this element is of prime importance. Here we report the measured value of the EA of astatine to be 2.41578(7) eV. This result is compared to state-of-the-art relativistic quantum mechanical calculations that incorporate both the Breit and the quantum electrodynamics (QED) corrections and the electron-electron correlation effects on the highest level that can be currently achieved for many-electron systems. The developed technique of laser-photodetachment spectroscopy of radioisotopes opens the path for future EA measurements of other radioelements such as polonium, and eventually super-heavy elements.

Relativistic Fock Space Coupled Cluster Method for Many-Electron Systems: Non-Perturbative Account for Connected Triple Excitations

The Fock space relativistic coupled cluster method (FS-RCC) is one of the most promising tools of electronic structure modeling for atomic and molecular systems containing heavy nuclei. Until recently, capabilities of the FS-RCC method were severely restricted by the fact that only single and double excitations in the exponential parametrization of the wave operator were considered. We report the design and the first computer implementation of FS-RCC schemes with full and simplified non-perturbative account for triple excitations in the cluster operator. Numerical stability of the new computational scheme and thus its applicability to a wide variety of molecular electronic states is ensured using the dynamic shift technique combined with the extrapolation to zero-shift limit. Pilot applications to atomic (Tl, Pb) and molecular (TlH) systems reported in the paper indicate that the breakthrough in accuracy and predictive power of the electronic structure calculations for heavy-element compounds can be achieved. Moreover, the described approach can provide a firm basis for high-precision modeling of heavy molecular systems with several open shells, including actinide compounds.

The periodic table and the physics that drives it

Mendeleev's introduction of the periodic table of elements is one of the most important milestones in the history of chemistry, as it brought order into the known chemical and physical behaviour of the elements. The periodic table can be seen as parallel to the Standard Model in particle physics, in which the elementary particles known today can be ordered according to their intrinsic properties. The underlying fundamental theory to describe the interactions between particles comes from quantum theory or, more specifically, from quantum field theory and its inherent symmetries. In the periodic table, the elements are placed into a certain period and group based on electronic configurations that originate from the Pauli and Aufbau principles for the electrons surrounding a positively charged nucleus. This order enables us to approximately predict the chemical and physical properties of elements. Apparent anomalies can arise from relativistic effects, partial-screening phenomena (of type lanthanide contraction) and the compact size of the first shell of every l-value. Further, ambiguities in electron configurations and the breakdown of assigning a dominant configuration, owing to configuration mixing and dense spectra for the heaviest elements in the periodic table. For the short-lived transactinides, the nuclear stability becomes an important factor in chemical studies. Nuclear stability, decay rates, spectra and reaction cross sections are also important for predicting the astrophysical origin of the elements, including the production of the heavy elements beyond iron in supernova explosions or neutron-star mergers. In this Perspective, we critically analyse the periodic table of elements and the current status of theoretical predictions and origins for the heaviest elements, which combine both quantum chemistry and physics.

Relativistic double-ionization equation-of-motion coupled-cluster method: Application to low-lying doubly ionized states

This article deals with the extension of the relativistic double-ionization equation-of-motion coupled-cluster (DI-EOMCC) method [H. Pathak et al. Phys. Rev. A 90, 010501(R) (2014)] for the molecular systems. The Dirac-Coulomb Hamiltonian with four-component spinors is considered to take care of the relativistic effects. The implemented method is employed to compute a few low-lying doubly ionized states of noble gas atoms (Ar, Kr, Xe, and Rn) and Cl₂, Br₂, HBr, and HI. Additionally, we presented results with two intermediate schemes in the four-component relativistic DI-EOMCC framework to understand the role of electron correlation. The computed double ionization spectra for the atomic systems are compared with the values from the non-relativistic DI-EOMCC method with spin-orbit coupling [Z. Wang et al. J. Chem. Phys. 142, 144109 (2015)] and the values from the National Institute of Science and Technology (NIST) database. Our atomic results are found to be in good agreement with the NIST values. Furthermore, the obtained results for the molecular systems agree well with the available experimental values.

An essay on periodic tables

After a compact history of the PT, from Döbereiner’s triads to the theoretical predictions up to element 172, a number of particular issues is discussed: Why may Z = 172 be a limit for stable electron shells? What are the expected stability limits of the nuclear isotopes? When are formally empty atomic orbitals used in molecular electronic structures? What is ‘Secondary Periodicity’? When do the elements (Ir, Pt, Au), at the end of a bond, simulate (N, O, I), respectively? Some new suggestions for alternative PTs are commented upon. As a local connection, Johan Gadolin’s 1794 analysis of the Ytterby mineral is mentioned.

QED effects on individual atomic orbital energies

Several issues, concerning QED corrections, that are important in precise atomic calculations are presented. The leading QED corrections, self-energy and vacuum polarization, to the orbital energy for selected atoms with 30 ≤ Z ≤ 118 have been calculated. The sum of QED and Breit contributions to the orbital energy is analyzed. It has been found that for ns subshells the Breit and QED contributions are of comparative size, but for np and nd subshells the Breit contribution takes a major part of the QED+Breit sum. It has also, been found that the Breit to leading QED contributions ratio for ns subshells is almost independent of Z. The Z-dependence of QED and Breit+QED contributions per subshell is shown. The fitting coefficients may be used to estimate QED effects on inner molecular orbitals. We present results of our calculations for QED contributions to orbital energy of valence ns-subshell for group 1 and 11 atoms and discuss about the reliability of these numbers by comparing them with experimental first ionization potential data.

Relativistic Normal Coupled-Cluster Theory for Accurate Determination of Electric Dipole Moments of Atoms: First Application to the ^{199}Hg Atom

Recent relativistic coupled-cluster (RCC) calculations of electric dipole moments (EDMs) of diamagnetic atoms due to parity and time-reversal violating (P,T-odd) interactions, which are essential ingredients for probing new physics beyond the standard model of particle interactions, differ substantially from the previous theoretical results. It is therefore necessary to perform an independent test of the validity of these results. In view of this, the normal coupled-cluster method has been extended to the relativistic regime [relativistic normal coupled-cluster (RNCC) method] to calculate the EDMs of atoms by simultaneously incorporating the electrostatic and P,T-odd interactions in order to overcome the shortcomings of the ordinary RCC method. This new relativistic method has been applied to ^{199}Hg, which currently has a lower EDM limit than that of any other system. The results of our RNCC and self-consistent RCC calculations of the EDM of this atom are found to be close. The discrepancies between these two results on the one hand and those of previous calculations on the other are elucidated. Furthermore, the electric dipole polarizability of this atom, which has computational similarities with the EDM, is evaluated and it is in very good agreement with its measured value.

Analytic Solutions to Coherent Control of the Dirac Equation

A simple framework for Dirac spinors is developed that parametrizes admissible quantum dynamics and also analytically constructs electromagnetic fields, obeying Maxwell's equations, which yield a desired evolution. In particular, we show how to achieve dispersionless rotation and translation of wave packets. Additionally, this formalism can handle control interactions beyond electromagnetic. This work reveals unexpected flexibility of the Dirac equation for control applications, which may open new prospects for quantum technologies.

Density functional theory modeling of C–Au chemical bond formation in gold implanted polyethylene

We have studied processes of gold ion implantation in polyethylene (PE) by theoretical chemistry methods.