\section{Introduction}

Since the birth of quantum mechanics, the question of determinism has been central.
While all previous science relied on determinism, in particular through the concept of reproducibility (the same causes always produce the same effects), quantum mechanics broke with this concept through the irreducible appearance of probabilities in its postulates (the same causes do not always produce the same effects).
In particular, this happens in a measurement process where identical measurements on identical systems in identical states may lead to different results.
More generally, each time decoherence takes place, a linear superposition of states randomly and irreversibly decoheres to a statistical mixture of states.

These phenomena always happen when the quantum system interacts with a macroscopic system (e.g.\ a readable apparatus) which displays uncontrollable internal or external (``environment'') degrees of freedom. A natural explanation of the apparent randomness could then be the very existence of these degrees of freedom, which would play the role of hidden variables characterizing the apparatus.
It was shown in reference \cite{sparenberg:13} that this type of hidden variables leads to a violation of Bell's inequalities, in agreement with standard quantum mechanics and in contrast with usual hidden variables characterizing the microscopic system.

A striking test model for this idea is the so-called Mott problem \cite{mott:1929}, i.e.\ the detection of an $\alpha$-particle emitted by a radioactive source in a cloud chamber (or any modern type of track chamber).
While the quantum description of the emission is a spherical wave, the actually detected tracks are linear (but randomly distributed with an isotropic probability distribution).
In an apparatus-hidden-variable approach, the direction of a particular track would be determined by the positions of the detector-gas atoms or molecules.
These move randomly and can be considered as fixed for a particular emission, given the different orders of magnitude of the $\alpha$-particle and
gas-atom energies.

This hypothesis was tested in reference \cite{sparenberg:13} and seemed to display a reduction of the spherical wave already through the interaction with a single atom. This explanation was however too simple: the observed effect is actually due to a misuse of the Born approximation, which is non unitary and cannot be used as such to evaluate absolute probability fluxes.
Here, to avoid this problem, we perform an exact calculation instead of using the Born approximation.
We also replace the complicated 3-dimensional problem by a simpler 1-dimensional version,
as recently proposed in reference \cite{carlone:15}.
There, the spherical wave is replaced by a linear superposition of two wave packets moving left and right
and the detector is schematically represented by two Dirac combs of equally-spaced ``atoms'',
symmetric with respect to the origin (see figure \ref{model}).
Each atom is described as a spin 1/2, which is in its down state initially and can be excited to its up ``ionized'' state at later times, through interaction with the $\alpha$ wave, hence leading to a possible detection.
\begin{figure}[ht]
\centering
\scalebox{0.7}{\input{DetectionModel.tikz}}
\caption{Strongly simplified scheme of an $\alpha$-particle detection process in a cloud chamber containing a radioactive source.}
\label{model}
\end{figure}

Below we first introduce the model, discussing some of its features useful for numerical calculations.
Next we revisit the results of reference \cite{carlone:15} using a time-dependent approach,
limiting ourselves to a single spin on each side of the detector.
This provides a clear representation of the configuration-space structure of this system,
which explains the absence of multiple (left {\em and} right) measurements,
but also makes appear a limitation of the one-dimensional two-sided model.
Then, with the need to allow for arbitrary positions of the spins in mind,
we develop a time-independent approach.
There, only a one-sided model is considered, which removes the above-mentioned limitation.
We discuss the numerical results obtained up to 8 spins.
These reveal first hints that the spin positions could make measurements more or less likely.
We finish with some conclusions and perspectives for future works.

\section{The model}

In reference \cite{carlone:15}, it is proposed to describe the evolution of a quantum particle of mass $m$
(e.g.\ an $\alpha$ particle) in a quantum environment (a mesh of $N$ non-interacting non-moving spins 1/2)
within a deterministic model based on the Schrödinger equation.
The particle is thus described by a pure state, $\Psi(x, t)$, which is a spinor made up of $2^N$ components.
As already noticed by Mott \cite{mott:1929} and as shown below, this complexity of the configuration space exponentially-increasing with the number of spins, i.e.\ with the size of the ``detector'', is the key element in the appearance of a classical-like behaviour for the quantum particle.

Each of the spinor component is a standard one-dimensional wave function $\psi_c(x,t)$, where $c=0, \dots, (2^N-1)$ represents a possible spin state of the mesh.
We found that a convenient way to index the spin states is to write $c$ in binary digits, 0 representing a down state and 1 representing an up state.
For instance, for a 6-spin case with 64 possible spin configurations, $c=53=110101_\mathrm{bin}$ represents a configuration where all spins are in the up state, except for the third and fifth ones.

The $\alpha$ particle only interacts with the spins at their locations $y_n$, $n=1, \dots, N$, where it can be either elastically scattered with a strength $\beta_n$ or inelastically scattered with a strength $\gamma_n$, loosing an energy $\epsilon_n$. Hence, in all generality the Hamiltonian reads (in units $\hbar^2/2m=1$)
\begin{equation}\label{H}
H =
\bigoplus_{n=1}^N \left(\begin{array}{cc}
-\frac{\partial^2}{\partial x^2} + \beta_n\,\delta(x-y_n) + \epsilon_n & \gamma_n\,\delta(x-y_n) \\
\gamma_n\,\delta(x-y_n) & -\frac{\partial^2}{\partial x^2} + \beta_n\,\delta(x-y_n)
\end{array}\right)\:.
\end{equation}
In the simplest case, all spins being identical, one has $\gamma_n=\gamma, \beta_n=\beta, \epsilon_n=\epsilon$.

Figure \ref{Hmat} schematically represents the structure of this Hamiltonian in the 6-spin case.
Due to the Dirac-peak interactions, each channel is only coupled to 5 other channels, which leads to a hollow fractal-like structure.
\begin{figure}[ht]
\centering
\includegraphics[scale=0.6]{HamiltonianMatrix64.pdf}
\caption{The schematic sparse filling of the Hamiltonian $64\times 64$ matrix for six spins.}
\label{Hmat}
\end{figure}

\section{Time-dependent approach and limitation of the 1D approach}

Let us first revisit the results obtained in reference \cite{carlone:15} using a time-dependent approach.
Space discretization is made using the same finite-difference scheme while time discretization is made using a time-exponential approach rather than a Crank-Nicolson scheme.

\begin{figure}[ht]
\centering
\scalebox{0.9}{\input{TwoSpinsSim.tikz}}
\caption{Time evolution of a double gaussian wave packet emitted at the origin and propagating in a ``two-atom symmetric detector''.
The probabilities of finding the propagating particle in the $ij_\mathrm{bin}$ part of the state space is given by $P_{ij}$ at times $t_0<t_1<t_2$.
\label{timedep}}
\end{figure}

Figure \ref{timedep} represents the time evolution of a double gaussian wave packet, starting from the origin at time $t=0$ and evolving in a ``one-dimensional detector'' made up of two spin-1/2 atoms, symmetrically positioned on the left side and on the right side of the ``radioactive source''.
At time $t_0$, the wave packets have not reached the atoms yet and the state has only a non-vanishing component in the $00_\mathrm{bin}$ spin configuration.
At time $t_1$, the left-propagating wave packet interacts with the left atom, populating the $10_\mathrm{bin}$ component of the state.
Symmetrically, the right-propagating wave packet excites the right-hand-side atom, populating the $01_\mathrm{bin}$ component.
At time $t_2$, the excited configurations are populated with new gaussian wave packets, with a probablity of 24\% each.
The probability of non excitation is 52\%.
Remarkably, the probability of exciting both atoms together, which would correspond to a simultaneous detection on the left- and right-hand sides
is practically 0 at time $t_2$.
It is in this sense that the state-space structure explains the absence of ``left {\em and} right'' measurements,
as stressed in reference \cite{carlone:15}.
This is the one-dimensional version of Mott's original argument \cite{mott:1929} which states that in a 3-dimensional cloud chamber two atoms can be simultaneously excited only if they are aligned with the radioactive source.

In the present one-dimensional system, the situation is a bit less simple though, as guessed from figure \ref{timedep}:
at later times, both atoms will be excited, populating configuration $11_\mathrm{bin}$,
because of the reflected wave packets in configurations $10_\mathrm{bin}$ and $01_\mathrm{bin}$.
This is a limitation of the one-dimensional model, which could have been better emphasized and numerically explored in reference \cite{carlone:15}.
Such a reflected wave also appears in a 3-dimensional scattering situation but with an extremely small probability,
making it negligible:
the high energy of the incident particle leads to a very thin scattering cone
(typically $10^{-4}$ rad angular aperture for a 1-MeV $\alpha$-particle incident on a 1-\AA-radius atom).
In the following, to avoid this limitation of the one-dimensional model,
we only consider a one-sided detector.
This simplifies the problem as the number of spins is smaller;
moreover, the two-sided system can always be reconstructed by linear superposition of two one-sided systems.

