THE SPINODAL CURVE OF THE SYSTEM WATER 1-BUTANOL 1-PROPANOL ACCORDING TO THE WHEELER WIDOM MODEL

Similar documents
ON THE DETERMINATION OF THE SPINODAL CURVE FOR THE SYSTEM WATER + CHLOROFORM + ACETIC ACID FROM THE MUTUAL DIFFUSION COEFFICIENTS

A THREE-COMPONENT MOLECULAR MODEL WITH BONDING THREE-BODY INTERACTIONS

HYDRODYNAMIC AND THERMODYNAMIC ASPECTS OF DIFFUSION COEFFICIENTS IN THE TERNARY SYSTEM WATER-CHLOROFORM-ACETIC ACID AT 25 C

Florin-Dorian Buzatu [ISI] Articles:

Extraction of Phenol from Industrial Water Using Different Solvents

Measurement and Calculation of Physico-Chemical Properties of Binary Mixtures Containing Xylene and 1- Alkanol

DETERMINATION AND PREDICTION OF THE ISOBARIC VAPOR- LIQUID-LIQUID EQUILIBRIUM DATA

Liquid-liquid equilibrium for the quaternary system of o-xylene(1)+water(2) +butyric acid(3)+ethyl acetate(4) at 25 o C and atmospheric pressure

The successes and limitations of this thermodynamic model in predicting the solubility of organic solutes in pure solvents

Vol. 114 (2008) ACTA PHYSICA POLONICA A No. 6 A

Outline of the Course

Prediction of phase equilibria in waterralcoholralkane systems

Liquid liquid equilibria of aqueous mixtures containing selected dibasic esters and/or methanol

Lecture 1 Chapter 12 Sections 1-2. Solutions Solubility

Introduction (1) where ij denotes the interaction energy due to attractive force between i and j molecules and given by; (2)

Chapter 4 Polymer solutions

A NEW AND A SIMPLE MODEL FOR SURFACE TENSION PREDICTION OF WATER AND ORGANIC LIQUID MIXTURES * R. TAHERY AND H. MODARRESS **

Development new correlations for NRTL Parameters in Polymer Solutions

MECH 6661 lecture 9/1 Dr. M. Medraj Mech. Eng. Dept. - Concordia University

4 Results of the static and dynamic light scattering measurements

Thermodynamic of polymer blends Assoc.Prof.Dr. Jatyuphorn Wootthikanokkhan

Densities and Viscosities of the Ternary Mixtures Water + Butyl Acetate + Methanol and Water + Ethyl Propionate + Methanol at 303.

CH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions

The Characteristics of a Soln

LIQUID-LIQUID EQUILIBRIUM FOR THE DESIGN OF EXTRACTION COLUMN

A theoretical study of Gemini surfactant phase behavior

ADSORPTION IN MICROPOROUS MATERIALS: ANALYTICAL EQUATIONS FOR TYPE I ISOTHERMS AT HIGH PRESSURE

Chap. 2. Polymers Introduction. - Polymers: synthetic materials <--> natural materials

Phase Transitions in Multicomponent Systems

Functional Groups SCH4C

Phase Separation in Ternary Systems Induced by Crosslinking*

Vapor-liquid equilibrium

= = 10.1 mol. Molar Enthalpies of Vaporization (at Boiling Point) Molar Enthalpy of Vaporization (kj/mol)

Lab 3: Solubility of Organic Compounds

Wetting Transitions at Fluid Interfaces and Related Topics

MODELING OF PHASE EQUILIBRIA FOR BINARY AND TERNARY MIXTURES OF CARBON DIOXIDE, HYDROGEN AND METHANOL

CHEM 254 EXPERIMENT 4. Partial Molar Properties of Solutions

Interpreting Phase Diagrams

Property Prediction in Reactive Solutions

Theory of Interfacial Tension of Partially Miscible Liquids


Chapter 11 section 6 and Chapter 8 Sections 1-4 from Atkins

Announcements Please contribute 1-item to the Tondo Children Project

Excess Heat Capacity Surfaces for Water-Alkanol Mixtures by the UNIQUAC Model

Chemistry 2.5 AS WORKBOOK. Working to Excellence Working to Excellence

Unit 10: Part 1: Polarity and Intermolecular Forces

Thermodynamic stability and transformation of pharmaceutical polymorphs*

CHEM 254 EXP 10 Chemical Equilibrium - Homogeneous and Heterogeneous Equilibrium

Models and Algorithms for Condensed Phase Equilibrium Data Regression in Ternary Systems: Limitations and Proposals

Modeling of the solubility of Naproxen and Trimethoprim in different solvents at different temperature

Overall: 75 ECTS: 7.0

Molecular interaction studies of acrylic esters with alcohols

Mixtures. Partial Molar Quantities

COMMENTS CONCERNING STUDY OF MOLECULAR INTERACTION IN BINARY LIQUID MIXTURES OF ETHYL ACETOACETATE WITH CHLOROFORM AND

General Chemistry (Third Quarter)

Journal of Chemical and Pharmaceutical Research

Solving mass transfer problems on the computer using Mathcad

SOLUTIONS. Heterogeneous Mixtures. Section 8.1: Solutions and Other Mixtures. Heterogeneous Mixtures (cont d) CHAPTER 8.

PHASE TRANSITIONS IN SOFT MATTER SYSTEMS

Determination of the nonsolvent polymer interaction parameter χ 13 in the casting solutions

Lecture 1 Solubility and Distribution Phenomena

