I'll use 3-D coordinates with the $z$-axis passing through the poles, so we start at $(\cos \theta, 0, \sin \theta)$. The optimal path $\gamma$ must be a geodesic, so $\gamma$ is an arc of a great circle through $(\cos \theta, 0, \sin \theta)$. Let $\nu$ be a vector in $\mathbb{R}^3$ othogonal to $\gamma$. Since $(\cos \theta, 0, \sin \theta) \in \gamma$, we have $(\cos \theta, 0, \sin \theta) \cdot \nu = 0$.

Let the end point of this geodesic have longitude $\psi$. Then $\gamma$ must be the shortest path to the line of longitude $(-\sin \psi, \cos \psi, 0)^{\perp}$. So $\gamma$ and $(-\sin \psi, \cos \psi, 0)^{\perp}$ must meet orthogonally, which means that $(-\sin \psi, \cos \psi, 0) \cdot \nu = 0$.

So $\nu$ is a multiple of $$(\cos \theta, 0, \sin \theta) \times (-\sin \psi, \cos \psi, 0) = (-\cos (\psi) \sin (\theta),-\sin (\psi) \sin (\theta),\cos (\psi) \cos (\theta )).$$

The point of intersection of the geodesic $\nu^{\perp}$ with the line of longitude $(-\sin \psi, \cos \psi, 0)^{\perp}$ is a scalar multiple of $$(-\sin \psi, \cos \psi, 0) \times \nu= (\cos ^2(\psi) \cos (\theta),\sin (\psi) \cos (\psi) \cos (\theta),\sin (\theta)).$$

One can compute the angle between $(\cos \theta, 0, \sin \theta)$ and $(\cos ^2(\psi) \cos (\theta),\sin (\psi) \cos (\psi) \cos (\theta),\sin (\theta))$ in order to find out how long the path is. I get the answer that the length of the path is $$\sin^{-1}(\cos \theta \sin \phi).$$ That's simple enough that there is probably a better derivation; I just slogged through trig identities to get it.

I'll use 3-D coordinates with the $z$-axis passing through the poles, so we start at $\alpha:=(\cos \theta, 0, \sin \theta)$. The optimal path $\gamma$ must be a geodesic, so $\gamma$ is an arc of a great circle through $\alpha$. Let $\nu$ be a vector in $\mathbb{R}^3$ orthogonal to $\gamma$. Since $\alpha \in \gamma$, we have $\alpha \cdot \nu = 0$.

Let the end point of this geodesic have longitude $\psi$. Then $\gamma$ must be the shortest path to the line of longitude $(-\sin \psi, \cos \psi, 0)^{\perp}$. Set $\lambda := (-\sin \psi, \cos \psi, 0)$. So $\gamma=\nu^{\perp}$ and $\lambda^{\perp}$ must meet orthogonally, which means that $\lambda \cdot \nu = 0$.

So $\nu$ is a multiple of $$\alpha \times \lambda = (-\cos (\psi) \sin (\theta),-\sin (\psi) \sin (\theta),\cos (\psi) \cos (\theta )).$$

The point of intersection of the geodesic $\nu^{\perp}$ with the line of longitude $\lambda^{\perp}$ is a scalar multiple of $$\lambda \times \nu=-(\lambda \times (\lambda \times \alpha))=(\cos ^2(\psi) \cos (\theta),\sin (\psi) \cos (\psi) \cos (\theta),\sin (\theta)).$$ Now, taking the cross product with $\lambda$ twice is just (up to scalar) orthogonal projection onto $\lambda^{\perp}$, so we can write this more simply as $\alpha - \langle \lambda, \alpha \rangle \lambda$. And the angle between $\alpha$ and the orthogonal projection of $\alpha$ onto $\lambda^{\perp}$ is $\sin^{-1} (\alpha \cdot \lambda)$ (using that $\alpha$ and $\lambda$ are both unit vectors), which is $$\sin^{-1} {\big(} (\sin \phi) (\cos \theta) {\big)}.$$

