Thermodynamic stability of fluid-fluid phase separation in binary athermal mixtures: the role of nonadditivity

We studied the thermodynamic stability of fluid-fluid phase separation in binary nonadditive mixtures of hard-spheres for moderate size ratios. We are interested in elucidating the role played by small amounts of nonadditivity in determining the stability of fluid-fluid phase separation with respect to the fluid-solid phase transition. The demixing curves are built in the framework of the modified-hypernetted chain and of the Rogers-Young integral equation theories through the calculation of the Gibbs free energy. We also evaluated fluid-fluid phase equilibria within a first-order thermodynamic perturbation theory applied to an effective one-component potential obtained by integrating out the degrees of freedom of the small spheres. A qualitative agreement emerges between the two different approaches. We also addressed the determination of the freezing line by applying the first-order thermodynamic perturbation theory to the effective interaction between large spheres. Our results suggest that for intermediate size ratios a modest amount of nonadditivity, smaller than earlier thought, can be sufficient to drive the fluid-fluid critical point into the thermodinamically stable region of the phase diagram. These findings could be significant for rare-gas mixtures in extreme pressure and temperature conditions, where nonadditivity is expected to be rather small.     

On the anatomy of the adsorption heat versus loading as a function of temperature and adsorbate for a graphitic surface

In this paper we review and classify the various patterns of isosteric heat versus loading for adsorption of gases on graphitised thermal carbon black at temperatures ranging from below the 3D triple point to temperatures above it, but less than the 3D critical point. We have identified the features of heat curve and highlighted the microscopic origin of these features. The patterns vary with temperature and with the relative strength of the fluid-fluid interaction and solid-fluid interaction. For simple adsorptives (by simple we meant there is no strong association between fluid particles), the heat curve is typified by fluid-fluid attraction and layering phenomena. For adsorptives showing strong association such as water, ammonia and methanol, the heat curve essentially begins below the condensation heat and then approaches it as loading is increased. This is mainly due to the strong hydrogen bonding in these fluids. A third group includes adsorptives such as benzene, where the heat curve is constant in the sub-monolayer coverage region (but is higher than the condensation heat) and then drops off to the condensation heat when higher layers are formed. The constant heat in the sub-monolayer region is due to the balance between the energy factor (from fluid-fluid interaction) and entropy factor (due to re-orientation of adsorbed molecules). Our proposed classification is supported by detailed GCMC simulations of various gases that have been reported in the literature, and we supplement these with new results for the adsorption of xenon on graphite to investigate in more detail the change in heat pattern with temperature. Xenon is chosen because of its high fluid-fluid interaction, allowing us to study the 2D-phase transition in the first as well as higher layers.     

Arrest of fluid demixing by nanoparticles: a computer simulation study

We use lattice Boltzmann simulations to investigate the formation of arrested structures upon demixing of a binary solvent containing neutrally wetting colloidal particles. Previous simulations for symmetric fluid quenches pointed to the formation of "bijels": bicontinuous interfacially jammed emulsion gels. These should be created when a glassy monolayer of particles forms at the fluid-fluid interface, arresting further demixing and rigidifying the structure. Experimental work has broadly confirmed this scenario, but it shows that bijels can also be formed in volumetrically asymmetric quenches. Here, we present new simulation results for such quenches, compare these to the symmetric case, and find a crossover to an arrested droplet phase at strong asymmetry. We then make extensive new analyses of the postarrest dynamics in our simulated bijel and droplet structures, on time scales comparable to the Brownian time for colloid motion. Our results suggest that, on these intermediate time scales, the effective activation barrier to ejection of particles from the fluid-fluid interface is smaller by at least 2 orders of magnitude than the corresponding barrier for an isolated particle on a flat interface.     

Phase behavior of short-range square-well model

Various Monte Carlo techniques are used to determine the complete phase diagrams of the square-well model for the attractive ranges lambda = 1.15 and lambda = 1.25. The results for the latter case are in agreement with earlier Monte Carlo simulations for the fluid-fluid coexistence curve and yield new results for the liquidus-solidus lines. Our results for lambda = 1.15 are new. We find that the fluid-fluid critical point is metastable for both cases, with the case lambda = 1.25 being just below the threshold value for metastability. We compare our results with prior studies and with experimental results for the gamma(II)-crystallin.     

Effect of surface-perturbed intermolecular interaction on adsorption of simple gases on a graphitized carbon surface

In this paper, we investigate the effect of the solid surface on the fluid-fluid intermolecular potential energy. This modified fluid-fluid interaction energy due to the inducement of a solid surface is used in the grand canonical Monte Carlo (GCMC) simulation of various noble gases, nitrogen, and methane on graphitized thermal carbon black. This effect is such that the effective interaction potential energy between two particles close to surface is less than the potential energy if the solid substrate is not present. With this modification the GCMC simulation results agree extremely well with the experimental data over a wide range of pressures while the simulation results with the unmodified potential energy give rise to a shoulder near the neighborhood of monolayer coverage and the significant overprediction of the second and higher layer coverages. In particular the unmodified GCMC results exhibit very sharp change in those higher layers while the experimental data have a much gradual change in the uptake. We will illustrate this theory with adsorption data of argon, xenon, neon, nitrogen, and methane on graphitized thermal carbon black.     

