Plot showing the perihelion and aphelion of 1000 planetesimal orbits at the last time step before they were scattered beyond the simulation domain and thus removed from the simulation, for one of the simulations in gas accretion scheme 1. The time evolution of the planet semimajor axes is included on the plot. The majority of the scattered planetesimals have an aphelion of ~50 au and a perihelion located ~1 au exterior of the planets, suggesting that they are not strongly scattered but rather are part of a scattered disk where eccentricities and aphelia increase gently with time.