Thin-thick surface phase coexistence and boundary tension of the square-well fluid on a weak attractive surface

Prewetting transition is studied for the square-well fluid of attractive-well diameter lambda(ff)sigma(ff)=1.5 in the presence of a homogeneous surface modeled by the square-well potential of attractive well from 0.8sigma(ff) to 1.8sigma(ff). We investigate surface phase coexistence of thin-thick film transition using grand-canonical transition matrix Monte Carlo (GC-TMMC) and histogram reweighting techniques. Molecular dynamics (MD) and GC-TMMC are utilized to predict the properties of the fluid for various surface fluid affinities. Occurrences of prewetting transition with the variation of surface affinity are observed for a domain of reduced temperature from T(*)=0.62 to 0.75. We have used MD and GC-TMMC+finite size scaling (FSS) simulations to calculate the boundary tension as a function of temperature as well as surface affinity. Boundary tensions via MD and GC-TMMC+FSS methods are in good agreement. The boundary tension increases with the decrease of wall-fluid affinity. Prewetting critical properties are calculated using rectilinear diameter approach and scaling analysis. We found that critical temperature and density increase with the decrease of wall-fluid affinity.     

Monte Carlo simulation of vapor-liquid equilibrium and critical asymmetry of square-well dimer fluid

The critical behavior of square-well dimer fluid has been investigated using grand canonical ensemble Monte Carlo simulations combined with a histogram reweighting technique, hyper-parallel tempering and finite-size scaling. The critical temperature and density obtained are T(c)*=1.5495±0.0009 and ρ(c)*=0.1473±0.0007, which are 2.5% lower and 5.2% higher than previous results. Coexistence curves both near to and far from the critical point were obtained. The vapor-liquid equilibrium data far from the critical point are consistent with previous results. Simulation results show that the contribution of |t|(1-α) to the coexistence diameter of square-well dimer fluid dominates the critical behavior and the contribution of |t|(2β) is larger than for a hard-core square-well fluid.     

Molecular simulation study of the effect of pressure on the vapor-liquid interface of the square-well fluid

We examine a model system to study the effect of pressure on the surface tension of a vapor-liquid interface. The system is a two-component mixture of spheres interacting with the square-well (A-A) and hard-sphere (B-B) potentials and with unlike (A-B) interactions ranging (for different cases) from hard sphere to strongly attractive square well. The bulk-phase and interfacial properties are measured by molecular dynamics simulation for coexisting vapor-liquid phases for various mixture compositions, pressures, and temperatures. The variation of the surface tension with pressure compares well to values given by surface-excess formulas derived from thermodynamic considerations. We find that surface tension increases with pressure only for the case of an inert solute (hard-sphere A-B interactions) and that the presence of A-B attractions strongly promotes a decrease of surface tension with pressure. An examination of density and composition profiles is made to explain these effects in terms of surface-adsorption arguments.     

The calculation of vapor-liquid coexistence curve of Morse fluid: application to iron

The vapor-liquid coexistence curve of Morse fluid was calculated within the integral equations approach. The critical point coordinates were estimated. The parameters of Morse potential, fitted for elastic constants in solid phase, were used here to apply the results of present calculations to the determination of iron binodal. The properties of copper and sodium were considered in an analogous way. The calculations of pair correlation functions and isobars at liquid phase have shown that only for sodium these potential parameters allow one to obtain agreement with the measurements data. For iron another parameters are necessary to get this agreement in liquid phase. However, they give rise to very low critical temperature and pressure with respect to the estimates of other authors. Consequently, one can suppose that Morse potential is possibly inapplicable to the calculation of high temperature properties of non-alkali metals in disordered phases.     

Mechanical stability of micropipet-aspirated giant vesicles with fluid phase coexistence