Process design using ionic liquids: Physical property modeling

Influence of the Temperature on the Liquid-Liquid- Solid Equilibria of the Water + Ethanol + 1-

3.012 PS 6 THERMODYANMICS SOLUTIONS Issued: Fall 2005 Due:

Chapter 2 Equilibria, Bubble Points, Dewpoints, Flash Calculations, and Activity Coefficients

Solution Experiment Collin College

Resolution of the ethanolamine salt of (±)mandelic acid by using the AS3PC method: Principle, application and results

Ch 12 and 13 Practice Problems

9.1 Mixtures and Solutions

Physical Pharmacy ( ) Unit 2 Phase Equilibria and the Phase Rule Solutions

Molecular Interaction Study of Binary Solutions of n-butyl Acetate and Isopropanol at Various Temperatures

Overview. Types of Solutions. Intermolecular forces in solution. Concentration terms. Colligative properties. Osmotic Pressure 2 / 46

THERMODYNAMIC CONSISTENCY TESTS FOR PHASE EQUILIBRIUM IN LIQUID SOLUTE+SUPERCRITICAL SOLVENT MIXTURES

Evaporation and Intermolecular Forces

States of matter Part 2

Ultrasonic Studies in Ternary Mixtures of Alcohols With NN DMA in N- Hexane at 303K

EXPERIMENT 7 - Distillation Separation of a Mixture

3.091 Introduction to Solid State Chemistry. Lecture Notes No. 9a BONDING AND SOLUTIONS

COSMO-RS Theory. The Basics

- intermolecular forces forces that exist between molecules

CORRELATION OF (LIQUID + LIQUID) EQUILIBRIUM OF SYSTEMS INCLUDING IONIC LIQUIDS

Status and results of group contribution methods

Some properties of water

Investigation of ternary systems

Chemistry General Chemistry II Spring 2006 Test #1

