Thanks both Mark and [email protected] 's advises. Listening to these
words of yours,is better than learning for ten years. Most people uses the
AMBER force field to run HEME simulation,I should focus on suitable force field
and associated correct parameters.thanks
> 1. Re: HEME-cysteine gromacs simulation
> 2. HEME-cysteine gromacs simulation ([email protected])
>
>Message: 1(Mark Abraham)
>Date: Wed, 26 Aug 2009 12:46:24 +1000
>From: Mark Abraham <[email protected]>
>Subject: Re: [gmx-users] HEME-cysteine gromacs simulation
>
>Visualization programs have to use heuristics to guess where bonds are
>and what names occur for which atoms (HG1 might be mercury or the first
>hydrogen on a gamma carbon, NE might be neon or an epsilon nitrogen,
>etc.). Only if you are able to read in a .tpr (such as with ngmx
>distributed inside GROMACS) can you get any visual confirmation of what
>GROMACS thinks you have described in your topology. There is no
>standardized widely-used mechanism for including such information in a
>PDB file, and so such usage of a file is unreliable.
>
>
>Read the .log files as well. GROMACS almost never segfaults (except when
>subjected to buggy MPI libraries), and when it does it will normally
>write some diagnostic information to stdout or the .log file.
>
>Sorry, I can't help there.
>
>> please help me how handle those problems,how can i go on the simulation?
>Mark
>
>Message: 2
>Date: Tue, 25 Aug 2009 23:10:17 -0400
>From: [email protected]
>Subject: [gmx-users] HEME-cysteine gromacs simulation
>You always need a correct topology. The main issue here is that you
>need to have correct parameters. Where did you get your heme and Fe
>parameters? Were you careful about the Fe state oxidation state? I
>suspect that most people use Amber because of their antechamber
>program, which seems like a brilliant idea even if it may overstretch
>it's own parameterization without letting you know. I have no specific
>advice for you here beyond saying that it is worth spending a month
>figuring out what parameters you should really be using (hint: why
>would you be using ffG43a1 if most others use Amber here? do you know
>something that they don't?). This will end up saving you time in the
>long run.
>
>Note that ffG43a1 is the GROMOS forcefield, not the GROMACS
>forcefield. It is unfortunate that many programs (Amber, Charmm,
>gromacs) have their own similarly named forcefield, but that does not
>mean that the forcefield must be used with the associated program.
>There are forcefields that are not associated with a program (OPLS)
>and programs that never developed their own force field (NAMD,
>Desmond, Tinker, LAMMPS), so it is perfectly ok for you to use the
>amber forcefield with the gromacs program.
>
>Bottom line: read, read, read.
>
>Chris.
_______________________________________________
gmx-users mailing list [email protected]
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [email protected].
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
