Step 1: Expand the vector triple product.
Using the identity $\vec a\times(\vec b\times\vec c)=(\vec a\cdot\vec c)\vec b-(\vec a\cdot\vec b)\vec c$, the left side becomes \[ (\vec a\cdot\vec c)\vec b-(\vec a\cdot\vec b)\vec c+(\vec a\cdot\vec b)\vec b. \]
Step 2: Group by $\vec b$ and $\vec c$.
This is $\big[(\vec a\cdot\vec c)+(\vec a\cdot\vec b)\big]\vec b-(\vec a\cdot\vec b)\vec c$, set equal to $(4-2\beta-\sin\alpha)\vec b+(\beta^2-1)\vec c$.
Step 3: Use the condition on $\vec a$.
From $(\vec c\cdot\vec c)\vec a=\vec c$ we get $\vec a=\dfrac{\vec c}{|\vec c|^2}$, hence $\vec a\cdot\vec c=1$ and $\vec a\cdot\vec b=\dfrac{\vec b\cdot\vec c}{|\vec c|^2}$, call this $\mu$.
Step 4: Compare coefficients.
Since $\vec b,\vec c$ are non-collinear we can match each: for $\vec c$, $-\mu=\beta^2-1$; for $\vec b$, $1+\mu=4-2\beta-\sin\alpha$.
Step 5: Add the two equations.
The $\mu$ terms cancel, giving $1=(4-2\beta-\sin\alpha)+(\beta^2-1)$, so $\sin\alpha=\beta^2-2\beta+2=(\beta-1)^2+1$.
Step 6: Solve using the range of sine.
Since $\sin\alpha\le1$ but $(\beta-1)^2+1\ge1$, equality forces $(\beta-1)^2=0$, so $\beta=1$ and $\sin\alpha=1$, giving $\alpha=2n\pi+\tfrac{\pi}{2}$. This is option 4.
\[ \boxed{\beta=1,\ \alpha=2n\pi+\tfrac{\pi}{2},\ n\in z} \]