4Currently only response to an infinitesimal field at zero field is implemented in CASTEP.