\section{Time-independent approach}

Having switched to a one-sided detector, a time-independent approach can be used, meaning that the Schrödinger equation is solved for a given energy $E=k_0^2$ defined with respect to the threshold energy of the $0\dots0_\mathrm{bin}$ ground channel.
Time-dependent solutions can always be built by linear combination of these single-energy solutions,
which can be seen as a basis for time-dependent wave-packet construction.
A time-independent approach allows to solve the problem analytically (for a very small number of spins),
symbolically (for a small number of spins) or numerically with algebraic methods (potentially for an arbitrary number of spins).
Here, we have chosen the analytic approach for a single spin, which allows to define a so-called ``spin entropy''.
Next, we solve the problem numerically and provide results for up to 8 spins.
Symbolic calculations are deferred to a future work.

\subsection{1 spin}

Let us start by solving analytically the single-spin case.
The solution of the stationary Schrödinger equation is expressed in terms of 4 coefficients defined in figure \ref{wf},
which allow to calculate reflection and transmission probabilities in the elastic and inelastic channels.
The continuity conditions on the wave-function components lead to the following expressions for these coefficients,
obtained by solving a four-dimensional algebraic system,
\begin{equation}\label{eq:One_spin_solution}
\cases{R_1 = T_1 = \frac{2\mathrm{i} k_0 \gamma}{(\beta-2\mathrm{i} k_0)(\beta-2\mathrm{i} k_1) - \gamma^2},               \\
R_0 = \frac{2\mathrm{i} k_1 \beta - \beta^2 + \gamma^2}{(\beta-2\mathrm{i} k_0)(\beta-2\mathrm{i} k_1) - \gamma^2}, \quad
T_0 = \frac{-2\mathrm{i} k_0 (\beta-2\mathrm{i} k_1)}{(\beta-2\mathrm{i} k_0)(\beta-2\mathrm{i} k_1) - \gamma^2} \:,
}
\end{equation}
where $k_1=\sqrt{k_0^2-\epsilon}$ is the wave number in the excited channel.
Figure \ref{wf} also presents a typical stationary solution.
The continuity of the solution, together with the discontinuity of its first derivative due to the Dirac-peak interaction,
can be seen.

\begin{figure}[ht]
\centering
\begin{tikzpicture}
\begin{axis}[
	width=0.6\textwidth,
	height=0.2\textheight,
	xmin = -3,
	xmax = 3.25,
	ymin = -1.5,
	ymax = 3,
	every axis x label/.style={at={(current axis.right of origin)}, above=1cm},
	every axis y label/.style={at={(current axis.north west)}, right},
	axis x line=center,
	axis y line=left,
	clip=false,
	xlabel=$x$,

	ylabel=$\psi_1(x)$,
	xtick={0},
	xticklabels={$0$},
	ytick=\empty
]
\draw[-latex, very thick] (axis cs:0,0) -- (axis cs:0,2);
\addplot[black, smooth] table[x index=0, y index=2] {OneSpinSty.dat};
\draw[->, black] (axis cs:-1,1) -- (axis cs:-2,1) node[pos=0.25, above]{$R_1\,e^{-i k_1 x}$};
\draw[->, black] (axis cs:1,1) -- (axis cs:2,1) node[pos=0.25, above]{$T_1\,e^{i k_1 x}$};
\node[black] at (axis cs:-3.2,0.5) {\rotatebox{90}{Excited state}};
\end{axis}
\end{tikzpicture}

\begin{tikzpicture}
\begin{axis}[
	width=0.6\textwidth,
	height=0.2\textheight,
	xmin = -3,
	xmax = 3.25,
	ymin = -1.5,
	ymax = 3,
	every axis x label/.style={at={(current axis.right of origin)}, above=1cm},
	every axis y label/.style={at={(current axis.north west)}, right},
	x label/.style={fill=white},
	axis x line=center,
	axis y line=left,
	clip=false,
	xlabel=$x$,

	ylabel=$\psi_0(x)$,
	xtick={0},
	xticklabels={$0$},
	ytick=\empty
]
\draw[-latex, very thick] (axis cs:0,2) -- (axis cs:0,0);
\addplot[black, smooth] table[x index=0, y index=1] {OneSpinSty.dat};
\draw[->, black] (axis cs:-2,2.2) -- (axis cs:-1,2.2) node[pos=0.25, above]{$e^{i k_0 x}$};
\draw[->, black] (axis cs:-1,1.1) -- (axis cs:-2,1.1) node[pos=0.25, above]{$R_0\,e^{-i k_0 x}$};
\draw[->, black] (axis cs:1,1.1) -- (axis cs:2,1.1) node[pos=0.25, above]{$T_0\,e^{i k_0 x}$};
\node[black] at (axis cs:-3.2,0.5) {\rotatebox{90}{Ground state}};
\end{axis}
\end{tikzpicture}
\caption{Stationary solution corresponding to a plane wave incident on a single spin 1/2 in its ground state. \label{wf}}
\end{figure}

Such solutions can be advantageously characterized by a ``spin entropy'', summarizing the splitting of the solution among the different spin states. For the single-spin case, this entropy reads
\begin{equation}
 S = -P_\mathrm{exc} \ln P_\mathrm{exc} - P_\mathrm{gnd} \ln P_\mathrm{gnd},
\end{equation}
where the excitation probability $P_\mathrm{exc}$ reads
\begin{equation}
 P_\mathrm{exc} = \frac{k_1}{k_0} \left(|R_1|^2+|T_1|^2 \right)
 = \frac{8 k_0 k_1 \gamma^2}{4\beta^2(k_0+k_1)^2 + (4 k_0 k_1-\beta^2 +\gamma^2)^2}
\end{equation}
and the ground-state probability $P_\mathrm{gnd} = |R_0|^2+|T_0|^2 = 1-P_\mathrm{exc}$.

The most interesting dependence of this excitation probability is on the coupling strength $\gamma$,
as seen in figure \ref{Pexc};
it displays a maximum value for $\gamma=\gamma_\mathrm{max}$
which can be interpreted as the optimal coupling strength for particle detection.
When $\epsilon=\beta=0$, the excitation probability simplifies to
\begin{equation}
 P_\mathrm{exc} = |R_1|^2+|T_1|^2 = \frac{8 k_0^2 \gamma^2}{(4 k_0^2+\gamma^2)^2} \quad (\epsilon=\beta=0),
\end{equation}
which leads to a maximal spin entropy of ln 2 for $\gamma_\mathrm{max}=2k_0$.
In the single-spin case, the maximum detection probability is thus 50\%.
For other values of $\epsilon$ and $\beta$, this probability is decreased but the general trend of the $\gamma$ dependence remains the same.

\begin{figure}
\centering
\begin{tikzpicture}
\begin{axis}[
	width=0.6\textwidth,
	height=0.2\textheight,
	xmin = 0,
	xmax = 20.5,
	ymin = 0,
	ymax = 1.1,
	every axis x label/.style={at={(current axis.right of origin)}, right},
	every axis y label/.style={at={(current axis.north west)}, right},
	axis x line=bottom,
	axis y line=left,
	xlabel=$\gamma$,
	ylabel=$P_\mathrm{exc}$,
	xtick={0, 5.3665},
	xticklabels={$0$, $\gamma_\mathrm{max}$},
	ytick={0.5,1},
	yticklabels={$P_\mathrm{max}\leq\frac{1}{2}$,$1$},
]
\draw[dashed, gray] (axis cs:0,0.5) -- (axis cs:5.3665,0.5) -- (axis cs:5.3665,0);
\addplot[black, smooth] table[x index=0, y index=1] {OneSpinExcProba.dat};
\fill[black] (axis cs:5.3665,0.5) circle[radius=1pt] node[right,above]{max};
\node[black] at (axis cs:1.5,0.15) {$\propto \gamma^2$};
\node[black] at (axis cs:18,0.15) {$\propto 1/\gamma^2$};
\end{axis}
\end{tikzpicture}
\caption{Excitation probability as a function of the coupling strength $\gamma$ for the single-spin case. \label{Pexc}}
\end{figure}

\subsection{$N$ spins}

Let us finally proceed to the general case of $N$ spins.
As a first approach we developed a purely numerical algorithm based on the construction and resolution of a $2N 2^{N}$-dimensional linear-algebra system obtained by combining the continuity equations satisfied by all the wave-function coefficients
$R_{ci}, T_{ci}$ with $c\in[0,2^N-1]$, $i\in[1, N]$.
No particular attempt was made to exploit the hollow structure of the Hamiltonian shown in figure \ref{Hmat}.
Hence, the algorithm is limited to $N<9$ on a standard personal computer.
In future works, we plan to improve this algorithm and to solve the problem both symbolically and numerically.

