Add cutoff energy. This term can be zero if Ecut is ignored.
Add virial.
Initial commit for calculating LJ potential and forces (not implemented yet). The virial calculation is not very clear!