Determination of the Ionic Association Constants of Na+ with CO32- and HCO3- Ions, in NaCl-NaHCO3-H2O Ternary Systems, at 25 °C.

The ionic association constants of sodium with carbonate ion (K1C') and acidic carbonate ions (K2C') were measured in NaCl-NaHCO3-H2O ternary systems to determine the distribution of sodium among the chemical species present in the growth medium of Chlorella homosphaera 424 algae. The mean activity coefficients of sodium chloride (in pure sodium chloride and in a mixture of electrolytes) were determined experimentally using two electrochemical cells, namely Ag, AgCl| KCl (3 M)|| NH4NO3 (1 M)| NaCl (mNaCl)| Na+-ISE and Ag, AgCl|KCl (3 M)|| NH4NO3 (1 M)| NaCl (mNaCl)| Cl--ISE. The studies carried out show that the values of the association constants of K1C' and K2C' do not depend on the composition of the medium, but only on the effective ionic strength. The experimentally obtained γNaCl0 values in the binary system are comparable to the mean activity coefficients values for NaCl, calculated using data from the literature, with -0.9 to 0.1% relative standard deviation. The obtained results show that the experimentally determined mean activity coefficient in the ternary system, γNaCl, is smaller than γNaCl0 in the binary system over the entire field of ionic strengths studied. The ternary system NaCl-NaHCO3-H2O obeys Harned's rule.

Modeling the Solubility of Phenolic Acids in Aqueous Media at 37 °C.

In this work, the solubility of vanillic, gallic, syringic, p-coumaric, ferulic and caffeic acids was determined at 37 °C under different conditions, namely pure water and two different ionic media, NaCl(aq) and NaClO4(aq), at different ionic strengths (i.e., 0.16, 0.50, 1.0, 2.0 and 3.0 M). The solubility in water of all the acids was found to be higher than that in both of the ionic media. Moreover, the solubility of hydroxycinnamic acids was lower than that of hydroxybenzoic acids. The activity coefficients of neutral species were calculated from these data; this knowledge is necessary when modeling the dependence of equilibrium constants on the ionic strength. Results obtained in this work can be useful for further studies regarding complex formation equilibria between these ligands and bioavailable metal cations.

Cluster-Based Thermodynamics of Interacting Dice in a Lattice.

In this paper, a model for two-component systems of six-sided dice in a simple cubic lattice is developed, based on a basic cluster approach previously proposed. The model represents a simplified picture of liquid mixtures of molecules with different interaction sites on their surfaces, where each interaction site can be assigned an individual energetic property to account for cooperative effects. Based on probabilities that characterize the sequential construction of the lattice using clusters, explicit expressions for the Shannon entropy, synonymously used as thermodynamic entropy, and the internal energy of the system are derived. The latter are used to formulate the Helmholtz free energy that is minimized to determine thermodynamic bulk properties of the system in equilibrium. The model is exemplarily applied to mixtures that contain distinct isomeric configurations of molecules, and the results are compared with the Monte-Carlo simulation results as a benchmark. The comparison shows that the model can be applied to distinguish between isomeric configurations, which suggests that it can be further developed towards an excess Gibbs-energy, respectively, activity coefficient model for chemical engineering applications.

Activity Coefficients for Liquid Organic Reactions: Towards a Better Understanding of True Kinetics with the Synthesis of Jasmin Aldehyde as Showcase.

The aldol condensation of benzaldehyde and heptanal is taken as an example of reversible liquid phase organic reactions to show that inclusion of activity coefficients reveal distinct differences in conversion and product distribution when different solvents methanol, ethanol, n-propanol, or n-butanol are used. The purpose of this work is to show a pronounced solvent effect for a given set of identical kinetic parameters, i.e., the same liquid phase kinetics can result in different conversion and yield values, depending on the choice of solvent. It was shown that subsequent parameter estimation without inclusion of the activity coefficients resulted in a pronounced deviation from the 'true' kinetics, up to a factor of 30. It is proposed that the usage of average activity coefficients gives already a significant improvement, resulting in acceptable parameter estimates.