PAPER No.6: PHYSICAL CHEMISTRY-II (Statistical

DESIGN AND CONTROL OF BUTYL ACRYLATE REACTIVE DISTILLATION COLUMN SYSTEM. I-Lung Chien and Kai-Luen Zeng

2. What property of water allows a needle to float on it without sinking? Answer: surface tension

Simulation of Ethanol Dehydration Using Cyclohexane as an Entrainer

Solutions (Concentration)

EPSRC Centre for Doctoral Training in Industrially Focused Mathematical Modelling

SUPPLEMENTAL HOMEWORK SOLUTIONS WEEK 10

Chapter 11 Solution Thermodynamics: Theory

Prediction of Vapor/Liquid Equilibrium Behavior from Quantum Mechanical Data

SPECTRAL INVESTIGATION OF TRIPHENYLFORMAZAN DERIVATIVES IN ULTRAVIOLET LIGHT

UNIT 3: CONCEPTS OF CHEMICAL BONDING. Chapter Chapter

An Introduction to Two Phase Molecular Dynamics Simulation

CH1010 Exam #2 Study Guide For reference see Chemistry: An Atoms-focused Approach by Gilbert, Kirss, and Foster

Reliable Computation of Binary Parameters in Activity Coefficient Models for Liquid-Liquid Equilibrium

EQUATION OF STATE DEVELOPMENT

Transcription:

CONDENSED MATTER THE SPINODAL CURVE OF THE SYSTEM WATER 1-BUTANOL 1-PROPANOL ACCORDING TO THE WHEELER WIDOM MODEL ROBERTO SARTORIO 1, CRISTINA STOICESCU 2 1 Department of Chemical Science University Federico II of Naples, Complesso di Monte Sant Angelo, Via Cinthia 80126, Naples, Italy E-mail: roberto.sartorio@unina.it, r.sartorio45@libero.it 2 Institute of Physical Chemistry Ilie Murgulescu, Thermodynamic Department, Romanian Academy, 202 Splaiul Independenţei, Bucharest, 060021, Romania Received June 12, 2014 The spinodal curve of the system water 1-butanol 1-propanol has been determined, from the literature Liquid Liquid Equilibrium data (LLE), by using the Mean Field Approximation (MFA) on the Generalized Wheeler Widom Huckaby Model (GWWHM). The results are compared with those obtained for the system water 1-pentanol 1-propanol. Key words: Ternary system, amphiphilic solvent, Wheeler-Widom model, spinodal curve. PACS: 64.70.Ja, 05.50.+q, 05.10.-a, 82.60.Lf, 82.70.Uv. 1. INTRODUCTION The Liquid Liquid Equilibrium data (LLE) of a ternary system presenting a miscibility gap between the two solutes can be used to obtain the spinodal curve of the system by mean of the Generalized Wheeler Widom Huckaby Model (GWWHM) in the Mean Field Approximation (MFA) [1 8]. The original Wheeler Widom model [1] considers a ternary system as constituted by three types of rod-like di-functional molecules, AA, BB, and AB, that occupy the bonds of a regular lattice with a single molecule for each bond. In this model, where AA and BB are the solutes and AB is the amphiphilic solvent, the interactions A B are forbidden (ε AB = ), while the interactions A A and B B have the same energy (ε AA = ε BB = 0). At microscopic level, because only the interactions AA AA, BB BB, AA AB, and AB BB may exist in solution, it follows that the sites of the lattice can be occupied only by A ends or only by B ends. The possibility to identify in the lattice sites A and sites B allows a direct analogy of this WW model with the Ising model, defined on the same lattice, where, at each site, the spin can be 1/2 or 1/2. At macroscopic level we can deduct for such a model the following characteristics: 1. The condition (ε AB = ) implies that the two solutes are completely immiscible and then the binodal curve insists on the entire interval of binary composition. RJP Rom. 60(Nos. Journ. Phys., 7-8), Vol. 1068 1086 60, Nos. 7-8, (2015) P. 1068 1086, (c) 2015 Bucharest, - v.1.3a*2015.9.2

2 Spinodal curve of water 1-butanol 1-propanol system. Wheeler-Widom model 1069 This means, by using a Cartesian representation of the system with x = (X 3 X 2 )/ 3 and y = X 1 where X 1, X 2, and X 3 are respectively the mole fractions of the solvent and of the two solutes, that the binodal curve is defined in the whole range 1/ 3 x +1/ 3. 2. The possibility of the solvent to interact with both the solutes, by setting respectively interactions AA AB and AB BB, increases the mutual solubility of the two solutes. Moreover the condition (ε AA = ε BB = 0) implies that the solvent, when added in small quantity to a mixture of the two solutes, is equally redistributed between the two ternary solutions in thermodynamic equilibrium. As a consequence of this the tie-lines are parallel to the x axis, the solubility gap is symmetrical and the plait point is coincident with the maximum of the binodal curve. For large concentration of the solvent the ternary system becomes homogeneous. 3. No effects of the temperature can be detected. Although the model catches the main aspects of the phase equilibrium in ternary systems, it is too simplistic to describe real systems that usually present asymmetric binodal curves, tie-lines with slopes different from zero and temperature dependence. Several attempts to modify and extend the original model are present in literature [2 8]. The up-dated model used in this paper (GWWHM) accounts for a honeycomb lattice where two- and three-body interactions between the molecular ends exist at the lattice sites. This model is equivalent to the standard Ising model on a 3 12 lattice and can be solved exactly to obtain the binodal curve [3, 4]; it allows describing systems with a non-zero solubility between the solutes, with asymmetric binodal curves, with tie-lines having slopes different from zero. This model can also account for binodal curve dependence on the temperature. Moreover it allows, by using the Mean Field Approximation, to obtain the spinodal curve of the system; in this case, for internal consistency, also the binodal curve must be determined within the same approximation. The most important aspect of this model is that both the binodal and the spinodal curve can be described, in function of the composition of the system, by the same constant parameters τ and that correspond respectively, in the Ising model, to the reduced temperature and to the asymmetry parameter. As a consequence of this, the calculation of the spinodal curve should be very simple. It would be sufficient to fit the experimental LLE data to the theoretical binodal GWWHM equation in order to obtain the values (τ, ); subsequently this parameters, inserted in the theoretical GWWHM spinodal equation, should allow obtaining the spinodal curve. To test the possibility to predict, in this way, the spinodal curve starting from the LLE data we considered at first the system water (W) chloroform (CHCl 3 ) acetic acid (AcH) [9 11]. Even if, by structural point of view, the species involved

1070 Roberto Sartorio, Cristina Stoicescu 3 in this system are very different from the rod-like molecules used by the models, the experimental solubility gap is very similar to that predicted by the original WW model. The reciprocal solubilities of water in pure chloroform and of chloroform in pure water are very small and then the binodal curve practically insists on the entire range of x. The binodal is almost symmetrical with a flat maximum, the plait point is almost coincident with the maximum and tie-lines are almost parallel to the x axis. Despite these characteristics, after several attempts it was clear that the fit, to the GWWHM binodal equation, of the experimental binodal points in the entire binodal range of x gives poor results. For this reason we worked out a local fit procedure. The results obtained with this procedure were very good; the obtained spinodal passes through the plait point and was in very good agreement with the literature spinodal point of the system water chloroform. The reliability of the results obtained with the GWWHM model has been checked by the use of the mutual diffusion coefficients of system. In fact two spinodal compositions have been determined by using the thermodynamic condition that the determinant of the matrix of the diffusion coefficients of the system must be zero on the spinodal curve [12]. The four mutual diffusion coefficients D ij have been determined, in the homogeneous region of the system, at two constant values of the mole fraction ratio X W /X CHCl3 and varying X AcH : by fitting the D ij values vs. X AcH the spinodal composition was chosen as that one for which D ij = 0 [13, 14]. The agreement between the GWWHM prediction of the spinodal curve and the values obtained by the use of the diffusion technique was excellent. Afterwards we considered systems constituted by water and two 1-alcohols with different chain length. In particular we took into consideration ternary system as HO-H HO-(CH 2 ) n -H HO-(CH 2 ) 3 -H with n = 4,5,6,..., for which LLE data are reported in literature [15 17]; in the following we will indicate this systems with the shorter notation HOH (2) HOC n (3) HOC 3 (1). In this system water and the 1-alcohol with the longer chain length are the solutes while 1-propanol is the solvent. The first terms of the n-alcohol series methanol, ethanol, and 1-propanol are completely miscible with water. For n 4 the binary systems HOH HOC n present solubility gap [18]. The ternary systems HOH HOC n HOC m with m < n do not present a solubility gap if n = 1,2,3 and m = 1,2,... On the contrary a I type solubility gap is present if n 4 independently on the value of m [19]. For both n 4 and m 4 a II type solubility gap should be present. The molecular structure of the species involved in these systems is more like to that hypothesized by the WW and GWWH models. The hydrocarbon tail of the solvent and that of the alcoholic solute are practically the same only differing in the number of C atoms; moreover it is proved by different techniques that the hydrocarbon tails of several classes of organic solutes (alcohols, diols, carboxylic acids, amines, etc.) interact each other by setting hydrophobic bonds also in dilute aque-

4 Spinodal curve of water 1-butanol 1-propanol system. Wheeler-Widom model 1071 Fig. 1 Experimental LLE data; experimental binodal points, experimental tie-line end points, - - - - experimental tie-line, experimental plait point (P), experimental binodal point of the binary system HOH (2) HOC 4 (3), representative experimental binodal (REB). ous solutions [20 24]. On the other hand the hydroxyl group of the solvent is very similar to that of water and then the hydrophilic end of the solvent can easily set interactions also with the other solute of the system. In conclusion the solvent molecule can interact with one end with a solute and with the other end with the other solute according to the GWWHM model main hypothesis. It is to observe that solute solute interactions are also possible via the OH groups present in their molecules. Regarding the general characteristics of the binodal curves of these systems we observe that the solubility of water in the long chain alcohol is always much larger than the solubility of the same alcohol in water; all the binodal curves are then shifted to one side of the Cartesian triangle graph. The curves are almost symmetric but the plait point, which is always on the left part of the curve, is distant from the maximum and the tie-lines close to the plait point are very steep. Overall, these systems are very well suited to evaluate the applicability of the GWWHM in describing the properties of real ternary systems with a miscibility gap. In this paper we present the results obtained applying the GWWHM to the LLE experimental data of the system HOH (2) HOC 4 (3) HOC 3 (1); the data are reported in Figure 1.

1072 Roberto Sartorio, Cristina Stoicescu 5 2. EXPERIMENTAL LLE DATA Under the heading LLE data two different types of experimental data are collected in literature: i) the composition of single binodal points and ii) the composition of the two conjugated phases in chemical equilibrium [25 27]. The compositions of single binodal point can be easily determined, with small errors, by the titration method. To a heterogeneous opalescent mixture of the two solutes the solvent is added drop by drop until the opalescence disappears; known the amount and the composition of the initial mixture and the quantity of solvent necessary to obtain a clear solution, the binodal composition can be calculated. The opposite procedure can also be used; to a homogeneous mixture of the solvent plus one solute the other solute is added drop by drop until the appearance of turbidity occurs. The composition of single binodal points are used to draw the binodal curve but do not provide complete information on phase equilibrium. The determination of the compositions of the two conjugated phases requires more elaborate experiments. A ternary mixture belonging to the heterogeneous region is first stirred for a rather extended period and then allowed to separate in the two conjugated solutions; the time required for the separation can vary from few hours to 1 or 2 days for aqueous solutions of hydrocarbons. After standing, samples are taken from the individual phases and analyzed. The analytical determination can be carried out using a combination of various physical and chemical properties as refractive index, density, etc. or by direct determination of at least two components of the solution for example by means of gas chromatography. If the physical chemical properties are employed it is necessary the use of calibration curves. Usually the measurements required for the calibration curves can not be carried out at compositions close to the binodal curve and thus it is very often necessary to dilute the conjugated solutions to reach the composition range of validity of the calibration. Some other methods to determine the composition of the tie-line end points use additional turbidimetric titrations and eventually the knowledge of the binodal curve equation [28, 29]. As can be seen a much longer and more complicated procedure is necessary to obtain the composition of the conjugated solutions; as a consequence their compositions are usually affected by a larger error in respect to that of the single binodal ones. Finally the plait point is obtained with a graphical method starting from the values of the conjugated compositions. Two conjugated compositions are connected by the tie-line ; using this segment as a basis, two opposite triangles are constructed with the other sides parallel to those of the concentration diagram. This construction is repeated for each couple of conjugate compositions i.e. for each tie-line ; the curve interpolating the vertices of these triangles outside the tie-line intercepts the binodal in the plait point (ref. 25 p. 278). It is evident that because the error on the

