Parity on gamma matrices I want to understand clearly why $ P \gamma^{\mu} P = \gamma^{\mu} $, where $ P $ is the parity operator. 
This result follow for example from pag. 66 of Peskin-Schroeder. 
The parity operator acts on Dirac fields in this way: $ P \psi (t, \textbf{x}) P = \gamma^0 \psi (t, -\textbf{x}) $, assuming no other phase factors. 
On the Dirac bilinear $ \bar{\psi} \gamma^{\mu} \psi (t, \textbf{x}) $, using the fact that $ P^2=1 $, you have  $P \bar{\psi} \gamma^{\mu} \psi (t, \textbf{x})P = P \bar{\psi} P P \gamma^{\mu} P P \psi (t, \textbf{x}) P =  \bar{\psi}\gamma^0 P \gamma^{\mu} P \gamma^0 \psi (t,- \textbf{x}) = \bar{\psi}\gamma^0 \gamma^{\mu}  \gamma^0 \psi (t,- \textbf{x}) = (-1)^{\mu}\bar{\psi} \gamma^{\mu} \psi (t, -\textbf{x}) $
With $(-1)^{\mu} =1 $ if $ \mu=0 $ and $(-1)^{\mu} =-1 $ if $ \mu=1,2,3 $.
So $ P \gamma^{\mu} P = \gamma^{\mu} $ follows because $ P$ and  $\gamma^{\mu} $ act on different spaces? Or there are other explanations? (Please, do all the steps in the answer)
To be clear: I use $ \gamma^0= \bigl(\begin{smallmatrix}
0&1\\ 1&0
\end{smallmatrix} \bigr)$ and $ \gamma^{i}= \bigl(\begin{smallmatrix}
0&\sigma^i\\ -\sigma^i&0
\end{smallmatrix} \bigr)$. 
 A: One of the ways for getting this result is the trivial consequence that parity operator must change the sign of integral values of 3-momentum and current while it leaves invariant full energy and charge values. For example, for energy density we have with $\Psi ' = \hat {P}\Psi$ ($\hat {P} = \hat {U}P_{\mathbf x \to -\mathbf x}$) the following result:
$$
E \to E' = \Psi{'}^{\dagger}(-\mathbf x, t) \left((\hat {\mathbf p} \cdot \hat {\alpha}) + \gamma_{0}m\right)\Psi{'}(-\mathbf x , t) = 
$$
$$
=\Psi^{\dagger}(\mathbf x , t)\hat {P}^{\dagger}\left((\hat {\mathbf p} \cdot\hat {\alpha} ) + \gamma_{0}m\right)\hat {P}\Psi(\mathbf x,  t) = 
$$
$$
=\Psi^{\dagger}(\mathbf x , t)\left(-\left(\hat {\mathbf p} \cdot \hat {P}^{+}\hat {\alpha}\hat {P} \right) + \hat {P}^{\dagger}\gamma_{0}\hat {P}m\right)\Psi(\mathbf x,  t) = 
$$
$$
=\Psi^{\dagger}(\mathbf x, t) \left((\hat {\mathbf p} \cdot\hat {\alpha} ) + \gamma_{0}m\right)\Psi(\mathbf x , t) \Rightarrow
$$
$$
\hat {U}^{\dagger}\hat {\alpha}\hat {U} = -\hat {\alpha}, \quad \hat {U}^{\dagger}\gamma_{0}\hat {U} = \gamma_{0} \Rightarrow \hat {U}^{\dagger}\gamma^{\mu}\hat {U} = \gamma_{\mu}, \qquad (1)
$$
so your equality is possible only in a form of $(1)$ ($\hat{P}^{\dagger}$ formally coincide with $\hat{P}$).
In the first line I write the expression for the energy density after transformation of the function of the state ($\Psi \to \Psi {'}$). In the second one I have used the expressions $\Psi{'} = \hat {P}\Psi , \Psi{'}^{+} = \Psi^{+}\hat {P}^{+}$; in the third one I have moved $\hat {P}^{\dagger}$ right and $\hat {P}$ left on ($\hat {P}^{\dagger}$ acts on $\hat {\mathbf p}$ by only changing its sign) and after that I have got bilinear forms $\hat {P}^{+}\hat {\alpha}\hat {P}, \hat {P}^{+}\gamma_{0}\hat {P}$. Finally, I have equated this result to the energy without inversion, because energy doesn't change (by its physical meaning) under spatial inversion and I have got $(2)$ by equating corresponding expressions from the last line and from the previous one: $-\left(\hat {\mathbf p} \cdot \hat {P}^{+}\hat {\alpha}\hat {P} \right)$ to $(\hat {\mathbf p} \cdot\hat {\alpha} )$ etc.
The last consequence can be obtained by the following way:
$$
\hat {U}^{\dagger}\gamma_{0}\hat {U} = \gamma_{0} \Rightarrow \hat {U}^{\dagger}\hat{\alpha}\hat {U} = \hat {U}^{\dagger}\gamma_{0}\mathbf {\gamma}\hat {U} = \gamma_{0}\hat {U}^{\dagger}\mathbf {\gamma}\hat {U} = -\gamma_{0} \gamma \Rightarrow
$$
$$
\hat {U}^{\dagger}\gamma_{0}\hat {U} = \gamma_{0}, \quad \hat {U}^{\dagger}\gamma \hat {U} = -\gamma .
$$
