No abstract

No abstract

In the last few years, much attention has been devoted to the control of the wettability properties of surfaces modified with functional groups. Molecular dynamics (MD) simulation is one of the powerful tools for microscopic analysis providing visual images and mean geometrical shapes of the contact line, e.g., of nanoscale droplets on solid surfaces, while profound understanding of wetting demands quantitative evaluation of the solid-liquid (SL) interfacial tension. In the present work, we examined the wetting of water on neutral and regular hydroxylated silica surfaces with five different area densities of OH groups ρ OH A , ranging from a nonhydroxylated surface to a fully hydroxylated one through two theoretical methods: thermodynamic integration (TI) and MD simulations of quasi-two-dimensional equilibrium droplets. For the former, the work of adhesion needed to quasi-statically strip the water film off the solid surface was computed by the phantom wall TI scheme to evaluate the SL interfacial free energy, whereas for the latter, the apparent contact angle θapp was calculated from the droplet density distribution. The theoretical contact angle θ YD and the apparent one θapp, both indicating the enhancement of wettability by an increase in ρ OH A , presented good quantitative agreement, especially for non-hydroxylated and highly hydroxylated surfaces. On partially hydroxylated surfaces, in which θ YD and θapp slightly deviated, the Brownian motion of the droplet was suppressed, possibly due to the pinning of the contact line around the hydroxyl groups. Relations between work of adhesion, interfacial energy, and entropy loss were also analyzed, and their influence on the wettability was discussed.

We have given theoretical expressions for the forces exerted on a so-called Wilhelmy plate, which we modeled as a quasi-2D flat and smooth solid plate immersed in a liquid pool of a simple liquid. All forces given by the theory, the local forces on the top, the contact line, and the bottom of the plate as well as the total force, showed an excellent agreement with the MD simulation results. The force expressions were derived by a purely mechanical approach, which is exact and ensures the force balance on the control volumes arbitrarily set in the system, and are valid as long as the solid–liquid (SL) and solid–vapor (SV) interactions can be described by mean-fields. In addition, we revealed that the local forces around the bottom and top of the solid plate can be related to the SL and SV interfacial tensions γSL and γSV, and this was verified through the comparison with the SL and SV works of adhesion obtained by the thermodynamic integration (TI). From these results, it has been confirmed that γSL and γSV as well as the liquid–vapor interfacial tension γLV can be extracted from a single equilibrium MD simulation without the computationally demanding calculation of the local stress distributions and the TI.

We describe a method to locate stationary points in the free-energy hypersurface of complex molecular systems using high-level correlated ab initio potentials. In this work, we assume a combined QM/MM description of the system although generalization to full ab initio potentials or other theoretical schemes is straightforward. The free-energy gradient (FEG) is obtained as the mean force acting on relevant nuclei using a dual level strategy. First, a statistical simulation is carried out using an appropriate, low-level quantum mechanical force-field. Free-energy perturbation (FEP) theory is then used to obtain the free-energy derivatives for the target, high-level quantum mechanical force-field. We show that this composite FEG-FEP approach is able to reproduce the results of a standard free-energy minimization procedure with high accuracy, while simultaneously allowing for a drastic reduction of both computational and wall-clock time. The method has been applied to study the structure of the water molecule in liquid water at the QCISD/aug-cc-pVTZ level of theory, using the sampling from QM/MM molecular dynamics simulations at the B3LYP/6-311+G(d,p) level. The obtained values for the geometrical parameters and for the dipole moment of the water molecule are within the experimental error, and they also display an excellent agreement when compared to other theoretical estimations. The developed methodology represents therefore an important step toward the accurate determination of the mechanism, kinetics, and thermodynamic properties of processes in solution, in enzymes, and in other disordered chemical systems using state-of-the-art ab initio potentials.