6 Spinodal curve of water 1-butanol 1-propanol system. Wheeler-Widom model 1073 conjugate compositions and the interpolating procedure, even in presence of large m, the error on the composition of the plait point is not negligible. Ultimately, the ternary LLE that can be found in literature can be summarized as: { (x exp 1. N experimental binodal points b,i, ) } yexp b,i ; i=1,2,...,n { (x exp 2. m experimental tie-line end points (L),k, ) ( yexp (L),k, x exp (R),k, (R),k) } yexp k=1,2,...,m where the subscripts (L) and (R) indicates that one end belongs, in respect to the plait point (PP) position, to the left branch of the binodal curve and the other to the right one; 3. The experimental plait point, ( x exp P, ) yexp P. In facing the phase properties of a ternary system, it is useful to add to these data also the experimental binodal points relative to the binary system containing the two solutes. All these data of interest are reported for the system in exam in Figure 1. 3. THE LOCAL FITTING PROCEDURE For the system water (W) chloroform (CHCl 3 ) acetic acid (AcH) the attempt to fit all together the experimental binodal points with the theoretical binodal equation of the GWWHM by using different couples of the parameters (τ, ) failed because the agreement between the experimental and theoretical binodal was poor or just limited to some parts of the curve. Considering that this system presents characteristics of the binodal curve similar to that predicted by the original WW model it was easy to foresee that this behavior would have been general for almost all the systems. Therefore a different strategy was necessary. Because a better agreement between the experimental and theoretical binodal curve was obtained considering small range of the variable, a local fit was considered the winning strategy to face the problem. The Local Fitting Procedure can be summarized as follows: { (x exp 1. The set of the N experimental binodal points, b,i, ) } yexp b,i are fitted to a polynomial of suitable degree obtaining the representative experimental i=1,2,...,n binodal curve (REB) with y exp b (x) equation. Generally in this fitting the 2 m experimental tie-line end points and the plait point compositions have not been used because these data have larger errors in respect to the single binodal points. 2. For some selected points of REB, ( x b,j, y b,j = y exp b (x b,j ) ) indicated as fitting points F j, the equation yb,j th (x) of the local theoretical binodal curve (thbj) that obeys to a predefined constraint is calculated. The choice of this constraint characterizes the fitting procedures; more details about this step are reported below.