Surface tension of the Widom-Rowlinson model

We consider the computation of the surface tension of the fluid-fluid interface for the Widom-Rowlinson [J. Chem. Phys. 52, 1670 (1970)] binary mixture from direct simulation of the inhomogeneous system. We make use of the standard mechanical route, in which the surface tension follows from the computation of the normal and tangential components of the pressure tensor of the system. In addition to the usual approach, which involves simulations of the inhomogeneous system in the canonical ensemble, we also consider the computation of the surface tension in an ensemble where the pressure perpendicular (normal) to the planar interface is kept fixed. Both approaches are seen to provide consistent values of the interfacial tension. The issue of the system-size dependence of the surface tension is addressed. In addition, simulations of the fluid-fluid coexistence properties of the mixture are performed in the semigrand canonical ensemble. Our results are compared with existing data of the Widom-Rowlinson mixture and are also examined in the light of the vapor-liquid equilibrium of the thermodynamically equivalent one-component penetrable sphere model.     

Many-body interactions between particles in a polydisperse polymer fluid

We use a continuum chain model and develop an analytical theory for the interaction between many spheres immersed in a fluid of ideal polydisperse polymers. Assuming local spherical symmetry of the polymer field about each particle, combined with a local approximation, compact expressions are derived for the many-body interaction between the spheres. We use a mean-field approximation to investigate the fluid-fluid phase diagram for the mixture.     

Structure, phase behavior, and inhomogeneous fluid properties of binary dendrimer mixtures

The effective pair potentials between different kinds of dendrimers in solution can be well approximated by appropriate Gaussian functions. We find that in binary dendrimer mixtures the range and strength of the effective interactions depend strongly upon the specific dendrimer architecture. We consider two different types of dendrimer mixtures, employing the Gaussian effective pair potentials, to determine the bulk fluid structure and phase behavior. Using a simple mean field density functional theory (DFT) we find good agreement between theory and simulation results for the bulk fluid structure. Depending on the mixture, we find bulk fluid-fluid phase separation (macrophase separation) or microphase separation, i.e., a transition to a state characterized by undamped periodic concentration fluctuations. We also determine the inhomogeneous fluid structure for confinement in spherical cavities. Again, we find good agreement between the DFT and simulation results. For the dendrimer mixture exhibiting microphase separation, we observe a rather striking pattern formation under confinement.     

Phase behavior of aqueous solutions containing dipolar proteins from second-order perturbation theory

Due to the interplay of Coulombic repulsion and attractive dipolar and van der Waals interactions, solutions of globular proteins display a rich variety of phase behavior featuring fluid-fluid and fluid-solid transitions that strongly depend on solution pH and salt concentration. Using a simple model for charge, dispersion and dipole-related contributions to the interprotein potential, we calculate phase diagrams for protein solutions within the framework of second-order perturbation theory. For each phase, we determine the Helmholtz energy as the sum of a hard-sphere reference term and a perturbation term that reflects both the electrostatic and dispersion interactions. Dipolar effects can induce fluid-fluid phase separation or crystallization even in the absence of any significant dispersion attraction. Because dissolved electrolytes screen the charge-charge repulsion more strongly than the dipolar attraction, the ionic strength dependence of the potential of mean force can feature a minimum at intermediate ionic strengths offering an explanation for the observed nonmonotonic dependence of the phase behavior on salt concentration. Inclusion of correlations between charge-dipole and dipole-dipole interactions is essential for a reliable calculation of phase diagrams for systems containing charged dipolar proteins and colloids.     

Stability of a fluid-fluid interface in a biconical pore segment

Pore networks that include biconical pore segments are frequently used to model two-phase flow. In this work, we describe in detail the displacement of a fluid-fluid interface in such a pore segment. We assume sharp edges in the throat, inlet, and outlet of the pore segment to be the limiting cases of round edges, the radii of which vanish. We account for interfacial and lineal tensions that cause nonconstant contact angles. For zero lineal tension, we provide analytical solutions for flow induced by changing infinitesimally slowly either capillary pressure or the volume of one fluid. In diverging and converging cones, the common line among the two fluids and the solid phase slides while it is pinned in the throat, inlet, and outlet. We observe hysteresis within the pore segment, and drainage entry pressures deviate from prior work.     

Acoustofluidics 16: acoustics streaming near liquid-gas interfaces: drops and bubbles

In this sixteenth part of the series on "Acoustofluidics-exploiting ultrasonic standing waves forces and acoustic streaming in microfluidic systems for cell and particle manipulation," we continue our discussion on the analytical aspects of the streaming phenomenon. In particular, the use of the singular perturbation technique for this class of problems is delineated with a set of examples where fluid-fluid interaction takes place. In this category, we focus on drops and bubbles, and deal specifically with the effect of interfacial mobility on the streaming flow.     

Phase diagram of the adhesive hard sphere fluid

