  We consider a discrete particle system of two species coupled through
nonlocal interactions driven by the one-dimensional Newtonian potential, with
repulsive self-interaction and attractive cross-interaction. After providing a
suitable existence theory in a finite-dimensional framework, we explore the
behaviour of the particle system in case of collisions and analyse the
behaviour of the solutions with initial data featuring particle clusters.
Subsequently, we prove that the empirical measure associated to the particle
system converges to the unique 2-Wasserstein gradient flow solution of a system
of two partial differential equations (PDEs) with nonlocal interaction terms in
a proper measure sense. The latter result uses uniform estimates of the
$L^m$-norms of a piecewise constant reconstruction of the density using the
particle trajectories.