Micropipet aspiration of phase-separated lipid bilayer vesicles can elucidate physicochemical aspects of membrane fluid phase coexistence. Recently, we investigated the composition dependence of line tension at the boundary between liquid-ordered and liquid-disordered phases of giant unilamellar vesicles obtained from ternary lipid mixtures using this approach. Here we examine mechanical equilibria and stability of dumbbell-shaped vesicles deformed by line tension. We present a relationship between the pipet aspiration pressure and the aspiration length in vesicles with two coexisting phases. Using a strikingly simple mechanical model for the free energy of the vesicle, we predict a relation that is in almost quantitative agreement with experiment. The model considers the vesicle free energy to be proportional to line tension and assumes that the vesicle volume, domain area fraction, and total area are conserved during aspiration. We also examine a mechanical instability encountered when releasing a vesicle from the pipet. We find that this releasing instability is observed within the framework of our model that predicts a change of the compressibility of a pipet-aspirated membrane cylinder from positive (i.e., stable) to negative (unstable) values, at the experimental instability. The model furthermore includes an aspiration instability that has also previously been experimentally described. Our method of studying micropipet-induced shape transitions in giant vesicles with fluid domains could be useful for investigating vesicle shape transitions modulated by bending stiffness and line tension.     

Surface tension of fluoroalkanes in a liquid phase

The parachor assigned to fluorine atoms in fluorinated alkanes was examined. The values vary systematically with the number of fluorine atoms in the molecule. Partially fluorinated alkanes show large values, while perfluoroalkanes give a smaller value which is constant for all the perfluoroalkanes examined in this work. The molecular structure seems to be important for the intermolecular interaction of this series of compounds, and may be responsible for the deviation from the parachor additivity rule.     

Effect of three-body interactions on the vapor-liquid phase equilibria of binary fluid mixtures

Gibbs-Duhem Monte Carlo simulations are reported for the vapor-liquid phase coexistence of binary argon+krypton mixtures at different temperatures. The calculations employ accurate two-body potentials in addition to contributions from three-body dispersion interactions resulting from third-order triple-dipole interactions. A comparison is made with experiment that illustrates the role of three-body interactions on the phase envelope. In all cases the simulations represent genuine predictions with input parameters obtained independently from sources other than phase equilibria data. Two-body interactions alone are insufficient to adequately describe vapor-liquid coexistence. In contrast, the addition of three-body interactions results in very good agreement with experiment. In addition to the exact calculation of three-body interactions, calculations are reported with an approximate formula for three-body interactions, which also yields good results.     

Direct determination of phase behavior of square-well fluids

We have combined Gibbs ensemble Monte Carlo simulations with the aggregation volume-biased method in conjunction with the Gibbs-Duhem method to provide the first direct estimates for the vapor-solid, vapor-liquid, and liquid-solid phase coexistences of square-well fluids with three different ranges of attraction. Our results are consistent with the previous simulations and verify the notion that the vapor-liquid coexistence behavior becomes metastable for cases where the attraction well becomes smaller than 1.25 times the particle diameter. In these cases no triple point is found.     

Surface confined ionic liquid as a stationary phase for HPLC

Trimethoxysilane "ionosilane" derivatives of room temperature ionic liquids based on alkylimidazolium bromides were synthesized for attachment to silica support material. The derivatives 1-methyl-3-(trimethoxysilylpropyl)imidazolium bromide and 1-butyl-3-(trimethoxysilylpropyl)imidazolium bromide were used to modify the surface of 3 microm diameter silica particles to act as the stationary phase for HPLC. The modified particles were characterized by thermogravimetric analysis (TGA) and (13)C and (29)Si NMR spectroscopies. The surface modification procedure rendered particles with a surface coverage of 0.84 micromol m(-2) for the alkylimidazolium bromide. The ionic liquid moiety was predominantly attached to the silica surface through two siloxane bonds of the ionosilane derivative (63%). Columns packed with the modified silica material were tested under HPLC conditions. Preliminary evaluation of the stationary phase for HPLC was performed using aromatic carboxylic acids as model compounds. The separation mechanism appears to involve multiple interactions including ion exchange, hydrophobic interaction, and other electrostatic interactions.     

Thermodynamics and phase transitions in a fluid confined by a harmonic trap

We study a fluid of interacting atoms confined by a three-dimensional anisotropic harmonic potential, similar to those produced by the magnetic traps used to confine cold atoms. We show that instead of the usual thermodynamic variables pressure and volume, no longer existing in this case, there appear "new" variables: the volume is replaced by (the inverse cube of) the geometric average of the oscillator frequencies of the trap, and the hydrostatic pressure is replaced by an intensive variable, conjugate to the previous one, and responsible for the mechanical equilibrium of the fluid in the trap. We discuss the origin and physical meaning of these new variables. With the aid of molecular dynamics simulations we show the emergence of novel liquid, vapor and solid-like phases in a classical fluid. In particular, we calculate the liquid-vapor-like coexistence curve and show evidence for the appearance of a critical point. These phase transitions should be observable in fluids of not-so-cold alkaline atoms.     

