@@ -0,0 +1,423 @@
+#! /usr/bin/env hoomd%0A%0Afrom hoomd_script import *%0A%0Ainit.create_random(N=10000, phi_p=0.2)%0Alj = pair.lj(r_cut=2.5)%0Alj.pair_coeff.set('A', 'A', epsilon=1.0, sigma=1.0, alpha=1.0)%0A%0Aall = group.all()%0Aintegrate.mode_standard(dt=0.005)%0Aintegrate.nvt(group=all, T=1.2, tau=0.5)%0A%0Aanalyze.log(filename = 'thermo.log', quantities = %5B'potential_energy', 'kinetic_energy', 'nvt_reservoir_energy','momentum'%5D, period=1000) %0A%0Arun(100e6)%0A

Test script for evaluating momentum conservation