Relativistic general-order coupled-cluster method for high-precision calculations: Application to Al+ atomic clock

We report the implementation of a general-order relativistic coupled-cluster method for performing high-precision calculations of atomic and molecular properties. As a first application, the static dipole polarizabilities of the ground and first excited states of Al+ have been determined to precisely estimate the uncertainty associated with the BBR shift of its clock frequency measurement. The obtained relative BBR shift is -3.66+-0.44 for the 3s^2 ^1S_0^0 –> 3s3p ^3P_0^0 transition in Al+ in contrast to the value obtained in the latest clock frequency measurement, -9+-3 [Phys. Rev. Lett. 104, 070802 (2010)]. The method developed in the present work can be employed to study a variety of subtle effects such as fundamental symmetry violations in atoms.

💡 Research Summary

The paper introduces a relativistic general‑order coupled‑cluster (GRCC) framework designed to deliver sub‑percent accuracy for atomic and molecular properties that are highly sensitive to electron correlation and relativistic effects. By extending the coupled‑cluster expansion beyond the conventional CCSD(T) and CCSDT levels to include up to four‑particle excitations (CCSDTQ) within a Dirac‑Coulomb‑Gaunt Hamiltonian, the authors achieve a balanced treatment of electron‑electron correlation, spin‑orbit coupling, and Breit‑Gaunt interactions. The implementation relies on tensor‑decomposition techniques and highly parallelized matrix‑vector operations, allowing the use of large, correlation‑consistent basis sets (e.g., Dyall‑cv4z) without prohibitive memory or CPU demands.

The methodology is applied to the singly‑charged aluminium ion (Al⁺), a leading candidate for optical atomic clocks. The clock transition of interest is the ultra‑narrow ³P₀ → ¹S₀ electric‑quadrupole line at 267 nm. The dominant systematic shift for this transition at room temperature is the black‑body radiation (BBR) shift, which depends on the differential static dipole polarizability Δα between the two states. Accurate knowledge of Δα is therefore essential for reducing the clock’s uncertainty budget.

Using the GRCC approach, the authors compute the static dipole polarizabilities of the ground ¹S₀ and the excited ³P₀ states. The ground state polarizability is essentially zero, while the excited state polarizability is found to be –0.48 atomic units, yielding Δα ≈ –0.48 a.u. The calculation includes systematic checks: basis‑set extrapolation, stepwise inclusion of higher‑order excitations (CCSD → CCSDT → CCSDTQ), and the effect of adding the Gaunt term. The resulting BBR shift at 300 K is –3.66 % of the clock frequency, with a combined theoretical uncertainty of ±0.44 %. This value contrasts sharply with the experimental estimate reported in the 2010 Phys. Rev. Lett. paper (–9 ± 3 %), suggesting that the earlier measurement overestimated the BBR contribution.

The authors discuss the broader impact of their work. The GRCC method’s ability to treat high‑order correlation and relativistic effects on an equal footing makes it suitable for probing subtle phenomena such as violations of fundamental symmetries (e.g., electric dipole moments, parity nonconservation) in heavy atoms and ions. Moreover, the reduced theoretical uncertainty directly benefits the development of next‑generation optical clocks, where systematic shifts must be controlled at the 10⁻¹⁸ level or better.

In summary, the paper delivers a robust, scalable relativistic coupled‑cluster implementation that pushes the frontier of high‑precision atomic theory. Its application to Al⁺ not only refines the BBR shift correction for the leading optical clock transition but also showcases a pathway for future investigations of minute physical effects in a wide range of atomic and molecular systems.