[gmx-users] qm/mm LA not found
sylwia.kacprzak at physchem.uni-freiburg.de
Mon Nov 7 10:04:57 CET 2016
thank you for your hint. There was no LA entry in the ff*nb.it file. I
have added it there.
But now I get a new error. I have a protein with cofactor. Part of the
cofactor is in QM part
and the rest in MM. For the cofactor I have a separate *itp file as for
the rest of protein.
By preparing the calculation I get and error for all MM atoms:
ERROR 1 [file DRGGMX_groups_v3_QMMM.itp, line 133]:
No default Bond types
example from *itp file
[ bonds ]
; ai aj fu c0, c1, ...
21 22 5 ; C9A N10
22 23 5 ; N10 C1' ;'
24 23 1 ; C2' C1'
24 25 1 ; C2' O2'
There is no error for atoms with bond type 5, but for all with bond type 1.
On 11/05/2016 03:40 PM, Groenhof, Gerrit wrote:
> Is the LA also defined for the force field you're using? Can you check whether there is a LA entry in the ff*nb.it of the force field? If not, add it there manually.
> Message: 2
> Date: Fri, 4 Nov 2016 16:10:03 +0100
> From: Sylwia Kacprzak <sylwia.kacprzak at physchem.uni-freiburg.de>
> To: gromacs.org_gmx-users at maillist.sys.kth.se
> Subject: [gmx-users] qm/mm LA not found
> Message-ID: <581CA4CB.2050002 at physchem.uni-freiburg.de>
> Content-Type: text/plain; charset=utf-8; format=flowed
> Dear all,
> I have set qm/mm calculations with Gromacs 5.1.4 and Orca.
> I treat first the full cofactor on a QM level. However when I cut the
> bond to take only part of it for QM and
> introduce link atom according to the manual with link atom (LA) I get
> following error.
> Program gmx grompp, VERSION 5.1.4
> Source code file:
> /home/widu/gromacs-5.1.4/src/gromacs/gmxpreprocess/toppush.c, line: 1345
> Fatal error:
> Atomtype LA not found
> Any ideas?
More information about the gromacs.org_gmx-users