Dear Jason,
Thank you very much for your kind reply.
>You don't give enough information here. How exactly did you modify
>NUMBER_EXCLUDED_ATOMS and EXCLUDED_ATOMS_LIST? Do you know what the numbers
>in each section mean (including the ordering)?
Sorry not to show how I modified topology file.
I showed the way in the E-mail in mailing list I sent to Mr. Simmerling,
please see it.
>you will probably have to adjust at least one of the pointer values
>in the first section (%FLAG POINTERS).
I changed it, but sander still doesn't recognize it...
Best Regards,
Akihiro Yokota
> Dear AMBER Users,
>
> I am studying the topology file of AMBER to modify it,
> especially for "excluded atoms" declaration. And I found
> a mysterious things below. Do you have any suggestions or
> comments to me?
>
> Let me explain you in a simple case.
>
> 1. An original pdb file have only 2 water molecules, separated in 4A.
> 2. I changed the topology file, to ignore the nonbonded interactions
> between 2 water molecules.
> -> 3 points are modified(NEXT number in 8th line,
> %FLAG_NUMBER_EXCLUDEDATOMS,
> and %FLAG_EXCLUDED_ATOMS_LIST)
3. I checked the interaction energies of modified one by using Anal module.
> Anal said that nonbonded interaction between 2 waters is zero. This is
> OK.
> 4. And optimization by Sander module is performed, but the results from the
> original and the modified are definitely the same. Why does my
> "excluded atoms" declaration not work in Sander?
> (This indicate that Anal module correctly recognizes the "excluded atoms"
> declaration in topology file, and the sander module does not.)
>
You don't give enough information here. How exactly did you modify
NUMBER_EXCLUDED_ATOMS and EXCLUDED_ATOMS_LIST? Do you know what the numbers
in each section mean (including the ordering)? Furthermore, if you want to
adjust the number of excluded atoms, you will probably have to adjust at
least one of the pointer values in the first section (%FLAG POINTERS).
Sander does recognize excluded atoms, which is why a simulation of a free
O2 will have no non-bonded energy contributions, so there has to be
something in your setup that is not quite correct. Keep in mind that the
prmtop is not intended to be modified by hand, and great care is needed when
you try this. In this case, it is very helpful to consult the source code.
For instance, look at the rdparm.f file, where it reads the
EXCLUDED_ATOMS_LIST (search for that phrase). It will only read in the
first "nnb" terms. Therefore, even if you added more terms to this list,
unless you also changed the nnb pointer (NEXT according to the amber file
formats, or the 11th pointer in the POINTERS block), sander will not use
them.
Hope this helps,
Jason
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Nov 09 2010 - 17:00:03 PST
