Table 3 AM1 (first row) and 6-311G*//3-21G (second row) standard
thermodynamic properties calculations, variation of the phenylaniline
dihedral angle, ∆φЊ, and absolute dipole moment changes, ∆µ, for
the tautomeric reaction of the N-salicylidene-p-X-aniline–methanol
hydrogen bond complexes
Table 2 AM1 Heat of formation and 6-311G*//3–21G MO energy
calculations for enol (e) and keto (k) N-salicylidene-p-X-aniline and
N-salicyl-p-X-aniline–methanol hydrogen bond complexes (hb)
AM1/kJ molϪ1
6-311G*//3-21Ga/au
a
∆GtЊau
∆HtЊau
T∆StЊau
∆φЊ(φketo
)
∆µ/debyeb
Ee
121.75
Ϫ116.98
Ϫ129.66
Ϫ12.63
Ϫ627.880314
Ϫ742.887668
Ϫ742.902068
Ϫ0.0144
Ee ϩ Emet
Eehb
∆Eehb = Eehb Ϫ (Ee ϩ Emet
-NMe2
-OMe
-Me
2.702
5.916
3.414
7.739
1.452
7.142
1.464
7.167
1.460
8.217
3.243
11.180
3.807
12.247
3.653
5.535
4.485
7.037
2.686
6.452
2.686
6.661
2.569
5.828
4.201
10.263
4.778
11.255
0.825 Ϫ7.4 (26.5)
Ϫ0.525 Ϫ4.0 (22.6)
1.148 Ϫ3.1 (28.8)
2.637
2.332
2.250
1.792
1.954
1.389
1.618
1.018
1.222
1.170
1.393
1.028
)
37.81 kJ molϪ1
Ϫ627.869401
Ϫ742.876755
Ϫ742.899528
Ϫ0.0228
Ek
Ek ϩ Emet
Ekhb
140.37
Ϫ98.32
Ϫ126.98
Ϫ28.66
Ϫ0.605
0.792
Ϫ1.045
0.794
Ϫ0.805
0.680
Ϫ0.996
0.961
2.0 (19.5)
4.1 (21.7)
9.4 (18.6)
5.4 (21.5)
10.3 (19.0)
12.8 (17.5)
18.7 (11.7)
12.7 (13.6)
19.1 (11.8)
∆Ekhb = Ekhb Ϫ (Ek ϩ Emet
)
-H
Ϫ59.79 kJ molϪ1
-COMe
-CN
a The ab initio calculation including the ZPE thermal corrections.
Ϫ2.247
1.108
-NO2
18.8 (12.8) Ϫ1.279
24.3 (9.3) Ϫ0.835
Ϫ2.554
a HF/3-21G calculations. b HF/6-311G*//3-21G calculations.
it is not possible to observe the visible band at 450 nm associ-
ated with the keto species, which is probably due to the fact that
in a non-polar solvent the energy between both species is about
20 kJ molϪ1 or more. The equilibrium should then be shifted to
the enol form and the absorption spectra should present only
the one band associated with this species. On the other hand,
in a protic solvent, the energy difference is reduced to 4 kJ
molϪ1, and a small fraction of the keto species is present in the
solution, and therefore a weak absorption band in the visible
region appears.
Thermodynamics calculations. For each compound under
study, we have calculated the total energy for the optimized
hydrogen bonded complex of the enol and keto species. In the
evaluation of the thermodynamics properties HЊ(hb), GЊ(hb)
and SЊ(hb), we have used the option “thermo”21 and “fre-
quency”19 in both AM1 and ab initio, calculations, respectively.
For the enol-imine→keto-amine tautomeric reaction, the
standard properties (∆PtЊau(hb)), were evaluated through ∆PtЊau-
(hb) = PЊ(khb) Ϫ PЊ(ehb), where PЊ represents any thermo-
dynamic property, HЊ, GЊ, SЊ.
In Table 3, we present the calculated ∆GtЊau(hb), ∆HtЊau(hb)
and T∆StЊau(hb) values for the tautomeric reaction at 298 K.
We observe that ∆GtЊau(hb) and ∆HtЊau(hb) values, calculated by
means of the 6-311G*//3-21G basis (second row), are system-
atically larger than the AM1 ones (first row). The greatest
difference between them was found in the nitro-substituted
compound.
In order to compare the calculated thermochemical proper-
ties with the experimental ones, we show in Figs. 4-B and 4-C
the calculated and experimental values against the Hammett
parameter. In Fig. 4-B we show the 6-311G*//3-21G ∆GtЊau(hb)
values. They are quite close to the experimental data with a
small negative deviation of 0.13 kJ molϪ1 (on average). On the
other hand, the AM1 results present a larger negative deviation
of Ϫ5.4 kJ molϪ1 (on average) with respect to the experi-
mental data, which is due to the positive value of the calculated
entropic contribution (see Table 3).
The AM1 and ab initio ∆HtЊau(hb) values are shown in the
Fig. 4-C. It can be seen that the AM1 calculations display a
negative deviation (1.9 kJ molϪ1 on average) and the ab initio
values predict a positive deviation (2.3 kJ molϪ1, on average),
with respect to the experimental values. Furthermore, we note
that the ab initio data follow the same trend with a good linear
correlation (r = 0.94). In contrast to this, the AM1 results do
not exhibit an important effect of the substituent, leading to a
poor linear correlation (r = 0.25).
Fig. 6 Optimized geometrical structures obtained from AM1 and
ab initio calculations.
respectively. Ee and Ek are the total energies of the optimized
geometry of the non-bonded enol and keto forms and Emet is
the energy of the optimized geometry of a methanol molecule.
To compare the AM1 and ab initio results, we have added, in
ab initio calculations, the thermal ZPE correction enthalpy,18,19
at 298.15 K to the total energy. This is a standard procedure,
and requires calculation of the force field. With the knowledge
of the force field and the moments of inertia even the cal-
culation of the partition function and the thermodynamic data
at room temperature is straightforward.20 The standard enthalpy
of the enol→keto tautomeric reaction is given by ∆HtЊau =
Ek Ϫ Ee for the non-bonded fragment and ∆HtЊau(hb) = Ekhb
Eehb for the H-bonded complex.
=
In Table 2 we display the AM1 and 6-311G*//3-21G
results of N-salicylidene-p-X-aniline and the N-salicylidene-p-
X-aniline–methanol hydrogen-bonded complexes. We note that
the formation energy of the keto–methanol complex (Ekhb
)
is higher than the formation energy of the enol–methanol
complex (Eehb). This higher stabilization energy of the keto
complex reduces the ∆HtЊau value, from 18.6 to 2.7 kJ molϪ1 in
AM1, and from 28.7 to 6.7 kJ molϪ1 in the ab initio calculation.
The energy values are in agreement with the absorption and
fluorescence excitation spectra shown in Fig. 2. In cyclohexane
J. Chem. Soc., Perkin Trans. 2, 2001, 1124–1129
1127