How to form the spin triplet - singlet states from two electrons with spin not an eigenstate of $ S_z $ (spins not along z-axis)) For a two electron system, we know that the total $ J^2 $ states (Triplet - Singlet) are related with the $\uparrow \downarrow $ , $\downarrow \uparrow$ , $\uparrow \uparrow$ , $\downarrow \downarrow $ states as:
$$
\[
\left(\begin{array}{c}
\left|11\right\rangle \\
\left|10\right\rangle \\
\left|1-1\right\rangle \\
\left|00\right\rangle 
\end{array}\right)=\left(\begin{array}{cccc}
1 & 0 & 0 & 0\\
0 & \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} & 0\\
0 & 0 & 0 & 1\\
0 & \frac{1}{\sqrt{2}} & -\frac{1}{\sqrt{2}} & 0
\end{array}\right)\left(\begin{array}{c}
\uparrow\uparrow\\
\uparrow\downarrow\\
\downarrow\uparrow\\
\downarrow\downarrow
\end{array}\right)
\] 
$$
But what if we want the singlet-triplet not in terms of 
$\uparrow = \left(\begin{array}{c}1\\0\end{array}\right)$ , 
$\downarrow =\left(\begin{array}{c}0\\1\end{array}\right)$ , but in terms of spin states that are at an angle $\theta$ from the $z$ axis?
For example when $\hat\eta $ is at $60^o$ from the $z$ axis the spin up and down states are 
\begin{eqnarray*}
\nearrow & = & \left(\begin{array}{c}\frac{\sqrt{3}}{2}\\
\frac{1}{2}\end{array}\right)=\frac{\sqrt{3}}{2}\uparrow+\frac{1}{2}\downarrow\\
\swarrow & = & \left(\begin{array}{c}-\frac{1}{2}\\\frac{\sqrt{3}}{2}
\end{array}\right)=-\frac{1}{2}\uparrow+\frac{\sqrt{3}}{2}\downarrow
\end{eqnarray*}
By solving for $\uparrow$ and $\downarrow$ in terms of $\nearrow$ and $\swarrow$ and rewriting the $\uparrow \downarrow $ , $\downarrow \uparrow$ , $\uparrow \uparrow$ , $\downarrow \downarrow $ states in terms of the 
$\nearrow \nearrow$ ,$\nearrow \swarrow$ , $\swarrow \nearrow$ , $\swarrow \swarrow$
we get (after some algebra): $$\left(\begin{array}{c}
\left|11\right\rangle \\
\left|10\right\rangle \\
\left|1-1\right\rangle \\
\left|00\right\rangle 
\end{array}\right)=\left(\begin{array}{cccc}
\frac{3}{4} & -\frac{\sqrt{3}}{4} & -\frac{\sqrt{3}}{4} & \frac{1}{4}\\
\sqrt{\frac{3}{8}} & \sqrt{\frac{1}{8}} & \sqrt{\frac{1}{8}} & -\sqrt{\frac{3}{8}}\\
\frac{1}{4} & \frac{\sqrt{3}}{4} & \frac{\sqrt{3}}{4} & \frac{3}{4}\\
0 & \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} & 0
\end{array}\right)\left(\begin{array}{c}
\nearrow\nearrow\\
\nearrow\swarrow\\
\swarrow\nearrow\\
\swarrow\swarrow
\end{array}\right)$$
What is the general method of writing the Triplet - Singlet in terms of two electron states with spin at an arbitrary axis $\hat\eta (\theta)$  ?
Thank you!
 A: You are looking for the Wigner d-functions. They relate angular momentum eigenstates through rotation. As you can see in the link the definition is 
$$d^{(j)}_{m,m'}(\theta) = \langle jm|e^{-i\theta J_y}|j m' \rangle$$
where $e^{-i\theta J_y}$ is a unitary rotation operator.
We'll have two sets of states: $|jm;0\rangle$ for the original basis and $|jm;\theta\rangle$ for eigenstates of momentum along the axis $\hat{\eta}(\theta)$. Let us assume that we've set up our basis so $\hat{\eta}$ lies in the $x-z$ plane. That way the two bases are related to each other by a rotation along the $y$ axis. The two sets of states are then related by a unitary rotation operator: $|jm;\theta\rangle = e^{-i\theta J_y}|jm;0\rangle $
You can look up the various $d$-functions here: http://pdg.lbl.gov/2002/clebrpp.pdf
Using $|jm;0\rangle$ as a complete set of states we can write
$$|jm;\theta\rangle = \sum_{m'} |jm';0\rangle\langle jm';0|jm;\theta\rangle$$
$$= \sum_{m'} |jm';0\rangle\langle jm';0|e^{-i\theta J_y}|jm;0\rangle $$
$$|jm;\theta\rangle= \sum_{m'}  d^{(j)}_{m',m}(\theta) |jm';0\rangle$$
So you want to perform two steps. First, rotate your total angular momentum state into the basis you'd like to compare to with the $d$-functions. Then second, use the Clebsch-Gordan coefficients to decompose the total-angular momentum states into tensor-product states. The final result is 
$$|jm;\theta\rangle= \sum_{j_1,j_2,m_1,m_2} \left(\sum_{m'}  d^{(j)}_{m',m}(\theta) \langle j_1 m_1;j_2 m_2|j m'\rangle \right) |j_1 m_1; j_2 m_2;\theta\rangle$$
I'll leave you to compare this expression to the matrix you wrote out in your question, but you'll find that plugging in $\theta = 60^\circ$ gives you exactly what you already have.
