- Short Report
- Open Access
A comparison of the chemical reactivity of naringenin calculated with the M06 family of density functionals
- Daniel Glossman-Mitnik^{1}Email author
https://doi.org/10.1186/1752-153X-7-155
© Glossman-Mitnik; licensee Chemistry Central Ltd. 2013
- Received: 11 July 2013
- Accepted: 22 August 2013
- Published: 16 September 2013
Abstract
Background
Chemicals generically referred to as flavonoids belong to the group of phenolic compounds and constitute an important group of secondary metabolites due to their applications as well as their biochemical properties. Flavonoids, which share a common benzo- γ-pyrone structure, constitute a kind of compound which are highly ubiquitous in the plant kingdom.
Findings
The M06 family of density functionals has been assessed for the calculation of the molecular structure and properties of the Naringenin flavonoid. The chemical reactivity descriptors have been calculated through Conceptual DFT. The active sites for nucleophilic and electrophilic attacks have been chosen by relating them to the Fukui function indices and the dual descriptor f^{(2)}(r). A comparison between the descriptors calculated through vertical energy values and those arising from the Koopmans’ theorem approximation have been performed in order to check for the validity of the last procedure.
Conclusions
The M06 family of density functionals (M06, M06L, M06-2X and M06-HF) used in the present work leads to the same qualitatively and quantitatively similar description of the chemistry and reactivity of the Naringenin molecule, yielding reasonable results. However, for the case of the M06-2X and M06-HF density functionals, which include a large portion of HF exchange, the calculations considering the validity of the Koopmans’ theorem lead to negative electron affinities.
Keywords
- Naringenin
- DFT
- M06 density functionals
- Conceptual DFT
- Chemical reactivity
Findings
Introduction
Chemicals generically referred to as flavonoids belong to the group of phenolic compounds and constitute an important group of secondary metabolites due to their applications as well as their biochemical properties. Flavonoids, which share a common benzo- γ-pyrone structure, constitute a kind of compound which are highly ubiquitous in the plant kingdom. Over 4, 000 different naturally occurring flavonoids have been discovered, and only in the case of flavones, a specific type of flavonoids, over 36, 000 different chemical structures are possible. Flavonoids are present in a wide variety of edible plant sources, such as fruits, vegetables, nuts, seeds, grains, tea and wine [1].
The knowledge of reactivity on a molecule is an essential concept; it is of a crucial interest because it allows to understand interactions that are operating during a reaction mechanism. In particular electrostatic interactions have been successfully explained by the use of the molecular electrostatic potential [2, 3].
On the other hand, there is no a unique tool to quantify and rationalize covalent interactions, however since 2005 a descriptor of local reactivity whose name is simply dual descriptor [4, 5], has allowed to rationalize reaction mechanisms in terms of overlapping nucleophilic regions with electrophilic regions in order to get a maximum stabilization thus leading to final products or intermediates; all those favorable nucleophilic–electrophilic interactions have been explained as a manifestation of the Principle of Maximum Hardness [6] in addition, chemical reactions have been understood in terms of the The Hard and Soft Acids and Bases Principle [7–10], principle that has been used even with the aim of replacing the use of the Molecular Orbital Theory to understand the whole Chemistry [11]. In fact the present work is a good chance to test the capability of the most recent reactivity descriptors coming from the Conceptual DFT [12–15], therefore the framework of this conceptual theory will be presented in the next section.
Theory and computational details
As mentioned above, DD allows one to obtain simultaneously the preferably sites for nucleophilic attacks (f^{(2)}(r) > 0) and the preferably sites for electrophilic attacks (f^{(2)}(r) < 0) into the system at point r. DD has demonstrated to be a robust tool to predict specific sites of nucleophilic and electrophilic attacks in a much more efficient way than the Fukui function by itself because dual descriptor is able to distinguish those sites of true nucleophilic and electrophilic behavior, in consequence some works have been published with the aim of remarking the powerfulness of f^{(2)}(r) and all those LRDs depending on DD [5, 9, 11, 16–19].
where densities of LUMO and HOMO are represented by ${\rho}_{{}_{\mathrm{L}}}\left(\mathbf{r}\right)$ and ${\rho}_{{}_{\mathrm{H}}}\left(\mathbf{r}\right)$, respectively.
Hence, when an interaction between two species is well described through the use of this LRD, it is said the reaction is controlled by frontier molecular orbitals (or frontier–controlled) under the assumption that remaining molecular orbitals do not participate during the reaction.
When ${f}_{k}^{\left(2\right)}>0$ the process is driven by a nucleophilic attack on atom k and then that atom acts an electrophilic species; conversely, when ${f}_{k}^{\left(2\right)}<0$ the process is driven by an electrophilic attack over atom k and therefore atom k acts as a nucleophilic species.
Settings and computational methods
All computational studies were performed with the Gaussian 09 [22] 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 MIDIY, which is the same basis set as MIDI! with a polarization function added to the hydrogen atoms. The MIDI! basis is a small double-zeta basis with polarization functions on N-F, Si-Cl, Br, and I [23–28].
For the calculation of the molecular structure and properties of the studied system, we have chosen the hybrid meta-GGA density functionals M06, M06L, M06-2X and M06HF [29], which consistently provide satisfactory results for several structural and thermodynamic properties [29–31]. All the calculations were performed in the presence of water as a solvent, by doing IEFPCM computations according to the SMD solvation model [32].
where χ is the electronegativity.
where ε_{ H } and ε_{ L } are the energies of the highest occupied and the lowest unoccupied molecular orbitals, HOMO and LUMO, respectively. However, within the context of density functional theory, the above inequalities are justified in light of the work of Perdew and Levy [33], where they commented on the significance of the highest occupied Kohn–Sham eigenvalue, and proved the ionization potential theorems for the exact Kohn–Sham density functional theory of a many–electron system. In addition the use of the energies of frontier molecular orbitals as an approximation to obtain I and A is supported by the Janak’s Theorem [34]. In particular, The negative of Hartree–Fock and Kohn–Sham HOMO orbital has been found to define upper and lower limits, respectively, for the experimental values of the first ionization potential [35] thus validating the use of energies of Kohn–Sham frontier molecular orbital to calculate reactivity descriptors coming from Conceptual DFT.
that is, the electron accepting power relative to the electron donating power.
Results and discussion
The molecular structure of Naringenin was pre-optimized by starting with the readily available PDB structure, and finding the most stable conformer by means of the Conformers module of Materials Studio through a random sampling with molecular mechanics techniques and a consideration of all the torsional angles. The structure of the resulting conformer was then optimized with the M06, M06L, M06-2X and M06-HF density functionals in conjunction with the MIDIY basis set.
The validity of the Koopmans’ theorem within the DFT approximation is controversial. However, it has been shown [35] that although the KS orbitals may differ in shape and energy from the HF orbitals, the combination of them produces Conceptual DFT reactivity descriptors that correlate quite well with the reactivity descriptors obtained through Hartree-Fock calculations. Thus, it is worth to calculate the electronegativity, global hardness and global electrophilicity for the studied systems using both approximations in order to verify the quality of the procedures.
HOMO and LUMO orbital energies (in eV), ionization potentials I and electron affinities A (in eV), and global electronegativity χ , total hardness η , and global electrophilicity ω of Naringenin calculated with the M06, M06L, M06-2X and M06-HF density functionals and the MIDIY basis set
Property | M06 | M06L | M06-2X | M06-HF |
---|---|---|---|---|
HOMO | -5.9155 | -4.7242 | -7.3281 | -9.2653 |
LUMO | -0.6052 | -1.4060 | 0.1298 | 1.4321 |
χ | 3.2604 | 3.0651 | 3.5992 | 3.9166 |
η | 2.6552 | 1.6591 | 3.7290 | 5.3487 |
ω | 2.0018 | 2.8313 | 1.7370 | 1.4339 |
I | 7.4000 | 6.9619 | 8.1214 | 8.8799 |
A | 0.9393 | 0.7328 | 0.9528 | 0.8997 |
χ | 4.1697 | 3.8474 | 4.5371 | 4.8898 |
η | 3.2304 | 3.1146 | 3.5843 | 3.9901 |
ω | 2.6911 | 2.3763 | 2.8716 | 2.9962 |
The condensed Fukui functions can also be employed to determine the reactivity of each atom in the molecule. The corresponding condensed functions are given by ${f}_{k}^{+}={q}_{k}(N+1)-{q}_{k}\left(N\right)$ (for nucleophilic attack), ${f}_{k}^{-}={q}_{k}\left(N\right)-{q}_{k}(N-1)$ (for electrophilic attack), and ${f}_{k}^{0}=\left[{q}_{k}\right(N+1)-{q}_{k}(N-1\left)\right]/2$ (for radical attack), where q_{ k } is the gross charge of atom k in the molecule.
with c_{ ai } being the LCAO coefficients and S_{ ab } the overlap matrix. The condensed Fukui functions are normalized, thus $\sum _{k}{f}_{k}=1$ and ${f}_{k}^{0}=[{f}_{k}^{+}+{f}_{k}^{-}]/2$.
The condensed Fukui functions have been calculated using the AOMix molecular analysis program [38, 39] starting from single-point energy calculations. We have presented, discussed and successfully applied the described procedure in our previous studies on different molecular systems [40–43].
The condensed dual descriptor has been defined as ${f}^{\left(2\right)}{\left(\mathbf{r}\right)}_{k}={f}_{k}^{+}-{f}_{k}^{-}$ [5, 9]. From the interpretation given to the Fukui function, one can note that the sign of the dual descriptor is very important to characterize the reactivity of a site within a molecule toward a nucleophilic or an electrophilic attack. That is, if f^{(2)}(r)_{ k } > 0, then the site is favored for a nucleophilic attack, whereas if f^{(2)}(r)_{ k } < 0, then the site may be favored for an electrophilic attack [5, 9, 44].
Electrophilic f^{ − } and nucleophilic f^{ + } condensed Fukui functions and f^{ (2) }(r) over the atoms of the Naringenin molecule calculated with the M06, M06L, M06-2X and M06-HF density functionals and the MIDIY basis set
M06 | M06L | M06-2X | M06-HF | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Atom | f^{−} | f^{−} | f^{(2)}(r) | f^{+} | f^{−} | f^{(2)}(r) | f^{+} | f^{−} | f^{(2)}(r) | f^{+} | f^{−} | f^{(2)}(r) |
1 O | 2.13 | 0.99 | 1.14 | 2.32 | 0.35 | 1.97 | 1.58 | 1.23 | 0.35 | 1.05 | 0.50 | 0.55 |
2 C | 11.39 | 0.30 | 11.09 | 10.86 | 0.29 | 10.57 | 11.68 | 0.28 | 11.48 | 11.58 | 0.06 | 11.52 |
3 C | 1.36 | 1.26 | 0.10 | 1.27 | 1.21 | 0.06 | 1.45 | 0.99 | 0.46 | 1.62 | 0.88 | 0.74 |
4 C | 0.61 | 0.73 | -0.12 | 0.51 | 0.32 | 0.19 | 0.82 | 1.14 | -0.32 | 1.38 | 0.23 | 1.15 |
5 C | 6.30 | 6.30 | 0.00 | 5.23 | 8.06 | -2.83 | 8.31 | 0.91 | 7.40 | 11.67 | 0.27 | 11.40 |
6 C | 0.36 | 9.08 | -8.72 | 0.30 | 10.14 | -9.84 | 0.34 | 3.45 | -3.11 | 0.33 | 2.46 | -2.13 |
7 C | 0.49 | 8.61 | -8.12 | 0.41 | 2.57 | -2.16 | 0.51 | 2.67 | -2.16 | 0.63 | 2.75 | -2.12 |
9 C | 17.38 | -0.01 | 17.39 | 16.22 | -0.01 | 16.23 | 18.71 | 0.03 | 18.68 | 16.60 | 0.03 | 16.57 |
10 C | 21.29 | 1.27 | 20.02 | 22.36 | 1.07 | 21.29 | 19.76 | 0.87 | 18.89 | 20.73 | 0.78 | 19.95 |
11 C | 9.32 | 0.21 | 9.11 | 9.57 | 0.15 | 9.42 | 9.10 | 0.31 | 8.79 | 9.04 | 0.09 | 8.95 |
14 C | 0.05 | 2.43 | -2.38 | 0.07 | 1.05 | -1.08 | 0.06 | 5.41 | -5.35 | 0.07 | 6.81 | -6.74 |
15 C | 0.59 | 1.70 | -1.11 | 0.42 | 0.81 | -0.39 | 0.55 | 3.29 | -2.74 | 0.64 | 4.18 | -3.54 |
16 O | 4.11 | 0.02 | 4.09 | 4.44 | 0.03 | 4.41 | 3.45 | 0.02 | 3.43 | 2.68 | 0.01 | 2.67 |
17 C | 1.02 | 1.03 | -0.01 | 0.69 | 0.54 | 0.15 | 1.53 | 0.84 | 0.69 | 2.24 | 0.11 | 2.13 |
18 O | 18.58 | 44.09 | -25.51 | 20.18 | 65.26 | -45.08 | 17.95 | 45.56 | -27.61 | 16.24 | 41.45 | -25.21 |
19 O | 2.31 | 0.98 | 1.33 | 2.65 | 1.10 | 1.55 | 1.81 | 0.40 | 1.41 | 1.28 | 0.06 | 1.22 |
20 C | 0.20 | 3.63 | -3.43 | 0.10 | 0.90 | -0.80 | 0.19 | 9.54 | -9.35 | 0.26 | 10.09 | -9.83 |
23 C | 0.05 | 3.72 | -3.67 | 0.02 | 0.91 | -0.89 | 0.03 | 10.56 | -10.53 | 0.03 | 10.99 | -10.96 |
28 C | 0.32 | 5.50 | -5.18 | 0.26 | 1.73 | -1.47 | 0.34 | 15.69 | -15.35 | 0.42 | 20.75 | -20.33 |
30 O | 0.08 | 7.15 | -7.07 | 0.08 | 2.19 | -2.11 | 0.06 | 15.62 | -15.56 | 0.05 | 12.24 | -12.19 |
It can be concluded from the analysis of the results on Table 2 that the M06, M06L, M06-2X and M06-HF density functionals predict that C10 will be the preferred site for nucleophilic attack. The four density functionals considered in this study display a large negative value of the condensed dual descriptor f^{(2)}(r) over O18, implying that this will be the preferred site for the electrophilic attack.
Electrodonating ( ω ^{ − } ) and electroaccepting ( ω ^{ + } ) powers and net electrophilicity Δ ω ^{ ± } of Naringenin calculated with the M06, M06L, M06-2X and M06-HF density functionals and the MIDIY basis set
Property | M06 | M06L | M06-2X | M06-HF |
---|---|---|---|---|
ω ^{−} | 3.9637 | 4.5712 | 4.0026 | 4.0609 |
ω ^{+} | 0.7035 | 1.5061 | 0.4035 | 0.1443 |
Δω^{±} | 4.6672 | 6.0773 | 4.4061 | 4.2052 |
ω ^{−} | 6.2176 | 4.6892 | 5.5882 | 5.9398 |
ω ^{+} | 1.2124 | 0.8419 | 1.0511 | 1.0500 |
Δω^{±} | 7.4300 | 5.5311 | 6.6393 | 6.9898 |
The results from Table 3 clearly indicate that Naringenin is an electrodonating molecule, with the same result predicted by all the four density functionals considered in this study. However, although the tendency is the same, the results for these descriptors are in poor agreement between those calculated assuming the validity of the Koopmans’ theorem, and those coming from energy differences.
Conclusions
From the whole of the results presented in this contribution it has been clearly demonstrated that the sites of interaction of the Naringenin molecule can be predicted by using DFT-based reactivity descriptors such as the hardness, softness, and electrophilicity, as well as Fukui function calculations. These descriptors were used in the characterization and successfully description of the preferred reactive sites and provide a firm explanation for the reactivity of the Naringenin molecule.
The M06 family of density functionals (M06, M06L, M06-2X and M06-HF) used in the present work leads to the same qualitatively and quantitatively similar description of the chemistry and reactivity of the Naringenin molecule, yielding reasonable results. However, for the case of the M06-2X and M06-HF density functionals, which include a large portion of HF exchange, the calculations considering the validity of the Koopmans’ theorem lead to negative electron affinities.
The calculated descriptors are in agreement with the known experimental facts about the chemical reactivity of the Naringenin molecule presented in the literature (with the exceptions mentioned on the paragraph above). Thus, this make us confidents that similar studies can be pursued with the same degree of accuracy on another flavonoids with analogue structures.
Declarations
Acknowledgements
This work has been partially supported by CIMAV, SC and Consejo Nacional de Ciencia y Tecnología (CONACYT, Mexico). DGM is a researcher of CONACYT and CIMAV.
Authors’ Affiliations
References
- Marín F, Frutos M, Pérez-Alvarez J, Martinez-Sánchez F, Río JD: Flavonoids as nutraceuticals: structural related antioxidant properties and their role on ascorbic acid preservation. Bioactive Natural Products, Volume 26, Part G of Studies in Natural Products Chemistry. Edited by: Rahman A. 2002, Oxford: Elsevier, 741-778.Google Scholar
- Politzer P, Murray J: The fundamental nature and role of the electrostatic potential in atoms and molecules. Theor Chem Acc. 2002, 108: 134-10.1007/s00214-002-0363-9.View ArticleGoogle Scholar
- Murray J, Politzer P: The electrostatic potential: an overview. WIREs Comput Mol Sci. 2011, 1: 153-163. 10.1002/wcms.19.View ArticleGoogle Scholar
- Morell C, Grand A, Toro-Labbé A: New dual descriptor for chemical reactivity. J Phys Chem A. 2005, 109: 205-212. 10.1021/jp046577a.View ArticleGoogle Scholar
- Morell C, Grand A, Toro-Labbé A: Theoretical support for using the Δf(r) descriptor. Chem Phys Lett. 2006, 425 (4–6): 342-346.View ArticleGoogle Scholar
- Pearson R: The principle of maximum hardness. Acc Chem Res. 1993, 26: 250-255. 10.1021/ar00029a004.View ArticleGoogle Scholar
- Pearson R: Hard and soft acids and bases. J Am Chem Soc. 1963, 85: 3533-3539. 10.1021/ja00905a001.View ArticleGoogle Scholar
- Pearson R: Recent advances in the concept of hard and soft acids and bases. J Chem Educ. 1987, 64: 561-10.1021/ed064p561.View ArticleGoogle Scholar
- Gázquez J: The hard and soft acids and bases principle. J Phys Chem A. 1997, 101: 4657-4659. 10.1021/jp970643+.View ArticleGoogle Scholar
- Ayers P, Parr R, Pearson R: Elucidating the hard/soft acid/base principle: a perspective based on half-reactions. J Chem Phys. 2006, 124: 194107-10.1063/1.2196882.View ArticleGoogle Scholar
- Cárdenas C, Rabi N, Ayers P, Morell C, Jaramillo P, Fuentealba P: Chemical reactivity descriptors for ambiphilic reagents: dual descriptor, local hypersoftness, and electrostatic potential. J Phys Chem A. 2009, 113: 8660-10.1021/jp902792n.View ArticleGoogle Scholar
- Chermette H: Density functional theory: a powerful tool for theoretical studies in coordination chemistry. Coord Chem Rev. 1998, 178–180: 699-701.View ArticleGoogle Scholar
- Chermette H: Chemical reactivity indexes in density functional theory. J Comput Chem. 1999, 20: 129-154. 10.1002/(SICI)1096-987X(19990115)20:1<129::AID-JCC13>3.0.CO;2-A.View ArticleGoogle Scholar
- Geerlings P, De Proft F, Langenaeker W: Conceptual density functional theory. Chem Rev. 2003, 103: 1793-1873. 10.1021/cr990029p.View ArticleGoogle Scholar
- Zevatskii Y, Samoilov D: Some modern methods for estimation of reactivity of organic compounds. Russ J Organic Chem. 2007, 43: 483-500. 10.1134/S107042800704001X.View ArticleGoogle Scholar
- Theoretical Aspects of Chemical Reactivity, Volume 19. Edited by: Toro-Labbé A. 2007, Amsterdam: Elsevier ScienceGoogle Scholar
- Ayers P, Morell C, De Proft F, Geerlings P: Understanding the Woodward-Hoffmann rules by using changes in electron density. Chem - Eur J. 2007, 13 (29): 8240-8247. 10.1002/chem.200700365.View ArticleGoogle Scholar
- Morell C, Ayers P, Grand A, Gutiérrez-Oliva S, Toro-Labbé A: Rationalization of Diels-Alder reactions through the use of the dual reactivity descriptor Δf(r). Phys Chem - Chem Phys. 2008, 10: 7239-10.1039/b810343g.View ArticleGoogle Scholar
- Morell C, Hocquet A, Grand A, Jamart-Grégoire B: A conceptual DFT study of hydrazino peptides: assessment of the nucleophilicity of the nitrogen atoms by means of the dual descriptor Δf(r). J Mol Struct: THEOCHEM. 2008, 849: 46-51. 10.1016/j.theochem.2007.10.014.View ArticleGoogle Scholar
- Fuentealba P, Parr RG: Higher-order derivatives in density-functional theory, specially the hardness derivative ∂η/∂N. J Chem Phys. 1991, 94 (8): 5559-5564. 10.1063/1.460491.View ArticleGoogle Scholar
- Parr RG, Yang W: Density Functional Theory of Atoms and Molecules. 1989, New York: Oxford University PressGoogle Scholar
- Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, Barone V, Mennucci B, Petersson GA, Nakatsuji H, Caricato M, Li X, Hratchian HP, Izmaylov AF, Bloino J, Zheng G, Sonnenberg JL, Hada M, Ehara M, Toyota K, Fukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Vreven T, Montgomery JA, et al: Gaussian 09 Revision A.1. 2009, Wallingford: Gaussian Inc.Google Scholar
- Huzinaga S, Andzelm J, Klobulowski M, Radzio-Audselm E, Sakai Y, Tatewaki H: Gaussian Basis Sets for Molecular Calculations. 1984, Amsterdam: ElsevierGoogle Scholar
- Easton R, Giesen D, Welch A, Cramer C, Truhlar D: The MIDI! basis set for quantum mechanical calculations of molecular geometries and partial charges. Theor Chem Acc. 1996, 93: 281-301. 10.1007/BF01127507.View ArticleGoogle Scholar
- Lewars E: Computational Chemistry - Introduction to the Theory and Applications of Molecular and Quantum Mechanics. 2003, Dordrecht: Kluwer Academic PublishersGoogle Scholar
- Young DC: Computational Chemistry - A Practical Guide for Applying Techniques to Real-World Problems. 2001, New York: John Wiley & SonsGoogle Scholar
- Jensen F: Introduction to Computational Chemistry, 2nd edition. 2007, Chichester: John Wiley & SonsGoogle Scholar
- Cramer CJ: Essentials of Computational Chemistry - Theories and Models, 2nd edition. 2004, Chichester: John Wiley & SonsGoogle Scholar
- Zhao Y, Truhlar DG: Density functionals with broad applicability in chemistry. Acc Chem Res. 2008, 41 (2): 157-167. 10.1021/ar700111a.View ArticleGoogle Scholar
- Zhao Y, Truhlar D: The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited States, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor Chem Acc. 2008, 120: 215-241. 10.1007/s00214-007-0310-x.View ArticleGoogle Scholar
- Zhao Y, Truhlar D: Applications and validations of the minnesota density functionals. Chem Phys Lett. 2011, 502: 1-13. 10.1016/j.cplett.2010.11.060.View ArticleGoogle Scholar
- Marenich A, Cramer C, Truhlar D: Universal solvation model based on solute electron density and a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J Phys Chem. 2009, 113: 6378-6396.View ArticleGoogle Scholar
- Perdew J, Burke K, Ersernhof M: Errata: generalized gradient approximation made simple. Phys Rev Lett. 1997, 78: 1396-View ArticleGoogle Scholar
- Janak J: Proof that ∂E/∂ni=ε in density functional theory. Phys Rev B. 1978, 18: 7165-7168. 10.1103/PhysRevB.18.7165.View ArticleGoogle Scholar
- Zevallos J, Toro-Labbé A: A theoretical analysis of the Kohn-Sham and Hartree-Fock orbitals and their use in the determination of electronic properties. J Chilean Chem Soc. 2003, 48: 39-47.View ArticleGoogle Scholar
- Gázquez JL, Cedillo A, Vela A: Electrodonating and electroaccepting powers. J Phys Chem A. 2007, 111 (10): 1966-1970. 10.1021/jp065459f.View ArticleGoogle Scholar
- Chattaraj PK, Chakraborty A, Giri S: Net electrophilicity. J Phys Chem A. 2009, 113 (37): 10068-10074. 10.1021/jp904674x.View ArticleGoogle Scholar
- Gorelsky S: AOMix program for molecular orbital analysis - version 6.5. 2011, [University of Ottawa, Ottawa, Canada]. http://www.sg-chem.net/,Google Scholar
- Gorelsky S, Lever A: Electronic structure and spectra of ruthenium diimine complexes by density functional theory and indo/s. comparison of the two methods. J Organometallic Chem. 2001, 635 (1–2): 187-196.View ArticleGoogle Scholar
- Ruiz-Anchondo T, Glossman-Mitnik D: Computational characterization of the β,β-carotene molecule. J Mol Struct: THEOCHEM. 2009, 913 (1–3): 215-220.View ArticleGoogle Scholar
- Glossman-Mitnik D: Computational Study of 3,4-Diphenyl-1,2,5-Thiadiazole 1-Oxide for organic photovoltaics. Int J Photoenergy. 2009, 2009: 1-7.View ArticleGoogle Scholar
- Glossman-Mitnik D: Computational molecular characterization of coumarin-102. J Mol Struct: THEOCHEM. 2009, 911 (1–3): 105-108.View ArticleGoogle Scholar
- Ruiz-Anchondo T, Flores-Holguín N, Glossman-Mitnik:: Natural carotenoids as nanomaterial precursors for molecular photovoltaics: a computational DFT study. Molecules. 2010, 15 (7): 4490-4510. 10.3390/molecules15074490.View ArticleGoogle Scholar
- Gázquez JL: Chemical reactivity concepts in density functional theory. Chemical Reactivity Theory: A Density Functional View. Edited by: Chattaraj PK. 2009, Boca Raton: CRC Press - Taylor & Francis Group, 7-21.Google Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.