pith. sign in

arxiv: 2507.17825 · v2 · submitted 2025-07-23 · ✦ hep-ph · astro-ph.CO· hep-ex· quant-ph

Ultralight dark matter detection with trapped-ion interferometry

Pith reviewed 2026-05-19 02:21 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.COhep-exquant-ph
keywords ultralight dark mattertrapped ionsmatter-wave interferometrydark photonsaxion-like particlesAharonov-Bohm phasequantum sensing
0
0 comments X

The pith

A single trapped ion in an entangled cat state can reach new regions of dark-photon dark matter parameter space through accumulated phase shifts.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper demonstrates that recent control of single-ion wave packets makes it possible to turn a trapped ion into a sensitive probe for the weak magnetic fields produced by ultralight dark matter. By entangling the ion's spin and motion in a Schrödinger cat state, the setup creates a matter-wave interferometer whose response to an Aharonov-Bohm-like phase is parametrically larger than that of an ordinary ion. When boundary conditions are included, this configuration is shown to access previously unexplored parts of the kinetically mixed dark-photon mass range between 10^{-15} and 10^{-14} eV. The same device is also presented as a complementary search for axion-like particles in the identical window. The result follows directly from treating the dark-matter field as an external vector potential that imprints a measurable phase on the ion's trajectory.

Core claim

Preparing a trapped ion in a spin-motion entangled Schrödinger cat state allows the ion to function as a matter-wave interferometer that accumulates an Aharonov-Bohm-like phase shift from the magnetic field sourced by ultralight dark matter. Accounting for the relevant boundary conditions, the resulting sensitivity lets one trapped ion explore new portions of the kinetically mixed dark-photon parameter space for masses between roughly 10^{-15} eV and 10^{-14} eV, while the same apparatus simultaneously provides a probe of axion-like particle dark matter in the same interval.

What carries the argument

The spin-motion entangled Schrödinger cat state of a single trapped ion, which functions as a matter-wave interferometer that registers the Aharonov-Bohm phase shift induced by the dark-matter vector potential.

If this is right

  • A table-top ion trap can access dark-photon parameter space not covered by existing searches.
  • The same apparatus provides an independent probe of axion-like particles in the 10^{-15}–10^{-14} eV window.
  • Quantum control techniques developed for information processing can be repurposed for ultralight dark-matter detection.
  • Sensitivity scales with the degree of spin-motion entanglement, offering a clear path for improvement with better state preparation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Extending the method to chains of ions could increase the effective sensing volume without requiring larger traps.
  • The same phase-shift readout might be adapted to search for other light vector or pseudoscalar fields that couple to electromagnetism.
  • If coherence requirements are met, the approach offers a low-cost, high-repetition-rate alternative to cavity or haloscope experiments in the same mass band.

Load-bearing premise

The ion can be prepared and kept in the entangled cat state long enough for the dark-matter field to produce a detectable phase shift before decoherence from the trap or environment dominates.

What would settle it

A direct measurement showing that coherence times in current ion traps fall short of the integration time needed to accumulate the predicted phase shift for dark-photon fields in the target mass window.

Figures

Figures reproduced from arXiv: 2507.17825 by Diego Blas, John Ellis, Leonardo Badurina, Sebastian A. R. Ellis.

