) and ⟨Y p i (t)⟩ = 0, where p, q ∈ {x, y, z} are Cartesian coordinates, δ is the Kronecker delta and F i is the deterministic force acting on the particle, originating from the gradient of potentials, given by Eq. 5, Eq. 6 and Eq. 7. We set the mass of each particle to m = 1, which in conjunction with σ = 1 and k B T = 1 defines the internal unit of time as τ = (mσ 2 /k B T ) 1/2 = 1 ≡ 1/γ = 100δt, where δt = 0.01τ is the integration time-step. Notably, the choice of the particle mass is arbitrary because it has no effect on the partition function and hence it does not affect any thermodynamic properties. It only affects dynamical properties of the system which we did not study here.