Abstract
The orbital eccentricity of massive black hole binaries (MBHBs) at binary formation shapes the stochastic gravitational-wave background (GWB) detectable by pulsar timing arrays (PTAs). Previous í µí±-body simulations show large run-to-run scatter in this quantity, dominated by Poisson noise, raising the question of whether physical substructure adds genuine astrophysical stochasticity. We test this with high-resolution re-simulations of a major merger from IllustrisTNG100-1, evolved with the Griffin í µí±-body code. A no-perturber control is compared with two matched suites in which í µí±“ target = 0.1 of the primary bulge mass is redistributed into equal-mass perturbers of 10 7 M ⊙ (í µí¼‡ p ≈ 3.2 × 10 −3) and 10 8 M ⊙ (í µí¼‡ p ≈ 3.2 × 10 −2), with four realisations per scenario. The control gives í µí¼Ž í µí±’ ≈ 0.11, consistent with the Poisson noise floor at this resolution. The 10 7 M ⊙ case gives í µí¼Ž í µí±’ ≈ 0.115, indistinguishable from the control in its scatter, whereas the 10 8 M ⊙ case gives í µí¼Ž í µí±’ ≈ 0.26, a factor of 2.4 above the floor, though statistically marginal given only four realisations. This excess scatter coincides with larger event-aligned residuals in orbital energy and angular momentum and stronger torque spikes, consistent with near-impulsive perturber–MBHB encounters. In binary–single scattering theory the transition is set by the perturber–MBHB mass ratio í µí¼‡ p : the 10 7 M ⊙ case remains diffusive while the 10 8 M ⊙ case approaches the near-impulsive regime. As the expected perturber population in massive ellipticals lies mostly below this regime, perturber-driven eccentricity randomisation is unlikely to matter for GWB-relevant MBHB mergers.