A study of square-well statistical associating fluid theory approximations

Six square-well (SW) statistical associating fluid theory (SAFT) models, fitted to the experimental saturated liquid volume and saturated vapor pressure for pure n-alkanes, are analyzed for predicting the coexisting densities, second virial coefficients, and binary phase equilibria. The models that result in low values of the segment energy and weak molecular weight dependence of the parameters are found to be more accurate for real fluids. The inclusion of the dimer structure in the SW chain term seems to produce no significant benefit for representing real substances.     

Surface tension of equilibrium spherical drops in the vapor phase

The properties of metastable and equilibrium drops that occur in the vapor phase and differ in the size and position of the dividing surface are compared. Using an equimolecular dividing surface, it was found that the total free energy and the total mass of the substance in drops of the same size over a broad temperature range differ by not more than 0.3%. A similar comparison for a dividing surface chosen from the equality condition of the moments of forces gives rather similar results. Using the surface where the maximum surface tension is attained as the dividing surface is impossible at low temperatures, because under these conditions, the notion of surface tension has no physical meaning. At high temperatures, the difference between the total mass of the substance for metastable and equilibrium drops does not exceed 0.6%. The calculations were carried out over a broad temperature range on the basis of the lattice-gas model in the quasichemical approximation.     

Surface tension and tolman length of spherical particulate in contact with fluid

The structure, surface tension and Tolman length of particulate-fluid interfaces were studied theoretically. Within the framework of density functional theory, the nonlocal, modified fundamental measure theory and direct correlation function from the first-order mean spherical approximation were incorporated. The theory accurately predicted the structure of fluid and the particulate-vapor surface tensions. The predictions of surface tensions for particulate-liquid interface and particulate in supercritical fluid are also reasonable. Especially, Tolman lengths for particulate-fluid interfaces were investigated systematically. The correct prediction of surface tension from Tolman length indicates that our analysis is reliable. Furthermore, Tolman length as a function of spherical particulate diameter, particulate-fluid interaction energy, and the properties of the fluid is fully discussed.     

Viscosity for a dense fluid of square-well rough spheres

Shear viscosity is calculated for a dense fluid of square-well rough spheres. The results are compared with calculated shear viscosities for the familiar hard-sphere, square-well and rough-sphere models.     

Density and chain conformation profiles of square-well chains confined in a slit by density-functional theory

Density and chain conformation profiles of square-well chains between two parallel walls were studied by using density-functional theory. The free energy of square-well chains is separated into two contributions: the hard-sphere repulsion and the attraction. The Heaviside function is used as the weighting function for both of the two parts. The equation of state of Hu et al. is used to calculate the excess free energy of the repulsive part. The equation of state of statistical associating fluid theory for chain molecules with attractive potentials of variable range [A. Gil-Villegas et al. J. Chem. Phys. 106, 4168 (1997)] is used to calculate the excess free energy of the attractive part. Because the wall is inaccessible to a mass center of a longer chain, there exists a sharp fall in the distribution of end-to-end distance near the wall as the chain length increases. When the average density of the system is not too low, the prediction of this work is in good agreement with computer simulation results for the density profiles and the chain conformation over a wide range of chain length, temperature, and attraction strength of the walls. However, when the average density and the temperature are very low, the prediction deviates to a certain degree from the computer simulation results for molecules with long chain length. A more accurate functional approximation is needed.     

Molecular dynamical calculations on the viscosity of a square-well fluid

The coefficient of viscosity for a square-well fluid is calculated by molecular dynamics as a function of the well-depth for densities up to the region of the fluid-solid phase transition. The inclusion of an attractive contribution in the intermolecular potential has a profound influence on the behaviour of the viscosity coeffient and is also responsible for the qualitative correspondence with real systems which has been found for densities above the critical one.     

Investigation of excess adsorption, solvation force, and plate-fluid interfacial tension for Lennard-Jones fluid confined in slit pores

