Application of a chemical equilibrium model to thermodynamic functions of transfer of alcohols from water to aqueous surfactants

In ternary aqueous solutions, hydrophobic solutes such as alcohols tend to aggregate with surfactants to form mixed micelles. These systems can be studied by meas of the functions of transfer of hydrophobic solutes from water to aqueous solutions of surfactant. These thermodynamic functions often go through extrema in the critical micellar concentration (CMC) region of the surfactant. A simple model based on interactions between surfactant and hydrophobic solute monomers, on the distribution of the hydrophobic solute between water and the micelles and on the shift in the CMC induced by the hydrophobic solute, can simulate the magnitude and trends of the transfer functions using parameters which are mostly derived from the binary systems. In order to check the model more quantitatively, volumes and heat capacities of transfer of alcohols from water to aqueous solutions of a nonionic surfactant, octyldimethylamine oxide, were measured. A quantitative agreement was achieved with three adjustable parameters. Good fits are also obtained for the transfers to the ionic surfactants, octylamine hydrobromide and sodium dodecylsulfate. When the equilibrium displacement contribution is small, the distribution constants and the partial molar properties of the alcohols in the micellar phase agree well with the parameters obtained with similar models.     

Volumes and heat capacities of some watersurfactant-alcohol ternary systems

The densities and volumetric heat capacities of urea and alcohols were measured in aqueous solutions of octylammonium bromide (OABr) and of OABr in aqueous urea and alcohol solutions. The alcohols studies were methanol, ethanol, 1-propanol, 2-propanol, n-butanol, t-butanol, n-pentanol, n-hexanol and 2-butoxyethanol (BE). In most experiments, the concentration of the reference solute was kept low, and volumes and heat capacities of transfer from water to the mixed solvent were calculated. A more complete study was made with the system BE-OABr-H₂O where both solutes were systematically changed. The observed trends in the thermodynamic functions can be explained through three effects: interactions between the reference solute and the cosolvent in the premiceller region of the surfactant or pre-aggregation region of the alcohol, a distribution of the reference solute between water and the micelle or microphase and an equilibrium displacement of the system, monomer-aggregate, in the vicinity of the reference solute.     

Micellar solubilization thermodynamics of acetophenone in surfactant mixtures: Case of benzyldimethyltetradecylammonium chloride and trimethyltetradecylammonium chloride

The thermodynamics of micellar solubilization of acetophenone in mixtures of two cationic surfactants [benzyldimethyltetradecylammonium chloride +trimethyltetradecylammonium chloride] has been derived from calorimetric measurements at controlled solute activity. The partition coefficient between micelles and water as well as the standard enthalpy and entropy of transfer between micelles and water were calculated. The results were compared to the case of benzylalcohol in the same cationic mixtures. For acetophenone, the variation of all thermodynamic transfer functions with micellar composition may be described by the regular solution formalism. The same conclusion has been achieved for most polar solutes in various surfactant mixtures: favorable interaction between unlike surfactants induces an unfavorable micellar solubilization. Exceptions should be found with the cases where solute solubilization induces profound micellar changes. It seems to be the case with some alcohols in the cationic surfactant mixtures studied.     

Activity of water in aqueous systems; a frequently neglected property

In this critical review, the significance of the term 'activity' is examined in the context of the properties of aqueous solutions. The dependence of the activity of water(l) at ambient pressure and 298.15 K on solute molality is examined for aqueous solutions containing neutral solutes, mixtures of neutral solutes and salts. Addition of a solute to water(l) always lowers its thermodynamic activity. For some solutes the stabilisation of water(l) is less than and for others more than in the case where the thermodynamic properties of the aqueous solution are ideal. In one approach this pattern is accounted for in terms of hydrate formation. Alternatively the pattern is analysed in terms of the dependence of practical osmotic coefficients on the composition of the aqueous solution and then in terms of solute-solute interactions. For salt solutions the dependence of the activity of water on salt molalities is compared with that predicted by the Debye-Hückel limiting law. The analysis is extended to consideration of the activities of water in binary aqueous mixtures. The dependence on mole fraction composition of the activity of water in binary aqueous mixtures is examined. Different experimental methods for determining the activity of water in aqueous solutions are critically reviewed. The role of water activity is noted in a biochemical context, with reference to the quality, stability and safety of food and finally with regard to health science.     

Hydration of aprotic donor solvents studied by means of FTIR spectroscopy

