Self-consistent field convergence for proteins: a comparison of full and localized-molecular-orbital schemes

Christian Wick, Matthias Hennemann, James Stewart, Tim Clark

Research output: Contribution to journalArticlepeer-review

74 Downloads (Pure)

Abstract

Proteins in the gas phase present an extreme (and unrealistic) challenge for self-consistent-field iteration schemes because their ionized groups are very strong electron donors or acceptors, depending on their formal charge. This means that gas-phase proteins have a very small band gap but that their frontier orbitals are localized compared to “normal” conjugated semiconductors. The frontier orbitals are thus likely to be separated in space so that they are close to, but not quite, orthogonal during the SCF iterations. We report full SCF calculations using the massively parallel EMPIRE code and linear scaling localized-molecular-orbital (LMO) calculations using Mopac2009. The LMO procedure can lead to artificially over-polarized wavefunctions in gas-phase proteins. The full SCF iteration procedure can be very slow to converge because many cycles are needed to overcome the over-polarization by inductive charge shifts. Example molecules have been constructed to demonstrate this behavior. The two approaches give identical results if solvent effects are included.
Original languageEnglish
Article number2159
Pages (from-to)2159
JournalJournal of Molecular Modeling
Volume20
Issue number2
DOIs
Publication statusPublished - Feb 2014

Keywords

  • Linear scaling,
  • LMO-SCF,
  • NDDO,
  • Proteins,
  • Self-consistent field

Fingerprint

Dive into the research topics of 'Self-consistent field convergence for proteins: a comparison of full and localized-molecular-orbital schemes'. Together they form a unique fingerprint.

Cite this