CC  Vol.6 No.3 , July 2018
Theoretical Study by Density Functional Theory Method (DFT) of Stability, Tautomerism, Reactivity and Prediction of Acidity of Quinolein-4-One Derivatives
A theoretical study of the reactivity of quinoline-4-one derivatives is undertaken in order to understand the involved mechanisms. The calculations were carried out in gas phase and in N, N-Dimethylformamide (DMF) solution. The Density Functional Theory (DFT) with B3LYP functional associated to 6-311G (d) and 6-311+G (d) bases is used to perform these calculations. The results of the thermodynamic parameters showed that there is an equilibrium relation between the different tautomers. This equilibrium can be used to explain the failure to obtain tetrahydroquinoline from 5,8-dimethoxy-quinolin-4-one. Reactivity analysis from Frontier Molecular Orbitals theory and Fukui function calculations revealed that ketone forms are less reactive than enol ones. The methoxyl substituent decreases the acidity of the nitrogen and oxygen atoms of quinolin-4-one while the bromine increases the acidity of the same sites. These results foresee that nitrogen deprotonation in the case of the brominated compound is easier than in the case of methoxylated ones.

1. Introduction

Life expectancy, as well as pollution of various kinds, associated with new technologies and different risk factors such as alcohol or tobacco, increase in the number of pathologies. Cancer is found among these pathologies [1] . Thus, many researchers focused themselves on new treatments finding. These researches help them to discover, the very interesting role of quinone compounds. They are often found in biological systems such as photosynthetic or enzymatic systems, and in chemistry (oxidation, conductive polymers, etc.) [2] [3] [4] . These molecules are present in many structures with high antitumoral potential [5] [6] [7] [8] [9] . According to the usefulness of the biological interest of this nucleus, several ways of synthesis have been invented in order to obtain it [6] [7] [8] [10] [11] [12] [13] . Quinones have the same basic motif which is cyclohexadiene type. They possess units having exocyclic double bonds, giving them a large number of physicochemical properties. The para-dimethoxyl tetrahydroquinoline ring is an essential element for the synthesis of these quinones. Thus, 5,8-dimethoxy-4,4-dimethyltetrahydroquinoline is presented as a key intermediate which is necessary for the synthesis of these quinones. Although that of tetrahydroquinoline is described, many synthetic ways have been envisaged with the purpose of its elaboration [11] [14] [15] [16] . Moreover, it has been shown that obtaining this key nucleus from 5,8-dimethoxy-quinolin-4-one by using the conditions of Reetz does not give the desired product [17] . These experimental studies have also shown that the reduction of the double bonds’ number of quinolin-4-one is made in order to obtain dihydroquinolin-4-one but it has not been observed [18] . A theoretical study of quinolein-4-one derivatives is therefore undertaken in order to try to understand and predict the stability and reactivity of these compounds. The calculations are carried out in the gaseous phase and in N, N-dimethylformamide (DMF) solution with the functional B3LYP of the DFT method associated to the following bases 6-311G (d) and 6-311+G (d). The experimental difficulties evoked in obtaining 5,8-dimethoxy-4,4-dimethyltetrahydroquinoline from quinoline lead us to envisage the theoretical calculation of enthalpies of formation, enthalpies of reaction and free reaction enthalpies. The equilibrium constants are determined. Currently, the method of functional density theory (DFT) has been accepted as a popular approach for the calculation of the structural characteristics and energies of molecules by the community [19] and for the efficiency and accuracy of the evaluation of molecular properties [20] . This work aims to determine theoretically the preferred sites of electrophilic and nucleophilic attack on the aromatic carbon atoms and heteroatoms contained in the benzimidazolyl nucleus by different quantum chemistry methods. The molecular structures of the tautomers of quinoline derivatives studied are shown in Figure 1 and Figure 2. Figure 1 shows the two studied tautomeric forms; that is, ketone and enol forms. As for Figure 2, it shows the ketone tautomers (Q-4-one, DMQ and DBrQ) and enol (Q-4-hydroxy, DMHQ and DBrHQ).

