On Fri, Apr 13, 2018 at 6:23 PM, Chuanxing Zhan <chuanxingzhan@...16...> wrote:
> Thanks for you information! Can you say anything specific about the second
> method for computing molecular g(r)? The question is, how to define this
> COM of molecules to particle and output to trajectory which can use to
> compute the molecular g(r)? I used the example command to compute the COM of
> molecule. If i want co continue, what command would be helpful, dump or
> restart ?
neither. like i wrote, this is a post-processing task. so you need to
write a little program to read the tmp.out file and compute the g(r)
from it yourself.
this is a frequent assignment in graduate level statistical mechanics
classes, and there plenty of examples on the web. you just need to
adjust it for this specific case.
> The example command for computing COM of molecule:
> compute cc1 all chunk/atom molecule
> compute myChunk all com/chunk cc1
> fix 1 all ave/time 100 1 100 c_myChunk[*] file tmp.out mode
> On Friday, 13 April 2018, 11:56:45 AM GMT-5, Axel Kohlmeyer
> <akohlmey@...24...> wrote:
> On Fri, Apr 13, 2018 at 12:40 PM, Chuanxing Zhan via lammps-users
> <email@example.com> wrote:
>> Dear lammps users,
>> Hello, everyone! I'm trying to compute the g(r) of COM of polymer, the
>> lammps command only compute the g(r) between the atoms. I read previously
>> thread in lammps forum, and what they mentioned is I need to build a
>> for calculating COM and defined the molecules. I'm wonder anyone have a
>> better way to do the molecular g(r) calculation. Any help would be
> please keep in mind, that a g(r) is a property which computes the
> relative probability of point particles in a given system versus the
> distribution of an ideal gas.
> so any application to molecular systems is of limited meaning.
> the easiest approach is usually to pick the atom that is closest to
> the center of mass, give it a different atom type number (but the same
> force field parameters) and then compute the g(r) for that type of
> atom. if you average over a large enough system and long enough
> trajectory, the differences to the g(r) of the COM should be
> a more involved approach would be to compute the COM(s) and output it
> to a file and then compute the g(r) in post-processing.
>> Thanks !
>> Check out the vibrant tech community on one of the world's most
>> engaging tech sites, Slashdot.org! http://sdm.link/slashdot
>> lammps-users mailing list
> Dr. Axel Kohlmeyer akohlmey@...24... http://goo.gl/1wk0
> College of Science & Technology, Temple University, Philadelphia PA, USA
> International Centre for Theoretical Physics, Trieste. Italy.
Dr. Axel Kohlmeyer akohlmey@...24... http://goo.gl/1wk0
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.