Acessibilidade / Reportar erro

THERMODYNAMIC PROPERTIES AND PITZER PARAMETER DETERMINATION FOR ORTHOPHOSPHORIC ACID FROM FREEZING POINT AND ISOPIESTIC MEASUREMENT DATA

Abstract

The freezing depression points of aqueous phosphoric acid have been determined using a method based on measurement of solution conductivity. The osmotic coefficient for freezing point lowering was calculated by the condition of solid-liquid equilibrium and by the equation proposed by Bjerrum (1918). A calculation method is presented in this study to evaluate the ion interaction parameters for the Pitzer model from freezing point and isopiestic molalities of aqueous phosphoric acid solution. The model of Chen, Bromely and the Three-Characteristic Parameter Correlation (TCPC) are also used to model this system. The isopiestic molality values from literature data are used to optimize these models. Using these parameters, we were able to estimate thermodynamic properties: osmotic coefficient, mean activity coefficient and water activity for the H3PO4-H2O system, and contribution to a better understanding of the thermodynamic behaviour of the aqueous system containing phosphoric acid.

Keywords
Phosphoric acid; Pitzer parameters; Water activity; Activity coefficient; Osmotic coefficient

INTRODUCTION

Phosphoric acid is one of the strategic elements in the agriculture and pharmaceutical domain. In the literature, there are several different studies predicting all the parameters of the H3PO4-H2O system, although they were only realized at 298.15 K and done on water vapor pressure.

The thermodynamic properties of phosphoric acid are of crucial importance in industrial processes such as crystallization, reactors and chemistry reactions. The thermodynamic model proposed by Pitzer was used in this study to predict thermodynamic properties based on freezing point depression experimental results. In general, the use of the Pitzer model is based on the determination of parameters from experimental data, for example: the determination of water activity in the ternary system KI-KNO3-H2O at 298.15 K using an electronic hydrometer by Galleguillos et al. (2003)Galleguillos, H. R., Graber, T. A., and Taboada, M. E., Activity of Water in the KI + KNO3 + H2O Ternary System at 298.15 K, Journal of Chemical & Engineering Data, 48(4), 851-855 (2003)., determination of activity coefficients of electrolytes in the ternary system NaCl-Na2SO4-H2O from potential difference measurements at 298.15, 303.15 and 308.15 K by Sirbu et al. (2011)Sirbu, F., Iulian, O., Ion, A. C. and Ion, I. Activity Coefficients of Electrolytes in the NaCl + Na2SO4 + H2O Ternary System from Potential Difference Measurements at (298.15, 303.15, and 308.15) K, Journal of Chemical & Engineering Data, 56(12), 4935-4943 (2011). and osmotic and activity coefficients of dilute aqueous solutions of unsymmetrical tetraalkylammonium Iodides at 298.15 K using the isopiestic method by Amado-Gonzalez et al. (2012)Amado-Gonzalez, E. and Blanco, L. H., Osmotic and Activity Coefficients of Dilute Aqueous Solutions of Unsymmetrical Tetraalkylammonium Iodides at 298.15 K, Journal of Chemical & Engineering Data, 57(4), 1044(2012)..

There are also other thermodynamic models widely used to determine the thermodynamic properties of electrolytic solutions. Bromely (1973)Bromley, L.A., Thermodynamic properties of b electrolytes in aqueous solutions, AlChE Journal, 19(2), 313-320 (1973). presented a method for calculating the osmotic coefficient and the activity coefficient. He presented the values of the parameters used in this equation at 298.15 K to determine the properties of various electrolyte systems. This model has been used in several fields, for example in the work of Vicum et al. (2003)Vicum. L., Mazzotti, M. and Baldyga, J., Applying a Thermodynamic Model to the Non Stoichiometric Precipitation of Barium Sulfate, Chemical Engineering & Technology, 26(3), 325-333 (2003)., where the Bromley model was used for the non-stoichiometric precipitation of barium sulfate, BaSO4. Chen et al. (1982)Chen, C. C., Britt, H. I., Boston, J. F., and Evans, L. B., Local composition model for excess Gibbs energy of electrolyte systems. Part I: Single solvent, single completely dissociated electrolyte systems. AIChE Journal, 28(4), 588-596 (1982). proposed a model to calculate the activity coefficient with long-range and short-range contributions. This model successfully correlated experimental data at 298.15 K, for example in the work of Taghi et al. (2003), who applied the Chen model to correlate experimental data of the osmotic coefficient and density of non-aqueous systems. Recently, the Three-Characteristic Parameter Correlation (TCPC) model proposed by Ge et al. (2007)Ge, X., Wang, X., Zhang, M., and Seetharaman, S., Correlation and Prediction of Activity and Osmotic Coefficients of Aqueous Electrolytes at 298.15 K by the Modified TCPC Model, Journal of Chemical & Engineering Data, 52(2), 538-547 (2007). was used to correlate experimental data of several aqueous and non-aqueous systems.

