Percolation transition of hydration water at hydrophilic surfaces

An analysis of water clustering is used to study the quasi-2D percolation transition of water adsorbed at planar hydrophilic surfaces. Above the critical temperature of the layering transition (quasi-2D liquid-vapor phase transition of adsorbed molecules) a percolation transition occurs at some threshold surface coverage, which increases with increasing temperature. The location of the percolation line is consistent with the existence of a percolation transition at the critical point. The percolation threshold at a planar surface is weakly sensitive to the size of the system when its lateral dimension increases from 80 to 150 A. The size distribution of the largest water cluster shows a specific two-peaks structure in a wide range of surface coverage : the lower- and higher-size peaks represent contributions from non-spanning and spanning clusters, respectively. The ratio of the average sizes of spanning and non-spanning largest clusters is about 1.8 for all studied planes. The two-peak structure becomes more pronounced with decreasing size of the planar surface and strongly enhances at spherical surfaces.Comment: 17 pages, 11 figure

On the surface critical behaviour in Ising strips: density-matrix renormalization-group study

Using the density-matrix renormalization-group method we study the surface critical behaviour of the magnetization in Ising strips in the subcritical region. Our results support the prediction that the surface magnetization in the two phases along the pseudo-coexistence curve also behaves as for the ordinary transition below the wetting temperature for the finite value of the surface field.Comment: 15 pages, 9 figure

Surface critical behavior of fluids: Lennard-Jones fluid near weakly attractive substrate

The phase behavior of fluids near weakly attractive substrates is studied by computer simulations of the coexistence curve of a Lennard-Jones (LJ) fluid confined in a slitlike pore. The temperature dependence of the density profiles of the LJ fluid is found to be very similar to the behavior of water near hydrophobic surfaces (Brovchenko et al. J.Phys.: Cond.Matt. v.16, 2004). A universal critical behavior of the local order parameter, defined as the difference between the local densities of the coexisting liquid and vapor phases at some distance z from the pore walls, Deltarho(z) = (rho_l(z) - rho_v(z))/2, is observed in a wide temperature range and found to be consistent with the surface critical behavior of the Ising model. Near the surface the dependence of the order parameter on the reduced temperature tau = (T_c - T)/T_c obeys a scaling law ~ tau^(beta_1) with a critical exponent beta_1 of about 0.8, corresponding to the ordinary surface transition. A crossover from bulk-like to surface-like critical behavior with increasing temperature occurs, when the correlation length is about half the distance to the surface. Relations between the ordinary and normal transitions in Ising systems and the surface critical behavior of fluids are discussed.Comment: 14 pages, 19 figures, submitted to PR

The marine kd and water/sediment interaction problem

The behavior of marine distribution coefficients is analyzed with the help of numerical experiments and ana-lytical solutions of equations describing kinetic models for uptake/release of radionuclides. The difficulties in measuring true kd in a marine environment perturbed by an external radionuclide source are highlighted. Differences between suspended matter and bed sediment kd are analyzed. The performances of different kinetic models (1-step/2step; single-layer/multi-layer) are studied in model/model and model/experiment compar-isons. Implications for the use of models to assess radioactive contamination after an emergency are given; as well as recommendations when kd data are compiled in order to create a useful databaseInternational Atomic Energy Agency (IAEA) KIOST PE9961

Electric field inside a "Rossky cavity" in uniformly polarized water

Electric field produced inside a solute by a uniformly polarized liquid is strongly affected by dipolar polarization of the liquid at the interface. We show, by numerical simulations, that the electric "cavity" field inside a hydrated non-polar solute does not follow the predictions of standard Maxwell's electrostatics of dielectrics. Instead, the field inside the solute tends, with increasing solute size, to the limit predicted by the Lorentz virtual cavity. The standard paradigm fails because of its reliance on the surface charge density at the dielectric interface determined by the boundary conditions of the Maxwell dielectric. The interface of a polar liquid instead carries a preferential in-plane orientation of the surface dipoles thus producing virtually no surface charge. The resulting boundary conditions for electrostatic problems differ from the traditional recipes, affecting the microscopic and macroscopic fields based on them. We show that relatively small differences in cavity fields propagate into significant differences in the dielectric constant of an ideal mixture. The slope of the dielectric increment of the mixture versus the solute concentration depends strongly on which polarization scenario at the interface is realized. A much steeper slope found in the case of Lorentz polarization also implies a higher free energy penalty for polarizing such mixtures.Comment: 9 pages, 8 figure

