Ab Initio Calculations of the Static and Dynamic Polarizability of BaOH

This paper presents high-precision relativistic coupled-cluster calculations of the static and dynamic polarizabilities for the ground and (010) vibrational states of BaOH, establishing a rigorous uncertainty quantification procedure and validating the results through comparison with experimental dipole moment data.

The Gist

Imagine a molecule, BaOH (Barium Monohydroxide), not as a static chemical formula, but as a tiny, squishy balloon floating in space. This balloon is made of three atoms: a heavy Barium, an Oxygen, and a Hydrogen.

Scientists are very interested in this specific balloon because they want to use it to hunt for a "ghost" in physics: the electron's electric dipole moment (eEDM). Think of the eEDM as a tiny, hidden imbalance in how an electron is shaped. Finding it would prove that our current understanding of the universe (the Standard Model) is incomplete and point us toward new, exciting physics.

However, to catch this ghost, scientists need to trap the BaOH balloon in a laser cage and hold it perfectly still. To do that, they need to know exactly how "squishy" the balloon is when you push on it with an electric field. This "squishiness" is called polarizability.

Here is a simple breakdown of what this paper did:

1. The Problem: We Need a Map for the Trap

To trap these molecules, scientists use lasers to create an "Optical Dipole Trap" (ODT)—essentially a bowl made of light. The depth of this bowl depends on how the molecule reacts to the laser light.

• If the molecule is too "stiff," the laser won't hold it.
• If it's too "squishy," it might get knocked out of the trap by the light itself.

The scientists needed a precise map of how BaOH reacts to different colors of light (frequencies), especially the specific infrared laser (1064 nm) they plan to use. Until now, they didn't have a reliable map.

2. The Solution: The Ultimate Simulation

Instead of building a physical trap and guessing, the team used a supercomputer to simulate the molecule from the ground up. They used a method called Coupled-Cluster Theory, which is like the "gold standard" of molecular modeling.

Think of it like this:

• The Ball: The Barium atom is a heavy, dense ball. Because it's so heavy, the laws of physics get a bit weird (relativity). The electrons zoom around so fast that they get heavier and squish closer to the nucleus. The team had to account for this "relativistic effect," or their map would be wrong.
• The Cloud: The electrons form a fuzzy cloud around the atoms. When an electric field (the laser) comes near, this cloud stretches and distorts. The team calculated exactly how much it stretches in different directions.

3. The "Uncertainty" Check: How Sure Are We?

In science, saying "the answer is 100" is useless if you don't know if it's actually 90 or 110. This paper is special because the authors didn't just give a number; they built a safety margin around it.

They tested their simulation by changing the "ingredients" of their calculation:

• Basis Sets: Imagine trying to draw a circle. Do you use 10 dots, 100 dots, or 1,000 dots? They tested different numbers of dots to see when the picture stopped changing.
• Relativity: They checked if ignoring the "fast electron" effects changed the result. (It did, significantly!)
• Vibrations: The molecule isn't a frozen statue; it wiggles and bends like a rubber band. They calculated how this wiggling changes the "squishiness."

By comparing all these variations, they created a "confidence interval." They can now say with high confidence: "The squishiness is X, give or take a tiny bit."

4. The Results: The Blueprint is Ready

The team produced a detailed blueprint for BaOH:

• Static Polarizability: How squishy it is when you push it gently (like a slow hand).
• Dynamic Polarizability: How squishy it is when you hit it with a fast-moving laser beam (the 1064 nm light).

They found that the molecule is much easier to squish from the side (perpendicular) than from the top (parallel). This is like a long, thin balloon: it's easy to squeeze the sides, but hard to squash the ends.

5. Why This Matters

This paper is the instruction manual for the next generation of physics experiments.

• For the Trap: Now that the NL-eEDM collaboration knows exactly how BaOH reacts to their specific laser, they can design a trap that is deep enough to hold the molecules but shallow enough not to knock them out.
• For Physics: With a better trap, they can hold the molecules longer and measure them more precisely. This increases their chances of spotting that "ghost" (the eEDM) and discovering new laws of the universe.

In a nutshell:
This paper is like an engineer running thousands of computer simulations to figure out exactly how much weight a specific bridge can hold before building it. They didn't just guess; they tested every variable (wind, material strength, vibration) to ensure the bridge (the experiment) won't collapse when they try to cross it (find new physics). They have successfully calculated the "squishiness" of the BaOH molecule, paving the way for a historic discovery.

Technical Summary

1. Problem Statement

The Barium Monohydroxide (BaOH) molecule is a leading candidate for next-generation experiments searching for the electron electric dipole moment (eEDM), a probe for physics beyond the Standard Model. To perform these experiments, BaOH molecules must be laser-cooled and trapped in optical dipole traps (ODTs).

• The Challenge: Designing efficient trapping and cooling schemes requires precise knowledge of the molecule's electric dipole polarizability (both static and dynamic).
• The Gap: While experimental data exists for some transitions, accurate theoretical values for the static and dynamic polarizability tensors of BaOH, including rigorous uncertainty quantification and vibrational corrections, were previously unavailable.
• Specific Needs: Researchers need values for the ground vibrational state $(000)$ and the bending excited state $(010)$, as the latter allows for the creation of near-zero g-factor spin states that suppress magnetic field sensitivity.

2. Methodology

