Analytic Solution for the Motion of Spinning Particles in Plane Gravitational Wave Spacetime

The interaction between spin and gravitational waves causes spinning bodies to deviate from their geodesics. In this work, we obtain the complete analytic solution of the Mathisson–Papapetrou–Dixon equations at linear order in the spin for a general plane gravitational wave with arbitrary polarization profiles. Our approach combines a parallel-transported tetrad with the translational Killing symmetries of plane wave spacetimes, yielding six conserved quantities that fully determine the momentum, spin evolution, and worldline. The resulting transverse and longitudinal motions are expressed in closed form as single integrals of the retarded time, providing a unified and model-independent framework for computing spin–curvature-induced deviations for realistic or theoretical gravitational-wave signals. This analytic solution offers a versatile tool for studying spin-dependent effects in gravitational memory, Penrose-limit geometries, and future high-precision space-based detectors.

💡 Research Summary

This paper presents a complete analytic solution of the Mathisson‑Papapetrou‑Dixon (MPD) equations for a spinning test particle moving in a general plane gravitational wave spacetime, retaining terms only to linear order in the particle’s spin (denoted O(s)). The authors adopt the Tulczyjew‑Dixon spin‑supplementary condition (SSC), which ensures that the spin magnitude is conserved and that the four‑momentum is timelike. Under this SSC, the linear‑spin approximation simplifies the MPD system: the momentum and velocity coincide (p^μ ≈ m u^μ), the two mass definitions become equal (m = M), and the spin evolves by Fermi‑Walker transport.

The background geometry is expressed in Rosen coordinates (u, v, y, z) as
ds² = 2 du dv + (1 – h₊(u)) dy² + (1 + h₊(u)) dz² – 2 h_×(u) dy dz,
where h₊(u) and h_×(u) are arbitrary functions describing the two polarization modes of the wave. This metric possesses three translational Killing vectors (∂_v, ∂_y, ∂_z). For any Killing vector ξ, the MPD dynamics conserves the quantity J_ξ = ξ·p – ½ S^{μν}∇_ν ξ_μ. Applying this to the three Killing vectors yields three conserved scalars: the longitudinal energy E = p_u and two transverse momenta J_α, J_β, which already determine the particle’s transverse momentum in the spin‑free (O(s⁰)) limit.

To capture the spin‑induced degrees of freedom, the authors construct a parallel‑transported orthonormal tetrad {e₀, e₁, e₂, e₃} along a fiducial geodesic. The first leg e₀ is the geodesic four‑velocity, while e₁ is built from the covariantly constant null vector l = ∂v. The transverse legs e₂, e₃ are initially defined in terms of the Rosen metric functions and then rotated by a single time‑dependent angle ψ(u) so that they satisfy parallel transport. The evolution equation for ψ(u) depends only on h₊, h× and their derivatives, and can be integrated once the initial angle is specified.

Decomposing the spin four‑vector in this tetrad,
s^μ = s_∥ e₁^μ + s_α e₂^μ + s_β e₃^μ,
the coefficients (s_∥, s_α, s_β) are constants of motion at O(s) because each tetrad leg is parallel‑transported. Together with the three Killing constants (E, J_α, J_β) they provide a full set of six conserved quantities that close the system.

Using these constants, the authors solve the linear system for the momentum components. The result can be written compactly as
p_u = E,
p_y = (1 + h₊) J_α + h_× J_β + Σ_A F^y_A(u) s^A,
p_z = h_× J_α + (1 – h₊) J_β + Σ_A F^z_A(u) s^A,
where the index A runs over {u, y, z} and the coefficient functions F^A_B(u) are explicit rational expressions involving h₊, h_×, their derivatives, and the scalar Δ = 1 – h₊² – h_ײ. These functions encode the entire influence of the wave profile on the particle’s momentum at linear order in spin.

Because p^μ = m u^μ at this order, the retarded time u serves as an affine parameter along the worldline (du/dτ = E/m). Consequently the transverse velocity components become simple first‑order ODEs:
dy/du = p_y/E, dz/du = p_z/E.
Since p_y(u) and p_z(u) are known functions of u, the transverse coordinates y(u) and z(u) are obtained by single quadratures. The longitudinal coordinate x(u) (or equivalently v(u)) follows from the relation between u, v, and the metric, yielding another single integral. Thus the full worldline (u, v(u), y(u), z(u)) is expressed analytically for any arbitrary wave profile.

The paper emphasizes several important implications. First, the solution is model‑independent: no specific waveform (e.g., sinusoidal, burst, or memory‑type) is assumed, making the result applicable to realistic astrophysical signals as well as theoretical constructs. Second, the six conserved quantities provide a transparent physical interpretation: E is the conserved longitudinal energy, J_α and J_β are the conserved transverse momenta, while s_∥, s_α, s_β encode the spin orientation relative to the wave propagation direction and the transverse polarization basis. Third, the analytic expressions enable precise estimates of spin‑induced corrections to particle trajectories, which could be relevant for future space‑based interferometers such as LISA, TianQin, and Taiji. For instance, a gyroscope or a small spacecraft with a known spin vector would experience a calculable deviation from geodesic motion when traversing a passing gravitational wave, potentially leaving a measurable imprint in the interferometer’s phase data.

Finally, the authors discuss connections to gravitational memory effects and Penrose limits. Since the worldline deviation integrates the wave’s history, the solution naturally captures permanent displacements (memory) that depend on the spin orientation. Moreover, because any spacetime admits a Penrose limit that locally reduces to a plane wave, the results can be used to approximate spin dynamics in more general backgrounds near null geodesics.

In summary, this work delivers the first closed‑form analytic solution of the linear‑spin MPD equations in a fully general plane gravitational wave spacetime. By exploiting the spacetime’s Killing symmetries and a parallel‑transported tetrad, the authors reduce the dynamics to six conserved quantities and express the particle’s momentum, spin, and worldline as elementary integrals of the wave’s retarded time. The framework provides a powerful tool for theoretical investigations of spin‑curvature coupling, gravitational memory, and for assessing potential observational signatures in upcoming high‑precision gravitational‑wave missions.