UNIQUAC

Last updated February 29, 2024

In statistical thermodynamics, UNIQUAC (a portmanteau of universal quasichemical) is an activity coefficient model used in description of phase equilibria.^[1]^[2] The model is a so-called lattice model and has been derived from a first order approximation of interacting molecule surfaces. The model is, however, not fully thermodynamically consistent due to its two-liquid mixture approach.^[2] In this approach the local concentration around one central molecule is assumed to be independent from the local composition around another type of molecule.

The UNIQUAC model can be considered a second generation activity coefficient because its expression for the excess Gibbs energy consists of an entropy term in addition to an enthalpy term. Earlier activity coefficient models such as the Wilson equation and the non-random two-liquid model (NRTL model) only consist of enthalpy terms.

Today the UNIQUAC model is frequently applied in the description of phase equilibria (i.e. liquid–solid, liquid–liquid or liquid–vapor equilibrium). The UNIQUAC model also serves as the basis of the development of the group contribution method UNIFAC,^[3] where molecules are subdivided into functional groups. In fact, UNIQUAC is equal to UNIFAC for mixtures of molecules, which are not subdivided; e.g. the binary systems water-methanol, methanol-acryonitrile and formaldehyde-DMF.

A more thermodynamically consistent form of UNIQUAC is given by the more recent COSMOSPACE and the equivalent GEQUAC model.^[4]

Equations

Like most local composition models, UNIQUAC splits excess Gibbs free energy into a combinatorial and a residual contribution:

G^{E}=(G^{E})^{C}+(G^{E})^{R}

The calculated activity coefficients of the i^th component then split likewise:

\ln \gamma _{i}=\ln \gamma _{i}^{C}+\ln \gamma _{i}^{R}

The first is an entropic term quantifying the deviation from ideal solubility as a result of differences in molecule shape. The latter is an enthalpic^{[nb 1]} correction caused by the change in interacting forces between different molecules upon mixing.

Combinatorial contribution

The combinatorial contribution accounts for shape differences between molecules and affects the entropy of the mixture and is based on the lattice theory. The Stavermann–Guggenheim equation is used to approximate this term from pure chemical parameters, using the relative Van der Waals volumes r_i and surface areas q_i^{[nb 2]} of the pure chemicals:

{\frac {G^{E}}{RT}}=\sum _{i}\,x_{i}\ln {V_{i}}+{\frac {z}{2}}q_{i}\,x_{i}\ln {\frac {F_{i}}{V_{i}}}

Differentiating yields the excess entropy γ^C,

\ln \gamma _{i}^{C}=(1-V_{i}+\ln V_{i})-{\frac {z}{2}}q_{i}\left(1-{\frac {V_{i}}{F_{i}}}+\ln {\frac {V_{i}}{F_{i}}}\right)

with the volume fraction per mixture mole fraction, V_i, for the i^th component given by:

V_{i}={\frac {r_{i}}{\sum _{j}x_{j}r_{j}}}

The surface area fraction per mixture molar fraction, F_i, for the i^th component is given by:

F_{i}={\frac {q_{i}}{\sum _{j}x_{j}q_{j}}}

The first three terms on the right hand side of the combinatorial term form the Flory–Huggins contribution, while the remaining term, the Guggenhem–Staverman correction, reduce this because connecting segments cannot be placed in all direction in space. This spatial correction shifts the result of the Flory–Huggins term about 5% towards an ideal solution. The coordination number, z, i.e. the number of close interacting molecules around a central molecule, is frequently set to 10. It is based on the coordination number of an methylene group in a long chain, which has in the approximation of a hexagonal close packing structure of spheres 10 intermolecular and 2 bonds ^{[nb 3]} .

In the case of infinite dilution for a binary mixture, the equations for the combinatorial contribution reduce to:

{\begin{cases}\ln \gamma _{1}^{C,\infty }=1-{\dfrac {r_{1}}{r_{2}}}+\ln {\dfrac {r_{1}}{r_{2}}}-{\dfrac {z}{2}}q_{1}\left(1-{\dfrac {r_{1}q_{2}}{r_{2}q_{1}}}+\ln {\dfrac {r_{1}q_{2}}{r_{2}q_{1}}}\right)\\\ln \gamma _{2}^{C,\infty }=1-{\dfrac {r_{2}}{r_{1}}}+\ln {\dfrac {r_{2}}{r_{1}}}-{\dfrac {z}{2}}q_{2}\left(1-{\dfrac {r_{2}q_{1}}{r_{1}q_{2}}}+\ln {\dfrac {r_{2}q_{1}}{r_{1}q_{2}}}\right)\end{cases}}

