Enhanced sensitivity to ultralight bosonic dark matter in the spectra of the linear radical SrOH Ivan Kozyryev,1, 2, ∗ Zack Lasner,1, 2,† and John M. Doyle1, 2 1Harvard-MIT Center for Ultracold Atoms, Cambridge, MA 02138 2Department of Physics, Harvard University, Cambridge, MA 02138 (Dated: January 22, 2021) Coupling between Standard Model particles and theoretically well-motivated ultralight dark matter (UDM) candidates can lead to time variation of fundamental constants, including the proton-to-electron mass ratio µ ≡ mp=me ≈ 1836. The presence of nearly-degenerate vibrational energy levels of different character in poly- atomic molecules can result in significantly enhanced relative energy shifts in molecular spectra originating from ¶t µ, relaxing experimental complexity required for high-sensitivity measurements. We analyze the am- plification of UDM effects in the spectrum of laser-cooled strontium monohydroxide (SrOH). SrOH was the first polyatomic molecule to be directly laser cooled to sub-millikelvin temperatures [Kozyryev et al., Phys. Rev. Lett. 118, 173201 (2017)], opening the possibility of long experimental coherence times and providing a 3 promising platform for suppressing systematic errors. Because of the high enhancement factors ( Qµ ≈ 10 ), measurements of the X˜ (200) $ X˜ 0310 rovibrational transitions of SrOH in the microwave regime can re- sult in ∼ 10−17 fractional uncertainty in dµ=µ with one day of integration, leading to significantly improved constraints for UDM coupling constants. We also detail how the use of more complex MOR-type radicals with additional vibrational modes arising from larger ligands R could lead to even greater enhancement factors, while still being susceptible to direct laser cooling. I. INTRODUCTION sensitivity [26–28]. Additionally, specific atomic transitions with enhanced sensitivities jQaj 1 have allowed measure- The quantum mechanical nature of dark matter remains a ments on a Dy beam to be competitive with atomic clock mystery despite significant experimental efforts [1–5]. Strin- limits [13, 29]. Exploring both Gµ and Ga is important as gent limits placed recently on the promising class of dark mat- these effects probe different underlying physical phenomena ter candidates, Weakly Interacting Massive Particles [3,6], as [18]. While the use of atomic clocks for probing dark-matter- well as the absence of signatures for supersymmetric partners induced oscillating, drifting and transient-in-time fundamen- at the Large Hadron Collider [1,7,8] and electron electric tal constants has been considered in depth [21, 30, 31], laser- dipole moment (EDM) experiments [9, 10], have motivated a cooled molecules have additional degrees of freedom that new generation of searches for other theoretically motivated could enable further breakthroughs in this area. dark matter candidates [11–17]. Bosonic ultralight dark mat- In molecular spectra, the energy scales for electronic, vibra- −1=2 ter (UDM) particles, like axions, axion-like particles (ALPs), tional and rotational transitions typically relate as 1 : µ : −1 dilatons, moduli, and relaxions [11], can form coherently os- µ [32]. Molecular transitions provide a system to study Gµ cillating classical fields f(r; t) = f cosw t − k ·r with the couplings without any contributions from Ga because vibra- 0 f f 1 oscillation frequency set by the mass of the dark matter par- tional transitions in molecules have Qµ = − 2 and Qa ≈ 0 µ ticle wf ' mf [18–20]. Coupling between UDM fields and [23]. Thus, isolating effects from variation in a model- ordinary matter can lead to variation in fundamental constants independent manner becomes possible [33]. Moreover, cer- X = a (fine-structure constant) and µ (proton-to-electron mass tain beyond the Standard Model theories predict larger µ vari- ratio) as [19, 21] ation dµ=µ = Rda=a with R ≈ 40 [23, 34, 35], further motivat- ing precision experiments in molecular spectroscopy. Molec- dX (t) n = GX f (r;t); (1) ular ions can also be used for such experiments and recent the- X oretical proposals consider using diatomic hetero- and homo- where the coupling strength is GX and n = 1(2) for lin- nuclear molecular ions to search for µ variation [36, 37]. In ear(quadratic) coupling. this paper, we propose to use laser-cooled samples of the neu- Transitions between different quantum levels with energy tral polyatomic radical SrOH that can be trapped at high den- separation 4E = ~w in atoms and molecules are dependent sities and low temperatures, allowing for large scalability and on the dimensionless constants dw = f (da;dµ) with [22, 23] enhanced sensitivity to UDM-induced µ variation. arXiv:1805.08185v2 [physics.atom-ph] 21 Jan 2021 dw da dµ = Qa + Qµ : (2) w a µ II. ENHANCED SENSITIVITY TO DARK MATTER WITH Sensitive probes of a variation due to UDM-induced ef- NEAR-DEGENERATE STATES fects have recently been explored with the use of ultrapre- −17 cise atomic clocks [18, 24, 25], reaching ¶t a=a ∼ 10 /yr As previously pointed out [22, 38–41], rovibrational spec- troscopy of diatomic and polyatomic molecules may provide significant enhancements in relative sensitivity to the variation ∗ [email protected] of µ with Qµ 1. An extensive list of enhancement factors † [email protected] calculated for diatomic and polyatomic molecules to µ varia- 2 To see how a large sensitivity to µ variation arises in the rovibrational spectrum of SrOH, we begin by following pre- vious treatments in Refs. [22, 23]. Consider two different energy levels Eg and Ee within the same electronic state with Eg < Ee. The energy difference is (assuming ~ = 1) w = Ee − Eg (3) (0330) with the change arising from µ variation given as (200) ¶E ¶E (0310) dw = e dµ − g dµ: (4) (0220) ¶µ ¶µ (100) (0200) Therefore, the fractional change in the level separation is (0110) dw 1 ¶Ee ¶Eg -1 = − dµ: (5) 527 cm w Ee − Eg ¶µ ¶µ 364 cm-1 ~ Equivalently, the relationship between the fractional changes X2Σ+ (000) in w and µ are related to each other as dw dµ Figure 1. Vibrational levels of SrOH in the electronic ground state. = Qµ (6) The dominant nature of the vibrational motion is indicated schemat- w µ l ically at the top. The states are labeled using v1v2v3 notation for with the proportionality constant Qµ also known as the dimen- the number of vibrational quanta in the Sr-O stretching (v1), Sr-O-H sionless enhancement factor defined as bending (v2) and O-H stretching (v3) vibrational motion. The num- ber of units of vibrational angular momentum present in the doubly- µ ¶Ee ¶Eg Qµ ≡ − (7) degenerate bending of linear triatomic molecules is denoted with the Ee − Eg ¶µ ¶µ superscript l. Nearly-degenerate excited stretching (200) and bend- or as more common in the literature ing 0310 vibrations are indicated with a red oval. 1 ¶E ¶E Q ≡ e − g : (8) µ w ¶(lnµ) ¶(lnµ) tion can be found in Refs. [22, 23], with large enhancements The absolute dependence of each energy level is calculated as of Q ∼ 300 and Q ∼ 700 estimated for CH OH and l-C H, µ µ 3 3 ¶Eg;e respectively. Other polyatomic molecules found in space like qg;e ≡ (9) ¶(lnµ) methanol [42], acetone [43] and ammonia [44] have been ana- −7 lyzed as well, leading to a stringent limit of dµ=µ . 10 from and has units of energy (usually cm-1). From Eq.8 one can the observations of astronomical methanol [23, 45]. While as- observe that a large enhancement factor Qµ will arise when trophysical observations place stringent time-variation limits two levels being probed are closely spaced (i.e. w ≈ 0) and −17 with ¶t µ=µ ∼ 10 /yr bounds due to large look-back times have different dependence on µ (i.e. qg 6= qe). (4t ∼ 7Gyr) [23], they have limited sensitivity to UDM- Generically, the interplay between harmonic and anhar- induced coherent oscillations since a linear drift, ¶t µ = dµ=4t, monic contributions (discussed in detail for SrOH below) must be assumed. to the difference in sensitivity coefficients, 4q, can lead to Here we analyze the enhancement factors for one of the enhancement factors Qµ significantly larger than unity. In simplest possible polyatomic molecules, the linear triatomic order to demonstrate the role of both harmonic and anhar- XYZ-type radical SrOH, and discover that enhancement fac- monic terms, we consider two vibrational levels separated 3 tors of Qµ ≈ 10−10 can be reached by probing rovibrational by 4E = 4Eharm + 4Eanharm. Using the dependence of vi- ˜ ˜ 1 brational constants on the proton-to-electron mass ratio (see transitions of the X (200) $ X 03 0 excitation spectrum in 1 the w ≈ 2p × 1 − 30 GHz transition frequency band. SrOH App.C) we calculate 4q = − 2 4Eharm − 4Eanharm for the µ was the first polyatomic molecule to be directly laser cooled sensitivity difference. Therefore, the absolute enhancement [46], and thus provides the additional significant advantages factor Qµ = 4q=4E becomes of low translational and internal temperatures, long experi- 1 DE Q = − 1 + anharm : (10) mental coherence times, and optical internal state prepara- µ 2 4E + 4E tion and efficient readout. Furthermore, the simple vibrational harm anharm structure of SrOH strongly limits the possibility of internal vi- For illustration, we consider three limiting cases, depending brational redistribution (IVR) or nonradiative transitions [47], on the relative contributions of 4Eharm and 4Eanharm: enabling highly sensitive laboratory measurements of both 1 4E = 0 : Q ! − dµ=µ and ¶t µ=µ in the frequency band of theoretical interest anharm µ 2 for promising UDM models. Figure1 shows the relevant vi- 4E = 0 : Q ! −1 (11) brational energy levels of SrOH in the ground electronic state harm µ X˜ . 4Eharm ∼ −4Eanharm : Qµ ! ±¥: 3 2 Therefore, a large enhancement factor is expected for a transi- −giili − B[v]N (N + 1); tion with anharmonic contributions comparable in magnitude to the harmonic oscillator energy difference and opposite in where l1(l2) = 0(1), d1(d2) = 1(2) and the sensitivity q[0;000] sign.