So, let's dig into this.
The question is ambigously asked.
1. Real-Valued interpretation
If the question is meant to read
\$x, y\sim \mathcal N(0,\sigma^2)\$: what is the distribution of \$A= \lvert\angle x-\angle y\rvert\$?
then the answer would be: since \$\angle x = 0 \, \forall x \in\mathbb R\$ (and the same for \$y\$), \$ a = \lvert\angle x-\angle y\rvert = \lvert0-0\rvert = 0\implies f_A(a) = \delta(a)\$.
The only choice that would be correct here would be "a)"; because the Dirac delta is the marginal case of the zero-mean normal distribution with zero variance.
2. Complex-value interpretation
I'm swinging towards this interpretation, because it makes a little more sense.
Likely, however, the question is meant to read
\$x, y\sim \mathcal{CN}(0,\sigma^2)\$: what is the distribution of \$A= \lvert\angle x-\angle y\rvert\$?
(which contradicts the literal text; when you say "normally distributed", you mean the real \$\mathcal N\$, not the circularly complex \$\mathcal {CN}\$)
You'd then realize that the phase of a circularly complex normal variable is uniform between \$0\$ and \$2\pi\$ or between \$-\pi\$ and \$\pi\$ (or, however you define your range of valid angles). It doesn't really matter; all it does in the end translate the PDF.
I'll pick \$\angle x,\angle y\sim \mathcal U(-\pi,\pi)\$ because it allows the PDF to be symmetric, which will be useful when we calculate the following:
\begin{align}
\tilde A &= \angle X - \angle Y& \tilde Y = -Y\\
&= \angle X + \angle\tilde Y\\
\implies\\
f_A(a) &= f_{\angle X} * f_{\angle\tilde Y}&*\text{ being the convolution}\\
&= f_{\angle X} * f_{\angle Y} & \text{due to symmetry of Y's PDF}\\
&= \text{rect}(\frac{u}{2\pi})*\text{rect}(\frac{u}{2\pi}) & \text{convolution of two $2\pi$ wide zero-centered rectangles}\\
\end{align}
That's a triangular distribution. If you look at \$f_{\lvert A \rvert}\$, it becomes a linear downslope. None of the options allow a triangular distribution.
3. Interpration as coordinates of point in plane
If the question is meant to read
\$x, y\sim \mathcal N(0,\sigma^2)\$: what is the distribution of \$A=\angle (x,y)\$?
You'd realize that \$(x,y)\sim\mathcal {CN}(0, \frac{\sigma^2}2)\$, and all these zero-mean circularly complex normally distributed variables have uniform phase over a range of \$2\pi$. That'd imply b) is right.