def METHOD_NAME(self):
    self.system.non_bonded_inter[0, 0].lennard_jones.set_params(
        epsilon=1.0, sigma=1.0, cutoff=2**(1.0 / 6.0), shift="auto")
    self.system.integrator.set_steepest_descent(
        f_max=1, gamma=0.01, max_displacement=0.01)
    self.system.integrator.run(100)
    self.system.integrator.set_vv()
    self.system.non_bonded_inter[0, 0].lennard_jones.set_params(
        epsilon=0.0, sigma=1.0, cutoff=2, shift="auto")