In the case of aqueous phosphoric acid solution, the work of Pitzer et al. (1976)Pitzer, K. S. and Silvester, L. F., Thermodynamics of Geothermal Brines I. Thermodynamic Properties of Vapor-Saturated NaCl(aq) Solutions From 0 - 300ºC, Lawrence Berkeley Laboratory Report, (1976). was an important effort to predict the thermodynamic parameters of phosphoric acid at 298.15 K up to 6 mol.kg-1 based on the statistical-mechanical equations. Their work employed isopiestic results coupled with electromotive force (EMF). Holmes et al. (1999)Holmes, H. F., and Mesmer, R. E., Isopiestic Studies of H3PO4(aq) at Elevated Temperatures, Journal of Solution Chemistry, 28(4), 327-340 (1999). determined experimentally the isopiestic molalities of aqueous solutions of phosphoric acid at high temperature in the range of 383.15 to 523.15 K. The recent study proposed by Yang et al. (2016)Yang, H., Zhao, Z., Zeng, D., and Yin, R., Isopiestic Measurements of Water Activity for the H2SO4-H3PO4-H2O System at 298.15 K, Journal of Solution Chemistry, 45(11), 1580-1587 (2016). presents experimental values of water activity at 298.15 K using the isopiestic method.

The objective of this work is to present new experimental values of freezing depression points up to 241.15 K by using the conductimetric method, thus combining the freezing point and isopiestic data to predict more general and precise thermodynamic properties for the phosphoric acid-water system up to 8 mol.kg-1, with a suggestion for adapting Pitzer's model to any temperature of phosphoric acid solutions, with comparison to the Chen model, Bromley model and Three-Characteristic Parameter Correlation model (TCPC).

EXPERIMENTAL METHODS

In the present study, the freezing point depression was determined experimentally by using a method based on measurements of conductivity, λ. In this method, the ice solubility curve is measured at fixed temperature by adding step-by-step phosphoric acid with a known concentration. The saturation by ice is achieved, which generates a variation in the slope of the curve of solution conductivity versus concentration. At this point the freezing temperature and molality are determined (Goundali et al., 2008Goundali, B. E. and Kaddami, M., The ternary system: H2O-Fe(NO3)3-Co(NO3)2 isotherms −15 and −25ºC, Journal of Alloys and Compounds, 460(1-2), 544-548(2008).; Goundali et al., 2007Goundali, B. E. and Kaddami, M., The ternary system H2O-Fe(NO3)3-Co(NO3)2 isotherms 0 and 15 ºC, Fluid Phase Equilibria, 260(2), 295-299 (2007).). For example, to determine freezing depression points of the system H2O-A, we start with the sample M, Figure 1-B. After each addition, the conductivity of the solution is measured until stability of its value. The curve (A) in Figure 1 (conductivity=f(volume of water)) obtained for the binary system H2O-A presents a break at phase change when an invariant equilibrium is analyzed. The operation is repeated until the total determination of the liquidus of the ice.

Figure 1
Conductivity study variation of the isothermal dilution of mixture M in the binary system H2O-A

Figure 2 shows the experimental setup used for the FPD determinations using the conductimetric method.

Figure 2
Experimental setup, (1) refrigerator; (2) conductimetric measuring zone; (3) temperature display; (4) thermometer; (5) conductivity meter; (6) capillary water-inlet; (7) thermostatic jacket; (8) aqueous solution; (9) magnetic bar; (10) closure cap.

RESULTS

Table 1 presents the ice saturation curve obtained for the H3PO4-H2O system in the range of temperature from 241.15 to 268.15 K with an estimated standard temperature uncertainty of ±0.1 K. The temperature control was performed by a VWR Model 1197P thermostat (USA). Water was distilled twice with conductivity on the order of 6 µS at 291.15 K. Two products were used to ensure the reliability of the experimental results, phosphoric acid with 85% concentration in water that presents 99.99% purity and phosphorus pentoxide with purity 99%. (Sigma-Aldrich) The concentration of phosphoric acid was verified using titration with sodium hydroxide.

Table 1
Experimental results for freezing point depression.