The authors employed a multi-faceted ab initio approach using Relativistic Coupled-Cluster (CC) theory to ensure high accuracy for the heavy Barium atom.

A. Electronic Structure Methods

• Static Polarizability ($\alpha(0)$):
  • Calculated using the Finite-Field (FF) method within the DIRAC23 software package.
  • Hamiltonian: Primarily used the Exact Two-Component (X2C) Hamiltonian to account for scalar relativistic effects, validated against the full four-component Dirac-Coulomb (DC) Hamiltonian.
  • Correlation: Used CCSD(T) (Coupled-Cluster with Single, Double, and perturbative Triple excitations).
  • Basis Sets: Employed large, uncontracted Dyall basis sets (s-, d-, and t-augmented up to quadruple-$\zeta$ quality).
  • Relativistic/QED Effects: Systematically included scalar relativistic effects, Breit interactions (estimated via atomic data), and Quantum Electrodynamics (QED) corrections (vacuum polarization and self-energy).
• Dynamic Polarizability ($\alpha(\omega)$):
  • Calculated using Linear-Response (LR) theory within the CFOUR software package at the CCSD level.
  • Approximation: Used Effective Core Potentials (ECP) for Barium to handle scalar relativistic effects computationally.
  • Frequency Range: Calculated from the static limit up to frequencies beyond the first and second electronic transitions ($A^2\Pi \leftarrow X^2\Sigma$ and $B^2\Sigma \leftarrow X^2\Sigma$).
  • Interpolation: Used a parametrized fit model to interpolate values at the specific laser wavelength of 1064 nm (282 THz), commonly used for ODTs.

B. Vibrational Corrections

• Calculations were performed for the equilibrium geometry and 10 displaced nuclear configurations along normal mode vectors.
• The Numerov–Cooley (NC) procedure was used to solve the nuclear Schrödinger equation, averaging the polarizability over the vibrational wavefunctions for the $(000)$ and $(010)$ states.
• This provided vibrational corrections to the electronic polarizabilities.

C. Uncertainty Quantification

A rigorous procedure was developed to estimate uncertainties by varying computational parameters:

• Basis set cardinality and augmentation.
• Active space size (number of correlated electrons).
• Level of electron correlation (CCSD vs. CCSD(T)).
• Relativistic and QED treatment.
• Uncertainties were combined using the Square Root of the Sum of Squares (SRSS) method, assuming mutual independence.

3. Key Contributions

1. Systematic Uncertainty Estimation: Unlike many previous studies, this work provides a comprehensive error budget for polarizability calculations in heavy polyatomic molecules, identifying basis set incompleteness and higher-order excitations as dominant error sources.
2. Validation of Approximations: The study validated the use of the X2C Hamiltonian and ECP approaches by comparing them against full four-component DC calculations, showing excellent agreement for polarizability properties.
3. Vibrational State Specificity: Provided distinct polarizability values for the ground state $(000)$ and the bending excited state $(010)$, which is critical for the specific "near-zero g-factor" experimental proposals.
4. Dynamic Polarizability at 1064 nm: Delivered precise values for the dynamic polarizability at the specific laser wavelength used in current and proposed trapping experiments, bridging the gap between static theory and experimental application.

4. Key Results

The final recommended values (in atomic units, a.u.) are summarized below:

Property	State	Static ($\omega=0$)	Dynamic ($\lambda=1064$ nm)
Dipole Moment ($\mu_z$)	$(000)$	0.583(11)	N/A
Parallel Polarizability ($\alpha_\parallel$)	$(000)$	200.8(24)	358(31)
	$(010)$	201.1(24)	360(31)
Perpendicular Polarizability ($\alpha_\perp$)	$(000)$	297(5)	715(29)
	$(010)$	297(5)	718(29)

• Dipole Moment: The calculated $\mu_z = 0.583(11)$ a.u. agrees well with the experimental value of $0.563(16)$ a.u., validating the theoretical framework.
• Anisotropy: The perpendicular component ($\alpha_\perp$) is significantly larger than the parallel component ($\alpha_\parallel$) due to the weaker binding of electrons perpendicular to the molecular axis.
• Vibrational Effects: The vibrational corrections are small (less than 1% change), but non-negligible for high-precision work.
• Dynamic Behavior: The dynamic polarizability increases significantly at 1064 nm compared to the static limit, particularly for the perpendicular component, due to proximity to electronic transitions.

5. Significance and Implications

• Experimental Design: The calculated values allow the NL-eEDM collaboration and other groups to accurately model the depth of optical dipole traps. For BaOH, the results suggest that a trap depth of $\sim 1$ mK can be achieved with $\sim 7$ W of laser power at a 25 $\mu$m waist.
• Scattering Rates: The dynamic polarizability data enables the estimation of residual photon scattering rates. The authors estimate that a 1 mK deep trap scatters a photon every second, while a 10 mK trap scatters 10 times per second. This informs the trade-off between trap depth and heating/decoherence.
• Benchmarking: The study serves as a benchmark for relativistic coupled-cluster methods in heavy polyatomic systems, demonstrating that high-accuracy predictions for response properties are achievable even with complex relativistic and correlation effects.
• Future Physics: By providing the necessary molecular constants, this work directly supports the feasibility of using BaOH in optical lattices to search for new physics (eEDM) with unprecedented sensitivity.