This pair of equations show that molecules of same shape, i.e. same r and q parameters, have $\gamma _{1}^{C,\infty }=\gamma _{2}^{C,\infty }=1$ .

Residual contribution

The residual, enthalpic term contains an empirical parameter, $\tau _{ij}$ , which is determined from the binary interaction energy parameters. The expression for the residual activity coefficient for molecule i is:

\ln \gamma _{i}^{R}=q_{i}\left(1-\ln {\frac {\sum _{j}q_{j}x_{j}\tau _{ji}}{\sum _{j}q_{j}x_{j}}}-\sum _{j}{\frac {q_{j}x_{j}\tau _{ij}}{\sum _{k}q_{k}x_{k}\tau _{kj}}}\right)

with

\tau _{ij}=e^{-\Delta u_{ij}/{RT}}

$\Delta u_{ii}$ [J/mol] is the binary interaction energy parameter. Theory defines $\Delta u_{ij}=u_{ij}-u_{ii}$ , and $\Delta u_{ji}=u_{ji}-u_{jj}$ , where $u_{ij}$ is the interaction energy between molecules $i$ and $j$ . The interaction energy parameters are usually determined from activity coefficients, vapor-liquid, liquid-liquid, or liquid-solid equilibrium data.

Usually $\Delta u_{ij}\neq \Delta u_{ji}$ , because the energies of evaporation (i.e. $\Delta u_{ii}$ ), are in many cases different, while the energy of interaction between molecule i and j is symmetric, and therefore $u_{ij}=u_{ji}$ . If the interactions between the j molecules and i molecules is the same as between molecules i and j, there is no excess energy of mixing, $\Delta u_{ij}=\Delta u_{ji}=0$ . And thus $\gamma _{i}^{R}=1$ .

Alternatively, in some process simulation software $\tau _{ij}$ can be expressed as follows :

\ln \tau _{ij}=A_{ij}+B_{ij}/T+C_{ij}\ln(T)+D_{ij}T+E_{ij}/T^{2}

.

The C, D, and E coefficients are primarily used in fitting liquid–liquid equilibria data (with D and E rarely used at that). The C coefficient is useful for vapor-liquid equilibria data as well. The use of such an expression ignores the fact that on a molecular level the energy, $\Delta u_{ij}$ , is temperature independent. It is a correction to repair the simplifications, which were applied in the derivation of the model.

Applications (phase equilibrium calculations)

Activity coefficients can be used to predict simple phase equilibria (vapour–liquid, liquid–liquid, solid–liquid), or to estimate other physical properties (e.g. viscosity of mixtures). Models such as UNIQUAC allow chemical engineers to predict the phase behavior of multicomponent chemical mixtures. They are commonly used in process simulation programs to calculate the mass balance in and around separation units.

Parameters determination

UNIQUAC requires two basic underlying parameters: relative surface and volume fractions are chemical constants, which must be known for all chemicals (q_i and r_i parameters, respectively). Empirical parameters between components that describes the intermolecular behaviour. These parameters must be known for all binary pairs in the mixture. In a quaternary mixture there are six such parameters (1–2,1–3,1–4,2–3,2–4,3–4) and the number rapidly increases with additional chemical components. The empirical parameters are obtained by a correlation process from experimental equilibrium compositions or activity coefficients, or from phase diagrams, from which the activity coefficients themselves can be calculated. An alternative is to obtain activity coefficients with a method such as UNIFAC, and the UNIFAC parameters can then be simplified by fitting to obtain the UNIQUAC parameters. This method allows for the more rapid calculation of activity coefficients, rather than direct usage of the more complex method.

Remark that the determination of parameters from LLE data can be difficult depending on the complexity of the studied system. For this reason it is necessary to confirm the consistency of the obtained parameters in the whole range of compositions (including binary subsystems, experimental and calculated lie-lines, Hessian matrix, etc.).^[5]^[6]