In Figure 3, our experimental freezing depression values and those of Lide (2003-2004) and Jones et al. (1907)Jones, H. C., Getman, F. H., Bassett, H. P. McMaster, L. and Uhler, H. S, Hydrates in Aqueous Solution: Evidence for the Existence of Hydrates in Solution. Their Approximate Composition and Certain Spectroscopic investigation bearing upon the hydrate problem (No. 60), Carnegie Institution of Washington: Baltimore, MD, (1907). are presented and compared. This figure shows that the data of Jones et al. (1907)Jones, H. C., Getman, F. H., Bassett, H. P. McMaster, L. and Uhler, H. S, Hydrates in Aqueous Solution: Evidence for the Existence of Hydrates in Solution. Their Approximate Composition and Certain Spectroscopic investigation bearing upon the hydrate problem (No. 60), Carnegie Institution of Washington: Baltimore, MD, (1907). have a large deviation and below 267.6 K the difference between the values is significant.

Figure 3
Experimental data for phosphoric acid freezing point depression.

The isopiestic molalities of aqueous phosphoric acid solutions have been measured by Holmes et al. (1999)Holmes, H. F., and Mesmer, R. E., Isopiestic Studies of H3PO4(aq) at Elevated Temperatures, Journal of Solution Chemistry, 28(4), 327-340 (1999). using NaCl for the reference solution in the range of temperature from 383.15 to 523.15K for molality up to 10 mol.kg-1. Table 2 presents the number of points used to model phosphoric acid's parameters:

Table 2
Number of data points used

In the presence of phosphoric acid in solution, the equilibrium between ice and water liquid is described by:

(1) H 2 O ( solid ) H 2 O ( liquid )

Below the temperature 273.15 K the condition for the equilibrium is:

(2) µ H 2 O l = µ H 2 O * s = µ ( T ) * l + RT ln ( a w )

In eq 2, µ designates the chemical potential and the asterisk (*) the pure phase, aw represents the activity of water. Following eq. 3, one can calculate directly values of water activity from the experimental data of freezing point based on the Gibbs-Duhem equation:

(3) ln ( a w ) = H * R 1 T f 1 273 . 15 C p R T f 273 . 15 T f ln T f 273 . 15

In the present calculations, the following values were used in this equation: the enthalpy of fusion for pure water ΔH*=6009.5J.mol1, the universal gas constant R=8.314J.mol1.K1 and ΔCp=37.87J.mol1.K (Osborne, 1939Osborne, N. S., Stimson, H. F., and Ginnings, D. C., Measurements of heat capacity and heat of vaporization of water in the range 0 to 100 ºC, J. Res. Natl. Bur. Stand, 23, 197-260 (1939).; Osborne et al., 1939Osborne, N. S., Heat of fusion of ice. A revision, J. Res. Natl. Bur. Stand, 23, 643-646 (1939).). The heat capacity value was used by Hassan et al. (2014)Hasan, M., Partanen, J. I., Vahteristo, K. P., and Louhi-kultanen, M., Determination of the Pitzer Interaction Parameters at 273.15 K from the Freezing-Point Data Available for NaCl and KCl Solutions, Industrial & Engineering Chemistry Research, 53(13), 5608-5616 (2014). to calculate Pitzer's parameters for NaCl and KCl from the freezing point and by Arrad et al. (2015)Arrad, M., Kaddami, M., Sippola, H., and Taskinen, P., Determination of Pitzer Parameters for Ferric Nitrate from Freezing Point and Solubility of Ice, Journal of Chemical & Engineering Data, 61(1), 674-678 (2015). in the case of ferric nitrate up to the temperature 248.15 K.

The osmotic coefficient ϕ is defined by the ratio πrealπideal 4, where π is the osmotic pressure:

(4) π = RT ϑ S ln ( a S )

where ϑs is the molar volume of solvent. Therefore, the expression of the osmotic coefficient becomes:

(5) ϕ = π ( real ) π ( ideal ) = ln ( a S real ) ln ( a S ideal )

In the case of an electrolyte solution, where one molecule of acid MX dissociates into ν ions, the ideal solution activity of the solvent (water in this case) is:

(6) ln a W ideal = M S 1000 Ѵ m MX

By using equations 5 and 6, the osmotic coefficient equation in the real solution becomes:

(7) ϕ = 1000 M S Ѵ m MX ln a W real

The activity coefficient for a real solution is calculated using Equation 3.

The other method to calculate the osmotic coefficient from the freezing point lowering was presented by Randal (1923)Randall, M., Methods of calculation of activity coefficients, Transactions of the Faraday Society, 23, 502-507 (1927).:

(8) j = 1 θ Ѵ δ m

A quantity ϕ identical in value with (1-j) was called the osmotic coefficient by Bjerrum (1918)Bjerrum, N., Die Dissoziation der starken Elektrolyte, Z. Electrochem., 24, 321-328 (1918).; the function j is a property of the solvent:

(9) ϕ = θ ѵ δ m

where ν is the number of molecules formed per molecule of solute, θ is the freezing point lowering, δ is a constant =1.858 (Randall, 1923) and m is the molality.

The application of equation 7 in our study is represented by:

(10) ϕ = 1000 v . M H 20 i = 1 n m i ln ( a w )

where mi is molality of species [mol.kg-1] in the solution, the molar mass of water MH2O=18.01528g.mol1 with ν=ν++ν. The results of application of equations 9 and 10 are presented in Table 1.

THERMODYNAMIC MODELING

Pitzer model

The Pitzer model was applied for the treatment of non-ideality and in the case of 1-1 electrolyte the osmotic coefficient ϕ of the solution and mean activity coefficient γ± of the electrolyte can be written as:

(11) ln γ MX ± = Z + Z f γ + m 2 v + v v B MX γ + m 2 2 v + v 3 2 v C MX γ

(12) ϕ 1 = Z + Z f ϕ + m 2 v + v v B MX ϕ + m 2 2 v + v 3 2 v C MX ϕ

where:

(13) f γ = A D H I 1 + b I + 2 b ln 1 + b I

(14) B MX γ = 2 β MX 0 + 2 β MX 1 a 2 I 1 1 + α I α 2 I 2 e a I

(15) C MX γ = 3 2 C MX ϕ

(16) f ϕ = A D H I 1 + b I

(17) B MX ϕ = 2 β MX 0 + 2 β MX 1 α 2 I e α I

(18) A D H = 1 3 2 π N 0 d 1 e 2 4 ε DKT 3 2

where ν+ and ν- are the numbers of M and X ions in the formula, Z+ and Z- are their respective charges, m is the conventional molality, AD-H is the Debye-Hückel constant, d1 is water density and ε is permittivity of vacuum, D the dielectric permeability of pure water at temperature T [K], N0 is Avogadro's number, K is the Boltzmann constant, the value of AD-H is 0.392 kg1/2.mol-1/2 at 298.15 K.

The Pitzer values of two parameters a and bare 2.0 and 1.2, respectively. I is the stoichiometric ionic strength of the solution [mol/kg].

Bromley model

The equation of the Bromley model for the calculation of the osmotic coefficient can be written as:

(19) 1 ϕ = 2 . 303 A γ Z + Z I 3 δ ρ I 2 . 303 0 . 06 + 0 . 6 B T Z + Z I 2 ψ aI 2 . 303 B ( T ) I 2

where:

(20) a = 1 . 5 Z + Z

(21) δ ρ I = 3 ρ I 1 + ρ I 1 1 + ρ I 2 ln 1 + ρ I

(22) ψ aI = 2 aI 1 + 2 aI 1 + aI 2 ln 1 + aI aI

where Aγ is the Debye-Hückel constant on a log e basis at temperature T [K], ρ is a constant related to the distance of closest approach; in this case the value of ρ =1.0, Iis the ionic strength and B(T) is Bromley's parameter.

Chen model

Chen et al. (1982)Chen, C. C., Britt, H. I., Boston, J. F., and Evans, L. B., Local composition model for excess Gibbs energy of electrolyte systems. Part I: Single solvent, single completely dissociated electrolyte systems. AIChE Journal, 28(4), 588-596 (1982). proposed a model for calculating the activity coefficient based on long-range and short-range contributions:

(23) ln f ± = ln f ± LR + ln f ± SR

where f± is the activity coefficient based on mole fraction, LR and SR refer to the long-range and short-range contributions.

The equation for short-range can be written as:

(24) ln f ± LR = 1 v v C ln f + LR + v a ln f LR

with:

(25) ln f + LR = x m 2 τ ca , m F cm x c F cm + x a F cm + x m 2 Z x a τ m , ca x m F ma x a + x m F ma 2 + Z + x m τ m , ca F mc x a + x m F ma Z + τ m , ca F cm τ ca , m

(26) ln f LR = x m 2 τ ca , m F am x c F cm + x a F am + x m 2 Z + x c τ m , ca x m F mc x a + x m F ma 2 + Z x m τ m , ca F ma x c + x m F ma Z τ m , ca F am τ ca , m

The mole fraction xi as a function of molality m is written as:

(27) x i = v i m vm + 1000 M H 2 O

where xc, xm, xa: are cation, solvent and anion mole fraction, τm,ca and τca,m are Chen parameters, Fmc=Fma=expατm,ca and Fcm=Fam=expατca,m.

The equation for long-range can be written as:

(28) ln f ± LR = 1000 M H 20 A D H 2 Z + Z ρ ln 1 + ρ I x + Z + Z I x 2 I x 3 2 1 + ρ I x

To make the data uniform, the mean activity coefficient is converted to the molality form using the relationship:

(29) ln ( γ MX ± ) = ln f ± ln 1 + 0 . 001 M H 20 vm

In our case, the determination of the Chen parameters is adjusted with the osmotic coefficient data. Pitzer et al. (1986)Pitzer, K.S. and Simonsont, J.M., Thermodynamics of Multicomponent, Miscible, Ionic Systems: Theory and Equations, Journal of Physical Chemistry, 90(13), 3005-3009 (1986). presented the long-range solvent equation:

(30) ln γ solvent PDH = 2 A x I x 3 / 2 1 + ρ I x 3 / 2

The short-rang contribution of solvent according to the Chen model is:

(31) ln f S SR = x cm τ ca , m + x am τ ca , m + Z + x c τ m , ca x a F mc x a + x m F ma 2 + Z x a τ m , ca x m F am x c x a + x m F ma 2 x c x m F cm τ ca , m x c F cm + x a F cm + x m 2 x a x m F am τ ca , m x c F cm + x a F cm + x m 2

with:

(32) x ij = x i F ij x a F aj + x c F cj + x m F jj

According to equations 30 and 31, the osmotic coefficient can be calculated by equation 10, using the expression lnas=lnγs+lnxs.

Three-Characteristic Parameter Correlation model

The following equation was proposed for calculation of the osmotic coefficient using the Three-Characteristic Parameter Correlation model:

(33) ln γ ± = Z + Z A D H I 1 + b I + 2 b ln 1 + b I + S T I 2 n v + v

The osmotic Coefficient can be calculated by integrating the equation for the mean activity coefficient:

(34) ϕ = 1 + 1 m 0 m m d ln γ ±

From eqs. 33 and 34, the osmotic coefficient can be written as:

(35) ϕ = 1 Z + Z A D H I 1 + b I + S T v + v 2 n 2 n 1 I 2 n

In this equation, b, S and n are the three characteristic parameters: S is the solvation parameter, b is the approach parameter representing the closest distance between ions and n is a distance parameter related to the distance between ion and solvent molecule. Using these parameters the osmotic coefficients of electrolytes in aqueous solution could be estimated.

For all thermodynamic models, the ionic strength of the solution [mol/kg], is given as:

(36) I = 1 2 i = 1 n m i Z i 2

The adjustable parameters βMX0, βMX1 and CMXϕ for Pitzer model, τm,ca and τca,m for the Chen model, B(T)for the Bromely model and b, S, n for the Three-Characteristic Parameter Correlation model are specific for the phosphoric acid-water system. They have been obtained from least-square fitting of the freezing point and isopiestic data.

RESULTS AND DISCUSSION

The Pitzer parameters for phosphoric acid-water were determined by using different types of data: solubility and isopiestic molalities. Pitzer's most extensive investigation was done for NaCl. As shown in the works of Sippola (2013)Sippola, H., Critical Evaluation of the Second Dissociation Constants for Aqueous Sulfuric Acid over a Wide Temperature Range, Journal of Chemical & Engineering Data, 58(11), 3009-3032 (2013). and Sippola et al. (2014)Sippola H., and Taskinen P., Thermodynamic Properties of Aqueous Sulfuric Acid, Journal of Chemical & Engineering Data, 59(8), 2389-2407 (2014) the variation of βMX0, βMX1 and CMXϕ was tested using a temperature dependency in the form =a+bT. In our work the variation of Pitzer's parameter values for phosphoric acid resulted in the following equations:

(37) β MX 0 = a 0 + b 0 T + c 0 ln ( T )

(38) β MX 1 = a 1 + b 1 T + c 1 T 2

(39) C MX 1 = a ϕ + b ϕ T + c ϕ ln ( T )

The parameter variation of the Bromley model as a function of the temperature is:

(40) B ( T ) = B 1 T 230 + B 2 T + B 3 + B 4 ln ( T )

The variations of the Chen parameters as a function of the temperature are:

(41) τ m , ca = a m , ca + b m , ca 1 T 1 298 . 15 + C m , ca 298 . 15 T T + ln T 298 + 15

(42) τ ca , m = a ca , m + b ca , m 1 T 1 298 . 15 + C ca , m 298 . 15 T T + ln T 298 . 15

In the case of the Three-Characteristic Parameter Correlation model nis a constant, the variation of S and b as a function of temperature is:

(43) S T = a S + b S T + C S ln T

(44) b T = a b + b b T + C b ln T

The parameter values were determined using the least-squares method with the objective function of eq 45 that minimizes the errors between the osmotic coefficients from the thermodynamic models and the experimental values:

(45) Ob = i = 1 n ϕ exp ϕ calc 2

The Deby-Hückel constant AD-H in the range of temperature 230 K to 298.15 K was calculated from the polynomial based on the tabulated values of Deby-Hückel constant presented by Pitzer et al. (1976)Pitzer, K. S. and Silvester, L. F., Thermodynamics of Geothermal Brines I. Thermodynamic Properties of Vapor-Saturated NaCl(aq) Solutions From 0 - 300ºC, Lawrence Berkeley Laboratory Report, (1976). as:

(46) A D H = 61 . 44534 exp T 273 . 15 273 . 15 + 2 . 864468 exp T 273 . 15 273 . 15 2 + 183 . 5379 ln T 273 . 15 0 . 6820223 T 273 . 15 + 0 . 0007875695 T 2 273 . 15 2 + 58 . 95788 273 . 15 T

For the range 298.15 to 523.5 K, the Deby-Hückel constant is given in the work of Holmes et al. (1999)Holmes, H. F., and Mesmer, R. E., Isopiestic Studies of H3PO4(aq) at Elevated Temperatures, Journal of Solution Chemistry, 28(4), 327-340 (1999)..

The standard deviation, σ, was defined as:

(47) σ T = i = 1 n ϕ exp ϕ calc 2 NP 1 / 2

where NP represents the number of experimental points at each temperature.

The values of ion-interaction parameters for the thermodynamic models are presented in Table 3:

Table 3
Values of thermodynamic models parameters for phosphoric acid in the range 241.15 to 523.15 K

The variation of βMX1 for the Pitzer model depends on a 1 and b 1 because c 1 is very small:

(48) β MX 1 = a 1 + b 1 T

Following eq 47, we calculated the standard deviation for each model. Table 4 presents the values for every temperature:

Table 4
Values of the standard deviation at high temperature

As can be seen from Table 4, the values of standard deviations for fitting the experimental data of osmotic coefficient for the Pitzer model are very low compared to the other models.

Figure 4 presents the objective function for the global treatment of each model.

Figure 4
Objective function for each model for the phosphoric acid-water system.

The comparison between experimental results of freezing depression curve and thermodynamic models is presented in Figure 5. According to Figure 5, the Pitzer model accurately predicts the experimental results, which is why their objective function value remains the lowest (Figure 4), whereas the Chen model diverges after 7.2 mol /kg.

The Pitzer model presents the lowest value of objective function in the case of the phosphoric acid-water system and was used to calculate the activity coefficient of water, osmotic coefficient and mean activity coefficient at different temperatures.

Figure 5 shows the calculated osmotic coefficient for the phosphoric acid-water system at different temperatures using the Pitzer parameters presented in Table 3.

Figure 5
A: Comparison between experimental results and thermodynamic models to predict the freezing depression curve; B: Calculated osmotic coefficient vs. ionic strength at 413.14 to 498.15 K using Pitzer model: acalculated osmotic coefficient by Pitzer model in this study, bexperimental values of osmotic coefficient.

It is important to calculate the mean activity coefficient γH3PO4± at temperatures below 273.15 K. Table 5 presents mean activity coefficient values at selected molalities:

Table 5
Mean activity coefficient at low temperatures using Pitzer model.

The mean coefficient calculated in Table 5 decreases as a function of temperature in high molality, but the values are almost equal in low molality. This confirms results presented earlier by Holmes et al. (1999)Holmes, H. F., and Mesmer, R. E., Isopiestic Studies of H3PO4(aq) at Elevated Temperatures, Journal of Solution Chemistry, 28(4), 327-340 (1999)..

Calculated water activity coefficients are usually reported in terms of the ionic strength. Figure 6 presents the values of aw at 298.15K calculated by equations 11 and 10 using the values of Pitzer parameters presented in Table 3:

Figure 6
Activity coefficient for water in the H3PO4-H2O system up to 8 mol.kg-1 at 298.15 K.

The water activity coefficients determined at 298.15 in this study using the values of Pitzer's model, Table 3, predict exactly the same values determined experimentally by Platford (1975)Platford, R. F., Thermodynamics of Aqueous Solutions of Orthophosphoric Acid from the Freezing Point to 298.15ºK, Journal of Solution Chemistry, 4(7), 591-598 (1975)., the isopiestic measurements of water activity in recent study by Yang et al. (2016)Yang, H., Zhao, Z., Zeng, D., and Yin, R., Isopiestic Measurements of Water Activity for the H2SO4-H3PO4-H2O System at 298.15 K, Journal of Solution Chemistry, 45(11), 1580-1587 (2016). and results presented by Elmore et al. (1946)Elmore, K. L., Mason, C. M., and Christensen, J. H., Activity of Orthophosphoric Acid in Aqueous Solution at 25º from Vapor Pressure Measurements, Journal of the American Chemical Society, 68(12), 2528-2532 (1946)., which confirms the Pitzer parameters presented in Table 3.

CONCLUSION

The present paper reports the results of an experimental freezing point study of aqueous phosphoric acid solutions using a conductimetric method. The model used in this work, that of Pitzer, was compared with 222 points of thermodynamic properties: activity coefficient, mean activity coefficient and osmotic coefficient in the range of temperature from 241.15 to 523.15 K up to 8 mol.kg-1 for H3PO4-H2O system. The experimental values related to water activity from the literature confirm the values calculated by the Pitzer model based on the experimental data of freezing point depression, which shows that the use of freezing point depression data for the modeling of electrolyte systems is very interesting.

The form of variation of the Pitzer parameters βMX0, βMX1 and CMX presented in this work produces better results compared to the Chen model, Bromley model and the Three-Characteristic Parameter Correlation model. According to the equations presented in this paper there are two approaches to calculate the osmotic coefficient value from freezing point depression. The first approach is based on the thermodynamic equilibrium condition, whereas the second approach is by direct calculation of the osmotic coefficient, based on the use of the equation proposed by Bjerrum (1918)Bjerrum, N., Die Dissoziation der starken Elektrolyte, Z. Electrochem., 24, 321-328 (1918).. Both methods resulted in almost the same values with the difference in the third digit after the decimal point, which shows that the Bjerrum equation is recommended to deduce the osmotic coefficient values from the freezing point depression.

NOMENCLAUTURE

  • aw  Activity of water
  • AD-H  Debye-Hückel parameter (kg12.mol12)
  • m  Molality (mol/kg)
  • I  Ionic strength of electrolyte
  • M  Molecular weight (g/mol)
  • T  Temperature (K)
  • FPD  Freezing point depression
  • TCPC  Three-Characteristic Parameter Correlation model
  • b  Debye-HUCKEL parameter(kg12.mol12)
    Greek letter
  • µ  Chemical potential
  • γ MX ±  Activity coefficient of element MX.
  • ϕ  Osmotic coefficient.
  • α  Universal parameter (kg12.mol12)
    Sub/Superscripts
  • f  freezing point
  • *  pure chemical compound

