Запишем второй закон Ньютона в проекции на линию, проходящую через электроны. Поскольку расстояние не меняется, ускорение равно $\frac{v^2}{r}$. Получаем в итоге уравнение:
$$
\frac{2mv^2}{d}=vBe-\frac{e^2}{4\pi\varepsilon_0 d^2}
$$
Решаем квадратное уравнение относительно $v$. Дискриминант этого уравнения равен $(Be)^2-\frac{2me^2}{\pi\varepsilon_0d^3}$.
Из дискриминанта видно, что решения для скорости существуют только для $d\geqslant2\sqrt[3]{\frac{m}{4\pi\varepsilon_0 B^2}}$. Это и будет минимальное возможное расстояние.
Так как дискриминант равен нулю: $$v_m = -\frac{b}{2a}= \frac{Be}{2m}d_{\text{min}}=\frac{e}{2}\sqrt[3]{\frac{B}{4\pi\varepsilon_0m^2}}$$
Центр масс движется под действием двух сил Лоренца: $$m\vec{a} = \left[\vec{v}_1\times \vec B\right]q+\left[\vec{v}_2\times \vec B\right]q = \left[\left(\vec{v}_1+\vec{v}_2\right)\times \vec B\right]q = 2\left[\vec{v}_{CM}\times \vec B\right]q$$
Так как других внешних сил нет,скорость центра масс поворачивается, но не меняется по модулю. Ускорение центра масс в таком случае равно $\vec\omega\times\vec v$
$$2m\left[\vec\omega_{CM}\times\vec v\right] \equiv 2\left[\vec{v}_{CM}\times \vec B\right]q$$
$$\omega_{CM} = \frac{Be}{m}$$
При движении, описанном в условии, кинетическая энергия сохраняется. Таким образом угловая скорость относительно центра масс не меняется. $$\omega = \frac{v_0}{d_1}$$
Скорость центра масс и относительная скорость равны, поэтому скорости электронов сравняются, когда угол между скоростью центра масс и относительной скоростью станет равным $\frac{\pi}{2}$.
$$t = \frac{\pi}{2}\frac{1}{\left|\frac{v_0}{d_1} - \frac{eB}{m}\right|}$$
$K=\frac{m}{2}(\omega^2r^2+v_r^2)$
Момент импульса меняется за счет составляющей силы Лоренца, отвечающей радиальной составляющей скорости.
$$\frac{dL}{dt} = v_rBer$$
Интегрируя выражение для $\frac{dL}{dt}$ получаем интеграл движения
$$J = L-\frac{eBr^2}{2}$$
Выражаем угловую скорость через $J$, затем подставляем в выражение для кинетической энергии. $\omega = \frac{1}{m}\left(\frac{J}{r^2}+\frac{eB}{2}\right)$.
К кинетической энергии надо добавить энергию кулоновского взаимодействия, приходящуюся на один электрон. $$E = \frac{mv_r^2}{2}+\frac{r^2}{2m}\left(\frac{J}{r^2}+\frac{eB}{2}\right)^2+\frac{e^2}{16\pi\varepsilon_0r}$$
Пределы $U(r)$ в точках $0$ и $+\infty$ равны $+\infty$, при этом функция непрерывна на отрезке $(0, +\infty)$
Вторая производная функции $U(r)$ имеет следующий вид:
$$\frac{Ar^4+Br+C}{r^4}; ~~A,~B,~C, ~r > 0$$
На интересующем нас отрезке она всегда больше нуля, оба предела бесконечные, следовательно есть минимум, причем единственный.
Так как мы будем оставаться в окрестности минимума, электроны будут совершать гармонические колебания вдоль радиуса.