1074 Roberto Sartorio, Cristina Stoicescu 7 Fig. 2 Exemplification of PPLFM (Plait Point Local Fitting Method); representative experimental binodal (REB), P experimental plait point, F j fitting point, - - - - local theoretical binodal (thbj), P j local theoretical plait point, - - - - F j P, - P j P. This step requires the use of three parameters M b,j,τ j, j where M b,j is a composition parameter dependent on x b,j. 3. By using the (τ j, j ) parameters determined in step 2) the equation ys,j th (x) of the local theoretical spinodal curve (thsj) is calculated. 4. By using the (M b,j,τ j, j ) parameters the equation ytl,j th (x) of the local theoretical tie-line, (thtlj) is also obtained. 5. The local spinodal point ( ) ( ) x s,j,y s,j, that corresponds to xb,j,y b,j, is then found as the intersection point of the two curves ys,j th (x) and yth tl,j (x). In this procedure the parameter M s,j is also determined. 6. At last the spinodal curve of the system is obtained by fitting the single spinodal points to a polynomial of suitable degree. In step 2) for determining the parameters τ j, j and M b,j three independent equations are necessary. The first two are intuitive and impose that the theoretical

8 Spinodal curve of water 1-butanol 1-propanol system. Wheeler-Widom model 1075 binodal has to pass through the selected point of REB. That is: x th b,α (M b,j;τ j, j ) = x b,j, (1) yb,α th b,j;τ j, j ) = y b,j, (2) where α = L,R depending if the binodal point belongs to the Right or to the Left branch of the curve (GWWHM has different equations for the two branches of the binodal). Obviously there are infinite theoretical binodal curves passing through the fitting points F j but only one of these curves is necessary and sufficient to obtain the values of the three fitting parameters. Therefore it is possible to use different constraints to select a single theoretical local binodal curve; obviously it is convenient that the choice of the constraint is such that the resulting (thbj) has at least one property as similar as possible to the corresponding property of the (REB). Until now three different procedures to obtain the spinodal curve from the LLE data according to the GWWHM model have been worked out. The first was indicated as GTLLFM (General Tie Line Local Fitting Method) [9, 14] because the third equation of step 2) imposes that the slope of the theoretical tie-line passing through the fitting point has to be as close as possible to the slope of the experimental tie-line passing through the same point. The second one was indicated as GSLFM (General Slope Local Fitting Methods) [10]; for this procedure the third equation imposes that the local slope of (thbj), calculated in F j, has to be as close as possible to the slope of (REB) calculated in the same point. Both these procedures have given poor results for the systems HOH (2) HOC n (3) HOC 3 (1) investigated until now. For this reason a new procedure has been proposed. In this Plait Point Local Fitting Method (PPLFM) the third equation of step 2) imposes, for any fitting point ( ) x b,j,y b,j, that the distance between the experimental plait point P, ( x exp ) P,yexp P, and the local theoretical plait point Pj, ( x th P j,yp th ) j, must be minimum. Because at P j it results M j = 0 this implies: [x P P j = th b,α (0;τ j, j ) x exp ] 2 [ P + y th b,α (0;τ j, j ) y exp ] 2 P = min. (3) It is more convenient to express this condition in dimensionless form: to obtain this is we divide equation 3 for the distance between the fitting point F j, ( x b,j,y b,j ), and the experimental plait point P : P P j P F j = [x th b,α (0;τ j, j ) x exp P [xb,j x exp P ] 2 [ + y th b,α (0;τ j, j ) y exp ] 2 P ] 2 + [ yb,j y exp P The elements involved in this procedure are illustrated in Figure 2. ] 2 (4)

1076 Roberto Sartorio, Cristina Stoicescu 9 Fig. 3 M vs. x ; local theoretical M j, experimental binodal points of the binary system HOH (2) HOC 4 (3). 4. RESULTS For the system in exam the following LLE data are present in literature [15]: { (x exp 1. N = 27 experimental binodal points b,i, ) } yexp b,i ; i=1,2,...,n 2. m { = 3 experimental tie-line end points, defined by the x and y coordinates: (x exp (L),k, ) ( yexp (L),k, x exp (R),k, (R),k) } yexp where the subscripts (L) and (R) k=1,2,...,m indicates that one end belongs, in respect to the plait point (PP) position, to the left branch of the binodal curve and the other to the right one; 3. the experimental plait point, ( x exp P, yexp P 4. the experimental binodal compositions relative to the binary system HOH (2) HOC 4 (3) [31]. The experimental data are reported in Figure 1. In the PPLFM only the binodal point and the plait point compositions are necessary [30]. The N = 27 + 2 binodal points have been fitted to a forth degree polynomial imposing that the curve has to pass through the experimental plait point and the binary binodal points; the interpolating function, REB, is reported in the figure. ) ;