Double dynamical regime of confined water

The Van Hove self correlation function of water confined in a silica pore is calculated from Molecular Dynamics trajectories upon supercooling. At long time in the $\alpha$ relaxation region we found that the behaviour of the real space time dependent correlators can be decomposed in a very slow, almost frozen, dynamics due to the bound water close to the substrate and a faster dynamics of the free water which resides far from the confining surface. For free water we confirm the evidences of an approach to a crossover mode coupling transition, previously found in Q space. In the short time region we found that the two dynamical regimes are overimposed and cannot be distinguished. This shows that the interplay between the slower and the faster dynamics emerges in going from early times to the $\alpha$ relaxation region, where a layer analysis of the dynamical properties can be performed.Comment: 6 pages with 9 figures. RevTeX. Accepted for pulbication in J. Phys. Cond. Mat

Simulation of fluid-solid coexistence in finite volumes: A method to study the properties of wall-attached crystalline nuclei

The Asakura-Oosawa model for colloid-polymer mixtures is studied by Monte Carlo simulations at densities inside the two-phase coexistence region of fluid and solid. Choosing a geometry where the system is confined between two flat walls, and a wall-colloid potential that leads to incomplete wetting of the crystal at the wall, conditions can be created where a single nanoscopic wall-attached crystalline cluster coexists with fluid in the remainder of the simulation box. Following related ideas that have been useful to study heterogeneous nucleation of liquid droplets at the vapor-liquid coexistence, we estimate the contact angles from observations of the crystalline clusters in thermal equilibrium. We find fair agreement with a prediction based on Young's equation, using estimates of interface and wall tension from the study of flat surfaces. It is shown that the pressure versus density curve of the finite system exhibits a loop, but the pressure maximum signifies the "droplet evaporation-condensation" transition and thus has nothing in common with a van der Waals-like loop. Preparing systems where the packing fraction is deep inside the two-phase coexistence region, the system spontaneously forms a "slab state", with two wall-attached crystalline domains separated by (flat) interfaces from liquid in full equilibrium with the crystal in between; analysis of such states allows a precise estimation of the bulk equilibrium properties at phase coexistence

Does Young's equation hold on the nanoscale? A Monte Carlo test for the binary Lennard-Jones fluid

When a phase-separated binary ($A+B$) mixture is exposed to a wall, that preferentially attracts one of the components, interfaces between A-rich and B-rich domains in general meet the wall making a contact angle $\theta$. Young's equation describes this angle in terms of a balance between the $A-B$ interfacial tension $\gamma_{AB}$ and the surface tensions $\gamma_{wA}$, $\gamma_{wB}$ between, respectively, the $A$- and $B$-rich phases and the wall, $\gamma _{AB} \cos \theta =\gamma_{wA}-\gamma_{wB}$. By Monte Carlo simulations of bridges, formed by one of the components in a binary Lennard-Jones liquid, connecting the two walls of a nanoscopic slit pore, $\theta$ is estimated from the inclination of the interfaces, as a function of the wall-fluid interaction strength. The information on the surface tensions $\gamma_{wA}$, $\gamma_{wB}$ are obtained independently from a new thermodynamic integration method, while $\gamma_{AB}$ is found from the finite-size scaling analysis of the concentration distribution function. We show that Young's equation describes the contact angles of the actual nanoscale interfaces for this model rather accurately and location of the (first order) wetting transition is estimated.Comment: 6 pages, 6 figure