These preliminary results already prove very interesting.
To discuss them, we introduce the $\bw$-spin excitation probabilities $P(\bw)$,
where $\bw$ is the Hamming weight, i.e.\ the number of 1-bits in a binary sequence representing an integer.
Following reference \cite{carlone:15},
we define the ground-state, one-spin excitation and track-detection probabilities in terms of these excitation probabilities as
\begin{equation}
\underbrace{P(\bw=0)}_{P_\mathrm{gnd}} +
\underbrace{P(\bw=1)}_{P_\mathrm{OS}} +
\underbrace{P(\bw\geq 2)}_{P_\mathrm{trk}} = 1.
\end{equation}
Table \ref{spin-exc} shows these probabilities for a regularly-spaced $N$-spin network.
Their general behavior can be interpreted through an entropic effect:
the no-spin excitation probability is close to 50\% for small $N$ and regularly decreases with increasing $N$,
as it becomes a more and more exceptional configuration among the $2^N$ possible configurations.
This shows that, contrary to the single-spin case, the probability of no excitation in the detector can become much smaller than 50\%.
The one-spin excitation, corresponding to $N$ possible configurations, stays roughly constant while the track-detection probability,
corresponding to $2^N-N-1$ possible configurations, regularly increases.
\begin{table}
\centerline{
\begin{tabular}{c|cccc}
$k_0=\pi, \gamma= 3 $          & $N=2$     &  $N=4$      & $N=6$       & $N=8$       \\
\hline
$P_\mathrm{gnd}$  & $0.60514$  & $0.41826$    & $0.24661$      & $0.15813$   \\
$P_\mathrm{OS}$ & $0.28992$  & $0.36631$    & $0.39652$      & $0.37044$   \\
$P_\mathrm{trk}$   & $0.10494$  & $0.21543$   & $0.35687$     & $0.47143$
\end{tabular} \hspace{1cm}}
\caption{\label{spin-exc}
Probabilities for no-spin ($P_\mathrm{gnd}$), one-spin ($P_\mathrm{OS}$) or more-than-one-spin ($P_\mathrm{trk}$) excitations
for a particle of energy $\pi^2$ incident on a regularly-spaced $N$-spin network of step 0.1
for $\epsilon=0.01$, $\beta=0.5$ and $\gamma=3$.}
\end{table}

However, these simple trends are not universal, as table \ref{Ptrk} shows.
There it is seen that the track-detection probability can strongly depend on the number of spins.
This seems to be the case in particular for large energies and
when the single-spin excitation is maximal, i.e.\ for $\gamma=\gamma_\mathrm{max}$
(see the surprisingly small value obtained for $k_0=400/3$, $\gamma=800/3$, $N=4$).
This suggests that a track detection could be a resonance phenomenon,
occurring for precise matches of the particle energy, atom positions and coupling strength.
This hypothesis will be systematically explored in future works.
\begin{table}
 \centerline{
\begin{tabular}{c|cccc}
 $P_\mathrm{trk}, k_0=400/3$ & $\gamma = 50$ & $\gamma = 100$ & $\gamma = 150$ & $\gamma_{\max} = 800/3$ \\
\hline
$N = 2$ & $0.0074691$   & $0.063126$     & $0.14733$      & $0.24907$            \\
$N = 3$ & $0.015572$    & $0.12502$      & $0.21374$      & $0.2812$              \\
$N = 4$ & $0.032093$    & $0.25924$      & $0.39707$      & $0.066758$           \\
$N = 6$ & $0.07459$     & $0.40374$      & $0.51122$      & $0.22914$             \\
$N = 7$ & $0.099126$    & $0.47012$      & $0.40409$      & $0.37756$             \\
$N = 8$ & $0.12356$     & $0.49174$      & $0.50721$      & $0.32307$
\end{tabular}}
\caption{Probabilities of ``track detection'', i.e.\ more-than-one-spin excitation,
for a particle of energy $k_0^2$ incident on a regularly-spaced $N$-spin network of step $0.05/N$
for $\epsilon=0.04$, $\beta=10^{-4}$ and different values of $\gamma$. \label{Ptrk}}
\end{table}

\section{Conclusions and perspectives}

In conclusion, the present paper revisits Mott's seminal measurement problem \cite{mott:1929} following the one-dimensional model proposed in reference \cite{carlone:15}.
Our approach differs from this last reference in two respects:
first, our aim is more ambitious as we aim at answering the question of why a particular track is chosen instead of just explaining why no multiple track is observed.
For that, we develop a model allowing to arbitrarily choose the atom positions,
as it is envisaged that these positions could be the apparatus hidden variables determining the measurement outcome.
Second, to avoid artificial reflections due to the one-dimensional character,
we limit ourselves to a one-sided detector.
This divides the number of spins by two and allows to use a time-independent approach which is more promising for large-scale calculations.

In the single-spin case, an analytical result was obtained, which shows that the most important parameter of the model is $\gamma$,
the coupling strength between the different spin states. The other parameters, $\epsilon$ and $\beta$, which represent the energy loss during the spin excitation and the elastic scattering on the spin, can probably be chosen to zero without much loss in the physical content of the model.
This conclusion should hold for the $N$-spin case, except if one intends to model the energy loss of the particle in the detector,
e.g.\ to calculate the detector stopping power.

For the $N$-spin case, first numerical results were obtained for an equally-spaced mesh.
They were limited to 8 spins for computational-cost reasons but were however sufficient to reveal that track-formation probabilities indeed depend on the detector microscopic state (the number of atoms in the present case).
We suspect this behavior to be due to a kind of resonance phenomenon \cite{barra:99},
allowing the incident wave to be transmitted through the spin network only for particular energy/position combinations.
In the near future, we plan to further study this phenomenon, combining symbolic and numeric calculations,
for randomly-spaced spin meshes.

On a longer term, a full 3-dimensional model should be set up.
Developing such a model with realistic orders of magnitudes for elastic-scattering and ionization probabilities could be very useful not only in the context of the quantum-measurement problem in a cloud-chamber apparatus.
It could provide a microscopic approach to existing macroscopic effective and non-linear models,
like the stopping-power formula by Bethe \cite{bethe:1930},
and replace these models when they are less reliable, in particular in low-energy situations.
On an even longer term, combining these microscopic and macroscopic models could be applied to the development of new track detectors,
like the active targets used in nuclear physics \cite{beceiro-novo:15} or the huge CLOUD cloud chamber at CERN \cite{kirkby:16},
which is dedicated to the study of cloud formation due to high-energy radiation present in the high atmosphere.

Regarding the measurement problem, our approach is similar and complimentary to the very complete detector description performed in reference \cite{allahverdyan:13}, where the relevant microscopic degrees of freedom could be explicitly taken into account.
There, the full thermodynamical process taking place during a measurement could be tackled, in addition to the preliminary switch from linear superposition to statistical mixture.
We hope a similar sophistication could be reached someday for track detectors too.
We suspect tracing out the microscopic degrees of freedom, i.e.\ the spin states,
could lead to the increase of the ``spin entropy'' for the scattered wave.
We conjecture the measured track could correspond to a minimal spin-entropy increase,
as the corresponding wave would be far reaching, with a higher probability of leading to a track formation.
Other waves would be stopped by an Anderson-localization-like phenomenon \cite{anderson:1958}
due to a spreading of the wave over the spin configuration space.
Note that the entropy of the total system, that would take into account the interaction between the atoms and the full vapor-condensation process,
would still be maximal, contrary to the spin entropy, since the track formation corresponds to a phase change with a strong entropy increase.

Let us finally stress that, once again, our philosophical starting point is more ambitious (or just more naive!) than the ones adopted in references \cite{carlone:15} or \cite{allahverdyan:13}, where no attempt is made at finding a deterministic explanation for a particular measurement outcome.
We consider that as long as this possibility has not been totally and definitely ruled out, it is an option very worth exploring as it could indeed bring quantum mechanics back into the world of deterministic science.

\ack{This text presents research results of the IAP program P7/12 initiated by the Belgian-state Federal Services for Scientific, Technical, and Cultural Affairs.
We thank Ruben Ceulemans (KU Leuven) for the check of our numerical results.}

\section*{References}

\bibliographystyle{jphysicsB}


\section{Introduction}

