Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Trying to run QM/MM simulation with calcium binding protein

  1. May 23, 2016 #1
    Hello everyone, I'm a graduate student in a biophysics lab. I'm trying to run a QM/MM simulation with a protein that binds 3 Ca2+ ions, and as such I am having a devil of a time getting the SCF to converge. I've tried DIIS, KDIIS, damping, and a few other methods to try to get convergence with not much luck. The couple of times it has converged, it was after more than 100 iterations, which is just too expensive. Does anyone have any suggestions of things I could try? It doesn't matter much, but I'm using Orca for the QM portion. Thanks a bunch.
     
  2. jcsd
  3. May 23, 2016 #2

    TeethWhitener

    User Avatar
    Science Advisor
    Gold Member

    Plot your energy after each SCF cycle. Are they oscillating or are they going down continually? If they're going down continually, then you probably just need more SCF cycles to reach convergence (I'm assuming the program terminates after a preset number of SCF cycles--which you should be able to adjust). If they're oscillating, you might need to loosen your convergence criteria, at least for preliminary runs. Sometimes if you start with a crappy initial guess for your molecular orbitals, it just takes a while and there's nothing you can do about it. But if you can get it to converge once, you can usually use these new MOs for a better initial guess. I'm sure there are ways to do this in Orca; it's a pretty robust program.

    Of course, all this assumes that your system is reasonably well-behaved. If you have a system which is inherently multi-configuration and you're running an HF-type calculation, your SCF runs will fail to converge because your model isn't a reasonable approximation of the system.
     
  4. May 23, 2016 #3

    TeethWhitener

    User Avatar
    Science Advisor
    Gold Member

    I just realized this made no sense when I read it to myself. What I meant was: 1) get the system to converge with an obscenely low level of theory/small basis set like HF/STO-3G or 3-21G. Then 2) use the MOs from this calculation as an initial guess for a production run with a bigger basis set and/or a higher level of theory.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted



Similar Discussions: Trying to run QM/MM simulation with calcium binding protein
  1. Denaturising Proteins? (Replies: 3)

  2. Protein folding (Replies: 6)

  3. Separate protein (Replies: 1)

  4. Protein fossil (Replies: 4)

Loading...