🔬 atomic physics

Atomic structure calculations for constraining new electron-electron forces

This paper presents a general atomic structure approach incorporating many-body corrections via the time-dependent Hartree-Fock method to constrain new electron-electron forces, establishing bounds on scalar-pseudoscalar and vector-axial vector couplings that rule out specific parameter regions and provide an updated Standard Model calculation for cesium parity non-conserving transitions.

Original authors: Cameron J. West, Benjamin M. Roberts

Published 2026-08-28

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Cameron J. West, Benjamin M. Roberts

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The universe is built on a foundation of known particles and forces, a framework physicists call the Standard Model. It explains how matter behaves and interacts with stunning precision, yet it leaves some of the cosmos's biggest mysteries unsolved. Why is there more matter than antimatter? What is dark matter made of? To find answers, scientists look for tiny cracks in the Standard Model's armor, searching for new particles or forces that might hide in plain sight. One promising place to look is within the atom itself, specifically in the way electrons interact with one another. While we understand how electrons behave around the nucleus, the subtle forces they exert on each other are incredibly difficult to measure. If a new, invisible particle were mediating a force between electrons, it would leave a faint fingerprint on the atom's energy or its shape. Detecting this requires not just sensitive instruments, but a way to calculate exactly how an atom should behave without any new physics, so that any deviation stands out clearly.

A team of researchers at the University of Queensland has developed a new, more powerful way to perform these calculations. They created a general method to predict how atoms would react if a new force existed between electrons, regardless of how heavy the particle carrying that force might be. Previous attempts to model these interactions had to rely on extreme assumptions: either the new particle was massless, or it was so heavy it acted like a direct contact point. The new approach bridges this gap, allowing for accurate predictions across the entire range of possible masses, including a middle ground that had been largely unexplored. By applying this method to complex atomic systems, the researchers were able to include the messy, collective behavior of many electrons, a factor often ignored in simpler models. This level of detail is crucial because the electrons in an atom do not act alone; they constantly influence one another, and ignoring this can lead to incorrect conclusions about whether a new force is present.

The researchers tested their method by looking for two specific types of hypothetical interactions. The first involves a force that would cause an atom to develop a tiny electric dipole moment, essentially making the atom slightly positive on one end and negative on the other. In the Standard Model, atoms do not have this property, so finding one would be a clear sign of new physics. The team calculated how this effect would appear in heavy atoms like thallium and cesium, accounting for the complex dance of electrons within them. They compared their theoretical predictions with the most precise experimental measurements available. The results showed that while their calculations were robust and included important corrections, the current experimental limits on these electric dipole moments are already quite tight. Their work confirmed that any new force of this type must be weaker than previously thought, or the particle carrying it must have a specific mass that makes it harder to detect in these particular atoms.

The second part of their investigation focused on a different kind of interaction that could cause an electron to jump between energy levels in a way that is normally forbidden. This is known as a parity-non-conserving transition, where the atom behaves differently depending on which way it is mirrored. The researchers focused on a specific jump in the cesium atom, from a lower energy state to a higher one, which has been measured with extreme precision. They calculated the contribution that the Standard Model itself makes to this jump, specifically the part caused by electrons interacting with each other via the weak nuclear force. Their updated calculation agreed with previous values but included more complete corrections for how the electrons move together. By comparing this refined theoretical value with the experimental measurement, they were able to set new limits on a hypothetical new particle that could be mediating this interaction.

The most significant finding emerged from this second analysis. The researchers were able to rule out a wide range of possibilities for a new type of particle that interacts with electrons but not with the heavier particles in the nucleus. For particles with a mass greater than about 10 million electron volts, they determined that the strength of the interaction between the new particle and electrons cannot exceed a specific threshold. This effectively closes a gap in our knowledge, eliminating a region of mass and interaction strength that had previously been unconstrained. It means that if such a particle exists, it must be even more elusive than scientists had hoped, or it does not exist in the form they were looking for.

This work does not just provide a list of numbers; it offers a refined tool for the future. The method the team developed can be applied to any atom and any type of new force, making it a versatile instrument for the ongoing search for physics beyond the Standard Model. By showing how to properly account for the collective behavior of electrons, they have removed a source of uncertainty that could have masked a discovery. While their specific results for the electric dipole moments did not reveal new physics, they sharpened the focus for future experiments. The new constraints on the heavy-mass region for the other interaction type provide a clear target for theorists and experimentalists alike. The study underscores that the path to discovering new particles often lies in the most subtle details of the familiar world, requiring calculations that are as precise as the measurements they aim to explain.

Technical Summary: Atomic Structure Calculations for Constraining New Electron-Electron Forces

Problem Statement
The Standard Model fails to account for phenomena such as dark matter and matter-antimatter asymmetry, motivating searches for new leptonic interactions. While precision atomic physics can probe electron-electron ($ee$) and electron-nucleus ( $eN$ ) interactions, investigations into the complex two-body $ee$ case, particularly those violating parity (P), remain scarce. Existing theoretical approaches for new $ee$ interactions typically rely on approximations of the mediator mass ( $m$ ), specifically the massless ( $m \to 0$ ) or contact ( $m \to \infty$ ) limits. These approximations are often invalid in the intermediate keV to MeV region, where constraints are sparse. Furthermore, previous calculations often lacked rigorous inclusion of many-body effects necessary for accurate predictions in heavy atoms.

