Quantum sensing of high-frequency gravitational waves with ion crystals

A detection method for high-frequency gravitational waves using two-dimensional ion crystals is investigated. Gravitational waves can resonantly excite the drumhead modes of the ion crystal, particularly the parity-odd modes. In the optical dipole force protocol, entanglement between the drumhead modes and the collective spins transfers the excitation of the drumhead modes to the rotation of the total spin. Furthermore, gravitational wave detection beyond the standard quantum limit becomes possible as a squeezed spin state is generated through this entanglement. The sensitivity gets better with a larger ions crystals as well as a larger number of the ions. Future realization of large ion crystals can significantly improve the sensitivity to gravitational waves in the 10 kHz to 10 MHz region.

💡 Research Summary

The paper proposes a tabletop quantum sensor for high‑frequency gravitational waves (GWs) in the 10 kHz–10 MHz band, based on a two‑dimensional crystal of trapped 9Be⁺ ions in a Penning trap. The ion crystal supports collective axial vibrations (“drumhead modes”). Because a GW couples to the crystal through a quadrupole‑like tidal field, it selectively excites parity‑odd (antisymmetric) drumhead modes while parity‑even (symmetric) modes remain essentially unperturbed. This selection rule provides a distinctive signature that can discriminate GW signals from other perturbations.

The authors describe the ion crystal’s equilibrium configuration, the harmonic approximation for the axial motion, and the resulting normal‑mode spectrum. The center‑of‑mass (COM) mode has the highest frequency ω₁ = ω_z, while the next two nearly degenerate modes (ω₂ ≈ ω₃) are the parity‑odd modes of interest. Each ion’s axial displacement z_i is expressed in terms of phonon ladder operators a_k with eigenvectors b_{ik}.

Detection relies on an optical dipole force (ODF) generated by two off‑axis laser beams. By adjusting the beams’ polarizations, frequencies, and relative phase, a spin‑dependent force F_i = F₀ z_i σ_z^i cos(ω_ODF t) is realized. A spatially varying phase pattern (implemented with deformable mirrors) allows the ODF to couple selectively to the desired parity‑odd mode. When the ODF frequency ω_ODF is tuned to resonance with ω₂ (or ω₃), the effective interaction reduces to H_eff ≈ g √N ( a + a† ) J_z, where J_z is the collective spin operator and g is the coupling strength proportional to the ODF amplitude and the mode shape.

An external GW of amplitude α drives the phonon mode via H_int = α a + α* a†. The sensing protocol is a Ramsey‑type sequence: (i) prepare the spins in a coherent state |N/2⟩_z and the phonon mode in vacuum; (ii) apply a π/2 rotation about y; (iii) apply the ODF for a duration τ; (iv) let the system evolve freely for time T − 2τ during which the GW imprints a displacement on the phonon; (v) apply the inverse ODF (implemented by a spin flip combined with the ODF); (vi) apply a final π/2 rotation about x and measure J_z via global fluorescence. The net effect is that the spin observable acquires a phase φ = 2 g τ (T − τ) √N Im α. The expectation value ⟨J_z⟩ ≈ (N/2) sin φ ≈ N g τ (T − τ) Im α for small φ, while the variance is ≈ N/4 in the ideal case. Consequently, the single‑shot sensitivity to the GW‑induced displacement η = |Im α| is δ η = 1/