Newer developments

UNIQUAC has been extended by several research groups. Some selected derivatives are: UNIFAC, a method which permits the volume, surface and in particular, the binary interaction parameters to be estimated. This eliminates the use of experimental data to calculate the UNIQUAC parameters,^[3] extensions for the estimation of activity coefficients for electrolytic mixtures,^[7] extensions for better describing the temperature dependence of activity coefficients,^[8] and solutions for specific molecular arrangements.^[9]

The DISQUAC model advances UNIFAC by replacing UNIFAC's semi-empirical group-contribution model with an extension of the consistent theory of Guggenheim's UNIQUAC. By adding a "dispersive" or "random-mixing physical" term, it better predicts mixtures of molecules with both polar and non-polar groups. However, separate calculation of the dispersive and quasi-chemical terms means the contact surfaces are not uniquely defined. The GEQUAC model advances DISQUAC slightly, by breaking polar groups into individual poles and merging the dispersive and quasi-chemical terms.

Notes

↑ Here it is assumed that the enthalpy change upon mixing can be assumed to be equal to the energy upon mixing, since the liquid excess molar volume is small and Δ H^ex=ΔU^ex+V^ex ΔP ≈ ΔU
↑ It is assumed that all molecules have the same coordination number as the methylene group of an alkane, which is the reference to calculate the relative volume and surface area.
↑ By setting q_I and r_i to the value an infinite long chain, infinite times the value of the methylene group, one finds with Eqn. B3 of the original paper the limiting value z=10.

Related Research Articles

In a chemical reaction, chemical equilibrium is the state in which both the reactants and products are present in concentrations which have no further tendency to change with time, so that there is no observable change in the properties of the system. This state results when the forward reaction proceeds at the same rate as the reverse reaction. The reaction rates of the forward and backward reactions are generally not zero, but they are equal. Thus, there are no net changes in the concentrations of the reactants and products. Such a state is known as dynamic equilibrium.

In physics and chemistry, an equation of state is a thermodynamic equation relating state variables, which describe the state of matter under a given set of physical conditions, such as pressure, volume, temperature, or internal energy. Most modern equations of state are formulated in the Helmholtz free energy. Equations of state are useful in describing the properties of pure substances and mixtures in liquids, gases, and solid states as well as the state of matter in the interior of stars.

In probability theory and statistics, the gamma distribution is a two-parameter family of continuous probability distributions. The exponential distribution, Erlang distribution, and chi-squared distribution are special cases of the gamma distribution. There are two equivalent parameterizations in common use:

With a shape parameter $k$ and a scale parameter $θ$
With a shape parameter $and an inverse scale parameter, called a rate parameter.$

The equilibrium constant of a chemical reaction is the value of its reaction quotient at chemical equilibrium, a state approached by a dynamic chemical system after sufficient time has elapsed at which its composition has no measurable tendency towards further change. For a given set of reaction conditions, the equilibrium constant is independent of the initial analytical concentrations of the reactant and product species in the mixture. Thus, given the initial composition of a system, known equilibrium constant values can be used to determine the composition of the system at equilibrium. However, reaction parameters like temperature, solvent, and ionic strength may all influence the value of the equilibrium constant.

In thermodynamics, an activity coefficient is a factor used to account for deviation of a mixture of chemical substances from ideal behaviour. In an ideal mixture, the microscopic interactions between each pair of chemical species are the same and, as a result, properties of the mixtures can be expressed directly in terms of simple concentrations or partial pressures of the substances present e.g. Raoult's law. Deviations from ideality are accommodated by modifying the concentration by an activity coefficient. Analogously, expressions involving gases can be adjusted for non-ideality by scaling partial pressures by a fugacity coefficient.

The Gibbs adsorption isotherm for multicomponent systems is an equation used to relate the changes in concentration of a component in contact with a surface with changes in the surface tension, which results in a corresponding change in surface energy. For a binary system, the Gibbs adsorption equation in terms of surface excess is:

The Kelvin equation describes the change in vapour pressure due to a curved liquid–vapor interface, such as the surface of a droplet. The vapor pressure at a convex curved surface is higher than that at a flat surface. The Kelvin equation is dependent upon thermodynamic principles and does not allude to special properties of materials. It is also used for determination of pore size distribution of a porous medium using adsorption porosimetry. The equation is named in honor of William Thomson, also known as Lord Kelvin.