The phase behavior of the Baxter adhesive hard sphere fluid has been determined using specialized Monte Carlo simulations. We give a detailed account of the techniques used and present data for the fluid-fluid coexistence curve as well as parametrized fits for the supercritical equation of state and the percolation threshold. These properties are compared with the existing results of Percus-Yevick theory for this system.     

Effect of solvent on the phase diagram of a simple anisotropic model of globular proteins

Anisotropic protein interactions play a role in globular protein phase transitions, including the shape of fluid-fluid coexistence curves and the formation of hemoglobin polymer fibers in sickle cell disease. Also, the solvent has been shown to play an important role in the phase behavior of some aqueous protein solutions, through the release or trapping of water molecules upon crystallization. Both anisotropy and solvent effects have been treated separately in earlier theoretical studies. Here we propose and analyze a simple, composite model that treats both anisotropy and solvent effects. We find that this model qualitatively explains some phase behavior that has been observed in protein solutions, including normal and retrograde solubility curves and lower critical points.     

A finite-size scaling study of a model of globular proteins

Grand canonical Monte Carlo simulations are used to explore the metastable fluid-fluid coexistence curve of the modified Lennard-Jones model of globular proteins of ten Wolde and Frenkel [Science, 277, 1975 (1997)]. Using both mixed-field finite-size scaling and histogram-reweighting methods, the joint distribution of density and energy fluctuations is analyzed at coexistence to accurately determine the critical-point parameters. The subcritical coexistence region is explored using the recently developed hyper parallel tempering Monte Carlo simulation method along with histogram reweighting to obtain the density distributions. The phase diagram for the metastable fluid-fluid coexistence curve is calculated in close proximity to the critical point, a region previously unattained by simulations.     

Tuning properties of silver particle monolayers via controlled adsorption-desorption processes

Nanoparticle self-assembly at fluid-fluid interfaces has been traditionally exploited in emulsification, encapsulation and oil recovery, and more recently in emerging applications including functional nanomaterials and biphasic catalysis. We provide a review of the literature focusing on the open challenges that still hamper the broader applicability of this potentially transformative technology, and we outline strategies to achieve improved control over interfacial self-assembly of nanoparticles. First, we discuss means to promote spontaneous adsorption by tuning the interfacial energies of the nanoparticles with the fluids using capping ligands, and the occurrence of energy barriers. We then examine the interactions between interfacial nanoparticles and how they affect the formation of equilibrium interfacial suspensions versus non-equilibrium two-dimensional phases, such as weakly attractive glasses and gels. Important differences with colloidal interactions in a bulk suspension arise due to the discontinuity in solvent properties at the interface. For instance, ligand brushes rearrange in asymmetric configurations, and thus play a significant role in determining interparticle interactions. Finally, we briefly discuss the link between interfacial microstructure and the dynamic response of particle-laden interfaces, including interfacial rheology and the fate of nanoparticle monolayers upon out-of-plane deformation.     

On the heat capacity of adsorbed phases using molecular simulation

The heat capacities of argon, ammonia, and methanol on carbon black at 87.3, 240, and 300 K, respectively, have been investigated. The carbon black surface has been modeled with and without carbonyl groups. Part of this investigation is a decomposition of the heat capacity into its contributions from the different interaction potentials of an adsorption system. All systems show a spectrum of heat capacity versus loading, and this behavior depends on the carbonyl configuration present on the surface. For methanol and ammonia the variation of the heat capacity between the two for the same carbonyl configurations is greater than the variation in the heat of adsorption. Heat capacities of methanol and ammonia are generally dominated by fluid-fluid interactions due to the strong association of fluid particles through hydrogen bonding. The difference in the heat capacity behavior of the two fluids is an indicator of their different clustering behaviors on the carbon black surface. The presence of carbonyl groups reduces the fluid-fluid contributions to the heat capacity. This is due to the compensation of fluid-fluid interactions with fluid-functional group interactions. At 87.3 K a first layer transition to a solidlike state is present for argon and results in a large peak in the heat capacity on a bare surface. The presence of functional groups greatly reduces this peak in the heat capacity by disrupting the packing of argon on the surface and preventing a transition to a solidlike state.     

Solvent mediated interactions close to fluid-fluid phase separation: microscopic treatment of bridging in a soft-core fluid

Using density functional theory we calculate the density profiles of a binary solvent adsorbed around a pair of big solute particles. All species interact via repulsive Gaussian potentials. The solvent exhibits fluid-fluid phase separation, and for thermodynamic states near to coexistence the big particles can be surrounded by a thick adsorbed "wetting" film of the coexisting solvent phase. On reducing the separation between the two big particles we find there can be a "bridging" transition as the wetting films join to form a fluid bridge. The effective (solvent mediated) potential between the two big particles becomes long ranged and strongly attractive in the bridged configuration. Within our mean-field treatment the bridging transition results in a discontinuity in the solvent mediated force. We demonstrate that accounting for the phenomenon of bridging requires the presence of a nonzero bridge function in the correlations between the solute particles when our model fluid is described within a full mixture theory based upon the Ornstein-Zernike equations.