I am trying to derive the weak-field Schwarzschild metric, but starting from the same form as Schwarzschild:
$ds^2=-(1+2\Phi(r))dt^2+(1-2\Psi(r))dr^2 +r^2 d\Omega^2$
which has $R=-2\partial_r^2 \Phi(r)$ and $R_{00}=\partial_r^2\Phi(r)$. When I use $T^0_0=\rho$, then I get to the equation
$T^0_0=R^0_0 -\frac{1}{2} \delta^0_0 R = -\partial_r^2\Phi(r)+\partial_r^2 \Phi(r)=0$
i.e. $\rho=0$. Is it impossible to derive taking a non-vacuum solution, or have I made a mistake somewhere?