Maximum eccentricity reached after , where N ∈ {10, 100, 300} are denoted by the colors. The initial conditions are m₁ = 1M_⊙, m₂ = 0.001M_⊙, m₃ = 10⁶M_⊙, a₁ = 10⁵ ℛ₁, i_tot = 65°, e₁ = 0.001, e₂ = 0.7, and a₂ was varied to simulate different ℛ. An eccentricity excitation around ℛ = 0.65 is evident. Simulations were performed using a regularized Hermite integrator with a time-step parameter of η = 0.003 for N = 10 and N = 100 and η = 0.002 for N = 300. Each data point is an independent simulation that took a wall-clock time of ∼20 min for N = 10, ∼3 hours for N = 100, and ∼12 hours for N = 300.