2. Material and Methods

2.1. Level of Calculation Theory

Theoretical calculations are made by using the software called Gaussian 03 [21] .

Figure 1. Keto-enol equilibrium in quinolin-4-one and atom numbers. The direction 1 is the equilibrium frome keto to enol and direction 2 is enol to keto.

Figure 2. Structures of tautomers of quinoline-4-one studied derivatives.

The different molecular structures have been optimized using the Functional Density Theory (DFT) [22] [23] [24] . The level of calculation used is B3LYP associated with bases 6-311G (d) and 6-311+G (d). The calculations are carried out gaseous phase and solution. N, N-Dimethylformamide (DMF) being one of the solvents used for the synthesis of these quinolines [18] , the calculations in solution are thus carried out in DMF. The Polarized Continuum Model (PCM) [25] [26] is taken as a solvation model. The optimization of the molecular geometries is followed by calculations of vibration frequencies at the same level of the theory in order to verify that the structures obtained are minimal.

2.2. Thermodynamic Parameters of Reaction

Knowledge of variations in energy contributions to internal energy which occurs between the products and the reagents helps to characterize energetically a chemical reaction. For a given energy parameter X, its variation is determined according to the following relation (1):

Δ r X = X ( products ) X ( reactants ) (1)

The thermodynamic determined quantities in this study are the reaction enthalpies ( Δ r H ), the reaction free enthalpies ( Δ r G ) and the equilibrium constants (Keq). These quantities are calculated by using relations (2) and (3) below:

Δ r G = Δ f G ( Products ) Δ f G ( Reactants ) (2)

Δ r H = Δ f H ( Products ) Δ f H ( Reactants ) (3)

The equilibrium constant denoted Keq characterizes the equilibrium state of a chemical reaction. The value of the equilibrium constant depends solely on the considered chemical reaction and the temperature. In this work, the equilibrium constants are determined at 298 K from the following relation:

K e q = exp ( Δ r G ( T ) / R T )

2.3. Reactivity Parameters Derived from Conceptual DFT

The Frontier Molecular Orbital (FMO) theory allows to understand the chemical reactivity. In particular, the prediction of the acidity of the different tautomers is made from the calculation of the deprotonation energies. The study of the reactivity of these molecules is carried out by taking into account the FMO [27] . This theory indicates the ease with which a molecule can be transformed or can react with other molecules. It gives a prediction of the excitation properties of a molecule. As a result, it constitutes quantum parameters for the determination of molecular reactivity [28] [29] . The smaller the energy gap between the HOMO and LUMO orbitals, the easier the molecule can be excited. Similarly a large gap HOMO-LUMO indicates a high stability of the studied molecule which implies its low reactivity during chemical reactions. The condensed Fukui functions are obtained using the Yang and Mortier procedure [30] [31] [32] based on the finite difference method. The different values of local descriptors are calculated from equations (4) and (5) below:

f k + = q k ( N + 1 ) q k ( N ) (4)

f k = q k ( N ) q k ( N 1 ) (5)

where: q k ( N ) : Electronic population of the atom k in the neutral form; q k ( N + 1 ) : Electronic population of the atom k in the anionic form; q k ( N 1 ) : Electronic population of the atom k in the cationic form.

We can gather all the information contained in the two Fukui functions in a single function called dual descriptor. This function is represented by expression (6) which is below.

Δ f = f k + f k (6)

where Δf is a function which represent the reactivity of the molecule, it is positive in the electrophilic areas and negative in the nucleophilic ones [33] .

The chemical potential μ measures the tendency of the electron cloud to escape from the molecule. It is a global property of the molecular system which is equal to the slope of the total electronic energy E. This energy is a function of the number of electrons N with an external potential V (r) constant. The studied parameter corresponds to the opposite of the electronegativity χ as defined by Pauling and Mulliken [34] [35] . This parameter is determined by the relation below (7)

μ = ( E N ) = χ (7)

μ can be expressed as a function of the ionization potential IP and the electronic affinity EA and it is calculated by using the following relation (8).