<span class="mw-page-title-main">UNIFAC</span> Liquid equilibrium model in statistical thermodynamics

In statistical thermodynamics, the UNIFAC method is a semi-empirical system for the prediction of non-electrolyte activity in non-ideal mixtures. UNIFAC uses the functional groups present on the molecules that make up the liquid mixture to calculate activity coefficients. By using interactions for each of the functional groups present on the molecules, as well as some binary interaction coefficients, the activity of each of the solutions can be calculated. This information can be used to obtain information on liquid equilibria, which is useful in many thermodynamic calculations, such as chemical reactor design, and distillation calculations.

The non-random two-liquid model is an activity coefficient model that correlates the activity coefficients $of a compound with its mole fractions in the liquid phase concerned. It is frequently applied in the field of chemical engineering to calculate phase equilibria. The concept of NRTL is based on the hypothesis of Wilson that the local concentration around a molecule is different from the bulk concentration. This difference is due to a difference between the interaction energy of the central molecule with the molecules of its own kind and that with the molecules of the other kind . The energy difference also introduces a non-randomness at the local molecular level. The NRTL model belongs to the so-called local-composition models. Other models of this type are the Wilson model, the UNIQUAC model, and the group contribution model UNIFAC. These local-composition models are not thermodynamically consistent for a one-fluid model for a real mixture due to the assumption that the local composition around molecule i is independent of the local composition around molecule j. This assumption is not true, as was shown by Flemr in 1976. However, they are consistent if a hypothetical two-liquid model is used.$

<span class="mw-page-title-main">Diffusion</span> Transport of dissolved species from the highest to the lowest concentration region

Diffusion is the net movement of anything generally from a region of higher concentration to a region of lower concentration. Diffusion is driven by a gradient in Gibbs free energy or chemical potential. It is possible to diffuse "uphill" from a region of lower concentration to a region of higher concentration, as in spinodal decomposition. Diffusion is a stochastic process due to the inherent randomness of the diffusing entity and can be used to model many real-life stochastic scenarios. Therefore, diffusion and the corresponding mathematical models are used in several fields beyond physics, such as statistics, probability theory, information theory, neural networks, finance, and marketing.

PSRK is an estimation method for the calculation of phase equilibria of mixtures of chemical components. The original goal for the development of this method was to enable the estimation of properties of mixtures containing supercritical components. This class of substances cannot be predicted with established models, for example UNIFAC.

Pitzer equations are important for the understanding of the behaviour of ions dissolved in natural waters such as rivers, lakes and sea-water. They were first described by physical chemist Kenneth Pitzer. The parameters of the Pitzer equations are linear combinations of parameters, of a virial expansion of the excess Gibbs free energy, which characterise interactions amongst ions and solvent. The derivation is thermodynamically rigorous at a given level of expansion. The parameters may be derived from various experimental data such as the osmotic coefficient, mixed ion activity coefficients, and salt solubility. They can be used to calculate mixed ion activity coefficients and water activities in solutions of high ionic strength for which the Debye–Hückel theory is no longer adequate. They are more rigorous than the equations of specific ion interaction theory, but Pitzer parameters are more difficult to determine experimentally than SIT parameters.

The Margules activity model is a simple thermodynamic model for the excess Gibbs free energy of a liquid mixture introduced in 1895 by Max Margules. After Lewis had introduced the concept of the activity coefficient, the model could be used to derive an expression for the activity coefficients $of a compound i in a liquid, a measure for the deviation from ideal solubility, also known as Raoult's law.$

Equilibrium chemistry is concerned with systems in chemical equilibrium. The unifying principle is that the free energy of a system at equilibrium is the minimum possible, so that the slope of the free energy with respect to the reaction coordinate is zero. This principle, applied to mixtures at equilibrium provides a definition of an equilibrium constant. Applications include acid–base, host–guest, metal–complex, solubility, partition, chromatography and redox equilibria.

