ACTA CHROMATOGRAPHICA, NO. 17, 2006 CALORIMETRY COMPARED WITH GAS CHROMATOGRAPHY: ARE COMPARABLE MOLAR HEAT CAPACITIES OBTAINED BY USE OF THESE TWO APPROACHES? K. Ciążyńska-Halarewicz, P. Korzenecki, M. Helbin, and T. Kowalska Institute of Chemistry, Silesian University, 9 Szkolna Street, 40-006 Katowice, Poland SUMMARY In this study molar heat capacities (Cp) were determined for aliphatic alcohols by using simple gas chromatographic data and mathematical models correctly derived from the laws and rules of physical chemistry and chromatography. The results obtained were compared with those derived by direct calorimetric measurement and with data calculated by use of the Kopp and Neumann rule. In this way an attempt was made to answer the question posed in the title of this paper–is the accuracy of the non-standard approach (even, if computationally rather complicated) comparable with that of the traditional measurement technique? From the results obtained it is clearly apparent that the results obtained by use of our novel approach are comparable with those originating from the other measurement and computational approaches. This paper is the latest in a series devoted to alternative methods for derivation of thermodynamic data from gas chromatographic results. INTRODUCTION In a previous publication [1] we dealt with estimation of the molar heat capacity (Cp) for three classes of compound, alkylbenzenes, aldehydes, and ketones. Instead of tedious and traditional calorimetric measurement of these values we proposed a novel approach which made use of several mathematical equations. These equations couple gas chromatographic data with physicochemical magnitudes and connect the theory of gas chromatography with the laws and rules of physical chemistry. They enable estimation of Cp values by indirect and relatively complex computation, yet from experiments much simpler than calorimetric determination. The - 72 - only data needed for these calculations are the retention times of the analytes, acquired under strictly defined chromatographic conditions. The starting point for derivation of the ultimate relationships enabling determination of molar heat capacity are given in Table I, where r is the relative retention, rG the non-reduced relative retention, k the retention factor, Vm the molar volume of the analyte, TB the boiling point of the analyte, ∆H° the enthalpy of vaporization of the analyte, recalculated per volume unit of a given congener, Tc the temperature of the column, and R the gas constant. A, B, and C are physicochemically meaningful fitting terms (in Table I we give the values of B and C only, because they alone are thermodynamically important; these terms are identical for all the relationships considered). Table I Basic equations used for derivation of relationships enabling calculation of molar heat capacity, Cp Equation Derivation V r = A exp BVm + C m TB Vm + const. b rG = A exp BVm + C TB a c V ln k = A + BVm + C m TB Fitting terms [2] [3] B= ∆H o RTc C=− ∆H o R [4] In all three equations the chromatographic property (r, rG, or ln k) depends on two physical magnitudes, the molar volume, Vm, and boiling point, TB, of the analyte. Thus the same procedure can be applied to each of these three equations, resulting in the dependences enabling estimation of Cp. Three simple relationships are needed: Vm = M d (1) where M and d are, respectively, the molar weight and density of the compound. - 73 - The molar weight can, however, be expressed as the ratio of the molar heat capacity (or the molar heat), Cp, to the specific heat capacity, c: M= Cp (2) c Combining eqs (1) and (2), we obtain: Vm = Cp (3) dc Now let us make use of the Kopp and Neumann rule. According to this rule the molar heat capacity of a compound in the condensed state is equal to the sum of the atomic heats of the individual atoms that constitute the compound. If we divide the molar heat capacity approximated in this way for a given compound by its molar weight, we obtain the specific heat capacity, which can be inserted in the denominator of eq. (3). The molar heat capacity in the numerator of the same expression is the physical magnitude we are trying to obtain. We should, perhaps, add that the Kopp and Neumann rule cannot be applied to the specific heat capacities of the elements in the sense that we cannot sum them to obtain the specific heat capacities of our analytes, the magnitude of which must be introduced in eq. (3). Table II Final relationships enabling estimation of the molar heat capacities of the chromatographed analytes Equation Fitting terms 1 1 r = A exp B + C dcTB dc 1 1 + const. II rG = A exp B + C dcTB dc I III ln k = A + B B= ∆H oCp RTc C=− ∆H oCp R 1 1 +C dc dcTB Further transformations consist in inserting the complete physical description of the fitting terms B and C and of eq. (3) into the relationships - 74 - given in Table I. The necessary rearrangements lead to the relationships presented in Table II. Here again the fitting terms B and C of each equation have the same physical description, thus the practical procedure for calculating the molar heat capacity will be the same in each case. These physical descriptions are: BRTc ∆H o from term B: Cp = from term C: Cp = − CR ∆H o The relationships listed in Table II are derived elsewhere [1] so the derivations will not be repeated here. It might, however, prove useful to recollect the main assumptions that connect all these models, because some of these assumptions will be referred to below. They focus mainly on the following issues: • The retention process. This consists in partitioning of the analyte between the stationary and mobile phases (in partition chromatography); the non-chromatographic physical analogy of the retention process can be alternate vaporization and condensation; • Transfer of the analyte from the low-polarity stationary phase to the gas phase; in non-chromatographic physical terms this can be described as vaporization of the pure liquid analyte; • Vaporization of a liquid. For the vast majority of liquids the entropy of their vaporization is approximately 85 ± 3 J mol–1 K–1 (Trouton’s Rule); • The processes of vaporization and condensation. These are the mutually opposite processes and according to the Lavoisier–Laplace law their thermodynamic characteristics should be the same absolute values but with opposite signs; • Intermolecular interactions. The assumed models do not take into consideration any intermolecular interactions and, therefore, systems with the weak, and therefore negligible, interactions should fit these models better than those in which the interactions are strong. It is apparent from previous results [1] that agreement between molar heat capacities available in the literature for alkylbenzenes, aldehydes, and ketones and those determined by our approach is fully satisfactory. The situation for the analytes investigated in this work, aliphatic alcohols, - 75 - depends on the results of the chromatographic experiment. Scrutiny of this agreement is the main objective of this study. EXPERIMENTAL As has already been mentioned, the analytes investigated in this study were aliphatic alcohols. n-Alkanes were used as reference analytes enabling determination of Kováts retention indices (PolyScience Corporation, Niles, IL, USA and J.T. Baker, Deventer, The Netherlands; purity grade: standards for GC). The analytes and the respective physical data needed for calculations are listed in Table III. The measurement conditions were the same as in our previous investigations. Experiments were performed with a Fisons Instruments (Rodano/Milan, Italy) GC 8000 series capillary gas chromatograph with oncolumn injection and flame-ionization detection. Isothermal chromatography was performed in the temperature range 323 to 423 K at 25-K intervals (five measurement temperatures in total). Analysis was performed on three different capillary columns (Agilent Technologies, Palo Alto, CA, USA, formerly J&W Scientific, Folsom, CA, USA) each of length 30 m, internal diameter 0.32 mm, and stationary phase film thickness 1 µm. Two low-polarity stationary phases, DB-1 (100% polydimethylsiloxane, polarity on the McReynolds scale 217) and DB-5 (95% polydimethylsiloxane + 5% phenyl, polarity 323) and one mediumpolarity stationary phase, DB-Wax (100% poly(ethylene glycol), polarity 2188) were used. Acquisition and immediate processing of the chromatographic data were achieved by use of the computer software Chromax (Artur Dzieniszewski, Pol-Lab, Warsaw, Poland). The software Statistica 5.1. PL (StatSoft, Tulsa, OK, USA) and Matlab 6.5.0. (The MathWorks, Natick, MA, USA) was used for statistical analysis of the data. Excel (Microsoft Corporation, Redmond, WA, USA) was used for other computations. EXPERIMENTAL AND COMPUTATIONAL PROCEDURE. RESULTS AND DISCUSSION Models I–III presented in this study were tested in the first instance to assess their ability to predict the chromatographic behaviour of the alcohols under the working condition used (an indispensable precondition) but also their usefulness for estimation of molar heat capacities for these - 76 - Table III Physicochemical properties of the aliphatic alcohols used in the gas chromatographic investigations Alcohol Methanol Ethanol 1-Propanol 1-Butanol 1-Pentanol 1-Hexanol 1-Heptanol 1-Octanol 1-Nonanol 1-Decanol 1-Undecanol 1-Dodecanol 2-Propanol 2-Butanol 2-Hexanol 2-Heptanol 2-Octanol 3-Hexanol 2-Methyl-1-propanol 2-Methyl-1-butanol 3-Methyl-1-butanol 2-Methyl-2-butanol 2-Methyl-1-pentanol 4-Methyl-2-pentanol 2-Methyl-3-pentanol 3-Methyl-3-pentanol 2,4-Dimethyl-3-pentanol Empirical formula CH3OH C2H5OH C3H7OH C4H9OH C5H11OH C6H13OH C7H15OH C8H17OH C9H19OH C10H21OH C11H23OH C12H25OH C3H7OH C4H9OH C6H13OH C7H15OH C8H17OH C6H13OH C4H9OH C5H11OH C5H11OH C5H11OH C6H13OH C6H13OH C6H13OH C6H13OH C7H15OH d20 Vm TB c 0.790 0.790 0.800 0.810 0.810 0.820 0.820 0.830 0.830 0.830 0.830 0.830 0.786 0.810 0.810 0.820 0.820 0.819 0.802 0.820 0.810 0.818 0.824 0.802 0.819 0.824 0.829 40.56 58.32 75.12 91.51 108.83 124.61 141.71 156.90 173.81 190.71 207.60 224.51 76.46 91.51 126.15 141.71 158.82 124.76 92.42 107.50 108.83 107.76 124.00 127.41 124.76 124.00 140.17 337.65 351.45 370.40 391.15 411.15 430.15 449.15 466.15 488.15 504.65 521.50 534.65 355.55 372.65 410.15 433.15 452.15 408.15 381.15 401.15 404.65 384.65 421.15 405.15 401.15 396.00 412.65 2.524 2.566 2.589 2.603 2.613 2.620 2.625 2.629 2.632 2.635 2.638 2.639 2.589 2.603 2.620 2.625 2.629 2.620 2.603 2.613 2.613 2.613 2.620 2.620 2.620 2.620 2.625 d20 is the density at 20°C (g cm–3), Vm is the molar volume (cm3 mol–1), TB is the boiling point (K), and c is the specific heat capacity (J g–1 K–1) compounds. The statistical correctness of the models, when used to predict retention data for the alcohols, enabled further acquisition of the thermodynamic data. As is clearly apparent from models I–III, to determine the molar heat capacity, Cp, knowledge of another thermodynamic magnitude is necessary, i.e. the enthalpy of vaporization of the compounds (more precisely, in each of these models we encounter the magnitude ∆H°, which must be multi- 77 - plied by the molar volume of the analyte to give its enthalpy of vaporization, ∆Hvap). Achieving the ultimate objective, the molar heat capacity, Cp, is rather complicated; let us present and comment on the entire data-acquisition procedure in detail. 1. Acquisition of the experimental data This stage is purely experimental and devoted to collecting the chromatographic data (i.e. the retention times of the analytes) under the working conditions given in the Experimental section. Here it must be stated that the group of 27 normal and branched aliphatic alcohols selected was in a sense a ‘risk’ group, because the models assumed are supposed to perform best for low-polarity analytes and some of the alcohols employed in this study (those the lowest molecular weight) are quite polar. Chromatography was performed at five different temperatures under isothermal conditions, using three stationary phases of different polarity in columns with identical dimensions, to avoid any effect of these on the results. The film thickness of the stationary phases (1 µm) on all the columns was relatively high, to prevent intermolecular interactions between the analytes and the capillary walls. To minimize intermolecular interactions between the analytes and the stationary phase, we selected two different low-polarity stationary phases and one medium-polarity stationary phase, deliberately avoiding stationary phases of high polarity. The retention times of the n-alkanes, needed to calculate the Kováts retention indices used in some equations, were determined analogously. 2. Experimental basis The retention data needed, i.e. the relative retention (r), the nonreduced relative retention (rG), the retention factor (k) (and the Kováts retention index, I, needed for determination of enthalpy of vaporization only) were derived from the chromatographic data collected. In that way a large database emerged. Determination of the numerical values of all these quantities was required by the retention models tested in this study. 3. Computational basis – search for statistically valid datasets Statistically valid datasets were built separately for each group of results collected under the given measurement conditions. This was done with aid of the LMS (least median of squares) procedure. The LMS method – very useful for elimination of results burdened with substantial ex- 78 - perimental error (i.e. outliers) and for preliminary determination of fitting and statistical data – can only be applied to linear equations, so non-linear equations must first be linearized. This applies to equations from which both molar heat capacity and the enthalpy of vaporization were derived. 4. Appropriate computational procedures Appropriate computational procedures performed with both sets of the equations enabled the fitting terms (A, B, and C) and statistical data characterizing the tested datasets to be obtained. For the linear equations the linear regression procedure was used; for the exponential equations the procedure of non-linear estimation (i.e. the quasi-Newton method) was used. 5. Statistical evaluation of the models After having constructed the computational basis (point 3) and performed the computational procedure appropriate for the type of the equation (point 4), the errors describing the quality of fits of the models to the experimental data (RMS, root-mean-square error) were determined and the predictive power of the models was assessed (by means of RMSCV, the root mean square of cross-validation). Last, but not least, the correlation coefficients, r, were also determined. This evaluation was performed for each equation (both those used to obtain the molar heat capacities and those used to determine the enthalpy of vaporization) according to the procedures described for each of the five working temperatures and for each of the three stationary phases. The retention data obtained differed by as much as three orders of magnitude, which might to some extent affect RMS and RMSCV values. Thus to compare the magnitudes of the errors the respective percentage values (%RMS, %RMSCV) had to be calculated. Lower values of the percentage errors indicated that the given model performs better for one dataset than for another (i.e. that for which error values were higher). The errors very seldom exceeded a few percent, and slightly higher values were obtained only very sporadically. Comparison of all the data revealed the best results (i.e. those with the lowest percentage errors) originated from the most polar (DB-Wax) column. The correlation coefficients (r) obtained were also very satisfactory, usually in excess of 0.99 (the lowest value was 0.94). Examples of results from statistical analysis are shown in Table IV. - 79 - Table IV Statistical evaluation of mathematical models I–III and of the values of the fitting terms B and C, with the respective errors (n is the number of valid cases) Eq. I II III Stationary phase DB-1 DB-5 DB-Wax DB-1 DB-5 DB-Wax DB-1 DB-5 DB-Wax Fitting terms B ± sB 78.13 ± 0.53 16.27 ± 0.35 80.04 ± 0.95 72.28 ± 1.39 11.69 ± 0.28 74.34 ± 5.61 37 ± 1.7 68 ± 0.6 14 ± 0.8 C ± sC −18263 ± 18 −15917 ± 24 −15285 ± 14 −16139 ± 21 −14049 ± 29 −13533 ± 28 −13551 ± 18 −16645 ± 52 −95236 ± 85 n 36 45 48 42 57 42 45 48 33 RMS % RMS 2.67 7.54 1.14 0.08 0.17 0.19 0.22 0.08 0.11 6.0 7.6 3.3 4.1 6.5 3.4 8.4 5.4 4.0 % RMSCV RMSCV 2.77 7.82 1.64 0.10 0.18 0.41 0.23 0.08 0.12 6.2 7.9 3.8 4.6 6.8 4.1 8.6 5.7 4.2 r 0.9987 0.9988 0.9998 0.9978 0.9982 0.9898 0.9812 0.9986 0.9938 In eqs. (I), (II), and (III) Tc = 348, 373, and 323 K, respectively To better illustrate the statistical data, predicted values were plotted against experimental results (Fig. 1) for one – although somewhat extended – example taken from Table IV. These plots show that agreement between experimental data and those predicted by regression analysis was excellent. The high quality statistical data confirmed that the indispensable precondition was fulfilled by the equations considered, which meant they could now be used for prediction of the retention of the test analytes in the chromatographic systems selected. Let us therefore proceed to the next stage of the computations, determination of the thermodynamic data. 6. Determination of thermodynamic data and their evaluation by comparison with those taken from the literature As already mentioned, the molar heat capacities (Cp) must be derived in two steps. First, numerical values of the enthalpy of vaporization, ∆Hvap, must be known for the compounds. These can be acquired by transformation of the fitting terms A, B, or C of eqs (a)–(i) given in Table V (the physical significance of these terms can be found elsewhere [2,5,6]). In models I–III there is a term related to the enthalpy of vaporization (∆H°); we therefore need to know the molar volumes of all the analytes investigated. When ∆H° was first introduced [4] it was defined as the enthalpy of vaporization calculated per the unit volume of a compound belonging to a given group of congeners (eq. 4), because of the method used for acquisition of the thermodynamic data: - 80 - ∆H vap = ∆H °×Vm (4) The magnitude of ∆H° and such physical constants as the analytes’ boiling points (TB), densities (d), and specific heat capacities (c) (all from reference sources) plus the gas constant (R) and the temperature of analysis (Tc) were used to determine the respective molar heat capacities from the regressionally fitted terms B or C of eqs (I)–(III). Observed 30 26 22 18 14 10 6 (a) 2 -2 0 4 8 12 16 20 Predicted Observed 30 26 22 18 14 10 6 (b) 2 -2 -2 2 6 10 14 18 22 26 30 Predicted 140 Observed 120 100 80 60 40 20 (c) 0 -20 -20 0 20 40 60 80 100 120 140 Predicted Fig. 1 Graphical evaluation of the predictive power of the models used. Plots of observed (experimental) values against predicted values (eq. (I), stationary phase: DB-Wax, dependent variable r) (a) Tc = 323 K, (b) Tc = 373 K, (c) Tc = 423 K - 81 - Table V Mathematical models used to determine enthalpies of vaporization (the fitting terms are described in detail elsewhere [2,5,6]) Model a V r = A exp BVm + C m TB b V rG = A exp BVm + C m + const. TB c V ln k = A + BVm + C m TB d R r = A exp BRm + C m TB e f R rG = A exp BRm + C m + const. TB R ln k = A + BRm + C m TB g rG = A exp(BTB ) + const. h I = A + BTB i I = A+ B TB The two-step determination of the thermodynamic magnitude has already been practised (eqs (h) and (i) [5,6]) for determination of the enthalpy of vaporization of the analytes (step 2) with the aid of the chemical potential of partitioning of one methylene group between the stationary and mobile phases (∆µp(-CH2-); step 1). It has been our previous experience that in the course of the two-step procedure overall determination error tends to increase and so we might now expect a similar effect of a marked error of determination with molar heat capacities also. Let us return to our current subject, however. The experimental enthalpy of vaporization values, i.e. those obtained for the test analytes by use - 82 - of eqs (a)–(i), are in very good agreement with data taken from the literature. The best results were obtained for DB-Wax, the most polar of the three stationary phases used. Most of the results from this stationary phase agree with the reference data within 10%. The results obtained from the low-polarity stationary phases were also very satisfactory and were easily within 25% of the reference data. (This issue will be discussed in more depth elsewhere [7].) To summarise, we obtained a vast number of enthalpy of vaporization data, ∆Hvap, from the nine different equations, eight of which had two thermodynamically relevant fitting terms (eqs (a)–(f), (h), and (i)) and eq. (g) had one such fitting term. All of these numerical values could be inserted in eqs (I)–(III) to derive the Cp values from either B or C. We tested all possible combinations and our assumptions about the quality of the results obtained in the two-step procedure proved true. In the other words, despite the abundance of the ∆Hvap data obtained, those furnishing acceptable values of molar heat capacity (Cp) proved much less numerous. These circumstances give rise the question: Would it not be more practical to use reference values of the enthalpy of vaporization rather than those calculated from our experiment? The answer is, however, no. First, reference data originating from direct (i.e. calorimetric) measurements are few and those originating from the other computational approaches are also burdened with approximation errors. Second, no criteria are available enabling assessment of the quality of thermodynamic results originating from the literature that might help us distinguish between ‘good’ and ‘wrong’ results. One argument in favour of our approach is that values for both enthalpy of vaporization and molar heat capacity were obtained from relationships that underwent an analogous statistical–computational procedure and were derived from the same database. Last, but not least, it is not only the two-step determination procedure which must be blamed, but also the many simplifications and approximations contained in the models given by eqs (a)–(i) and (I)–(III). There is yet one more reason for the moderate success of our comparisons – the availability of reference data. Molar heat capacity values for primary alcohols are relatively readily available in the literature whereas for the secondary and the ternary alcohols most of the reference data were calculated by use of the highly approximate Kopp and Neumann rule (Table VI). This procedure consists in adding the atomic heats for all the atoms of the elements present in the structure of a given compound. In our experiment we were making use of the atomic heats for three elements only, car- 83 - bon (graphite; 8.53 J mol–1 K–1), hydrogen (14.415 J mol–1 K–1), and oxygen (14.68 J mol–1 K–1) [12]. Table VI Literature values of the molar heat capacities (Cp, kJ mol−1) of the alcohols Numerical values of Cp taken from literature [kJ mol–1] Alcohol a) b) c) d) e) f) g) h) i) Methanol Ethanol 1-Propanol 1-Butanol 1-Pentanol 1-Hexanol 1-Heptanol 1-Octanol 1-Nonanol 1-Decanol 1-Undecanol 1-Dodecanol 2-Propanol 2-Butanol 2-Hexanol 2-Heptanol 2-Octanol 3-Hexanol 2-Methyl-1-propanol 2-Methyl-1-butanol 3-Methyl-1-butanol 2-Methyl-2-butanol 2-Methyl-1-pentanol 4-Methyl-1-pentanol 2-Methyl-3-pentanol 3-Methyl-3-pentanol 2,4-Dimethyl-3-pentanol 81.06 112.41 143.64 177.15 208.03 240.49 272.49 304.74 336.48 368.47 400.47 432.46 156.26 81.92 115.90 146.88 177.70 207.40 241.32 270.80 305.55 335.56 367.14 398.72 430.31 85.73 121.79 157.38 194.13 227.21 262.90 295.40 329.50 363.70 398.50 441.00 467.90 93.71 137.66 179.50 219.46 255.45 295.20 332.40 368.20 404.70 436.50 481.00 507.00 100.22 149.95 195.40 235.40 272.87 313.00 353.00 390.00 427.70 458.90 502.00 529.60 112.00 179.20 218.00 255.40 293.00 333.00 373.00 410.00 448.20 483.10 526.00 554.30 121.80 191.71 230.40 265.70 303.00 342.00 381.00 417.00 454.20 491.80 536.00 562.80 73.15 106.37 139.59 172.81 206.03 240.10 272.10 305.20 334.40 370.60 407.20 438.30 80.87 118.23 155.59 192.95 230.31 267.67 305.03 342.39 379.75 417.11 454.47 491.83 155.59 192.95 267.67 305.03 342.39 267.67 192.95 230.31 230.31 230.31 267.67 267.67 267.67 267.67 305.03 198.03 181.59 181.05 247.30 247.63 293.38 a) measurement at the boiling point, p = 1 atm [8]; b) [9]; c) 320 K [10]; d) 350 K [10]; e) 370 K [10]; f) 400 K [10]; g) 420 K [10]; h) 298.15 K, HP-5 [11]; i) Kopp and Neumann rule Comparison of the results obtained in this study with those from the literature and with other reference data (calculated by use of the Kopp and Neumann rule) enables evaluation of the models tested and, consequently, evaluation of the sensitivity and correctness of our approach. Let us discuss our final results. Here we present a small fraction of the data only, because of their gargantuan abundance (all the results obtained are available elsewhere [13]). The data shown in Tables VII–IX alone, however, enable assessment of the quality of the thermodynamic results. It - 84 - is readily apparent that a substantial number of these results can be described as very good or good, i.e. deviating by ±10 or ±25%, respectively, from the reference values (i.e. from those taken from literature or calculated by use of the Kopp and Neumann rule). These data are indicated in bold characters. Another batch of the results can be described as satisfactory, because their deviation from the reference data is within the range ±50% (these data are written in italics). From eqs (a)–(c) one can obtain average values of the enthalpy of vaporization only, hence single molar heat capacity values also, for the whole set of the congeners chromatographed under the given working conditions. Another insight into the data is provided in Fig. 2. From this figure it is apparent that use of fitting term B from eq. (I) in combination with enthalpy of vaporization values derived from eq. (i) (term A) results in unexpectedly good prediction of the molar heat capacity values, despite the complexity of this two-step approach. The vast majority of the results obtained deviate from the reference data by ±50% and at the lowest measurement temperature the deviation is even less (±25% or ±10% only). The dependence of the quality of predicted molar heat capacity values on stationary phase polarity is readily apparent from Fig. 3. This comparison was made for analyses run at the same temperature and for computations with the same equations. From the data in Fig. 3a it can be concluded that the performance of the mathematical models tested does not depend on the type of stationary phase used. For each stationary phase the number of results falling outside the range of ±50% is very low and most of the results fall well within ±25% or ±10%. To summarise, it can be stated that the best results (i.e. the molar heat capacity values closest to the reference data) were obtained from use of the medium-polarity DB-Wax stationary phase. This may seem somewhat striking, because one of the basic assumptions made when deriving the equations of the model was that the models – taking no account of any intermolecular interactions within the chromatographic system – should perform best for systems lacking specific intermolecular interactions. How can we explain this seeming contradiction? Among the model assumptions we find another that refers to transfer of the analyte from the liquid stationary phase to the gaseous mobile phase – it was assumed that this transfer resembles vaporization of a liquid from its bulk phase. Because our test analytes are alcohols and the most polar stationary phase applied is 100% poly(ethylene glycol) (and thus the analytes and stationary phase have hydroxyl groups in common), the assumption that the physical natu- 85 - Table VII Molar heat capacity values (Cp) calculated from the fitting terms B or C of eq. (III) by use of experimental gas chromatographic data (DB-Wax, Tc = 323K) and implementation of the enthalpy of vaporization values derived from eqs (a)–(i) (terms A, B, or C) Analyte Term, eq. (III) Tc = 323K Methanol Ethanol 1-Propanol 1-Butanol 1-Pentanol 1-Hexanol 1-Heptanol 1-Octanol 1-Nonanol 1-Decanol 1-Undecanol 1-Dodecanol 2-Propanol 2-Butanol 2-Hexanol 2-Heptanol 2-Octanol 3-Hexanol 2-Methyl-1-propanol 2-Methyl-1-butanol 3-Methyl-1-butanol 2-Methyl-2-butanol 2-Methyl-1-pentanol 4-Methyl-1-pentanol 2-Methyl-3-pentanol 3-Methyl-3-pentanol 2,4-Dimethyl-3-pentanol a/B B a/C C 87.82 189.01 B C 72.84 156.76 Cp values (kJ mol–1) calculated by use of the fitting terms of eq: b/B b/C c/B c/C B C B C B C B C 79.19 170.43 66.02 142.10 - 86 - 86.16 185.44 73.79 158.81 d/B B 96.45 87.49 83.00 80.17 78.32 77.02 76.02 75.26 74.60 74.04 73.63 73.63 84.37 80.58 77.74 76.52 76.14 77.49 80.85 78.20 79.12 79.17 77.00 78.15 77.16 77.00 75.89 d/C C 207.57 188.29 178.62 172.54 168.56 165.76 163.60 161.97 160.56 159.35 158.46 158.46 181.58 173.42 167.31 164.69 163.87 166.77 174.01 168.31 170.28 170.39 165.72 168.20 166.07 165.72 163.33 B 80.16 72.71 68.98 66.63 65.09 64.01 63.18 62.54 62.00 61.53 61.19 61.19 70.12 66.97 64.61 63.60 63.28 64.40 67.19 64.99 65.76 65.80 64.00 64.95 64.13 64.00 63.07 C 172.51 156.49 148.45 143.39 140.09 137.76 135.96 134.61 133.44 132.43 131.70 131.70 150.90 144.13 139.04 136.87 136.19 138.60 144.61 139.88 141.52 141.61 137.73 139.79 138.02 137.73 135.74 Analyte Term, eq. (III) Tc = 323K Methanol Ethanol 1-Propanol 1-Butanol 1-Pentanol 1-Hexanol 1-Heptanol 1-Octanol 1-Nonanol 1-Decanol 1-Undecanol 1-Dodecanol 2-Propanol 2-Butanol 2-Hexanol 2-Heptanol 2-Octanol 3-Hexanol 2-Methyl-1-propanol 2-Methyl-1-butanol 3-Methyl-1-butanol 2-Methyl-2-butanol 2-Methyl-1-pentanol 4-Methyl-1-pentanol 2-Methyl-3-pentanol 3-Methyl-3-pentanol 2,4-Dimethyl-3-pentanol e/B e/C B C B C 96.31 87.37 82.88 80.06 78.21 76.91 75.91 75.15 74.50 73.93 73.53 73.53 84.25 80.47 77.63 76.41 76.03 77.38 80.74 78.09 79.01 79.06 76.89 78.04 77.06 76.89 75.78 207.28 188.03 178.37 172.30 168.32 165.52 163.37 161.74 160.34 159.12 158.24 158.24 181.32 173.18 167.07 164.45 163.64 166.54 173.76 168.07 170.04 170.15 165.49 167.96 165.84 165.49 163.10 80.08 72.65 68.91 66.57 65.03 63.95 63.12 62.49 61.95 61.48 61.14 61.14 70.05 66.91 64.55 63.54 63.22 64.34 67.13 64.94 65.70 65.74 63.94 64.89 64.07 63.94 63.01 172.35 156.35 148.32 143.27 139.96 137.63 135.84 134.48 133.32 132.31 131.58 131.58 150.77 144.00 138.92 136.75 136.07 138.48 144.48 139.75 141.39 141.48 137.61 139.66 137.89 137.61 135.62 Cp values (kJ mol–1) calculated by use of the fitting terms of eq: f/B f/C g/B h/A h/B B C B C B C B C B C 105.90 96.06 91.13 88.02 86.00 84.57 83.46 82.63 81.91 81.29 80.84 80.84 92.63 88.48 85.35 84.02 83.60 85.08 88.77 85.87 86.87 86.93 84.55 85.81 84.73 84.55 83.33 227.91 206.74 196.12 189.45 185.08 182.00 179.63 177.83 176.29 174.96 173.99 173.99 199.37 190.42 183.70 180.82 179.93 183.11 191.05 184.80 186.97 187.09 181.96 184.68 182.34 181.96 179.33 91.23 82.76 78.51 75.83 74.09 72.85 71.90 71.19 70.57 70.03 69.65 69.65 79.81 76.22 73.53 72.38 72.02 73.30 76.48 73.97 74.84 74.89 72.84 73.93 72.99 72.84 71.79 196.35 178.11 168.96 163.21 159.44 156.79 154.75 153.20 151.88 150.73 149.89 149.89 171.75 164.04 158.26 155.78 155.01 157.75 164.59 159.21 161.07 161.18 156.76 159.10 157.09 156.76 154.50 34.29 47.36 57.88 66.77 75.55 82.68 90.05 96.07 101.63 107.86 113.62 119.85 61.38 70.09 87.79 93.38 100.25 87.24 69.21 76.49 76.76 79.96 84.04 89.76 88.77 89.37 96.95 73.79 101.93 124.58 143.71 162.59 177.95 193.81 206.75 218.72 232.13 244.53 257.94 132.10 150.84 188.93 200.96 215.76 187.76 148.95 164.61 165.21 172.09 180.86 193.17 191.04 192.35 208.66 23.52 31.80 38.22 44.27 50.61 68.45 82.25 95.29 34.52 47.69 58.66 69.07 74.29 102.64 126.25 148.65 41.76 89.87 62.58 134.68 47.35 101.90 72.27 155.55 58.96 126.89 92.96 200.07 62.16 47.39 51.80 52.41 58.81 58.82 62.38 60.12 60.90 68.18 133.79 101.99 111.48 112.79 126.57 126.60 134.25 129.40 131.06 146.73 96.06 72.56 80.91 81.87 88.16 92.32 96.94 93.91 94.54 106.17 206.73 156.15 174.12 176.21 189.73 198.69 208.64 202.12 203.47 228.49 i/A i/B B C B C 30.27 43.52 56.05 68.28 81.21 92.98 105.74 117.08 129.70 142.31 154.91 167.53 57.05 68.28 94.13 105.74 118.51 93.10 68.96 80.22 81.21 80.41 92.53 95.07 93.10 92.53 104.59 65.14 93.66 120.64 146.96 174.78 200.12 227.58 251.97 279.13 306.27 333.39 360.55 122.79 146.96 202.59 227.58 255.06 200.36 148.42 172.64 174.78 173.06 199.14 204.61 200.36 199.14 225.11 42.80 61.54 79.26 96.56 114.83 131.48 149.53 165.56 183.40 201.23 219.05 236.90 80.68 96.56 133.11 149.53 167.58 131.64 97.52 113.43 114.83 113.70 130.84 134.44 131.64 130.84 147.90 92.11 132.44 170.59 207.81 247.14 282.98 321.81 356.31 394.71 433.08 471.44 509.84 173.63 207.81 286.47 321.81 360.67 283.32 209.88 244.12 247.14 244.71 281.59 289.34 283.32 281.59 318.31 a/B means the enthalpy of vaporization was determined from eq. (a) using term B; values in bold are experimental values within ±10 and ±25% of values taken from the literature; values in italics are the experimental values within ±50% of values taken from the literature - 87 - Table VIII Molar heat capacity values (Cp) calculated from the fitting terms B or C of eq. (III) by use of experimental gas chromatographic data (DB-5, Tc = 373K) and implementation of the enthalpy of vaporization values derived from eqs (a)–(i) (terms A, B, or C) Analyte Term, eq. (III) Tc = 323K Methanol Ethanol 1-Propanol 1-Butanol 1-Pentanol 1-Hexanol 1-Heptanol 1-Octanol 1-Nonanol 1-Decanol 1-Undecanol 1-Dodecanol 2-Propanol 2-Butanol 2-Hexanol 2-Heptanol 2-Octanol 3-Hexanol 2-Methyl-1-propanol 2-Methyl-1-butanol 3-Methyl-1-butanol 2-Methyl-2-butanol 2-Methyl-1-pentanol 4-Methyl-1-pentanol 2-Methyl-3-pentanol 3-Methyl-3-pentanol 2,4-Dimethyl-3-pentanol a/B B a/C C B C Cp values (kJ mol–1) calculated by use of the fitting terms of eq: b/B b/C c/B c/C B C B C B C B C d/B B 201.73 182.99 173.59 167.68 163.81 161.09 158.99 157.40 156.04 154.86 154.00 154.00 176.46 150.29 321.28 135.49 289.64 150.81 322.40 136.24 291.26 171.50 366.63 197.30 421.77 168.54 162.59 160.05 159.25 162.08 169.11 163.57 165.49 165.59 161.06 163.46 161.39 161.06 158.73 - 88 - d/C C B C 431.24 391.19 371.10 358.46 350.20 344.37 339.89 336.49 333.58 331.05 329.22 329.22 377.23 360.30 347.59 342.15 340.45 346.48 361.51 349.67 353.77 354.00 344.30 349.45 345.02 344.30 339.33 183.42 166.38 157.84 152.46 148.95 146.47 144.56 143.12 141.88 140.80 140.02 140.02 160.45 153.24 147.84 145.52 144.80 147.37 153.76 148.72 150.47 150.57 146.44 148.63 146.75 146.44 144.32 392.10 355.68 337.42 325.93 318.41 313.12 309.04 305.95 303.30 301.00 299.34 299.34 342.99 327.60 316.04 311.09 309.55 315.04 328.70 317.93 321.66 321.87 313.05 317.73 313.71 313.05 308.53 Analyte Term, eq. (III) Tc = 323K Methanol Ethanol 1-Propanol 1-Butanol 1-Pentanol 1-Hexanol 1-Heptanol 1-Octanol 1-Nonanol 1-Decanol 1-Undecanol 1-Dodecanol 2-Propanol 2-Butanol 2-Hexanol 2-Heptanol 2-Octanol 3-Hexanol 2-Methyl-1-propanol 2-Methyl-1-butanol 3-Methyl-1-butanol 2-Methyl-2-butanol 2-Methyl-1-pentanol 4-Methyl-1-pentanol 2-Methyl-3-pentanol 3-Methyl-3-pentanol 2,4-Dimethyl-3-pentanol e/B e/C B C B C 203.63 184.72 175.23 169.27 165.36 162.61 160.49 158.89 157.52 156.32 155.46 155.46 178.13 170.13 164.13 161.56 160.76 163.61 170.70 165.11 167.05 167.16 162.58 165.01 162.92 162.58 160.23 435.32 394.89 374.61 361.85 353.51 347.63 343.10 339.67 336.73 334.18 332.33 332.33 380.80 363.70 350.87 345.38 343.67 349.76 364.92 352.98 357.12 357.35 347.55 352.75 348.29 347.55 342.53 186.80 169.45 160.74 155.27 151.69 149.17 147.22 145.75 144.49 143.40 142.60 142.60 163.40 156.07 150.56 148.20 147.47 150.08 156.59 151.46 153.24 153.34 149.14 151.37 149.45 149.14 146.98 399.33 362.24 343.63 331.93 324.28 318.88 314.73 311.59 308.89 306.55 304.85 304.85 349.31 333.63 321.86 316.82 315.25 320.84 334.75 323.79 327.59 327.80 318.82 323.58 319.49 318.82 314.21 Cp values (kJ mol–1) calculated by use of the fitting terms of eq: f/B f/C g/B h/A h/B B C B C B C B C B C 229.45 208.14 197.45 190.72 186.33 183.23 180.84 179.04 177.48 176.14 175.16 175.16 200.71 191.70 184.94 182.04 181.14 184.35 192.35 186.05 188.23 188.35 183.19 185.93 183.58 183.19 180.54 490.51 444.95 422.10 407.72 398.32 391.70 386.60 382.73 379.42 376.55 374.46 374.46 429.07 409.81 395.36 389.17 387.24 394.10 411.19 397.72 402.39 402.65 391.62 397.47 392.44 391.62 385.96 265.19 240.56 228.21 220.43 215.35 211.77 209.01 206.92 205.13 203.58 202.45 202.45 231.98 221.56 213.75 210.40 209.36 213.07 222.31 215.03 217.55 217.69 211.73 214.89 212.17 211.73 208.67 566.91 514.26 487.85 471.24 460.37 452.71 446.82 442.35 438.52 435.20 432.79 432.79 495.91 473.65 456.94 449.79 447.56 455.49 475.24 459.68 465.07 465.37 452.62 459.38 453.57 452.62 446.08 44.62 61.63 75.34 86.91 98.33 107.59 117.18 125.02 132.25 140.36 147.85 155.98 79.88 91.20 114.23 121.52 130.47 113.53 90.06 99.53 99.89 104.05 109.36 116.81 115.52 116.31 126.16 95.38 131.75 161.05 185.78 210.20 230.00 250.50 267.27 282.71 300.06 316.07 333.44 170.76 194.97 244.20 259.78 278.92 242.71 192.53 212.78 213.55 222.44 233.78 249.72 246.94 248.64 269.70 41.32 55.56 67.61 78.11 88.65 96.56 105.10 111.54 118.57 125.48 131.18 88.34 118.77 144.54 166.97 189.50 206.42 224.68 238.45 253.47 268.26 280.43 43.97 60.13 74.19 86.76 99.49 109.54 120.29 128.74 137.95 146.99 154.80 94.00 128.55 158.60 185.48 212.68 234.16 257.15 275.21 294.90 314.23 330.92 70.73 80.35 100.86 108.27 116.17 96.09 151.20 171.78 215.61 231.46 248.34 205.41 77.10 88.66 113.64 123.17 133.33 109.12 164.82 189.54 242.94 263.31 285.02 233.28 88.37 84.27 92.68 97.57 103.76 100.89 101.02 110.19 188.91 180.15 198.13 208.58 221.81 215.67 215.96 235.55 98.99 95.60 102.77 110.32 116.47 113.40 113.38 124.60 211.61 204.38 219.71 235.85 248.99 242.43 242.38 266.37 i/A i/B B C B C 34.94 50.24 64.71 78.82 93.74 107.34 122.07 135.15 149.72 164.27 178.82 193.39 65.86 78.82 108.66 122.07 136.80 107.47 79.61 92.60 93.74 92.82 106.81 109.75 107.47 106.81 120.74 74.69 107.39 138.33 168.51 200.40 229.46 260.95 288.92 320.06 351.18 382.28 413.42 140.79 168.51 232.29 260.95 292.45 229.74 170.18 197.95 200.40 198.43 228.34 234.61 229.74 228.34 258.11 43.95 63.19 81.39 99.15 117.92 135.01 153.54 170.00 188.32 206.63 224.93 243.26 82.84 99.15 136.68 153.54 172.08 135.18 100.14 116.48 117.92 116.76 134.35 138.05 135.18 134.35 151.87 93.95 135.08 174.00 211.96 252.08 288.63 328.24 363.42 402.59 441.73 480.85 520.02 177.10 211.96 292.20 328.24 367.87 288.98 214.07 249.00 252.08 249.60 287.22 295.11 288.98 287.22 324.67 a/B means the enthalpy of vaporization was determined from eq. (a) using term B; values in bold are experimental values within ±10 and ±25% of values taken from the literature; values in italics are the experimental values within ±50% of values taken from the literature - 89 - Table IX Molar heat capacity values (Cp) calculated from the fitting terms B or C of eq. (I) by use of experimental gas chromatographic data (DB-1, Tc = 373K) and implementation of the enthalpy of vaporization values derived from eqs (a)–(i) (terms A, B, or C) Analyte Term, eq. (III) Tc = 323K Methanol Ethanol 1-Propanol 1-Butanol 1-Pentanol 1-Hexanol 1-Heptanol 1-Octanol 1-Nonanol 1-Decanol 1-Undecanol 1-Dodecanol 2-Propanol 2-Butanol 2-Hexanol 2-Heptanol 2-Octanol 3-Hexanol 2-Methyl-1-propanol 2-Methyl-1-butanol 3-Methyl-1-butanol 2-Methyl-2-butanol 2-Methyl-1-pentanol 4-Methyl-1-pentanol 2-Methyl-3-pentanol 3-Methyl-3-pentanol 2,4-Dimethyl-3-pentanol a/B B a/C C B C Cp values (kJ mol–1) calculated by use of the fitting terms of eq: b/B b/C c/B c/C B C B C B C B C d/B B 325.31 295.09 279.94 270.40 264.17 259.78 256.39 253.83 251.63 249.73 248.34 248.34 284.56 239.46 440.06 240.55 442.06 238.93 439.09 244.34 449.03 233.04 428.26 280.99 516.37 271.79 262.20 258.10 256.82 261.37 272.70 263.77 266.87 267.04 259.72 263.60 260.27 259.72 255.97 - 90 - d/C C B C 597.82 542.29 514.45 496.92 485.47 477.39 471.17 466.47 462.43 458.92 456.38 456.38 522.95 499.47 481.85 474.31 471.95 480.32 501.15 484.74 490.42 490.74 477.29 484.43 478.30 477.29 470.40 339.37 307.85 292.04 282.09 275.59 271.00 267.47 264.80 262.51 260.52 259.08 259.08 296.86 283.54 273.54 269.25 267.92 272.67 284.49 275.17 278.40 278.58 270.95 275.00 271.52 270.95 267.03 623.66 565.73 536.68 518.40 506.45 498.03 491.54 486.63 482.42 478.76 476.11 476.11 545.55 521.06 502.68 494.81 492.35 501.08 522.81 505.69 511.62 511.96 497.92 505.37 498.97 497.92 490.73 Analyte Term, eq. (III) Tc = 323K Methanol Ethanol 1-Propanol 1-Butanol 1-Pentanol 1-Hexanol 1-Heptanol 1-Octanol 1-Nonanol 1-Decanol 1-Undecanol 1-Dodecanol 2-Propanol 2-Butanol 2-Hexanol 2-Heptanol 2-Octanol 3-Hexanol 2-Methyl-1-propanol 2-Methyl-1-butanol 3-Methyl-1-butanol 2-Methyl-2-butanol 2-Methyl-1-pentanol 4-Methyl-1-pentanol 2-Methyl-3-pentanol 3-Methyl-3-pentanol 2,4-Dimethyl-3-pentanol e/B e/C B C B C 325.60 295.36 280.19 270.65 264.41 260.01 256.62 254.06 251.86 249.95 248.57 248.57 284.82 272.04 262.44 258.33 257.05 261.61 272.95 264.01 267.11 267.28 259.96 263.84 260.50 259.96 256.20 598.36 542.78 514.91 497.37 485.91 477.82 471.60 466.89 462.85 459.34 456.80 456.80 523.42 499.92 482.29 474.74 472.38 480.75 501.60 485.18 490.87 491.19 477.72 484.87 478.73 477.72 470.82 347.20 314.95 298.77 288.60 281.94 277.26 273.64 270.91 268.56 266.53 265.05 265.05 303.71 290.08 279.84 275.46 274.10 278.96 291.05 281.52 284.82 285.01 277.20 281.34 277.78 277.20 273.19 638.05 578.78 549.06 530.36 518.13 509.51 502.88 497.86 493.54 489.80 487.09 487.09 558.13 533.08 514.27 506.22 503.71 512.64 534.87 517.35 523.42 523.76 509.41 517.02 510.48 509.41 502.05 Cp values (kJ mol–1) calculated by use of the fitting terms of eq: f/B f/C g/B h/A h/B B C B C B C B C B C 309.38 280.65 266.23 257.17 251.24 247.06 243.84 241.41 239.31 237.50 236.19 236.19 270.63 258.48 249.37 245.46 244.24 248.57 259.35 250.86 253.80 253.97 247.01 250.70 247.53 247.01 243.44 568.55 515.75 489.26 472.60 461.70 454.02 448.11 443.63 439.79 436.46 434.04 434.04 497.35 475.02 458.26 451.09 448.85 456.81 476.61 461.01 466.41 466.72 453.93 460.71 454.88 453.93 447.37 376.73 341.74 324.19 313.15 305.93 300.84 296.92 293.95 291.41 289.20 287.60 287.60 329.55 314.75 303.65 298.90 297.41 302.68 315.81 305.47 309.05 309.25 300.78 305.27 301.41 300.78 296.43 692.32 628.02 595.76 575.47 562.20 552.85 545.65 540.20 535.52 531.47 528.52 528.52 605.61 578.42 558.02 549.28 546.56 556.24 580.36 561.36 567.94 568.32 552.74 561.00 553.90 552.74 544.75 73.78 101.93 124.57 143.69 162.59 177.93 193.80 206.73 218.71 232.13 244.52 257.90 132.09 150.83 188.93 200.94 215.76 187.72 148.94 164.58 165.19 172.06 180.85 193.14 191.01 192.34 208.63 135.59 187.32 228.93 264.07 298.80 326.98 356.15 379.91 401.92 426.58 449.35 473.94 242.75 277.18 347.20 369.26 396.50 344.98 273.71 302.46 303.58 316.19 332.35 354.93 351.02 353.46 383.40 63.76 86.54 105.58 123.00 139.41 152.25 165.81 175.94 188.21 117.18 159.03 194.03 226.04 256.20 279.79 304.71 323.33 345.87 69.95 96.43 119.32 140.55 161.16 177.96 195.71 209.56 225.83 128.54 177.21 219.27 258.29 296.17 327.04 359.65 385.11 415.01 112.36 127.82 160.08 171.80 183.79 156.90 145.36 141.72 140.46 152.51 154.27 163.40 158.59 159.79 173.59 206.49 234.90 294.18 315.72 337.74 288.33 267.13 260.45 258.12 280.27 283.51 300.29 291.43 293.65 319.01 125.50 144.67 185.47 201.15 217.42 182.15 159.42 162.70 162.05 172.02 179.57 188.85 183.67 184.44 202.24 230.63 265.86 340.84 369.65 399.56 334.74 292.97 298.99 297.81 316.11 329.99 347.04 337.52 338.95 371.66 i/A i/B B C B C 58.50 84.11 108.34 131.98 156.96 179.72 204.38 226.28 250.67 275.05 299.40 323.79 110.27 131.98 181.94 204.38 229.05 179.93 133.29 155.04 156.96 155.41 178.84 183.75 179.93 178.84 202.16 107.50 154.57 199.10 242.54 288.44 330.26 375.59 415.84 460.66 505.45 550.22 595.04 202.65 242.54 334.35 375.59 420.93 330.66 244.95 284.92 288.44 285.60 328.65 337.68 330.66 328.65 371.50 76.44 109.91 141.57 172.46 205.11 234.85 267.07 295.70 327.57 359.42 391.25 423.12 144.10 172.46 237.75 267.07 299.32 235.13 174.18 202.60 205.11 203.09 233.70 240.12 235.13 233.70 264.17 140.48 201.99 260.17 316.94 376.93 431.58 490.80 543.41 601.98 660.51 719.01 777.57 264.81 316.94 436.91 490.80 550.06 432.10 320.09 372.32 376.93 373.22 429.47 441.28 432.10 429.47 485.47 a/B means the enthalpy of vaporization was determined from eq. (a) using term B; values in bold are experimental values within ±10 and ±25% of values taken from the literature; values in italics are the experimental values within ±50% of values taken from the literature - 91 - (c) - 92 - 2-Meth-3-pent 4-Meth-2-pent 2-Meth-1-pent 2-Meth-1-but 3-Meth-1-but 2-Meth-2-but 2-Meth-1-prop 3-Hex 2-Meth-3-pent 4-Meth-2-pent 2-Meth-1-pent 2-Meth-1-but 3-Meth-1-but 2-Meth-2-but 2-Meth-1-prop 3-Hex 2-Oct 2-Hept 2-Hex 2-But 2-Prop 1-Dodec 1-Undec 1-Dec 1-Non 1-Oct 1-Hept 1-Hex 1-Pent 0 2-Oct 2-Hept 2-Hex 2-But 2-Prop 1-Dodec 1-Undec 1-Dec 1-Non 1-Oct 1-Hept 1-Hex 1-Pent 1-But 1-Prop Eth 1-But 1-Prop 2-Meth-3-pent 4-Meth-2-pent 2-Meth-1-pent 2-Meth-1-but 3-Meth-1-but 2-Meth-2-but 2-Meth-1-prop 3-Hex 2-Oct 2-Hept 2-Hex 2-But 2-Prop 1-Dodec 1-Undec 1-Dec 1-Non 1-Oct 1-Hept 1-Hex 1-Pent 1-But 1-Prop Eth Meth (a) Eth Meth molar heat capacity [J/(mol K)] (b) Meth molar heat capacity [J/(mol K)] molar heat capacity [J/(mol K)] 700 600 500 400 300 200 100 0 0 Exp. T = 323K Ref. +/- 10% Ref. +/- 25% Ref. +/- 50% 700 600 500 400 300 200 100 Exp. T = 373K Ref. +/- 10% Ref. +/- 25% Ref. +/- 50% 700 600 500 400 300 200 100 Exp. T = 423K Ref. +/- 10% Ref. +/- 25% Ref. +/- 50% Fig. 2 Molar heat capacity (Cp) values calculated from experimental gas chromatographic data by use of fitting term B from eq. (I), with the indicated ranges of the literature data (DB-Wax) (The enthalpy of vaporization values originated from eq. (i), term A). (a) Tc = 323 K, (b) Tc = 373 K, (c) Tc = 423 K (c) - 93 - 2-Meth-3-pent 4-Meth-2-pent 2-Meth-1-pent 2-Meth-1-but 3-Meth-1-but 2-Meth-2-but 2-Meth-1-prop 3-Hex 1-But 2-Meth-3-pent 4-Meth-2-pent 2-Meth-1-pent 2-Meth-1-but 3-Meth-1-but 2-Meth-2-but 2-Meth-1-prop 3-Hex 2-Oct 2-Hept 2-Hex 2-But 2-Prop 1-Dodec 1-Undec 1-Dec 1-Non 1-Oct 1-Hept 1-Hex 1-Pent 0 2-Oct 2-Hept 2-Hex 2-But 2-Prop 1-Dodec 1-Undec 1-Dec 1-Non 1-Oct 1-Hept 1-Hex 1-Pent 1-But 1-Prop Eth Eth 1-Prop 2-Meth-3-pent 4-Meth-2-pent 2-Meth-1-pent 2-Meth-1-but 3-Meth-1-but 2-Meth-2-but 2-Meth-1-prop 3-Hex 2-Oct 2-Hept 2-Hex 2-But 2-Prop 1-Dodec 1-Undec 1-Dec 1-Non 1-Oct 1-Hept 1-Hex 1-Pent 1-But 1-Prop Eth Meth (a) Meth molar heat capacity [J/(mol K)] (b) Meth molar heat capacity [J/(mol K)] molar heat capacity [J/(mol K)] 700 600 500 400 300 200 100 0 0 Exp. DB-1 Ref. +/- 10% Ref. +/- 25% Ref. +/- 50% 700 600 500 400 300 200 100 Exp. DB-5 Ref. +/- 10% Ref. +/- 25% Ref. +/- 50% 700 600 500 400 300 200 100 Exp. DB-Wax Ref. +/- 10% Ref. +/- 25% Ref. +/- 50% Fig 3 Molar heat capacity (Cp) values calculated from experimental gas chromatographic data by use of fitting term B from eq. (III), with the indicated ranges of the literature data (Tc = 348 K) (the enthalpy of vaporization values originated from eq. (i), term A). (a) DB-1, (b) DB-5, (c) DB-Wax re of vaporization was similar to the process of retention seems to predominate. The least polar stationary phase, DB-1, furnishes retention data that result in molar heat capacity values which deviate most from the reference data. This is another confirmation that Nature apparently favours intermolecular interactions among similar species. Models I and III proved to be those enabling prediction of the most accurate values of molar heat capacity. The good performance was independent of the stationary phase used and the fitting term B for both performed better than the fitting term C (probably because its full physical description includes the working temperature, Tc, which does not appear in C). The worst results were obtained from model II; irrespective of the stationary phase used the thermodynamic data obtained from its use deviated most from the reference values. An analogous ranking can be applied to the equations used for determination of the enthalpies of vaporization of the test analytes. Irrespective of the stationary phase used, eqs (g)–(i) proved to be the best performing. For DB-Wax eqs (d)–(f) also performed very well. With this series of relationships term C afforded better results than B (irrespective of the stationary phase and model used). An important relationship was noticed between the enthalpies of vaporization used for the computations performed in step 2 and the resulting molar heat capacities – the less the enthalpies of vaporization deviated from the reference data, the better were the heat capacities. The opposite was also noticed – enthalpies of vaporization that deviated strongly from the reference data usually resulted in poor molar heat capacity values. Very sporadically the opposite behaviour was observed and can be regarded as a meaningless effect of the numbers. CONCLUSIONS The mathematical models used proved physicochemically sound and performed well statistically. Their usefulness for prediction of retention data and thermodynamic values to some extent fulfilled our expectations. The models used performed best for experimental data derived from the most polar stationary phase (i.e. DB-Wax) and the best coincidence of our results and reference data was observed for the lowest working temperatures used in this study. Models (I) and (III) proved best for determination of molar heat capacities, and B proved to be the best fitting term. - 94 - The enthalpies of vaporization derived from eqs (g)–(i) resulted in the best numerical values of molar heat capacity for the alcohols studied (i.e. the values closest to the reference data). Now term C was found to perform better than term B). Combination of capillary gas chromatography with mathematical models has once again proved an effective means of acquisition the thermodynamic data. This time the approach was used to determine the molar heat capacities of aliphatic alcohols, relatively polar analytes. This approach is certainly a good alternative to traditional, well established calorimetric methods. REFERENCES [1] K. Ciążyńska-Halarewicz and T. Kowalska, Acta Chromatogr., 16, 119 (2006) [2] K. Ciążyńska-Halarewicz and T. Kowalska, Acta Chromatogr., 15, 97 (2005) [3] K. Ciążyńska-Halarewicz, E. Borucka, and T. Kowalska, Acta Chromatogr., 12, 49 (2002) [4] T. Kowalska, Acta Chromatogr., 11, 7 (2001) [5] K. Ciążyńska-Halarewicz and T. Kowalska, Acta Chromatogr., 13, 81 (2003) [6] K. Ciążyńska-Halarewicz and T. Kowalska, J. Chromatogr. Sci., 41, 467 (2003) [7] K. Ciążyńska-Halarewicz, M. Helbin, P. Korzenecki, and T. Kowalska, Mathematical models of solute retention in gas chromatography as sources of thermodynamic data. Part IV. Alcohols as the test analytes (paper in preparation) [8] CRC Handbook of Chemistry and Physics, 81st Edition, CRC Press, 2000 [9] E.S. Domalski and E.D. Hearing, J. Phys. Chem. Ref. Data, 19, 881 (1990) [10] M. Zábranský, V. Růžička, and V. Majer, J. Phys. Chem. Ref. Data, 19, 719 (1990) [11] G. Défayes, K.S. Reddy, A. Dallos, and E. sz. Kováts, J. Chromatogr. A, 699, 131 (1995) [12] Poradnik fizykochemiczny, WNT, Warszawa, 1974 - 95 - [13] P. Korzenecki, An alternative approach to determination of the thermodynamic data for the homologues and congeners from the retention data derived by means of the capillary gas chromatography, combined with the semi-empirical mathematical models. The heat capacity of aliphatic alcohols (in Polish), MSc Thesis, The University of Silesia, Katowice, Poland, 2006 - 96 -
© Copyright 2026 Paperzz