μ = IP + EA 2 = χ (8)

The hardness ɳ η is determined from the first derivative of the chemical potential. And it must be noted that the softness S is the inverse of the hardness [36] [37] . They are given by the following relation (9)

η = ( μ N ) = ( 2 E E 2 ) = 1 S (9)

According to the theory of acids and bases, developed by Pearson [38] , these quantities can be expressed as a function of both the Ionization Potentials (IP) and Electronic Affinity (EA). It permits to determine the relation (10)

η = 1 S = ( IP EA ) / 2 (10)

In order to characterize other electronic behaviors of the systems, we have interested ourselves to other descriptors such as the electrophile index ω which evaluate the ability of a molecule to generate an electron transfer. It is calculated according to the following formula [39] .

ω = μ 2 2 η = χ 2 2 η (11)

These descriptors can be calculated in a very simple way in the framework of Koopmans’ approximation [40] according to which:


E HOMO and E LUMO are the energies of the HOMO and LUMO orbitals.

2.4. Parameters of Tautomers Acidity Prediction

De protonation energies are widely used for the prediction of Brönsted acidity in gaseous and aqueous phases [41] [42] [43] [44] . The deprotonation reaction is defined by the following relationship (12).

A H A + H + (12)

Δ G 0 298 = Δ H 0 298 T Δ S 0 298 (13)

The contribution of entropy is given by the expression (14):

T Δ S 0 298 = T [ S ( A ) + S ( H + ) S ( A H ) ] (14)

For a temperature of T = 298 K at standard pressure, the second term T S ( H + ) = 32.5 Kj .mol -1 gives the expression (15) [42] [43] .

Δ G 0 298 = Δ H 0 298 T [ S ( A ) S ( A H ) ] 32.5 (15)

which leads to the following relation (16):

Δ G 0 298 = G ( A ) G ( A H ) 32.5 (16)

Previous theoretical work on the prediction of acidity has shown that gradient-corrected functional and hybrid functionals such as B3LYP associated with bases 6-311+G(d) and 6-311+G(3df, 3pd) lead to results in good agreement with the experimental data [45] [46] [47] [48] . This motivates our choice of the level of calculation implemented in this work.

3. Results and Discussion

The ketone-enol tautomeric equilibrium in the case of quinolin-4-one and the tautomeric structures which were the subject of this study are shown respectively in Figure 1 and Figure 2.

The performed calculations highlight the following parameters:

1) Thermodynamic parameters (Table 1);

2) The energies of boundaries molecular orbital and energetic gaps (Table 2);

3) Reactivity indices (Table 3);

4) The functions of Fukui (Table 4);

5) The deprotonation energies (Table 5).

3.1. Thermodynamic Parameters and Relative Stability

The calculated thermodynamic parameters are free enthalpies and enthalpies of formation. From these thermodynamic quantities, the reaction enthalpies and the equilibrium constants are determined. The order of relative stability of the tautomers is established from the enthalpy energies of formation of the different compounds. The values of the thermodynamic quantities calculated in gaseous phase and in DMF solution are gathered in Table 1.

Keq1 and Keq2 are respectively equilibrium constants in directions 1 and 2

Table 1. Enthalpy of formation, reaction enthalpy and equilibrium constants of tautomers of quinoline derivatives calculated at the B3LYP/6-311G (d) theory level.

Table 2. Energies of HOMO, LUMO boundary orbitals and energy gaps of tautomers of quinolin-4-one derivatives calculated at B3LYP/6-311G (d).

(Figure 1). The values of the equilibrium constants mentioned in Table 1 are approximately equal to 1 which indicates that the reaction is in its equilibrium domain. Experimental data have shown that the probable existence of a tautomer relationship could explain the failure of tetrahydroquinoline to be obtained from 5,8-dimethoxy-quinolin-4-one by Reetz conditions [17] . Equilibrium constancy calculations confirm this hypothesis. They highlight a tautomeric equilibrium relationship between different tautomers. For two given tautomers, the results show that the ketone forms (Q-4-one, DMQ, DBrQ) have the lowest formation enthalpy values (Table 1). These are therefore more stable than the enol forms. In the case of 5,8-dimethoxy-quinolin-4-one, the enthalpies of formation of tautomers have very similar values. Both ketone and enol forms have a very similar stability. The reaction enthalpies are all positive which means that the equilibrium reaction is an endothermic one.

