General constraints for invariance of magnetic properties in a gauge transformation are analyzed. Sum rules relative to the transformation from Coulomb to Landau gauges are examined in particular. Numerical tests for hydrogen fluoride, water, ammonia, and methane molecule have been carried out in large basis set calculations, using random-phase approximation. The conditions for invariance are severe conditions for accuracy of variational molecular wave functions.