Re: [AMBER] masking in the QMMM input

From: Gustavo Seabra <gustavo.seabra.gmail.com>
Date: Fri, 15 Apr 2011 09:45:26 -0300

On Friday, April 15, 2011, Jason Swails <jason.swails.gmail.com> wrote:
> Hello,
>
> On Thu, Apr 14, 2011 at 9:25 PM, setyanto md <stwahyudi.md.gmail.com> wrote:
>
>> Dear amber user and developer
>>
>> I am doin qmmm simulation. I'd like to mask atom rather than residue. Here
>> is list of atoms from my pdb file:
>> :
>> ATOM    935  CG  ASP    61      -4.452  12.837  -8.642  1.00  0.00
>> ATOM    936  OD1 ASP    61      -5.589  12.636  -8.154  1.00  0.00
>> ATOM    937  OD2 ASP    61      -3.962  13.964  -8.866  1.00  0.00
>>
>> I've been read in page 91 manual of Amber 10. we can use iqmatoms to mask
>> our QM region.
>>
>> So I do like this :
>>
>> iqmatoms='.935,936,937',
>>
> No. This is an integer list. Just put "iqmatoms=935,936,937,"

Just to add to Jason's answer, note that the atom numbers that go here
are *not* the atom numbers of the PDB file, but the sequential atom
number from the amber input (prmtop, inpcrd), which is not necessarily
the same as in the PDB file. You can get the correct atom numbers by
opening the amber input in a program such as VMD and checking the
numbering.

Gustavo.

-- 
Gustavo Seabra
Professor Adjunto
Departamento de Química Fundamental
Universidade Federal de Pernambuco
Fone: +55-81-2126-7417
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Apr 15 2011 - 06:00:03 PDT
Custom Search