\section{Introduction}
Mathematical description of dynamics of aggregating, swarming or flocking particles plays a significant role in modelling of various physical, biological and sociological phenomena.
 From the mathematical point of view such models take the form of transport-type PDE's or systems of ODE's, where the individuals are subjected to a force generated by their nonlocal interactions. Depending on the nature of the interactions, the individuals' behavior may differ for instance they can aggregate, align their velocities or disperse. There is a wide range of applications of such models that include such seemingly unrelated phenomena like distribution of goods, reaching a consensus among individuals or emergence of common languages in primitive cultures (see \cite{lang, cons, cons2, goods}). Among such models, our main subject of interest is the Cucker--Smale (in short, CS) flocking model introduced in \cite{cuc1} with the purpose of describing motion of self-propelled agents with tendency to flock. The CS dynamical system reads as follows:
\begin{align}
\begin{aligned}\label{cspart}
\left\{
\begin{array}{lcr}
\displaystyle \dot{x}_i=v_i, \quad i=1,\cdots, N,  \quad t > 0,\\[2mm]
\displaystyle \dot{v}_i=\frac{1}{N}\sum_{j=1}^N\psi(|x_i-x_j|)(v_j-v_i),
\end{array}
\right.
\end{aligned}
\end{align}
subject to the initial data
\bq\label{ini_cspart}
(x_i,v_i)(0) =: (x_{i0}, v_{i0}), \quad i =1,\cdots, N
\eq

Here $N$ is the number of particles, while $x_i(t)$ and $v_i(t)$ denote the position and velocity of $i$-th particle at the time $t$, respectively. The function $\psi$ is referred to as
{\it the communication weight} and it is nonnegative and nonincreasing. The state of the art for CS model (and more generally, for kinetic models of interacting particles) is rich; it includes
results in various directions, such as time asymptotics (see e.g. \cite{hakalaru, car}), pattern formation (see e.g. \cite{hajekaka, top}), models with additional deterministic (see e.g. \cite{car3, hahaki})
 or stochastic (see e.g. \cite{cuc4, aha2, choi, dua1, halele}) forces, aggregation with leaders (see e.g. \cite{cuc3, shen}), analysis of CS model with singular weight $\psi$ (see e.g. \cite{haliu, ahn1, jpe, jps}),
 passage from particle to the kinetic description for CS and similar models (see e.g. \cite{can, CH16, haliu, hatad, rec, deg1, deg3, mp}) and, particularly interesting from the point of view of this paper, collision avoidance
(see e.g. \cite{cuc2, park}). We refer to \cite{CCP, CHL} for recent surveys.

In \cite{haliu}, the CS model with a singular communication weight of the form
\begin{align}\label{psi}
\psi(s)=s^{-\alpha}\ \ \ {\rm for}\ \ \ \alpha>0
\end{align}
is considered and the time asymptotics behavior of solutions for the particle system \eqref{cspart} are provided. The question of existence, uniqueness and regularity of solutions was being consecutively answered after that work. In \cite{ahn1}, the particles' tendency to avoid collision is proved for $\alpha\geq 1$. To be more precise, the authors provide a set of initial configurations leading to no finite-time collision between particles and, as a consequence, the existence, uniqueness, and regularity of solutions to \eqref{cspart} are established for such initial data. On the other hand, existence, regularity, uniqueness of solutions to \eqref{cspart} in case of $\alpha\in(0,1)$ are investigated in \cite{jpe} and \cite{jps}. These results established a dichotomy: for $\alpha\in(0,1)$ the particles can collide and stick together and existence of solutions is obtained thanks to the low singularity of $\psi$ (particularly due to integrability of $\psi$ at $0$); on the other hand for $\alpha\in[1,\infty)$, the particles exhibit a tendency to avoid collisions.

The issue of collision-avoiding plays a significant role in quantitative analysis and applications of the model. Since system \eqref{cspart} with weight \eqref{psi} is singular only at times at which particles collide, knowledge that the particles do not collide enables us to immediately deduce existence, uniqueness and regularity of solutions. From the point of view of applications collision-avoiding is important in any situation when the described phenomena involves flocking of agents that naturally avoid collisions such as birds, fish or robots. In fact many models of flocking, such as in \cite{can}, that include effects of aggregation and alignment also include a repulsion effect to ensure the lack of collisions between particles.

The goal of this paper is to refine the results of \cite{ahn1} proving that for $\alpha\in[1,\infty)$ the CS particles indeed cannot collide (regardless of the initial configuration, as long as $x_{i0}\neq x_{j0}$ for $i\neq j$) and thus the solutions to \eqref{cspart} globally exist in time and are smooth and unique.

Further, in case of $\alpha\geq 2$ we improve our results by introducing a uniform with respect to the number of particles non-collisioness estimate. We prove the estimate for the following generalization of the CS particle system. Given $\delta \geq 0$ we introduce
\begin{align}
\begin{aligned}\label{dcspart}
\left\{
\begin{array}{lcr}
\displaystyle \dot{x}_i=v_i, \quad i=1,\cdots, N,  \quad t > 0,\\[2mm]
\displaystyle \dot{v}_i=\frac{1}{N}\sum_{j=1}^N\psi(|x_i-x_j|-\delta)(v_j-v_i),
\end{array}
\right.
\end{aligned}
\end{align}
subject to the initial data \eqref{ini_cspart} with $|x_{i0}-x_{j0}|>\delta$ for $i\neq j$.
In case of $\alpha>2$, the uniformly estimated quantity that determines $\delta$-distance between particles is given by
\begin{equation}\label{es}
\frac{1}{N^2}\sum_{i,j=1}^N(|x_i-x_j|-\delta)^{2-\alpha}\leq C \quad \mbox{for}\quad \alpha>2.
\end{equation}
Note that if $\delta = 0$ then \eqref{dcspart} reduces to \eqref{cspart} and \eqref{es} serves as a uniform estimate of non-collisioness. Similarly for $\delta>0$ it implies that
the distances between particles are always greater than $\delta$ (provided that $|x_{i0}-x_{j0}|>\delta$ for $i\neq j$). This result can be viewed in the following way. We expand the set of singular points of $\psi$ from $\{0\}$ to $[0,\delta]$ and prove that the property that the solutions do not enter the singular set (in the sense that their distances do not enter this set) is preserved. It is also interesting from the point of view of applications, since it shows that the alignment kernel alone can be used to establish a minimal distance between particles (for which usually a separate repulsion kernel is utilized).

Another motivation of \eqref{es} comes from the fact that regularity of solutions to \eqref{cspart} can be controlled by the distances between particles. Thus estimates similar to \eqref{es} can be used in the passage from the particle system to the kinetic CS equation by mean-field limit as done for example in \cite{mp}.

The paper is organized as follows. In Section \ref{noncoll} we present the proof of non-collisioness of solutions to \eqref{cspart} with $\alpha\geq 1$ from which we deduce the global existence, uniqueness and regularity of solutions to the particle system \eqref{cspart}. In Section \ref{estim}, we show the $N$-independent non-collisioness estimate for $\alpha\geq 2$.

\section{A global existence theory for the particle system}\label{noncoll}
In this section, we analyze the dynamics of the CS model with a singular weight given by \eqref{psi}. Our goal is to prove that for $\alpha\geq 1$ the particle system
\eqref{cspart}-\eqref{ini_cspart} does not allow any collisions between particles and as a consequence admits a unique smooth solution globally in time. For this, we
follow the ideas of \cite{haliu} establishing a system of locally dissipative differential inequalities (SDDI) for the quantities that control collisions between the particles.

From now on, $x=(x_1,\cdots,x_N)$, where $x_i=(x_{i,1},\cdots,x_{i,d})$ with $d \geq 1$ denotes the position of the particles, while $v = \dot{x}$ is their velocity.
Throughout the paper $C$ denotes a generic positive constant that may change from line to line even in the same inequality.

Let us be more precise about the meaning of collision between particles.

\begin{defi}
We say that the $i$-th and $j$-th particles collide at $t_0$ if and only if $x_i(t_0)=x_j(t_0)$ and we say that they stick together if they collide and $v_i(t_0)=v_j(t_0)$.
\end{defi}

The following proposition contains useful properties of solutions to the CS model showing the uniform bound estimate of the velocity in time. We refer to \cite{carchoha, jpe} for its proof.
\begin{prop}\label{bounds}
Let $(x,v)$ be a solution to (\ref{cspart}) on an arbitrary finite interval $[0,T]$, with singular communication weight given by (\ref{psi}). Then we have
\bq\label{vM}
\sup_{1 \leq i \leq N}|v_i(t)|\leq \sup_{1 \leq i \leq N}|v_{i0}| \quad \mbox{and} \quad \sup_{1 \leq i \leq N}|x_i(t)|\leq \sup_{1 \leq i \leq N}|x_{i0}| +
\sup_{1 \leq i \leq N}|v_{i0}|\,t \nonumber \quad \mbox{for } t \geq 0.
\eq
\end{prop}
For notational simplicity, we set $M := \sup_{1 \leq i \leq N}|v_{i0}|$ and $R=R(T):= \sup_{1 \leq i \leq N}|x_{i0}|$ + TM.

In order to introduce a general strategy of the proofs, we begin with a simplified, one--dimensional case that serves as a good starting point and introduces the main ideas
behind our argumentation in higher dimensions.
\begin{prop}\label{d1}
Let $d=1$ and $\alpha\geq 1$. Suppose that the initial data $(x_0,v_0)$ satisfy
\[
x_{i0}\neq x_{j0} \quad \mbox{for } 1 \leq i \neq j \leq N.
\]
Then the system \eqref{cspart}-\eqref{ini_cspart} admits a unique smooth solution. Moreover, the trajectories of this solution do not collide in finite time, i.e.,
\[
x_i(t)\neq x_j(t) \quad \mbox{for } 1 \leq i \neq j \leq N, \quad t \geq 0.
\]
\end{prop}
\begin{proof}
Suppose  $\alpha\geq 1$ and that we have an arbitrary non--collision initial data $(x_0,v_0)$. Then since $\psi(|x_i-x_j|)$ is regular as long as $x_i\neq x_j$,
then there exists a unique smooth solution to \eqref{cspart} on the time interval $[0,t_0)$, where $t_0$ is a supposed time of the first collision between any number of particles. Suppose that $l$th particle collides with some other particles at $t_0$ and denote by $[l]$ the set of all indices $i$ such that $i$th particle collides with $l$th particle at $t_0$.
We will show that $t_0 < \infty$ actually does not exist, i.e., $t_0 = \infty$. Since $x(t)$ is Lipschitz continuous on $[0,t_0)$, it can be continuously extended to $t_0$, which means that $|x_i(t)-x_j(t)|\to 0$ as $t\to t_0$ if $i,j\in[l]$. Suppose without a loss of generality that $x_i(t)$ with $i\in[l]$ is the furthest to the right particle
from $[l]$, i.e., we have $x_j(t) < x_i(t)$ for all $j\in[l], j\neq i$ on $[0,t_0)$. For $t\in[0,t_0)$, let
\begin{align*}
r(t):= \sum_{j\in[l]}(v_i(t)-v_j(t)).
\end{align*}
Then it follows from $(\ref{cspart})_2$ that
\begin{align}\label{d1eq1}
\begin{aligned}
\frac{d}{dt}r
&= \frac{1}{N}\sum_{j\in[l]}\sum_{k=1}^N\psi(|x_k-x_i|)(v_k-v_i)
-\frac{1}{N}\sum_{j\in[l]}\sum_{k=1}^N\psi(|x_k-x_j|)(v_k-v_j)\\
&= \frac{1}{N}\sum_{j,k\in[l]}\psi(|x_k-x_i|)(v_k-v_i)
-\frac{1}{N}\sum_{j,k\in[l]}\psi(|x_k-x_j|)(v_k-v_j)\\
&\quad +\frac{1}{N}\sum_{\stackrel{j\in[l]}{k\notin[l]}} \left(\psi(|x_k-x_i|)(v_k-v_i)-\psi(|x_k-x_j|)(v_k-v_j)\right)\\
&=: I_1 + I_2 + I_3,
\end{aligned}
\end{align}
where by substituting indices $j$ and $k$ from the set $[l]$, we obtain
\begin{align*}
I_2 = \frac{1}{N}\sum_{j,k\in[l]}\psi(|x_k-x_j|)(v_j-v_k) = -I_2=0
\end{align*}
and $I_3$ is integrable on $[0,t_0)$. On the other hand, since $x_i(t)>x_k(t)$ for $k\in[l]$, we find
\begin{align*}
\sum_{k\in[l]}\left(\Psi(|x_k-x_i|)\right)' = \sum_{k\in[l]}\left(\Psi(x_i-x_k)\right)' = \sum_{k\in[l]}\psi(x_i-x_k)(v_i-v_k) = I_1.
\end{align*}
Here $\Psi$ is the primitive of $\psi$, i.e.,
\bq\label{def_pri}
\Psi(s) = \left\{ \begin{array}{ll}
 \ln(s) & \textrm{if $\alpha = 1$,}\\[1mm]
\displaystyle \frac{1}{1-\alpha}s^{1-\alpha} & \textrm{if $\alpha > 1$}.
  \end{array} \right.
\eq
Integrating \eqref{d1eq1} over $[0,t]$ with $t < t_0$ it yields
\[
r(t)-r(0)  = \sum_{k\in[l]}\Psi(|x_{0,k}-x_{0,i}|)-\sum_{k\in[l]}\Psi(|x_k(t)-x_i(t)|) + \int_0^t I_3\,ds,
\]
and by the bounds from Proposition \ref{bounds}, this implies that the function
\[
t\mapsto \sum_{k\in[l]}\Psi(|x_k(t)-x_i(t)|)
\]
is bounded on the time interval $[0,t_0)$. However, this cannot be the case since $|x_k(t)-x_i(t)|\to 0$ as $t\to t_0$ for all $k\in[l]$ and $\Psi$ has a singularity at $0$.
Thus this contradicts the supposed existence of a finite time of collision. Hence the unique smooth solution that was previously assumed to exist on the time interval $[0,t_0)$
 can be actually prolonged up to an arbitrary finite time $T$. This completes the proof.
\end{proof}

We are now in a position to provide the proof of non-collisioness for the multi-dimensional case.

\begin{theo}\label{thm_main1} Let $d \geq 1$ and $\alpha \geq 1$. Suppose that the initial data $(x_0,v_0)$ are non-collisional, i.e. they satisfy
\[
x_{i0}\neq x_{j0} \quad \mbox{for } 1 \leq i \neq j \leq N.
\]
Then the system \eqref{cspart}-\eqref{ini_cspart} admits a unique smooth solution. Moreover, the trajectories of this solution are also non-collisional, i.e.,
\[
x_i(t)\neq x_j(t) \quad \mbox{for } 1 \leq i \neq j \leq N, \quad t \geq 0.
\]

\end{theo}

\begin{proof}[Proof of Theorem \ref{thm_main1}]
At the very beginning of the proof, let us reformulate the thesis in a more suitable way. Let us fix $\alpha\geq 1$ and $T>0$. We need to prove that on $[0,T]$
there exists a unique solution to (\ref{cspart}) and that collisions of particles are impossible. However, since at $t=0$ the particles have distinct
positions and the communication weight $\psi$ is singular only in a neighborhood of $0$, local existence of a unique smooth solution is standard. In fact, there
are two possibilities: the particles do not collide in $[0,T]$ and the local existence can be extended up to $[0,T]$ or there exists $t_0\in(0,T]$, the first time of
 collision of any particles and we only know that the solution exists and is unique and smooth on $[0,t_0)$. Thus let us suppose that such $t_0$ exists. Then, by its definition,
 there exists an index $l=1,...,N$ such that the $l$th particle collides with some other particles. Let us denote by $[l]$ the set of all those indices $j \in \{1,\cdots, N\}$ that
the $j$th particle collides with $l$th particle, i.e.,
\begin{align}
&|x_l(t)-x_j(t)|\to 0\ \ \ {\rm as}\ \ \  t\to t_0\ \ \ {\rm for\ all}\ \ \ j\in[l],\nonumber\\
&|x_l(t)-x_j(t)|\geq\delta>0\ \ \ {\rm in}\ \ \ [0,t_0)\ \ \ {\rm for\ all}\ \ \ j\notin[l]\ \ \ {\rm and\ some}\
 \ \ \delta>0.\label{delta}
\end{align}
Thus by the definition of $t_0$ there exists at least one set $[l]$ such that $\lt|[l]\rt| > 1$. Let us fix one of such sets denoting it simply by $[l]$. Let us also  denote
\begin{align*}
\|x\|_{[l]}(t):= \sqrt{\sum_{i,j\in[l]}|x_i(t)-x_j(t)|^2} \quad \mbox{and} \quad \|v\|_{[l]}(t):= \sqrt{\sum_{i,j\in[l]}|v_i(t)-v_j(t)|^2},
\end{align*}
where the sum is taken over all $i,j\in[l]$. Then we get
\begin{align}\label{contr}
\|x\|_{[l]}(t)\to 0\ \ \  {\rm as}\ \ \  t\nearrow t_0.
\end{align}
In particular if we show that it is impossible for $\|x\|_{[l]}$ to converge to $0$ as $t\nearrow t_0$, then it will contradict the assumption that $t_0$ is the first time of collision of any particles. We have
\begin{align*}
\frac{d}{dt}\|x\|_{[l]}^2 = 2\sum_{i,j\in[l]}(x_i-x_j)\cdot(v_i-v_j)\leq 2\sqrt{\sum_{i,j\in[l]}|x_i-x_j|^2}\sqrt{\sum_{i,j\in[l]}|v_i-v_j|^2}= 2\|x\|_{[l]}\|v\|_{[l]},
\end{align*}
which implies that
\begin{align}\label{rprime}
\left|\frac{d}{dt}\|x\|_{[l]}\right|\leq \|v\|_{[l]}.
\end{align}
On the other hand, we find from $(\ref{cspart})_2$ that
\begin{align}\label{dv}
\frac{d}{dt}\|v\|_{[l]}^2 &= 2\sum_{i,j\in[l]}(v_i-v_j)\cdot\left[\frac{1}{N} \sum_{k=1}^N\psi(|x_k-x_i|)(v_k-v_i)- \frac{1}{N}\sum_{k=1}^N\psi(|x_k-x_j|)(v_k-v_j)\right]\\
&=\frac{2}{N}\left(\sum_{i,j,k\in[l]}+\sum_{\stackrel{i,j\in[l]}{k\notin[l]}}\right) \left[\psi(|x_i-x_k|)(v_i-v_j)\cdot(v_k-v_i)-\psi(|x_j-x_k|)(v_i-v_j)\cdot(v_k-v_j)\right]\nonumber\\
&=: J_1 + J_2.\nonumber
\end{align}
The estimate of $J_1$ follows by its antisymmetry, which we explain below.
By substituting indices $k$ and $i$ in the first term of $J_1$, we obtain
\begin{align}\label{symm}
&\frac{2}{N}\sum_{i,j,k\in[l]}\psi(|x_k-x_i|)(v_i-v_j)\cdot(v_k-v_i) \\
&\quad =\frac{1}{N}\sum_{i,j,k\in[l]}\psi(|x_k-x_i|)(v_i-v_j)\cdot(v_k-v_i) +\frac{1}{N}\sum_{i,j,k\in[l]}\psi(|x_k-x_i|)(v_k-v_j)\cdot(v_i-v_k)\nonumber\\
&\quad = \frac{1}{N}\sum_{i,j,k\in[l]}\psi(|x_k-x_i|)(v_i-v_k)\cdot(v_k-v_i)\nonumber\\
&\quad =-\frac{|[l]|}{N}\sum_{i,j\in[l]}\psi(|x_i-x_j|)|v_i-v_j|^2,\nonumber
\end{align}
and similarly
\begin{align*}
- \frac{2}{N}\sum_{i,j,k\in[l]}^N\psi(|x_k-x_j|)(v_i-v_j)\cdot(v_k-v_j) = -\frac{|[l]|}{N}\sum_{i,j\in[l]}\psi(|x_i-x_j|)|v_i-v_j|^2.
\end{align*}
This yields
\begin{align*}
J_1 \leq -\frac{2|[l]|}{N}\sum_{i,j\in[l]}\psi(|x_i-x_j|)|v_i-v_j|^2.
\end{align*}
Since we have $|x_i-x_j|\leq  \|x\|_{[l]}$ for all $i,j\in[l]$, by monotonicity of $\psi$ we obtain
\begin{align*}
J_1 \leq-2c_0\psi(\|x\|_{[l]})\sum_{i,j\in[l]}|v_i-v_j|^2 = -2c_0\psi(\|x\|_{[l]})\|v\|_{[l]}^2,
\end{align*}
where $c_0 := |[l]|/N$. The estimate of $J_2$ follows by the fact that $|x_k-x_i|$ is separated from $0$ for $i\in[l]$ and $k\notin[l]$, which is written explicitly in (\ref{delta}). We have
\begin{align*}
J_2 &= \frac{2}{N}\sum_{\stackrel{i,j\in[l]}{k\notin[l]}}\psi(|x_i-x_k|)(v_i-v_j)\cdot(v_j-v_i) + \frac{2}{N}\sum_{\stackrel{i,j\in[l]}{k\notin[l]}}\left(\psi(|x_i-x_k|)-\psi(|x_j-x_k|)\right)(v_i-v_j)\cdot(v_k-v_j) \\
&= \underbrace{-\frac{2}{N}\sum_{\stackrel{i,j\in[l]}{k\notin[l]}} \psi(|x_i-x_k|)|v_i-v_j|^2}_{\leq 0} + \frac{2}{N}\sum_{\stackrel{i,j\in[l]}{k\notin[l]}}\left(\psi(|x_i-x_k|)-\psi(|x_j-x_k|)\right)(v_i-v_j)\cdot(v_k-v_j) \\
&\leq \frac{2L(\delta)}{N}\sum_{\stackrel{i,j\in[l]}{k\notin[l]}}|(v_i-v_j)\cdot(v_k-v_j)||x_i-x_j|,
\end{align*}
where $L(\delta)$ is the Lipschitz constant of $\psi$ in the interval $(\delta,\infty)$. Moreover by Proposition \ref{bounds} velocity $|v_k-v_j|$ is bounded by $2M$ and ultimately by H\"older's inequality
\[
J_2 \leq \frac{4ML(\delta)}{N}\sum_{\stackrel{i,j\in[l]}{k\notin[l]}}|v_i - v_j||x_i - x_j| = \frac{4ML(\delta)(N - |[l]|)}{N}\sum_{i,j \in [l]}|v_i - v_j||x_i - x_j| \leq 2c_1 \|v\|_{[l]}\|x\|_{[l]},
\]
where $c_1$ is a positive constant given by
\[
c_1 := \frac{2ML(\delta)(N- |[l]|)}{N}.
\]
Combining all of the above estimates, we have
\[
\frac{d}{dt}\|v\|_{[l]}^2\leq -2c_0\psi\lt(\|x\|_{[l]})\rt)\|v\|_{[l]}^2 + 2c_1\|v\|_{[l]}\|x\|_{[l]}.
\]
Note that if $\|v\|_{[l]} \neq 0$, then we get $\frac{d}{dt}\|v\|_{[l]}^2 = 2\|v\|_{[l]}\frac{d}{dt}\|v\|_{[l]}$. On the other hand, if $\|v\|_{[l]} \equiv 0$ on an open sub interval of $(s,t_0)$, then it is clear $\frac{d}{dt}\|v\|_{[l]} \equiv 0 \leq -c_0\psi\lt(\|x\|_{[l]})\rt)\|v\|_{[l]} + c_1\|x\|_{[l]} = c_1\|x\|_{[l]}$ on that interval. Thus we obtain
\[
\frac{d}{dt}\|v\|_{[l]}\leq -c_0\psi\lt(\|x\|_{[l]})\rt)\|v\|_{[l]} + c_1\|x\|_{[l]} \quad a.e. \mbox{ on } (s,t_0).
\]
Applying Gronwall's inequality to the above differential inequality together with the continuity of $\|v\|_{[l]}$ on the time interval $(s,t_0)$ yields
\begin{align}\label{rbis}
\|v\|_{[l]}(t)\leq \left(c_1\int_s^t\|x\|_{[l]}(\tau) e^{c_0\int_s^\tau\psi(\|x\|_{[l]}(\sigma))d\sigma}d\tau +\|v\|_{[l]}(s)\right)\exp\lt(-c_0\int_s^t\psi(\|x\|_{[l]}(\tau))\,d\tau\rt).
\end{align}