Experimental Evidence for Azeotrope Formation from Protic Ionic Liquids.

Herein, we present experimental evidence that protic ionic liquids (PILs), derived from 1 : 1 liquid mixtures of the organic superbases 1,5-diazabicyclo[4.3.0]non-5-ene (DBN) and 1,8-diazabicyclo[5.4.0]undec-7-ene (DBU) with carboxylic acids, form azeotropic mixtures with acid/base molar fractions different from 1 : 1. The ability of the carboxylic acids to form strong hydrogen bonds with the PIL ion pair leads to an azeotropic composition richer in the acid component. The results show that the azeotropic composition is ruled by the extent of acid-base equilibrium and the relative volatility of the neutral species in the PIL medium. The PILs show marked negative deviations from Raoult's Law with the stronger superbase (DBU) leading to an azeotropic composition closer to the equimolar 1 : 1 ratio.

Modeling food matrix effects on chemical reactivity: Challenges and perspectives.

The same chemical reaction may be different in terms of its position of the equilibrium (i.e., thermodynamics) and its kinetics when studied in different foods. The diversity in the chemical composition of food and in its structural organization at macro-, meso-, and microscopic levels, that is, the food matrix, is responsible for this difference. In this viewpoint paper, the multiple, and interconnected ways the food matrix can affect chemical reactivity are summarized. Moreover, mechanistic and empirical approaches to explain and predict the effect of food matrix on chemical reactivity are described. Mechanistic models aim to quantify the effect of food matrix based on a detailed understanding of the chemical and physical phenomena occurring in food. Their applicability is limited at the moment to very simple food systems. Empirical modeling based on machine learning combined with data-mining techniques may represent an alternative, useful option to predict the effect of the food matrix on chemical reactivity and to identify chemical and physical properties to be further tested. In such a way the mechanistic understanding of the effect of the food matrix on chemical reactions can be improved.

Predicting the Solvent Effect on Esterification Kinetics.

It is well known that solvents influence reaction kinetics. The classical concentration-based kinetic modeling is unable to describe these effects. In this work, the reaction kinetics was studied for the esterifications of acetic acid and propionic acid with ethanol at 303.15 K. It was found that the reactant ratio as well as the applied solvents (acetonitrile, tetrahydrofurane, dimethylformamide) significantly affect the reaction rate. The thermodynamic model PC-SAFT was applied to account for the interactions between the reacting species and the solvents via activity coefficients. This allowed the identification of solvent-independent kinetic constants and the prediction of the solvent effect on reaction kinetics in almost quantitative agreement with experimental data. The presented approach shows the importance of taking into account thermodynamic non-idealities and significantly reduces experimental effort for finding the best solvent candidate for a given target reaction.

Solid-liquid equilibria of triacylglycerols and vitamin E mixtures.

Oils and fats are important ingredients for food and pharmaceutical industries. Their main compounds, such as triacylglycerols (TAG), are responsible for determining their physical properties during food storage and consumption. Lipid-rich foods are also sources of minority compounds, which is the case of vitamin E, mainly represented by (±)-α-tocopherol. These compounds can interact with the main lipid molecules in food formulation leading to modification on lipids' physicochemical properties during processes, storage, as well as during digestion, possibly altering their nutritional functionalities, which is the case of vitamin E antioxidant abilities, but also their solubility in the systems. In this case, the study of the phase-behavior between (±)-α-tocopherol and lipid compounds can elucidate these physicochemical changings. Therefore, this work was aimed at determining the solid-liquid equilibrium (SLE) of binary mixtures of TAG (tripalmitin, triolein and tristearin) and (±)-α-tocopherol including the complete description of their phase diagrams. Melting data were evaluated by Differential Scanning Calorimetry, Microscopy, X-Ray Diffraction, and thermodynamic modeling by using Margules, UNIFAC, and COSMO-SAC models. Experimental results showed that systems presented a monotectic-like behavior, with a significant decreasing in TAG melting temperature by the addition of (±)-α-tocopherol. This high affinity and attractive strengths between these molecules were also observed by thermodynamic modeling, whose absolute deviations were below 2 %. Micrographs and X-ray diffraction evidenced the possible formation of solid solutions. Both behaviors are interesting by avoiding phase separation on food in solid and liquid phases, possibly improving the antioxidant role the (±)-α-tocopherol in lipid-base systems.

