Electrostatic and induction effects in the solubility of water in alkanes

Evaluation of Osmotic Virial Coefficients via Restricted Gibbs Ensemble Simulations, with Support from Gas-Phase Mixture Coefficients

We present a method for computing osmotic virial coefficients in explicit solvent via simulation in a restricted Gibbs ensemble. Two equivalent phases are simulated at once, each in a separate box at constant volume and temperature and each in equilibrium with a solvent reservoir. For osmotic coefficient BN, a total of N solutes are individually exchanged back and forth between the boxes, and the average distribution of solute numbers between the boxes provides the key information needed to compute BN. Separately, expressions are developed for BN as a series in solvent reservoir density ρ₁, with the coefficients of the series expressed in terms of the usual gas-phase mixture coefficients B_ij. Normally, the B_ij are defined for an infinite volume, but we suggest that the observed dependence of B_ij on system size L can be used to estimate L dependence of the BN, allowing them to be computed accurately at L → ∞ while simulating much smaller system sizes than otherwise possible. The methods for N = 2 and 3 are demonstrated for two-component mixtures of size-asymmetric additive hard spheres. The proposed methods are demonstrated to have greater precision than established techniques, for a given amount of computational effort. The ρ₁ series for BN when applied by itself is (for this noncondensing model) found to be the most efficient in computing accurate osmotic coefficients for the solvent densities considered here.

Miscibility at the immiscible liquid/liquid interface: A molecular dynamics study of thermodynamics and mechanism

Molecular dynamics simulations are used to study the dissolution of water into an adjacent, immiscible organic liquid phase. Equilibrium thermodynamic and structural properties are calculated during the transfer of water molecule(s) across the interface using umbrella sampling. The net free energy of transfer agrees reasonably well with experimental solubility values. We find that water molecules "prefer" to transfer into the adjacent phase one-at-a-time, without co-transfer of the hydration shell, as in the case of evaporation. To study the dynamics and mechanism of transfer of water to liquid nitrobenzene, we collected over 400 independent dissolution events. Analysis of these trajectories suggests that the transfer of water is facilitated by interfacial protrusions of the water phase into the organic phase, where one water molecule at the tip of the protrusion enters the organic phase by the breakup of a single hydrogen bond.