Received date: October 10, 2016; Accepted date: October 26, 2016; Published date: November 15, 2016
Citation: Frau J, Munoz F, Mitnik DG. Validation of the Koopmans Theorem in DFT by Means of the Calculation of the Conceptual DFT Descriptors of Three Fluorescent DNA Staining Dyes. Chem Inform. 2016, 2:2. doi:10.21767/2470-6937.100021
The fulfillment of the "Koopmans' theorem in DFT" is verified by means of the calculation of several global descriptors arising from Conceptual DFT. They have been calculated through a ΔSCF procedure, and by means of the HOMO and LUMO frontier orbitals. The latest Minnesota family of density functionals has been considered and three fluorescent DNA staining dyes have been studied: Hoechst 33258, Hoechst 33342 and Hoechst 34580. On the basis of the obtained Conceptual DFT indices, a series of descriptors have been devised in order to determine the accuracy of each model chemistry considered here in the verification of the mentioned theorem. It is shown that the only density functionals that fulfill this task are those denoted as range-separated hybrids (RSH), while the local density functionals are not useful at all.
Computational chemistry; Molecular modeling; DNA staining dyes; Conceptual DFT; Chemical reactivity theory
Conceptual Density Functional Theory (DFT) or Chemical Reactivity Theory (as it is also known) is a powerful tool for the prediction, analysis and interpretation of the outcome of chemical reactions [1-4].
Following the pioneering work of Parr and others , a useful number of concepts have been derived from the analysis of the density of any molecular system through DFT. These concepts that allow a researcher to make qualitative predictions about the chemical reactivity of a given system can also be quantified and are collectively known as Conceptual DFT Descriptors.
In order to obtain quantitative values of the Conceptual DFT Descriptors, it is necessary to resort to the Kohn-Sham theory trough calculations of the molecular density, the energy of the system, and the orbital energies in particular, those related to the frontier orbitals, that is, HOMO and LUMO [5-10].
The usual way to proceed implies as a first step the choice of model chemistry for the study of the molecular system or chemical reaction of interest. Model chemistry is a combination of a density functional, a basis set, and an implicit solvent model that one consider that can be adequate for the problem under study. There is a plethora of information in the literature about how to choose this model chemistry and one generally follows the experience of previous researchers and his/her own work.
Although the foundations of DFT have established that a universal density functional must exist, and that all of the properties of the system can be obtained through calculations with this functional, in practice one needs to resort to some of the approximate density functional that have been developed during the last thirty years. Due to the fact that these are approximate functional (that is, not a universal functionals), many of them are good for predicting some properties and others are good for another properties. Sometimes, you can find density functionals that are excellent for describing the properties of a given molecular system with a particular functional group, but it is necessary to resort to other density functionals for a different functional group that you want to include in the molecular system under study.
When one is dealing with the study of the chemical reactivity, that is, a process that involve the transference of electrons, it is usual to perform calculations not only of the ground state, but also for open systems like the radical cation and radical anion. These systems are often difficult to converge giving trustworthy results, specially if diffuse functions must be included in the basis set [5-10]. For this reason, it is convenient to have a method that can give all information that one needs directly from the results of the calculation of the ground state of the molecular system under study. In particular, one may want to obtain the ionization potential (I) and electron affinity (A) of the system avoiding the calculation of the radical’s anion and cation. Indeed, the link for this s given by the so-called Koopmans' theorem [7-10], that states that within Hartree-Fock (HF) theory, the I can be approximated by minus the energy of the HOMO, that is, I=-εH. By extension, it is considered that the A can be approximated by minus the energy of the LUMO, that is, A=-εL.
However, the validity of the Koopmans' theorem in DFT is controversial and the problem has been identified with the difference between the fundamental band gap and the HOMO-LUMO gap that is called the derivative discontinuity. Notwithstanding, it has been mentioned recently  that an exact physical meaning can be assigned to the Kohn-Sham (KS) HOMO using "the KS analog of Koopmans' theorem in Hartree- Fock theory", which states that for the exact theory, the KS HOMO is equal to and opposite of the ionization potential, εH=-I [12-15]. Due to the mentioned problem of the discontinuity, a similar Koopmans' theorem that relates the LUMO energy to the electron affinity does not exist. Thus, it has been proposed to circumvent the problem, to consider that the I of the N+1 electron system (the anion) is the same that the A of N electron system . By considering range-separated hybrids (RSH) functional [16-18], where the repulsive Coulomb potential is split into a long-range (LR) and short-range (SR) term, e.g., via r-1=r-1erf(ϒr)+r-1 erfc(ϒr), with ϒ the range-separation parameter, Kronik et al.  showed that with a judicious choice of this last parameter, the validity of the Koopmans' theorem could be enforced.
This ϒ tuning technique can be used for the improvement of the description of the properties predicted by these density functionals. This is a consequence of the better fulfillment of the Koopmans' theorem that leads to better agreement of the orbital energies with the I and A. For example, Lima et al.  have recently presented an improved description of the optical properties of carotenoids by tuning some long-range corrected functionals.
This means that the goodness of a given density functional can be estimated can be estimated by checking how well it follows the "Koopmans' theorem in DFT" that makes it behave closer to the exact density functional, and this will be crucial for a good calculation of the Conceptual DFT descriptors that predict and explain the chemical reactivity of molecular systems. However, the ϒ tuning procedure for the RSH density functionals is system dependent and that implies that different density functionals are going to be used for the calculation of the descriptors for the different molecular systems. Thus, it will be interesting to study other RSH density functional where the parameter is fixed by constructions, although other parameters have been fitted to reproduce some molecular properties. In particular, we are going to consider several density functionals that have shown great accuracy across a broad spectrum of databases in chemistry and physics .
The aim of this work is to conduct a comparative study of the performance of the latest Minnesota family of density functionals for the calculation of the Conceptual DFT descriptors of three fluorescent DNA staining dyes: Hoechst 33258 (HD33258) , Hoechst 33342 (HD33342)  and Hoechst 34580 (HD34580)  whose molecular structures are shown in Figures 1-3. These dyes have important applications in detecting nucleic acids, cancer cells, human papilloma virus (HPV), nucleic acid quantification, detecting abnormal gametes, apoptosis, polynucleotides, proteins and peptides, Alzheimer's disease, treating of cancer and detecting malaria infected red blood cells .
Where χ is the electronegativity.
The global η hardness can be seen as the resistance to charge transfer:
Where εH and εL are the energies of the highest occupied and the lowest unoccupied molecular orbitals, HOMO and LUMO, respectively.
The electrophilicity index ω has been defined as:
The electrodonating (ω-) and electroaccepting (ω+) powers have been definedp as :
It follows that a larger ω+ value corresponds to a better capability of accepting charge, whereas a smaller value of ω- value of a system makes it a better electron donor. In order to compare ω+ with ω-, the following definition of net electrophilicity has been proposed :
That is, the electroaccepting power relative to the electrodonating power.
All computational studies were performed with the Gaussian 09  series of programs with density functional methods as implemented in the computational package. The equilibrium geometries of the molecules were determined by means of the gradient technique. The force constants and vibrational frequencies were determined by computing analytical frequencies on the stationary points obtained after the optimization to check if there were true minima. The basis set used in this work was Def2SVP for geometry optimization and frequencies while Def2TZVP was considered for the calculation of the electronic properties [26,27].
For the calculation of the molecular structure and properties of the studied systems, we have chosen several density functionals from the Minnesota density functionals family, which consistently provide satisfactory results for several structural and thermodynamic properties : M11, which is a is a rangeseparated hybrid meta- GGA , M11L, which is a dual-range local meta-GGA , MN12L, which is a nonseparable local meta-NGA , MN12SX, which is a range-separated hybrid nonseparable meta-NGA , N12, which is a nonseparable gradient approximation , N12SX, which is a range-separated hybrid nonseparable gradient approximation , SOGGA11, which is a GGA density functional  and SOGGA11X, which is a hybrid GGA density functional . In these functionals, GGA stands for generalized gradient approximation (in which the density functional depends on the up and down spin densities and their reduced gradient) and NGA stands for nonseparable gradient approximation (in which the density functional depends on the up/down spin densities and their reduced gradient, and also adopts a nonseparable form). All the calculations were performed in the presence of water as a solvent, by doing IEFPCM computations according to the SMD solvation model .
The molecular structures of HD33258, HD33342 and HD34580 were pre-optimized by starting with the readily available MOL structures, and finding the most stable conformers by means of the Avogadro 1.2.0 program [36,37] through a random sampling with molecular mechanics techniques and a consideration of all the torsional angles. The structures of the resulting conformers were then reoptimized with the M11, M11L, MN12L, MN2SX, N12, N12SX, SOGGA11 and SOGGA11X density functionals in conjunction with the Def2SVP basis set and the SMD solvation model, using water as a solvent.
The HOMO and LUMO orbital energies (in eV), ionization potentials I and electron a nities A (in eV), and global electronegativity χ, total hardness η, global electrophilicity ω, electrodonating power, (ω-), electroaccepting power (ω+), and net electrophilicity Δω± of the HD33258, HD33342 and HD34580 molecules calculated with the M11, M11L, MN12L, MN12SX, N12, N12SX, SOGGA11, and SOGGA11X density functionals and the Def2TZVP basis set using water as solvent simulated with the SMD parameterization of the IEF-PCM model are presented in Tables 1-3, respectively. The upper part of the tables shows the results derived assuming the validity of Koopmans' theorem in DFT (hence the subscript K) and the lower part shows the results derived from the calculated vertical I and A.
Table 1. HOMO and LUMO orbital energies (in eV), ionization potentials I and electron affinities A (in eV), and global electronegativity χ , total hardness η global elec-trophilicity ω, electroaccepting power (ω+), and net electrophilicity Δ ω± of HD33258 calculated with the M11, M11L, MN12L, MN12SX, N12, N12SX, SOGGA11 and SOGGA11X density functionals and the Def2TZVP basis set using water as as solvent simulated with the SMD parametrization of the IEF-PCM model. The upper part of the table shows the results derived assuming the validity of Koopmans' theorem and the lower part shows the results derived from the calculated vertical I and A.
Table 2. HOMO and LUMO orbital energies (in eV), ionization potentials I and electron affinities A (in eV), and global electronegativity χ, total hardness η global elec-trophilicity ω, electroaccepting power (ω+), and net electrophilicity Δ ω± of HD33258 calculated with the M11, M11L, MN12L, MN12SX, N12, N12SX, SOGGA11 and SOGGA11X density functionals and the Def2TZVP basis set using water as as solvent simulated with the SMD parametrization of the IEF-PCM model. The upper part of the table shows the results derived assuming the validity of Koopmans' theorem and the lower part shows the results derived from the calculated vertical I and A.
Table 3. HOMO and LUMO orbital energies (in eV), ionization potentials I and electron affinities A (in eV), and global electronegativity χ, total hardness η global elec-trophilicity ω, electroaccepting power (ω+), and net electrophilicity Δ ω± of HD33258 calculated with the M11, M11L, MN12L, MN12SX, N12, N12SX, SOGGA11 and SOGGA11X density functionals and the Def2TZVP basis set using water as as solvent simulated with the SMD parametrization of the IEF-PCM model. The upper part of the table shows the results derived assuming the validity of Koopmans' theorem and the lower part shows the results derived from the calculated vertical I and A.
Inspired from previous works on this subject [11,19], and with the object of analyzing our results in order to verify the fulfillment of the "Koopmans' theorem in DFT", we have designed several descriptors that relate the results obtained through the HOMO and LUMO calculations with those obtained by means of the vertical I and A with a ΔSCF procedure. However, it must be stressed that it is not our intention to perform a gap-fitting by minimizing a descriptor by choosing optimal range-separation parameter γ, but to check if the density functionals considered in this study, in which, some of the contain a fixed range-separation parameter γ, obbey the "Koopmans' theorem in DFT". As a matter fact, there is no range-separation parameter γ in our designed descriptors. Moreover, we have considered A as minus the energy of the LUMO of the neutral system instead of considering A as minus the energy of the HOMO of the N+1 electron system, as it was in the mentioned works [11,19].
The first three descriptors are related to the simplest fulfillment of the Koop-mans' theorem by relating εH with -I, εL with -A, and the behavior of them in the description of the band gap:
Next, we consider four other descriptors that analyze how well the studied density functionals are useful for the prediction of the electronegativity χ, the global hardness η and the global electrophilicity ω, and for a combination of these Conceptual DFT descriptors, just considering the energies of the HOMO and LUMO or the vertical I and A:
Where D1 stands for the first group of Conceptual DFT descriptors. Finally, we designed other four descriptors to verify the goodness of the studied density functionals for the prediction of the electroaccepting power ω+, the electrodonating power ω-, the net electrophilicity Δω±, and for a combination of these Conceptual DFT descriptors, just considering the energies of the HOMO and LUMO or the vertical I and A:
Where D2 stands for the first group of Conceptual DFT descriptors. The results of the calculations of JI, JA, JGap, Jχ, Jη, Jω, JD1, Jω+,Jω-, JΔω± and JD2 for the HD33258, HD33342 and HD34580 molecules are displayed in Tables 4-6, respectively.
Table 4. Descriptors JI, JA, JGap, Jχ, Jη, Jω, JD1, Jω+,Jω-, JΔω± and JD2 for the HD33258 molecule calculated from the results of Table 1.
As can be seen from Tables 1-3, and the results presented in Tables 4-6, the "Koopman's theorem in DFT" holds with great accuracy for the MN12SX and N12SX density functionals, which are a range-separated hybrid meta-NGA and a range-separated hybrid NGA density functionals, respectively. Indeed, the values of JI, JA and JGap are not exactly zero. However, their values can be favorably compared with the results presented for these quantities in the work of Lima et al. , where the minima has been obtained by choosing a parameter that enforces that behavior.
Table 5. Descriptors JI, JA, JGap, Jχ, Jη, Jω, JD1, Jω+,Jω-, JΔω± and JD2 for the HD33342 molecule calculated from the results of Table 2.
Table 6. Descriptors JI, JA, JGap, Jχ, Jη, Jω, JD1, Jω+,Jω-, JΔω± and JD2 for the HD34580 molecule calculated from the results of Table 3.
It is interesting to see that the same density functionals also fulfill the "Koop-mans' theorem in DFT" for the other descriptors, namely Jχ, Jη, Jω and JD1 as well as for Jω-,Jω+, JΔω±, and JD2. These results are very important, because they show that it is not enough to rely only in JI, JA and JGap. For example, if we consider only Jχ, for all of the density functionals considered, the values are very close to zero. As for the other descriptors, only the MN12SX and N12SX density functionals show this behavior. That means that the results for Jχ are due to a fortituous cancellation of errors.
The usual GGA (SOGGA11) and hybrid-GGA (SOGGA11X) are not good for the fulfillment of the "Koopmans' theorem in DFT", and the same conclusion is valid for the local functionals M11L, MN12L and N12.
An important fact is that although the range-separated hybrid NGA and range-separated hybrid meta-NGA density functionals can be useful for the calculation of the Conceptual DFT descriptors, it is not the same for the range-separated hybrid GGA (M11) density functional. An inspection of Tables 1-3 shows that this is due to the fact that this functional describes inadequately the energy of the LUMO, leading to negative values of A, which are in contradiction with the ΔSCF results.
From the whole of the results presented in this contribution it has been clearly demonstrated that the chemical reactivity of the HD33258, HD33342 and HD34580 molecules can be predicted by using DFT-based reactivity descriptors such as the electronegativity, global hardness, global electrophilicity, electrodonating and electroaccepting powers, and net electrophilicity.
The Minnesota family of density functionals (M11, M11L, MN12L, MN12SX, N12, N12SX, SOGGA11 and SOGGA11X) have been tested for the fulfillment of the "Koopmans' theorem in DFT" by comparison of the HOMO- and LUMO-derived values with those obtained through a ΔSCF procedure. It has been shown that the range-separated hybrid meta-NGA density functional (MN12SX) and the range-separated hybrid NGA density functional (N12SX) are the best for the accomplishment of this objective. As such, they are a good alternative to those density functionals whose behavior have been tuned through a gap-fitting procedure and a good prospect for their usefulness in the description of the chemical reactivity of molecular systems of larger size.
This work has been partially supported by CIMAV, SC and Consejo Nacional de Ciencia y Tecnolog a (CONACYT, Mexico) through Grant 219566/2014 for Basic Science Research and Grant 265217/2016 for a Foreign Sabbatical Leave. Daniel Glossman- Mitnik conducted this work while a Sabbatical Fellow at the University of the Balearic Islands from which support is gratefully acknowl-edged. This work was cofunded by the Ministerio de Economia y Competitividad (MINECO) and the European Fund for Regional Development (FEDER) (CTQ2014-55835-R).
All Published work is licensed under a Creative Commons Attribution 4.0 International License
Copyright © 2018 All rights reserved. iMedPub LTD Last revised : September 23, 2018