Emergence of a thermoplasmonic attractor. Trajectories x=x(τ), y=y(τ) (upper row), projections x=x(z), y=y(z) (second row) and ż=ż(x), y=y(x) (third row) at increasing interparticle energy parameter ψ0. Transition from instability to a thermoplasmonic attractor occurs near ψ0c≈19.447. The trajectories eventually reach a subset of the state space from which, over long times, they no longer escape. Bottom row: corresponding fast Fourier transforms showing a cascade of uneven spectral lines and a sharp drop in the low-frequency tail. Quantities x, y, z, τ are the same of Figure 2 and obey the dynamic system presented in Section S7 (Supporting Information) for an extended clustering dynamics, ż≡ϱ̈ being the second time derivative of the cluster density (Equations (7.3)). Other parameters and computational details are in Section S7.2 (Supporting Information).