Electric field dependent g factors of RaOCH$_3$ molecule

This paper presents a method for calculating electric-field-dependent $g$-factors in symmetric top molecules and applies it to RaOCH$_3$ to identify $K$-doublet levels with small $g$-factor differences, which is crucial for enhancing electron electric dipole moment (eEDM) experiments via laser cooling.

The Gist

Imagine you are a detective trying to catch a ghost. This ghost is the electron electric dipole moment (eEDM). In the world of physics, finding this "ghost" would prove that our current understanding of the universe (the Standard Model) is incomplete and open the door to new, exotic physics.

To catch this ghost, scientists use tiny, spinning tops called molecules. Specifically, they are looking at a molecule called RaOCH₃ (Radium-Methyl-Oxide). Think of this molecule as a tiny, spinning ice skater.

Here is the story of what Alexander Petrov did in this paper, explained simply:

1. The Problem: The Magnetic "Noise"

To find the electron ghost, scientists spin these molecules and look at how they react to electric and magnetic fields.

• The Goal: They want to see a tiny wobble caused by the electron ghost.
• The Trouble: The Earth's magnetic field (and stray magnets in the lab) creates a much louder "noise" that drowns out the ghost. It's like trying to hear a whisper in a rock concert.

2. The Clever Trick: The "Twin" Strategy

Scientists use a special trick involving K-doublets. Imagine the spinning skater has a "twin" spinning in the exact opposite direction.

• The ghost makes the twins wobble in opposite directions (one left, one right).
• The magnetic noise makes them wobble in the same direction (both left, or both right).

If you measure both twins and subtract one result from the other, the magnetic noise cancels out (Left minus Left = Zero), but the ghost signal doubles (Left minus Right = Big Wobble). This is a brilliant way to silence the noise.

3. The New Challenge: The Twins Aren't Perfect

Here is the catch: The twins are almost identical, but not quite. They have slightly different "magnetic personalities" (called g-factors).

• If the twins react slightly differently to the magnetic field, the "noise" doesn't cancel out perfectly. A little bit of noise leaks through, and the experiment gets messy.
• The size of this leak depends on how strong the electric field is. Scientists need to know exactly how much the twins differ at different electric field strengths to plan their experiment perfectly.

4. What This Paper Did: The "Twin Calculator"

Before this paper, scientists knew how to do this math for simple, two-atom molecules (like a dumbbell). But RaOCH₃ is a symmetric top molecule (it looks more like a spinning top with a handle). No one had done the complex math for this shape yet.

Alexander Petrov wrote a new computer program (a "calculator") to figure out:

1. How the "magnetic personalities" (g-factors) of the RaOCH₃ twins change as you turn up the electric field.
2. Why they differ.

5. The Results: Good News for the Ghost Hunters

The calculations revealed some very encouraging things:

• The Twins are Very Similar: The difference in their magnetic personalities is incredibly tiny (about one-millionth of their total magnetism).
• The Electric Field Helps: As you increase the electric field (which is needed to make the molecules spin nicely), the difference between the twins stays very small and predictable.
• Why? The paper explains that the "rules of the quantum world" (specifically the Pauli Exclusion Principle, which acts like a strict bouncer at a club) prevent certain confusing interactions that usually mess things up.

The Bottom Line

This paper is like a user manual for a new, high-tech tool. It tells experimentalists: "If you use RaOCH₃ molecules and set your electric field to about 500 mV/cm, the magnetic noise will be suppressed almost perfectly."

Because the "leak" of magnetic noise is so small, this molecule is a superior candidate for catching the electron ghost. It promises to make the experiment much more sensitive, potentially helping us discover new laws of physics that have been hiding from us for decades.

Technical Summary

1. Problem Statement

The search for the electron electric dipole moment (eEDM) is a primary method for testing physics beyond the Standard Model. Experiments utilizing symmetric top molecules (like RaOCH$_3$) offer significant advantages, particularly when combined with laser cooling, which can enhance statistical sensitivity by orders of magnitude.

However, a major source of systematic error in these experiments is the control of magnetic fields. In eEDM experiments, the energy splitting between levels with opposite angular momentum projections is measured. While the eEDM contribution to this splitting has opposite signs in the two components of a "doublet" (K-doublet in symmetric tops), the Zeeman (magnetic) contribution has the same sign. By subtracting the splittings of the two doublet components, systematic errors related to magnetic field imperfections can be suppressed.

The efficiency of this suppression depends on the ratio $\Delta g / g$, where $\Delta g$ is the difference in magnetic $g$-factors between the two K-doublet components. A smaller $\Delta g$ leads to better suppression of magnetic systematic errors. While $g$-factor differences have been studied for diatomic and linear triatomic molecules, there was a lack of data regarding the magnitude and specific contributions to $\Delta g$ for symmetric top molecules (specifically RaOCH$_3$) under experimental electric fields.

