Let $z_1,\ldots,z_n$ be $n\ge 1$ distinct points of $\mathbb R^2$. Define the potential function $U: \mathbb R^2 \to\mathbb R$ by

$$U(x):=\sum_{1\le i\le n} \log(|x-z_i|),$$

where $|\cdot|$ denotes the Euclidean norm. Denote by $F$ be the negative gradient of $U$, i.e. 

$$F(x):=-\nabla U(x)=-\sum_{1\le i\le n} \frac{x-z_i}{|x-z_i|^2}.$$

For any $x\in\mathbb R^2$ that is different from $z_1,\ldots,z_n$, define the integral curve $Y\equiv Y_x$ by $Y(0)=x$ and

$$\frac{dY(t)}{dt}=F(Y(t)),\quad \forall t\ge 0.$$

By standard results about flows on vector fields, the curve $Y$ can be defined over some maximal domain $[0,\tau\equiv \tau_x)$ with $0<\tau\le \infty$. Can we prove 

$$\lim_{t\uparrow\tau} Y(t) \mbox{ exists and belongs to } \{z_1,\ldots,z_n\}?$$

If so, the space $\mathbb R^2$ is divided into a partition $V_1,\ldots, V_n$ defined by

$$V_i:=\big\{x\in \mathbb R^2: \lim_{t\uparrow\tau_x} Y_x(t)=z_i\big\}.$$