Solid-liquid equilibrium of free form of oil contaminants (3-MCPD and glycidol) in lipidic systems.

Esters of 3-monochloropropane-1,2-diol (3-MCPD) and glycidol are reported vegetable fats and oils contaminants formed during processing. During digestion, esterified forms are hydrolyzed making the free forms, with high toxicity to human health, possibly available in the digestive system, which could depend, among other factors, on the phase condition stablished in the gastrointestinal tract between contaminants and lipids. Therefore, this work was aimed at evaluating the solid-liquid equilibrium (SLE) of binary and pseudobinary mixtures of fatty acids (palmitic, oleic and stearic), triacylglycerols (tripalmitin, triolein and tristearin) and a partial acylglycerols, with 3-MCPD and glycidol, by differential scanning calorimetry (DSC) and thermodynamic modeling with Margules 2 and 3 suffixes, UNIFAC and ideal models. Melting properties of the contaminants were determined by DSC thermograms, with microstructure micrographs and compared to some predictive group contribution models (GC). 14 complete SLE phase diagrams could be determined. Results showed that the lipid compounds had high affinity with the free contaminants with probably the formation of liquid crystalline structures and/or solid solutions. The Margules 3 suffixes model showed the best fit with the lowest average relative deviations, no higher than 1%, although the parameters probably incorporated the non-ideality promoted by the crystalline phase behavior. The formation of mesophases and/or solid solutions with very high melting temperature probably means a difficulty of separating the contaminants from the lipid matrix in which, on the other hand, could avoid separation process or affect (probably decreasing) their absorption during digestion.

A New Model Algorithm for Estimating the Inhalation Exposure Resulting from the Spraying of (Semi)-Volatile Binary Liquid Mixtures (SprayEva).

The spraying of liquid multicomponent mixtures is common in many professional and industrial settings. Typical examples are cleaning agents, additives, coatings, and biocidal products. In all of these examples, hazardous substances can be released in the form of aerosols or vapours. For occupational and consumer risk assessment in regulatory contexts, it is therefore important to know the exposure which results from the amount of chemicals in the surrounding air. In this research, a mechanistic mass balance model has been developed that covers the spraying of (semi)-volatile substances, taking into account combined exposure to spray mist, evaporation from droplets, and evaporation from surfaces as well as the nonideal behaviour of components in liquids and backpressure effects. For wall-spraying scenarios, an impaction module has been developed that quantifies the amount of overspray and the amount of material that lands on the wall. Mechanistically, the model is based on the assumption that continuous spraying can be approximated by a number of sequentially released spray pulses, each characterized by a certain droplet size, where the total aerosol exposure is obtained by summation over all release pulses. The corresponding system of differential equations is solved numerically using an extended Euler algorithm that is based on a discretisation of time and space. Since workers typically apply the product continuously, the treated area and the corresponding evaporating surface area grows over time. Time-dependent concentration gradients within the sprayed liquid films that may result from different volatilities of the components are therefore addressed by the proposed model. A worked example is presented to illustrate the calculated exposure for a scenario where aqueous solutions of H2O2 are sprayed onto surfaces as a biocidal product. The results reveal that exposure to H2O2 aerosol reaches relevant concentrations only during the spraying phase. Evaporation from sprayed surfaces takes place over much longer time periods, where backpressure effects caused by large emission sources can influence the shape of the concentration time curves significantly. The influence of the activity coefficients is not so pronounced. To test the plausibility of the developed model algorithm, a comparison of model estimates of SprayExpo, SprayEva, and ConsExpo with measured data is performed. Although the comparison is based on a limited number (N = 19) of measurement data, the results are nevertheless regarded as supportive and acceptable for the plausibility and predictive power of SprayEva.

Estimating Inhalation Exposure Resulting from Evaporation of Volatile Multicomponent Mixtures Using Different Modelling Approaches.

In many professional and industrial settings, liquid multicomponent mixtures are used as solvents, additives, coatings, biocidal products, etc. Since, in all of these examples, hazardous liquids can evaporate in the form of vapours, for risk assessments it is important to know the amount of chemicals in the surrounding air. Although several models are available in legal contexts, the current implementations seem to be unable to correctly simulate concentration changes that actually occur in volatile mixtures and in particular in thin films. In this research, the estimation of evaporation rates is based on models that take into account non-ideal behaviour of components in liquids and backpressure effects as well. The corresponding system of differential equations is solved numerically using an extended Euler algorithm that is based on a discretisation of time and space. Regarding air dispersion of volatile components, the model builds upon one-box and two-box mass balance models, because there is some evidence that these models, when selected and applied appropriately, can predict occupational exposures with sufficient precision. That way, numerical solutions for a wide variety of exposure scenarios with instantaneous and continuous/intermittent application, even considering "moving worker situations", can be obtained. A number of example calculations have been carried out on scenarios where binary aqueous solutions of hydrogen peroxide or glutaraldehyde are applied as a biocidal product to surfaces by wiping. The results reveal that backpressure effects caused by large emission sources as well as deviations from liquid-phase ideality can influence the shape of the concentration time curves significantly. The results also provide some evidence that near-/far-field models should be used to avoid underestimation of exposure in large rooms when small/medium areas are applied. However, the near-field/far-field model should not be used to estimate peak exposure assuming instantaneous application, because then the models tend to overestimate peak exposure significantly. Although the example calculations are restricted to aqueous binary mixtures, the proposed approach is general and can be used for arbitrary liquid multicomponent mixtures, as long as backpressure effects and liquid-phase non-idealities are addressed adequately.

Extremely Low Vapor-Pressure Data as Access to PC-SAFT Parameter Estimation for Ionic Liquids and Modeling of Precursor Solubility in Ionic Liquids.

Precursor solubility is a crucial factor in industrial applications, dominating the outcome of reactions and purification steps. The outcome and success of thermodynamic modelling of this industrially important property with equations of states, such as Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT), vastly depends on the quality of the pure-component parameters. The pure-component parameters for low-volatile compounds such as ionic liquids (ILs) have been commonly estimated using mixture properties, e. g. the osmotic pressure of aqueous solutions. This leads to parameters that depend on the solvent, and transferability to other mixtures often causes poor modeling results. Mixture-independent experimental properties would be a more suitable basis for the parameter estimation offering a way to universal parameter sets. Model parameters for ILs are available in the literature [10.1016/j.fluid.2012.05.029], but they were estimated using pure-IL density data. The present work focuses on a step towards a more universal estimation strategy that includes new experimental vapor-pressure data of the pure IL. ILs exhibit an almost negligible vapor pressure in magnitude of usually 10-5  Pa even at elevated temperatures. In this work, such vapor-pressure data of a series of 1-ethyl-3-methyl-imidazolium-based [C2 mim]-ILs with various IL-anions (e. g. tetrafluoroborate [BF4 ]- , hexafluorophosphate [PF6 ]- , bis(trifluoromethylsulfonyl)imide [NTf2 ]- ) were experimentally determined and subsequently used for PC-SAFT parameter estimation. The so-determined parameters were used to predict experimental molecular precursor solubility in ILs and infinitely diluted activity coefficients of various solvents in ILs. The parameters were further compared to modeling results using classical parametrization methods (use of liquid-density data only for the molecular PC-SAFT and the ion-based electrolyte PC-SAFT). As a result, the modeled precursor solubilities using the new approach are much more precise than using the classical parametrization methods, and required binary parameters were found to be much smaller (if needed). In sum, including the pure-component vapor-pressure data of ILs opens the door towards parameter estimation that is not biased by mixture data. This procedure might be suitable also for polymers and for all kind of ionic species but needs extension to ion-specific parametrization in the long term.

Solubility and thermodynamic analysis of famotidine polymorphs in pure solvents.

The present study investigates the solubility of famotidine polymorphs forms A and B between 298.15 K and 348.15 K in a range of pure solvents: water, methanol, ethanol, isopropanol and acetonitrile. Empirical and semi-empirical models have been fitted to solubility data determined experimentally by a gravimetric method. The solid phases have been characterized by FTIR and Raman spectroscopy, SEM and PXRD. In addition, heat capacities and melting data determined by DSC have been used to estimate the fusion thermodynamics and the activity of the solid phases as a function of temperature. The relationship between the famotidine polymorphs is monotropic, with form A being the stable polymorph. For both polymorphs, in terms of mass ratio, the solubility in the studied solvents decreases in the order methanol > water > ethanol > acetonitrile > isopropanol. The activity coefficient at saturation in all the solutions exceeds unity, showing a positive deviation with respect to ideality, which translates into solubilities significantly lower than the ideal values. Among the alcohols, a consistent correlation is observed between the polarity and the order of solubility.The Hildebrand solubility parameter is also well correlated with the order of solubilities in the studied solvents, with a higher solubility in more polar solvents, revealing the importance of the hydrogen bonding of the sulfamoyl group oxygens.

Analytical modeling of micelle growth. 3. Electrostatic free energy of ionic wormlike micelles - Effects of activity coefficients and spatially confined electric double layers.

HYPOTHESES: To correctly predict the aggregation number and size of wormlike micelles from ionic surfactants, the molecular-thermodynamic theory has to calculate the free energy per molecule in the micelle with accuracy better than 0.01 kT, which is a serious challenge. The problem could be solved if the effects of mutual confinement of micelle counterion atmospheres, as well as the effects of counterion binding, surface curvature and ionic interactions in the electric double layer (EDL), are accurately described. THEORY: The electric field is calculated using an appropriate cell model, which takes into account the aforementioned effects. Expressions for the activity coefficients have been used, which vary across the EDL and describe the electrostatic, hard sphere, and specific interactions between the ions. New approach for fast numerical calculation of the electrostatic free energy is developed. FINDINGS: The numerical results demonstrate the variation of quantities characterizing the EDL of cylindrical and spherical micelles with the rise of electrolyte concentration. The effect of activity coefficients leads to higher values of the free energy per surfactant molecule in the micelle as compared with the case of neglected ionic interactions. The results are essential for the correct prediction of the size of wormlike micelles from ionic surfactants. This study can be extended to mixed micelles of ionic and nonionic surfactants for interpretation of the observed synergistic effects.

Determination of physicochemical properties of ionic liquids by gas chromatography.

Over the years room temperature ionic liquids have gained attention as solvents with favorable environmental and technical features. Both chromatographic and conventional methods afford suitable tools for the study of their physicochemical properties. Use of gas chromatography compared to conventional methods for the measurement of physicochemical properties of ionic liquids have several advantages; very low sample concentrations, high accuracy, faster measurements, use of wider temperature range and the possibility to determine physicochemical properties of impure samples. Also, general purpose gas chromatography instruments are widely available in most laboratories thus alleviating the need to purchase more specific instruments for less common physiochemical measurements. Some of the main types of physicochemical properties of ionic liquids accessible using gas chromatography include gas-liquid partition constants, infinite dilution activity coefficients, partial molar quantities, solubility parameters, system constants of the solvation parameter model, thermal stability, transport properties, and catalytic and other surface properties.

Activity coefficients of binary methanol alcohol mixtures from cluster weighting.

The hydrogen bond network of different small alcohols is investigated via cluster analysis. Methanol/alcohol mixtures are studied with increasing chain length and branching of the molecule. Those changes can play an important role in different fields, including solvent and metal extraction. The extended tight binding method GFN2-xTB allows the evaluation and geometry optimization of thousands of clusters built via a genetic algorithm. Interaction energies and geometries are evaluated and discussed for the neat systems. Thermodynamic properties, such as vaporization enthalpies and activity coefficients, are calculated with the binary quantum cluster equilibrium (bQCE) approach using our in-house code peacemaker 2.8. Combined distribution functions of the distances against the angles of the hydrogen bonds are evaluated for neat and mixed clusters and weighted by the equilibrium populations achieved from bQCE calculations.

Solubility and thermodynamic analysis of ketoprofen in organic solvents.

