Search

We propose diatomic molecules built from gold and carbon-group atoms as promising candidates for optical cycling and precision measurements. We show that this class of molecules (AuX, X = C, Si, Ge, Sn, Pb) features laser-accessible electronic transitions with nearly diagonal Franck-Condon factors. The $^2\Pi_{1/2}$ ground states can be easily polarized in the laboratory frame and have near-zero magnetic moments, valuable features for quantum science and precision measurement applications. The sensitivities of AuX molecules to the electron electric dipole moment (EDM) are found to be favorable, with effective electric fields of 10-30 GV/cm. Together, these features imply that AuX molecules may enable significantly improved searches for time-reversal symmetry violation., Comment: 10 pages, 3 figures

Relativistic exact two-component coupled-cluster calculations of molecular sensitivity factors for nuclear Schiff moments (NSMs) are reported. We focus on molecules containing heavy nuclei, especially octupole-deformed nuclei. Analytic relativistic coupled-cluster gradient techniques are used and serve as useful tools for identifying candidate molecules that sensitively probe for physics beyond the Standard Model in the hadronic sector. Notably, these tools enable straightforward ``black-box'' calculations. Two competing chemical mechanisms that contribute to the NSM are analyzed, illuminating the physics of ligand effects on NSM sensitivity factors.

The vibrational branching ratios from the lowest excited electronic state for $\textrm{SrOCH}_3$, $\textrm{SrNH}_2$, and $\textrm{SrSH}$ are measured at the $< 0.1\%$ level. Spectra are obtained by driving the $\tilde{X} - \tilde{A}$ transitions and dispersing the fluorescence on a grating spectrometer. We also perform $\textit{ab initio}$ calculations for the energies of vibrational levels relevant for laser cooling, as well as branching ratios to support the interpretations of all molecular spectra. Symmetry group analysis is applied in conjunction with our data to study rotational closure in these molecules. These analyses indicate favorable prospects for laser cooling $\textrm{SrNH}_2$ and other similar alkaline-earth(-like) amides for future beyond the Standard Model physics searches using polyatomic molecules with long-lived parity doublets.

A new computational framework for spinor-based relativistic exact two-component (X2C) calculations is developed using contracted basis sets with a spin-orbit contraction scheme. Generally contracted j-adapted basis sets using primitive functions in the correlation-consistent basis sets are constructed for the X2C Hamiltonian with atomic mean-field spin-orbit integrals (the X2CAMF scheme). The contraction coefficients are taken from atomic X2CAMF Hartree-Fock spinors, hereby following the simple concept of linear combination of atomic orbitals (LCAOs). Benchmark calculations of spin-orbit splittings, equilibrium bond lengths, and harmonic vibrational frequencies demonstrate the accuracy and efficacy of the j-adapted spin-orbit contraction scheme.

Combinatorial optimization (CO) on graphs is a classic topic that has been extensively studied across many scientific and industrial fields. Recently, solving CO problems on graphs through learning methods has attracted great attention. Advanced deep learning methods, e.g., graph neural networks (GNNs), have been used to effectively assist the process of solving COs. However, current frameworks based on GNNs are mainly designed for certain CO problems, thereby failing to consider their transferable and generalizable abilities among different COs on graphs. Moreover, simply using original graphs to model COs only captures the direct correlations among objects, which does not consider the mathematical logicality and properties of COs. In this paper, we propose a unified pre-training and adaptation framework for COs on graphs with the help of the maximum satisfiability (Max-SAT) problem. We first use Max-SAT to bridge different COs on graphs since they can be converted to Max-SAT problems represented by standard formulas and clauses with logical information. Then, we further design a pre-training and domain adaptation framework to extract the transferable and generalizable features so that different COs can benefit from them. In the pre-training stage, Max-SAT instances are generated to initialize the parameters of the model. In the fine-tuning stage, instances from CO and Max-SAT problems are used for adaptation so that the transferable ability can be further improved. Numerical experiments on several datasets show that features extracted by our framework exhibit superior transferability and Max-SAT can boost the ability to solve COs on graphs.

