Physical Review D, cilt.110, sa.4, 2024 (SCI-Expanded)
We simulated spin-spin interactions of N-bodies in linearized general relativity (GR) and linearized massive gravity of the Fierz-Pauli type (mGR). It was noted earlier that there is a discrete difference between the spin-spin interaction potential in GR and mGR for a 2-body system, akin to the van Dam-Veltman-Zakharov discontinuity in the static Newton's potential. Specifically, at large distances, GR favors antiparallel spin orientation with total spin pointing along the interaction axis, while mGR favors parallel spin orientation with total spin perpendicular to the axis between the sources. For an N-body system, a simulation in mGR hitherto has not been done and one would like to know the total spin of the system in both theories. Here we remedy this. In the simulations of GR, we observed that the total spin tends to decrease from a random initial configuration, while for mGR with a large distance, the total spin increases.