Ana C Gómez Marigliano^{*} and Ruth N Medina Naessens  
Departamento de Física, Facultad de Ciencias Exactas y Tecnología, Universidad Nacional de Tucumán, Avenida Independencia 1800, 4000San Miguel de Tucumán, Argentina  
Corresponding Author :  Ana C. Gómez Marigliano Departamento de Física Facultad de Ciencias Exactas y Tecnología Universidad Nacional de Tucumán Avenida Independencia 1800 4000San Miguel de Tucumán, Argentina, Email: agomezmarigliano@herrera.unt.edu.ar 
Received April 20, 2015; Accepted May 13, 2015; Published May 15, 2015  
Citation: Marigliano ACG, Naessens RNM (2015) Density, Viscosity, VaporLiquid Equilibrium, Excess Molar Enthalpy and IR Spectrocopy of [Chloroform + Diisopropyl Ether (DIPE)]. J Phys Chem Biophys 5:179. doi: 10.4172/21610398.1000179  
Copyright: © 2015 Marigliano ACG, et al. This is an openaccess article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.  
Related article at Pubmed, Scholar Google 
Visit for more related articles at Journal of Physical Chemistry & Biophysics
Density, viscosity and refractive index measurements in the T=[288.15 to 303.15] K range of pure chloroform and diisopropyl ether, as well as of the binary system [x1 chloroform+(1x1) diisopropyl ether (DIPE)] over the whole concentration range at T=298.15 K were made. The experimental results for the pure components were fitted to empirical equations. Calculated values are in agreement with the experimental ones. Data of the binary mixture were further used to calculate the excess molar volume, refractive index deviations and viscosity deviations. Excess molar enthalpy at T=(298 ±1) K and vaporliquid equilibrium measurements at T=(313.15 ±0.05) K were also measured for the binary system. Activity coefficients and the excess molar Gibbs energy are calculated. This binary system shows strong negative deviations from ideality and exhibits a minimum pressure azeotrope, whose coordinates are: P=(32.90 ±0.05) kPa and x1=(0.402 ±0.002). Excess and deviation properties were fitted to the RedlichKister polynomial relation. The optimized structures and vibrational frequencies for DIPE, chloroform and their 1:1 complex were calculated by means of density functional theory (DFT) techniques using the B3LYP functional combined with the 631G(d,p) basis set. Raman and FTinfrared spectra of pure chloroform (1) and DIPE (2), as well as of the binary system were recorded.. The recorded Raman and FTIR spectra confirm the existence of these type of hydrogenbonded complexes. The ideal associated solution model was used to calculate the equilibrium constants, which compares well with results obtained from spectroscopic data and the enthalpy of formation of the hydrogen bond.
Keywords 
Association; Intermolecular complexes; Vibrational spectra; Equilibrium constants; (Vapour+liquid) equilibrium; Excess enthalpy 
Introduction 
Experimental properties of systems consisting of pairs of liquids are of great practical importance because they are used both to understand the macroscopic behavior and infer structures and interactions as for the design of equipment and methods. 
In recent times ethers have received particular attention because they are a possible additive to gasoline and also are used as solvents or reactants in the development of new products. It is for this reason that the knowledge of their properties as well as the mixtures containing is very important to optimize processes in which they are involved. 
On the other hand, excess thermodynamic functions and deviations of nonthermodynamic ones of binary liquid mixtures are fundamental to understand the interactions between molecules in this type of mixtures and used as a qualitative and quantitative guide to predict the extent and form of complex formation. 
As we assume that the system presents molecular interactions with the formation of hydrogen bond, then perform: 
1) Experimental determinations of: 
a) Density ρ , refractive index nD and viscosity η in the T=[288.15 to 303.15] K range of pure chloroform and DIPE, as well as for the binary system constituted by these two in the whole mole fraction range at T=(298.15 ±0.01) K. 
b) Excess enthalpy at T=(298 ±1) K and vaporliquid equilibrium (VLE) at T=(313.15 ±0.05) K for this binary system. 
From the experimental result for the binary system found empirical equations for density, refractive index and viscosity of both pure components as a function of temperature, as well as for the binary system as a function of composition and were calculated excess functions and deviations from ideality, excess molar volume V^{E}, refractive index deviation Δn_{D}, viscosity deviations from the ideal behavior Δη, activity coefficients γ_{i}, and excess molar Gibbs energy G^{E}. 
The excess thermodynamic properties and nonthermodynamic ones were fitted to a RedlichKister [1] type equation. 
2) The optimized structures and vibrational frequencies for DIPE, chloroform and their 1:1 complexe (as shown in Figure 1) were calculated by means of density functional theory (DFT) techniques using the B3LYP functional combined with the 631G(d,p) basis set. 
3) Raman and FTinfrared spectroscopic studies were carried out for chloroform, DIPE and the binary system in order to confirm the class and stoichiometry of nmer and to calculate the equilibrium constants, absorptivities of the monomer bands and enthalpy of formation of the hydrogen bond. 
4) We made used of the simplest and more convenient thermodynamic model: the “ideal associated solution model”. We calculate the equilibrium constant and the enthalpy of formation of the hydrogen bond.In the literature for the pure components we have found only, a single job on mixing volumes [2] and two on enthalpy of mixing [3,4]. 
Experimental Section 
Materials. 
Chloroform (p.a.) and DIPE (p.a.) were supplied by Cicarelli (Argentina) and MERCK (Germany), respectively. To verify the purity of reagents gas chromatography and IR spectroscopy was performed. Infrared spectra of pure components and mixtures were recorded on Perkin Elmer FTIR Spectrometer – Spectrum RXI, using a KBr sealed cell. A Hewlett Packard 6890 gas chromatograph with an automatic injector (Agilent G2613A) directly connected to a ChemStation (HP G2070AA) with a FID detector. The chemical purities were 0.998 and 0.995 (in mass fraction), respectively. The pure components were stored over 0.3 nm molecular sieves to prevent water absorption and their water contents were periodically checked by Karl Fischer titration using an automatic Mettler DL18 Karl Fischer titrator. These water contents never surpassed 20 ppm. 
Experimental Procedure 
The same methodology previously applied for mixtures preparation was used in this work [5]. All the weighing were preformed on an electronic balance (Mettler Toledo AG245) accurate to ±0.1 mg. The uncertainty in the mole fractions is estimated to be lower than ±1.10^{4}. 
Density and viscosity were simultaneously measured with an Anton Paar Stabinger viscometer (SVM 3000/G2) [6] calibrated by the manufacturer. This instrument is equipped with its own Peltier effect control system, which maintains temperature at ±0.01 K. The uncertainties in this work were ±0.3% of the measured viscosity value, ±0.1 kg.m^{3} for density, and ±0.01 K for temperature. 
Refractive Index was measured with a Leica AR600 refractometer. The temperature was maintained at ±0.01 K with a thermostat Haake D^{8}. 
Vapor pressure data were obtained with a commercial equipment (Labodest, model 602S) available from Fischer Technology, Germany [6]. The observed uncertainties are: ±0.05 K for temperature and ±0.05 kPa for pressure. 
Compositions of both conjugated phases for VLE data were determined by density measurements using the abovementioned equipment. The uncertainty in the mole fractions is estimated to be ±0.001. 
Excess enthalpy measurements at atmospheric pressure and at T=(298 ±1) K of the binary system were performed with an adiabatic calorimeter previously described and using the same methodology [7]. The uncertainty in the mole fraction is estimated to be lower than ±0.001 while, for enthalpy of mixing measurements, it is estimated to be lower than 4% of the measured value. 
FTRaman Spectra: Raman absorption spectra at a resolution of 4 cm^{1} were recorded using a Thermo Scientific , DX Smart Raman Spectrometer both for pure components and the binary mixtures (3.5 M and 6.2 M; in DIPE and chloroform, respectively, where M denotes concentrations in mol.L^{1}). 
FTInfrared Spectra: Infrared spectra of pure components and their binary solutions were recorded on Perkin Elmer FTIR Spectrometer – Spectrum RXI, using a KBr sealed cell with 0.500 mm of path, provided for Perkin Elmer. For this purpose, several highly diluted ternary mixtures were prepared mixing dissimilar masses of both pure, and diluted in CCl_{4}. The concentrations of these ternary solutions were in the range 0.013–0.41 M for both solutes. The infrared spectra of pure components and their binary mixtures were recorded against pure CCl_{4} as a reference in order to eliminate the absorption bands of the solvent. 
Results and Discussion 
Macroscopic properties 
The experimental results of the density, viscosity and refractive index of the pure compounds in the T=(273.15 to 318.15) K range, can be seen in Table 1 together with those of literature [5,811]. 
Experimental results of density, refractive index and viscosity for the [x_{1} chloroform+(1x_{1}) DIPE] binary system in the whole composition range at (298.15 ±0.01) K are listed in Table 2. 
At least three measurements of these properties were made at each temperature and for each solution or pure component. The experimental results reported in Tables 1 and 2 are averages values. 
Polynomial equations were used to correlate the density and viscosity of pure components as a function of temperature and the same properties as a function of mole fraction for the [x_{1} chloroform+(1x_{1}) DIPE] binary system: 
(1) 
(2) 
where a_{i} and b_{i} are parameters, T is the Kelvin temperature and x1 is the mole fraction of chloroform . Table 3 shows the equations obtained for the pure components as a function of temperature. Table 4 shows the same properties for the binary system as a function of composition together with their standard deviations σ. 
Equation 1 and 2 were fitted using a nonlinear regression method based on the LevenbergMarquardt algorithm [12]. with all points equally weighted. 
The excess molar volume, refractive index deviations, viscosity deviation, and excess molar Gibbs energy, were calculated with the following expressions: 
(3) 
(4) 
(5) 
(6) 
where x_{1} is the mole fraction of chloroform, M_{1}, M_{2}, ρ_{1}, ρ_{2}, n_{D1}, n_{D2}, η_{1}, η_{2}, γ_{1} and γ_{2} are the molar masses, densities, refractive indexes, viscosities and activity coefficients the of chloroform and DIPE; respectively, ρ, nD and are the density, refractive index and viscosity of the mixture. The activity coefficients in the liquid phase are related to the vaporliquid equilibrium by: 
(7) 
Here, yi is the mole fraction of component i in the vapor phase, P and are the total pressure and the pure component vapor pressure, respectively, and V_{i} is the molar volume of component i. The value of δ is given by 
(8) 
where Bij the second virial, estimated at T=313.15 K by the Hayden and O'Connell method [13], are: B_{11}=1041.947, B_{22}=1687.182, and B_{12}=1969.413, all expressed in cm^{3} mol^{1}. 
Table 5 shows VLE data for the binary system [x_{1} chloroform+(1x_{1}) DIPE] at T=(313.15 ±0.05) K. These are thermodynamically consistent according to the pointtopoint Van Ness et al. Test [14], as modified by Fredenslund et al. [15], with an average absolute deviation Δy=0.01, and average absolute deviation in pressure ΔP=0.4 kPa, using a treeparameter Legendre polynomial. 
Table 6 shows the excess molar enthalpy H^{E} of the binary system [x_{1} chloroform+(1x_{1}) DIPE] as a function of the mole fraction of chloroform at T=(298 ±1) K. 
Excess thermodynamic functions and deviations from ideality were fitted by means of a RedlichKister type equation1. 
(9) 
where YE represents either V^{E}, Δn_{D}, Δη, H^{E} or G^{E}. The corresponding equations are shown in Table 4, together with the standard deviation σ, defined as: 
(10) 
where N and p are the experimental point and parameter numbers, respectively. 
In all correlation equations the number of constants and the appropriate number of significant digits are selected taking into account the experimental uncertainty reported above for each measured property, obtaining a good agreement with experimental data within experimental error. 
Figures 2 and 3 show the excess molar volume and viscosity deviations plotted against the mole fraction of chloroform at 298.15 K, respectively, while Figure 4 shows the total vapor pressure P, liquidphase x_{1}, and vaporphase y_{1} mole fraction measurements, plotted against the mole fraction of chloroform at 313.15 K. 
Figure 5 and 6 show the excess molar Gibbs energy GE and excess molar enthalpy HE, against x1 at 313.15 K and 298.15 K respectively. As shown G^{E}, H^{E}, and V^{E} are negative over the entire concentration range, with a minimum close to x_{1}=0.5. 
The concentration dependence of these excess properties indicates that the interactions between unlike molecules in the mixture are stronger than between like molecules in the pure liquids and association forces dominate the mixing process. 
It is known that systems with strong interactions between molecules exhibit positive viscosity deviations [1617]. Figure 3 shows that the viscosity deviation is positive over the whole composition range. On the other hand, Δn_{D} close to or greater than 0.004 are indicative of complex formation [1820]. This system presents Δn_{D} values higher than 0.04. Figure 4 shows that this binary system has negative deviations from ideal behavior and Figure 5 shows that this binary system exhibits an exothermic behavior. This behavior was already observed by other researchers [3,4]. 
From this result we can conclude that association is the dominant process in this system. 
Theoretical calculations 
The optimized structures and vibrational frequencies for DIPE, chloroform and their 1:1 complex were calculated by means of density functional theory (DFT) techniques using the B3LYP functional [21,22] combined with the 631G(d,p) basis set. Although MP2 [23] method is preferred to the B3LYP one regarding binding energies, the same tendency can be observed in B3LYP calculations with the advantage of lowering the computational cost [24]. All calculations were made for the isolated species using the Gaussian 03 set of programs [25]. 
For both levels of theory, the geometry was fully optimized with C_{3V} symmetry for chloroform and without symmetry constraints for DIPE and their complex. The frequency calculation confirmed that all the structures are global minimum on the potential surface. 
Theoretical frequencies, infrared intensities and Raman scattering were used to simulate the experimental spectra. The vibrational frequencies of DIPE, chloroform and their complex calculated with the B3LYP/631G(d,p) combination were multiplied by the factor 0.9613 [26] to better reproduce the experimental data. No anharmonicity corrections of these frequencies were made because of lack of necessary experimental data. These calculated values appear in Table 7 along with the frequencies observed in the infrared and Raman spectra. 
Spectroscopic studies 
IR spectra for the pure components and their mixtures in different concentrations, diluted in a solvent "inert" (from the point of view of interaction by hydrogen bonding) and Raman spectra for the pure components and one mixture were performed. 
Figure 7 shows the Raman spectra at T=293 K of pure chloroform and DIPE, together with one binary mixture with a concentration (3.5M and 6.2M in DIPE and chloroform, respectively). It is observed that there are significant differences between the spectra of the pure components and the mixture. 
In the CH stretching region around 3000 cm^{1} the pure chloroform has a very intense band at 3017 cm^{1}, assigned to the CH stretching vibrational mode of free molecules (monomers), which remarkably decreases when this chemical is mixed with DIPE and a new band appears at 2992 cm^{1} which is attributed to the hydrogenbonded species. Pure DIPE presents several bands around 2900 cm^{1}, which are assigned to symmetric and asymmetric stretching modes of the methyl groups. The shape of the band is not significantly affected, but the intensity changes due to dilution 
From the above observations, we conclude that the CH group of the chloroform molecule is directly involved in the formation of a hydrogenbonded complex, as occurs for the system chloroformx butanone [27]. 
In the COC stretching zone are two bands at approximately 750 cm^{1} (strong) and 1000 cm^{1} (very weak) corresponding to the symmetric and antisymmetric stretching respectibely for pure DIPE. The shape of the COC stretching band is significantly modified when DIPE is mixed with chloroform. As we can see the change in shape suggests the presence of new bands underneath it. Using the same methodology previously applied in other studies curve fitting is performed in both bands [27]. The most appropriate set corresponds to two bands in each, whose maximums are located at 791 cm^{1} 794 cm^{1} and 1007 cm 11003 cm^{1}, which were assigned to the symmetric and antisymmetric stretching to the COC the monomer and dimer respectively. 
Infrared spectra were recorded at highly dilution of both components and the different mixtures, using CCl_{4} as solvent. IR spectra show a similar behavior to that observed for Raman. 
When comparing the results of the theoretical calculations with IR and Raman spectra, both frequencies and intensities show obvious similarities. This is conclusive evidence of AB complex formation between unlike molecules in this binary system. Therefore the chloroform (A) and DIPE (B) are in chemical equilibrium with the complex (AB), represented by reactions (11). 
A+B → AB (11) 
The equilibrium constants for these reactions in terms of molar concentrations K_{D}, may be written as 
(12) 
where K_{D} is the association constant, and C_{i} the concentrations of the species present. 
Macroscopic molar concentrations C_{A} is 
(13) 
where is the monomer concentration of chloroform and C_{AB} is the concentration of dimers formed in the mixing process. Using Beer's law A= a_{0}.L.C and (12) in (13) 
(14) 
Where are the absorbances and are the absorptivities of the monomeric species corresponding to the bands under study, and L is the optical path. Absorptivities (of the monomer bands) are obtained from the spectrum of the pure components diluted in CCl_{4} . All spectra are always performed under the same conditions of temperature and pressure. 
Using the same process of curve fitting the CH stretching bands of chloroform and CO stretch of DIPE in the IR spectra of the pure components and mixtures could be calculated absorptivities and the equilibrium constant. 
The absorptivity value for CO stretch of DIPE is (17000 ±1000) dm^{2}.mol^{1} and the CH stretching bands of chloroform is the same result we obtained in a previous work [27]. The K_{D} values at T=293.15 K, evaluated through eq (12), is (2.9 ±0.1) L mol^{1}. 
On the other hand, K_{D}, in term of molar concentration is related to K, in term of mole fractions, by 
(15) 
where P, R, and T are the normal pressure (1 atm), universal gas constant ( 0.08206 L.atm.mol^{1}.K^{1}), and absolute temperature ( 293 K), respectively, while Δn=1 for dimer formation, that leads to: 
The ideal associated solution model was also applied, the simplest case, as this system. The total heat of formation (h_{T, f}) for this model and for this situation (unassociated pure compounds and in the mixture with associations in a 1:1) is obtained from the expression; 
(16) 
The heat of formation (h_{T, f}) values at T=293.15 K, evaluated through equation (16), 
Comparing the thermodynamic results with those calculated by the model of ideal associated solutions from spectroscopic data, we can conclude that, within experimental error, the model is applicable as shown in Figure 5. 
Conclusions 
Empirical equations of density, refractive index and viscosity were deduced, to estimate these properties both for pure and mixed components within the reported standard deviations. Calculated values are in agreement with the experimental results, within the experimental error 
Isothermal vaporliquid equilibrium data of the binary system were obtained at T=313.15 K, showing negative deviations from Raoult's law and a minimum pressure azeotrope. 
From the measured properties, excess functions and deviations from ideality are calculated and it can be seen that all excess properties and deviations from ideality lead to the conclusion that there is an association between different molecules and it can be assumed an 1:1 stoichiometry. 
The IR and Raman spectrum show agreement with the theoretical calculations from the proposed interaction model. There is an association by hydrogen bonding with AB stoichiometric relationship. This conclusion was obtained analyzing the IR and Raman spectrum of chloroform, DIPE, and mixtures together with quantum chemical calculations with DFT procedures. 
Furthermore, the ideal associated solution model is able to represent the experimental excess molar enthalpy, as we can see in Figure 5 that leads to equilibrium constants values well comparable with that calculated from infrared spectroscopy. 
Acknowledgements 
Financial support from the Consejo de Investigaciones de la Universidad Nacional de Tucumán of Argentina, through Grant 26/E516, is gratefully acknowledged. As well as the doctoral fellowship granted by CONICET to R.N.Medina Naessens. 
References 

Table 1  Table 2  Table 3  Table 4 
Table 5  Table 6  Table 7 
Figure 1  Figure 2  Figure 3  Figure 4 
Figure 5  Figure 6  Figure 7  Figure 8 