  We show how to describe the coupling of electrons to non-uniform magnetic
fields in the framework of the widely used norm-conserving pseudopotential
appro ximation for electronic structure calculations. Our derivation applies to
magnetic fields that are smooth on the scale of the core region. The method is
validated by application to the calculation of the magnetic susceptibility of
molecules. Our results are compared with high quality all electron quantum
chemical results, and another recently proposed formalism.
