{Prediction of excess molar volumes of binary mixtures by Prigogine-Flory-Patterson (PFP) and Extended Real Association Solution (ERAS) models} J. Serb. Chem. Soc. 82 (12) 1379–1390 (2017) UDC 544.038.2+544.344.2’032:53.02:544.3 JSCS–5048 Original scientific paper 1379 Prediction of excess molar volumes of binary mixtures by Prigogine–Flory–Patterson (PFP) and extended real association solution (ERAS) models• IVONA R. RADOVIĆ#, NIKOLA D. GROZDANIĆ#, BOJAN D. DJORDJEVIĆ#, SLOBODAN P. ŠERBANOVIĆ# and MIRJANA LJ. KIJEVČANIN*# Faculty of Technology and Metallurgy, University of Belgrade, Karnegijeva 4, 11120 Belgrade, Serbia (Received 17 August, revised 31 August, accepted 8 September 2017) Abstract: In this paper, the Prigogine–Flory–Patterson (PFP) theoretical model and the extended real association solution (ERAS) model for the prediction of the excess molar volumes of 17 binary sets of alcohols (methanol, ethanol, pro- pan-1-ol, butan-1-ol, butan-2-ol and pentan-1-ol) with benzene, chlorobenzene, acetonitrile, hexan-1-amine and n-heptane were tested at 298.15 K and atmo- spheric pressure. The binary sets contained 333 experimental data points for 11 different compounds. The ERAS model demonstrated better correlative charac- teristics for binary mixtures comparing to PFP model. The ERAS model is successfully applied on all binary systems except ethanol + benzene and pro- pan-1-ol + benzene. PFP model gave less satisfactory results for all systems except butan-1-ol + hexan-1-amine, for which obtained results were acceptable. Keywords: binary mixtures; excess molar volume; prediction model; PFP; ERAS. INTRODUCTION Design of process and industrial equipment is based on the knowledge of the thermodynamic properties of fluids. The knowledge of thermodynamic and trans- port properties of pure components and their mixtures is very important for understanding intermolecular interactions in binary mixtures. Solvent mixtures play an important role in chemical industries and their physiochemical properties often show large deviations from ideal behaviour. Data of this properties are important for various engineering calculations and represent the starting points for the designing, analysis and optimization of technological processes. In addi- tion to the experimental determination of thermodynamic parameters of pure sub- * Corresponding author. E-mail: mirjana@tmf.bg.ac.rs # Serbian Chemical Society member. • This paper is dedicated to Prof. Bojan D. Djordjević on the occasion of his 80th birthday. https://doi.org/10.2298/JSC170817103R ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. 1380 RADOVIĆ et al. stances and their mixtures, it is necessary to develop and test a reliable thermos- dynamic model to replace, in some cases, very expensive experiments. Prigogine–Flory–Patterson (PFP) model1–7 is the theoretical approach that can describe and predict the non-ideal behaviour of multicomponent systems. On the other hand, the original extended real associated solution (ERAS) model, introduced by Heintz,8 combines the association effects arising from hydrogen bonding with free volume effects from the differences of van der Waals interact- ions between unlike molecules in a mixture. The aim of this work was to test the applicability of these two models to correlate the excess molar volume data of 17 binary systems containing the widely applied organic compounds, alcohols (methanol, ethanol, propan-1-ol, butan-1-ol, butan-2-ol and pentan-1-ol) and ben- zene, chlorobenzene, acetonitrile, hexan-1-amine and n-heptane.9–13 The analyzed systems are divided into five groups: alcohols + aromates, alcohols + n-alkane, alcohols + amine, alcohols + acetonitrile and amine + n-alkane. Selected systems show diverse deviations from ideal behaviour, e.g., highly negative, positive or S-shaped excess molar volumes vs. mixture composition curve. In this way, the applicability of the model for a specific type of non-ideal behaviour is confirmed. The modelling was performed at 298.15 K and atmospheric pressure. MODELLING In this paper predictive capabilities of PFP and ERAS models were tested for density calculations and excess molar volumes for large number of binary systems based on the pre- viously measured experimental data. The selected systems are characterized by different forms of non-ideal behaviour and the results of this work can show the applicability of both models for different groups of compounds and mixtures. The PFP model is designed to analyze the thermodynamic functions of the binary sys- tems. This theory is used for interpretation of the measurement results of some excess pro- perties of the binary systems formed by molecular species that differ in size and shape. Pure empirical application of the Flory model4 can provide some interesting correlations of excess molar properties for complex mixtures. According to PFP model, the excess molar volume contains three contributions: interactional contribution, a free volume contribution and a pres- sure contribution. The expression for excess molar volume which separates the three contri- butions is given as: ( ) ( ) ( )( ) ( ) ( )( ) ( )( ) ( )( ) 2 1/3 2/3 1/3 1 2 12 1 2 1 2 * * 1/3 * 1/31 1 2 2 1 * * 1 2 1 2 1 2 * * 1 2 2 1 1 14 / 9 / 1 4 / 3 / 1 4 / 3 / 1 E V V V V VV x V x V V P V V V V P P P P ψ θ χ ψ ψ ψ ψ ψ ψ − − − = − + + − − − − + +           (1) where VE is excess molar volume, x mole fraction, V* characteristic volume, P* characteristic pressure, ψ molecular contact energy fraction, θ molecular surface fraction and χ12 parameter. Reduced volume of binary mixture is defined as:14 1 1 2 2V V VΨ Ψ= +   (2) ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. PFP AND ERAS MODELS 1381 The characteristic parameters V* and P* are obtained from thermal expansion coefficient (1 / )( / )PV V Tα = ∂ ∂ and isothermal compressibility, (1 / )( / )V V P τκ = ∂ ∂ : ( )( ) ( ) 3 1 4 / 3 , 1, 2 1 i i i T V i T α α  + = =  +   (3) * /i i iV V V= (4) 2 * , 1, 2i iii ii P TV P i P α κ = = =   (5) The molecular contact energy fraction is calculated by: ( ) * 1 1 1 2 * * 1 1 2 2 1 P P P ϕ Ψ Ψ ϕ ϕ = − = + (6) ( ) * 1 1 1 2 * * 1 1 2 2 1 x V x V x V ϕ ϕ= − = + (7) where ϕ is a hard-core volume fraction. The molecular surface fraction is calculated by: 1 1 11 2 1 1 2 2 1 2 1 2 1 ( / ) S S S S S ϕ ϕ θ θ ϕ ϕ ϕ ϕ = − = = + + (8) where Si is the molecular fraction surface/volume ratio for the components determined by Bondi’s method:15 1/3* 1 1 * 2 2 S V S V −   =     (9) , ln ρ α ∂ = − ∂ p xT (10) 2 21 / ( ) p T u C α κ ρ ρ = + (11) The values of the parameters for the pure liquid components and the mixture are obtained using the Flory theory.4 The cross-interaction parameter required for the calculation of VE using PFP theory was derived from values of excess molar volumes at equimolar com- position for each mixture. On the other hand, ERAS model combines real association solution model16–19 with Flory equation of state. The excess function is written as: E E Echem physV V V= + (12) where VEchem is the chemical contribution, mainly due to the association reactions, and V E phys represents the physical contribution, which results from the physical interactions between molecules: ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. 1382 RADOVIĆ et al. ( ) ( ) ( ) ( ) E * 0 * 0 M A A A A1 B B B B11 1chem B1 A A1* A AB AB 1 B B B1 AB B1 1 1 A B A V V x K v x K v K x K v V V K K Δ φ φ Δ φ φ φ φ Δ φ φ− = − + − + − +  − +   (13) ( )( )E * *A A B B M A A B Bphys V x V x V V V Vφ φ= + − −   (14) where Ki are the association constants (i is A, B or AB) and it is assumed to be independent of the chain length. Their temperature dependence is given by: * i i 0 0 1 1 exp h K K R T T Δ    = − −        , i = A, B or AB (15) where K0 is the equilibrium constant at the atmospheric pressure and standard temperature T = = 298.15 K and ∆hi * the enthalpy change of reaction, which corresponds to the hydrogen bonding energy. These reactions are also characterized by the volume change ∆vi *, related to the formation of the linear chains. ERAS model is used to correlate the experimental values for the investigated mixtures and needs pure component properties to obtain the characteristic parameters V* and P* of each compound. It also needs data for density, ρ, isobaric thermal expansivity, α, and isothermal compressibility, κ, at atmospheric pressure and T = 298.15 K. RESULTS AND DISCUSSION The goal of this research was to check if PFP and ERAS models are capable to correlate the excess molar volumes of different binary mixtures containing alcohols, aromates, n-alkanes, amines and nitriles. Table I shows the literature density and the molar mass data for all constituents of investigated binary systems at atmospheric pressure and T = 298.15 K. TABLE I. Densities, ρ, and molar masses, M, of the pure compounds at 298.15 K and atmospheric pressure Compound ρ / g cm-3 M / g mol-1 Methanol9 0.78665 32.04 Ethanol9 0.78525 46.07 Propan-1-ol10 0.799692 60.09 Butan-1-ol11 0.805762 74.12 Butan-2-ol12 0.802528 74.12 Pentan-1-ol13 0.810968 88.15 Acetonitrile9 0.77669 41.05 Hexan-1-amine11 0.760128 101.19 Benzene10 0.873582 78.11 Chlorobenzene13 1.101038 112.56 n-Heptane11 0.679494 100.02 The values of parameters for the pure liquid components required for both PFP and ERAS models were obtained from the literature.20–26 Table S-I of the Supplementary material to this paper lists the literature for pure component data, ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. PFP AND ERAS MODELS 1383 self-association constant, K, thermal expansion coefficient, α, isothermal com- pressibility κ, reduced pressure P*, reduced volume V*, molar enthalpy change of association, ∆h*, molar volume change of association, ∆v*, and surface volume ratio S at atmospheric pressure and T = 298.15 K. The PFP parameter, χ12, in Eq. (1), can be calculated at given temperature by fitting the experimental excess molar volumes of binary mixture over the entire composition range. The values of the interaction parameter χ12 for the studied mixture have been calculated using the equimolar VE values and they are listed in Table S-II (Supplementary material) together with the equimolar values for all contributions to VE (interactional, free volume and P* effect contribution). The parameters KAB, ∆vAB*, and XAB in the ERAS model were correlated and listed in Table S-III (Supplementary material). The quality of the fit for the applied models was estimated by the standard deviation, σ and average percent deviation, D: ( ) 0.5 2 exp 1 N E E cal i V V N npar σ =   −    =  −       (16) E E exp cal E exp1 100 N i V V D N V= − =  (17) where npar is the number of parameters. The obtained results with PFP and ERAS models, for all systems, are listed in Table II. The results obtained for the selected binary systems, representing each group of analyzed binary mixtures, are presented in Figs. 1–7. The systems selected for the analyses are divided into five groups: alcohols + aromates, alcohols + n-alkane, alcohols + amine, alcohols + acetonitrile and amine + n-alkane. Selected systems show different deviations from ideal behaviour, e.g., negative, positive or S-shape (positive and negative) VE vs. x1 curves. In all cases the ERAS model was significantly superior comparing to PFP model. The systems of alcohols with aromates were the largest group in this res- earch, i.e., five systems with benzene and six systems with chlorobenzene were analyzed. The systems of analyzed alcohols with chlorobenzene exhibit S-shape VE vs. x1 curve. As it was already discussed in the literature10 the addition of large amounts of chlorobenzene to a pure alcohol induces the de-polymerization of normal alcohols, resulting in volume expansion and positive VE values. The con- traction in volume is characteristic for higher alcohol concentrations due to the ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. 1384 RADOVIĆ et al. formation of OH–Cl atom hydrogen bonded complexes, between an alcohol and chlorobenzene. Structural effects can also contribute to closer geometrical pack- ing and volume contraction, which is the prevailing factor in the mixture with methanol, since methanol is the smallest molecule. The systems with benzene (Figs. 1 and 2) also show sigmoidal VE vs. x1 curves, but the negative effects are less pronounced comparing to the equivalent systems with chlorobenzene. The possible interactions between alcohols and benzene aromatic ring are the OH-π type27, which are much weaker comparing to OH-Cl interactions between alco- hols and chlorobenzene. It is obvious that PFP model is not capable to correlate VE experimental data of those systems, with errors exceeding 100 % (Fig. 2). On the other hand ERAS model was rather successful in some cases. For the alcohol + chlorobenzene the systems percentage deviations are ranged from about 7 to 15 %, while for the alcohol + benzene systems the deviations were higher especially for the systems containing ethanol or 1-propanol (see Fig. 2). Although the errors for the ethanol + benzene system are high (10.97 %), Fig. 1 confirms that ERAS model fits experimental point rather well, except in the dilute regions. TABLE II. Average relative deviation, D and the standard deviation, σ obtained with PFP and ERAS models Mixture PFP ERAS D / % σ / cm3 mol-1 D / % σ / cm3 mol-1 Methanol + benzene27 103.25 60.097 10.97 0.0017 Ethanol + benzene28 203.20 51.452 217.45 0.0276 Propan-1-ol + benzene9 1120.70 37.7 181.05 0.0203 Butan-1-ol + benzene29 56.65 21.7 12.65 0.0096 Butan-2-ol + benzene31 18.34 9.55 6.71 0.0194 Methanol + chlorobenzene27 865.59 13.7 7.28 0.0030 Ethanol + chlorobenzene27 202.61 17.6 14.02 0.0051 Propan-1-ol + chlorobenzne12 156.60 26.1 8.53 0.0062 Butan-1-ol + chlorobenzene12 249.09 32.8 15.18 0.0061 Butan-2-ol + chlorobenzene12 229.09 36.5 6.72 0.0030 Pentan-1-ol + chlorobenzene12 158.62 37.9 14.17 0.0063 Butan-1-ol + n-heptane30 45.45 20.9 9.87 0.0078 Butan-2-ol + n-heptane32 20.18 10.2 5.36 0.0160 Butan-1-ol + hexan-1-amine11 2.66 1.41 1.95 0.0129 Methanol + acetonitrile7 26.11 15.4 3.03 0.0026 Ethanol + acetonitrile27 240.09 39.0 29.06 0.0037 Hexan-1-amine + n-heptane9 20.90 8.95 3.16 0.0042 Two systems containing an alcohol (butan-1-ol and butan-2-ol) and n-hept- ane and one system of an amine (hexan-1-amine) and n-heptane were also ana- lyzed. All systems exhibit positive VE vs. x1 trend. In these mixtures the unlike molecules are not capable to make heteroassociates due to hydrogen bonding. Con- ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. PFP AND ERAS MODELS 1385 0.0 0.2 0.4 0.6 0.8 -0.02 -0.01 0.00 0.01 0.02 0.03 10 6 V E /m 3 m ol -1 x 1 0.0 0.2 0.4 0.6 0.8 -0.03 -0.02 -0.01 0.00 0.01 0.02 0.03 0.04 10 6 V E /m 3 m ol -1 x 1 a) b) Fig. 1. The experimental excess molar volumes compared with the correlated values via PFP theory and ERAS model for binary mixture of methanol + benzene at 298.15K; a) PFP theory and b) ERAS model (- - -) physical contribution, (– ·–) chemical contribution. 0.0 0.2 0.4 0.6 0.8 1.0 -0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 0.08 10 6 V E / m 3 m ol =1 x 1 0.0 0.2 0.4 0.6 0.8 1.0 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 10 6 V E /m 3 m ol -1 x 1 a) b) Fig. 2. The experimental excess molar volumes compared with the correlated values via PFP theory and ERAS model for binary mixture of ethanol + benzene at 298.15K; a) PFP theory and b) ERAS model (- - -) physical contribution, (– ·–) chemical contribution. trary to that, both alcohol and amine molecules are making homoassociates. It is assumed that n-heptane added to a pure alcohol/amine causes disruption of alcohol or amine homoassociates and volume expansion, i.e., positive VE values. Increase in the volume in butan-2-ol + n-heptane mixture is higher, compared to butan-1-ol + n-heptane system, while both VE curves are slightly skewed towards mixtures rich in n-heptane. The same trend, the VE–x1 shape, is also characteristic for the hexan-1-amine + n-heptane system. Both models work better for the sys- tem butan-2-ol + n-heptane, since the VE curves are more symmetric. The average percentage deviation for the ERAS model and butan-1-ol + n-heptane system ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. 1386 RADOVIĆ et al. goes up to 10 %, though ERAS model curve follows experimental data very well (see Figs. 3 and 4). 0.0 0.2 0.4 0.6 0.8 0.00 0.05 0.10 0.15 0.20 10 6 V E / m 3 m ol -1 x 1 0.0 0.2 0.4 0.6 0.8 -0.05 0.00 0.05 0.10 0.15 0.20 0.25 1 06 V E /m 3 m ol -1 x 1 a) b) Fig. 3. The experimental excess molar volumes compared with the correlated values via PFP theory and ERAS model for binary mixture of butan-1-ol + n-heptane at 298.15K; a) PFP theory and b) ERAS model(- - -) physical contribution, (– ·–) chemical contribution. 0.0 0.2 0.4 0.6 0.8 0.0 0.1 0.2 0.3 0.4 0.5 10 6 V E /m 3 m ol -1 x 1 0.0 0.2 0.4 0.6 0.8 0.0 0.1 0.2 0.3 0.4 0.5 10 6 V E /m 3 m ol -1 x 1 a) b) Fig. 4. The experimental excess molar volumes compared with the correlated values via PFP theory and ERAS model for binary mixture of butan-2-ol + n-heptane at 298.15K; a) PFP theory and b) ERAS model(- - -) physical contribution, (– ·–) chemical contribution. In Figs. 5 and 6 two systems of methanol and ethanol with acetonitrile are shown. The excess molar volumes for the system with methanol show small negative values, while the VE–x1 curve for the system with ethanol is S-shaped, i.e., positive for the mixtures rich in acetonitrile and negative for higher ethanol concentrations. From shown figures it is obvious that ERAS model succeeds to follow the S-shaped VE curve, although with much higher average relative devi- ation, comparing to the system with ethanol. On the other hand PFP model com- ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. PFP AND ERAS MODELS 1387 pletely fails to describe the complex VE behaviour of ethanol + acetronitrile sys- tem. For the system with methanol the situation is better, however the PFP model still does not follow the asymmetric VE vs. x1 trend. 0.0 0.2 0.4 0.6 0.8 -0.16 -0.14 -0.12 -0.10 -0.08 -0.06 -0.04 10 6 V E / m 3 m ol -1 x 1 0.0 0.2 0.4 0.6 0.8 -0.16 -0.14 -0.12 -0.10 -0.08 -0.06 -0.04 -0.02 0.00 10 6 V E /m 3 m ol -1 x 1 a) b) Fig. 5. The experimental excess molar volumes compared with the correlated values via PFP theory and ERAS model for binary mixture of methanol + acetonitrile at 298.15K; a) PFP theory and b) ERAS model(- - -) physical contribution, (– ·–) chemical contribution. 0.0 0.2 0.4 0.6 0.8 1.0 -0.08 -0.06 -0.04 -0.02 0.00 0.02 10 6 V E /m 3 m ol -1 x 1 0.0 0.2 0.4 0.6 0.8 -0.10 -0.08 -0.06 -0.04 -0.02 0.00 0.02 0.04 10 6 V E /m 3 m ol -1 x 1 a) b) Fig. 6. The experimental excess molar volumes compared with the correlated values via PFP theory and ERAS model for binary mixture of ethanol + acetonitrile at 298.15K; a) PFP theory and b) ERAS model (- - -) physical contribution, (– ·–) chemical contribution. The mixtures of alcohols and amines show highly negative excess molar volume data with the symmetric VE vs. x1 curve,9 due to the hydrogen bonding and hetero-association between alcohol and amine molecules. This type of non- ideal behaviour can be successfully fitted by both PFP and ERAS models, yet the ERAS model fits minimal region slightly better than PFP model. From the all mentioned above it is obvious that except for the alcohol + amine mixtures, the ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. 1388 RADOVIĆ et al. PFP model was not able to describe the different VE–x1 curve shapes. The best results were achieved for the hexan-1-amine + n-heptane system (see Fig. 7). 0.0 0.2 0.4 0.6 0.8 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 10 6 V E /m 3 m ol -1 x 1 0.0 0.2 0.4 0.6 0.8 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 10 6 V E /m 3 m ol -1 x 1 a) b) Fig. 7. The experimental excess molar volumes compared with the correlated values via PFP theory and ERAS model for binary mixture of hexan-1-amine + n-heptane at 298.15K; a) PFP theory and b) ERAS model (- - -) physical contribution, (– ·–) chemical contribution. CONCLUSIONS In this paper, two theoretical models, the Prigogine–Flory–Patterson (PFP) and the extended real association solution (ERAS) for the excess molar volume correlation were tested on 17 binary mixture systems divided into five groups: alcohols + aromates, alcohols + n-alkane, alcohols + amine, alcohols + aceto- nitrile and amine + n-alkane. The excess molar volumes for binary mixtures are described qualitatively for all ranges of mole fractions by the PFP and the ERAS model. Both models can be applied on the systems that have only positive or negative excess molar volume values over the whole range of mole fractions, where the ERAS model gives better results. Also, the ERAS model was much better comparing to the PFP for systems that have both, positive and negative excess molar volume values, over the whole range of mole fractions. The PFP theory is able to describe the systems with only negative or only positive excess molar volume values over the whole range of mole fractions. However, if the VE– –x1 curve is asymmetric and tends to shifts towards higher or lower concentration than an equimolar point, the errors are getting higher. The comparison of cal- culated and experimental excess properties of binary systems shows that the ERAS model is able to describe all the excess functions ranging from positive values to strong negative values, as well as the S-shaped VE–x1 curves, which cannot be completed for the PFP model. ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. PFP AND ERAS MODELS 1389 SUPPLEMENTARY MATERIAL Tables S-I–S-III are available electronically at the pages of journal website: http:// //www.shd.org.rs/JSCS/, or from the corresponding author on request. Acknowledgement. This work was financially supported by the Research Fund of Ministry of Education, Science and Environmental Protection, Serbia, the Faculty of Tech- nology and Metallurgy, University of Belgrade, Serbia (project No 172063). И З В О Д ОДРЕЂИВАЊЕ ДОПУНСКЕ МОЛАРНЕ ЗАПРЕМИНЕ БИНАРНИХ СИСТЕМА ПОМОЋУ PRIGOGINE–FLORY–PATTERSON (PFP) И EXTENDED REAL ASSOCIATION SOLUTION (ERAS) МОДЕЛА ИВОНА Р. РАДОВИЋ, НИКОЛА Д. ГРОЗДАНИЋ, БОЈАН Д. ЂОРЂЕВИЋ, СЛОБОДАН П. ШЕРБАНОВИЋ и МИРЈАНА Љ. КИЈЕВЧАНИН Технолошко–металуршки факултет, Универзитет у Београду, Карнегијева 4, 11120 Београд У овом раду је извршено израчунавање допунских моларних запремина 17 бинар- них система алкохола (метанол, етанол, пропан-1-ол, бутан-1-ол, бутан-2-ол и пентан-1- -ол) са бензеном, хлорбензеном, ацетонитрилом, хексан-1-амином и n-хептаном на температури од 298,15 K и на атмосферском притиску, коришћењем теоријских модела: Prigogine–Flory–Patterson (PFP) и Extended Real Association Solution (ERAS). Бинарни системи садржали су 333 експерименталне тачке за 11 различитих једињења. ERAS модел показао је боље корелативне карактеристике за бинарне смеше у поређењу са PFP моделом. ERAS модел успешно је примењен на свим бинарним системима, осим на смеше етанол + бензен и пропан-1-ол + бензен. PFP модел генерално даје лошије резултате за све системе, осим за систем бутан-1-ол + хексан-1-амин, за који су доби- јени прихватљиви резултати. (Примљено 17. августа, ревидирано 31. августа, прихваћено 8. септембра 2017) REFERENCES 1. I. Prigogine, A. Bellemans, V. Mathot, The Molecular Theory of Solutions, North Holland, Amsterdam, 1957 2. P. J. Flory, R. A. Orwoll, A. Vrij, J. Am. Chem. Soc. 86 (1964) 3507 3. P. J. Flory, R. A. Orwoll, A. Vrij,J. Am. Chem. Soc. 86 (1964) 3515 4. P. J.Flory, J. Am. Chem. Soc. 87 (1965) 1833 5. A. Abe, P. J.Flory, J. Am. Chem. Soc. 87 (1965) 1838 6. R. A. Orwoll, P. J.Flory, J. Am. Chem. Soc. 89 (1967)6814 7. D. Patterson, G. Delmas, Discuss. Faraday Soc. 49 (1970) 98 8. A. Heintz, Ber. Bunsenges. Phys. Chem. 89 (1985) 172 9. I. R. Grgurić, A. Ž. Tasić, B. D. Djordjević, M. Lj. Kijevčanin, S. P. Šerbanović, J. Serb. Chem. Soc. 67 (2002) 581 10. M. Lj. Kijevčanin, I. M. Purić, I. R. Radović, B. D. Djordjević, S. P. Šerbanović, J. Chem. Eng. Data 52 (2007) 2067 11. I. R. Radović, M. Lj. Kijevčanin, S. P. Šerbanović, B. D. Djordjević, Fluid Phase Equilib. 298 (2010) 117 12. I. R. Radović, M. Lj. Kijevčanin, E. M. Djordjević, B. D. Djordjević, S. P. Šerbanović, Fluid Phase Equilib. 263 (2008) 205 13. I. R. Radović, M. Lj. Kijevčanin, A. Ž. Tasić, B. D. Djordjević, S. P. Šerbanović, J. Serb. Chem. Soc. 74 (2009) 1303 ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. 1390 RADOVIĆ et al. 14. Y. Sun, L. Su, H.Wang, J. Chem. Thermodynamics 41 (2009) 1154 15. A. Bondi, Physical Properties of Molecular Crystals, Liquids and Glasses, Wiley, New York, 1968 16. C. B. Kretschmer, R. Wiebe, J. Chem. Phys. 22 (1954) 1697 17. H. Renon, J.M. Prausnitz, Chem. Eng. Sci. 22 (1967) 299 18. H. V. Kehiaian, Bull. Acad. Pol. Sci. 16 (1968) 165 19. H. V. Kehiaian, A. J. Treszczanowicz, Bull. Acad. Pol. Sci. 16 (1968) 171 20. M. Bender, A. Heintz, Fluid Phase Equilib. 89 (1993) 197 21. E. N. Rezanova, K. Kammerer, R. N. Lichtenthaler, J. Chem. Thermodynamics 32 (2000) 1569 22. A. Heintz, P. K. Naicker, S. P. Verevkin, R. Pfestorf, Ber. Bunsen – Ges. Phys. Chem. 102 (1998) 953 23. R. B. Torres, A. Z. Francesconi, P. L. Volpe, Fluid Phase Equilib. 210 (2003) 287 24. U. Domanska, M. Gloskowska, Fluid Phase Equilib. 216 (2004) 135 25. B. Orge, M. Iglesias, G. Marino, M. Domınguez, M. M. Pineiro, J. Tojo, Fluid Phase Equilib. 170 (2000) 151 26. S. P. Ijardar, N. I. Malek, S. L. Oswal, Ind. J. Chem., A 50 (2011) 1709 27. S. P. Šerbanović, M. Lj. Kijevčanin, I. R. Radović, B. D. Djordjević, Fluid Phase Equilib. 239 (2006) 69 28. I. R. Grgurić, S. P. Šerbanović, M. Lj. Kijevčanin, A. Ž. Tasić, B. D. Djordjević, Thermochim. Acta 412 (2004) 25 29. J. D. Smiljanić, M. Lj. Kijevčanin, B. D. Djordjević, D. K. Grozdanić, S. P. Šerbanović, Int. J. Thermophys. 29 (2008) 586 30. M. Lj. Kijevčanin, I. R. Radović, S. P. Šerbanović, A. Ž. Tasić, B. D. Djordjević, Thermochim. Acta 496 (2009) 71 31. 31. J. D. Smiljanić, M. Lj. Kijevčanin, B. D. Djordjević, D. K. Grozdanić, S. P. Šerbanović, J. Chem. Eng. Data 53 (2008) 1965 32. M. Lj. Kijevčanin, I. R. Radović, S. P. Šerbanović, E. M. Živković, B. D. Djordjević, J. Chem. Eng. Data 55 (2010) 1739. ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. 1391 A tribute to Professor Bojan Djordjević in honour of his 80th birthday Journal of the Serbian Chemical Society wishes to draw attention to Professor Djordjević’s rich legacy and superior contribution to the field of ther- modynamics in general and to JSCS in particular. Prof. Dr. Bojan D. Djordjević was born in 1937 in Belgrade where he graduated from the Faculty of Technology, University of Belgrade in 1961 and commenced his career as Research Fellow in IHTM, Belgrade. He stayed in IHTM for the period from 1961 to 1964. He received his PhD in 1970 from the Institute of Nuclear Sciences “Vinča“ and the Fac- ulty of Technology where in 1971 he was promoted to the position of the Assistant Professor and in 1976 to the position of the Associate Professor. Since 1982 he had been Full Professor at the Faculty of Technology and Metallurgy, where he had been extensively involved in curricula development and teaching of various subjects of thermodynamics, including but not limited to chemical engineering thermodynamics, energy of process industry, and mole- cular thermodynamics. During his fruitful carrier he also mentored hundreds of graduate students supporting them selflessly to commence their own academic paths. Further, Professor Djordjević lectured different subjects of thermodyn- amics at many other university centers (Novi Sad, Podgorica, Priština, The Military–Technical Academy, Belgrade, Higher Technical School, Kosovska Mitrovica, etc.). He is the author of 27 text-books and monographs including Thermodynamics and Thermotechnics, Thermodynamics of Fluids, Chemical Engineering Thermodynamics, and Heat Transfer. Professor Djordjević was a foundational member of the Laboratory for chemical engineering parameters. In two terms he was the Head of Fundamental Technical Science Department and, also, led a national research projects and sub- projects. He contributed to over 50 technical studies and projects, many of which have been applied in chemical industry (Mineloprema, Pančevo; Oil Refinery, ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. 1392 Beolgrade; Oil Refinery, Pančevo; HIP, Pančevo; Miloje Zakić, Kruševac, Železara, Smederevo, etc.) and he also established cooperation with Imperial Chemical Industries, United Kingdom and University of Berkley, USA. His research has incorporated a wide range of fundamental areas, with a focus on equations of state, vapour–liquid equilibria, excess properties of solution and thermo-hydraulic design of heat equipment. Professor Djordjevic has also been very active in the Yugoslav Association of Chemists and Technologists. He is a member of Serbian Scientific Society, Serbian Chemical Society and American Chemical Society. He received a number of national and international acknow- ledgments, awards, plaques and certificates: Men of Achievement for Distin- guished Achievement, Cambridge, England, 1982.; 5.000 Personalities of the Word, American Biographical Institute; The International Who’ S Who of Intel- lectuals, Cambridge, England, 1983; Certificate of Merit for Distinguished Ser- vice to the Community, Cambridge, England, 1984.; Who’s Who in the World 1980/1981; Who’s Who in the World 1982/1983; Medal for outstanding and long-lasting contribution to science (for development of chemical engineering thermodynamics), Serbian Chemical Society, 2000; Honorary diploma on the 25th anniversary of the Faculty of Technology, University of Novi Sad, for a significant contribution to the work and development of the Faculty of Tech- nology, University of Novi Sad, 1984; Certificate as a sign of acknowledgment and gratitude for the great contribution to the work and development of the Serbian Chemical Society on the occasion of 100th anniversary 1997. Professor Djordjević authored and contributed to over 100 scientific articles, more than hounded scientific papers presented on international and national con- gresses and symposia and over 30 professional papers. His articles were pub- lished in international journals, such as: Journal of Chemical Thermodynamic, Fluid Phase Equilibria, Journal of Chemical Engineering Data, Chemical Eng- ineering Science, AIChE Journal, Canadian Journal of Chemical Engineering, Industrial & Engineering Chemistry Process Design and Development, Indus- trial & Engineering Chemistry Research, Thermochimica Acta, Physics and Chemistry of Liquids, Chemie Ingenieur Technik, Korean Journal of Chemical Engineering, Journal of Applied Polymer Science, Hungarian Journal of Ind- ustrial Chemistry, Chemische Technik, Journal of Chemical Education, etc. Being internationally recognized and respected, Professor Djordjević was cited over 500 times in the international journals, university textbooks and mono- graphs. He was a referee in the leading international and national journals and congresses and was a member of Editorial Boards of the Journal of the Serbian Chemical Society and Chemical Industry. Professor Djordjević has made an exceptional contribution to the field of thermodynamics throughout his career, and has also been graciously involved in supporting the Journal of the Serbian Chemical Society. ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. 1393 Being a member of Editorial Board of the Journal of the Serbian Chemical Society (more than 30 years) and publishing more than 40 scientific papers in this journal*, professor Djordjević provided an outstanding contribution to Journal of the Serbian Chemical Society development and reputation as well. It is therefore great pleasure to honour him in this tribute, on the celebration of his 80th birthday. We join in wishing him well in the future and thanking him for all the work and the high standards he set for all of us. Mirjana Kijevčanin Ivona Radović Slobodan Šerbanović * The paper JSCS–5048 in this issue is dedicated to Professor Djordjević on the occasion of his jubilee. ________________________________________________________________________________________________________________________Available on line at www.shd.org.rs/JSCS/ (CC) 2017 SCS. << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice