On 10/31/13 7:03 PM, Ehsan Sadeghi wrote:
Thanks Justin,


Here is the steps that I followed:

pdb2gmx -f eth.pdb -o eth.gro -p eth.top


pdb2gmx -f ionomer-4.pdb -o ionomer4.gro -p ionomer4.top


editconf -f ionomer4.gro -o box.gro -bt dodecahedron -d 0.5


genbox -cp box.gro -cs spc216.gro -ci eth.gro -nmol 1000 -o solvated.gro

I got the box with ionomer in ethanol and water solution. However, I only have 
400 of eth and 423 water molecule. The box wasn't large enough?


Follow the advice at the link I posted before. Trying to add ethanol while simultaneously solvating is a bad idea. If genbox cannot add 1000 ethanol molecules without trying to do anything else, then yes, the box is simply too small. At this point, we can't say because it could simply be that simultaneous solvation might have been competing for available box volume.

What should I do with .top file? I have the .top file for ionomer and .top file 
for one eth molecule.


It is simple to covert between .top and .itp and #include them however you see fit. Probably the conceptually easiest methhod is to consider ethanol just like any other solvent (i.e. water) and #include it within the .top for the ionomer, which is just like any other solute.

-Justin

--
==================================================

Justin A. Lemkul, Ph.D.
Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441

==================================================
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to