I'll use 3-D coordinates with the $z$-axis passing through the poles, so we start at $(\cos \theta, 0, \sin \theta)$. The optimal path $\gamma$ must be a geodesic, so $\gamma$ is an arc of a great circle through $(\cos \theta, 0, \sin \theta)$. Let $\nu$ be a vector in $\mathbb{R}^3$ othogonal to $\gamma$. Since $(\cos \theta, 0, \sin \theta) \in \gamma$, we have $(\cos \theta, 0, \sin \theta) \cdot \nu = 0$.

Let the end point of this geodesic have longitude $\psi$. Then $\gamma$ must be the shortest path to the line of longitude $(-\sin \psi, \cos \psi, 0)^{\perp}$. So $\gamma$ and $(-\sin \psi, \cos \psi, 0)^{\perp}$ must meet orthogonally, which means that $(-\sin \psi, \cos \psi, 0) \cdot \nu = 0$.

So $\nu$ is a multiple of $$(\cos \theta, 0, \sin \theta) \times (-\sin \psi, \cos \psi, 0) = (-\cos (\psi) \sin (\theta),-\sin (\psi) \sin (\theta),\cos (\psi) \cos (\theta )).$$

The point of intersection of the geodesic $\nu^{\perp}$ with the line of longitude $(-\sin \psi, \cos \psi, 0)^{\perp}$ is a scalar multiple of $$(-\sin \psi, \cos \psi, 0) \times \nu= (\cos ^2(\psi) \cos (\theta),\sin (\psi) \cos (\psi) \cos (\theta),\sin (\theta)).$$

One can compute the angle between $(\cos \theta, 0, \sin \theta)$ and $(\cos ^2(\psi) \cos (\theta),\sin (\psi) \cos (\psi) \cos (\theta),\sin (\theta))$, in order to find out how long the path is, but it doesn't seem to simplify in any nice way.

I'll use 3-D coordinates with the $z$-axis passing through the poles, so we start at $(\cos \theta, 0, \sin \theta)$. The optimal path $\gamma$ must be a geodesic, so $\gamma$ is an arc of a great circle through $(\cos \theta, 0, \sin \theta)$. Let $\nu$ be a vector in $\mathbb{R}^3$ othogonal to $\gamma$. Since $(\cos \theta, 0, \sin \theta) \in \gamma$, we have $(\cos \theta, 0, \sin \theta) \cdot \nu = 0$.

Let the end point of this geodesic have longitude $\psi$. Then $\gamma$ must be the shortest path to the line of longitude $(-\sin \psi, \cos \psi, 0)^{\perp}$. So $\gamma$ and $(-\sin \psi, \cos \psi, 0)^{\perp}$ must meet orthogonally, which means that $(-\sin \psi, \cos \psi, 0) \cdot \nu = 0$.

So $\nu$ is a multiple of $$(\cos \theta, 0, \sin \theta) \times (-\sin \psi, \cos \psi, 0) = (-\cos (\psi) \sin (\theta),-\sin (\psi) \sin (\theta),\cos (\psi) \cos (\theta )).$$

The point of intersection of the geodesic $\nu^{\perp}$ with the line of longitude $(-\sin \psi, \cos \psi, 0)^{\perp}$ is a scalar multiple of $$(-\sin \psi, \cos \psi, 0) \times \nu= (\cos ^2(\psi) \cos (\theta),\sin (\psi) \cos (\psi) \cos (\theta),\sin (\theta)).$$

One can compute the angle between $(\cos \theta, 0, \sin \theta)$ and $(\cos ^2(\psi) \cos (\theta),\sin (\psi) \cos (\psi) \cos (\theta),\sin (\theta))$ in order to find out how long the path is. I get the answer that the length of the path is $$\sin^{-1}(\cos \theta \sin \phi).$$ That's simple enough that there is probably a better derivation; I just slogged through trig identities to get it.