I am working with a system that has 13 chlorophylls.

I am trying to do a QM/MM Md calculation with orca, where I only want to set one chlorophyll without the tail part in the qm region. By using gaussview and vmd I could select the appropriate atoms for this region, but orca shows an error which says " MM1 atom 2705 (bonded to QM atom 2593) is a hydrogen atom. Cannot cut boundaries with hydrogen atoms, Please include atom 2705 in the QM region"

however there shouldn't be covalent bond between these atoms.

if you visualize these two atoms the distance between them is almost 9 angstrom.

I included the input and the coordinates that I used.

Thanks for any help anyone can give!

More Zahra Tohidi Nafe's questions See All
Similar questions and discussions