Antw: Re: [gmx-users] total charge of the system
Emanuel Peter
Emanuel.Peter at chemie.uni-regensburg.de
Tue Jan 25 11:03:13 CET 2011
Hello,
You have to look into your ions.itp which is included in your top-file by
#include ions.itp.
There all the types of ions have to be defined.
The atom-types which you can see in the ions.itp, you will
find in the ff"your-forcefield".itp where your atomtypes are defined.
On top of your .top file the Atomtype-itp-file is included.
All these files are normally placed in /usr/share/gromacs/top/.
But you also can place them into your current directory.
Bests,
Emanuel
>>> Mark Abraham 25.01.11 10.48 Uhr >>>
On 01/25/11, ahmet yıldırım wrote:Dear Mark,
I looked at gromacs mail list but I could not find a proper solution .What
should I add to the .top file? Please look at the following reconstructed1 .top
and reconstructed1 .top files
I have error as the following reconstructed1 .top file:
Fatal error:
moleculetype CU1 is redefined
I have error as the following reconstructed2 .top file:
Fatal error:
No such moleculetype CL-
I don't have any knowledge of the context, so can't answer. It looks to me like
you are mixing copies of ions.itp from multiple sources. Don't. Use the one for
the force field you are targetting. pdb2gmx generated the right invocation -
all you should have to do is use that by generating correctly-named ions. See
http://www.gromacs.org/Documentation/Gromacs_Utilities/genion
Mark
Original .top file:
; Include topology for ions
#include "gromos43a1.ff/ions.itp"
[ system ]
; Name
GP41 MPER-DERIVED PEPTIDE; ANTI-HIV-1 ANTIBODY 2F5 LIGHT CHAIN; ANTI-HIV-1
ANTIBODY 2F5 HEAVY CHAIN in water
[ molecules ]
; Compound #mols
Protein_chain_P 1
Protein_chain_L 1
Protein_chain_H 1
SOL 10
SOL 127
SOL 157
SOL 41779
reconstructed1 .top file
; Include topology for ions
#include "gromos43a1.ff/ions.itp"
#include "ions.itp"
[ system ]
; Name
GP41 MPER-DERIVED PEPTIDE; ANTI-HIV-1 ANTIBODY 2F5 LIGHT CHAIN; ANTI-HIV-1
ANTIBODY 2F5 HEAVY CHAIN in water
[ molecules ]
; Compound #mols
Protein_chain_P 1
Protein_chain_L 1
Protein_chain_H 1
SOL 10
SOL 127
SOL 157
SOL 41771
CL- 8
reconstructed2 .top file
; Include topology for ions
#include "gromos43a1.ff/ions.itp"
[ system ]
; Name
GP41 MPER-DERIVED PEPTIDE; ANTI-HIV-1 ANTIBODY 2F5 LIGHT CHAIN; ANTI-HIV-1
ANTIBODY 2F5 HEAVY CHAIN in water
[ molecules ]
; Compound #mols
Protein_chain_P 1
Protein_chain_L 1
Protein_chain_H 1
SOL 10
SOL 127
SOL 157
SOL 41771
CL- 8
25 Ocak 2011 10:49 tarihinde Mark Abraham <mark.abraham at anu.edu.au> yazdı:
On 01/25/11, ahmet yıldırım <ahmedo047 at gmail.com> wrote:Hi,
In my simulation, total charge of the system is a noninteger number (System has
non-zero total charge: 8.000004e+00). I neutralized it with 8 chlorine atoms.
Then, grompp -f em.mdp -p topol.top -c solvated.gro -o em.tpr
Fatal error:
moleculetype CU1 is redefined
is it some thing wrong?
ions.itp defines molecule types for ions. Molecule types cannot be redefined.
When you #included ions.itp GROMACS thought you were doing illegal
redefinitions. Look back in the .top to find the original definitions, and then
take suitable action.
Mark
Below is the first and final version of the .top file:
First topol.top File
....
[ molecules ]
; Compound #mols
Protein_chain_P 1
Protein_chain_L 1
Protein_chain_H 1
SOL 10
SOL 127
SOL 157
SOL 41779
Final topol.top File
#include "ions.itp"
[ molecules ]
; Compound #mols
Protein_chain_P 1
Protein_chain_L 1
Protein_chain_H 1
SOL 10
SOL 127
SOL 157
SOL 41771
CL- 8
--
Ahmet YILDIRIM
--
gmx-users mailing list gmx-users at 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-request at gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
--
Ahmet YILDIRIM
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20110125/7d5ba124/attachment.html>
More information about the gromacs.org_gmx-users
mailing list