Search

Closed-shell atoms and molecules such as Hg or TlF provide some of the best low-energy tests of hadronic $\mathcal{CP}$-violation which is considered to be a necessary ingredient to explain the observed excess of matter over antimatter in our universe. $\mathcal{CP}$-violation is, however, expected to be strongly enhanced in octupole deformed nuclei such as $^{225}$Ra. Recently, closed-shell radium-containing symmetric-top molecular ions were cooled sympathetically in a Coulomb crystal [M. Fan et al., Phys. Rev. Lett. 126, 023002 (2021)] and shown to be well-suited for precision spectroscopy in the search for fundamental physics [P. Yu and N. R. Hutzler, Phys. Rev. Lett. 126, 023003 (2021)]. In closed-shell molecules hadronic $\mathcal{CP}$-violation contributes to a net electric dipole moment (EDM) that violates parity and time-reversal symmetry ($\mathcal{P,T}$), which is the target of measurements. To interpret experiments, it is indispensable to know the electronic structure enhancement parameters for the various sources of $\mathcal{P,T}$-violation which contribute to the net $\mathcal{P,T}$-odd EDM. In this paper we employ relativistic Hartree--Fock and density functional theory calculations to determine relevant parameters for interpretation of possible EDM measurements in RaOCH$_3^+$, RaSH$^+$, RaCH$_3^+$, RaCN$^+$, and RaNC$^+$ and perform accurate relativistic coupled cluster calculations of the Schiff moment enhancement in RaSH$^+$ to gauge the quality of the density functional theory approach. Finally, we project to bounds on various fundamental $\mathcal{P,T}$-odd parameters that could be achievable from an experiment with RaOCH$_3^+$ in the near future and asses the complementarity of this experiment to experiments with Hg and TlF., Comment: Updated references, corrected Wm values for RaCH3+, RaCN+ and RaNC+ in Table VII and in the supplement

In the last decades, material discovery has been a very active research field driven by the need to find new materials for many different applications. This has also included materials with heavy elements, beyond the stable isotopes of lead, as most actinides exhibit unique properties that make them useful in various applications. Furthermore, new heavy elements beyond actinides, collectively referred to as super-heavy elements (SHEs), have been synthesized, filling previously empty space of Mendeleev periodic table. Their chemical bonding behavior, of academic interest at present, would also benefit of state-of-the-art modeling approaches. In particular, in order to perform first-principles calculations with planewave basis sets, one needs corresponding pseudopotentials. In this work, we present a series of scalar- and fully-relativistic optimized norm-conserving Vanderbilt pseudopotentials (ONCVPs) for thirty-four actinides and super-heavy elements, for three different exchange-correlation functionals (PBE, PBEsol and LDA). The scalar-relativistic version of these ONCVPs is tested by comparing equations of states for crystals, obtained with \textsc{abinit} 9.6, with those obtained by all-electron zeroth-order regular approximation (ZORA) calculations, without spin-orbit coupling, performed with the Amsterdam Modeling Suite \textsc{band} code. $\Delta$-Gauge and $\Delta_1$-Gauge indicators are used to validate these pseudopotentials. This work is a contribution to the PseudoDojo project, in which pseudopotentials for the whole periodic table are developed and systematically tested. The pseudopotential files are available on the PseudoDojo web-interface pseudo-dojo.org in psp8 and UPF2 formats, both suitable for \textsc{abinit}, the latter being also suitable for Quantum ESPRESSO.

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

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.

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