The behavior of the resolvent at low energies has implications for many kinds
of asymptotics, including for the scattering matrix and phase, for the
Dirichlet-to-Neumann map, and for wave evolution. In this paper we present a
robust method, based in part on resolvent identity arguments following Vodev
and boundary pairing arguments following Melrose, for deriving such expansions,
and implement it in detail for compactly supported perturbations of the
Laplacian on . We obtain precise results for general self-adjoint
black box perturbations, in the sense of Sj\"ostrand--Zworski, and also for
some non-self-adjoint ones. The most important terms are the most singular
ones, and we compute them in detail, relating them to spaces of zero
eigenvalues and resonances.