Re: QM/MM Question

From: Gerard Rowe (GerardR_at_usca.edu)
Date: Mon Nov 05 2018 - 16:19:08 CST

Unless you have a huge number of atoms, going over 250 iterations isn't useful. If that doesn't work, you may have to walk a long road of trying different convergence aids to achieve a good wavefunction.

________________________________
From: McGuire, Kelly <mcg05004_at_byui.edu>
Sent: Monday, November 5, 2018 2:46:39 PM
To: Gerard Rowe; namd-l_at_ks.uiuc.edu
Subject: Re: QM/MM Question

Excellent, that worked, thanks! In your opinion, what is the maxIter you would try before looking at starting geometry problems or changing things like basis set, level of theory, etc..?

Kelly L. McGuire

PhD Scholar

Biophysics

Department of Physiology and Developmental Biology

Brigham Young University

LSB 3050

Provo, UT 84602

________________________________
From: Gerard Rowe <GerardR_at_usca.edu>
Sent: Monday, November 5, 2018 12:41:11 PM
To: McGuire, Kelly; namd-l_at_ks.uiuc.edu
Subject: Re: QM/MM Question

That should definitely work. Make sure you use something like: qmConfigLine "%%scf maxiter 500 end"

The two percent signs are necessary for the input file to be constructed properly.

________________________________
From: McGuire, Kelly <mcg05004_at_byui.edu>
Sent: Monday, November 5, 2018 1:14:02 PM
To: Gerard Rowe; namd-l_at_ks.uiuc.edu
Subject: Re: QM/MM Question

  My calculations won't converge after 125 iterations. I tried increasing the iterations:

%scf

MaxIter 500

end

    But it says that MaxIter is unrecognized...I followed the manual...is there a different way to increase max iterations for version 4.0.1?

Kelly L. McGuire

PhD Scholar

Biophysics

Department of Physiology and Developmental Biology

Brigham Young University

LSB 3050

Provo, UT 84602

________________________________
From: Gerard Rowe <GerardR_at_usca.edu>
Sent: Monday, November 5, 2018 6:32:06 AM
To: namd-l_at_ks.uiuc.edu; McGuire, Kelly
Subject: Re: QM/MM Question

I don't think Orca has somehting like that. Gaussian has a feature they call ONIOM that can achieve what you describe. However, for bonded atoms that are the focus of the calculation, it's not a good idea to use different levels of theory because the boundary atoms won't behave correctly. Orca does have some computational shortcuts you can use to speed up the calculation, though. If you use a pure DFT functional like BP86 or M06L, you can take advantage of the RI approximation for a massive speedup with little loss in accuracy. You also want to use as small of a basis set as possible, which is probably going to be "sv" or "def2-svp".

See this page for a quick tutorial.

https://sites.google.com/site/orcainputlibrary/dft6770436761403461&sdata=Gbn5iiVuCBoPE4aav8oO1WedRhLm%2BBHeVhc5LQVBcKw%3D&reserved=0>

The Orca manual is much more extensive and explains the keywords better.

-Gerard

________________________________
From:
owner-namd-l_at_ks.uiuc.edu <owner-namd-l_at_ks.uiuc.edu> on behalf of McGuire, Kelly <mcg05004_at_byui.edu>
Sent: Sunday, November 4, 2018 1:59:51 PM
To: namd-l_at_ks.uiuc.edu
Subject: namd-l: QM/MM Question

  A bit of an ignorant question, but does anyone know if it's possible with the NAMD QM/MM to use different theories with different atoms in the QM region? I ask because I have a divalent copper ion, 12 water molecules, and four histidine sidechains in my QM region. MP2 or HF is taking a very long time to complete just one step. ORCA doesn't have PM7. I am using MOPAC's PM7, but I'm not sure when we'll be able to open mopac output files in VMD to view the orbitals. So, with ORCA QM/MM, can I maybe use PM3 for the waters and histidines, and assign copper with MP2 or HF in the same run? Or would that be bad to use different levels of theory in the same region? Thanks!

Kelly L. McGuire

PhD Scholar

Biophysics

Department of Physiology and Developmental Biology

Brigham Young University

LSB 3050

Provo, UT 84602

This archive was generated by hypermail 2.1.6 : Wed Sep 18 2019 - 23:20:27 CDT