[gmx-users] QM/MM error: no such bonded interactions with 5 atoms

[Mark Abraham]

Hi,

CHARMM27 uses CMAP dihedrals that have 5 atoms. How that relates to your
topology isn't clear on the information you've provided. In what context
did you get the error message?

Mark

On Thu, Mar 10, 2016 at 12:58 PM Roman Zeiss <romanzeiss at gmx.net> wrote:

> Dear Gromacs users,
>
> I am trying to set up a QM/MM simulation on Gromacs 5.0.4 with Charmm27. I
> get the error "no such bonded interactions with 5 atoms" when I try to run
> the simulation. What I did so far:
>
> 1. I modified the force field files atomtypes.atp and ffnonbonded.itp to
> include the atom type "LA" and also added a residue file "la.rtp".
>
> 2. In my topology file for the Protein I have (after adding the LA atoms
> to that chain at position 1763 and 1764):
> [ virtual_sites2 ]
> ; LA QMatom MMatom func length
> 1763    30  32  1   0.65
> 1764    1747    1745    1   0.65
> [ constraints ]
> ; QMatom MMatom func length
> 30  32  2   0.153
> 1747    1745    2   0.153
>
> 3. Then I created an index.ndx with a group [ QMatoms ].
>
> 4. For all bonds between QMatoms I changed the bond type from 1 to 5.
>
> Do I have to change anything more like dihedrals or something? Maybe
> delete bonds/dihedrals containing the connecting atoms? But nothing so far
> changed the error message.
>
> I appreciate any kind of advice very much!
> Thanky you and best regards,
> Roman
>
>
>
>
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-request at gromacs.org.