We consider Hardy spaces of functions harmonic on smooth domains in Euclidean spaces of dimension greater than two with respect to the Laplacian perturbed by lower order terms. We deal with the gradient and Schrödinger perturbations under appropriate Kato conditions. In this context we show the usual correspondence between the harmonic Hardy spaces and the $L^{p}$ spaces (or the space of finite measures if p = 1) on the boundary. To this end we prove the uniform comparability of the respective harmonic measures for a class of approximating domains and the relative Fatou theorem for harmonic functions of the perturbed operator.