The estimates \eqref{rprime} and \eqref{rbis} enable us to finally approach the conclusion of the proof. Let us recall $\Psi$ is the primitive of $\psi$ given in \eqref{def_pri}.
Then, by (\ref{rprime}) in the interval $(s,t_0)$, we have
\begin{align*}
\left|\Psi(\|x\|_{[l]}(t))\right| &= \left|\int_s^t\frac{d}{dt}\Psi(\|x\|_{[l]}(\tau))\,d\tau
 +\Psi(\|x\|_{[l]}(s))\right|\\
&=\left|\int_s^t\psi(\|x\|_{[l]}(\tau))\left(\frac{d}{dt}\|x\|_{[l]}\right)(\tau)\,d\tau
+ \Psi(\|x\|_{[l]}(s))\right|\\
&\leq \int_s^t\psi(\|x\|_{[l]}(\tau))\|v\|_{[l]}(\tau)\,d\tau +|\Psi(\|x\|_{[l]}(s))|.
\end{align*}
We apply (\ref{rbis}) to obtain
\begin{align}
\left|\Psi(\|x\|_{[l]}(t))\right|
&\leq \!\!\int_s^t\!\!\psi(\|x\|_{[l]}(\tau))\left(c_1\int_s^\tau\|x\|_{[l]}(\sigma) e^{c_0\int_s^\sigma\psi(\|x\|_{[l]}(\rho))d\rho}d\sigma +\|v\|_{[l]}(s)\right) e^{-c_0\int_s^\tau\psi(\|x\|_{[l]}(\sigma))d\sigma}d\tau\nonumber\\
&\quad +|\Psi(\|x\|_{[l]}(s))| \cr
&=: A + |\Psi(\|x\|_{[l]}(s))|\label{a}.
\end{align}
In order to estimate $A$ in the above inequality let us simplify the notation by taking $a:=\psi(\|x\|_{[l]})$. Then we have
\begin{align*}
A &= c_1\int_s^ta(\tau)\int_s^\tau\|x\|_{[l]}(\sigma) e^{c_0\int_s^\sigma a(\rho)d\rho}d\sigma\ e^{-c_0\int_s^\tau a(\sigma)d\sigma}\,d\tau\\
&\quad + \|v\|_{[l]}(s)\int_s^ta(\tau)e^{-c_0\int_s^\tau a(\sigma)d\sigma}\,d\tau =: A_1 + A_2.
\end{align*}
By Proposition \ref{bounds} there exists a constant $c_3=c_3(T) = \max\{M,R\}$ such that
\[
c_1\|x\|_{[l]}\leq c_3\quad \mbox{and}\quad \|v\|_{[l]}\leq c_3 \quad\mbox{on}\quad [0,t_0).
\]
Therefore we have
\begin{align*}
A_1\leq c_3\int_s^t\left(\int_s^\tau e^{c_0\int_s^\sigma a(\rho)d\rho}d\sigma\right)\left(a(\tau) e^{-c_0\int_s^\tau a(\sigma)d\sigma}\right)\,d\tau
\end{align*}
and noting that
\begin{align*}
\int_s^t a(\tau)e^{-c_0\int_s^\tau a(\sigma)\,d\sigma}d\tau = -\frac{1}{c_0}e^{-c_0\int_s^t a(\tau)\,d\tau}
\end{align*}
by integration by parts we obtain
\begin{align*}
A_1 &\leq -\frac{c_3}{c_0}\int_s^te^{c_0\int_s^\tau a(\sigma)d\sigma}d\tau \ e^{-c_0\int_s^ta(\tau)d\tau} + \frac{c_3}{c_0}\int_s^t e^{c_0\int_s^\tau a(\sigma)d\sigma} e^{-c_0\int_s^\tau a(\sigma)\,d\sigma}\\
&\leq \frac{c_3}{c_0}T.
\end{align*}
Similar calculation reveals that
\begin{align*}
A_2\leq \frac{c_3}{c_0}\left(1-e^{-c_0\int_s^t a(\tau)d\tau}\right)\leq \frac{c_3}{c_0}
\end{align*}
and altogether
\begin{align*}
A\leq \frac{c_3}{c_0}(T+1).
\end{align*}
We apply the above estimate of $A$ in \eqref{a} to obtain
\begin{align*}
\left|\Psi(\|x\|_{[l]}(t))\right| \leq \frac{c_3}{c_0}(T+1) +|\Psi(\|x\|_{[l]}(s))|
\end{align*}
and by Proposition \ref{bounds} the right-hand side of the above inequality is bounded regardless of the choice of $s\in[0,t_0)$. Thus $|\Psi(\|x\|_{[l]})|$ is also bounded in $(s,t_0)$,
 which subsequently implies that $\|x\|_{[l]}$ is separated from $0$ in $(s,t_0)$. To be more precise, we have
