 Solution

Retain the homogeneous lapse function $\bar{N}(t)$, and define the physical Hubble parameter $H=\dot{a}/(\bar{N}a)$. At first order the scalar metric perturbations give $\delta g_{00}=-2{\bar{N}}^2\Phi$, $\delta g_{0i}=a^2{B}_{,i}$, and $\delta g_{ij}=-2a^2(\Psi {\delta }_{ij}+E_{,ij})$. Applying the passive coordinate displacement component by component gives the scalar gauge transformations with a background lapseFor example the $0i$ transformation is $\delta {\tilde{g}}_{0i}=a^2{\partial }_{i}B-a^2{\partial }_i\dot{\lambda }+{\bar{N}}^2{\partial }_i{\xi }^0$. The spatial isotropic part acquires $-2a\dot{a}{\xi }^0{\delta }_{ij}$, explaining the positive sign in the transformation of $\Psi$.

Preserving synchronous gauge requires both transformed $\Phi$ and $B$ to remain zero. The first condition is ${\partial }_t(\bar{N}{\xi }^0)=0$, and the second is $\dot{\lambda }={\bar{N}}^2{\xi }^0/a^2$. Integrating them givesThus fixing the lapse function and shift vector perturbations does not exhaust the scalar coordinate freedom. The functions $C,D$ describe residual synchronous-gauge freedom; an integration-origin change in the integral is absorbed into $D$.

To move from synchronous gauge to Newtonian gauge, choose $\lambda =-E_S$ to set $\tilde{E}=0$. Requiring $\tilde{B}=0$ then givesA scalar energy density transforms as $\delta {\rho }_N=\delta {\rho }_S-\dot{\bar{\rho }}{\xi }^0$. Dividing by the background energy density, which changes the contrast only at second order if perturbed in that denominator, yields the synchronous-to-Newtonian density transformation with a lapseThe last form uses background stress-energy conservation, $\dot{\bar{\rho }}=-3\bar{N}H(\bar{\rho }+\bar{P})$. With cosmic time $\bar{N}=1$ this is ${\delta }_N={\delta }_S-3(1+w)a^{2}H{\dot{E}}_S$; with conformal time $\bar{N}=a$ it is ${\delta }_N={\delta }_S-3(1+w)\mathcal{H}{E}_S^{\prime }$, where $\mathcal{H}=a^{\prime }/a$. Keeping these clock conventions distinct is essential to the factors of $a$.