Statistical Thermodynamics of Solutions

A simple theory for interfacial properties of dilute solutions

Recent studies suggest that cosolute mixtures may exert significant non-additive effects upon protein stability. The corresponding liquid–vapor interfaces may provide useful insight into these non-additive effects. Accordingly, in this work, we relate the interfacial properties of dilute multicomponent solutions to the interactions between solutes. We first derive a simple model for the surface excess of solutes in terms of thermodynamic observables. We then develop a lattice-based statistical mechanical perturbation theory to derive these observables from microscopic interactions. Rather than adopting a random mixing approximation, this dilute solution theory (DST) exactly treats solute–solute interactions to lowest order in perturbation theory. Although it cannot treat concentrated solutions, Monte Carlo (MC) simulations demonstrate that DST describes the interactions in dilute solutions with much greater accuracy than regular solution theory. Importantly, DST emphasizes a fundamental distinction between the “intrinsic” and “effective” preferences of solutes for interfaces. DST predicts that three classes of solutes can be distinguished by their intrinsic preference for interfaces. While the surface preference of strong depletants is relatively insensitive to interactions, the surface preference of strong surfactants can be modulated by interactions at the interface. Moreover, DST predicts that the surface preference of weak depletants and weak surfactants can be qualitatively inverted by interactions in the bulk. We also demonstrate that DST can be extended to treat surface polarization effects and to model experimental data. MC simulations validate the accuracy of DST predictions for lattice systems that correspond to molar concentrations.

Entropy determination for mixtures in the adiabatic grand-isobaric ensemble

The entropy change that occurs upon mixing two fluids has remained an intriguing topic since the dawn of statistical mechanics. In this work, we generalize the grand-isobaric ensemble to mixtures and develop a Monte Carlo algorithm for the rapid determination of entropy in these systems. A key advantage of adiabatic ensembles is the direct connection they provide with entropy. Here, we show how the entropy of a binary mixture A–B can be readily obtained in the adiabatic grand-isobaric ( μA, μB, P, R) ensemble, in which μA and μB denote the chemical potential of components A and B, respectively, P is the pressure, and R is the heat (Ray) function, that corresponds to the total energy of the system. This, in turn, allows for the evaluation of the entropy of mixing and the Gibbs free energy of mixing. We also demonstrate that our approach performs very well both on systems modeled with simple potentials and with complex many-body force fields. Finally, this approach provides a direct route to the determination of the thermodynamic properties of mixing and allows for the efficient detection of departures from ideal behavior in mixtures.

Conformal Sites Theory for Adsorbed Films on Energetically Heterogeneous Surfaces

We present a conformal sites theory for a solid substrate whose surface is both geometrically and energetically heterogeneous and that interacts with an adsorbed film. The theory is based on a perturbation expansion for the grand potential of a real system with a rough surface about that of a reference system with an ideal reference surface, thus mapping the real system onto a much simpler interfacial system. The expansion is in powers of the intermolecular potential parameters, and leads to mixing rules for the potential parameters of the reference system. Grand canonical Monte Carlo simulations for the adsorption of argon at 87.3 K, carbon dioxide at 273 K, and water vapor at 298 K on heterogeneous carbon surfaces are investigated to explore the limits of applicability of the theory. Simulation results indicate that the theory works well with typical asymmetry of the potential parameters in the force field. However, care should be taken when applying the theory to strongly associating fluids and in the low-pressure region where the active surface sites play an important role. The conformal sites theory can be used to predict the adsorption properties and to characterize the solid substrate by taking advantage of the corresponding states principle. Other possible applications are also discussed.

The rehabilitation of irreversible processes and dissipative structures' 50th anniversary

In 2017, Ilya Prigogine would have been 100 years of age. As for any human being, this centenary is a notable event. For him, as a scientist, 2017 was also above all the 50th anniversary of dissipative structures It was indeed in 1967 that for the first time he used this denomination at the occasion of an important scientific event and in publications. The attribution of this qualification for self-organized behaviours of matter only possible far from equilibrium coincided with the outcome of a research effort of more than 25 years. Centred in thermodynamics and statistical physics on the role played by irreversible processes in the physical evolution of matter, the aim of this research is clear from the outset of his scientific career. With visionary personal intuition and iron-willed determination, it was pursued. The road to success had been long and sinuous, but finally it led to what he called the rehabilitation of irreversible processes The progresses that stand out as major landmarks of this endeavour that imposed a U-turn with respect to conceptions of classical physics deeply rooted since the nineteenth century will be described. This article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (Part 1)'.