\begin{align*}
\|x\|_{[l]}(t)\geq \left|\Psi^{-1}\left(\frac{c_3}{c_0}(T+1) + |\Psi(\|x\|_{[l]}(s))|\right)\right|>0.
\end{align*}
This contradicts \eqref{contr}, and consequently it proves that the assumption of existence of $t_0$, the first time of collision between any number of the particles was false.
\end{proof}

\begin{rem} The estimates of derivatives of $\|x\|_{[l]}$ and $\|v\|_{[l]}$, in principle, are similar to (and to some degree were inspired by) the SDDI derived in \cite{haliu}. In fact, we
divided the particles into two groups: one of those particles that collide with each other at $t_0$ and the second group of those that do not collide with the first group. Next, for the first group
we basically derived the SDDI similarly to \cite{haliu} with some additional terms that originated from the interaction between the groups.
\end{rem}

\begin{rem}As a direct application of \cite[Theorem 3.1]{ahn1}, we also have the flocking estimates: Consider the CS particle system \eqref{cspart}-\eqref{ini_cspart}. Suppose that the initial
 configurations $(x_0,v_0)$ satisfy
\[
x_{i0}\neq x_{j0} \quad \mbox{for } 1 \leq i \neq j \leq N \quad \mbox{and} \quad \|v_0 - v_c(0)\|_{\ell^\infty} < \frac12 \int_{2\|x_0 - x_c(0)\|_{\ell^\infty}}^\infty \psi(s)\,ds,
\]
where $x_c(t)$ and $v_c(t)$ denote the average quantities of the position and velocity, i.e.,
\[
x_c(t) = \frac1N\sum_{i=1}^N x_i(t) \quad \mbox{and} \quad v_c(t) = \frac1N\sum_{i=1}^N v_i(t),
\]
respectively. Then there exist positive constants $c_1 > c_0(t) > 0$ such that
\[
c_0(t)\leq |x_i(t) - x_j(t)| \leq c_1 \,, i\neq j\quad \mbox{and} \quad \|v(t) - v_c(0)\|_{\ell^\infty} \leq \|v_0 - v_c(0)\|_{\ell^\infty} e^{-\psi(2c_1) t},
\]
for $t \geq 0.$ Note that $c_0(t)$ might go to zero as $t\to\infty$ but $c_1$ is uniform in time. This implies that the communication rate is bounded below uniformly in time giving the exponential convergence to the mean velocity, see \cite{ahn1,car} for details. In particular, for the critical case, $\alpha = 1$, if we assume
\[
x_{i0}\neq x_{j0} \quad \mbox{for } 1 \leq i \neq j \leq N,
\]
then there is no collision for all time and we have exponential flocking estimate because $\int^\infty \psi(s)\,ds = \infty$.
\end{rem}

\section{Uniform estimate of minimal interparticle distance}\label{estim}

In this part, we consider system \eqref{dcspart} denoting the singular communication weights $\psi_\delta(s) = \psi(s-\delta)$, where $\delta \geq 0$ is a control parameter to
make particles stay away from each other. We subject system \eqref{dcspart} to the initial data
\bq\label{dini_cspart}
(x_i,v_i)(0) =: (x_{i0}, v_{i0}), \quad i =1,\cdots, N.
\eq
It is worth mentioning that the strategy used in Section \ref{noncoll} cannot be used for the system \eqref{dcspart}. Thus we propose a new argument based on an energy estimate.
For $\delta=0$, Theorem \ref{thm_d} implies the lack of collisions between particles and on top of that, it provides a uniform with respect to $N$ estimate of non-collisioness for system \eqref{cspart}.

For a fixed $\delta>0$ we introduce a function $\ml^\beta(t)$ which determines the $\delta$-distance between particles:
\[
\ml^\beta(t) := \frac{1}{N^2}\sum_{i,j=1}^N\left(|x_i(t) - x_j(t)|-\delta\right)^{-\beta} \quad \mbox{with} \quad \beta > 0.
\]
Note that there exists a $\beta > 0$ such that $\ml^\beta(t) < \infty$ for $t \in [0,T]$ if and only if the distances between particles are no less than $\delta$ for $t \in [0,T]$. Thus the
solution to the $\delta$-CS system \eqref{dcspart} is well-defined provided that we start from admissible initial data i.e., $|x_{i0}-x_{j0}|>\delta$ for $i\neq j$. Notice that the definition of $\ml^\beta$ and its property of  determining the $\delta$-distance between particles is valid for all $\beta>0$. However in the proof of Theorem \ref{thm_d} we take $\beta = \alpha -2$.