3.2. Chemical Reactivity

3.2.1. Theory of Frontier Molecular Orbitals

Table 2 contains the energies of the HOMO LUMO boundary orbitals and the HOMO-LUMO energy gaps.

The energy gap values decrease in the DMF solvent (Table 2) compared with those obtained in the gas phase. The molecules are therefore more reactive in solution. The obtained values of the ketone tautomers are greater than those obtained with enols in the two areas. The ketone forms are less polarizable than tautomers ones. And they have low chemical reactivity with high kinetic stability compared to enol tautomers. The substituents on quinoline cause a decrease in the energy gap. The smallest energy gap is obtained with bromine. This substituent therefore increases the reactivity of quinolin-4-one. At the enol form level, the quinoline-4-one methoxylated has the lowest value of ΔE which means that the methoxyl group leads to the most reactive enol form. We note however that among two tautomers, the different energy gaps are close to each other. This denotes that the tautomers reactivity is very closed to each other. This can explain the equivalence of N- and O-benzylation yielding rates in the benzylation reaction of quinolin-4-one [17] .

3.2.2. Global Indices of Reactivity

The study of the global reactivity of tautomers is based on the computation of global indices deduced from electronic properties. Global indices of the reactivity of studied tautomers are shown in Table 3.

According to these parameters, the chemical reactivity varies with the structure of the tautomers. Global softness (S) values of ketonic tautomers are lower than enol ones, indicating that ketone tautomers are soft. Also, we note that ketones have their electronegativity values (χ) lower than those of enols; except (Q-4-one) in DMF (3.823 eV). They are therefore good donors of electrons. In terms of ketone forms, the brominated derivative has the greatest value of electronegativity. It is therefore the best electron acceptor compared to the methoxylated derivatives. These data show that the bromine atom contributes to increase this index of reactivity. Concerning the methoxyl substituent it causes a decrease in the value of the same parameter. In addition, electrophile index

Table 3. Ionization potential (IP, eV), electron affinity (EA, eV), chemical potential (μ, eV), chemical hardness (η, eV), softness (S, eV-1) and overall electrophile index (ω, eV). of tautomers and some quinolin-4-one derivatives calculated at B3LYP/6-311G (d).

values of ketonic tautomers are lower than those of enol indicating that ketonic tautomers are less electrophilic.

3.2.3. Fukui Functions

The calculation of Fukui’s functions was performed using Hirschfeld populations [49] . The results are summarized in Table 4. This calculation relates to all the atoms of the quinolin-4-one ring (Figure 1).

The calculated Fukui functions indicate positive values of Δf at the C2, C3, C4, C5, C7 and C8 carbon atoms in all molecules (Table 4). These atoms are electrophilic sites. A nucleophilic attack in these molecules will happen on C4 and C5 carbons which have the highest values of Δf. Failure of C4 carbon dimethylation and reduction of the C2 = C3 double bond of the quinoline ring may be due to the fact that these atoms are nucleophilic sites. They are therefore favorable to an electrophilic attack. The N1, C8, C9 and O11 atoms have a negative value of Δf. These atoms are nucleophilic sites. They are susceptible to electrophilic attack. The nitrogen atom N1 has the highest negative value of Δf which means that an electrophilic attack will happen preferentially on this atom. The probable sites of nucleophilic and electrophilic attack are the same in all the molecules. The substitution does not modify the reactivity sites of these molecules.

3.3. Prediction of Acidity of Tautomers of Quinoline-4-one Derivatives

Gibbs free energies of deprotonation ( Δ G 0 298 ) are calculated in the gas phase at B3LYP/6-311+G (d). The values of this thermodynamic quantity are grouped together in Table 5.

