[gmx-users] QM/MM error: no such bonded interactions with 5 atoms
mark.j.abraham at gmail.com
Thu Mar 10 14:37:36 CET 2016
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?
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,
> Gromacs Users mailing list
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> * 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.
More information about the gromacs.org_gmx-users