Rotation of second quantized operator in Fock space I'm quite familiar with rotation in quantum/classical mechanics. I know rotation for an operator $O$ or state $|\psi \rangle$ acts like:
$$O \rightarrow R O R^{-1} \\ |\psi \rangle \rightarrow R |\psi \rangle $$
However, I don't understand how to apply this to the second quantized operator for example $c_{i \sigma} c_{j \sigma}^\dagger$ where $\sigma$ is spin index. I mean naively I can rotate this like an operator and so acts like: 
$$c_{i \sigma} c_{j \sigma}^\dagger \rightarrow R c_{i \sigma} R^{-1} R c_{j \sigma}^\dagger R^{-1}$$ where $R \in SU(2)$, but then $c$ is in a different space from $R$. Maybe I'm supposed to do this? 
$$R c_{i \sigma} R^{-1} = c_{i R\sigma R^{-1}}$$
And if the above is true, then I have another confusion: in quantum mechanics, say two particle state, the only rotation invariant state is the singlet $\frac{1}{\sqrt{2}} |\uparrow  \downarrow \rangle - | \downarrow \uparrow  \rangle$ with $S(S+1) = 0$, and so analogously I think this operator is also rotational invariant
$$c_{i \uparrow} c_{j \downarrow} -c_{i \downarrow} c_{j \uparrow}  $$
but this is not the only rotational invariant operator? Because it seems like $c_{i \sigma} c_{j \sigma}^\dagger$ is also invariant.. Is my thinking right? 
In summary, I may be confused about how rotation acts in Fock space (for second quantized operator) versus Hilbert space (for spin state). 
 A: Let's focus, like your question, only on rotations in spin-space (the generalization to real-space is straight-forward). Then we first choose a basis of operators in 2nd quantization, let's say $c_{i, \uparrow}, c_{i, \downarrow}$ and their conjugates, which maintain the canonical anti-commutation relations
$$ \{ c_{i,\sigma}, c^{\dagger}_{j,\sigma '} \} = \delta_{i,j} \delta_{\sigma, \sigma '}$$
A rotation in spin-space will be a linear transformation within this basis
$$ R_{S_i} c_{i,\sigma} R^{-1}_{S_i} = a_r c_{i,\sigma} + b_r c_{i, \bar{\sigma}} $$
where $\bar{\sigma}$ is the opposite spin to $\sigma$. How to carry out this rotation is straight-forward from the construction of the spin operators in second quantization 
$$ S_{i}^{\alpha} = \frac{\hbar}{2} \sum_{\lambda, \lambda '} c^{\dagger}_{i,\lambda} \sigma^{\alpha}_{\lambda, \lambda '}c_{i,\lambda '}$$
with $\alpha = x,y,z$ for different spin operations at the point $i$, and $\sigma^{\alpha}$ is the respective Pauli matrix. From this, you can explicitly calculate how does each operator transform, as $R_{S_i} = e^{-i \sum_\alpha \theta_\alpha S^{\alpha}/\hbar}$.
The singlet you wrote is the only rotationaly invariant state for two spins. The point is that in second quantization we don't restrict ourselves from the beginning to a fixed number of particles, and allow them to change. We can construct, however, operators that are invariant under local spin rotations. For example $c^{\dagger}_{i,\uparrow}c_{i, \uparrow} + c^{\dagger}_{i,\downarrow}c_{i, \downarrow}$ is such an operator, as it just counts the number of particles at site $i$, and their rotation is immaterial to it.
A: Based on yu-v's answer, consider a simple rotation around z-axis by /2 and so let
$$R_S = e^{-i \pi/4 c^\dagger_\alpha \sigma_{\alpha \beta} c_\beta}$$
The answer we expect is $S_x\rightarrow S_y$ or $1/2 (c^\dagger_\uparrow c_\downarrow + c^\dagger_\downarrow c_\uparrow ) \rightarrow 1/2 (-i c^\dagger_\uparrow c_\downarrow + i c^\dagger_\downarrow c_\uparrow )$
To carry out this calculation, we need to expand the exponential as
$$R_S = 1 -i \pi/4 c^\dagger_\alpha \sigma_{\alpha \beta} c_\beta + (-i \pi/4 c^\dagger_\alpha \sigma_{\alpha \beta} c_\beta)^2 +...\\
 = 1 - i \pi/4 (c^\dagger_\uparrow c_\uparrow - c^\dagger_\downarrow c_\downarrow )  - \pi^2/32 (c^\dagger_\uparrow c_\uparrow + c^\dagger_\downarrow c_\downarrow+ 2 c^\dagger_\downarrow c^\dagger_\uparrow c_\downarrow c_\uparrow )+...$$
I see to be a faithful representation of the spin operator, we need the single occupancy constraint so the quarctic term should be zero, but there are extra quadratic terms that are higher order with coefficient of $\pi^n$ in n-th order expansion and I'm not sure how to simplify it.