A Cholesky decomposition (CD)-based implementation of relativistic two-component coupled-cluster (CC) and equation-of-motion CC (EOM-CC) methods using an exact two-component Hamiltonian augmented with atomic-mean-field spin-orbit integrals (the X2CAMF scheme) is reported. The present CD-based implementation of X2CAMF-CC and EOM-CC methods employs atomic-orbital-based algorithms to avoid the construction of two-electron integrals and intermediates involving three and four virtual indices. Our CD-based implementation extends the applicability of X2CAMF-CC and EOM-CC methods to medium-sized molecules with the possibility to correlate around 1000 spinors. Benchmark calculations for uranium-containing small molecules have been performed to assess the dependence of the CC results on the Cholesky threshold. A Cholesky threshold of $10^{-4}$ is shown to be sufficient to maintain chemical accuracy. Example calculations to illustrate the capability of the CD-based relativistic CC methods are reported for the bond-dissociation energy of the uranium hexafluoride molecule, UF$_6$, with up to quadruple-zeta basis sets, and the lowest excitation energy in solvated uranyl ion [UO$_2^{2+}$(H$_2$O)$_{12}$].

A study of the intensity-borrowing mechanisms important to optical cycling transitions in laser-coolable polyatomic molecules arising from non-adiabatic coupling, contributions beyond the Franck-Condon approximation, and Fermi resonances is reported. It has been shown to be necessary to include non-adiabatic coupling to obtain computational accuracy that is sufficient to be useful for laser cooling of molecules. The predicted vibronic branching ratios using perturbation theory based on the non-adiabatic mechanisms have been demonstrated to agree well with those obtained from variational discrete variable representation calculations for representative molecules including CaOH, SrOH, and YbOH. The electron-correlation and basis-set effects on the calculated transition properties, including the vibronic coupling constants, the spin-orbit coupling matrix elements, and the transition dipole moments, and on the calculated branching ratios have been thoroughly studied. The vibronic branching ratios predicted using the present methodologies demonstrate that RaOH is a promising radioactive molecule candidate for laser cooling.

Molecular predissociation is the spontaneous, nonradiative bond breaking process that can occur upon excitation. In the context of laser cooling, predissociation is an unwanted consequence of molecular structure that limits the ability to scatter a large number of photons required to reach the ultracold regime. Unlike rovibrational branching, predissociation is irreversible since the fragments fly apart with high kinetic energy. Of particular interest is the simple diatomic molecule, CaH, for which the two lowest electronically excited states used in laser cooling lie above the dissociation threshold of the ground potential. In this work, we present measurements and calculations that quantify the predissociation probabilities affecting the cooling cycle. The results allow us to design a laser cooling scheme that will enable the creation of an ultracold and optically trapped cloud of CaH molecules. In addition, we use the results to propose a two-photon pathway to controlled dissociation of the molecules, in order to gain access to their ultracold fragments, including hydrogen., Comment: 16 pages, 4 figures

We study the internal state dynamics of optically trapped polyatomic molecules subject to room temperature blackbody radiation. Using rate equations that account for radiative decay and blackbody excitation between rovibrational levels of the electronic ground state, we model the microscopic behavior of the molecules' thermalization with their environment. As an application of the model, we describe in detail the procedure used to determine the blackbody and radiative lifetimes of low-lying vibrational states in ultracold CaOH molecules, the values of which were reported in previous work [Hallas et al., arXiv:2208.13762]. Ab initio calculations are performed and are found to agree with the measured values. Vibrational state lifetimes for several other laser-coolable molecules, including SrOH and YbOH, are also calculated., Comment: arXiv admin note: substantial text overlap with arXiv:2208.13762

