Opportunities for fundamental physics research with radioactive molecules.

Molecules containing short-lived, radioactive nuclei are uniquely positioned to enable a wide range of scientific discoveries in the areas of fundamental symmetries, astrophysics, nuclear structure, and chemistry. Recent advances in the ability to create, cool, and control complex molecules down to the quantum level, along with recent and upcoming advances in radioactive species production at several facilities around the world, create a compelling opportunity to coordinate and combine these efforts to bring precision measurement and control to molecules containing extreme nuclei. In this manuscript, we review the scientific case for studying radioactive molecules, discuss recent atomic, molecular, nuclear, astrophysical, and chemical advances which provide the foundation for their study, describe the facilities where these species are and will be produced, and provide an outlook for the future of this nascent field.

Nuclear Charge Radii of Silicon Isotopes.

The nuclear charge radius of ^{32}Si was determined using collinear laser spectroscopy. The experimental result was confronted with ab initio nuclear lattice effective field theory, valence-space in-medium similarity renormalization group, and mean field calculations, highlighting important achievements and challenges of modern many-body methods. The charge radius of ^{32}Si completes the radii of the mirror pair ^{32}Ar-^{32}Si, whose difference was correlated to the slope L of the symmetry energy in the nuclear equation of state. Our result suggests L≤60 MeV, which agrees with complementary observables.

Systematic study and uncertainty evaluation of P, T-odd molecular enhancement factors in BaF.

A measurement of the magnitude of the electric dipole moment of the electron (eEDM) larger than that predicted by the Standard Model (SM) of particle physics is expected to have a huge impact on the search for physics beyond the SM. Polar diatomic molecules containing heavy elements experience enhanced sensitivity to parity (P) and time-reversal (T)-violating phenomena, such as the eEDM and the scalar-pseudoscalar (S-PS) interaction between the nucleons and the electrons, and are thus promising candidates for measurements. The NL-eEDM collaboration is preparing an experiment to measure the eEDM and S-PS interaction in a slow beam of cold BaF molecules [P. Aggarwal et al., Eur. Phys. J. D 72, 197 (2018)]. Accurate knowledge of the electronic structure parameters, Wd and Ws, connecting the eEDM and the S-PS interaction to the measurable energy shifts is crucial for the interpretation of these measurements. In this work, we use the finite field relativistic coupled cluster approach to calculate the Wd and Ws parameters in the ground state of the BaF molecule. Special attention was paid to providing a reliable theoretical uncertainty estimate based on investigations of the basis set, electron correlation, relativistic effects, and geometry. Our recommended values of the two parameters, including conservative uncertainty estimates, are 3.13 ±0.12×1024Hzecm for Wd and 8.29 ± 0.12 kHz for Ws.

Hyperfine Structure Constants on the Relativistic Coupled Cluster Level with Associated Uncertainties.

Accurate predictions of hyperfine structure (HFS) constants are important in many areas of chemistry and physics, from the determination of nuclear electric and magnetic moments to benchmarking of new theoretical methods. We present a detailed investigation of the performance of the relativistic coupled cluster method for calculating HFS constants within the finite-field scheme. The two selected test systems are 133Cs and 137BaF. Special attention has been paid to construct a theoretical uncertainty estimate based on investigations on basis set, electron correlation and relativistic effects. The largest contribution to the uncertainty estimate comes from higher order correlation contributions. Our conservative uncertainty estimate for the calculated HFS constants is ∼5.5%, while the actual deviation of our results from experimental values is <1% in all cases.

Enhanced P,T-violating nuclear magnetic quadrupole moment effects in laser-coolable molecules.

Nuclear magnetic quadrupole moments (MQMs), such as intrinsic electric dipole moments of elementary particles, violate both parity and time-reversal symmetry and, therefore, probe physics beyond the standard model. We report on accurate relativistic coupled cluster calculations of the nuclear MQM interaction constants in BaF, YbF, BaOH, and YbOH. We elaborate on estimates of the uncertainty of our results. The implications of experiments searching for nonzero nuclear MQMs are discussed.

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.

First Ionization Potentials of Fm, Md, No, and Lr: Verification of Filling-Up of 5f Electrons and Confirmation of the Actinide Series.

