Re: Minimization Question

From: Gerard Rowe (
Date: Mon Nov 12 2018 - 13:08:00 CST

I had similar issues with PM7 QM/MM calculations with my model system. I don't know that it's particularly robust when dealing with open shell systems with unusual metal geometries. I frequently had SCF failures with a high-spin sextet system, and things weren't much better when I simplified to a low-spin doublet. As far as restraints go, you might want to inspect the force constants produced by the PM7 calculation. If they are unreasonably high, they can overcome the restraint placed on the atoms unless your restraint barrier multiplier is also unreasonably high.

Unfortunately, my only recourse was to fall back to DFT. Is it possible to substitute your copper with zinc in the semiempirical calculations? They're similar as long as you don't need to monitor redox behavior.

From: <> on behalf of McGuire, Kelly <>
Sent: Sunday, November 11, 2018 8:46:19 PM
Subject: namd-l: Minimization Question

  Usually, during a normal MM simulation, I run the minimization for 1,000 to 2,000 steps. Now, with a QM/MM simulation, using PM7, MOPAC, and GPUs, I can minimize for 1,000 steps in 25 minutes. However, looking at the minimization trajectory, in the QM region some of the sidechain atoms move drastically, causing their bond to stretch really far. I tried restraining the whole sidechain, but the it still happens.

1) Should I not use restraints on the atoms in the QM region?

2) Should I use fewer steps during minimization?

3) Could it be PM7, maybe use a non-semi-empirical theory for minimization?

Same questions for the annealing and equilibration steps...

Kelly L. McGuire

PhD Scholar


Department of Physiology and Developmental Biology

Brigham Young University

LSB 3050

Provo, UT 84602

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2018 - 23:21:32 CST