The paper attempts to explain the mutual influence of nonpolar and electron-donor groups on solute hydration, the problem of big importance for biological aqueous systems. Aprotic organic solvents have been used as model solutes, differing in electron-donating power. Hydration of acetonitrile, acetone, 2-butanone, and triethylamine has been studied by HDO and (partially) H2O spectra. The quantitative version of difference spectra method has been applied to determine solute-affected water spectra. Analysis of the data suggests that solvent-water interaction via the donor center of the solute is averaged between water-water interactions around the solute. Such behavior can be simply explained by the model of solute rotating in a cavity of water structure, which is formed by clathratelike hydrogen-bonded water network. On the basis of the band shape of solute-affected HDO spectra and the corresponding distribution of intermolecular distances, the criterion for hydrophobic type hydration has been proposed. From that point of view, all the studied solutes could be treated as hydrophobic ones. The limiting band position and the corresponding intermolecular distance of affected water, gained with increasing electron-donating power of solutes, has been inferred from the data obtained. These observations are important for interpretation of vibrational spectra of water as well as for volumetric measurements of solutions. The simple model of hydration, proposed to better justify the results, connects the values obtained from the methods providing microscopic and macroscopic characteristics of the system studied.     

Volumes of mixtures of alcohols in water at 25°C

Many studies support the idea that alcohols in water undergo microphase transitions which are, in many respects, similar to micellization. To investigate the interactions in these systems even further, the volumes of transfer of normal alcohols of intermediate chain length, kept near infinite dilution, were measured from water to aqueous solutions of 2-propanol and 2-butoxyethanol. These results were compared with the volumes of transfer of the same alcohols to aqueous solutions of octyldimethylamine oxide, a well-characterized non-ionic surfactant. The trends observed are all very similar, exhibiting in many cases a maximum in the transition region. This tends to confirm the formation of mixed aggregates in aqueous mixtures of alcohols, but, in a general way, it is also shown that the magnitude of an extremum in the functions of transfer is related to the relative hydrophobicities of the present solutes, the extremum appearing in most cases only when the transferred solute is more hydrophobic than the main solute.     

Characterization of solvation properties of lipid bilayer membranes in liposome electrokinetic chromatography

The nature of solute interactions with biomembrane-like liposomes, made of naturally occurring phospholipids and cholesterol, was characterized using electrokinetic chromatography (EKC). Liposomes were used as a pseudo-stationary phase in EKC that provided sites of interactions for uncharged solutes. The retention factors of uncharged solutes in liposome EKC are directly proportional to their liposome-water partition coefficients. Linear solvation energy relationship (LSER) models were developed to unravel the contributions from various types of interactions for solute partitioning into liposomes. Size and hydrogen bond acceptor strength of solutes are the main factors that determine partitioning into lipid bilayers. This falls within the general behavior of solute partitioning from an aqueous into organic phases such as octanol and micelles. However, there exist subtle differences in the solvation properties of liposomes as compared to those of octanol and various micellar pseudo-phases such as aggregates of sodium dodecyl sulfate (SDS), sodium cholate (SC), and tetradecylammonium bromide (TTAB). Among these phases, the SDS micelles are the least similar to the liposomes, while octanol, SC, and TTAB micelles exhibit closer solvation properties. Subsequently, higher correlations are observed between partitioning into liposomes and the latter three phases than that into SDS.     

Comparison of retention of aromatic hydrocarbons with polar groups in binary reversed-phase high-performance liquid chromatography systems

The retention of aromatic hydrocarbons with polar groups has been correlated as log k1 versus log k2 for reversed-phase high-performance liquid chromatography systems with different binary aqueous mobile phases containing methanol, acetonitrile or tetrahydrofuran as modifiers. Distinct changes in separation selectivity have been observed between tetrahydrofuran and acetonitrile or methanol systems. Methanol and acetonitrile systems show lower diversity of separation selectivity. The changes in retention and selectivity of aromatic hydrocarbons with various polar groups between any two chromatographic systems with binary aqueous eluents (tetrahydrofuran vs. acetonitrile, tetrahydrofuran vs. methanol and methanol vs. acetonitrile) have been interpreted in terms of molecular interactions of the solute with especially one component of the stationary phase region, i.e. extracted modifier, and stationary phase ordering. The ordering of the stationary phase region caused by modifier type influences the chromatographic selectivity of solutes with different molecular shape.     

Effects of nonpolar solutes on the thermodynamic response functions of aqueous mixtures

We investigate the effect of adding nonpolar solutes at atmospheric pressure on water's temperature of maximum density, isothermal compressibility, and isobaric heat capacity, using a statistical mechanical model of water solutions [H. S. Ashbaugh, T. M. Truskett, and P. G. Debenedetti, J. Chem. Phys. 116, 2907 (2002)]. We find that the temperature of maximum density increases with solute hydrophobicity, as characterized by its size, and decreases with its van der Waals attractive parameter a, in agreement with experiment. We predict similar trends for the addition of solutes on the isothermal compressibility and isobaric heat capacity: solute hydrophobicity causes an upward shift in water's anomalies, whereas dispersive interactions as measured by the solute's van der Waals attractive parameter shift the anomalies to lower temperatures. The locus along which the competing contributions of solute size sigma and interaction strength a to the shift in water's response functions balance each other obeys the scaling relationship sigma6 approximately a.     

Thermodynamic properties of ternary non-electrolyte solutions: Part 2. Prediction of infinite dilution activity coefficients and solubilities based on the various nibs and microscopic partition models

The limitations and applications of the various nearly ideal binary solvent (NIBS) and microscopic partition (MP) models for predicting the thermochemical properties of solutes dissolved in binary solvent mixtures are re-examined using published solute solubility and infinite dilution activity coefficient data for 48 systems. Expressions derived from the basic NIBS and extended NIBS models provide very reasonable predictions for anthracene and carbazole solubilities in systems containing both non-specific and specific interactions. For many of the systems considered, deviations between experimental values and NIBS predictions are of the order of 6% or less. In comparison, the two expressions derived from the MP model grossly overpredict the observed solubilities, with the average absolute deviations for several of the carbazole systems being 30% or more.     

Thermodynamics of mixed micellar systems of sodium octanoate and ethoxylated alcohols in water at 25°C

Thermodynamic volume and compressibility properties of binary aqueous systems of sodium octanoate (C₈Na) and polyoxyethylene butyl ether compounds (ethoxylated alcohols), with one to three oxyethylene groups (C₄EO_X), and ternary systems of these compounds, have been determined as a function of surfactant and alcohol concentrations at 25°C from density and ultrasonic velocity measurements. Values of standard state apparent molar volume and apparent molar adiabatic compressibility properties of transfer of ethoxylated alcohols from water to aqueous C₈Na solutions were obtained. The observed trends in the properties of the aqueous binary C₈Na systems and of the transfer functions for C₄EO_x, at a given low concentration of 0.05m, were analyzed using theoretical models. Good agreement between the simulated results and the experimental data was achieved. The results obtained for the distribution coefficients of the alcohols and the thermodynamic properties of both components of the mixed micelle show that the hydrophobicity of the ethoxylated alcohols is enhanced by introducing more oxyethylene groups into the alcohol. The distribution behavior of these compounds between the C₈Na micelles and the aqueous phase also depends on the difference between the hydrophobicity of the surfactant and the alcohols.     

CORRELATION BETWEEN INTERFACIAL TENSIONS AND MICELLAR STRUCTURES

It is shown that results of surface and interfacial tension measurements can be used to predict the type of micelles and of liquid crystalline phases which are formed in binary and ternary surfactant solutions. In particular it is possible to predict the position of l.c. cubic phases in ternary systems consisting of surfactant, hydrocarbon and water. Data to demonstrate the conclusions were obtained on the surfactants Alkyltrimethylammoniumbromides, Alkyldimethylaminoxides and Alkyldimethylphosphinoxides. It was found that the interfacial tension of a dilute micellar solution against a reference hydrocarbon is a most sensitive and indicative parameter for the prediction of the different structures. Large changes of the interfacial tension were observed for the three systems having the same hydrocarbon chainlength. The value of the interfacial tension directly reflects also the amount of hydrocarbon which can be solubilized in the micellar solution. Interfacial tensions larger than 1mN/m are indicative of globular micelles while interfacial tensions between 0.1 and 1 mN/m indicate the formation of rods. Values below 0.1 mN/m indicate disclike micelles or lamellar phases.

The interfacial tension depends somewhat on the kind of hydrocarbon which is used for the measurements. It is observed that for several surfactant solutions the interfacial tension passes through a shallow minimum when the chainlength of the hydrocarbon is increased from six to sixteen.     

Fluid-phase behavior of binary mixtures in which one component can have two critical points

