I need to define the atoms of a ligand which would be in the QM region since the ligand is big. Therefore I defined them as
define res2 sele segid HETB .and. (type C7 .or. type C7M .or. type HC10 .or. type HC9 .or. type HC11 .or. type C8 .or. type C8M .or. type HC12 .or. type HC14 .or. type HC13 .or. type C6 .or. type HC8 .or. type C5X .or. type C9A .or. type C9 .or. type HC7 .or. type N10 .or. type C1' .or. type HC5 .or. type HC6 .or. type N5) end ! QM3 region of substrate define qmsub2 sele segid HETB .and. (resid 2 .and. res2) show end
As it was executed,
HARMM> define res2 sele segid HETB .and. (type C7 .or. type C7M .or. type HC10 .or. type HC9 .or. type HC11 .or. type C8 .or. type C8M .or. type HC12 .or. type HC14 .or. type HC13 .or. type C6 .or. type HC
***** LEVEL 0 WARNING FROM ***** ***** ENDselection terminator missing ****************************************** BOMLEV ( -1) IS NOT REACHED. WRNLEV IS 5
calculation was terminated giving the above error.
Output file was attached for more information.
Is it correct the way I defined atoms? or am I missing something?