Experimental Data of Fluid Phase Equilibria- Correlation and Prediction Models: A Review

Elements and Chemical Bonds Contribution Estimation of Activity Coefficients in Nonideal Liquid Mixtures

Based on the contribution of elements and chemical bonds, the UNICAC (Universal Quasi-Chemical elements and chemical bonds Activity Coefficient) method was proposed to estimate the activity coefficients of nonelectrolyte liquid mixtures. The UNICAC method defined 10 elements and 33 chemical bonds as contribution groups. The calculation of activity coefficients was divided into the combination term and the residual term. The combination term represents molecular size differences, and the residual term describes the interaction between molecules. The interaction energy parameters of 43 groups were regressed simultaneously with the experimental data of the vapor–liquid equilibrium of 1085 binary systems. According to the molecular structural information of compounds, the UNICAC method can accurately predict the activity coefficients of nonelectrolyte liquid mixtures. The vapor–liquid equilibrium of 16 groups of the binary system, which were not included in the parameters regress, was predicted using UNICAC. The average relative error of vapor composition was 1.53%. Compared with UNIFAC (2003), UNIFAC (Lyngby), UNIFAC (Dortmund), and ASOG (2011), the UNICAC model employs fewer parameters, provides a broader scope of application, and receives more precise predicted results of the vapor–liquid equilibrium. The UNICAC method would play an important reference role in the design of the chemical separation process.

Comparison of CP-PC-SAFT and SAFT-VR-Mie in Predicting Phase Equilibria of Binary Systems Comprising Gases and 1-Alkyl-3-methylimidazolium Ionic Liquids

This study compares performances of the Critical Point-based revision of Perturbed-Chain SAFT (CP-PC-SAFT) and the SAFT of Variable Range and Mie Potential (SAFT-VR-Mie) in predicting the available data on VLE, LLVE, critical loci and saturated phase densities of systems comprising CO, O2, CH4, H2S, SO2, propane, the refrigerants R22, R23, R114, R124, R125, R125, R134a, and R1234ze(E) and ionic liquids (ILs) with 1-alkyl-3-methylimidazolium ([Cnmim]+) cations and bis(trifluoromethanesulfonyl)imide ([NTf2]-), tetrafluoroborate ([BF4]-) and hexafluorophosphate ([PF6]-) anions. Both models were implemented in the entirely predictive manner with k12 = 0. The fundamental Global Phase Diagram considerations of the IL systems are discussed. It is demonstrated that despite a number of quantitative inaccuracies, both models are capable of reproducing the regularities characteristic for the considered systems, which makes them suitable for preliminary estimation of selectivity of the ILs in separating various gases.

Prediction of the Synergistic Glass Transition Temperature of Coamorphous Molecular Glasses Using Activity Coefficient Models

The glass transition temperature (T_g) of a binary miscible mixture of molecular glasses, termed a coamorphous glass, is often synergistically increased over that expected for an athermal mixture due to the strong interactions between the two components. This synergistic interaction is particularly important for the formulation of coamorphous pharmaceuticals since the molecular interactions and resulting T_g strongly impact stability against crystallization, dissolution kinetics, and bioavailability. Current models that describe the composition dependence of T_g for binary systems, including the Gordon-Taylor, Fox, Kwei, and Braun-Kovacs equations, fail to describe the behavior of coamorphous pharmaceuticals using parameters consistent with experimental ΔC_P and Δα. Here, we develop a robust thermodynamic approach extending the Couchman and Karasz method through the use of activity coefficient models, including the two-parameter Margules, non-random-two-liquid (NRTL), and three-suffix Redlich-Kister models. We find that the models, using experimental values of ΔC_P and fitting parameters related to the binary interactions, successfully describe observed synergistic elevations and inflections in the T_g versus composition response of coamorphous pharmaceuticals. Moreover, the predictions from the NRTL model are improved when the association-NRTL version of that model is used. Results are reported and discussed for four different coamorphous systems: indomethacin-glibenclamide, indomethacin-arginine, acetaminophen-indomethacin, and fenretinide-cholic acid.

Molecular Assembling in Mixtures of Hydrophilic 1-Butyl-1-Methylpyrrolidinium Dicyanamide Ionic Liquid and Water

The infrared absorbance spectrum of the ionic liquid 1-butyl-1-methylpyrrolidinium dicyanamide, mixed with water at two different concentrations, was measured between 160 and 300 K in the mid infrared range. Both mixtures do not crystallize on cooling; however, remarkably, the one with an ionic liquid (IL):water composition of 1:3 displays a cold crystallization process on heating in a restricted temperature range between 240 and 250 K. A portion of the water participates to the cold crystallization. On the contrary, with an IL:water composition of 1:6.6 no crystallization takes place. Upon water addition the vibration frequencies of the anion and of some lines of the cation are blue shifted, while the absorption lines of water are red shifted. These facts are interpreted as the evidence of the occurrence of the hydrogen bonding of water, as the hydrogen bonding acceptor with respect to the anion (anion∙∙∙O-H bonds develop) and as hydrogen donor for the cation (C-H∙∙∙O bonds can form). Microscopic inhomogeneities in the samples and their evolution with temperature are discussed.

Assessment of the SM12, SM8, and SMD Solvation Models for Predicting Limiting Activity Coefficients at 298.15 K

The SMx (x = 12, 8, or D) universal solvent models are implicit solvent models which using electronic structure calculations can compute solvation free energies at 298.15 K. While solvation free energy is an important thermophysical property, within the thermodynamic modeling of phase equilibrium, limiting (or infinite dilution) activity coefficients are preferred since they may be used to parameterize excess Gibbs free energy models to model phase equilibrium. Conveniently, the two quantities are related. Therefore the present study was performed to assess the ability to use the SMx universal solvent models to predict limiting activity coefficients. Two methods of calculating the limiting activity coefficient where compared: (1) the solvation free energy and self-solvation free energy were both predicted and (2) the self-solvation free energy was computed using readily available vapor pressure data. Overall the first method is preferred as it results in a cancellation of errors, specifically for the case in which water is a solute. The SM12 model was compared to both the Universal Quasichemical Functional-group Activity Coefficients (UNIFAC) and modified separation of cohesive energy density (MOSCED) models. MOSCED was the highest performer, yet had the smallest available compound inventory. UNIFAC and SM12 exhibited comparable performance. Therefore further exploration and research should be conducted into the viability of using the SMx models for phase equilibrium calculations.

Expanding the Solubility Parameter Method MOSCED to Pyridinium-, Quinolinium-, Pyrrolidinium-, Piperidinium-, Bicyclic-, Morpholinium-, Ammonium-, Phosphonium-, and Sulfonium-Based Ionic Liquids

MOSCED (modified separation of cohesive energy density) is a solubility parameter method that offers an improved treatment of association interactions. Solubility parameter methods are well known for their ability to both make quantitative predictions and offer a qualitative description of the underlying molecular-level driving forces, lending themselves to intuitive solvent selection and design. Currently, MOSCED parameters are available for 130 organic solvents, water, and 33 imidazolium-based room temperature ionic liquids (ILs). In this work, we expand MOSCED to cover 66 additional ILs containing the pyridinium, quinolinium, pyrrolidinium, piperidinium, bicyclic, morpholinium, ammonium, phosphonium, and sulfonium cations using 10,052 experimental limiting activity coefficients. The resulting parameters may readily be used to predict the phase behavior in mixtures involving ILs.