 Linearized cosmological Hamiltonian constraint

For a spatially flat Friedmann-Lemaitre-Robertson-Walker metric, negative extrinsic curvature, and $K=-3H+\kappa$, the Hamiltonian constraint gives the displayed scalar perturbation equation. Indeed $K^2=9H^2-6H\kappa$ and $K^i{}_j{K}^j{}_i=3H^2-2H\kappa$ to first order. The spatial Ricci scalar is $4\Delta \Phi$; subtract the background equation $6H^2=16\pi G\bar{\rho }$ and divide by four.