No Arabic abstract
We introduce an accurate and efficient method for characterizing surface wetting and interfacial properties, such as the contact angle made by a liquid droplet on a solid surface, and the vapor-liquid surface tension of a fluid. The method makes use of molecular simulations in conjunction with the indirect umbrella sampling technique to systematically wet the surface and estimate the corresponding free energy. To illustrate the method, we study the wetting of a family of Lennard-Jones surfaces by water. We estimate contact angles for surfaces with a wide range of attractions for water by using our method and also by using droplet shapes. Notably, as surface-water attractions are increased, our method is able to capture the transition from partial to complete wetting. Finally, the method is straightforward to implement and computationally efficient, providing accurate contact angle estimates in roughly 5 nanoseconds of simulation time.
Within mean-field theory we study wetting of elastic substrates. Our analysis is based on a grand canonical free energy functional of the fluid number density and of the substrate displacement field. The substrate is described in terms of the linear theory of elasticity, parametrized by two Lame coefficients. The fluid contribution is of the van der Waals type. Two potentials characterize the interparticle interactions in the system. The long-ranged attraction between the fluid particles is described by a potential $w(r)$, and $v(r)$ characterizes the substrate-fluid interaction. By integrating out the elastic degrees of freedom we obtain an effective theory for the fluid number density alone. Its structure is similar to the one for wetting of an inert substrate. However, the potential $w(r)$ is replaced by an effective potential which, in addition to $w(r)$, contains a term bilinear in $v(r)$. We discuss the corresponding wetting transitions in terms of an effective interface potential $omega(ell)$, where $ell$ denotes the thickness of the wetting layer. We show that in the case of algebraically decaying interactions the elasticity of the substrate may suppress critical wetting transitions, and may even turn them first order.
A mean-field theory is presented which describes the basic observations of recent experiments revealing rich wetting behaviour of n-alkane/methanol mixtures at the liquid-vapour interface. The theory, qualitative and in part heuristic, is based on a microscopic lattice-gas model from which a Cahn-Landau approach is distilled. Besides the physics associated with the short-range components of the intermolecular interactions, effects of the long-range tails of the net van der Waals forces between interfaces are also taken into account. Including weak long-range forces which favour wetting in the theory does not visibly alter the critical wetting transition for the nonane/methanol mixture, in contrast with the generic expectation of first-order wetting for such systems, but in good agreement with experiment. For decane/methanol weak long-range forces bring the transition very close to the prewetting critical point, leading to an adsorption behaviour closely reminiscent of short-range tricritical wetting, observed experimentally for alkane chain length between 9.6 and 10. Finally, for undecane/methanol the transition is clearly of first order. First-order wetting is also seen in the experiment.
Water density fluctuations are an important statistical mechanical observable that is related to many-body correlations, as well as hydrophobic hydration and interactions. Local water density fluctuations at a solid-water surface have also been proposed as a measure of its hydrophobicity. These fluctuations can be quantified by calculating the probability, $P_v(N)$, of observing $N$ waters in a probe volume of interest $v$. When $v$ is large, calculating $P_v(N)$ using molecular dynamics simulations is challenging, as the probability of observing very few waters is exponentially small, and the standard procedure for overcoming this problem (umbrella sampling in $N$) leads to undesirable impulsive forces. Patel et al. [J. Phys. Chem. B, 114, 1632 (2010)] have recently developed an indirect umbrella sampling (INDUS) method, that samples a coarse-grained particle number to obtain $P_v(N)$ in cuboidal volumes. Here, we present and demonstrate an extension of that approach to other basic shapes, like spheres and cylinders, as well as to collections of such volumes. We further describe the implementation of INDUS in the NPT ensemble and calculate $P_v(N)$ distributions over a broad range of pressures. Our method may be of particular interest in characterizing the hydrophobicity of interfaces of proteins, nanotubes and related systems.
We report on the effect of intermolecular forces on the fluctuations of supported liquid films. Using an optically-induced thermal gradient, we form nanometer-thin films of wetting liquids on glass substrates, where van der Waals forces are balanced by thermocapillary forces. We show that the fluctuation dynamics of the film interface is strongly modified by intermolecular forces at lower frequencies. Data spanning three frequency decades are in excellent agreement with theoretical predictions accounting for van der Waals forces. Our results emphasize the relevance of intermolecular forces on thermal fluctuations when fluids are confined at the nanoscale.
The computational study of conformational transitions in RNA and proteins with atomistic molecular dynamics often requires suitable enhanced sampling techniques. We here introduce a novel method where concurrent metadynamics are integrated in a Hamiltonian replica-exchange scheme. The ladder of replicas is built with different strength of the bias potential exploiting the tunability of well-tempered metadynamics. Using this method, free-energy barriers of individual collective variables are significantly reduced compared with simple force-field scaling. The introduced methodology is flexible and allows adaptive bias potentials to be self-consistently constructed for a large number of simple collective variables, such as distances and dihedral angles. The method is tested on alanine dipeptide and applied to the difficult problem of conformational sampling in a tetranucleotide.