Topological and thermodynamic investigations of binary mixtures: Molar excess volumes, molar excess enthalpies and isentropic compressibility changes of mixing
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Fluid Phase Equilibria 235 (2005) 42–49
Topological and thermodynamic investigations of molecularinteractions in binary mixtures: Molar excess volumes
and molar excess enthalpies
Sanjeev Makena,∗, Bal Raj Deshwalb, Renu Chadhac, Anud,Krishan Chander Singhd, Hwayong Kime, Jin-Won Parka
a Department of Chemical Engineering, Yonsei University, Seoul 120-749, Republic of Koreab Korea Institute of Energy Research, Daejon 305 600, Republic of Korea
c Department of Pharmaceutical Science, Punjab University, Chandigarh, Indiad Department of Chemistry, M. D. University, Rohtak 124001, Haryana, India
e School of Chemical Engineering, Seoul National University, Seoul 151-744, Republic of Korea
Received 18 November 2003; received in revised form 8 June 2005; accepted 22 June 2005
Molar excess volumes and molar excess enthalpies of butyl acetate (i) with cyclohexane or benzene or toluene oro-, m- or p-xylene (j)inary mixtures have been measured dilatometrically and calorimetrically over the entire composition range at 308.15 K. The obsave also been analyzed in terms of graph theoretical approach. The analysis ofVE data by graph theoretical approach suggests thatcetate in pure state exists as associated entity and (i + j) mixtures are characterized by the presence of (i–j) molecular entity. It has furtheeen observed thatVE andHE values calculated by this approach agree well with the corresponding experimental values. The preolecular entity is further confirmed by IR study of (i + j) mixture.2005 Elsevier B.V. All rights reserved.
A number of experimental as well as theoretical stud-es on thermodynamic properties for alkyl ester withlkanes [1–9], aromatic hydrocarbons[8–17], alcohols
18,19], chloroalkane[20] and acetonitrile[16,21]have beeneported in literature. Among them, systemic data for alkylster + aromatic polar solvent are relatively rare[10–17].lkyl esters are characterized by dipole–dipole interac-
ions in the pure state[1–3]. The degree and strength ofipole–dipole interactions decreases with the increasing sizef alkyl group in the esters. The additions of inert solvents
ike alkane generally break the orientation order of pure
alkyl esters to give the positive value of excess thermnamic functions like excess molar volume,VE, excess molaenthalpy,HE, and excess molar Gibb’s free energy,GE [4,5].However, these values become very less and eventive in the mixtures of alkyl esters + aromatic hydrocarb[10,11,14,15]. Specific interactions of dipole–induce diptype are postulated to account for such behavior[14,15].
Moreover, it has been revealed[22–26]that that graph theoretical approach, based on molecular connectivity param[27] of third degree, could be of great use not only in euating theirVE andHE but also in understanding the natof molecular interactions between the components of bmixtures. These considerations prompted us to carrysystemic study on the thermodynamic properties of aesters with aromatic hydrocarbons. In first paper ofseries, we report the measuredVE andHE data for buty
S. Maken et al. / Fluid Phase Equilibria 235 (2005) 42–49 43
Table 1Densities,ρ, and refractive indices,nD, of the pure components at 308.15 K
Substance ρ (kg m−3) nD
This work Literature This work Literature
Butyl acetate 865.45 865.4a 1.3908 1.3910a
Cyclohexane 764.43 764.46b 1.4205 1.4201c
Benzene 862.93 862.95b 1.4921 1.49170d
Toluene 852.80 852.85b 1.4885 1.48870d
o-Xylene 867.41 867.38b 1.5031 1.50295e,f
m-Xylene 851.53 851.57b 1.4946 1.49464e,f
p-Xylene 847.91 847.87b 1.4879 1.4881d
a Ref. [21].b Ref. [31].c Ref. [32].d Ref. [33].e Ref. [28].f At 298.15 K.
acetate + cyclohexane or benzene or toluene oro-, m- or p-xylene at 308.15 K and their interpretation in terms of graphtheoretical approach.
2. Experimental
Butyl acetate (BA) (Merck) was dried over anhydrousmagnesium sulphate and then fractionally distilled[28]. Themiddle fraction of distilled BA was then dried over type0.3 nm molecular sieves (Merck) in an amber colored bottlefor several days before use. Cyclohexane, benzene, tolueneand xylenes (Merck) were purified with standard procedure[28,29]. The purities of the purified samples were checked bymeasuring their densities and refractive indices at 308.15 Kas described earlier[30] and these compared well with theliterature values[21,28,31–33]as shown inTable 1.
Molar excess volumes,VE, for the binary mixtures havebeen measured dilatometrically at 308.15 K in the mannerdescribed elsewhere[34]. The temperature of the thermostatwas controlled within±0.01 K. The uncertainties in the mea-suredVE values are±0.5%.
Molar excess enthalpies (HE) were determined with a heatflux calorimeter (model C-80, Setaram, France) at 308.15 Kas described elsewhere[35]. The uncertainties in the mea-s E −1
wherexi is the mole fraction of BA andXn (n= 0–3) areadjustable parameters. These parameters were evaluated bymethod of least squares and are given along with standarddeviations,σ(X), in Table 4. TheVE values for BA + benzeneat 308.15 K have been reported by Oswal[14]. Our experi-mental curve for BA + benzene system is in agreement withthat reported by him.VE values ofo-,m- orp-xylene systemsat 303.15 K have also been reported by Ramachandran et al.[15]. Our experimental values for these systems are lowerthan their values at each mole fraction because reported val-ues are at lower temperature. However, the general shape ofthe curves is same. We are unaware of any previously pub-lishedHE data for the present systems at 308.15 K with whichto compare our results. TheHE values for all the binary sys-tems are negative over the entire composition range except
F(3r
uredH values are about±0.1 J mol .
. Results and discussion
The VE and HE for BA + cyclohexane or benzeneoluene oro-, m-, p-xylene as a function of mole fractioi , at 308.15 K are reported inTables 2 and 3and shown inigs. 1 and 2. The results were fitted to the Redlick-Kisquation
Vn, σ(VE) are in cm3 mol−1 andHn, σ(HE) are in J mol−1.
for cyclohexane andm-xylene systems. At equimolar com-position,HE values vary in the order:
cyclohexane> m-xylene> benzene> o-xylene
> p-oxylene> toluene.
On the other hand, whileVE values for BA + toluene oro-,m-orp-xylene systems are negative over the whole compositionrange, these are positive for BA + cyclohexane or benzenesystems. At equimolar composition,VE values follow thesequence:
cyclohexane> benzene> toluene≈ o-xylene
> m-xylene> p-xylene.
At the simplest qualitative level, the observedHE andVE
values may be attributed to the resultant of two opposingeffects. The positive contribution toHE andVE values arisesfrom the breaking and stretching of dipole–dipole interac-tions in self-associated BA and negative contribution arisesdue to the formation of attractive interactions between polarcarboxyl group of BA and�-electrons of aromatic hydro-carbons. As theHE and VE values for cyclohexane + BAmixture are highly positive, it indicates the dissociation ofdipole–dipole interactions and absence of attractive interac-tions. When cyclohexane is replaced by benzene,HE andVE E es alsot and� isi nsityo omesm es ofH oft ausef nd itse the
presence of steric hindrance between the two methyl groupsof xylenes and the alkyl groups of BA, which restrict theproper orientation of these molecules and obstruct the acetatetoward the ring thus making interactions weaker. Among thexylenes,p-xylene is symmetrical molecule, thus it offers leaststeric hindrance. TheHE andVE values for this system aremost negative[12]. The substitution of methyl group inmeta-position seems to be most sterically hindered and thusHE isslightly positive for this system.
3.1. Conceptual aspects of graph theoretical approachand results
According to mathematical discipline of graph theory[36], if the atoms in a structural formula of a molecule arerepresented by vertices and bonds joining them by edges,then the resulting graph describes the totality of informationcontained in that molecule[27,36–39]. Consequently, ifδv
m,δvn, etc. represent the degrees ofmth andnth vertices of the
graph of a molecule (shown inFig. 3 around atoms), thenconnectivity parameters of third degree,3ξi , is defined[39]by Eq.(2) and shown inFig. 3:
3ξi =∑
m<n<o<p
(δvmδv
nδvoδ
vp)−0.5 (2)
wm ,Zb
δ
F ard-i d bye
tso s
values decrease considerably andH becomes negativhowing the breaking of dipole–dipole interactions andhe formation of interactions among polar ester groups-electrons on benzene ring[12]. When a methyl group
ntroduced in benzene (as in toluene), the electron def �-electron cloud increases and these interactions becore stronger and this should lead to more negative valuE andVE than benzene + BA mixture. The introduction
wo methyl groups in benzene (as in xylenes) should curther enhancement of these attractive interactions ahould further yield more negative values ofHE andVE. How-ver, it is not observed experimentally. It may be due to
hereδvm, etc. reflect explicitly the valency ofmth vertex in
olecular graph ofi and is related[39] to maximum valencym, and number of H-atoms,hm, attached tomth, etc. vertexy the relation:
vm = Zm − hm (3)
urther, Kier[27] has suggested that the information regng effect of branching in the molecules can be obtainevaluation of3ξi of molecules.
SinceVE of an (i + j) mixture reflects interactional effecn the packing of molecules and as (3ξi)−1 of i determine
46 S. Maken et al. / Fluid Phase Equilibria 235 (2005) 42–49
Fig. 3. Structures of various components of binary mixture along with3�′ values.
[23] the effectiveness with which the molecular surface ofiinteracts with that of anotheri molecule, the interactional partof the molar volume of purei has been taken to (3ξi)−1. Theideal interactional molar volume of an (i + j) mixture would
then be proportional to∑
xi(3ξi)−1
, wherexi and (3ξi) denotemole fraction and connectivity parameter of third degree,respectively, ofi. The interactional molar volume of an (i + j)mixture should also be proportional to 1/(3ξ)m of the mix-ture. If (3ξ)m of the real mixture in the molecular graph isexpressed by
∑{x1(3ξi)m}, where (3ξi)m denotes3ξ of i in
the mixture and if the proportionality constantαij is assumeto be the same for the mixture and its pure components then,
VE may be expressed by[23]
VE = αij
[{∑xi(
3ξi)m}−1 −
∑ xi
3ξi
](4)
whereαij is the constant characteristic of the (i + j) mixtureand can be evaluated using equimolar experimentalVE value.
As the degree of association ofi andj is not known in themixture and pure state, we regarded (3ξi) and (3ξi)m (i = i orj) as adjustable parameters and evaluated them by employingVE data to Eq.(4). Only those values of (3ξi) and (3ξi)m wereretained that best reproduced theVE data. Various (3ξi) or(3ξi)m (i = i or j) parameters along withαij are recorded in
S. Maken et al. / Fluid Phase Equilibria 235 (2005) 42–49 47
Table 5Comparison of smoothenedVE andHE values with their corresponding values calculated from graph theory along with the various3ξ values, interaction energyparametersχ′
αij andσ(VE) are in cm3 mol−1; χ∗ij , χ12 andσ(HE) are in J mol−1.
Table 5along withVE values obtained from Eq.(4) for thevarious binary mixtures as a function ofxi and are also com-pared with their experimental values. Examination ofTable 5reveals thatVE values agree well with their correspondingexperimental values and thus (3ξi) or (3ξi)m (i = i or j) valuescan be relied upon to extract information about the state ofpure components as well as their mixture.
A number of structures (I–VIII ) were then assumed fori andj and their (3ξ′
i) values were calculated from structuralconsiderations of Eq.(2). These values were then comparedwith corresponding values, (3ξi) obtained from Eq.(4). Anystructure or combination of structures that yielded (3ξ′
i) valueswhich compared well with the corresponding (3ξi) obtained
from Eq. (4) were considered as representative structure ofthat component. It was assumed that BA exists as molec-ular entitiesI and II . The calculated values of (3ξ′
i) forthese molecular entities were 0.889 and 1.880, respectively.The (3ξi) values for BA in BA (i) + cyclohexane or benzeneor toluene oro-, m- or p-xylene (j) mixtures from Eq.(4)are 1.261, 1.261, 1.261, 1.661, 1.261 and 1.666, respec-tively (Table 5). It suggests that BA exists as a mixture ofmonomer and dimer (average (3ξ′
i) = 1.384). Further, (3ξ′i)
values of 0.666, 0.840, 1.426, 1.174, 1.250 and 1.500 forbenzene, toluene,o-, m-, p-xylene and cyclohexane (molec-ular entitiesIII –VIII in Fig. 3) suggest that they exist asmonomers.
48 S. Maken et al. / Fluid Phase Equilibria 235 (2005) 42–49
Fig. 4. IR spectra of butyl acetate,p-xylene and their equimolar mixture.
(3ξi)m values were next evaluated to see the postulationof BA in aromatic hydrocarbons (Fig. 3). For this purpose, itwas assumed that if BA (i) in benzene or toluene or xyleneis assumed to exist as molecular entityIX (3ξ′
i = 1.20), then(3ξi)m values of 1.261–1.666 (Table 5) for these mixturessuggest that BA (i) + aromatic hydrocarbons (j) mixtures con-tain molecular entityIX . For calculating (3ξ′
i)m value, it isassumed that molecular entityIX is characterized by inter-actions between�-electrons spilling over O–C–O of BA anddelocalized�-electrons cloud of aromatic hydrocarbon. Thepresence of molecular entityIX in the studied mixtures thensuggests that addition of BA to aromatic hydrocarbons shouldhave influenced symmetric vibrations of carbonyl group ofBA and also the ring vibrations of aromatic hydrocarbons.In view of this, we analyzed IR spectra of pure BA and itsequimolar mixture withp-xylene. It was observed from theIR spectra (Fig. 4) that BA andp-xylene in their pure stateshowed characteristic vibrations at 1744 cm−1 (C O stretch-ing), 1069 and 1034 cm−1 (C O stretching) and 1516, 1460and 1382 cm−1 (ring vibrations)[41]; while their mixtureshowed characteristic vibrations at 1686, 1103, 1063 and1541, 1470, 1370 cm−1, respectively. Thus, IR spectra pro-vide an additional support to the presence of molecular entityIX (Fig. 3) in these mixtures.
To study the energetic of the interactions operatingb umedt
((
(
np e top
�
whereSj is the surface fraction ofj involved in (i–j) contactand is defined[42,43]by:
Sj = xjVj∑xiVi
(6)
So that Eq.(5) reduces to
�Ha = xixjχijVj∑xiVi
(7)
Further, ifχii andχ12 are molecular interaction parametersfor (i–i) and specific interactions betweeni andj components,then enthalpy change due to processes (b) and (c) are[42–44]given by:
�Hb = x2i xjχiiVj∑
xiVi
(8)
�Hc = xix2jχ12Vj∑xiVi
(9)
The total enthalpy change due to processes (a–c) is given by:
HE =[
xixjVj∑xiVi
][χij + xiχii + xjχ12] (10)
SinceVj /Vi = (3ξi /3ξj), Eq.(10) then reduces to:
H
FE
H
E
ppwx
ap hisg n ofEc
H
E
w( lyuH yi
nsb
etween the components of these mixtures, it was asshat (i + j) mixtures formation involves:
a) the formation of unlike contact betweeni andj;b) unlike contact then leads to depolymerization of (i) to
yield monomer;c) monomer ofi andj then formi–j molecular entity.
Consequently, ifχij is the molar enthalpy interactioarameter of (i–j) contacts, then change in enthalpy durocess (a) is expressed[42,43]by:
Ha = xiχijSj (5)
E =[
xixj(3ξi/3ξj)
xi + xj(3ξi/3ξj)
][χij + xiχii + xjχ12] (11)
or these mixtures, if it is assumed thatχij = χii = χ∗ij, then
q.(11) reduces to:
E =[
xixj(3ξi/3ξj)
xi + xj(3ξi/3ξj)
][(1 + xi)χ
∗ij + xjχ12] (12)
q.(12)contains two unknown parametersχ∗ij andχ12. These
arameters were calculated employingHE data at two com-ositions (xi = 0.4 and 0.5) for various (i + j) mixtures andere subsequently used to evaluateHE at other values ofi . Such values ofHE along with parametersχ∗
ij and χ12
re recorded inTable 5and calculatedHE found to com-are well with their corresponding experimental values. Tives additional support to assumptions made in derivatioq.(10)–(12). For BA (i) + cyclohexane mixture (j), χ12 = 0,onsequently Eq.(12) for this mixture reduces to:
E =[
xixj(3ξi/3ξj)
xi + xj(3ξi/3ξj)
][(1 + xi)χ
∗ij] (13)
q. (13) contains only one unknown parameterχ∗ij and it
as calculated employingHE data only at one compositionxi = 0.5) for various (i + j) mixtures and were subsequentsed to evaluateHE at other values ofxi . Such values ofE were tabulated inTable 5and comparison was not ver
mpressive.Values ofχ12, a measure of extent of specific interactio
etweeni andj for BA + benzene or toluene orp-xylene are
S. Maken et al. / Fluid Phase Equilibria 235 (2005) 42–49 49
−0.17,−9.44 and−7.18 J mol−1 which again indicate theincrease in interactions with increase in electron donatingtendency due to addition of methyl groups to benzene, slightlyless negative (−7.18 J mol−1) value forp-xylene may be dueto steric hindrance offered by its two methyl groups.
List of symbolsBA butyl acetateGE molar excess Gibb’s free energyhm number of H-atomsHE molar excess enthalpyHn adjustable parameters of Eq.(1)Sj surface fraction ofjth componentVE molar excess volumeVi , Vj molar volume ofith andjth componentsVn adjustable parameters of Eq.(1)xi , xj mole fraction ofith andjth componentsZm maximum valency
Greek lettersαij constant characteristic of (i + j) mixture in Eq.(4)δvm degree ofmth vertex of the graph of a molecule
σ standard deviation3ξi , (3ξi)m, (3ξ′
i) connectivity parameters of third degree ofith component in pure state, mixture from Eqs.(4)
χ er-
A
nksK ita-t AjitM ov-e tateC
R
004)
.
her-
004)
ieu,
2001)
[7] M. Pintos, R. Bravo, M.C. Baluja, M.I. Paz Andrade, G. Roux-Desgranges, J.P.E. Grolier, Can. J. Chem. 66 (1988) 1179–1186.
[8] E. Jimenez, L. Romanı, E. Wilhelm, G. Roux-Desgranges, J.P.E.Grolier, J. Chem. Thermodyn. 26 (1994) 817–827.
[9] J.-P.E. Grolier, D. Ballet, A. Viallard, J. Chem. Thermodyn. 6 (1974)895–908.
[10] S. Delcros, E. Jimenez, L. Romanı, A.H. Roux, J.P.E. Grolier, H.V.Kehiaian, Fluid Phase Equilib. 111 (1995) 71–88.
[11] G. Figurski, U. Emmerling, D. Nipprasch, H.V. Kehiaian, FluidPhase Equilib. 137 (1997) 173–183.
[12] J.M. Resa, C. Gonzalez, R.G. Concha, M. Iglesias, Phys. Chem.Liqs. 42 (2004) 493–520.
[13] D. Borchardt, M. Opel, H.-J. Rademacher, M.T. Ratzsch, Z. Phys.Chem. (Leipzig) 262 (1981) 479–484.
[14] S. Oswal, J. Chem. Eng. Data 66 (1988) 111–115.[15] D. Ramachandran, K. Rambabu, K.M. Krishan, P. Venkateshwarlu,
G.K. Raman, J. Chem. Eng. Data 40 (1995) 815–817.[16] T. Ohta, T. Kinoshita, I. Nagata, J. Chem. Eng. Data 28 (1983)
36–39.[17] J. Gmehling, B. Meents, Int. DATA Ser. Data Mixtures Sel. Ser. A
20 (1992) 177–184.[18] A. Blanco, M.J. Ortega, J. Chem. Eng. Data 43 (1998) 638–645.[19] J.M. Resa, C. Gonzalez, M. Juez, S.O. de Landaluce, Fluid Phase
Equilib. 217 (2004) 175–180.[20] G. Avedis, J.P.E. Grolier, J. Fernandez, H.V. Kehiaian, J. Solution
Chem. 23 (1994) 135–152.[21] T.N. Aminabhavi, K. Banerjee, J. Chem. Eng. Data 43 (1998)
514–518.[22] P.P. Singh, Thermochim. Acta 66 (1983) 37–73.[23] P.P. Singh, S. Maken, Pure Appl. Chem. 66 (1994) 449–454.[24] P.P. Singh, S. Maken, R. Malik, Indian J. Chem. 36A (1997) 46–51.[ 9.[[ .),
New
[ urth
[ d.,
[ Han,
[ 9.[ 98)
[ ata
[[ . 29
[ 9.[ mic
[[ ug
[ Aca-
[[[ 989)
and(2)ii, χij, χ
∗ij, χ12 molecular interaction parameters for diffent binary interactions in Eqs.(10)–(12)
cknowledgement
Bal Raj Deshwal and Sanjeev Maken (S.M.) thaorean Federation of Science and Technology for inv
ion under Brain Pool program, and S.M. also thanks. Sharon (IAS), Commissioner, Technical Education, Grnment of Haryana, for award of study leave from C.R. Sollege of Engineering, Haryana, India.
eferences
[1] J. Ortega, F. Espiau, F.J. Toledo, J. Chem. Thermodyn. 36 (2193–209.
[2] M. Chaar, J. Ortega, F.J. Toledo-Marante, C. Gonzalez, J. ChemThermodyn. 33 (2001) 689–711.
[3] F.J. Toledo-Marante, J. Ortega, M. Chaar, M. Vidal, J. Chem. Tmodyn. 32 (2000) 1013–1036.
[4] J. Ortega, F. Espiau, R. Dieppa, Fluid Phase Equilib. 215 (2175–186.
[5] H.V. Kehiaian, R. Bravo, M.P. Barral, M.I. Paz Andrade, R. GuFluid Phase Equilib. 17 (1984) 187–216.
[6] J. Ortega, C. Gonzalez, S. Galvan, J. Chem. Eng. Data 46 (904–912.