10 Spinodal curve of water 1-butanol 1-propanol system. Wheeler-Widom model 1077 Fig. 4 M vs. x ; local theoretical M j, experimental plait point (P), linear interpolation. In applying the PPLFM procedure we used a large set of fitting points that, because the plait point position is known, have a definite belonging either to the physical left branch or to the physical right branch of the binodal curve. In particular we have used 109 points for the left branch and 238 for the right branch of the binodal. For each ( x b,j,y b,j ) binodal fitting point we have, as result of the procedure, the values of parameters τ j, j, of the variable M j, and the composition of the corresponding spinodal point ( x s,j,y s,j ). We can also calculate the composition of the theoretical plait point P j that minimizes equation 4. Because of the use of the local fitting procedure the parameters τ j and j can change on the binodal composition and their variation on x has not any constrain and is not very significant. Differently the new variable M j has defined limits: in fact for each branch the following limits hold 0 M M max M = M (α) y (α) = 0 (binary limit). { M = 0 theoretical plait point; max While the M (α) max do not assume definite values the condition M = 0 at the experimental plait point must be fulfilled. We can use the respect of this condition to evaluate the reliability of our results. b

1078 Roberto Sartorio, Cristina Stoicescu 11 Fig. 5 Local theoretical plait point coordinates; representative experimental binodal (REB), experimental plait point (P), fitting point F j in the range 0.494 x Fj 0.468, corresponding P j in the range 0.468 x Pj 0.468, fitting point F j in the range 0.468 x Fj 0.429, corresponding P j in the range 0.471 x Pj 0.468. In Figure 3 the M j values are reported as a function of x; the definition range of M j is also indicated with two vertical line passing through the binodal points of the binary system HOH (2) HOC 4 (3). As it can be observed in the graph in the right branch of the curve there are at least two discontinuities; similar behavior is observed for the trend of the parameters τ j and j when reported versus x. However the curve is strictly decreasing in the left branch and strictly increasing in the right one as expected. In Figure 4 is reported a zoom of the previous graph for x values in the neighborhood of the plait point; the experimental composition ( x P,y P ) is also reported. Interpolating, with a straight line, separately the last nine points relative to the left branch and the first nine relative to the right one we obtain an excellent agreement between the x value that correspond for both the branches at M = 0, x = 0.4680, and the experimental composition of the plait point ( x P = 0.4679, y P = 0.900 ). It is worth noting that for the other two procedures the condition M = 0 at the experimental PP is not fulfilled. It is also interesting to analyze how the P j position changes as a function of x along the binodal curve. In Figures 5 and 6 the P j coordinates are reported for some selected set of fitting points of the binodal curve. In particular in Figure 5 the two

12 Spinodal curve of water 1-butanol 1-propanol system. Wheeler-Widom model 1079 Fig. 6 Local theoretical plait point coordinates; representative experimental binodal (REB), experimental plait point (P), fitting point F j in the range 0.300 x Fj 0.193, corresponding P j in the range 0.474 x Pj 0.454, fitting point F j in the range 0.050 x Fj 0.024 corresponding P j in the range 0.469 x Pj 0.468. selected portions of the REB belong one to the left branch and the other to the right branch of the curve and both include the experimental plait point. In Figure 6 the two selected portions belong to the right branch of REB. It is interesting to observe that, while for the left branch of binodal the local theoretical plait points are always located under the REB, for almost all the right branch, excluded a very limited portion of REB, the theoretical plait points are outside the binodal curve. Obviously there is not any constrain that imposes P j to be internal to the REB but we have never found for the other examined systems a similar behavior. This anomaly seems to be related to the high asymmetry of the binodal curve with the plait point very far from the maximum of the curve; this aspect will be analyzed in details in another paper. This unexpected behavior of the theoretical plait point also influences the spinodal point compositions. The spinodal ternary points are reported in Figures 7, 8 and 9 as well as the spinodal binary points of the system HOH (2) HOC 4 (3); these last data have been computed using the NRTL (Non Random Two Liquids) model and the corresponding constant present in literature [18].

1080 Roberto Sartorio, Cristina Stoicescu 13 Fig. 7 Local theoretical spinodal point coordinates, Left branch; - - - - representative experimental binodal (REB), experimental plait point (P), spinodal point of the binary system HOH (2) HOC 4 (3), ternary spinodal points, interpolated spinodal curve. To analyze the spinodal curve we consider three different ranges of x: the first corresponds to the left branch of REB, the second and the third to two consecutive ranges of the right branch. We can make the following observations: 1. Regarding the left branch of the binodal curve we observe that the corresponding spinodal points have a reasonable trend. They tend to the binodal curve approaching the plait point and moreover they are in a fairly good agreement with the binary spinodal point; the deviation observed at low concentration HOC 3 are probably due to the inability of the GWWHM to gives good results close to the binary limit. Then it is rather easy to draw this section of the spinodal curve just interpolating the spinodal points, the plait point, and the binary spinodal point, see Figure 7. 2. The first range of the right branch is defined by 0.4679 x 0.35 where x = 0.4679 corresponds to the plait point abscissa. In this range the spinodal points are outside the binodal curve, see Figure 8; this has no physical sense because the spinodal points must be by definition inside the binodal curve. It is easy to deduct that this unphysical result derives from the previously illustrated behavior of the theoretical plait point. Because all the obtained spinodal points

14 Spinodal curve of water 1-butanol 1-propanol system. Wheeler-Widom model 1081 Fig. 8 Local theoretical spinodal point coordinates, Right branch; - - - - representative experimental binodal (REB), experimental plait point (P), ternary spinodal points, interpolated spinodal curve. must be rejected as unphysical we could not draw the spinodal curve in this range of x; however looking at the general trend of the spinodal points it is reasonable to predict that in this range the spinodal curve is practically coincident with the binodal one. 3. In this range of composition the spinodal points are inside the binodal curve and then they are significant by physical point of view. However we remind that also these spinodal points, as those of the previous range, derive form theoretical spinodal curves whose plate points were outside the experimental binodal curve. This anomaly seems to influence the trend of the spinodal points; in fact they result not to be in good agreement with the binary spinodal point, see Figure 9. Ultimately the fitting of the spinodal points does not permit to draw a significant spinodal curve. We could guess a reasonable spinodal curve using, in the polynomial fitting, only the set of spinodal points close to the maximum and imposing that the polynomial has to pass through the binary binodal point.

1082 Roberto Sartorio, Cristina Stoicescu 15 Fig. 9 Local theoretical spinodal point coordinates, Right branch; - - - - representative experimental binodal (REB), spinodal point of the binary system HOH (2) HOC 4 (3), ternary spinodal points, interpolated spinodal curve. 5. CONCLUSIONS In this paper we have presented the result of the analysis of the LLE data for the system HOH (2) HOC 4 (3) HOC 3 (1) using the GWWHM. The binodal curve of the examined system is, in the triangle diagram, completely shifted towards the pure water corner accounting for the low solubility of 1-butanol in pure water and the much higher solubility of water in pure 1-butanol. The binodal curve is highly asymmetric and the tie lines have a high slope, see Figure 1. The PPLFM has given very good results for the left branch of the binodal, while, for the right branch, has given unphysical result for a long piece of the binodal close to the plait point and almost poor results for the remaining part. However, with the help of the spinodal points of the binary system system HOH (2) HOC 4 (3), a reasonable spinodal curve can be drawn, see Figure 10. If we compare this figure with Figure 11 where the binodal and spinodal curves are reported for the companion system HOH (2) HOC 5 (3) HOC 3 (1) [30] we observe some similarities. In both the systems the spinodal curve is practically coincident with the binodal one for a long section of the right branch just after the plait point; a preliminary

16 Spinodal curve of water 1-butanol 1-propanol system. Wheeler-Widom model 1083 Fig. 10 System HOH (2) HOC 4 (3) HOC 3 (1); - - - - representative experimental binodal (REB), experimental plait point (P), spinodal point of the binary system HOH (2) HOC 4 (3), ternary spinodal points, interpolated spinodal curve. analysis of the LLE data for the system HOH (2) HOC 6 (3) HOC 3 (1) confirms this behavior also for this system. The water rich metastable region is small while the alcohol rich metastable region is very large. This is very interesting because it should makes possible to run diffusion experiments in the right metastable region; it is worth reminding that the knowledge of the diffusion coefficients in these regions is a fundamental parameter in any phase separation model. Finally it is to observe that the y s,max is much closer to the y b,max in the case of system HOH (2) HOC 5 (3) HOC 3 (1) than in case of HOH (2) HOC 4 (3) HOC 3 (1); this trend seems to be confirmed by the preliminary analysis of the system HOH (2) HOC 6 (3) HOC 3 (1). In the next future this analysis will be extended to systems containing as solutes 1-alcohols with n = 7,8,9,10 to see the effect of much longer hydrophobic chains on the shape of the spinodal and on the extension of the metastable area.

1084 Roberto Sartorio, Cristina Stoicescu 17 Fig. 11 System HOH (2) HOC 5 (3) HOC 3 (1); - - - - representative experimental binodal (REB), experimental plait point (P), spinodal point of the binary system HOH (2) HOC 5 (3), ternary spinodal points, interpolated spinodal curve. Acknowledgements. We want to thank to Dr. Radu P. Lungu for his help to obtain the data corresponding to the GWWHM, and also for some important advices concerning the theoretical properties of this model. REFERENCES 1. J. C. Wheeler and B. Widom, Phase transitions and critical points in a model three-component system, J. Am. Chem. Soc. 90, 3064-3071 (1968). 2. D. A. Huckaby and M. Shinmi, Exact solution of a three-component system on the honeycomb lattice, J. Stat. Phys. 48, 135-144 (1986). 3. D. L. Strout, D. A. Huckaby and F. Y. Wu, An exactly solvable model ternary solution with threebody interactions, Physica A 173, 60-71 (1991). 4. F. D. Buzatu and D. A. Huckaby, An exactly solvable model ternary solution with three-body interactions, Physica A 299, 427-440 (2001). 5. F. D. Buzatu, D. Buzatu and J. G. Albright, Spinodal curve of a model ternary solution, J. Solution Chem. 30, 969-983 (2001). 6. F. D. Buzatu, R. P. Lungu and D. A. Huckaby, An exactly solvable model for a ternary solution with three-body interactions and orientationally dependent bonding, J. Chem. Phys. 121, 6195-6206 (2004).

18 Spinodal curve of water 1-butanol 1-propanol system. Wheeler-Widom model 1085 7. R. P. Lungu, D. A. Huckaby and F. D. Buzatu, Phase separation in an exactly solvable model binary solution with three-body interactions and intermolecular bonding, Phys. Rev. E 73, 021508 (14 pages) (2006). 8. R. P. Lungu and D. A. Huckaby, The microscopic structure of an exactly solvable model binary solution that exhibits two closed loops in the phase diagram, J. Chem. Phys. 129, 034505 (11 pages) (2008). 9. F. D. Buzatu, R. P. Lungu, D. Buzatu, R. Sartorio and L. Paduano, Spinodal composition of the system water + chloroform + acetic acid at 25 0 C, J. Solution Chem. 38 (4), 4003-415 (2009). 10. R. P. Lungu, R. Sartorio and F. D. Buzatu, New Method for Theoretical Spinodals Corresponding to Ternary Solutions with an Amphiphile Component, J. Solution Chem. 40 (10), 1687-1700 (2011). 11. R. P. Lungu, R. Sartorio and F. D. Buzatu, Theoretical spinodal for ternary solution with an amphiphile component using a generalized Wheeler-Widom model, Romanian Journal of Physics 56(9-10), 1069-1079 (2011). 12. V. Vitagliano, R. Sartorio, S. Scala and D. Spaduzzi, Diffusion in a ternary system and critical mixing point, J. Solution Chem. 7 (8), 605-622 (1978). 13. D. Buzatu, F. D. Buzatu, L. Paduano and R. Sartorio, Diffusion coefficients for the ternary system water + chloroform + acetic acid at 25 0 C, J. Solution Chem. 36 (11-12), 1373-1384 (2007). 14. D. Buzatu, F. D. Buzatu, R. P. Lungu, L. Paduano and R. Sartorio, On the determination of the spinodal curve for the system water + chloroform + acetic acid from the mutual diffusion coefficients Romanian Journal of Physics 55 (3-4), 342-351 (2010). 15. C. Stoicescu, O. Iulian and R. Isopescu, Liquid-Liquid Phase Equilibria of 1-Propanol + Water + n-alcohol Ternary Systems at 298.15 K and Atmospheric Pressure J. Chem. Eng. Data 56 (7), 3214-3221, (2011). 16. C. Stoicescu, O. Iulian and F. Sirbu, Liquid + liquid equilibrium data for the ternary mixtures of 1-propanol + water with 1-butanol, 1-hexanol, 1-octanol or 1-decanol at 294.15 K Rev. Roum. Chim. 53 (5), 363-367, (2008). 17. C. Stoicescu, O. Iulian and R. Isopescu, Liquid-Liquid Phase Equilibria of 1-Propanol + Water + n-alcohol Ternary Systems at 298.15 K and Atmospheric Pressure Rev. Roum. Chim. 56 (5), 553-560, (2011). 18. J. M. Sǿrensen and W. Arlt, Liquid-Liquid Equilibrium Data Collection, DECHEMA Chemistry Data Series, Vol. 1 Part 1 (1979), Frankfurt: DECHEMA. 19. J. M. Sǿrensen and W. Arlt, Liquid-Liquid Equilibrium Data Collection, DECHEMA Chemistry Data Series, Vol. 5 Parts 2,3 (1979), Frankfurt: DECHEMA. 20. C. Cascella, G. Castronuovo, V. Elia, R. Sartorio and S. Wurzburger, Hydrophobic interactions of alkanols. A calorimetric study in water at 25 0 C Journal of the Chemical Society, Faraday Transactions Trans. I 85 (10), 3289-3299 (1989). 21. S. Wurzburger, R. Sartorio, V. Elia and C. Cascella, Volumetric properties of aqueous solutions of alcohols and diols at 25 0 C Journal of the Chemical Society, Faraday Transactions 86 (23), 3891-3895 (1990). 22. L. Ambrosone, S. Andini, G. Castronuiovo, V. Elia and G. Guarino, Empirical correlations between thermodynamic and spectroscopic properties of aqueous solutions of aklkan-n,m diols. Excess enthalpies and spin lattice relaxation time at 298.15 K Journal of the Chemical Society, Faraday Transactions 87, 2989-2993 (1991). 23. G. Castronuovo, V. Elia and F. Velleca, On the role of the functional group in determining the strength of hydrophobic interactions of aqueous solutions of α amino acid derivatives J. Solution Chem. 25, 971-982 (1996). 24. G. Castronuovo, V. Elia and F. Velleca, Hydrophilic groups determine preferential configurations

1086 Roberto Sartorio, Cristina Stoicescu 19 in aqueous solutions. A calorimetric study of monocarboxylic acid and monoalkylamines at 298.15 K Thermochim, Acta 291, 21-26 (1997). 25. J. P. Novàk, J. Matou s and J. Pick, Liquid - Liquid Equilibria (Academia Prague 1987). 26. J. M. Prausnitz, R. N. Lichtenhaler and E. Gomez de Azevedo, Molecular Thermodynamics of Fluid Phase Equilibria (Prentice - Hall, Inc., New Jersey 1986). 27. P. S. Drzaic, Liquid Crystal Dispersion (World Scientific, Singapore 1995). 28. A. S. Smith, Analysis of ternary mixtures containing nonconsolute components, Ind. Eng Chem. 37 (2), 185-187 (1945). 29. D. M. T. Newsham and S. B. Ng Liquid-liquid equilibria I mixtures of water n-propanol and n- butanol J. Chem. Eng. Data 17 (2), 205-207 (1972). 30. R. P. Lungu and R. Sartorio Spinodal Curve for the System Water + 1-Pentanol + 1-Propanol at 25 0 C J. Solution Chem. 43 (1), 109-125 (2014). 31. J. M. Sǿrensen and W. Arlt, Liquid-Liquid Equilibrium Data Collection, DECHEMA Chemistry Data Series, Vol. 5 Parts 1 (1979), Frankfurt: DECHEMA.