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 (Rudberg 2012 J. Phys.: Condens. Matter 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.S Online supplementary data available from stacks.iop.org/JPhysCM/25/152101/mmedia

Predicting solvatochromic shifts and colours of a solvated organic dye: The example of nile red

“…For the TDDFT calculation, an effective system of nile red together with an extended solvation shell is cut out of the much larger classical MD periodic box, where the extent of the solvation shell is defined through a cutoff radius R shell centered on all heavy atoms in nile red and solvent molecules are retained for the TDDFT system if their centre of mass lies within the volume V (R shell ) defined by the combined volume of the spheres on all heavy atoms. To avoid a spurious closure of the band gap [49] and reduce spurious charge transfer states in polar solvents that have been reported for both (semi)-local and hybrid exchange-correlation functionals [9,11], the entire TDDFT system containing the solute and the explicit solvation shell is placed in an implicit solvent model [50] with 0 set to the relative static permittivity of the solvent in question. This guarantees that the explicit solvent molecules on the edge of the explicit solvent volume are correctly screened and the long range continuum effects of the solvent are accounted for.…”

Universal tight binding model for chemical reactions in solution and at surfaces. II. Water

“…6,55 In our TB model, the dipole moment of the water monomer is fitted to 1.86 D, and the average dipole moment in the liquid at 300 K is then found to be 2.40 D. The dipole in the liquid is smaller than predicted using DFT, and it is smaller than was claimed in Ref. 6 due to an error in extracting the induced dipole moments from the trajectories.…”

“…If we define "band gap" as the energy difference between highest unoccupied and lowest occupied states, 55 averaged over an NPT MD simulation, then we can examine this quantity as a function of temperature. The result is shown in Fig.…”

Advanced Potential Energy Surfaces for Molecular Simulation

“…This is of practical significance, because it allows chemical accuracy to be obtained with QM regions smaller than those used in fixed point-charge electrostatic embedding schemes, and is vastly superior to neglecting embedding entirely (i.e. performing DFT on a truncated system), with the latter suffering from additional issues beyond its high computational cost 92 . To illustrate this point, Figure 1 presents a test case of a single diphenylhydramine solute in 330…”