<div>Hi again,<br><br>I have the following:<br><br>I have not changed anything
in ffG53a6.atp and .rtp<br>As suggested by gmx users, I have the following in
the ions,itp file<br><br>[ moleculetype ]<br>
; molname nrexcl<br>
OH 1<br>
<br>
[ atoms ]<br>
; id at type res nr residu name at name cg nr
charge<br>
1 OH 1 OH OH 1 0<br>
1 OW 1 OW OW 1 -1<br><br><br>I executed the
following commands:<br><br>(a) pdb2gmx -f his8.pdb -ignh -ter -his -p topol.top
-o his8.gro<br>(b) editconf -f his8.gro -box 5 5 5 -o his8box.gro<br>(c) genbox
-cp his8box.gro -ci 1methanol.gro -nmol 1800 -o his8_meoh.gro -p
topol.top<br>(d) grompp -f ions.mdp -c his8_meoh.gro -p topol.top -o ions.tpr
-p topol.top<br>(e) genion -s ions.tpr -o added_ions.gro -p topol.top -pname NA
-nname OH -nn 8<br>(f) genbox -cp added_ions.gro -cs
tip4p.gro -o his8_meohi_wat.gro -p topol.top<br> <br>Commands
(a)-(f), all goes smooth.<br><br>When I do the following<br> grompp -f
minim.mdp -c his8_meohi_wat.gro -p topol.top -o em.tpr<br><br>I run into error
which says moleculetype OH not found.<br><br><br>Added the following
specifically in .top file <br> [ moleculetype ]<br>
; molname nrexcl<br>
OH 1<br>
<br>
[atoms]<br>
1
OH 1 OH
OH 1
0.21
1.008 ;<br>
2
OW 1 OH
OW 1
-1.21
15.99994 ;<br>
[ bonds ]<br>
OW HW1 gb_38<br><br>Going back
to the command<br>grompp -f minim.mdp -c his8_meohi_wat.gro -p topol.top -o
em.tpr<br><br>runs into different error as mentioned earlier : Fatal error:
number of coordinates in coordinate file
(5his8_meohi_wat.gro, 12487) does not match topology (topol.top,
12519)<br><br>The .top file includes tip4p.itp, methanol.itp and
ions.itp.<br><br>So I am totally confused as to where I am going wrong.
<br>Please guide.<br><br>Thanks,<br>Shivangi<br><br><br><br>Shivangi
Nangia<br>Garrison Group<br>Graduate Student<br>201A Chemistry
Building<br>University Park PA 16802<br><br>>
<br><br><br><br><br><br>Hello,<br>><br>> Thanks a lot for the
suggestion.<br>><br>> I think I made some progress but I am not there
yet.<br>><br>> I added the following in the ions.itp file<br>><br>>
[ moleculetype ]<br>> ; molname nrexcl<br>>
OH 1<br>><br>> [ atoms ]<br>> ; id at
type res nr residu name at name cg nr charge<br>>
1 OH 1 OH OH 1 0<br>> 1 OW
1 OW OW 1 -1<br>><br>><br>> Still I ran
into the same error that it could not find moleculetype OH.<br>><br>>
Than I specifically added the following to the .top file<br>><br>> [
moleculetype ]<br>> ; molname nrexcl<br>>
OH 1<br>><br>>
[atoms]<br>>
1
OH 1 OH
OH 1
0.21 1.008
;<br>>
2
OW 1 OH
OW 1
-1.21 15.99994 ;<br>> [ bonds
]<br>> OW HW1
gb_38<br>><br>> This helped.<br>> However, there is new error I am
running into which says:<br>> Fatal error: number of coordinates in
coordinate file (5his8_meohi_wat.gro, 12487) does not match topology
(topol.top, 12519)<br>><br>> I am trying to add 8 hydroxide ions in my
methanol-water system.<br>><br><br>And how are you attempting this?
What command(s) did you give? What subsequent modifications to your
topology did you make?<br><br>The generic solution to your problem can be found
here:<br><br>http://www.gromacs.org/Documentation/Errors#Number_of_coordinates_in_coordinate_file_does_not_match_topology<br><br>> The confusion I have is:<br>> The .gro file has OH but is gromacs really accounting OH as 2 atoms?<br><br>Gromacs assigns parameters to the molecules it finds based on the entries in the respective [moleculetypes]. If you tell grompp or any other program that the [moleculetype] "OH" is comprised of two atoms, then there's no reason to believe anything is wrong.<br><br>> Also, the difference (12519-12487) is greater than 8.<br>><br><br>Without some context as to what you're doing and what you've done up to this point, this information does not lead to any useful solution. The number of molecules (either methanol or OH) in your topology does not match the contents of your coordinate file.<br><br>-Justin<br><br>> Kindly help.<br>><br>> Thanks,<br>> Shivangi<br>><br>><br>><br>> Shivangi Nangia<br>!
> Garrison Group<br>> Graduate Student<br>> 201A Chemistry Building<br>> University Park PA 16802<br>><br>> On 22/03/2011 7:44 AM, SHIVANGI NANGIA wrote:<br>> > Hello,<br>> ><br>> > I am new gromacs user.<br>> > I am trying to set up a water-methanol system with few hyrdroxide ions.<br>> ><br>> > Since hyrdroxide ion does not already exist in gromacs database, I defined the following in the .rtp file ( I am using ffG53a6)<br>> ><br>> > [ OH ]<br>> > [ atoms ]<br>> > OW OW -1.41000 0<br>> > HW1 H 0.41000 0<br>> > [ bonds ]<br>> > OW HW1 gb_38<br>> > [ angles ]<br>>&!
nbsp; > ; ai aj a!
k &
nbsp; gromos typpe<br>> > [ impropers ]<br>> > ; ai aj ak al gromos type<br>> > [ dihedrals ]<br>> > ; ai aj ak al gromos type<br>><br>> That's only useful for pdb2gmx when generating .top files.<br>><br>> > I also included OH in the .atp file ( I realise that hydroxide ion is NOT an atom, its a diatomic anion, I was just trying!)<br>><br>> Don't. The .atp defines *atom* types.<br>><br>> > I ran into the error: No such moleculetype OH<br>> ><br>> > I further included information in the ions.itp as follows:<br>> ><br>> > [ moleculetype ]<br>> > ; molname nrexcl<br>> > OH 1<br>> ><br>> > [ atoms ]<br>> !
> ; id at type res nr residu name at name cg nr charge<br>> > 1 OH 1 OH OH 1 -1<br>><br>> Use two atoms, like you did for your .rtp entry.<br>><br>> > The genion adds the "OH" in the .gro file and also the topol.top but when I try to minimize the system using<br>> > grompp -f minim.mdp -c 5his8_meohi_wat.gro -p topol.top -o em.tpr<br>> > I am still running into the same error, No such moleculetype OH.<br>><br>> Your .top still does not have a functional [ moleculetype ] OH. You're probably not using the #include "ions.itp" mechanism correctly.<br>><br>> Mark<br>><br><br>-- ========================================<br><br>Justin A. Lemkul<br>Ph.D. Candidate<br>ICTAS Doctoral Scholar<br>MILES-IGERT Trainee<br>Department of Biochemistry<br>Virginia Tech<br>Blacksburg, VA<br>jalemkul[at]vt.edu | (540) 231-9080<br>http://www.beva!
nlab.biochem.vt.edu/Pages/Personal/justin<br><br>=============!
========
===================<br>-- gmx-users mailing list gmx-users@gromacs.org<br>http://lists.gromacs.org/mailman/listinfo/gmx-users<br>Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting!<br>Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-request@gromacs.org.<br>Can't post? Read http://www.gromacs.org/Support/Mailing_Lists<br><br><br></div>