Differencing estimator in two-way fixed effects model Consider a simple linear panel regression model with two-way fixed effects:
$$Y_{it}=\alpha_i+\mu_t+X_{it}'\beta+u_{it}, \quad i=1,...,n, \ t=1,...,T,$$
where $n$ and $T$ represent the total units in individual and time dimensions, respectively. Variables $\alpha_i$ and $\mu_t $ are unobserved individual and time effects, respectively, which are correlated with (time variant) regressors $X_{it}\in \mathcal{R}^p$. The random noise $u_{it}$ satisfies $E(u_{it}|\alpha_i,\mu_t,X_{it})=0$. When both $n$ and $T$ are large, the least square dummy variable approach becomes infeasible for a consistent estimator of $\beta$. One popular method is the demean approach that wipes out $\alpha_i$ and $\mu_t $.
Question: If we instead use differencing estimator, can the estimator for $\beta$ be still consistent? Namely, we perform differencing on both sides of the model above to have
$$\ddot{\Delta}Y_{it}=\ddot{\Delta}X_{it}'\beta+\ddot{\Delta}u_{it}, \quad i=2,...,n, \ t=2,...,T,$$
where $\ddot{\Delta}Y_{it}=Y_{it}-Y_{i,t-1}-Y_{i-1,t}+Y_{i-1,t-1}$ and similarly for $\ddot{\Delta}X_{it} $ and $\ddot{\Delta}u_{it}$. I am sure that the $\beta$ can be estimated through OLS if we assume $E(\ddot{\Delta}u_{it}|\alpha_i,\mu_t,X_{it})=0.$
However, in such differencing process, $Y_{i-1,t}$ is not uniquely defined, because the order of $i=1,...,n$ can be arbitrary. For instance, if $i=3$ stands for the U.S., then $i-1=2$ can be any other country in a dataset, like China or Canada. I am not sure if this is a problem for the differencing estimator. I found no closed reference from the literature, so any help would be appreciated.
Thanks.

When you difference you throw away your first observation, I think you acknowledge this by letting your differenced model start with T=2,3,...,T. However you also set i=2,3,...,N , this is unusual, you throw away your first observation from a time perspective, this is because you do a time difference operation delta on your data Y. The set of entities i in your differenced data should be the exact same as those you start out with, the time differencing has no impact on your entities i. Hope this makes sense. Maybe I misunderstand you, because the way you phrase your differenced model is unusual to me, typically we represent the differenced model as $\Delta Y_{i,t}=\Delta X_{i,t}'\beta+ \Delta u_{i,t}$.
