Wetting behavior of associating binary mixtures at attractive walls: a lattice Monte Carlo study

The Lattice Model of Particles with Orientation-Dependent Interactions at Solid Surfaces: Wetting Scenarios

Wetting phenomena in a lattice model of particles having two chemically different halves (A and B) and being in contact with solid substrates have been studied with Monte Carlo methods. The energy of the interaction between a pair of neighboring particles has been assumed to depend on the degree to which the AA, AB and BB regions face each other. In this work, we have assumed that uAA=−1.0 and considered three series of systems with uAB=uBB, uAB=0 and uBB=0. The phase behavior of bulk systems has been determined. In particular, it has been shown that at sufficiently low temperatures the bulk systems order into the superantiferromagnetic (SAF) phase, or into the antiferromagnetic (AF) phase, depending on the magnitudes of AA, AB and BB interaction energies, uAA, uAB and uBB. The SAF structure occurs whenever ϵ=uAA+uBB−2uAB is lower than zero and the AF structure is stable when ϵ is greater than zero. The wetting behavior has been demonstrated to depend strongly on the structure of the bulk condensed phase, the interactions between fluid particles and the strength of the surface potential. In all series, we have found the dewetting transition, resulting from the limited stability of different ordered structures of surface phases. However, in the systems that exhibit the gas–liquid transition in the bulk, the reentrant wetting transition has been observed at sufficiently high temperatures. The mechanism of dewetting and reentrant wetting transitions has been determined. Moreover, we have also demonstrated, how the dewetting transition in the series with uAB=0 is affected by the wall selectivity, i.e., when the interaction between the parts A and B of fluid particles and the solid is different.

Molecular Dynamics Study of Wetting and Adsorption of Binary Mixtures of the Lennard-Jones Truncated and Shifted Fluid on a Planar Wall

The wetting of surfaces is strongly influenced by adsorbate layers. Therefore, in this work, sessile drops and their interaction with adsorbate layers on surfaces were investigated by molecular dynamics simulations. Binary fluid model mixtures were considered. The two components of the fluid mixture have the same pure component parameters, but one component has a stronger and the other a weaker affinity to the surface. Furthermore, the unlike interactions between both components were varied. All interactions were described by the Lennard-Jones truncated and shifted potential with a cutoff radius of 2.5σ. The simulations were carried out at constant temperature for mixtures of different compositions. The parameters were varied systematically and chosen such that cases with partial wetting as well as cases with total wetting were obtained and the relation between the varied molecular parameters and the phenomenological behavior was elucidated. Data on the contact angle as well as on the mole fraction and thickness of the adsorbate layer were obtained, accompanied by information on liquid and gaseous bulk phases and the corresponding phase equilibrium. Also, the influence of the adsorbate layer on the wetting was studied: for a sufficiently thick adsorbate layer, the wall's influence on the wetting vanishes, which is then only determined by the adsorbate layer.

Highly non-additive symmetric mixtures at a wall

This paper discusses the results of the grand canonical ensemble Monte Carlo simulation of the wetting behavior of non-additive symmetric mixtures at non-selective walls. We have focused on the mixtures that exhibit closed immiscibility loops in the bulk, and the results obtained have demonstrated that such systems show a rather complex wetting behavior. In particular, such mixtures may exhibit a complete wetting at temperatures below the bulk demixing point, and an incomplete wetting at higher temperatures. Such a situation occurs when the adsorbed film remains mixed up to the bulk coexistence. However, close to the bulk tricritical point, being the onset of the continuous demixing transition (λ-line), a second wetting transition takes place. On the other hand, in the systems in which the adsorbed films undergo the demixing transition, a complete wetting may occur below as well as at and above the bulk demixing transition temperature. It has also been demonstrated that the wetting behavior depends strongly on the strength of the surface-fluid interaction.

Histogram analysis as a method for determining the line tension of a three-phase contact region by Monte Carlo simulations

A method is proposed for determining the line tension, which is the main physical characteristic of a three-phase contact region, by Monte Carlo (MC) simulations. The key idea of the proposed method is that if a three-phase equilibrium involves a three-phase contact region, the probability distribution of states of a system as a function of two order parameters depends not only on the surface tension, but also on the line tension. This probability distribution can be obtained as a normalized histogram by appropriate MC simulations, so one can use the combination of histogram analysis and finite-size scaling to study the properties of a three phase contact region. Every histogram and results extracted therefrom will depend on the size of the simulated system. Carrying out MC simulations for a series of system sizes and extrapolating the results, obtained from the corresponding series of histograms, to infinite size, one can determine the line tension of the three phase contact region and the interfacial tensions of all three interfaces (and hence the contact angles) in an infinite system. To illustrate the proposed method, it is applied to the three-dimensional ternary fluid mixture, in which molecular pairs of like species do not interact whereas those of unlike species interact as hard spheres. The simulated results are in agreement with expectations.