All, I have a system of nitric acid modelled as NO3- and HO3+ plus SPC water. 
After setting the system up using genbox I then energy minimise. Althouth the 
input conf.gro file looks fine (NO3- labelled 'NO3', HO3+ labelled 'HO3' and 
water labelled 'SOL') the confout file has all of the H3O renamed to 'SOL'. 
This makes it difficult to analyse my system if I can't distinguish between 
water and hydronium.

To construct the system I started with a single NO3- and a single HO3+ in 
separate gro files. I used the following commands in this order:

genbox -cs NO3.gro -nmol 124 -box 4 4 4 -o NO3conf.gro
genbox -cp NO3conf.gro -ci hydronium.gro -nmol 124 -o NO3+H3Oconf.gro
genbox -cp NO3+HO3conf.gro -cs -o conf.gro

This results in the following order in the gro file: NO3, HO3, SOL (I manually 
edit the topology file so that the number of water molecules matches those 
added by genbox).
This is the same as the order in the topology file. Residues in the topology 
file have the same name as in the original gro files.

I should point out that I'm using my own force field file, based on information 
on the literature, with the parameters for SPC water taken from the gromacs 
AMBER force field (all the literature data is based on AMBER). I don't see that 
this would affect residue naming though.

I'm using gromacs v 4.5.4

If anyone can shed some light on this, or suggest a work around, I'd appreciate 
it.

Laura
--
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