Studying vapor-liquid transition using a generalized ensemble.

Homogeneous vapor-liquid nucleation is studied using the generalized Replica Exchange Method (gREM). The generalized ensemble allows the study of unstable states that cannot directly be studied in the canonical ensemble. Along with replica exchange, this allows for efficient sampling of the multiple states in a single simulation. Statistical Temperature Weighted Histogram Analysis Method is used for postprocessing to get a continuous free energy curve from bulk vapor to bulk liquid. gREM allows the study of planar, cylindrical, and spherical interfaces in a single simulation. The excess Gibbs free energy for the formation of a spherical liquid droplet in vapor for a Lennard-Jones system is calculated from the free energy curve and compared against the umbrella sampling results. The nucleation free energy barrier obtained from gREM is then used to calculate the nucleation rate without relying on any classification scheme for separating the vapor and liquid.

Electrostatic and induction effects in the solubility of water in alkanes.

Experiments show that at 298 K and 1 atm pressure, the transfer free energy, µex, of water from its vapor to liquid normal alkanes CnH2n+2 (n=5…12) is negative. Earlier it was found that with the united-atom TraPPE model for alkanes and the SPC/E model for water, one had to artificially enhance the attractive alkane-water cross interaction to capture this behavior. Here we revisit the calculation of µex using the polarizable AMOEBA and the non-polarizable Charmm General (CGenFF) forcefields. We test both the AMOEBA03 and AMOEBA14 water models; the former has been validated with the AMOEBA alkane model while the latter is a revision of AMOEBA03 to better describe liquid water. We calculate µex using the test particle method. With CGenFF, µex is positive and the error relative to experiments is about 1.5 kBT. With AMOEBA, µex is negative and deviations relative to experiments are between 0.25 kBT (AMOEBA14) and 0.5 kBT (AMOEBA03). Quantum chemical calculations in a continuum solvent suggest that zero point effects may account for some of the deviation. Forcefield limitations notwithstanding, electrostatic and induction effects, commonly ignored in consideration of water-alkane interactions, appear to be decisive in the solubility of water in alkanes.

Competition between Intra- and Intermolecular Association of Chain Molecules with Water-like Solvent.

Fluid properties and phase behavior of systems such as glycol ethers, carboxylic acids, and proteins are affected by the competition between intra- and intermolecular hydrogen bonding. Here we study this competition by extending Wertheim's first-order thermodynamic perturbation theory to include intramolecular hydrogen bonding in chain molecules in the presence of an explicit water-like solvent. The theory derived here is found to be in good agreement with molecular simulation. It is shown that intramolecular association is most important for shorter chains at low temperature, low density, and high chain concentration. The theory is also extended into a density functional theory formalism to study the effect of intramolecular association on the structuring of the different segments of the molecules close to a hydrophobic surface. Intramolecular association is found to be enhanced close to the surface, with the total density of the system having the most effect on structuring close to the surface.

Isolating the non-polar contributions to the intermolecular potential for water-alkane interactions.

Intermolecular potential models for water and alkanes describe pure component properties fairly well, but fail to reproduce properties of water-alkane mixtures. Understanding interactions between water and non-polar molecules like alkanes is important not only for the hydrocarbon industry but has implications to biological processes as well. Although non-polar solutes in water have been widely studied, much less work has focused on water in non-polar solvents. In this study we calculate the solubility of water in different alkanes (methane to dodecane) at ambient conditions where the water content in alkanes is very low so that the non-polar water-alkane interactions determine solubility. Only the alkane-rich phase is simulated since the fugacity of water in the water rich phase is calculated from an accurate equation of state. Using the SPC/E model for water and TraPPE model for alkanes along with Lorentz-Berthelot mixing rules for the cross parameters produces a water solubility that is an order of magnitude lower than the experimental value. It is found that an effective water Lennard-Jones energy ε(W)/k = 220 K is required to match the experimental water solubility in TraPPE alkanes. This number is much higher than used in most simulation water models (SPC/E-ε(W)/k = 78.2 K). It is surprising that the interaction energy obtained here is also higher than the water-alkane interaction energy predicted by studies on solubility of alkanes in water. The reason for this high water-alkane interaction energy is not completely understood. Some factors that might contribute to the large interaction energy, such as polarizability of alkanes, octupole moment of methane, and clustering of water at low concentrations in alkanes, are examined. It is found that, though important, these factors do not completely explain the anomalously strong attraction between alkanes and water observed experimentally.

Hydrophobic and hydrophilic interactions in aqueous mixtures of alcohols at a hydrophobic surface.

Aqueous solutions of alcohols are interesting because of their anomalous behavior that is believed to be due to the molecular structuring of water and alcohol around each other in solution. The interfacial structuring and properties are significant for application in alcohol purification processes and biomolecular structure. Here we study aqueous mixtures of short alcohols (methanol, ethanol, 1-propanol, and 2-propanol) at a hydrophobic surface using interfacial statistical associating fluid theory which is a perturbation density functional theory. The addition of a small amount of alcohol decreases the interfacial tension of water drastically. This trend in interfacial tension can be explained by the structure of water and alcohol next to the surface. The hydrophobic group of an added alcohol preferentially goes to the surface preserving the structure of water in the bulk. For a given bulk alcohol concentration, water mixed with the different alcohols has different interfacial tensions with propanol having a lower interfacial tension than methanol and ethanol. 2-propanol is not as effective in decreasing the interfacial tension as 1-propanol because it partitions poorly to the surface due to its larger excluded volume. But for a given surface alcohol mole fraction, all the alcohol mixtures give similar values for interfacial tension. For separation of alcohol from water, methods that take advantage of the high surface mole fraction of alcohol have advantages compared to separation using the vapor in equilibrium with a water-alcohol liquid.

Wertheim's association theory applied to one site patchy colloids: beyond the single bonding condition.

We apply Wertheim's theory to develop an equation of state for one site patchy colloids where the patch can bond multiple times. We allow for the possibility of ring formation without the introduction of empirical parameters and show that for moderate patch coverage the infinite series of chain graphs is well represented by the first two terms. The theory is found to be in excellent agreement with new NVT and NPT Monte Carlo simulations. The approach described here can easily be converted to the form of a density functional theory to describe inhomogeneous patchy colloid systems.