Electrostatic considerations affecting the calculated HOMO-LUMO gap in protein molecules

Abstract

A detailed study of energy differences between the highest occupied and lowest unoccupied molecular orbitals (HOMO-LUMO gaps) in protein systems and water clusters is presented. Recent work questioning the applicability of Kohn-Sham density-functional theory to proteins and large water clusters (E. Rudberg, J. Phys.: Condens. Mat. 2012, 24, 072202) has demonstrated vanishing HOMO-LUMO gaps for these systems, which is generally attributed to the treatment of exchange in the functional used. The present work shows that the vanishing gap is, in fact, an electrostatic artefact of the method used to prepare the system. Practical solutions for ensuring the gap is maintained when the system size is increased are demonstrated. This work has important implications for the use of large-scale density-functional theory in biomolecular systems, particularly in the simulation of photoemission, optical absorption and electronic transport, all of which depend critically on differences between energies of molecular orbitals.

13 pages, 4 figures

References in corpus (8)

• Nearsightedness of Electronic Matter
• O(N) methods in electronic structure calculations
• Minimal parameter implicit solvent model for ab initio electronic structure calculations
• Calculating optical absorption spectra for large systems using linear-scaling density-functional theory
• Difficulties in applying pure Kohn-Sham density functional theory electronic structure methods to protein molecules
• Importance of many body effects in the kernel of hemoglobin for ligand binding
• Ligand Discrimination in Myoglobin from Linear-Scaling DFT+U
• Fermi level pinning can determine polarity in semiconductor nanorods

Cited by in corpus (8)

• How large should the QM region be in QM/MM calculations? The case of catechol O-methyltransferase
• Challenges in Large Scale Quantum Mechanical Calculations
• Combining the Ensemble and Franck-Condon Approaches for Spectral Shapes of Molecules in Solution
• Large scale quantum chemistry with Tensor Processing Units
• Polarization Saturation in Multi-layered Interfacial Ferroelectrics
• Supercell convergence of charge-transfer energies in pentacene molecular crystals from constrained DFT
• Revisiting Artifacts of Kohn-Sham Density Functionals for Biosimulation
• Linear scalability of density functional theory calculations without imposing electron localization