Characterization of the inner-shell decay processes in molecules containing heavy elements is key to understanding x-ray damage of molecules and materials and for medical applications with Auger-electron-emitting radionuclides. The 1s hole states of heavy atoms can be produced by absorption of tunable x-rays and the resulting vacancy decays characterized by recording emitted photons, electrons, and ions. The 1s hole states in heavy elements have large x-ray fluorescence yields that transfer the hole to intermediate electron shells that then decay by sequential Auger-electron transitions that increase the ion's charge state until the final state is reached. In molecules the charge is spread across the atomic sites, resulting in dissociation to energetic atomic ions. We have used x-ray/ion coincidence spectroscopy to measure charge states and energies of I$^{q+}$ and Br$^{q'+}$ atomic ions following 1s ionization at the I and Br \textit{K}-edges of IBr. We present the charge states and kinetic energies of the two correlated fragment ions associated with core-excited states produced during the various steps of the cascades. To understand the dynamics leading to the ion data, we develop a computational model that combines Monte-Carlo/Molecular Dynamics simulations with a classical over-the-barrier model to track inner-shell cascades and redistribution of electrons in valence orbitals and nuclear motion of fragments.

We report optical trapping of a polyatomic molecule, calcium monohydroxide (CaOH). CaOH molecules from a magneto-optical trap are sub-Doppler laser cooled to $20(3)~\mu\text{K}$ in free space and loaded into an optical dipole trap. We attain an in-trap molecule number density of $3(1) \times 10^9~\text{cm}^{-3}$ at a temperature of $57(8)~\mu$K. Trapped CaOH molecules are optically pumped into an excited vibrational bending mode, whose $\ell$-type parity doublet structure is a potential resource for a wide range of proposed quantum science applications with polyatomic molecules. We measure the spontaneous, radiative lifetime of this bending mode state to be $\sim$$0.7~\text{s}$.

