Search

An efficient approach to calculate approximate pure-state and transition reduced density matrices in the framework of the multireference relativistic Fock-space coupled cluster (FS CC) theory is proposed. The method is based on the effective operator formalism and consists of the direct substitution of the FS CC Ansatz for a wave operator into the effective operator expression with the subsequent truncation of expansion at the terms quadratic in cluster amplitudes. The final density matrix is defined by active-space density matrices of different ranks "dressed" with contributions from cluster operators. The method gives a connected expression for pure-state density matrices, provided that the intermediate normalization condition is fulfilled. Moreover, under some additional assumptions, the connectivity can also be ensured for calculated transition property matrix elements and natural transition spinors. The developed technique allows for fast and accurate calculations of one-particle reduced density matrices for a wide range of electronic states. A pilot application of the new technique to construct averaged atomic natural orbital (ANO) basis sets for fully relativistic electronic structure calculations is presented.

The calculations of the first and the second ionization potentials of lawrencium and lutetium and the electron affinity of lawrencium are performed within the relativistic coupled cluster framework. These results are corrected by including the contributions of extrapolation to the complete basis set limit and higher order contributions due to relativity and electron correlation. The excellent agreement between our predictions of the ionization potentials of Lu and Lr and experimental values supports the accuracy of our predictions of the second ionization potential and the electron affinity of Lr.

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 \leq 60$\,MeV, which agrees with complementary observables.

The tiny-core generalized (Gatchina) relativistic pseudopotential (GRPP) model provides an accurate approximation for many-electron Hamiltonians of molecules containing heavy atoms, ensuring a proper description of the effects of non-Coulombian electron-electron interactions, electronic self-energy and vacuum polarization. Combining this model with electron correlation treatment in the frames of the intermediate Hamiltonian Fock space coupled cluster theory employing incomplete main model spaces, one obtains a reliable and economical tool for excited state modeling. The performance of this method is assessed in applications to \textit{ab initio} modeling of excited electronic states of the thorium monoxide molecule with term energies below 20000 cm$^{-1}$. Radiative lifetimes of excited states are estimated using truncated expansions of effective and metric operators in powers of cluster amplitudes.

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.

A simple procedure to incorporate one-loop quantum electrodynamic (QED) corrections into the generalized (Gatchina) nonlocal shape-consistent relativistic pseudopotential model is described. The pseudopotentials for Lu, Tl, and Ra replacing only inner core shells (with principal quantum numbers $n\le 3$ for the two former elements and $n\le 4$ for the latter one) are derived from the solutions of reference atomic SCF problems with the Dirac-Coulomb-Breit Hamiltonian to which the model Lamb shift operator added. QED contributions to atomic valence excitation energies evaluated at the SCF level are demonstrated to exceed the errors introduced by the pseudopotential approximation itself by an order of magnitude. Pilot applications of the new model to calculations of excitation energies of two-valence-electron atomic systems using the intermediate-Hamiltonian relativistic Fock space coupled cluster method reformulated here for incomplete main model spaces are reported. Implications for high-accuracy molecular excited state calculations are discussed.

High-precision tests of fundamental symmetries are looking for the parity- (P), time-reversal- (T) violating electric dipole moment of the electron (eEDM) as proof of physics beyond the Standard Model. Particularly, in polyatomic molecules, the complex vibrational and rotational structure gives the possibility to reach high enhancement of the P, T-odd effects in moderate electric fields. Additionally, it is possible to increase the statistical sensitivity by using laser cooling. In this work, we calculate the P, T-odd electronic structure parameters $W_\mathrm{d}$ and $W_\mathrm{s}$ for the promising candidates BaCH$_3$ and YbCH$_3$ for the interpretation of future experiments. We employ high-accuracy relativistic coupled cluster methods and systematically evaluate the uncertainties of our computational approach. Compared to other Ba- and Yb-containing molecules, BaCH$_3$ and YbCH$_3$ exhibit larger $W_\mathrm{d}$ and $W_\mathrm{s}$ associated to increased covalent character of the M--C bond. The calculated values are $3.22\pm 0.11 \times 10^{24}\frac{h\text{Hz}}{e\text{cm}}$ and $13.80\pm 0.35 \times 10^{24}\frac{h\text{Hz}}{e\text{cm}}$ for $W_\mathrm{d}$, and $8.42\pm0.29$~$h$kHz and $45.35\pm1.15$~$h$kHz for $W_\mathrm{s}$, in BaCH$_3$ and YbCH$_3$, respectively. The robust, accurate, and cost-effective computational scheme reported in this work makes our results suitable for extracting the relevant fundamental properties from future measurements and also can be used to explore other polyatomic molecules sensitive to various violations of fundamental symmetries.

