The wave function of a free Dirac particle moving with momentum $\vec{p}=(p,0,0)^T$ is given in the rest frame and the laboratory frame as
$$\Psi_r=N_r\left(\begin{array}{c}1\\0\\0\\0\end{array}\right)e^{-\frac{i}{\hbar}mc^2t}\qquad\Psi_l = N_l\sqrt{\frac{E+mc^2}{2mc^2}}\left(\begin{array}{c}1\\0\\0\\\frac{pc}{E+mc^2}\end{array}\right)e^{-ip_\mu x^\mu/\hbar},$$
where $N_r$ and $N_l$ are normalization constants which are
$$N_i =\left\{\begin{array}{l}1/\sqrt{V},\quad i=r\\\sqrt{\frac{mc^2}{EV}},\quad i=l\end{array}\right.$$
and which I calculated using $\int\mathop{d^3x}\Psi_i^\dagger\Psi_i=1$.
Now, I would like to calculate the spin expectation values $\left<\sigma_z\right> = \frac{1}{N}\int\mathop{d^3x}\Psi^\dagger\Sigma_z\Psi$ in the rest and the laboratory frame, where $\Sigma_z = \frac{\hbar}{2}\left(\begin{array}{cc}\sigma_z & 0 \\ 0 & \sigma_z\end{array}\right)$ is the spin operator and $\sigma_z$ the third Pauli matrix. For the rest frame, I get
$$\left<\sigma_z\right>_r = \frac{\hbar\sqrt{V}}{2}\int\mathop{d^3x}(1,0,0,0)e^{i mc^2t/\hbar}\left(\begin{array}{cc}\sigma_z & 0 \\ 0 & \sigma_z\end{array}\right)\left(\begin{array}{c}1\\0\\0\\0\end{array}\right)e^{-imc^2t/\hbar}=\frac{V^{3/2}\hbar}{2}$$
and analogously in the lab frame
$$\left<\sigma_z\right>_l = \frac{\hbar}{2}\sqrt{\frac{E}{mc^2}}V^{3/2}\frac{mc^2}{E+mc^2}.$$
I'm quite sure that these results cannot be true. Especially in the rest frame, I would expect an expectation value which is independent of the volume $V$. So it would be nice if someone could tell me how to do the calculation correctly.