The Villarceau circles can be parametrized by
$$\vec{x}=\left(\begin{array}{ccc}
\cos\phi & \sin\phi & \\
-\sin\phi & \cos\phi & \\
&&1
\end{array}
\right)
\left(
\begin{array}{c}
\rho\sin\theta\\
r+R\cos\theta\\
r\sin\theta
\end{array}
\right)\;,
$$
where $\theta$ is the angle along the circle, $\phi$ is the orientation of the circle around the $z$-axis and $\rho=\sqrt{R^2-r^2}$. At $\phi=0$, the derivatives with respect to $\theta$ and $\phi$ are
$$\frac{\partial\vec{x}}{\partial\theta}=
\left(
\begin{array}{c}
\rho\cos\theta\\
-R\sin\theta\\
r\cos\theta
\end{array}
\right)\;,
\;\;\;
\frac{\partial\vec{x}}{\partial\phi}=
\left(
\begin{array}{c}
r+R\cos\theta\\
-\rho\sin\theta\\
0
\end{array}
\right)\;.
$$
We want to know where changing $\phi$ will be least effective in getting us away from the circle. This is determined on the one hand by the magnitude of the derivative with respect to $\phi$ and on the other hand by the angle between the two derivatives; the rate at which we move away from the circle as we change $\phi$ is the product of that magnitude and the sine of that angle:
$$
\begin{eqnarray}
\frac{\mathrm ds}{\mathrm d\phi}
&=&
\left|\frac{\partial\vec{x}}{\partial\phi}\right|
\sqrt{1-\left(\frac{\frac{\partial\vec{x}}{\partial\phi}\cdot\frac{\partial\vec{x}}{\partial\theta}}{\left|\frac{\partial\vec{x}}{\partial\phi}\right|\left|\frac{\partial\vec{x}}{\partial\theta}\right|}\right)^2}
\\
&=&
\sqrt{\left|\frac{\partial\vec{x}}{\partial\phi}\right|^2-\left(\frac{\frac{\partial\vec{x}}{\partial\phi}\cdot\frac{\partial\vec{x}}{\partial\theta}}{\left|\frac{\partial\vec{x}}{\partial\theta}\right|}\right)^2}
\\
&=&
\sqrt{(r+R\cos\theta)^2+\rho^2\sin^2\theta-\left(\rho(r\cos\theta+R)/R\right)^2}\;.
\end{eqnarray}
$$
We need to find the minimum of this expression with respect to $\theta$. Setting the derivative of the radicand with respect to $\cos\theta$ to zero yields
$$R(r+R\cos\theta)-\rho^2\cos\theta-r\rho^2(r\cos\theta+R)/R^2=0\;,$$
which simplifies to
$$\cos\theta=-\frac{R}{r}\;.$$
Since $R> r$, this has no solution, so the extrema occur on the boundary at $\cos\theta=\pm1$. At these points, $|\partial\vec{x}/\partial\theta|$ is parallel to the projection of $|\partial\vec{x}/\partial\phi|$ into the cutting plane, so the result is the same as would have been obtained by using the angle between $|\partial\vec{x}/\partial\phi|$ and the plane instead:
$$
\frac{\mathrm ds}{\mathrm d\phi}=(R\pm r)\frac{r}{R}\;\;\text{for}\;\;\cos\theta=\pm1\;.
$$
The factor $R\pm r$ is $|\partial\vec{x}/\partial\phi|$, and the factor $r/R$ is the sine of the angle between $\partial\vec{x}/\partial\phi$ and the plane. The lesser of these is the inner one at $\cos\theta=-1$, and the angle by which you need to rotate to make two tori with minor radii $r'$ fit is thus (to first order for $r'\ll r$)
$$\Delta\phi=\frac{2r'}{\mathrm ds/\mathrm d\phi}=\frac{2r'R}{(R-r)r}\;.$$