I'm confused about the number of degrees of freedom in General Relatity. There are two ways to count it. However, they are contradictory. For simplicity, we consider vacuum solution.
First, $G_{\mu\nu}=0$ gives $10$ equations and $g_{\mu\nu}$ have $10$ degrees of freedom(d.o.f). While $\nabla^\mu G_{\mu\nu}=0$ are $4$ identities, only $6$ of original $10$ ($G_{\mu\nu}=0$) equations are independent. So now there are $6$ independent equations and $10$ degrees of freedom. Because there is coordinate transformation freedom, $4$ of $10$ are gauge freedom. Given $4$ coordinate conditions, only $6$ are physical degrees of freedom and there are $6$ independent equations, so it is well-defined.
Second, if we consider Hamiltonian formalism, we need ADM decomposition where $N$ and $N^i$ are Lagrangian mutipliers and can be given arbitrarily. So there are $6$ degrees of freedom which are $g_{ij}$. While $N$ and $N^i$ give $4$ constraints, only $6-4=2$ are physical degrees of freedom.
Therefore there is contradiction about the number of d.o.f between Einstein Field Equations which give $6$ and Hamiltionian formalism which gives $2$.
So I have following questions:
1）How to reconcile above contradiction? How much physical degrees of freedom are in GR?
2) There is a saying that massless spin-2 particle have two d.o.f. Is this d.o.f same as the physical d.o.f in $g_{\mu\nu}$?
3) We always say that because of the freedom of coordiante transformation (or gauge freedom) ,$x^\mu \xrightarrow{} x^\mu+\xi^\mu(x)$ , we can decrease $4$ d.o.g in $g_{\mu\nu}$. However, in electrodynamics, $A^\mu$ also has gauge freedom $A^\mu \xrightarrow{} A^\mu + \partial^\mu \Lambda(x) $ which can decrease $2$ d.o.f in $A^\mu$. I know how to derive these. I just want to know why in GR $4$ functions freedom $\xi^\mu(x)$ can decrease 4 d.o.f , while in electrodynamics $1$ function freedom $\Lambda(x)$ can decrease $2$ d.o.f .
Remark: I heard someone saying that in first case, $g_{\mu\nu}$ have $10$ d.o.f, gauge freedom decrease $4$ d.o.f, and Bianchi's identities $\nabla^\mu G_{\mu\nu}=0$ are $4$ constraints which decrease $4$ d.o.f of $g_{\mu\nu}$ , so $10-4-4=2$ physical d.o.f are left. I think it is not right, because identity is different from constriant. Identity decreases the number of idependent equations, while constriant decreases the number of d.o.f. Because given any $g_{\mu\nu}$, $\nabla^\mu G_{\mu\nu}=0$ are always right, $\nabla^\mu G_{\mu\nu}=0$ will not have any constraint on $g_{\mu\nu}$.
Best Answer
The key point in all of this is that general relativity is a gauge theory, and, as the saying goes, "the gauge always hits twice" (apparently attributed to Claudio Teitelboim). What this means is that (1) you have an arbitrary freedom in defining your evolution, corresponding to the ability to make gauge transformations, and (2) some of the evolution equations will be constraints. This second fact means that you are not allowed to choose arbitrary initial data for your theory; rather, the initial data that you pick is subject to the constraints, which arise since your action is gauge invariant.
It's usually easiest to start with vacuum electrodynamics. There the equations of motion read $$\partial^\mu(\partial_\mu A_\nu - \partial_\nu A_\mu)=0.$$ Not all of these equations are second order in time; just look at the $\nu=0$ component: $$\partial_t^2 A_0 - \nabla^2A_0 -\partial_t(\partial_t A_0 - \nabla\cdot\vec{A}) = 0 \\ \implies\partial_t\nabla\cdot\vec{A}-\nabla^2A_0 = 0.$$
This is basically the $\nabla\cdot \vec{E} = 0$ vacuum Maxwell equation (i.e. Coulomb gauge with $\nabla\cdot\vec{A}=0$ and $\vec{E} = -\nabla A_0$). This is a constraint on your initial data, because you are not allowed to make an arbitrary choice for $(A_0, \vec{A})$ and $(\partial_t A_0, \partial_t \vec{A})$; rather, they need to satisfy this constraint. So this cuts down the number of initial conditions from 4 to 3. Then the gauge transformation $A_\mu \mapsto A_\mu + \partial_\mu \lambda$ allows you to cut off another piece of initial data, by imposing a gauge fixing condition (i.e. $\nabla\cdot\vec{A}=0$). This brings us to 2 degrees of freedom.
For general relativity, you now have 4 gauge freedoms generated by diffeomorphisms described by a vector $\xi^\mu$. So applying the maxim, we should expect to cut down $4\times2=8$ degrees of freedom. In fact the Bianchi identity tells where to look for the constraints. Let's expand it out a bit: $$0=\nabla_\mu G^{\mu\nu} = \partial_0 G^{0\nu}+\partial_i G^{i\nu} + \Gamma^\mu_{\mu\alpha}G^{\alpha \nu}+ \Gamma^{\nu}_{\mu\alpha}G^{\mu\alpha}.$$ This tells us that the first time derivative ($\partial_0$) of $G^{0\nu}$ is related to spatial derivatives of $G^{i\nu}$ as well as terms with no derivatives of $G^{\mu\alpha}$. The important thing here is that this is an identity, so it holds even if you don't impose the vacuum Einstein equations $G^{\mu\nu}=0$. The tensor $G^{\mu\nu}$ has two derivatives of the metric in it. But if $G^{0\nu}$ had two time derivatives appearing, there would be no way to satisfy the Bianchi identity because no other term in the identity has three time derivatives acting on the metric. This means $G^{0\nu}$ are not evolution equations--they involve only one time derivative of the dynamical variables, and thus are initial value constraints. So that kills 4 degrees of freedom, and you kill 4 more from gauge fixing. This is how you get the $10-4-4=2$ degrees of freedom in general relativity.
And in regards to your second question, yes general relativity describes the two degrees of freedom of a massless spin-2 particle.