The grouped data in Table 5 show an increase in the free enthalpy of the deportation reaction with the methoxyl group. This energy decreases with bromine.

Table 4. Calculated Fukui functions in gas phase at B3LYP/6-311 + G (d) ketone forms.

Table 5. Free enthalpies of formation and deprotonation of quinoline-4-one tautomers calculated in the gas phase at B3LYP/6-311+G (d) level.

According to these data, we can establish the increasing order of the studied molecules by taking into account the free energy of deprotonation ( Δ G 0 298 ) at 298 K: This order is: ΔG[DBrQ (O11-H)] < ΔG[DBrQ (N1-H)] < ΔG[Q-4-one (O11-H)] < ΔG[Q-4-one (N1-H)] < ΔG[DMQ (O11-H)] < ΔG[DMQ (N1-H)]. Bromine contributes to reducing the free energies of deprotonation of the nitrogen and oxygen atoms while the methoxyl group increases the same thermodynamic parameter on the same atoms. These results indicate that a deprotonation of the N1-H nitrogen in the case of brominated quinolin-4-one will be easier than that of the nitrogen (N1-H) of the methoxylated derivative. Work has shown that there is an inverse relationship between the amplitude of ΔG and the strength of the acid. The lower the value of ΔG, the more acidic the site is [41] [42] [43] [44] . In this way, the nitrogen atoms N1-H and oxygen O11-H in the brominated compounds are more acidic. These sites are less acidic in the methoxylated derivatives. Through the analysis of these free energies of deprotonation, we can establish the following decreasing order of studied molecules according to acidity: DBrQ (O11-H) > DBrQ (N1-H) > Q-4-one (O11-H)> Q-4-one (N1-H)] > DMQ (O11-H)> DMQ (N1-H). For a given quinoline, the deprotonation energy of the hydroxyl (O11-H) of the enol form is lower than that of the nitrogen (N1-H) of the ketone form.

4. Conclusion

The results of this theoretical study on the quinoline-4-one derivatives have shown that the ketone forms are therefore more stable than the enol ones. Experimental data have shown that the probable existence of a tautomerism can explain the failure to obtain tetrahydroquinoline from 5,8-dimethoxy-quinolin-4-one. Equilibrium constancy calculations confirmed this hypothesis. They highlight equilibrium relationship between different tautomers which is an endothermic reaction. Reactivity analysis from boundary molecular orbital theory has also shown that bromine (Br) increases the reactivity of quinoline. On the other hand, brominated quinolin-4-one is more reactive than its methoxylated analogues. Calculation of Fukui functions showed that C2 and C3 carbon atoms are electrophilic sites. This fact is a possible explanation for the failure of the reduction of the double bond in quinolin-4-one. These results indicate that deprotonation of N1-H nitrogen in the case of brominated quinoline would be easier than that of nitrogen (N1-H) of the methoxylated derivatives. The results predict a higher acidity of the N1-H and O11-H sites in the case of the brominated compounds. A theoretical study of the reaction mechanism of the N1-H nitrogen deprotonation reaction of these molecules is envisaged in order to bring more precision to this first conclusion.


The authors would like to thank Lamoussa Ouattara for his contribution to the realization of this work.

Cite this paper
Bédé, A. , Assoma, A. , Yapo, K. , Koné, M. , Koné, S. , Koné, M. , N’Guessan, B. and Bamba, E. (2018) Theoretical Study by Density Functional Theory Method (DFT) of Stability, Tautomerism, Reactivity and Prediction of Acidity of Quinolein-4-One Derivatives. Computational Chemistry, 6, 57-70. doi: 10.4236/cc.2018.63005.

[1]   Belot, A., Grosclaude, P., Bossard, N., Jougla, E. and Benhamou, E. (2008) Cancer Incidence and Mortality in France over the period 1980-2005. Revue d’épidémiologie et de santé publique, 56, 159-175.