We investigate theoretically the binary fluid-phase behavior of mixtures in which one water-like component can have two critical points. We consider three equal-sized nonpolar solutes that differ in the strength of their dispersive interactions (a1 < a2 < a3, where a denotes the van der Waals attractive parameter). In each case, we compare the phase behavior predicted using two sets of parameters for water: one giving rise to a pure component low-temperature liquid-liquid transition terminating at a critical point (two-critical-point parameter set), and one in which no such second critical point exists (singularity-free parameter set). Regardless of the parameter values used, we find five mixture critical lines. Using the two-critical-point parameter set, we find that a critical line originates at water's second critical point for aqueous mixtures involving solutes 1, 2, or 3. For mixtures involving solutes 1 or 2, this line extends towards low pressures and high temperatures as the solute mole fraction increases, and is closely related to the critical line originating at water's ordinary vapor-liquid critical point: these two critical lines are loci of upper and lower consolute points corresponding to the same liquid-liquid transition. In mixtures involving solute 2, the critical locus emanating from water's second critical point is shifted to higher temperatures compared to mixtures involving solute 1, and extends up to T approximately 310 K at moderate pressures (ca. 200 bars). This suggests the possibility of an experimentally accessible manifestation of the existence of a second critical point in water. For binary mixtures involving solutes 1 or 2, changing the water parameters from the two critical points to the singularity-free case causes the disappearance of a lower consolute point at moderate pressures. For binary mixtures involving solute 3, the differences between two-critical-point and singularity-free behaviors occur only in the experimentally difficult-to-probe low-temperature and high-pressure region.     

Inferences on dynamic water structure in aqueous solutions from diffusion measurements

The dynamic, hydrogen-bonded structure of water can be profoundly affected by addition of solutes as reflected by the resulting solute—solvent interactions. Measurements of diffusion coefficients for the solute species and for water are a useful probe for studying those interactions and changes in the dynamic water structure.

The analysis of transport properties of electrolyte solutions is presently approached both theoretically and experimentally by the use of generalised transport coefficients such as Onsager phenomenological coefficients, velocity correlation coefficients or the closely related distinct diffusion coefficients of Friedman. However, two of these approaches generally use a solvent-fixed frame of reference and thereby exclude the information available from studying the diffusion of the solvent.

The velocity correlation coefficient approach is used here to examine hydrogen-bonding effects in binary solution containing water as one component. We discuss also examples where intra-diffusion data for water in aqueous solutions give important insights into the dynamic structure of aqueous solutions.     

Correlation and comparison of the infinite dilution activity coefficients in aqueous and organic mixtures from a modified excess Gibbs energy model

Activity coefficients at infinite dilution (ln γ^∞) of aqueous systems were calculated using a modified excess Gibbs energy model. More than 95 binary systems with 15 solute families from nonpolar alkanes to polar alcohols and acids were employed in this study. Based on the local composition lattice model developed by Aranovich and Donohue, a modified excess Gibbs energy equation (m-AD model) was derived in this study. With two generalized parameters for each homologous series of solutes, this modified model yields satisfactory results for the limiting activity coefficients. The overall absolute average deviation (AAD) of ln γ^∞ for all aqueous systems investigated in this study is 2% for the m-AD model, and the corresponding AAD in γ^∞ unit is 7%. The calculated infinite dilution activity coefficients from the m-AD model are comparable to those from the MOSCED, SPACE, PDD, LSER or the modified UNIFAC model. The m-AD model shows lower peak deviation than those from other methods. Satisfactory generalized correlation results are also observed for organic solvents other than water. With the generalized parameters, the m-AD model satisfactorily predicts the limiting activity coefficients for other solutes not included in the correlation.     

A molecular dynamics study of the dielectric properties of aqueous solutions of alanine and alanine dipeptide

Molecular dynamics simulations were used to compute the frequency-dependent dielectric susceptibility of aqueous solutions of alanine and alanine dipeptide. We studied four alanine solutions, ranging in concentration from 0.13-0.55 mol/liter, and two solutions of alanine dipeptide (0.13 and 0.27 mol/liter). In accord with experiment we find a strong dielectric increment for both solutes, whose molecular origin is shown to be the zwitterionic nature of the solutes. The dynamic properties were analyzed based on a dielectric component analysis into solute, a first hydration shell, and all remaining (bulk) waters. The results of this three component decomposition were interpreted directly, as well as by uniting the solute and hydration shell component to a "suprasolute" component. In both approaches three contributions to the frequency-dependent dielectric properties can be discerned. The quantitatively largest and fastest component arises from bulk water [i.e., water not influenced by the solute(s)]. The interaction between waters surrounding the solute(s) (the hydration shell) and bulk water molecules leads to a relaxation process occurring on an intermediate time scale. The slowest relaxation process originates from the solute(s) and the interaction of the solute(s) with the first hydration shell and bulk water. The primary importance of the hydration shell is the exchange of shell and bulk waters; the self-contribution from bound water molecules is comparatively small. While in the alanine solutions the solute-water cross-terms are more important than the solute self-term, the solute contribution is larger in the dipeptide solutions. In the latter systems a much clearer separation of time scales between water and alanine dipeptide related properties is observed. The similarities and differences of the dielectric properties of the amino acid/peptide solutions studied in this work and of solutions of mono- and disaccharides and of the protein ubiquitin are discussed.     