2. Methodology

The authors developed a theoretical method to calculate the electric-field-dependent $g$-factors for the first excited rotational level ($N=1$) of the RaOCH$_3$ molecule.

• Hamiltonian Construction: The total Hamiltonian ($\hat{H}$) was constructed in the molecular reference frame, comprising:
  • $\hat{H}_{mol}$: The molecular rotational Hamiltonian.
  • $\hat{H}_{hfs}$: Hyperfine interaction with hydrogen nuclei.
  • $\hat{H}_{S}$: Stark interaction with the external electric field.
  • $\hat{H}_{Z}$: Zeeman interaction with the external magnetic field.
• Basis Set and Symmetry: The calculation utilized a basis set of electronic-rotational-vibrational-nuclear spin wavefunctions. Special attention was paid to the Pauli principle for the hydrogen nuclei (spin $I=1/2$), which creates a strong correlation between the rotational quantum number $K$ and the total nuclear spin $I$.
  • $K=0$ corresponds to $I=3/2$.
  • $|K|=1$ corresponds to $I=1/2$.
  • Only the $|K|=1$ states form K-doublets relevant to the eEDM search.
• Numerical Diagonalization: The Hamiltonian was numerically diagonalized to obtain wavefunctions, energy levels, and $g$-factors. The calculation included electronic matrix elements for orbital angular momentum ($\hat{L}$) and spin ($\hat{S}$), derived from previous ab initio studies.
• Perturbation Analysis: The authors analyzed the specific perturbations (hyperfine and rotational interactions) that cause the slight difference in $g$-factors between the upper ($g_u$) and lower ($g_l$) K-doublet components.

3. Key Contributions

• First Calculation for Symmetric Tops: This work provides the first calculation of the $g$-factor difference ($\Delta g$) between K-doublet components in a symmetric top molecule (RaOCH$_3$).
• Identification of Perturbation Mechanisms: The paper identifies the specific physical mechanisms driving $\Delta g$:
  • For the $F=1$ states, the difference arises primarily from the mixing of $|N=1, J=3/2\rangle$ and $|N=2, J=3/2\rangle$ states via molecular interaction, which is sensitive to the precise values of electronic matrix elements.
  • For the $F=2$ states (the most promising for experiments), the difference is significantly smaller, arising from subtle differences in how positive and negative parity states are perturbed by higher rotational levels ($N=2$) and hyperfine interactions.
• Electric Field Dependence: The study explicitly maps how $\Delta g$ evolves as a function of the laboratory electric field ($E$), a critical parameter for experimental planning.

4. Results

The calculations were performed for the $N=1$ rotational level, focusing on Zeeman sublevels with $M_F = 1$ and $M_F = 2$.

• Magnitude of $\Delta g$:
  • For the $N=1, J=3/2, F=1$ levels, $\Delta g \approx 4.7 \times 10^{-4}$.
  • For the $N=1, J=3/2, F=2$ levels (the target for eEDM experiments), $\Delta g$ is extremely small. At an electric field of $E = 500$ mV/cm, $\Delta g \approx -2.3 \times 10^{-7}$.
• Ratio $\Delta g / g$:
  • For the $F=2$ state, the ratio $\Delta g / g \sim 10^{-6}$.
  • This is significantly more favorable (smaller) than values found in other systems:
    • ThO and HfF$^+$ (diatomic/linear): $\Delta g / g \sim 10^{-3}$.
    • YbOH (linear triatomic): $\Delta g / g \le 4 \times 10^{-5}$.
• Field Dependence:
  • The $g$-factor difference for $M_F=2$ increases almost linearly with the electric field.
  • The value is stable and independent of initial small differences because it is not the result of the cancellation of large numbers (unlike the field-free case).
  • The required electric field to polarize the molecule ($E \approx 250-500$ mV/cm) is sufficient to maximize the eEDM signal while keeping $\Delta g$ minimal.

5. Significance

• Experimental Viability: The extremely small ratio of $\Delta g / g \sim 10^{-6}$ for RaOCH$_3$ suggests that this molecule is exceptionally robust against systematic errors caused by magnetic field inhomogeneities. This makes it a superior candidate for next-generation eEDM experiments compared to current leading candidates like ThO or HfF$^+$.
• Laser Cooling Compatibility: Unlike diatomic molecules with $\Omega$-doubling (which have complex electronic structures), symmetric top molecules like RaOCH$_3$ in their ground electronic state are amenable to laser cooling. The combination of laser cooling (improving statistical sensitivity) and the tiny $\Delta g$ (suppressing systematic errors) positions RaOCH$_3$ as a highly promising system for future precision tests of fundamental physics.
• Theoretical Framework: The developed method for calculating $g$-factors in symmetric tops fills a critical gap in the literature, providing a template for analyzing other symmetric top molecules in the search for new physics.