[2]   Kumar, A., Chelvam, V., Sakkarapalayam, M., Li, G.P., Sanchez-Cruz, N.S., Pinero, P.S., Low, P.S. and Alegria, A.E. (2016) Synthesis and Evaluation of Folate-Conjugated Phenanthraquinones for Tumor-Targeted Oxidative Chemotherapy. Open Journal of Medicinal Chemistry, 6, 1-17.

[3]   Apel, K. and Hirt, H. (2004) Reactive Oxygen Species: Metabolism, Oxidative Stress, and Signal Transduction. Annual Review of Plant Biology, 55, 373-399.

[4]   Fieser, L.F. (1929) The Reduction Potentials of Various Phenanthrenequinones. Journal of the American Chemical Society, 51, 3101-3111.

[5]   Shinkre, B.A., Raisch, K.P., Fan, L. and Velu, S.E. (2007) Analogs of the Marine Alkaloid Makaluvamines: Synthesis, Topoisomerase II Inhibition, and Anticancer Activity. Bioorganic & Medicinal Chemistry Letters, 17, 2890-2893.

[6]   Shinkre, B.A., Raisch, K.P., Fan, L. and Velu, S.E. (2008) Synthesis and Antiproliferative Activity of Benzyl and Phenethyl Analogs of Makaluvamines. Bioorganic & Medicinal Chemistry, 16, 2541-2549.

[7]   Ezell, S.J., Li, H., Xu, H., Zhang, X., Gurpinar, E., Zhang, X., Rayburn, E.R., Sommers, C.I., Yang, X., Velu, S.E., Wang, W. and Zhang, R. (2010) Preclinical Pharmacology of BA-TPQ, a Novel Synthetic Iminoquinone Anticancer Agent. Marine Drugs, 8, 2129-2141.

[8]   Hoang, H., LaBarbera, D., Mohammed, K., Ireland, C. and Skibo, E. (2007) Synthesis and Biological Evaluation of Imidazoquinoxalinones, Imidazole Analogues of Pyrroloiminoquinonemarine Natural Products. Journal of Medicinal Chemistry, 50, 4561-4571.

[9]   Barrows, L.R., Radisky, D.C., Copp, B.R., Swaffar, D.S., Kramer, R.A., Warters, R.L. and Ireland, C.M. (1993) Makaluvamines, Marine Natural Products, Are Active Anti-Cancer Agents and DNA Topo II Inhibitors. Anti-Cancer Drug Design, 8, 333-347.

[10]   Delfourne, E. (2008) Analogues of Marine Pyrroloiminoquinone Alkaloids: Synthesis and Antitumor Properties. Anticancer Agents. Journal of Medicinal Chemistry, 8, 910-916.

[11]   Hajra, S., Maji, B. and Karmakar, A. (2005) Lewis Acid Catalyzed Intramolecular Halo-Arylation of Tethered Alkenes Using N-halosuccinimide (NXS) as the Halogen Source: A General Method for the Synthesis of Chromanones, Chromans, Quinolones, Tetrahydroquinolines and Tetralins. Tetrahedron Letters, 46, 8599-8603.

[12]   Antunes, E.M., Beukes, D.R., Kelly, M., Samaai, T., Barrows, L.R., Marshall, K.M., Sincich, C. and Davies-Coleman, M.T. (2004) Cytotoxic Pyrroloiminoquinones from Four New Species of South African Latrunculid Sponges. Journal of Natural Products, 67, 1268-1276.

[13]   Monks, A., Scuderio, D.A., Skehan, P., Shoemaker, R.H., Paull, K.D., Hose, V.C., Langley, J., Cronise, P., Vaigro-Wolff, A., Grey-Goodrich, M., Campball, H. and Boyd, M.R. (1991) Feasibility of a High-Flux Anticancer Drug Screen Using a Diverse Panel of Cultured Human Tumor Cell Lines. Journal of the National Cancer Institute, 83, 757-766.

[14]   Clerici, A. and Porta, O. (1990) Arylative Amination of Aldehydes Promoted by Aqueous Titanium Trichloride. Tetrahedron Letters, 31, 2069-2072.