The excess Helmholtz free energy functional is formulated in terms of a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)] for a short ranged repulsion and a first-order mean-spherical approximation theory [Y. P. Tang, J. Chem. Phys. 118, 4140 (2003)] for a long ranged attraction. Within the framework of the density functional theory, the density profile, excess adsorption, solvation force, and plate-fluid interfacial tension of a Lennard-Jones fluid confined in slit pores are predicted, and the results agree well with the simulation data. The phase equilibria inside the slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal, and the plate-fluid interfacial tensions at equilibrium states are predicted consequently.     

Molecular association of heteronuclear vibrating square-well dumbbells in liquid-vapor phase equilibrium

Molecular aggregates are formed by heteronuclear vibrating square-well dumbbells. In a recent article [G. A. Chapela and J. Alejandre, J. Chem. Phys., 132(10), 104704 (2010)], it is shown that heteronuclear vibrating square-well dumbbells with a diameter ratio between particles of 1/2 and interacting potential ratio of 4 form micelles of different sizes and shapes which manifest themselves in both the liquid and vapor phases, up to and above the critical point. This means that micellization and phase separation are present simultaneously in this simple model. These systems present a maximum in the critical temperature when plotted against the potential well depth of the second particle ε(2). In the same publication, it was speculated that the formation of micelles was responsible for the appearance of the maximum. A thorough study on this phenomena is presented here and it is found that there is a threshold on the size of the second particle and its corresponding depth of interaction potential, where the micelles are formed. If the diameter and well depth of the second particle are small enough for the first and deep enough for the second, micelles are formed. For σ(2)/σ(1) between 0.25 and 0.65 and ε(2)/ε(1) larger than 5.7, micelles are formed up to and above the critical temperature. Outside these ranges micelles appear only at temperatures lower than the critical point. There is a strong temperature dependence on the formation and persistence of the aggregates. For the deepest wells and large enough second particles, a gel interconnected aggregate is obtained. In this work, the micelles are formed at temperatures as low as the triple point and as high as the critical point and, in some cases, persist well above it. The presence of these maxima in critical temperatures T(c) when plotted against ε(2) as follows. At lower values of ε(2), an increase of T(c) is obtained as is expected by the increase of the attractive volume as indicated by the principle of corresponding states. As ε(2) increases further, the formation of molecular aggregates produce a saturation effect of the deepening of the potential well by encapsulating the particles of the second kind inside the micelles, so the resulting T(c) represents a new poly disperse system of molecular aggregates and not the original heteronuclear vibrating square-well dumbbells. The surface tension is also analyzed for these systems, and it is shown that decreases with increasing attraction due to the formation of molecular aggregates.     

Molecular simulation of the shear viscosity and the self-diffusion coefficient of mercury along the vapor-liquid coexistence curve

In earlier work [G. Raabe and R. J. Sadus, J. Chem. Phys. 119, 6691 (2003)] we reported that the combination of an accurate two-body ab initio potential with an empirically determined multibody contribution enables the prediction of the phase coexistence properties, the heats of vaporization, and the pair distribution functions of mercury with reasonable accuracy. In this work we present molecular dynamics simulation results for the shear viscosity and self-diffusion coefficient of mercury along the vapor-liquid coexistence curve using our empirical effective potential. The comparison with experiment and calculations based on a modified Enskog theory shows that our multibody contribution yields reliable predictions of the self-diffusion coefficient at all densities. Good results are also obtained for the shear viscosity of mercury at low to moderate densities. Increasing deviations between the simulation and experimental viscosity data at high densities suggest that not only a temperature-dependent but also a density-dependent multibody contribution is necessary to account for the effect of intermolecular interactions in liquid metals. An analysis of our simulation data near the critical point yields a critical exponent of beta = 0.39, which is identical to the value obtained from the analysis of the experimental saturation densities.     

Phase equilibria and interfacial tension of fluids confined in narrow pores

Correlation between phase behaviors of a Lennard-Jones fluid in and outside a pore is examined over wide thermodynamic conditions by grand canonical Monte Carlo simulations. A pressure tensor component of the confined fluid, a variable controllable in simulation but usually uncontrollable in experiment, is related with the pressure of a bulk homogeneous system in equilibrium with the confined system. Effects of the pore dimensionality, size, and attractive potential on the correlations between thermodynamic properties of the confined and bulk systems are clarified. A fluid-wall interfacial tension defined as an excess grand potential is evaluated as a function of the pore size. It is found that the tension decreases linearly with the inverse of the pore diameter or width.