MOSCED is a thermodynamic model for the estimation of limiting activity coefficients. From a historical point of view MOSCED can be regarded as an improved modification of the Hansen method and the Hildebrand solubility model by adding higher interaction term such as polarity, induction and separation of hydrogen bonding terms. This allows the prediction of polar and associative compounds, which most solubility parameter models have been found to do poorly. In addition to making quantitative prediction, MOSCED can be used to understand fundamental molecular level interaction for intuitive solvent selection and formulation.

The Van Laar equation is a thermodynamic activity model, which was developed by Johannes van Laar in 1910-1913, to describe phase equilibria of liquid mixtures. The equation was derived from the Van der Waals equation. The original van der Waals parameters didn't give good description of vapor-liquid equilibria of phases, which forced the user to fit the parameters to experimental results. Because of this, the model lost the connection to molecular properties, and therefore it has to be regarded as an empirical model to correlate experimental results.

In chemical thermodynamics, excess properties are properties of mixtures which quantify the non-ideal behavior of real mixtures. They are defined as the difference between the value of the property in a real mixture and the value that would exist in an ideal solution under the same conditions. The most frequently used excess properties are the excess volume, excess enthalpy, and excess chemical potential. The excess volume, internal energy, and enthalpy are identical to the corresponding mixing properties; that is,

VTPR is an estimation method for the calculation of phase equilibria of mixtures of chemical components. The original goal for the development of this method was to enable the estimation of properties of mixtures which contain supercritical components. These class of substances couldn't be predicted with established models like UNIFAC.

COSMO-RS is a quantum chemistry based equilibrium thermodynamics method with the purpose of predicting chemical potentials µ in liquids. It processes the screening charge density σ on the surface of molecules to calculate the chemical potential µ of each species in solution. Perhaps in dilute solution a constant potential must be considered. As an initial step a quantum chemical COSMO calculation for all molecules is performed and the results are stored in a database. In a separate step COSMO-RS uses the stored COSMO results to calculate the chemical potential of the molecules in a liquid solvent or mixture. The resulting chemical potentials are the basis for other thermodynamic equilibrium properties such as activity coefficients, solubility, partition coefficients, vapor pressure and free energy of solvation. The method was developed to provide a general prediction method with no need for system specific adjustment.

The shear viscosity of a fluid is a material property that describes the friction between internal neighboring fluid surfaces flowing with different fluid velocities. This friction is the effect of (linear) momentum exchange caused by molecules with sufficient energy to move between these fluid sheets due to fluctuations in their motion. The viscosity is not a material constant, but a material property that depends on temperature, pressure, fluid mixture composition, local velocity variations. This functional relationship is described by a mathematical viscosity model called a constitutive equation which is usually far more complex than the defining equation of shear viscosity. One such complicating feature is the relation between the viscosity model for a pure fluid and the model for a fluid mixture which is called mixing rules. When scientists and engineers use new arguments or theories to develop a new viscosity model, instead of improving the reigning model, it may lead to the first model in a new class of models. This article will display one or two representative models for different classes of viscosity models, and these classes are:

References

↑ Abrams, Denis S.; Prausnitz, John M. (1975). "Statistical thermodynamics of liquid mixtures: A new expression for the excess Gibbs energy of partly or completely miscible systems". AIChE Journal. 21 (1): 116–128. doi:10.1002/aic.690210115. ISSN 0001-1541.
1 2 Maurer, G.; Prausnitz, J.M. (1978). "On the derivation and extension of the uniquac equation". Fluid Phase Equilibria. 2 (2): 91–99. doi:10.1016/0378-3812(78)85002-X. ISSN 0378-3812.
1 2 Fredenslund, Aage; Jones, Russell L.; Prausnitz, John M. (1975). "Group-contribution estimation of activity coefficients in nonideal liquid mixtures". AIChE Journal. 21 (6): 1086–1099. doi:10.1002/aic.690210607. ISSN 0001-1541.
↑ Egner, K.; Gaube, J.; Pfennig, A. (1997). "GEQUAC, an excess Gibbs energy model for simultaneous description of associating and non-associating liquid mixtures". Berichte der Bunsengesellschaft für Physikalische Chemie. 101 (2): 209–218. doi:10.1002/bbpc.19971010208. ISSN 0005-9021.
↑ Marcilla, Antonio; Reyes-Labarta, Juan A.; Olaya, M.Mar (2017). "Should we trust all the published LLE correlation parameters in phase equilibria? Necessity of their Assessment Prior to Publication". Fluid Phase Equilibria. 433: 243–252. doi:10.1016/j.fluid.2016.11.009. hdl: 10045/66521 .
↑ Graphical User Interface, (GUI). "Topological Analysis of the Gibbs Energy Function (Liquid-Liquid Equilibrium Correlation Data. Including a Thermodinamic Review and Tie-lines/Hessian matrix analysis)". University of Alicante (Reyes-Labarta et al. 2015-17). hdl:10045/51725.{{cite journal}}: |last1= has generic name (help); Cite journal requires |journal= (help)
↑ "The Extended UNIQUAC model".
↑ Wisniewska-Goclowska B., Malanowski S.K., “A new modification of the UNIQUAC equation including temperature dependent parameters”, Fluid Phase Equilib., 180, 103–113, 2001
↑ Andreas Klamt, Gerard J. P. Krooshof, Ross Taylor “COSMOSPACE: Alternative to conventional activity-coefficient models”, AIChE J., 48(10), 2332–2349,2004

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[5] Here it is assumed that the enthalpy change upon mixing can be assumed to be equal to the energy upon mixing, since the liquid excess molar volume is small and Δ H^ex=ΔU^ex+V^ex ΔP ≈ ΔU

