[gmx-users] building stacks with a desired system density
Moeed
lecielll at googlemail.com
Thu May 13 04:41:40 CEST 2010
Dear gmx users,
I tried to generate a stack of molecules using
*genconf_d -f Hexane.gro -nbox 4.0 8.0 8.0 -dist 1.1 0.55 0.55 -o
Hexane-stack.gro*
I ran the MD simulation successfully but now I want to change the number of
molecules so that I get the desired system density. To do so I played a bit
with -dist and -box options but I find it not easy to reach desired
density. I need to have nearly 7 molecules per nm^3.
Regarding genbox program in the literature I see that it is recommended for
solvation of a solute in water. (The reason I want to try this is because it
seems I can directly enter the number of molecules with -nmol option)
1-I am sorry for asking ridiculous questions but does it make sense if I use
genbox and I ask the program solvate a single hexane molecule with hexane
molecules?
genbox –cp Hexane.gro –cs HExane.gro –p Hexane.top –o solvated.pdb
(although this command line is not working I want to know if there is any
way to modify the command to get it working or my idea is totally wrong!)
2- Initialiy one single hexane molecule is in a box of 3*3*3 nm^3
(Hexane.gro). What is the best approach for doing such calculations. Can I
take any arbitrary box size and then based on the box size put a specific
No. of particles in it so that I get density?
Thank you,
Hexane.gro:
**
Glycine aRginine prOline Methionine Alanine Cystine Serine
20
1HEX C1 1 1.786 1.639 1.501
1HEX H11 2 1.881 1.608 1.500
1HEX H12 3 1.768 1.696 1.421
1HEX H13 4 1.770 1.692 1.584
1HEX C2 5 1.693 1.518 1.499
1HEX H21 6 1.711 1.461 1.579
1HEX H22 7 1.710 1.465 1.416
1HEX C3 8 1.547 1.561 1.501
1HEX H31 9 1.528 1.617 1.420
1HEX H32 10 1.530 1.615 1.584
1HEX C4 11 1.454 1.440 1.500
1HEX H41 12 1.472 1.384 1.581
1HEX H42 13 1.472 1.386 1.418
1HEX C5 14 1.307 1.483 1.501
1HEX H51 15 1.289 1.539 1.420
1HEX H52 16 1.289 1.537 1.583
1HEX C6 17 1.214 1.361 1.500
1HEX H61 18 1.119 1.392 1.501
1HEX H62 19 1.231 1.307 1.418
1HEX H63 20 1.232 1.305 1.581
3.00000 3.00000 3.00000 0.00000 0.00000 0.00000 0.00000
0.00000 0.00000
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20100512/7155255f/attachment.html>
More information about the gromacs.org_gmx-users
mailing list