Theoretical calculations based on the Dirac--Coulomb--Breit relativistic coupled cluster method have been carried out for the electron affinities and ionization potentials of the superheavy element nihonium (Nh) and its lighter homologues In and Tl. The In and Tl calculations are in agreement with measurement within uncertainties. For Nh, where experiment is yet unknown, we predict the ionization potential of 7.569(48) eV and electron affinity of 0.776(30) eV.

We present high accuracy relativistic coupled cluster calculations of the first and second ionisation potentials and the electron affinity of the heaviest element in the Periodic Table, Og. The results were extrapolated to the basis set limit and augmented with the higher order excitations (up to perturbative quadruples), the Breit contribution, and the QED self energy and vacuum polarisation corrections. We have performed an extensive investigation of the effect of the various computational parameters on the calculated properties, which allowed us to assign realistic uncertainties on our predictions. Similar study on the lighter homologue of Og, Rn, yields excellent agreement with experiment for the first ionisation potential and a reliable prediction for the second ionisation potential.

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-\textit{e}EDM collaboration is preparing an experiment to measure the eEDM and S-PS interaction in a slow beam of cold BaF molecules [Eur. Phys. J. D, 72, 197 (2018)]. Accurate knowledge of the electronic structure parameters, $W_d$ and $W_s$, 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 $W_d$ and $W_s$ 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 $\pm$ $0.12 \times 10^{24}\frac{\text{Hz}}{e\cdot \text{cm}}$ for $W_d$ and 8.29 $\pm$ 0.12 kHz for $W_s$.

Measurements of nuclear spin-dependent parity-violating (NSD-PV) effects provide an excellent opportunity to test nuclear models and to search for physics beyond the Standard Model. Molecules possess closely-spaced states with opposite parity which may be easily tuned to degeneracy to greatly enhance the observed parity-violating effects. A high-sensitivity measurement of NSD-PV effects using light triatomic molecules is in preparation [E. B. Norrgard, et al., Commun. Phys. 2, 77 (2019)]. Importantly, by comparing these measurements in light nuclei with prior and ongoing measurements in heavier systems, the contribution to NSD-PV from $Z^0$-boson exchange between the electrons and the nuclei may be separated from the contribution of the nuclear anapole moment. Furthermore, light triatomic molecules offer the possibility to search for new particles, such as the postulated $Z^{\prime}$ boson. In this work, we detail a sensitive measurement scheme and present high-accuracy molecular and nuclear calculations needed for interpretation of NSD-PV experiments on triatomic molecules composed of light elements Be, Mg, N, and C. The ab initio nuclear structure calculations, performed within the No-Core Shell Model (NCSM) provide a reliable prediction of the magnitude of different contributions to the NSD-PV effects in the four nuclei. These results differ significantly from the predictions of the standard single-particle model and highlight the importance of including many-body effects in such calculations. In order to extract the NSD-PV contributions from measurements, a parity-violating interaction parameter $W_{\text{PV}}$, which depends on molecular structure, needs to be known with high accuracy. We have calculated these parameters for the triatomic molecules of interest using the relativistic coupled-cluster approach., Comment: 16 pages, 4 figures, and 6 tables

The four-component relativistic Fock space coupled cluster method is used to describe the magnetic hyperfine interaction in low-lying electronic states of the KCs molecule. Both diagonal and off-diagonal matrix elements as functions of the internuclear separation $R$ are calculated within the finite-field scheme. The resulting matrix elements exhibit very weak dependence on $R$ for the separations exceeding 8 \AA , whereas in the vicinity of the ground-state equilibrium the deviation of molecular HFS matrix elements from the atomic values reaches 15\%. The dependence of the computed HFS couplings on the level of core correlation treatment is discussed., Comment: 3 figures

