polarisation virial implemented