This work discusses a protocol for constructing highly accurate potential energy curves (PECs) for the lowest two states of Rb$_{2}^+$, i.e. $X\,{}^2{\Sigma}{_g^+}$ and $(1) {}^2\Sigma{_u^+}$, using an additivity scheme based on coupled-cluster theory. The approach exploits the findings of our previous work [J. Schnabel, L. Cheng and A. K\"ohn, J. Chem. Phys. 155, 124101 (2021)] to avoid the unphysical repulsive long-range barrier occurring for symmetric molecular ions when perturbative estimates of higher-order cluster operators are employed. Furthermore, care was taken to reproduce the physically correct exchange splitting of the $X {}^2{\Sigma}{_g^+}$ and $(1) {}^2{\Sigma}{_u^+}$ PECs. The accuracy of our computational approach is benchmarked for ionization energies of Rb and for spectroscopic constants as well as vibrational levels of the $a {}^3{\Sigma}{_u^+}$ triplet state of Rb\textsubscript{2}. We study high-level correlation contributions, high-level relativistic effects and inner-shell correlation contributions and find very good agreement with experimental reference values for the atomic ionization potential and the binding energy of Rb$_{2}$ in the $a\,{}^3{\Sigma}{_u^+}$ triplet state. Our final best estimate for the binding energy of the Rb$_{2}^+$ $X {}^2{\Sigma}{_g^+}$ state including zero-point vibrational contributions is $D_0 = 6179\,\mathrm{cm}^{-1}$ with an estimated error bound of $\mathcal{O}(\pm 30\,\mathrm{cm}^{-1})$. This value is smaller than the experimentally inferred lower bond of $D_0\ge 6307.5\,\mathrm{cm}^{-1}$ [Bellos et al., Phys. Rev. A 87, 012508 (2013)] and will require further investigation. For the $(1) {}^2{\Sigma}{_u^+}$ state a shallow potential with $D_0 = 78.4\,\mathrm{cm}^{-1}$ and an error bound of $\pm 9\,\mathrm{cm}^{-1}$ is computed.

The vibrational branching ratios of SrOH for radiative decay to the ground electronic state, $X^{2}\Sigma^{+}$, from the first two electronically excited states, $A^{2}\Pi$ and $B^{2}\Sigma^{+}$, are determined experimentally at the $\sim10^{-5}$ level. The observed small branching ratios enable the design of a full, practical laser-cooling scheme, including magneto-optical trapping and sub-Doppler laser cooling, with $>10^4$ photon scatters per molecule. Ab initio calculations sensitive to weak vibronic transitions are performed to facilitate the experimental measurement and analysis, and show good agreement with experiment., Comment: 6 pages, 3 figures, 3 tables

Static electric dipole moments of nondegenerate systems probe mass scales for physics beyond the Standard Model well beyond those reached directly at high energy colliders. Discrimination between different physics models, however, requires complementary searches in atomic-molecular-and-optical, nuclear and particle physics. In this report, we discuss the current status and prospects in the near future for a compelling suite of such experiments, along with developments needed in the encompassing theoretical framework., Comment: Contribution to Snowmass 2021; updated with community edits and endorsements

An excellent candidate molecule for the measurement of the electron's electric dipole moment (eEDM) is thorium monofluoride (ThF$^+$) because the eEDM-sensitive state, $^3\Delta_1$, is the electronic ground state, and thus is immune to decoherence from spontaneous decay. We perform spectroscopy on $X\,^3\Delta_1$ to extract three spectroscopic constants crucial to the eEDM experiment: the hyperfine coupling constant, the molecular frame electric dipole moment, and the magnetic $g$-factor. To understand the impact of thermal blackbody radiation on the vibrational ground state, we study the lifetime of the first excited vibrational manifold of $X\,^3\Delta_1$. We perform ab initio calculations, compare them to our results, and discuss prospects for using ThF$^+$ in a new eEDM experiment at JILA., Comment: 9 pages, 8 figures

The YbF molecule is a sensitive system for measuring the electron's electric dipole moment. The precision of this measurement can be improved by direct laser cooling of the molecules to ultracold temperature. However, low-lying electronic states arising from excitation of a 4f electron may hinder laser cooling. One set of these "4f hole" states lies below the $A^2\Pi_{1/2}$ excited state used for laser cooling, and radiative decay to these intermediate levels, even with branching ratios as small as $10^{-5}$, can be a hindrance. Other 4f hole states lie very close to the $A^2\Pi_{1/2}$ state, and a perturbation results in states of mixed character that are involved in the laser cooling cycle. This perturbation may enhance the loss of molecules to states outside of the laser cooling cycle. We model the perturbation of the $A^2\Pi_{1/2}$ state to determine the strength of the coupling between the states, the de-perturbed potential energy curves, and the radiative branching ratios to various vibrational levels of the ground state, $X ^{2}\Sigma^+$. We use electronic structure calculations to characterise the 4f hole states and the strengths of transitions between these states and the $A^2\Pi_{1/2}$ and $X ^{2}\Sigma^+$ states. We identify a leak out of the cooling cycle with a branching ratio of roughly $5 \times 10^{-4}$, dominated by the contribution of the ground state configuration in a 4f hole state. Finally, we assess the impact of these results for laser cooling of YbF and molecules with similar structure., Comment: 11 pages, 3 figures, minor changes from previous version

The development of a quadratic unitary coupled-cluster singles and doubles (qUCCSD) based self-consistent polarization propagator method is reported. We present a simple strategy for truncating the commutator expansion of the UCC transformed Hamiltonian $\bar{H}$. The qUCCSD method for the electronic ground-state includes up to double commutators for the amplitude equations and up to cubic commutators for the energy expression. The qUCCSD excited-state eigenvalue equations include up to double commutators for the singles-singles block of $\bar{H}$, single commutators for the singles-doubles and doubles-singles blocks, and the bare Hamiltonian for the doubles-doubles block. Benchmark qUCCSD calculations of the ground-state properties and excitation energies for representative molecules demonstrate significant improvement of the accuracy and robustness over the previous UCC3 scheme derived using M{\o}ller-Plesset perturbation theory.

We reveal limitations of several standard coupled-cluster (CC) methods with perturbation-theorybased noniterative or approximate iterative treatments of triple excitations when applied to thedetermination of highly accurate potential energy curves (PECs) of ionic dimers, such as the X $^2\Sigma^+_g$ electronic ground state of Rb$_2^+$. Such computations are of current interest for the understanding ofion-atom interactions in the ultracold regime. We demonstrate that these coupled-cluster methodslead to an unphysical long-range barrier for the Rb$_2^+$ system. The barrier is small but clearly spoilsthe entire long-range behavior of the PEC. The effect is also found for other X$_2^+$ systems, like X = Li, Na, and K. Calculations using a new flexible framework for obtaining leading perturbativetriples corrections derived using an analytic CC singles and doubles (CCSD) energy derivativeformulation demonstrate that the origin of this problem lies in the use of $T_3$ amplitudes obtained from approximate CC singles, doubles and triples (CCSDT) amplitude equations. It is shown thatthe unphysical barrier is related to a symmetry instability of the underlying Hartree-Fock mean-fieldsolution. Physically meaningful perturbative corrections in the long-range tail of the PEC mayinstead be obtained using symmetry-broken reference determinants.

We report a generally applicable computational and experimental approach to determine vibronic branching ratios in linear polyatomic molecules to the $10^{-5}$ level, including for nominally symmetry forbidden transitions. These methods are demonstrated in CaOH and YbOH, showing approximately two orders of magnitude improved sensitivity compared with the previous state of the art. Knowledge of branching ratios at this level is needed for the successful deep laser cooling of a broad range of molecular species.

We develop an analytic-gradient based method for relativistic coupled-cluster calculations of effective electric field, $\mathcal{E}_{\text{eff}}$, with improved efficiency and robustness over the previous state of the art. The enhanced capability to calculate this time-reversal symmetry violation sensitivity parameter enables efficient screening of candidate molecules for the electron electric dipole moment (eEDM) search. As examples, the |$\mathcal{E}_{\text{eff}}$| values of metal methoxides including BaOCH$_3$, YbOCH$_3$, and RaOCH$_3$ are shown to be as large as those of the corresponding fluorides and hydroxides, which supports the recent proposal of using these symmetric-top molecules to improve the sensitivity of eEDM measurements. The computational results also show that molecules containing late actinide elements, NoF, NoOH, LrO, and LrOH$^+$, exhibit particularly large |$\mathcal{E}_{\text{eff}}$| values of around 200 GV/cm.

Medium resolution ($\Delta \nu$ ~ 3 GHz) laser-induced fluorescence (LIF) excitation spectra of a rotationally cold sample of YbOH in the 17300-17950 cm$^{-1}$ range have been recorded using two-dimensional (excitation and dispersed fluorescence) spectroscopy. High resolution ($\Delta \lambda$ ~ 0.65 nm) dispersed laser induced fluorescence (DLIF) spectra and radiative decay curves of numerous bands detected in the medium resolution LIF excitation spectra were recorded. The vibronic energy levels of the $\tilde{X} \, ^2\Sigma^+$ state were predicted using a discrete variable representation approach and compared with observations. The radiative decay curves were analyzed to produce fluorescence lifetimes. DLIF spectra resulting from high resolution ($\Delta \nu$ < 10 MHz) LIF excitation of individual low-rotational lines in the $\tilde{A} \, ^2\Pi_{1/2}(0,0,0) - \tilde{X} \, ^2\Sigma^+(0,0,0)$, $\tilde{A} \, ^2\Pi_{1/2}(1,0,0) - \tilde{X} \, ^2\Sigma^+(0,0,0)$, $[17.73]\Omega=0.5(0,0,0) - \tilde{X} \, ^2\Sigma^+(0,0,0)$ bands were also recorded. The DLIF spectra were analyzed to determine branching ratios which were combined with radiative lifetimes to obtain transition dipole moments. The implications for laser cooling and trapping of YbOH are discussed., Comment: 62 pages, 14 figures

An extensive analysis has been carried out of the performance of standard families of basis sets with the hierarchy of coupled cluster methods CC2, CCSD, CC3 and CCSDT in computing selected Oxygen, Carbon and Nitrogen K-edge (vertical) core excitation and ionization energies within a core-valence separated scheme in the molecules water, ammonia, and carbon monoxide. Complete basis set limits for the excitation energies have been estimated via different basis set extrapolation schemes. The importance of scalar relativistic effects has been established within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e)., Comment: Manuscript submitted to Adv. Quantum Chemistry