The solubility of the racemic solid phase of ketoprofen (KTP) in methanol, ethanol, isopropanol, butanol, acetonitrile, ethyl acetate, 1,4-dioxane and toluene has been determined between 273 and 303 K by a gravimetric method. FTIR and Raman spectroscopy, SEM and PXRD, have been used to characterise the solid phase. The melting data and heat capacity of solid and melt have been determined by DSC, and used to estimate fusion thermodynamics and the activity of the solid phase as functions of temperature. Empirical and semi-empirical models have been fitted to experimental solubility data. The solution activity coefficients reveal positive deviation from ideality in all solvents except for in dioxane, and very close to ideality in methanol. The solubility is fairly high in the alcohols but decrease with increasing hydrocarbon chain. Generally and due to the presence of the carboxylic acid group, KTP is more readily dissolved in polar protic solvents, followed in order by polar aprotic and non-polar solvents. However, the highest solubility is found in dioxane, classified as a non-polar solvent, but notably though the molecule having two strong hydrogen bond accepting functionalities, and no hydrogen bond donation capability.

Thermodynamic Activity-Based Solvent Design for Bioreactions.

To improve the kinetics of enzyme-catalyzed reactions, cosolvents are commonly added to reaction mixtures. The search for a good cosolvent is still empirical and experimentally based. We discuss a thermodynamic activity-based approach that improves biocatalytic processes by predicting cosolvent influences on Michaelis constants, ultimately reducing time and cost.

Activity-based analysis of potentiometric pH titrations.

The discrepancy between concentrations and activities is a predicament well known to the analytical chemist. Because of the difficulty of determining activity coefficients, the standard technique for quantitative equilibrium studies is to work under a particular 'constant ionic strength' by adding an excess of an inert salt. Under such conditions, activity coefficients are approximately constant and can be taken into the equilibrium constants which are defined for the chosen ionic strength (I). Here we propose a fundamentally different approach. Throughout the numerical analysis of the titration data, activity coefficients for all individual species are approximated by well-known equations based on the work of Debye-Hückel. The computational analysis of the measurements strictly obeys the law of mass conservation and obeys the law of mass action only approximately. The main novelty is that now the addition of inert salts is no longer required and measurements are done at minimal I. Consequently, the thermodynamic equilibrium constants are now determined much more robustly based on experiments taken at low I. The approach has been tested and validated with the two very well investigated 3-protic phosphoric and citric acids. In summary: the technique of artificially keeping ionic strength constant has been replaced by improved computational analysis.

Study of the influence of surfactants on the activity coefficients and mass transfer coefficients of methanol in aqueous mixtures by reversed-flow gas chromatography.

This work focuses on the influences of surfactants on the activity coefficients, γ, of methanol in binary mixtures with water, as well as on the mass transfer coefficients, kc, for the evaporation of methanol, which is a ubiquitous component in the troposphere, from mixtures of methanol with water at various surfactant's and methanol's concentrations. The technique used is the Reversed-Flow Gas Chromatography (R.F.G.C.), a version of Inverse Gas Chromatography, which allows determining both parameters by performing only one experiment for the kc parameter and two experiments for the γ parameter. The kc and γ values decrease in the presence of the three surfactants used (CTAB, SDS, TRITON X-100) at all methanol's and surfactant's concentrations. The decrease in the methanol's molar fraction, at constant number of surfactant films leads to a decrease in the kc and γ values, while the decrease in the surfactant's concentration, at constant methanol's molar fraction leads to an increase in both the kc and γ parameters. Mass transfer coefficients for the evaporation of methanol at the surfactant films, are also calculated which are approximately between 4 and 5 orders of magnitude larger than the corresponding mass transfer coefficients at the liquid films. Finally, thicknesses of the boundary layer of methanol in the mixtures of methanol with water were determined. The quantities found are compared with those given in the literature or calculated theoretically using various empirical equations. The precision of the R.F.G.C. method for measuring γ and kc parameters is approximately high (94.3-98.0%), showing that R.F.G.C. can be used with success not only for the thermodynamic study of solutions, but also for the interphase transport.