Lattice model for water-solute mixtures

A lattice model for the study of mixtures of associating liquids is proposed. Solvent and solute are modeled by adapting the associating lattice gas (ALG) model. The nature of interaction of solute/solvent is controlled by tuning the energy interactions between the patches of ALG model. We have studied three set of parameters, resulting in, hydrophilic, inert, and hydrophobic interactions. Extensive Monte Carlo simulations were carried out, and the behavior of pure components and the excess properties of the mixtures have been studied. The pure components, water (solvent) and solute, have quite similar phase diagrams, presenting gas, low density liquid, and high density liquid phases. In the case of solute, the regions of coexistence are substantially reduced when compared with both the water and the standard ALG models. A numerical procedure has been developed in order to attain series of results at constant pressure from simulations of the lattice gas model in the grand canonical ensemble. The excess properties of the mixtures, volume and enthalpy as the function of the solute fraction, have been studied for different interaction parameters of the model. Our model is able to reproduce qualitatively well the excess volume and enthalpy for different aqueous solutions. For the hydrophilic case, we show that the model is able to reproduce the excess volume and enthalpy of mixtures of small alcohols and amines. The inert case reproduces the behavior of large alcohols such as propanol, butanol, and pentanol. For the last case (hydrophobic), the excess properties reproduce the behavior of ionic liquids in aqueous solution.

Solvation theory to provide a molecular interpretation of the hydrophobic entropy loss of noble-gas hydration

An equation for the chemical potential of a dilute aqueous solution of noble gases is derived in terms of energies, force and torque magnitudes, and solute and water coordination numbers, quantities which are all measured from an equilibrium molecular dynamics simulation. Also derived are equations for the Gibbs free energy, enthalpy and entropy of hydration for the Henry's law process, the Ostwald process, and a third proposed process going from an arbitrary concentration in the gas phase to the equivalent mole fraction in aqueous solution which has simpler expressions for the enthalpy and entropy changes. Good agreement with experimental hydration free energies is obtained in the TIP4P and SPC/E water models although the solute's force field appears to affect the enthalpies and entropies obtained. In contrast to other methods, the approach gives a complete breakdown of the entropy for every degree of freedom and makes possible a direct structural interpretation of the well-known entropy loss accompanying the hydrophobic hydration of small non-polar molecules under ambient conditions. The noble-gas solutes experience only a small reduction in their vibrational entropy, with larger solutes experiencing a greater loss. The vibrational and librational entropy components of water actually increase but only marginally, negating any idea of water confinement. The term that contributes the most to the hydrophobic entropy loss is found to be water's orientational term which quantifies the number of orientational minima per water molecule and how many ways the whole hydrogen-bond network can form. These findings help resolve contradictory deductions from experiments that water structure around non-polar solutes is similar to bulk water in some ways but different in others. That the entropy loss lies in water's rotational entropy contrasts with other claims that it largely lies in water's translational entropy, but this apparent discrepancy arises because of different coordinate definitions and reference frames used to define the entropy terms.

The second-order theory of conformal solutions

This paper describes a theoretical contribution to the statistical thermodynamics of mixtures of spherical molecules. The second-order perturbation free energy of a conformal solution is obtained by a rigorous Taylor-series expansion of the configuration integral in powers of the differences between intermolecular energy and size parameters, about an ideal unperturbed reference solution. Unlike the first-order terms, those of the second order contain statistical functions of the reference solution which cannot, in general, be related to its thermodynamic properties. All but one of these functions are concerned with departures from a random molecular distribution, and have been called molecular fluctuation integrals ; the remaining function can be related exactly to thermodynamic properties for the Lennard-Jones form of the intermolecular potential. The expressions for the molecular fluctuation integrals implied by the full random mixing approximation and by the semi-random mixing approximation of the cell theory, are derived and compared with the correct expressions given by the cell theory. The role of the Taylor series expansion as a critique of solution theories is discussed.