Compute potential energy of a pair style not used in integration of equation of motion

Dear Lammps Users,

In my work, I need to run MD according to one Hamiltonian, and sample the average energy according to a different Hamiltonian. In other words, I need to compute the potential energy of a pair style that was NOT used in the integration of motion. I'm wondering whether there is a simple way to achieve this.

Thanks,

Feng

Dear Lammps Users,

In my work, I need to run MD according to one Hamiltonian, and sample
the average energy according to a different Hamiltonian. In other words,
I need to compute the potential energy of a pair style that was NOT used
in the integration of motion. I'm wondering whether there is a simple
way to achieve this.

please check out the rerun command: http://lammps.sandia.gov/doc/rerun.html

axel.

Checkout the rerun command which can do exactly that :slight_smile:

Anders

Sounds like you might be doing free energy perturbation, in which case, it makes sense to also check out the USER-FEP package.