I have $n$ random points $u_1, u_2, \ldots, u_n$ on the unit circle. Points are uniformly distributed.
I want to think about the points as associated angles $\theta_1, \theta_2, \ldots, \theta_n$ that are uniformly picked in $(-\pi, \pi]$.
Suppose, we rotate the circle, thus $\theta_i$ is a function of the angle $\theta$ of the new axis in the original frame.
For example, if $\theta_1 = \dfrac{\pi}{2}$ and we rotate the circle $\dfrac{\pi}{6}$ in the positive direction, then $\theta_1\left(\dfrac{\pi}{6}\right) = \dfrac{2}{3}\pi$.
Angles are adjusted to always stay in $(-\pi, \pi]$ range. See the animation below
Denote $S_{\theta} \triangleq \theta_1(\theta) + \theta_2(\theta) + \ldots + \theta_n(\theta)$.
Does there exist a rotation $\theta$ such that $S_{\theta} = 0$ ?
My gut feeling, is that yes. But I do not see a direction to prove the claim.
What I have done, is that rotating a circle clock-wise creates $S_{\theta}$ monotonically decreasing, except, when a point "crosses a negative X axis", at that point $2\pi$ are added.
But I don't see why $S_{\theta}$ changes the sign. It could stay negative or positive always.

Yes, such a value exists.
First, let's understand how any of the $\theta_i(\theta)$ functions looks like. Let's choose the direction of rotation such that $\theta_i(\theta)$ increases, unless it encounters the discontinuity jumping between $-\pi$ and $\pi$. In that case, $\theta_i(\theta)$ is just a linear function, increasing with slope $1$. If it reaches the value $\pi$, it will jump to $-\pi$ and continue from there linearly with slope 1 until again the value reaches $\pi$, which happens at an increment of $\theta$ of exactly $2\pi$.
So it's easy to see understand that $\theta_i(\theta)$ is $2\pi$-periodic (rotating the circle by $\theta+2\pi$ leaves the new axis in exactly the same position as rotating it by $\theta$). For a starting value of e.g. $\theta_i(0)=\frac{\pi}3$, the function looks like this:
The key insight now is that, no matter the starting value of $\theta_i(0)$, the Rieman integral of $\theta_i(\theta)$ over $[0,2\pi]$ is always zero. One can do this by hand calculation, but an easier way to see this is to realize that because the function is $2\pi$-periodic, the integral doesn't change when one moves the integration interval around, as long as the length stays at $2\pi$.
In this case, we can choose the interval to be $[-\theta_i(0),2\pi-\theta_i(0)]$, which means the function will always look like this in this interval, starting from $0$ up to $\pi$, the jumping down to $-\pi$ and climb back up to $0$.
The integral over the positive part of the function cancels the integral over the negative part exactly.
We are almost done now!
Since $S(\theta)$ (I'll use that instead of $S_\theta$ to stress the function character) equals $\sum_{i=1}^n \theta_i(\theta)$, we get
$$ \int_0^{2\pi}S(\theta) = \int_0^{2\pi}\sum_{i=1}^n \theta_i(\theta) = \sum_{i=1}^n\int_0^{2\pi}\theta_i(\theta) = 0.$$
In order for $S(\theta)$'s integral to be $0$, there must be positive and negative function values (if all were $0$, we would be done already). So there exist $x_1, x_2$ with $S(x_1) < 0, S(x_2) > 0$. Since $S(\theta)$ is of course also $2\pi$-periodic, we can assume $x_1 < x_2$.
$S(\theta)$ is the sum of piecewise increasing functions, so $S(\theta)$ is as well. $S(x_1)$ is negative, from here $S(\theta)$ increases (linearly, so continuosly) until it reaches a point of discontinuity, where it jumps down by $2\pi$. So if it reached function value $0$ during the linear increase, we found what we were looking for. If not, it was still negative when the jump happened, so it is still negative after the jump down, when the function will again start to rise.
But we know that the function will reach a positve value at point $x_2$, so at some point in a linerly rising part of $S(\theta)$ it will have to obtain the value $0$, which concludes the proof.