We report the first ionization potentials (IP1) of the heavy actinides, fermium (Fm, atomic number Z = 100), mendelevium (Md, Z = 101), nobelium (No, Z = 102), and lawrencium (Lr, Z = 103), determined using a method based on a surface ionization process coupled to an online mass separation technique in an atom-at-a-time regime. The measured IP1 values agree well with those predicted by state-of-the-art relativistic calculations performed alongside the present measurements. Similar to the well-established behavior for the lanthanides, the IP1 values of the heavy actinides up to No increase with filling up the 5f orbital, while that of Lr is the lowest among the actinides. These results clearly demonstrate that the 5f orbital is fully filled at No with the [Rn]5f147s2 configuration and that Lr has a weakly bound electron outside the No core. In analogy to the lanthanide series, the present results unequivocally verify that the actinide series ends with Lr.

Electronic structure of three-dimensional isotropic quantum dots by four-component relativistic coupled cluster methods.

Quantum dots with three-dimensional isotropic harmonic confining potentials and up to 60 electrons are studied. The Dirac-Coulomb Hamiltonian serves as a framework, so that relativistic effects are included, and electron correlation is treated at a high level by the Fock-space coupled cluster method, with single and double excitations summed to all orders. Large basis sets composed of spherical Gaussian functions are used. Energies of ground and excited states are calculated. The orbital order is 1s, 2p, 3d, 3s, 4f, 4p, 5g, ... , and closed-shell structures appear for 2, 8, 18, 20, 34, 40, and 58 electrons. Relativistic effects are negligible for low strengths of the harmonic potential and increase rapidly for stronger potentials. Breit contributions, coming from the lowest order relativistic correction to the interelectronic repulsion terms, are also studied. Correlation effects are significant for these systems, in particular for weak confining potentials and for small systems, where they constitute up to 6% of the total energies. Their relative weight goes down (although they increase in absolute value) for larger systems or confining potentials. Planned applications to quantum dots with impurities are discussed briefly.

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.

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.

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.; Wahlgren, U. J. Chem. Phys. 2007, 126, 214302], a sizable difference between CASPT2 and linear response coupled cluster (LRCC) was observed both with and without the perturbative inclusion of spin-orbit coupling. This poses a serious question as to which of the two would be more reliable for investigating molecules containing actinides. In this paper we address this question by comparing CASPT2 and LRCC to a method known to accurately describe the spectra of actinide-containing molecules: the four-component intermediate Hamiltonian Fock-space coupled cluster (IHFSCC) method, where electron correlation and spin-orbit coupling are treated on an equal footing. Our results indicate that for UO22+ there is little difference between treatments of spin-orbit coupling, making electron correlation the main cause of discrepancies. We have found IHFSCC and LRCC to be the most alike in the overall description of excited states, even though individual LRCC energies are blue-shifted in comparison to IHFSCC due to the difference in the parametrization of the excited states' wave functions. CASPT2, on the other hand, shows good agreement with IHFSCC for individual frequencies but significantly less so for the spectrum as a whole, due to the difference in the degree of correlation recovered in both cases.

On the performance of two-component energy-consistent pseudopotentials in atomic Fock-space coupled cluster calculations.

The four-component atomic intermediate-Hamiltonian Fock-space coupled cluster (IHFSCC) code of Landau et al. [J. Chem. Phys. 115, 6862 (2001)] has been adapted to two-component calculations with relativistic pseudopotentials of the energy-consistent variety. Recently adjusted energy-consistent pseudopotentials for group 11 and 12 transition elements as well as group 13 and 14 post-d main group elements, which were fitted to atomic valence spectra from four-component multiconfiguration Dirac-Hartree-Fock calculations, are tested in IHFSCC calculations for ionization potentials, electron affinities, and excitation energies of a variety of atoms and ions. Where comparison is possible, the deviations from experimental data are in good agreement with those found in previously published IHFSCC all-electron calculations: experimental data are usually reproduced within a few hundred wavenumbers.

Nuclear quadrupole moment of 197Au from high-accuracy atomic calculations.

The electric field gradient (EFG) at the gold nucleus is calculated using a finite field approach, to make the extraction of the nuclear quadrupole moment Q from experimental nuclear quadrupole coupling constants possible. The four-component Dirac-Coulomb Hamiltonian serves as the framework, 51 of the 79 electrons are correlated by the relativistic Fock-space coupled cluster method with single and double excitations, and the contribution of the Gaunt term, the main part of the Breit interaction, is evaluated. Large basis sets (up to 26s22p18d12f8g5h uncontracted Gaussians) are employed. Energy splittings of the 2D5/2 and 2D3/2 levels, rather than level shifts, are used to extract the EFG constants, as the former remain linear with Q up to 10(-5) a.u., whereas the latter display significant nonlinearity even at Q=10(-8) a.u. Larger Q values lead to larger energy changes and better precision. Excellent agreement (0.1%) is obtained between Q values derived from 2D5/2 and 2D3/2 data. Systematic errors connected with neglecting triple and higher excitations, truncating the basis and orbital active space, and approximating the Gaunt contribution are evaluated. The final value of Q(197Au) is 521(7) mb. It is lower than the muonic 547(16) mb and agrees within error bounds with the recent value of 510(15) mb obtained from molecular calculations.