Methodology
The authors present a general theoretical framework for calculating new $ee$ interactions using atomic structure theory, valid for any mediator mass. The methodology involves:

Interaction Potential: Deriving a new expression for the interaction potential based on an effective Lagrangian for scalar ( $S=0$ ) and vector ( $S=1$ ) bosons. This potential takes the form of a Yukawa potential, $V_{0}^{\alpha\beta} \propto \frac{e^{-\mu r_{jk}}}{r_{jk}}$ , which is valid across all mass regimes and reduces to limiting values appropriately.
Many-Body Corrections:
- Time-Dependent Hartree-Fock (TDHF): To account for first-order corrections induced by new P-odd interactions (which cannot be directly included in the static Dirac-Hartree-Fock Hamiltonian without breaking symmetry), the authors apply TDHF equations. This is equivalent to the Random Phase Approximation (RPA) and yields corrections to the DHF potential and wavefunctions.
- All-Order Correlation: Significant higher-order electron correlation effects are included using the all-order correlation potential method.
- Extension to Multi-Valence: The framework is extended to multi-valence atoms using the Combined Configuration Interaction and Many-Body Perturbation Theory (CI+MBPT) method.
Matrix Element Calculation:
- For P-even interactions, energy shifts are calculated via expectation values.
- For P-odd interactions (Scalar-Pseudoscalar [SP] and Vector-Axial Vector [VA]), the method calculates off-diagonal matrix elements mixing states of opposite parity. This involves a sum-over-states approach to determine induced atomic Electric Dipole Moments (EDMs) or parity non-conserving (PNC) transition amplitudes.
- In the contact limit ( $\mu \to \infty$ ), a Fierz identity is utilized to simplify spatial terms for VA interactions, allowing the neglect of spatial contributions in favor of time-like terms.

Key Contributions

General Mass Validity: The derivation of interaction potentials valid for all mediator masses, bridging the gap between massless and contact limits and enabling accurate calculations in the previously unconstrained intermediate mass region.
Rigorous Many-Body Treatment: The inclusion of TDHF and all-order correlation corrections specifically for new $ee$ interactions, addressing a deficiency in previous calculations where such effects were often omitted or approximated.
Updated Standard Model Baseline: An updated calculation of the Standard Model $ee$ contribution (via $Z^0$ exchange) to the $6s \to 7s$ PNC transition amplitude in Cesium ( $^{133}\text{Cs}$ ), incorporating improved many-body contributions.

Results

Scalar-Pseudoscalar (SP) Interactions:
- Calculations of atomic EDMs for single-valence systems ( $^{173}\text{Yb}^+$ , $^{205}\text{Tl}$ , alkali metals) and the divalent $^{199}\text{Hg}$ system were performed.
- Results show a strong dependence on atomic number $Z$ , with heavier systems being more affected.
- A sign change in the EDM is observed in the intermediate mass range (near the electron mass), creating regions of insensitivity to experimental constraints.
- Bounds on the coupling strength $|g_S^e g_P^e|$ were extracted by comparing calculations with experimental EDM limits for $^{133}\text{Cs}$ and $^{205}\text{Tl}$ . While these bounds are less stringent than existing molecular and diamagnetic limits, they account for many-body effects and reveal insensitivity zones.
Vector-Axial Vector (VA) Interactions:
- New bounds were extracted for a leptophilic vector boson by comparing theoretical $ee$ contributions to experimental measurements of the $6s \to 7s$ PNC transition in $^{133}\text{Cs}$ .
- The study rules out any coupling-mass ratio $|g_V^e g_A^e|/m_X^2$ greater than $10^{-11} \text{ MeV}^{-2}$ for mediator masses $m_X \gtrsim 10 \text{ MeV}$ .
- These represent the first direct bounds obtained outside the ultralight mass region for this specific interaction type.

Significance and Claims
The paper claims to provide a robust improvement over existing methods by removing the reliance on mass-limit approximations and rigorously incorporating many-body corrections. The authors emphasize that their approach allows for the investigation of the intermediate mass range where previous constraints were non-existent or weak.

Specifically, the work:

Demonstrates that many-body effects (particularly in trivalent systems like Thallium) significantly enhance sensitivity to new physics, suggesting that atoms with lower experimental uncertainties and specific electronic configurations are optimal for future probes.
Provides the first direct constraints on leptophilic vector bosons in the high-mass sector ( $m_X > 10 \text{ MeV}$ ).
Offers an updated, more precise estimate of the Standard Model $ee$ contribution to the Cesium PNC transition, which is crucial for isolating new physics signals.
Identifies that while current atomic EDM bounds are competitive with molecular limits, renewed interest in atomic measurements is warranted given the potential for improved experimental precision.

The authors note that spatial contributions to the VA interaction (Breit-like terms) were neglected in the contact limit approximation and are reserved for future work, as are calculations for P-even interactions which induce energy shifts in light systems.

Technical Summary: Atomic Structure Calculations for Constraining New Electron-Electron Forces

More like this