Dear LAMMPS users ,
Greetings .
I want to calculate radial distribution function from the center of mass for fullerene ( C60 ) .
Is it possible with LAMMPS ?
Can I use " compute com " command ?
http://lammps.sandia.gov/doc/compute_com.html
Thanks for your attention .
Regards,
Saeed .
I forgot to say that I have fullerene inside of water molecules and I want to calculate RDF from center of mass of fullerene to oxygen atom of water molecules .
2015-04-01 11:49 GMT+02:00 SaeeD PourasaD <[email protected]...>:
I forgot to say that I have fullerene inside of water molecules and I want
to calculate RDF from center of mass of fullerene to oxygen atom of water
molecules .
Dear LAMMPS users ,
Greetings .
I want to calculate radial distribution function from the center of mass
for fullerene ( C60 ) .
There's the compute rdf command, but it only calculated the "real" rdf
between all pairs in the specified groups. I think you're best off by
either writing your own compute (based on the existing compute rdf), our do
it in post-processing of dump files.
Kristof
saaed,
you are running the C60 as rigid objects, yes?
then you could add an atom with a tiny mass (say 1.0e-10) with zero
potential energy in the center of each C60 and then compute the g(r)
relative to this reference particle. otherwise you need to post
process.
axel.
Hi .
Greetings.
Thanks alot for you replay.
I have no idea about post-processing so i will check it .
About Axel’s suggestion , that is very nice idea i will do that ^^
Regards,
Saeed.
A conceptual note. The comptue rdf command
calculates an RDF between pairs of atoms.
There is no command in LAMMPS to compute
RDFs between arbitrary point(s) and atoms.
I suppose one could be written to infer points
from the COM or geometric center of a molecule
and compute an RDF between those points and
all other atoms. The new compute com/chunk
would be a starting point, since it shows how
to calculate the COM of each chunk and chunks
can be molecules.
So a new compute rdf/chunk would be possible.
Steve
Hi .
Dear LAMMPS users ,
I have a C60 inside of plenty of water molecules.
I want to calculate COM radial distribution function. I used Axel’s suggestion ;I add an atom with a tiny mass (say 1.0e-10) with zero
potential energy in the center of C60 .
Unfortunately the result is not OK , and I get large amount for RDF .
Do you have any suggestion for solving this problem ?
a part of my data file :
4 atom types
3 bond types
3 angle types
1 dihedral types
1 improper types
-50.0 40 xlo xhi
-50.0 40 ylo yhi
-50.0 40 zlo zhi
Masses
1 12.011 #CA “Aromatic C” 6 partial charge=-0.115
2 15.9994
3 1.008
4 1.0e-10 #this is inside of C60 for COM_RDF !
and a part of my input file :
…
group C60_COM type 1 4
…
NVT
fix 1 SPC shake 1.0e-4 500 1000 a 2 b 2
fix 2 SPC nvt temp 298.0 298.0 100.0
fix 4 C60_COM rigid/nvt molecule temp 298.0 298.0 100.0
time step 0.05
run 2000000
…
NPT
fix 5 SPC npt temp 300.0 300.0 100.0 iso 1.0 1.0 1000.0 dilate all
fix 6 C60_COM rigid/nvt molecule temp 300.0 300.0 100.0
compute myRDF all rdf 250 4 3
fix 7 all ave/time 2 500 5000 c_myRDF file COM.rdf mode vector
time step 1.0
run 1000000
Thanks in advance .
Best ,
Saeed.
Dear Axel ,
Do you have any idea about this ? I can not understand the reason of this result .