Monte Carlo simulation strategies to compute interfacial and bulk properties of binary fluid mixtures

Coupled Monte Carlo and Molecular Dynamics Simulations on Interfacial Properties of Antifouling Polymer Membranes

We use molecular simulation to study the wetting behavior of antifouling polymer-tethered membranes. We obtain the interfacial properties (e.g., contact angle) of water at various temperatures for five polymer membranes, including a base polysulfone (PSF) membrane and four other PSF membranes grafted with antifouling polymers (two poly(ethylene glycol) (PEG) tethers and two zwitterionic tethers). We implement a coupled Monte Carlo (MC)/molecular dynamics (MD) approach to determine the interface potentials of water on the membrane surfaces in an efficient manner. Within this method, short MC and MD simulations are performed in cycles to collect the surface excess free energy of a thin water film on polymer membrane surfaces. Simulation results show that the grafting of zwitterionic tethers provides a more significant enhancement in the hydrophilicity of the PSF membrane than that of the PEG tethers. Water completely wets the surface of zwitterionic polymer membranes.

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.

Effect of Carboxylic Acid on the Wetting Properties of a Model Water-Octane-Silica System

Monte Carlo simulations are employed to determine the effects of acetic acid on the wetting properties of a model water-octane-silica system. We first compute the bulk liquid-vapor saturation properties of pure acetic acid and subsequently explore the bulk liquid-liquid saturation properties of the ternary water-octane-acid system. We introduce an expanded ensemble approach to compute the coexistence properties of the ternary system. An interface potential approach is then used to capture the evolution of the wetting properties of the water-octane-silica system upon the addition of acetic acid. We track the change in the octane-water interfacial tension and the contact angle of water droplet on a silica substrate in a mother octane fluid over a range of acetic acid activities. The structure of the fluid, including the partitioning of acetic acid within the interfacial system, is also considered at several state points. We observe that acetic acid has a strong tendency to adsorb at the octane-water interface, resulting in a reduction in the octane-water interfacial tension. The response of the contact angle is more sensitive to the temperature and the hydrophilicity of the silica substrate.

Structure and Interfacial Tension of a Hard-Rod Fluid in Planar Confinement

The structural properties and interfacial tension of a fluid of rodlike hard-spherocylinder particles in contact with hard structureless flat walls are studied by means of Monte Carlo simulation. The calculated surface tension between the rod fluid and the substrate is characterized by a nonmonotonic trend as a function of the bulk concentration (density) over the range of isotropic bulk concentrations. As suggested by earlier theoretical studies, a surface-ordering scenario is confirmed by our simulations: the local orientational order close to the wall changes from uniaxial to biaxial nematic when the bulk concentration reaches about 85% of the value at the onset of the isotropic-nematic phase transition. The surface ordering coincides with a wetting transition whereby the hard wall is wetted by a nematic film. Accurate values of the fluid-solid surface tension, the adsorption, and the average particle-wall contact distance are reported (over a broad range of densities into the dense nematic region for the first time), which can serve as a useful benchmark for future theoretical and experimental studies on confined rod fluids. The simulation data are supplemented with predictions from second-virial density functional theory, which are in good qualitative agreement with the simulation results.