INTERCONVERSION EQUILIBRIUM OF PARA–META ISOMERIZATION
449
of di-tert-butylbenzenes with the meta interaction of culated by the Ambrose–Walton method [24] with the
the tert-butyl substituents in TBB molecules.
critical parameters (Table 6) and normal boiling points
(NBP) (Table 2).
o
In this work, ∆ H (g, í ) of reaction (4) serves as
r
m
av.
These systems exhaust the variety of alkylbenzenes
for which not only their chemical para–meta isomeriza-
tion equilibrium have been studied, but also experimen-
tal data on the vapor pressure of the reactants or, at
least, on the NBP of the compounds are available. It is
obvious that, as the structure of the reactant molecules
becomes more complex, it becomes less likely to find
the complete set of experimentally obtained informa-
tion for the description of chemical equilibria in both
liquid- and gas-phase transformations. In this case, an
approach with a minimum of initial information must
inevitably be employed for the prediction of equilib-
rium in the transformations in question.
an important supplement to the data discussed above
and allows the general conclusion that the enthalpy of
para–meta isomerization and, hence, the value of the
enthalpic effect of interaction between 1,3-substituents,
which stabilizes the molecule, depends on the effective
size of the substituents.
An additional support for this conclusion is pro-
vided by the information that is presented in Figs. 1–4
and relates to the contribution of reactions (1)–(4) to
o
m
o
m
∆rS (g, í ). In fact, ∆ H (g, í) of the reaction is
i
r
o
directly related to ∆ S (g, í) at close values of K
r
m
p
(
Table 1) for transformations of the same type. From the
We have shown that a fairly accurate prediction of
the equilibrium is possible for the transformations of
the given type, which includes the calculation (quan-
tum-chemical and by statistical thermodynamics) of
o
data presented in Table 3, it follows that the enthalpy
o
o
∆r H (g, í) and the entropy ∆ S (g, í) of reactions
m
r
m
–1
–1
(1)−(4) regularly change, from 5.5 to 3.1 J mol K and
from –0.47 to 1.23 kJ/mol, respectively, upon passing contributions to ∆
from xylenes to tert-butylbenzenes.
The ratio between the individual contributions to the methods verified for this group of compounds.
S
(g, í) of the reaction of interest
r
m
and the estimation of the NBP of the reactants using the
entropy of the reaction changes in this case (Figs. 1–4).
An increase in the contribution of the vibrational term
The potential of this approach is illustrated by the
material presented in Table 7, which collates the results
of the calculation of the liquid-phase equilibrium con-
stants with the experimental data for the selected sys-
tems. Note that there are no reliable ebulliometric data
on the saturated vapor pressure and normal boiling
point for any component of reactions (6)–(9).
o
leads to a decrease in ∆ S (g, í).
r
m
The entropies of the compounds were calculated by
means of the statistical thermodynamics procedure
using the force-field molecular mechanics method
(Allinger MM2 modernized field) to determine the
As a source of information on the NBP of com-
pounds, we used data obtained using GLC [32]. The
critical parameters were calculated on the basis of the
Randic molecular connectivity indices using the
method described in [30], and the critical pressures
were calculated by the Lydersen method [31]. The val-
ues of the critical properties and the NBP of the com-
pounds used in the calculation of the vapor pressure via
the Ambrose–Walton method [24] are given in Table 6.
geometry of molecules, and rotation barriers, which are
necessary for the calculation of the entropy components
due to both the rotation of the molecule as a whole and
the internal rotation of the groups. The vibration term
of the entropy was calculated by means of the density
functional theory (DFT) with the use of the B3LYP/6-
3
11++G92d,2p) basis set. The application of the
molecular-mechanics method to the calculation of the
rotational terms is caused by the necessity to reduce the
computation time in the matching of the geometric
The procedure used to calculate the equilibrium
parameters used in the entropy terms for the rotation of constants corresponds to the sequence of the presenta-
the molecule as a whole and the internal rotation.
tion of the material in Table 7. Using the statistical ther-
modynamics routine, the rotational (of the molecule as
a whole and groups in the molecule) and vibrational
The coefficients for the equation
o
or
ir
vib
ln∆rS = B(1000/T) + Cln(1000/í) + D(1000/í) (2)
(S + S + S ) contributions to the entropy of the com-
i, m
pounds were calculated; the contribution to the entropy
change (S + S + S ) –(S + S + S ) of the reaction
of the temperature dependence (in the interval 273–598 K),
or
ir
vib
or
ir
vib
fin
in
o
i, m
where ∆rS are the contributions to the entropy due to was determined; the enthalpies of the reactions were
calculated according to Eq. (3) on the basis of the
external and internal rotation or the vibrational term, are
presented in Table 5.
o
obtained data; and, ∆ S (g, í), Kp (calc.), and Kx (calc.)
—
r
from the calculated vapor pressure—were subsequently
found and compared with the experimental values.
A chemical experiment on isopropyltoluenes was
performed in [6, 7] in the presence of AlCl in the liquid
3
In this procedure, only the form of the equation
phase for quite a broad temperature range (293–473 K).
However, this range is covered by experimental
o
or
ir
vib
or
ir
∆
H (g, í) = 0.3253 × [(S + S + + S ) –(S + S +
r
fin
vib
S ) ] – 0.5379 (3) needs to be discussed.
−ídata only for 383–423 K [17]; moreover, these data
in
are not very precise (Table 2). In accordance with the
At present, we consider the linear form to be the
recommendations [1], we used the vapor pressures cal- most appropriate; nonetheless, it should be mentioned
PETROLEUM CHEMISTRY Vol. 48 No. 6 2008