The molecular ion ThF$^+$ is the species to be used in the next generation of search for the electron's Electric Dipole Moment (eEDM) at JILA. The measurement requires creating molecular ions in the eEDM sensitive state, the rovibronic ground state $^3\Delta_1$, $v^+=0$, $J^+=1$. Survey spectroscopy of neutral ThF is required to identify an appropriate intermediate state for a Resonance Enhanced Multi-Photon Ionization (REMPI) scheme that will create ions in the required state. We perform broadband survey spectroscopy (from 13000 to 44000~cm$^{-1}$) of ThF using both Laser Induced Fluorescence (LIF) and $1+1'$ REMPI spectroscopy. We observe and assign 345 previously unreported vibronic bands of ThF. We demonstrate 30\% efficiency in the production of ThF$^+$ ions in the eEDM sensitive state using the $\Omega = 3/2$ [32.85] intermediate state. In addition, we propose a method to increase the aforementioned efficiency to $\sim$100\% by using vibrational autoionization via core-nonpenetrating Rydberg states, and discuss theoretical and experimental challenges. Finally, we also report 83 vibronic bands of an impurity species, ThO., Comment: 49 pages, 7 figures

Photoionization cross sections and partial ion yields of Xe and XeF$_2$ from Xe 3d$_{5/2}$, Xe 3d$_{3/2}$, and F 1s subshells in the 660--740 eV range are compared to explore effects of the F ligands. The Xe 3d - $\epsilon$f continuum shape resonances dominate the photoionization cross sections of both the atom and molecule, but prominent resonances appear in the XeF$_2$ cross section due to nominal excitation of Xe 3d and F 1s electrons to the lowest unoccupied molecular orbital (LUMO), a delocalized anti-bonding MO. The subshell ionization thresholds, the LUMO resonance energies and their oscillator strengths are calculated by relativistic coupled-cluster methods. Several charge states and fragment ions are produced from the atom and molecule due to alternative decay pathways from the inner-shell holes. Total and partial ion yields vary in response to the shape resonances and LUMO resonances. Previous calculations and measurements of atomic Xe 3d core-hole decay channels and our calculated results for XeF$_2$ guide interpretations of the molecular ion products. The partial ion yields of XeF$_2$ are dominated by Xe 3d core-hole decays, but distinct ion products are measured at the F 1s - LUMO resonance. Xe 3d core-hole decays from XeF$_2$ produce lower charge states in comparison with atomic Xe, and energetic F ions are produced by Coulomb explosions of the molecular ions. The measurements support a model of molecular core-hole decay that begins with a localized hole, stepwise Auger electron emission spreads charge across neighboring atoms, and the system fragments energetically.

We study the inverse boundary crossing problem for diffusions. Given a diffusion process $X_t$, and a survival distribution $p$ on $[0,\infty)$, we demonstrate that there exists a boundary $b(t)$ such that $p(t)=\mathbb{P}[\tau >t]$, where $\tau$ is the first hitting time of $X_t$ to the boundary $b(t)$. The approach taken is analytic, based on solving a parabolic variational inequality to find $b$. Existence and uniqueness of the solution to this variational inequality were proven in earlier work. In this paper, we demonstrate that the resulting boundary $b$ does indeed have $p$ as its boundary crossing distribution. Since little is known regarding the regularity of $b$ arising from the variational inequality, this requires a detailed study of the problem of computing the boundary crossing distribution of $X_t$ to a rough boundary. Results regarding the formulation of this problem in terms of weak solutions to the corresponding Kolmogorov forward equation are presented., Comment: Published in at http://dx.doi.org/10.1214/10-AAP714 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)