[15]   Kitahara, Y., Tamura, F., Nishimura, M. and Kubo, A. (1998) Synthesis of Meridine, Cystodamine, and Related Compounds Including Iminoquinolinequinone Structure. Tetrahedron, 54, 8421-8432.

[16]   Marcos, A., Pedregal, C. and Avendano, C. (1995) Synthesis of 1,5-and 1,8-Diazaanth-raquinones by Reaction of Aminoquinolinequinones with β-dielectrophiles. Tetrahedron, 51, 1763-1774.

[17]   Reetz, M.T. and Seitz, T. (1987) Regio-und Stereoselektive Carbosulfenylierung von Olefinen. Angewandte Chemie, 99, 1081-1082.

[18]   Bouclé, S. (2010) Synthèse d’analogues d’alcaloides marins à potentiel anti-tumoral. Thèse de Doctorat, Université Francois-Rabelais de Tours.

[19]   Kurt, M., Sertbakan, T.R. and Ozduran, M. (2008) Spectrochim, an Experimental and Theoretical Study of Molecular Structure and Vibrational Spectra of 3-and 4-Pyridineboronic Acid Molecules by Density Functional Theory Calculations. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 70, 664-673.

[20]   Ravikumar, C., Joe, I.H. and Jayakumar, V.S. (2008) Charge Transfer Interactions and Nonlinear Optical Properties of Push-Pull Chromophore Benzaldehyde Phenylhydrazone: A Vibrational Approach. Chemical Physics Letters, 460, 552-558.

[21]   Frisch, M.J., Trucks, G.W., Schlegel, H.B., et al. (2003) Gaussian 03, Revision A.1. Gaussian Inc., Pittsburgh.

[22]   Kohn, W. and Sham, L. (1965) Self-Consistent Equations Including Exchange and Correlation Effects. Physical Review, 140, 1133-1138.

[23]   Becke, A.D. (1997) Density-Functional Thermochemistry. V. Systematic Optimization of Exchange-Correlation Functionals. The Journal of Chemical Physics, 107, 8554-8560.

[24]   Bauernschmitt, R. and Ahlrichs, R. (1996) Stability Analysis for Solutions of the Closed Shell Kohn-Sham Equation. The Journal of Chemical Physics, 104, 9047-9052.

[25]   Canses, E., Mennucci, B. and Tomasi, J. (1997) A New Integral Equation Formalism for the Polarizable Continuum Model: Theoretical Background and Applications to Isotropic and Anisotropic Dielectrics. The Journal of Chemical Physics, 107, 3032-3041.

[26]   Tomasi, J. and Persico, M. (1994) Molecular Interactions in Solution: An Overview of Methods Based on Continuous Distributions of the Solvent. Chemical Reviews, 94, 2027-2094.

[27]   Fujimoto, H. and Fukui, K. (1972) Molecular Orbital Theory of Chemical Reactions. Advances in Quantum Chemistry, 6, 177-201.

[28]   Belletête, M., Morin, J.F., Leclerc, M. and Durocher, G. (2005) A Theoretical, Spectroscopic, and Photophysical Study of 2,7-Carbazolenevinylene-Based Conjugated Derivatives. The Journal of Physical Chemistry A, 109, 6953-6959.

[29]   Jun-ichi, A. (1999) Reduced HOMO-LUMO Gap as an Index of Kinetic Stability for Polycyclic Aromatic Hydrocarbons. The Journal of Physical Chemistry A, 103, 7487-7495.

[30]   Morell, C., Grand, A. and Toro-Labbé, A. (2005) A New Dual Descriptor for Chemical Reactivity. The Journal of Physical Chemistry A, 109, 205-212.

[31]   Bultinck, P., Clarisse, D., Ayers, P.W. and Carbo-Dorca, R. (2011) The Fukui Matrix: A Simple Approach to the Analysis of the Fukui Function and Its Positive Character. Physical Chemistry Chemical Physics, 13, 6110-6115.

[32]   Melin, J., Ayers, P.W. and Ortiz, J.V. (2007) Removing Electrons Can Increase the Electron Density: A Computational Study of Negative Fukui Functions. The Journal of Physical Chemistry A, 111, 10017-10019.

