No Arabic abstract
In our first paper, we showed how a non-local effective Hamiltionian for short-ranged wetting may be derived from an underlying Landau-Ginzburg-Wilson model. Here, we combine the Greens function method with standard perturbation theory to determine the general diagrammatic form of the binding potential functional beyond the double-parabola approximation for the Landau-Ginzburg-Wilson bulk potential. The main influence of cubic and quartic interactions is simply to alter the coefficients of the double parabola-like zig-zag diagrams and also to introduce curvature and tube-interaction corrections (also represented diagrammatically), which are of minor importance. Non-locality generates effective long-ranged many-body interfacial interactions due to the reflection of tube-like fluctuations from the wall. Alternative wall boundary conditions (with a surface field and enhancement) and the diagrammatic description of tricritical wetting are also discussed.
Coulomb interactions are present in a wide variety of all-atom force fields. Spherical truncations of these interactions permit fast simulations but are problematic due to their incorrect thermodynamics. Herein we demonstrate that simple analytical corrections for the thermodynamics of uniform truncated systems are possible. In particular results for the SPC/E water model treated with spherically-truncated Coulomb interactions suggested by local molecular field theory [Proc. Nat. Acad. Sci. USA 105, 19136 (2008)] are presented. We extend results developed by Chandler [J. Chem. Phys. 65, 2925 (1976)] so that we may treat the thermodynamics of mixtures of flexible charged and uncharged molecules simulated with spherical truncations. We show that the energy and pressure of spherically-truncated bulk SPC/E water are easily corrected using exact second-moment-like conditions on long-ranged structure. Furthermore, applying the pressure correction as an external pressure removes the density errors observed by other research groups in NPT simulations of spherically-truncated bulk species.
We show that spherical truncations of the 1/r interactions in models for water and acetonitrile yield very accurate results in bulk simulations for all site-site pair correlation functions as well as dipole-dipole correlation functions. This good performance in bulk simulations contrasts with the generally poor results found with the use of such truncations in nonuniform molecular systems. We argue that Local Molecular Field (LMF) theory provides a general theoretical framework that gives the necessary corrections to simple truncations in most nonuniform environments and explains the accuracy of spherical truncations in uniform environments by showing that these corrections are very small. LMF theory is derived from the exact Yvon-Born-Green (YBG) hierarchy by making physically-motivated and well-founded approximations. New and technically interesting derivations of both the YBG hierarchy and LMF theory for a variety of site-site molecular models are presented in appendices. The main paper focuses on understanding the accuracy of these spherical truncations in uniform systems both phenomenologically and quantitatively using LMF theory.
Work fluctuations and work probability distributions are fundamentally different in systems with short- ranged versus long-ranged correlations. Specifically, in systems with long-ranged correlations the work distribution is extraordinarily broad compared to systems with shortranged correlations. This difference profoundly affects the possible applicability of fluctuation theorems like the Jarzynski fluctuation theorem. The Heisenberg ferromagnet , well below its Curie temperature, is a system with long-ranged correlations in very low magnetic fields due to the presence of Goldstone modes. As the magnetic field is increased the correlations gradually become short-ranged. Hence, such a ferromagnet is an ideal system for elucidating the changes of the work probability distribution as one goes from a domain with long-ranged correlations to a domain with short-ranged correlations by tuning the magnetic field. A quantitative analysis of this crossover behaviour of the work probability distribution and the associated fluctuations is presented.
Spherical truncations of Coulomb interactions in standard models for water permit efficient molecular simulations and can give remarkably accurate results for the structure of the uniform liquid. However truncations are known to produce significant errors in nonuniform systems, particularly for electrostatic properties. Local molecular field (LMF) theory corrects such truncations by use of an effective or restructured electrostatic potential that accounts for effects of the remaining long-ranged interactions through a density-weighted mean field average and satisfies a modified Poissons equation defined with a Gaussian-smoothed charge density. We apply LMF theory to three simple molecular systems that exhibit different aspects of the failure of a naive application of spherical truncations -- water confined between hydrophobic walls, water confined between atomically-corrugated hydrophilic walls, and water confined between hydrophobic walls with an applied electric field. Spherical truncations of 1/r fail spectacularly for the final system in particular, and LMF theory corrects the failings for all three. Further, LMF theory provides a more intuitive way to understand the balance between local hydrogen bonding and longer-ranged electrostatics in molecular simulations involving water.
The excess adsorption $Gamma $ in two-dimensional Ising strips $(infty times L)$ subject to identical boundary fields, at both one-dimensional surfaces decaying in the orthogonal direction $j$ as $-h_1j^{-p}$, is studied for various values of $p$ and along various thermodynamic paths below the critical point by means of the density-matrix renormalization-group method. The crossover behavior between the complete wetting and critical adsorption regimes, occurring in semi-infinite systems, are strongly influenced by confinement effects. Along isotherms $T=const$ the asymptotic power law dependences on the external bulk field, which characterize these two regimes, are undercut by capillary condensation. Along the pseudo first-order phase coexistence line of the strips, which varies with temperature, we find a broad crossover regime where both the thickness of the wetting film and $Gamma$ increase as function of the reduced temperature $tau$ but do not follow any power law. Above the wetting temperature the order parameter profiles are not slab-like but exhibit wide interfacial variations and pronounced tails. Inter alia, our explicit calculations demonstrate that, contrary to opposite claims by Kroll and Lipowsky [Phys. Rev. B {bf 28}, 5273 (1983)], for $p=2$ critical wetting transitions do exist and we determine the corresponding wetting phase diagram in the $(h_1,T)$ plane.