Re: QM/MM Question

From: Gerard Rowe (GerardR_at_usca.edu)
Date: Mon Nov 05 2018 - 13:41:11 CST

Next message: McGuire, Kelly: "Re: QM/MM Question"
Previous message: Giacomo Fiorin: "Re: Error with D-amino acids mid-sequence using NAMD 2.12/2.13b2 and CHARMm36"
In reply to: McGuire, Kelly: "Re: QM/MM Question"
Next in thread: McGuire, Kelly: "Re: QM/MM Question"
Reply: McGuire, Kelly: "Re: QM/MM Question"
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

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/dft

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!