REFERENCES

  • Arrad, M., Kaddami, M., Sippola, H., and Taskinen, P., Determination of Pitzer Parameters for Ferric Nitrate from Freezing Point and Solubility of Ice, Journal of Chemical & Engineering Data, 61(1), 674-678 (2015).
  • Bjerrum, N., Die Dissoziation der starken Elektrolyte, Z. Electrochem., 24, 321-328 (1918).
  • Bromley, L.A., Thermodynamic properties of b electrolytes in aqueous solutions, AlChE Journal, 19(2), 313-320 (1973).
  • Chen, C. C., Britt, H. I., Boston, J. F., and Evans, L. B., Local composition model for excess Gibbs energy of electrolyte systems. Part I: Single solvent, single completely dissociated electrolyte systems. AIChE Journal, 28(4), 588-596 (1982).
  • Lide, D.R., CRC, Handbook of Chemistry and Physics, (2003-2004).
  • Goundali, B. E. and Kaddami, M., The ternary system H2O-Fe(NO3)3-Co(NO3)2 isotherms 0 and 15 ºC, Fluid Phase Equilibria, 260(2), 295-299 (2007).
  • Goundali, B. E. and Kaddami, M., The ternary system: H2O-Fe(NO3)3-Co(NO3)2 isotherms −15 and −25ºC, Journal of Alloys and Compounds, 460(1-2), 544-548(2008).
  • Amado-Gonzalez, E. and Blanco, L. H., Osmotic and Activity Coefficients of Dilute Aqueous Solutions of Unsymmetrical Tetraalkylammonium Iodides at 298.15 K, Journal of Chemical & Engineering Data, 57(4), 1044(2012).
  • Elmore, K. L., Mason, C. M., and Christensen, J. H., Activity of Orthophosphoric Acid in Aqueous Solution at 25º from Vapor Pressure Measurements, Journal of the American Chemical Society, 68(12), 2528-2532 (1946).
  • Sirbu, F., Iulian, O., Ion, A. C. and Ion, I. Activity Coefficients of Electrolytes in the NaCl + Na2SO4 + H2O Ternary System from Potential Difference Measurements at (298.15, 303.15, and 308.15) K, Journal of Chemical & Engineering Data, 56(12), 4935-4943 (2011).
  • Hasan, M., Partanen, J. I., Vahteristo, K. P., and Louhi-kultanen, M., Determination of the Pitzer Interaction Parameters at 273.15 K from the Freezing-Point Data Available for NaCl and KCl Solutions, Industrial & Engineering Chemistry Research, 53(13), 5608-5616 (2014).
  • Galleguillos, H. R., Graber, T. A., and Taboada, M. E., Activity of Water in the KI + KNO3 + H2O Ternary System at 298.15 K, Journal of Chemical & Engineering Data, 48(4), 851-855 (2003).
  • Holmes, H. F., and Mesmer, R. E., Isopiestic Studies of H3PO4(aq) at Elevated Temperatures, Journal of Solution Chemistry, 28(4), 327-340 (1999).
  • Jones, H. C., Getman, F. H., Bassett, H. P. McMaster, L. and Uhler, H. S, Hydrates in Aqueous Solution: Evidence for the Existence of Hydrates in Solution. Their Approximate Composition and Certain Spectroscopic investigation bearing upon the hydrate problem (No. 60), Carnegie Institution of Washington: Baltimore, MD, (1907).
  • Pitzer, K.S. and Simonsont, J.M., Thermodynamics of Multicomponent, Miscible, Ionic Systems: Theory and Equations, Journal of Physical Chemistry, 90(13), 3005-3009 (1986).
  • Vicum. L., Mazzotti, M. and Baldyga, J., Applying a Thermodynamic Model to the Non Stoichiometric Precipitation of Barium Sulfate, Chemical Engineering & Technology, 26(3), 325-333 (2003).
  • Zafarani-Moattar, M.T. and Sardroodi, J.J., Isopiestic and Volumetric Study of CuCl2 + Ethanol at 25 ºC, Journal of Chemical & Engineering Data, 48(2), 308-313 (2003).
  • Osborne, N. S., Heat of fusion of ice. A revision, J. Res. Natl. Bur. Stand, 23, 643-646 (1939).
  • Osborne, N. S., Stimson, H. F., and Ginnings, D. C., Measurements of heat capacity and heat of vaporization of water in the range 0 to 100 ºC, J. Res. Natl. Bur. Stand, 23, 197-260 (1939).
  • Pitzer, K. S., and Silvester, L. F., Thermodynamics of electrolytes. VI. Weak electrolytes including H3PO4, Journal of Solution Chemistry, 5, 269-277 (1976).
  • Pitzer, K. S. and Silvester, L. F., Thermodynamics of Geothermal Brines I. Thermodynamic Properties of Vapor-Saturated NaCl(aq) Solutions From 0 - 300ºC, Lawrence Berkeley Laboratory Report, (1976).
  • Platford, R. F., Thermodynamics of Aqueous Solutions of Orthophosphoric Acid from the Freezing Point to 298.15ºK, Journal of Solution Chemistry, 4(7), 591-598 (1975).
  • Randall, M., Methods of calculation of activity coefficients, Transactions of the Faraday Society, 23, 502-507 (1927).
  • Sippola H., and Taskinen P., Thermodynamic Properties of Aqueous Sulfuric Acid, Journal of Chemical & Engineering Data, 59(8), 2389-2407 (2014)
  • Sippola, H., Critical Evaluation of the Second Dissociation Constants for Aqueous Sulfuric Acid over a Wide Temperature Range, Journal of Chemical & Engineering Data, 58(11), 3009-3032 (2013).
  • Ge, X., Wang, X., Zhang, M., and Seetharaman, S., Correlation and Prediction of Activity and Osmotic Coefficients of Aqueous Electrolytes at 298.15 K by the Modified TCPC Model, Journal of Chemical & Engineering Data, 52(2), 538-547 (2007).
  • Yang, H., Zhao, Z., Zeng, D., and Yin, R., Isopiestic Measurements of Water Activity for the H2SO4-H3PO4-H2O System at 298.15 K, Journal of Solution Chemistry, 45(11), 1580-1587 (2016).

Publication Dates

  • Publication in this collection
    Jul-Sep 2018

History

  • Received
    14 Aug 2017
  • Reviewed
    04 Oct 2017
  • Accepted
    05 Oct 2017
Brazilian Society of Chemical Engineering Rua Líbero Badaró, 152 , 11. and., 01008-903 São Paulo SP Brazil, Tel.: +55 11 3107-8747, Fax.: +55 11 3104-4649, Fax: +55 11 3104-4649 - São Paulo - SP - Brazil
E-mail: rgiudici@usp.br