Ermod: Fast and versatile computation software for solvation free energy with approximate theory of solutions

ERmod is a software package to efficiently and approximately compute the solvation free energy using the method of energy representation. Molecular simulation is to be conducted at two condensed‐phase systems of the solution of interest and the reference solvent with test‐particle insertion of the solute. The subprogram ermod in ERmod then provides a set of energy distribution functions from the simulation trajectories, and another subprogram slvfe determines the solvation free energy from the distribution functions through an approximate functional. This article describes the design and implementation of ERmod, and illustrates its performance in solvent water for two organic solutes and two protein solutes. Actually, the free‐energy computation with ERmod is not restricted to the solvation in homogeneous medium such as fluid and polymer and can treat the binding into weakly ordered system with nano‐inhomogeneity such as micelle and lipid membrane. ERmod is available on web at http://sourceforge.net/projects/ermod . © 2014 Wiley Periodicals, Inc.     

Ions in water: the microscopic structure of concentrated hydroxide solutions

Neutron-diffraction data on aqueous solutions of hydroxides, at solute concentrations ranging from 1 solute per 12 water molecules to 1 solute per 3 water molecules, are analyzed by means of a Monte Carlo simulation (empirical potential structure refinement), in order to determine the hydration shell of the OH- in the presence of the smaller alkali metal ions. It is demonstrated that the symmetry argument between H+ and OH- cannot be used, at least in the liquid phase at such high concentrations, for determining the hydroxide hydration shell. Water molecules in the hydration shell of K+ orient their dipole moment at about 45 degrees from the K+-water oxygen director, instead of radially as in the case of the Li+ and Na+ hydration shells. The K+-water oxygen radial distribution function shows a shallower first minimum compared to the other cation-water oxygen functions. The influence of the solutes on the water-water radial distribution functions is shown to have an effect on the water structure equivalent to an increase in the pressure of the water, depending on both ion concentration and ionic radius. The changes of the water structure in the presence of charged solutes and the differences among the hydration shells of the different cations are used to present a qualitative explanation of the observed cation mobility.     

Comparison of the retention characteristics of different pseudostationary phases for microemulsion and micellar electrokinetic chromatography of betamethasone and derivatives

Five electrokinetic chromatography systems were compared concerning retention behavior and lipophilicity. Comparison was based on capacity (retention) factors of some steroidal drugs, and on log P(OW) values derived by the aid of reference substances. In all systems the aqueous buffer consisted of phosphate (20 mM, pH 7.5). Two systems had micelles, three systems microdroplets as negatively charged pseudostationary phases. The micelles were formed by sodium dodecyl sulfate (SDS) and sodium cholate, respectively. One microemulsion consisted (as usual) from octane as oil, butanol as cosurfactant and SDS as charged tenside. Two microemulsions were made from biosurfactants (phosphatidylcholine, isopropylmyristate) to better simulate biopartitioning of the drugs. Even for noncharged analytes a change in migration sequence and thus in log P(OW) was observed for the systems consisting of the biosurfactants, compared to the others. For the former systems, log P(OW) derived from the capacity factors agree for all analytes with those obtained from calculation by computer software based on the structure of the drugs, and with experimental data directly obtained from octanol/water partitioning.     

Diffusion Coefficients of Three Organic Solutes in Aqueous Sodium Dodecyl Sulfate Solutions

Tracer diffusion coefficients of phenol, toluene, and benzoic acid in aqueous solutions of sodium dodecyl sulfate (SDS) were measured by the Taylor dispersion technique. In addition, the viscosities and densities of the SDS solutions were measured. For phenol and toluene, the effect of micelle formation on the diffusion coefficient is pronounced. When the SDS concentration is below the critical micelle concentration (cmc), the diffusion coefficients are almost independent of the SDS concentration. However, above the cmc there is a rapid decrease in the diffusion coefficients, and the apparent diffusion coefficients of the two solutes are the weighted average of free solute diffusion and the micelle diffusion. A model is presented to describe the diffusion behavior of the two solutes in aqueous micellar solutions of SDS. The interaction between the two solutes and the micelles has been investigated and the fraction of each solute that is solubilized by the micelles is estimated from the measured apparent diffusion coefficient. For benzoic acid, the diffusion coefficient is dependent on the joint contribution of the benzoic acid molecules that are solubilized by the micelles as well as the corresponding benzoate ions. The effect of micelle formation on the diffusion coefficient of benzoic acid is not as pronounced as for phenol and toluene. Copyright 2000 Academic Press.