Hydration of an apolar solute in a two-dimensional waterlike lattice fluid

In a previous work, we investigated a two-dimensional lattice-fluid model, displaying some waterlike thermodynamic anomalies. The model, defined on a triangular lattice, is now extended to aqueous solutions with apolar species. Water molecules are of the "Mercedes Benz" type, i.e., they possess a D3 (equilateral triangle) symmetry, with three equivalent bonding arms. Bond formation depends both on orientation and local density. The insertion of inert molecules displays typical signatures of hydrophobic hydration: large positive transfer free energy, large negative transfer entropy (at low temperature), strong temperature dependence of the transfer enthalpy and entropy, i.e., large (positive) transfer heat capacity. Model properties are derived by a generalized first order approximation on a triangle cluster.Comment: 9 pages, 5 figures, 1 table; submitted to Phys. Rev.

Effect of hydrophobic solutes on the liquid-liquid critical point

Jagla ramp particles, interacting through a ramp potential with two characteristic length scales, are known to show in their bulk phase thermodynamic and dynamic anomalies, similar to what is found in water. Jagla particles also exhibit a line of phase transitions separating a low density liquid phase and a high density liquid phase, terminating in a liquid-liquid critical point in a region of the phase diagram that can be studied by simulations. Employing molecular dynamics computer simulations, we study the thermodynamics and the dynamics of solutions of hard spheres (HS) in a solvent formed by Jagla ramp particles. We consider the cases of HS mole fraction x = 0.10, 0.15 and 0.20, and also the case x = 0.50 (a 1:1 mixture of HS and Jagla particles). We find a liquid-liquid critical point, up to the highest HS mole fraction; its position shifts to higher pressures and lower temperatures upon increasing x. We also find that the diffusion coefficient anomalies appear to be preserved for all the mole fractions studied.Comment: 8 pages, 7 figures, 1 table. In press (Phys. Rev. E

A lattice model of hydrophobic interactions

Hydrogen bonding is modeled in terms of virtual exchange of protons between water molecules. A simple lattice model is analyzed, using ideas and techniques from the theory of correlated electrons in metals. Reasonable parameters reproduce observed magnitudes and temperature dependence of the hydrophobic interaction between substitutional impurities and water within this lattice.Comment: 7 pages, 3 figures. To appear in Europhysics Letter

Solvent-free coarse-grained lipid model for large-scale simulations

A coarse-grained molecular model, which consists of a spherical particle and an orientation vector, is proposed to simulate lipid membrane on a large length scale. The solvent is implicitly represented by an effective attractive interaction between particles. A bilayer structure is formed by orientation-dependent (tilt and bending) potentials. In this model, the membrane properties (bending rigidity, line tension of membrane edge, area compression modulus, lateral diffusion coefficient, and flip-flop rate) can be varied over broad ranges. The stability of the bilayer membrane is investigated via droplet-vesicle transition. The rupture of the bilayer and worm-like micelle formation can be induced by an increase in the spontaneous curvature of the monolayer membrane.Comment: 13 pages, 19 figure

Micellar Aggregates of Gemini Surfactants: Monte Carlo Simulation of a Microscopic Model

We propose a "microscopic" model of gemini surfactants in aqueous solution. Carrying out extensive Monte Carlo simulations, we study the variation of the critical micellar concentration (CMC) of these model gemini surfactants with the variation of the (a) length of the spacer connecting the two hydrophilic heads, (b) length of the hydrophobic tail and (c) the bending rigidity of the hydrocarbon chains forming the spacer and the tail; some of the trends of variation are counter-intuitive but are in excellent agreement with the available experimental results. Our simulations also elucidate the dependence of the shapes of the micellar aggregates and the magnitude of the CMC on the geometrical shape and size of the surfactant molecules and the electrical charge on the hydrophilic heads

Adding Salt to an Aqueous Solution of t-Butanol: Is Hydrophobic Association Enhanced or Reduced?

Recent neutron scattering experiments on aqueous salt solutions of amphiphilic t-butanol by Bowron and Finney [Phys. Rev. Lett. {\bf 89}, 215508 (2002); J. Chem. Phys. {\bf 118}, 8357 (2003)] suggest the formation of t-butanol pairs, bridged by a chloride ion via ${O}-{H}...{Cl}^-$ hydrogen-bonds, and leading to a reduced number of intermolecular hydrophobic butanol-butanol contacts. Here we present a joint experimental/theoretical study on the same system, using a combination of molecular dynamics simulations and nuclear magnetic relaxation measurements. Both theory and experiment clearly support the more intuitive scenario of an enhanced number of hydrophobic contacts in the presence of the salt, as it would be expected for purely hydrophobic solutes [J. Phys. Chem. B {\bf 107}, 612 (2003)]. Although our conclusions arrive at a structurally completely distinct scenario, the molecular dynamics simulation results are within the experimental errorbars of the Bowron and Finney work.Comment: 15 pages twocolumn revtex, 11 figure

Simulations of a lattice model of two-headed linear amphiphiles: influence of amphiphile asymmetry

Using a 2D lattice model, we conduct Monte Carlo simulations of micellar aggregation of linear-chain amphiphiles having two solvophilic head groups. In the context of this simple model, we quantify how the amphiphile architecture influences the critical micelle concentration (CMC), with a particular focus on the role of the asymmetry of the amphiphile structure. Accordingly, we study all possible arrangements of the head groups along amphiphile chains of fixed length $N=12$ and 16 molecular units. This set of idealized amphiphile architectures approximates many cases of symmetric and asymmetric gemini surfactants, double-headed surfactants and boloform surfactants. Consistent with earlier results, we find that the number of spacer units $s$ separating the heads has a significant influence on the CMC, with the CMC increasing with $s$ for $s<N/2$. In comparison, the influence of the asymmetry of the chain architecture on the CMC is much weaker, as is also found experimentally.Comment: 30 pages, 17 fgure

Origin of Native Driving Force in Protein Folding

We derive an expression with four adjustable parameters that reproduces well the 20x20 Miyazawa-Jernigan potential matrix extracted from known protein structures. The numerical values of the parameters can be approximately computed from the surface tension of water, water-screened dipole interactions between residues and water and among residues, and average exposures of residues in folded proteins.Comment: LaTeX file, Postscript file; 4 pages, 1 figure (mij.eps), 2 table

Induced Crystallization of Polyelectrolyte-Surfactant Complexes at the Gas-Water Interface

Synchrotron-X-ray and surface tension studies of a strong polyelectrolyte (PE) in the semi-dilute regime (~ 0.1M monomer-charges) with varying surfactant concentrations show that minute surfactant concentrations induce the formation of a PE-surfactant complex at the gas/solution interface. X-ray reflectivity and grazing angle X-ray diffraction (GIXD) provide detailed information of the top most layer, where it is found that the surfactant forms a two-dimensional liquid-like monolayer, with a noticeable disruption of the structure of water at the interface. With the addition of salt (NaCl) columnar-crystals with distorted-hexagonal symmetry are formed.Comment: 4 pages, 5 eps figure

Solvent-induced micelle formation in a hydrophobic interaction model

We investigate the aggregation of amphiphilic molecules by adapting the two-state Muller-Lee-Graziano model for water, in which a solvent-induced hydrophobic interaction is included implicitly. We study the formation of various types of micelle as a function of the distribution of hydrophobic regions at the molecular surface. Successive substitution of non-polar surfaces by polar ones demonstrates the influence of hydrophobicity on the upper and lower critical solution temperatures. Aggregates of lipid molecules, described by a refinement of the model in which a hydrophobic tail of variable length interacts with different numbers of water molecules, are stabilized as the length of the tail increases. We demonstrate that the essential features of micelle formation are primarily solvent-induced, and are explained within a model which focuses only on the alteration of water structure in the vicinity of the hydrophobic surface regions of amphiphiles in solution.Comment: 11 pages, 10 figures; some rearrangement of introduction and discussion sections, streamlining of formalism and general compression; to appear in Phys. Rev.