High-accuracy calculation of nuclear quadrupole moments of atomic halogens.

Electric field gradients at the nuclei of halogen atoms are calculated using a finite field approach. The four-component Dirac-Coulomb Hamiltonian serves as the framework, all electrons are correlated by the relativistic Fock-space coupled cluster method with single and double excitations, and the Gaunt term, the main part of the Breit interaction, is included. Large basis sets (e.g., 28s24p21d9f4g2h Gaussian-type functions for I) are used. Combined with experimental nuclear quadrupole coupling constants, accurate estimates of the nuclear quadrupole moments are obtained. The calculated values are Q(35Cl)=-81.1(1.2) mb, Q(79Br)=302(5) mb, and Q(127I)=-680(10) mb. Currently accepted reference values [Pyykko, Mol. Phys. 99, 1617 (2001)] are -81.65(80), 313(3), and -710(10) mb, respectively. Our values are lower for the heavier halogens, corroborating the recent work of van Stralen and Visscher [Mol. Phys. 101, 2115 (2003)], who obtained Q(127I)=-696(12) mb in a series of molecular calculations.

A Fock space coupled cluster study on the electronic structure of the UO(2), UO(2) (+), U(4+), and U(5+) species.

The ground and excited states of the UO(2) molecule have been studied using a Dirac-Coulomb intermediate Hamiltonian Fock-space coupled cluster approach (DC-IHFSCC). This method is unique in describing dynamic and nondynamic correlation energies at relatively low computational cost. Spin-orbit coupling effects have been fully included by utilizing the four-component Dirac-Coulomb Hamiltonian from the outset. Complementary calculations on the ionized systems UO(2) (+) and UO(2) (2+) as well as on the ions U(4+) and U(5+) were performed to assess the accuracy of this method. The latter calculations improve upon previously published theoretical work. Our calculations confirm the assignment of the ground state of the UO(2) molecule as a (3)Phi(2u) state that arises from the 5f(1)7s(1) configuration. The first state from the 5f(2) configuration is found above 10,000 cm(-1), whereas the first state from the 5f(1)6d(1) configuration is found at 5,047 cm(-1).

Extrapolated intermediate Hamiltonian coupled-cluster approach: theory and pilot application to electron affinities of alkali atoms.

The intermediate Hamiltonian (IH) coupled-cluster method makes possible the use of very large model spaces in coupled-cluster calculations without running into intruder states. This is achieved at the cost of approximating some of the IH matrix elements, which are not taken at their rigorous effective Hamiltonian (EH) value. The extrapolated intermediate Hamiltonian (XIH) approach proposed here uses a parametrized IH and extrapolates it to the full EH, with model spaces larger by several orders of magnitude than those possible in EH coupled-cluster methods. The flexibility and resistance to intruders of the IH approach are thus combined with the accuracy of full EH. Various extrapolation schemes are described. A pilot application to the electron affinities (EAs) of alkali atoms is presented, where converged EH results are obtained by XIH for model spaces of approximately 20,000 determinants; direct EH calculations converge only for a one-dimensional model space. Including quantum electrodynamic effects, the average XIH error for the EAs is 0.6 meV and the largest error is 1.6 meV. A new reference estimate for the EA of Fr is proposed at 486+/-2 meV.

Mixed-sector intermediate Hamiltonian Fock-space coupled cluster approach.

An alternative formulation of the intermediate Hamiltonian Fock-space coupled cluster scheme developed before is presented. The methodological and computational advantages of the new formulation include the possibility of using a model space with determinants belonging to different Fock-space sectors. This extends the scope of application of the multireference coupled cluster method, and makes possible the use of quasiclosed shells (e.g., p2, d4) as reference states. Representative applications are described, including electron affinities of group-14 atoms, ionization potentials of group-15 elements, and ionization potentials and excitation energies of silver and gold. Excellent agreement with experiment (a few hundredths of an electronvolt) is obtained, with significant improvement (by a factor of 5-10 for p3 states) over Fock-space coupled cluster results. Many states not reachable by the Fock-space approach can now be studied.