We introduce the maximal life-span $T(x_0)$ of the initial datum $x_0$:
\[
T(x_0):= \sup\lt\{ s \in \R_+: \exists \mbox{ solution $(x(t),v(t))$ for the system \eqref{dcspart} in a time-interval $[0,s)$}  \rt\}.
\]
\begin{theo}\label{thm_d}Suppose that $\alpha \in [2,\infty)$ and the initial data $x_0$ satisfy
\begin{equation}\label{assume_basicd}
|x_{i0} - x_{j0}|>\delta \quad \mbox{for any} \quad 1 \leq i \neq j \leq N.
\end{equation}
Then there exists a global smooth solution $(x(t),v(t))$ to the system \eqref{dcspart}, i.e., $T(x_0) = \infty$. Moreover for $t\geq 0$ we have
\begin{align}
&\lt|\frac{1}{N^2}\sum_{i,j = 1}^N\log \lt(|x_i(t) - x_j(t)|-\delta\rt)\rt| \cr
&\qquad \quad \leq \lt|\frac{1}{N^2}\sum_{i,j = 1}^N\log \lt(|x_i(0) - x_j(0)|-\delta\rt)\rt| + \frac{T_0}{2} + \frac{1}{2N}\sum_{i=1}^N |v_i(0)|^2 \quad \mbox{if} \quad \alpha = 2.\label{estim1d}\\
&\ml^{\alpha - 2}(t) \leq \ml^{\alpha - 2}(0)e^{Ct} + Ce^{Ct}\frac{1}{N}\sum_{i=1}^N |v_i(0)|^2 \hspace{4.3cm} \mbox{if} \quad \alpha>2,\label{estim2d}
\end{align}
where $C$ is a constant depending only on $\alpha$.
\end{theo}
\begin{proof}
Throughout the proof we will denote for simplicity $T_0:=T(x_0)$.
Clearly, if the distances between particles are bigger than $\delta$ for any finite time then the solution can be prolonged indefinitely and $T_0=\infty$. Moreover if \eqref{estim1d}
or \eqref{estim2d} holds then $\inf_{t\geq 0}\min_{i,j=1,...,N}|x_i(t)-x_j(t)|>\delta$ and thus also $T_0=\infty$. Therefore in order to finish the proof it suffices to show \eqref{estim1d} and \eqref{estim2d}.

First we provide the energy dissipation of the system \eqref{dcspart}. By $\eqref{dcspart}_2$ and a symmetry argument, we have
\[
\frac{d}{dt} \frac{1}{N}\sum_{i=1}^N |v_i(t)|^2 + \frac{1}{N^2}\sum_{i,j=1}^N\frac{|v_i(t) - v_j(t)|^2}{(|x_i(t) - x_j(t)|-\delta)^\alpha} = 0 \quad \mbox{for} \quad t \in [0,T_0),
\]
and this yields
\begin{equation}\label{est_energy}
\frac{1}{N}\sum_{i=1}^N |v_i(t)|^2 + \int_0^t \frac{1}{N^2}\sum_{i,j=1}^N\frac{|v_i(s) - v_j(s)|^2}{(|x_i(s) - x_j(s)|-\delta)^\alpha} ds =  \frac{1}{N}\sum_{i=1}^N|v_i(0)|^2 \quad \mbox{for} \quad t \in [0,T_0).
\end{equation}
We divide the proof into two cases: $\alpha = 2$ and $\alpha > 2$.  (i) $\alpha = 2$: We estimate
$$\begin{aligned}
\lt|\frac{d}{dt} \frac{1}{N^2}\sum_{i,j = 1}^N\log \lt(|x_i(t) - x_j(t)|-\delta\rt)\rt| &= \lt|\frac{1}{N^2}\sum_{i,j=1}^N \frac{(x_i(t) - x_j(t))\cdot (v_i(t) - v_j(t))}{|x_i(t) - x_j(t)|}\frac{1}{(|x_i(t)-x_j(t)|-\delta)}\rt|\cr
&\leq \frac{1}{N^2}\sum_{i,j=1}^N \frac{|v_i(t) - v_j(t)|}{|x_i(t) - x_j(t)|-\delta}\cr
&\leq \frac{1}{2} + \frac{1}{2N^2}\sum_{i,j=1}^N \frac{|v_i(t) - v_j(t)|^2}{(|x_i(t) - x_j(t)|-\delta)^2},
\end{aligned}$$
for $t \in [0,T_0)$. We use the energy estimate \eqref{est_energy} to have
$$\begin{aligned}
&\lt|\frac{1}{N^2}\sum_{i,j = 1}^N\log \lt(|x_i(t) - x_j(t)|-\delta\rt)\rt| \cr
&\quad \leq \lt|\frac{1}{N^2}\sum_{i,j = 1}^N\log \lt(|x_i(0) - x_j(0)|-\delta\rt)\rt| + \frac{t}{2} + \frac{1}{2}\int_0^t \frac{1}{N^2}\sum_{i,j=1}^N\frac{|v_i(s) - v_j(s)|^2}{(|x_i(s) - x_j(s)|-\delta)^2} ds\cr
&\quad \leq \lt|\frac{1}{N^2}\sum_{i,j = 1}^N\log \lt(|x_i(0) - x_j(0)|-\delta\rt)\rt| + \frac{T_0}{2} + \frac{1}{2N}\sum_{i=1}^N |v_i(0)|^2 \quad \mbox{for} \quad t \in [0,T_0).
\end{aligned}$$
This yields \eqref{estim1d} in case of $\alpha=2$.

(ii) $\alpha > 2$: It is a straightforward to get
$$\begin{aligned}
\frac{d \ml^\beta(t)}{dt} &= -\beta\frac{1}{N^2}\sum_{i,j=1}^N (|x_i(t) - x_j(t)|-\delta)^{-\beta - 1}\frac{(x_i(t) - x_j(t))\cdot (v_i(t) - v_j(t))}{|x_i(t)-x_j(t)|}\cr
&\leq C\frac{1}{N^2}\sum_{i,j = 1}^N (|x_i(t) - x_j(t)|-\delta)^{-\beta - 1}|v_i(t) - v_j(t)|\cr
&\leq C\frac{1}{N^2}\sum_{i,j = 1}^N\frac{1}{(|x_i(t) - x_j(t)|-\delta)^\beta} + C\frac{1}{N^2}\sum_{i,j = 1}^N\frac{|v_i(t) - v_j(t)|^2}{(|x_i(t) - x_j(t)|-\delta)^{\beta +2}}\cr
&= C\ml^\beta(t) + C\frac{1}{N^2}\sum_{i,j = 1}^N\frac{|v_i(t) - v_j(t)|^2}{(|x_i(t) - x_j(t)|-\delta)^{\beta +2}} \quad \mbox{for} \quad t \in [0,T_0),
\end{aligned}$$
where we used the Young's inequality. Then, by using Gronwall's inequality, we obtain
\[
\ml^\beta(t) \leq \ml^\beta(0)e^{Ct} + Ce^{Ct}\int_0^{t} \frac{1}{N^2}\sum_{i,j = 1}^N\frac{|v_i(s) - v_j(s)|^2}{(|x_i(s) - x_j(s)|-\delta)^{\beta +2}} ds \quad \mbox{for} \quad t \in [0,T_0).
\]
We choose $\alpha = \beta + 2$ together with \eqref{est_energy} to deduce
\[
\ml^{\alpha - 2}(t) \leq \ml^{\alpha - 2}(0)e^{Ct} + Ce^{Ct}\frac{1}{N}\sum_{i=1}^N |v_i(0)|^2 \quad \mbox{for} \quad t \in [0,T_0)
\]
and the proof of \eqref{estim2d} is finished.
Hence if $\alpha \geq 2$, then there is no collision between particles  for all time and $T_0=\infty$.
\end{proof}
\begin{rem}Similarly as before, we also have the flocking estimates for the system \eqref{dcspart}. Suppose that the initial configurations $(x_0,v_0)$ satisfy
\[
\min_{1 \leq i \neq j \leq N}|x_{i0} - x_{j0}| > \delta,
 \quad \|x_0 - x_c(0)\|_{\ell^\infty} > \delta, \quad \mbox{and} \quad \|v_0 - v_c(0)\|_{\ell^\infty} < \frac12 \int_{2\|x_0 - x_c(0)\|_{\ell^\infty}}^\infty \psi(s - \delta)\,ds,
\]
Then if $\alpha \geq 2$ there exist positive constants $c_1 > c_0 > \delta > 0$ such that
\[
\|x(t) - x_c(0)\|_{\ell^\infty} \in [c_0,c_1] \quad \mbox{and} \quad \|v(t) - v_c(0)\|_{\ell^\infty} \leq \|v_0 - v_c(0)\|_{\ell^\infty} e^{-\psi(2c_1) t},
\]
for $t \geq 0.$
\end{rem}

\section*{Acknowledgments}
\small{JAC was partially supported by the Royal Society via a Wolfson Research Merit Award. YPC was supported by the ERC-Starting grant HDSPCONTR ``High-Dimensional Sparse Optimal Control''.  JAC and YPC were partially supported by EPSRC grant EP/K008404/1. YPC is also supported by the Alexander Humboldt Foundation through the Humboldt Research Fellowship for Postdoctoral Researchers. JP was supported by the Polish NCN grant PRELUDIUM 2013/09/N/ST1/04113}

\bibliographystyle{abbrv}