[33]   Padmanabhan, J., Parthasarathi, R., Elango, M., Subramanian, V., Krishnamoorthy, B.S., Gutierrez-Oliva, S., Toro-Labbe, A., Roy, D.R. and Chattaraj, P.K. (2007) Multiphilic Descriptor for Chemical Reactivity and Selectivity. The Journal of Physical Chemistry, 111, 9130-9138.

[34]   Parr, R.G., Donnelly, R.A., Levy, M. and Palke, W.E. (1978) Electronegativity: The Density Functional Viewpoint. The Journal of Chemical Physics, 68, 3801-3807.

[35]   Mulliken, R.S. (1934) A New Electroaffinity Scale: Together with Data on Valence States and on Ionization Potential and Electron Affinities. The Journal of Chemical Physics, 2, 782-793.

[36]   Zhou, Z. and Parr, R.G. (1989) New Measures of Aromaticity—Absolute Hardness and Relative Hardness. Journal of the American Chemical Society, 111, 7371-7379.

[37]   Parr, R.G. and Gazquez, J.L. (1993) Hardness Functional. The Journal of Physical Chemistry, 97, 3939-3940.

[38]   Pearson, R.G. (1987) Recent Advances in the Concept of Hard and Soft Acids and Bases. Journal of Chemical Education, 64, 561-567.

[39]   Parr, R.G., Szentpály, L.V. and Liu, S. (1999) Electrophilicity Index. Journal of the American Chemical Society, 121, 1922-1924.

[40]   Koopmans, T. (1934) über die zuordnung von wellen funktionen und eigenwerten zu den einzelnen elektronen eines atom. Physica, 1, 104-113.

[41]   Lim, C., Bashford, D. and Karplus, M. (1991) Absolute pKa Calculations with Continuum Dielectric Methods. The Journal of Physical Chemistry, 95, 5610-5620.

[42]   Topol, I.A., Tawa, G.J., Burt, S.K. and Rashin, A.A. (1997) Calculation of Absolute and Relative Acidities of Substituted Imidazoles in Aqueous Solvent. The Journal of Physical Chemistry A, 101, 10075-10081.

[43]   Remko, M. (2003) Theoretical Study of Molecular Structure and Gas-Phase Acidity of Some Biologically Active Sulfonamides. The Journal of Physical Chemistry A, 107, 720-725.

[44]   Assoma, A.B., Bede, A.L., Kone, M. and N’Guessan, Y.T. (2010) Theoretical Study of Stability, Tautomerism, Equilibrium, Constants (pkT), Activation Energies and Acidity of 6-Thioxanthine in Gas and Aqueous Phase by the Ab Initio Method and Functional Density Theory Calculations. European Journal of Scientific Research, 44, 337-354.

[45]   Klicic, J.J., Friesner, R.A., Shi-Yi, L. and Wayne, C.G. (2002) Accurate Prediction of Acidity Constants in Aqueous Solution via Density Functional Theory and Self-Consistent Reaction Field Methods. The Journal of Physical Chemistry A, 106, 1327-1335.

[46]   Keith, J.A. and Carter, E.A. (2012) Quantum Chemical Benchmarking, Validation, and Prediction of Acidity Constants for Substituted Pyridinium Ions and Pyridinyl Radicals. Journal of Chemical Theory and Computation, 8, 3187-3206.

[47]   Remko, M. and von der Lieth, C.W. (2004) Theoretical Study of Gas-Phase Acidity, pKa, Lipophilicity, and Solubility of Some Biologically Active Sulfonamides. Bioorganic & Medicinal Chemistry, 12, 5395-5403.

[48]   Vianello, R. and Maksi?, Z.B. (2005) High Acidity of Polycyano Azatriquinanes— Theoretical Prediction by the DFT Calculations. Tetrahedron Letters, 46, 3711-3713.

[49]   Roy, R.K. (1999) On Non-Negativity of Fukui Function Indices. The Journal of Chemical Physics, 110, 8236-8245.