Modern challenges arising in the fields of theoretical and experimental physics require new powerful tools for high-precision electronic structure modelling; one of the most perspective tools is the relativistic Fock space coupled cluster method (FS-RCC). Here we present a new extensible implementation of the FS-RCC method designed for modern parallel computers. The underlying theoretical model, algorithms and data structures are discussed. The performance and scaling features of the implementation are analyzed. The software developed allows to achieve a completely new level of accuracy for prediction of properties of atoms and molecules containing heavy and superheavy nuclei.

One of the most important properties influencing the chemical behavior of an element is the energy released with the addition of an extra electron to the neutral atom, referred to as the electron affinity (EA). Among the remaining elements with unknown EA is astatine, the purely radioactive element 85. Astatine is the heaviest naturally occurring halogen and its isotope $^{211}$At is remarkably well suited for targeted radionuclide therapy of cancer. With the At$^-$ anion being involved in many aspects of current astatine labelling protocols, the knowledge of the electron affinity of this element is of prime importance. In addition, the EA can be used to deduce other concepts such as the electronegativity, thereby further improving the understanding of astatine's chemistry. Here, we report the first measurement of the EA for astatine to be 2.41578(7)eV. This result is compared to state-of-the-art relativistic quantum mechanical calculations, which require incorporation of the electron-electron correlation effects on the highest possible level. 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, which are produced at a one-atom-at-a-time rate.

DIRAC is a freely distributed general-purpose program system for 1-, 2- and 4-component relativistic molecular calculations at the level of Hartree--Fock, Kohn--Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, coupled cluster and electron propagator 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 the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding, and frozen density embedding models. DIRAC was one of the earliest codes for relativistic molecular calculations and remains a reference in its field., Comment: The following article has been submitted to The Journal of Chemical Physics. After it is published, it will be found at https://aip.scitation.org/toc/jcp/current

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 withing the finite-field scheme. The two selected test systems are $^{133}$Cs and $^{137}$BaF. 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 $\sim$ 5.5\%, while the actual deviation of our results from experimental values was $<1$\% in all cases.

Nuclear magnetic quadrupole moments (MQMs), like intrinsic electric dipole moments of elementary particles, violate both parity and time-reversal symmetry and therefore probe physics beyond the Standard Model of particle physics. 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.

Polyatomic polar molecules are promising systems for future experiments that search for violation of time-reversal and parity symmetries due to their advantageous electronic and vibrational structure, which allows laser cooling, full polarisation of the molecule, and reduction of systematic effects [I. Kozyryev and N.R. Hutzler, Phys, Rev. Lett. {\bf 119}, 133002 (2017)]. In this work we investigate the enhancement factor of the electric dipole moment of the electron ($E_\text{eff}$) in the triatomic monohydroxide molecules BaOH and YbOH within the high-accuracy relativistic coupled cluster method. The recommended $E_\text{eff}$ values of the two systems are 6.65 $\pm$ 0.15 GV/cm and 23.4 $\pm$ 1.0 GV/cm, respectively. We compare our results with similar calculations for the isoelectronic diatomic molecules BaF and YbF, which are currently used in experimental search for $P,T$-odd effects in molecules. The $E_\text{eff}$ values prove to be very close, within about 1.5 $\%$ difference in magnitude between the diatomic and the triatomic compounds. Thus, BaOH and YbOH have a similar enhancement of the electron electric dipole moment, while benefiting from experimental advantages, and can serve as excellent candidates for next-generation experiments.

We present high accuracy relativistic coupled cluster calculations of the P-odd interaction coefficient $W_A$ describing the nuclear anapole moment effect on the molecular electronic structure. The molecule under study, BaF, is considered a promising candidate for the measurement of the nuclear anapole moment, and the preparation for the experiment is now underway [Altunas et al., Phys. Rev. Lett. 120, 142501 (2018)]. Influence of various computational parameters (size of the basis set, treatment of relativistic effects, and treatment of electron correlation) on the calculated $W_A$ coefficient is investigated and a recommended value of 147.7 Hz with an estimated uncertainty of 1.5% is proposed., Comment: 9 pages, 3 figures, 3 tables, minor changes, the published version