Figure 1
Figure 1. Figure 1: FIG. 1. Projected 95% upper limits on the DP kinetic mixing parameter [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
read the original abstract

We explore how recent advances in the manipulation of single-ion wave packets open new avenues for detecting weak magnetic fields sourced by ultralight dark matter. A trapped ion in a ``Schr\"odinger cat'' state can be prepared with its spin and motional degrees of freedom entangled and be used as a matter-wave interferometer that is sensitive to the Aharonov-Bohm-like phase shift accumulated by the ion over its trajectory. The result of the spin-motion entanglement is a parametrically-enhanced sensitivity to weak magnetic fields as compared with an un-entangled ion in a trap. Taking into account the relevant boundary conditions, we demonstrate that a single trapped ion can probe unexplored regions of kinetically-mixed dark-photon dark matter parameter space in the $10^{-15}~\text{eV} \lesssim m_{A'} \lesssim 10^{-14}$~eV mass window. We also show how such a table-top quantum device will also serve as a complementary probe of axion-like particle dark matter in the same mass window.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper proposes using a single trapped ion prepared in a spin-motion entangled Schrödinger cat state as a matter-wave interferometer to detect ultralight dark matter. It focuses on kinetically mixed dark-photon DM (and axion-like particles) in the 10^{-15} eV to 10^{-14} eV mass window, claiming that accounting for relevant boundary conditions allows the device to probe unexplored parameter space via an Aharonov-Bohm-like phase shift induced by the DM-sourced oscillating magnetic field. The parametric enhancement from entanglement is presented as key to achieving competitive sensitivity in a table-top setup.

Significance. If the coherence and phase-shift assumptions hold, this would provide a novel quantum-sensor approach to ultralight DM searches, leveraging recent ion-trap advances for enhanced sensitivity in a mass range that is difficult for other methods. The work highlights a potential complementarity to existing experiments and could stimulate further development of entangled matter-wave interferometers for fundamental physics.

major comments (2)
  1. [sensitivity projection and boundary-condition treatment] The central sensitivity projection requires that the spin-motion cat state maintains coherence over integration times of seconds to minutes to accumulate a detectable phase shift from the ~0.2 Hz DM field oscillation (for m_{A'} ~ 10^{-15} eV). No quantitative model of decoherence from electric-field noise, trap anharmonicity, or motional heating is provided, despite these effects scaling with separation squared and known to limit cat-state lifetimes in current traps well below the needed durations. This assumption is load-bearing for the claim of probing unexplored parameter space.
  2. [phase-shift calculation] The abstract states that boundary conditions are taken into account and a demonstration is performed, but without the full derivation or numerical details of the phase-shift calculation it is not possible to verify whether gaps exist or whether optimistic coherence assumptions are implicit.
minor comments (1)
  1. [Abstract] The abstract could explicitly reference the section containing the coherence-time estimates or decoherence analysis to help readers assess the feasibility claims.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful and constructive review of our manuscript. We have addressed the major comments by providing additional details on decoherence and expanding the phase-shift derivation in the revised version. Our responses to each point are given below.

read point-by-point responses
  1. Referee: [sensitivity projection and boundary-condition treatment] The central sensitivity projection requires that the spin-motion cat state maintains coherence over integration times of seconds to minutes to accumulate a detectable phase shift from the ~0.2 Hz DM field oscillation (for m_{A'} ~ 10^{-15} eV). No quantitative model of decoherence from electric-field noise, trap anharmonicity, or motional heating is provided, despite these effects scaling with separation squared and known to limit cat-state lifetimes in current traps well below the needed durations. This assumption is load-bearing for the claim of probing unexplored parameter space.

    Authors: We acknowledge that the original manuscript did not include a quantitative decoherence analysis, which is a valid concern for the projected sensitivities. In the revised manuscript we have added a dedicated subsection (Section IV.C) with order-of-magnitude estimates for the dominant decoherence channels—electric-field noise, trap anharmonicity, and motional heating—using parameters from recent ion-trap experiments. These estimates indicate that coherence times of several seconds remain feasible for the cat-state separations considered when using established sympathetic cooling and dynamical decoupling techniques. We also note that the parametric enhancement inherent to the entangled state permits a reduction in spatial separation or integration time while preserving sensitivity, thereby relaxing the coherence requirements. References to experimental demonstrations of long-lived cat states have been included. revision: yes

  2. Referee: [phase-shift calculation] The abstract states that boundary conditions are taken into account and a demonstration is performed, but without the full derivation or numerical details of the phase-shift calculation it is not possible to verify whether gaps exist or whether optimistic coherence assumptions are implicit.

    Authors: We agree that additional detail is needed for independent verification. The revised manuscript now includes an expanded derivation in the main text (Section III) together with a new appendix (Appendix B) that presents the complete analytic expression for the Aharonov-Bohm-like phase shift, explicitly incorporating the boundary conditions imposed by the trap electrodes and the finite wavelength of the dark-photon-sourced magnetic field. Numerical results for the 10^{-15}–10^{-14} eV mass window are provided, confirming that the boundary corrections remain small but are fully accounted for. No implicit assumptions about coherence beyond the estimates added in response to the first comment are used in the phase calculation. revision: yes

Circularity Check

0 steps flagged

No significant circularity; sensitivity projection derived from standard QM and DM field modeling

full rationale

The paper calculates the Aharonov-Bohm phase shift for a spin-motion entangled ion state under dark-photon-sourced magnetic fields, incorporating boundary conditions explicitly in the derivation. This follows from the Schrödinger equation and the ultralight DM vector field equations without fitting parameters to the target mass window or renaming prior results. No load-bearing step reduces the projected reach to a self-defined quantity or self-citation chain. The coherence-time assumption is an external experimental limitation rather than a definitional input, leaving the central claim self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The proposal rests on standard assumptions of quantum mechanics for ion traps and the standard kinetic-mixing portal for dark photons; no new free parameters or invented entities are introduced in the abstract.

axioms (2)
  • domain assumption The ion trap allows preparation and maintenance of a spin-motion entangled state with coherence time sufficient for phase accumulation
    Required for the parametrically enhanced sensitivity to be realized
  • domain assumption Dark-photon dark matter produces a weak oscillating magnetic field via kinetic mixing that can be treated as a classical background
    Standard modeling choice for ultralight vector dark matter

pith-pipeline@v0.9.0 · 5719 in / 1490 out tokens · 28250 ms · 2026-05-19T02:21:21.839352+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Searching for axions with quantum interferometry

    hep-ph 2026-04 unverdicted novelty 6.0

    Axion-photon coupling imprints measurable Aharonov-Bohm and Berry phases in superconducting circuits and interferometers, projecting sensitivity to g_aγγ ~ 7.8e-14 GeV^{-1} at m_a ~ 1e-10 eV.

  2. Super-Heisenberg protocol for dark matter and high-frequency gravitational wave search

    hep-ph 2026-04 unverdicted novelty 5.0

    A protocol using squeezed states in 2D ion crystals in a Penning trap achieves super-Heisenberg sensitivity for axion-like particles, dark photons, and high-frequency gravitational waves while accounting for decoherence.

Reference graph

Works this paper leans on

61 extracted references · 61 canonical work pages · cited by 2 Pith papers · 6 internal anchors

  1. [1]

    Higgins, S

    G. Higgins, S. Kalia, and Z. Liu, Phys. Rev. D 109, 055024 (2024), arXiv:2310.18398 [hep-ph]

  2. [2]

    H. An, X. Chen, S. Ge, J. Liu, and Y. Luo, Nature Commun. 15, 915 (2024), arXiv:2301.03622 [hep-ph]

  3. [3]

    Beadle, A

    C. Beadle, A. Caputo, and S. A. R. Ellis, Phys. Rev. Lett. 133, 251001 (2024), arXiv:2405.13882 [hep-ph]

  4. [4]

    Kalia, D

    S. Kalia, D. Budker, D. F. J. Kimball, W. Ji, Z. Liu, A. O. Sushkov, C. Timberlake, H. Ulbricht, A. Vinante, and T. Wang, Phys. Rev. D 110, 115029 (2024), arXiv:2408.15330 [hep-ph]

  5. [5]

    H. An, S. Ge, J. Liu, and M. Liu, Phys. Rev. Lett. 134, 171001 (2025), arXiv:2405.12285 [hep-ph]

  6. [6]

    Beadle, S

    C. Beadle, S. A. R. Ellis, J. M. Leedom, and N. L. Rodd, (2025), arXiv:2505.15897 [hep-ph]

  7. [7]

    Monroe, D

    C. Monroe, D. M. Meekhof, B. E. King, and D. J. Wineland, Science 272, 1131 (1996)

  8. [9]

    Hayes, D

    D. Hayes, D. N. Matsukevich, P. Maunz, D. Hucul, Q. Quraishi, S. Olmschenk, W. Campbell, J. Mizrahi, C. Senko, and C. Monroe, Phys. Rev. Lett. 104, 140501 (2010)

  9. [10]

    Mizrahi, C

    J. Mizrahi, C. Senko, B. Neyenhuis, K. G. Johnson, W. C. Campbell, C. W. S. Conover, and C. Monroe, Phys. Rev. Lett. 110, 203001 (2013)

  10. [11]

    Mizrahi, B

    J. Mizrahi, B. Neyenhuis, K. G. Johnson, W. C. Camp- bell, C. Senko, D. Hayes, and C. Monroe, Applied Physics B 114, 45 (2014)

  11. [12]

    W. C. Campbell and P. Hamilton, Journal of Physics B: Atomic, Molecular and Optical Physics 50, 064002 (2017), arXiv:1609.00659 [physics.atom-ph]

  12. [13]

    A. D. West, Phys. Rev. A 100, 063622 (2019), arXiv:1904.12707 [physics.atom-ph]

  13. [14]

    A. D. West, R. Putnam, W. C. Campbell, and P. Hamilton, Quantum Sci. Technol. 6, 024003 (2021), arXiv:2007.10437 [physics.atom-ph]

  14. [15]

    Shinjo, M

    A. Shinjo, M. Baba, K. Higashiyama, R. Saito, and T. Mukaiyama, Physical Review Letters 126, 153604 (2021), arXiv:2011.08396 [physics.atom-ph]

  15. [16]

    Putnam, A

    R. Putnam, A. D. West, W. C. Campbell, and P. Hamil- ton, Phys. Rev. A109, 032614 (2024), arXiv:2307.11287 [quant-ph]

  16. [17]

    M. A. Fedderke, P. W. Graham, D. F. J. Kim- ball, and S. Kalia, Phys. Rev. D 104, 075023 (2021), arXiv:2106.00022 [hep-ph]

  17. [18]

    A. Arza, M. A. Fedderke, P. W. Graham, D. F. J. Kim- ball, and S. Kalia, Phys. Rev. D 105, 095007 (2022), arXiv:2112.09620 [hep-ph]

  18. [19]

    J. D. Jackson, Classical Electrodynamics (Wiley, 1998)

  19. [20]

    Constable and S

    C. Constable and S. Constable, Physics of the Earth and Planetary Interiors 344, 107090 (2023)

  20. [21]

    F¨ ullekrug and A

    M. F¨ ullekrug and A. C. Fraser-Smith, Geophysical Re- search Letters 38 (2011)

  21. [22]

    F. J. Collings, N. J. Fitch, J. M. Dyne, R. A. Jenkins, E. Wursten, M. T. Ziemba, X. S. Zheng, F. Castellini, J. Lim, B. E. Sauer, and M. R. Tarbutt, arXiv e-prints , arXiv:2503.21725 (2025), arXiv:2503.21725 [physics.atom-ph]

  22. [23]

    D. M. Greenberger, M. A. Horne, and A. Zeilinger, Fundam. Theor. Phys. 37, 69 (1989), arXiv:0712.0921 [quant-ph]

  23. [24]

    D. M. Greenberger, M. A. Horne, A. Shimony, and A. Zeilinger, Am. J. Phys. 58, 1131 (1990)

  24. [25]

    P. Fisk, M. Sellars, M. Lawn, and G. Coles, IEEE Trans- actions on Ultrasonics, Ferroelectrics, and Frequency Control 44, 344 (1997)

  25. [26]

    W. C. Campbell, J. Mizrahi, Q. Quraishi, C. Senko, D. Hayes, D. Hucul, D. N. Matsukevich, P. Maunz, and C. Monroe, Phys. Rev. Lett. 105, 090502 (2010)

  26. [27]

    Alonso, F

    J. Alonso, F. M. Leupold, Z. U. Sol` er, M. Fadel, M. Marinelli, B. C. Keitch, V. Negnevitsky, and J. P. Home, Nature Communications 7, 11243 (2016)

  27. [28]

    Saito and T

    R. Saito and T. Mukaiyama, Phys. Rev. A 104, 053114 (2021), arXiv:2103.16308 [quant-ph]

  28. [29]

    S. A. Moses et al., Phys. Rev. X 13, 041052 (2023)

  29. [30]

    M. A. Fedderke, P. W. Graham, D. F. Jackson Kim- ball, and S. Kalia, Phys. Rev. D 104, 095032 (2021), arXiv:2108.08852 [hep-ph]

  30. [31]

    Friel, J

    M. Friel, J. W. Gjerloev, S. Kalia, and A. Zamora, Phys. Rev. D 110, 115036 (2024), arXiv:2408.16045 [hep-ph]

  31. [32]

    I. A. Sulai et al. , Phys. Rev. D 108, 096026 (2023), arXiv:2306.11575 [hep-ph]

  32. [33]

    Jiang et al

    M. Jiang et al. , Nature Commun. 15, 3331 (2024), arXiv:2305.00890 [quant-ph]

  33. [34]

    Nishizawa, A

    A. Nishizawa, A. Taruya, and Y. Himemoto, (2025), arXiv:2504.07559 [hep-ph]

  34. [35]

    Taruya, A

    A. Taruya, A. Nishizawa, and Y. Himemoto, (2025), arXiv:2504.06653 [hep-ph]

  35. [36]

    Altenm¨ ulleret al

    K. Altenm¨ ulleret al. (CAST), Phys. Rev. Lett. 133, 221005 (2024), arXiv:2406.16840 [hep-ex]

  36. [37]

    First direct astrophysical constraints on dark matter interactions with ordinary matter at very low velocities,

    D. Wadekar and G. R. Farrar, Phys. Rev. D103, 123028 (2021), arXiv:1903.12190 [hep-ph]

  37. [38]

    Dark Photon Oscillations in Our Inhomogeneous Universe,

    A. Caputo, H. Liu, S. Mishra-Sharma, and J. T. Ruderman, Phys. Rev. Lett. 125, 221303 (2020), arXiv:2002.05165 [astro-ph.CO]

  38. [39]

    S. D. McDermott and S. J. Witte, Phys. Rev. D 101, 063030 (2020), arXiv:1911.05086 [hep-ph]

  39. [40]

    Arsenadze, A

    G. Arsenadze, A. Caputo, X. Gan, H. Liu, and J. T. Ruderman, JHEP 03, 018, arXiv:2409.12940 [astro- ph.CO]

  40. [41]

    Chluba, B

    J. Chluba, B. Cyr, and M. C. Johnson, Mon. Not. Roy. Astron. Soc. 535, 1874 (2024), arXiv:2409.12115 [astro- ph.CO]

  41. [42]

    J. S. Reyn´ es, J. H. Matthews, C. S. Reynolds, H. R. Russell, R. N. Smith, and M. C. D. Marsh, Mon. Not. Roy. Astron. Soc. 510, 1264 (2021), arXiv:2109.03261 [astro-ph.HE]

  42. [43]

    M. C. D. Marsh, H. R. Russell, A. C. Fabian, B. P. McNamara, P. Nulsen, and C. S. Reynolds, JCAP 12, 036, arXiv:1703.07354 [hep-ph]

  43. [44]

    C. S. Reynolds, M. C. D. Marsh, H. R. Russell, A. C. Fabian, R. Smith, F. Tombesi, and S. Veilleux, Astro- phys. J. 890, 59 (2020), arXiv:1907.05475 [hep-ph]

  44. [45]

    Dessert, J

    C. Dessert, J. W. Foster, and B. R. Safdi, Phys. Rev. Lett. 125, 261102 (2020), arXiv:2008.03305 [hep-ph]

  45. [46]

    Hoof and L

    S. Hoof and L. Schulz, JCAP 03, 054, arXiv:2212.09764 [hep-ph]

  46. [47]

    Ning and B

    O. Ning and B. R. Safdi, Phys. Rev. Lett. 134, 171003 (2025), arXiv:2404.14476 [hep-ph]

  47. [48]

    T. Monz, P. Schindler, J. T. Barreiro, M. Chwalla, D. Nigg, W. A. Coish, M. Harlander, W. H¨ ansel, M. Hennrich, and R. Blatt, Phys. Rev. Lett. 106, 130506 (2011)

  48. [49]

    Y. Zhao, R. Zhang, W. Chen, X.-B. Wang, and J. Hu, npj Quantum Information 7, 24 (2021)

  49. [50]

    R. G. Barrera, G. A. Estevez, and J. Giraldo, European 8 Journal of Physics 6, 287 (1985)

  50. [51]

    Solutions to Axion Electrodynamics in Various Geometries

    J. Ouellet and Z. Bogorad, Phys. Rev. D 99, 055010 (2019), arXiv:1809.10709 [hep-ph]

  51. [52]

    Alken et al., Earth, Planets and Space 73, 49 (2021)

    P. Alken et al., Earth, Planets and Space 73, 49 (2021)

  52. [53]

    Hoburg, IEEE Transactions on Electromagnetic Compatibility 37, 574 (1995)

    J. Hoburg, IEEE Transactions on Electromagnetic Compatibility 37, 574 (1995)

  53. [54]

    S. S. Wilks, Annals Math. Statist. 9, 60 (1938)

  54. [55]

    Asymptotic formulae for likelihood-based tests of new physics

    G. Cowan, K. Cranmer, E. Gross, and O. Vitells, Eur. Phys. J. C 71, 1554 (2011), [Erratum: Eur.Phys.J.C 73, 2501 (2013)], arXiv:1007.1727 [physics.data-an]

  55. [56]

    J. W. Foster, N. L. Rodd, and B. R. Safdi, Phys. Rev. D 97, 123006 (2018), arXiv:1711.10489 [astro-ph.CO]

  56. [57]

    Badurina, A

    L. Badurina, A. Beniwal, and C. McCabe, Phys. Rev. D 108, 083016 (2023), arXiv:2306.16477 [hep-ph]

  57. [58]

    Variable Stars: which Nyquist Frequency ?

    L. Eyer and P. Bartholdi, Astron. Astrophys. Suppl. Ser. 135, 1 (1999), arXiv:astro-ph/9808176

  58. [59]

    Koen, Monthly Notices of the Royal Astronomical Society 371, 1390 (2006)

    C. Koen, Monthly Notices of the Royal Astronomical Society 371, 1390 (2006)

  59. [60]

    PUWEs/Zp2j/l3snwl3adYbzEM0A=

    L. Mandel and E. Wolf,Optical Coherence and Quantum Optics (Cambridge Univ. Pr., 1995). 9 Ultralight dark matter detection with trapped-ion interferometry Supplemental Material Leonardo Badurina, Diego Blas, John Ellis, Sebastian A. R. Ellis In the following, we provide additional details in support of the results presented in the main body of this work. ...

  60. [61]

    When evaluating the probability P[↑] with these P -representations, we find P[↑] = 1 2 1 + exp − 2(2N ηsin(βz/2))2 (1 + ¯nx + ¯ny + (¯ny − ¯nx) cosβz) cos 4N η yd sin βz x0

    Thermal Phonons If there is a non-zero thermal occupation of phonon modes ¯n in the trap, the Glauber-Sudarshan P -representation that must be integrated over is P (αi) = e−|αi|2/¯ni π¯ni , (S66) where the occupation numbers ¯ni in the x and y directions may differ. When evaluating the probability P[↑] with these P -representations, we find P[↑] = 1 2 1 +...

  61. [62]

    In this case, the trap potential is modelled as V (x, y) = κxx2 + κyy2 + κxyxy + γxx3 + γyy3 + γxyx2y + γyxxy2 + λxx4 + λyy4 +

    Semi-Classical Treatment of Imperfections As discussed in the main text, a serious possible source of error is anharmonicity in the trap potential. In this case, the trap potential is modelled as V (x, y) = κxx2 + κyy2 + κxyxy + γxx3 + γyy3 + γxyx2y + γyxxy2 + λxx4 + λyy4 + . . . , (S68) where the κij, γ i(j), λ i coefficients are treated phenomenological...