Discrete-to-continuum limits of interacting particle systems in one dimension with collisions

Patrick van Meurs
2024-03-20
Abstract:We study a class of interacting particle systems in which $n$ signed particles move on the real line. At close range particles with the same sign repel and particles with opposite sign attract each other. The repulsion and attraction are described by the same singular interaction force $f$. Particles of opposite sign may collide in finite time. Upon collision, pairs of colliding particles are removed from the system. In a recent paper by Peletier, Požár and the author, one particular particle system of this type was studied; that in which $f(x) = \frac1x$ is the Coulomb force. Global well-posedness of this particle system was shown and a discrete-to-continuum limit (i.e. $n \to \infty$) to a nonlocal PDE for the signed particle density was established. Both results rely on innovative use of techniques in ODE theory and viscosity solutions. In the present paper we extend these results to a large class of particle systems in order to cover many new applications. Motivated by these applications, we consider the presence of an external force $g$, consider interaction forces $f$ with a large range of singularities and allow $f$ to scale with $n$. To handle this class of $f$ we develop several new proof techniques in addition to those used for the Coulomb force.
Analysis of PDEs
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the dynamic behavior of the signed particle system in one - dimensional space during collisions and its limit problem from discrete to continuous. Specifically, the paper studies a class of signed particle systems. These particles move on the real line and have different charge signs. When particles with opposite charges collide, they are removed from the system. The main objectives of the paper are: 1. **Well - posedness of the system**: Prove the well - posedness of this type of particle system, that is, there is a unique solution, and the solution is stable with respect to perturbations of the initial conditions. The work in this part mainly focuses on dealing with the singularity problem on the right side of the system equations during particle collisions and the situation of multiple particles colliding simultaneously. 2. **Limit from discrete to continuous**: Study the limit behavior of the particle system when the number of particles \(n\) tends to infinity. Specifically, the author hopes to prove that the empirical measure of the particle system converges to a non - local partial differential equation (PDE), which describes the evolution of the signed particle density. To achieve the above objectives, the paper introduces some new techniques and methods. In particular, when dealing with the collision problem, the "annihilation rule" is proposed, that is, when a pair of particles with opposite charges collide, they are removed from the system, and the remaining particles continue to evolve according to the original dynamic equations. In addition, the paper also considers the external potential field and different types of interaction potential energies, thus expanding the application range of previous work. ### Formula summary 1. **Basic equations of the particle system**: \[ \frac{dx_i}{dt}=\frac{1}{n}\sum_{j\neq i}b_ib_jV'_{\alpha_n}(x_i - x_j)-b_iU'(x_i),\quad t\in(0,T),\quad i = 1,\ldots,n \] where \(V_{\alpha_n}(x)=\alpha_nV(\alpha_nx)\) is the scaled form of the interaction potential energy, and \(U\) is the external potential field. 2. **Limit PDE**: - When \(\alpha_n\rightarrow\alpha\): \[ \partial_t\kappa=\partial_x(|\kappa|(V'_{\alpha}*\kappa)+|\kappa|U') \] - When \(1\ll\alpha_n\ll n\): \[ \partial_t\kappa=\partial_x(f_2(\kappa)\partial_x\kappa+|\kappa|U') \] where \(f_2(\kappa)=\|V\|_{L^1(\mathbb{R})}|\kappa|\). - When \(\frac{\alpha_n}{n}\rightarrow\beta\): \[ \partial_t\kappa=\partial_x(f_3(\kappa)\partial_x\kappa+|\kappa|U') \] where \[ f_3(\kappa)=\begin{cases} \frac{\beta^3}{\kappa^2}\Psi\left(\frac{\beta}{\kappa}\right)&\text{if }\kappa\neq0\\ 0&\text{if }\kappa = 0 \end{cases} \] and \[ \Psi(x)