Dear all, Is there any way to modify the LJ potential as a time dependent value, e.g. V(r)=4*e(t)*((sigma/r)^12-(sigma/r)^6), here e(t) is time-dependent, for instance, e(t)=e0*cos(w*t)? If the code needs to be modified, which C file should I find? Xinghua