From the dump file of the simulation and by the “rerun” command I want to calculate the bond orientational order of atoms by the “orienorder/atom” command or the “RDF” or the “MSD” of the system!
you can use the pair_style zero (https://docs.lammps.org/pair_zero.html) since you don’t need to reevaluate the potential to compute these three properties.
HTH
Evangelos
Στις Δευ 18 Απρ 2022 στις 8:49 π.μ., ο/η Lamm Gro <[email protected]> έγραψε:
What the command “pair_style zero” does is to set all non-bonded interactions equal to zero and use the specified cutoff to build the neighbor lists and acquire ghost atoms
This is fine for computing the msd, rdf and orientation parameters in a post-processing step. These quantities do not need the value of the potential energy.
Also, what you observe for the pe values is absolutely expected and documented: since the non-bonded interactions are zero, the pe/atom is zero and the sum is also zero.
Personally, I do not have any example file to deal with simulations using the “fix external” command in a post-processing step as you now ask for.
Στις Τρί 19 Απρ 2022 στις 8:23 π.μ., ο/η Lamm Gro <[email protected]> έγραψε:
It is difficult to make specific suggestions without having any information about how fix external is used (there are two modes) and whether it is being called from C or from python.