[6] It is assumed that all molecules have the same coordination number as the methylene group of an alkane, which is the reference to calculate the relative volume and surface area.

[7] By setting q_I and r_i to the value an infinite long chain, infinite times the value of the methylene group, one finds with Eqn. B3 of the original paper the limiting value z=10.

[AbramsPrausnitz1975-1] Abrams, Denis S.; Prausnitz, John M. (1975). "Statistical thermodynamics of liquid mixtures: A new expression for the excess Gibbs energy of partly or completely miscible systems". AIChE Journal. 21 (1): 116–128. doi:10.1002/aic.690210115. ISSN 0001-1541.

[:1-2] 1 2 Maurer, G.; Prausnitz, J.M. (1978). "On the derivation and extension of the uniquac equation". Fluid Phase Equilibria. 2 (2): 91–99. doi:10.1016/0378-3812(78)85002-X. ISSN 0378-3812.

[:0-3] 1 2 Fredenslund, Aage; Jones, Russell L.; Prausnitz, John M. (1975). "Group-contribution estimation of activity coefficients in nonideal liquid mixtures". AIChE Journal. 21 (6): 1086–1099. doi:10.1002/aic.690210607. ISSN 0001-1541.

[EgnerGaube1997-4] Egner, K.; Gaube, J.; Pfennig, A. (1997). "GEQUAC, an excess Gibbs energy model for simultaneous description of associating and non-associating liquid mixtures". Berichte der Bunsengesellschaft für Physikalische Chemie. 101 (2): 209–218. doi:10.1002/bbpc.19971010208. ISSN 0005-9021.

[8] Marcilla, Antonio; Reyes-Labarta, Juan A.; Olaya, M.Mar (2017). "Should we trust all the published LLE correlation parameters in phase equilibria? Necessity of their Assessment Prior to Publication". Fluid Phase Equilibria. 433: 243–252. doi:10.1016/j.fluid.2016.11.009. hdl: 10045/66521 .

[9] Graphical User Interface, (GUI). "Topological Analysis of the Gibbs Energy Function (Liquid-Liquid Equilibrium Correlation Data. Including a Thermodinamic Review and Tie-lines/Hessian matrix analysis)". University of Alicante (Reyes-Labarta et al. 2015-17). hdl:10045/51725.{{cite journal}}: |last1= has generic name (help); Cite journal requires |journal= (help)

[10] "The Extended UNIQUAC model".

[NewExtention-11] Wisniewska-Goclowska B., Malanowski S.K., “A new modification of the UNIQUAC equation including temperature dependent parameters”, Fluid Phase Equilib., 180, 103–113, 2001

[12] Andreas Klamt, Gerard J. P. Krooshof, Ross Taylor “COSMOSPACE: Alternative to conventional activity-coefficient models”, AIChE J., 48(10), 2332–2349,2004

[1]

[2]

[3]

[4]

[nb 1]

[nb 2]

[nb 3]

[5]

[6]

[7]

[8]

[9]