Frédéric Leroy

We investigated the solid-liquid work of adhesion of water on a model silica surface by molecular dynamics simulations, where a methodology previously developed to determine the work of adhesion through thermodynamic integration was extended to a system with long-range electrostatic interactions between solid and liquid. In agreement with previous studies, the work of adhesion increased when the magnitude of the surface polarity was increased. On the other hand, we found that when comparing two systems with and without solid-liquid electrostatic interactions, which were set to have approximately the same total solid-liquid interfacial energy, former had a significantly smaller work of adhesion and a broader distribution in the interfacial energies, which has not been previously reported in detail. This was explained by the entropy contribution to the adhesion free energy; i.e., the former with a broader energy distribution had a larger interfacial entropy than the latter. While the entropy contribution to the work of adhesion has already been known, as a work of adhesion itself is free energy, these results indicate that, contrary to common belief, wetting behavior such as the contact angle is not only governed by the interfacial energy but also significantly affected by the interfacial entropy. Finally, a new interpretation of interfacial entropy in the context of solid-liquid energy variance was offered, from which a fast way to qualitatively estimate the work of adhesion was also presented.

The high interfacial tension between two immiscible liquids can provide the necessary driving force for the self-assembly of nanoparticles at the interface. Particularly, the interface between water and oily liquids (hydrocarbon chains) has been exploited to prepare networks of highly interconnected graphene sheets of only a few layers thick that are well suited for industrial applications. Studying such complex systems through particle-based simulations could greatly enhance the understanding of the various driving forces in action and could possibly give more control over the self-assembly process. However, the interaction potentials used in particle-based simulations are typically derived by reproducing bulk properties and are therefore not suitable to describe systems dominated by interfaces. To address this issue, we introduce a methodology to derive solid-liquid interaction potentials that yield an accurate representation of the balance between interfacial interactions at atomistic and coarse-grained resolutions. Our approach is validated through its ability to lead to the adsorption of graphene nanoflakes at the interface between water and n-hexane. The development of accurate coarse-grained potentials that our approach enables makes it possible to perform large scale simulations to study the assembly of graphene nanoparticles at the interface between immiscible liquids. Our methodology is illustrated through a simulation of many graphene nanoflakes adsorbing at the interface.

Coarse-grained models have increasingly been used in large-scale particle-based simulations. However, due to their lack of degrees of freedom, it is a priori unlikely that they straightforwardly represent thermal properties with the same accuracy as their atomistic counterparts. We take a first step in addressing the impact of liquid coarse-graining on interfacial heat conduction by showing that an atomistic and a coarse-grained model of water may yield similar values of the Kapitza conductance on few-layer graphene with interactions ranging from hydrophobic to mildly hydrophilic. By design the water models employed yield similar liquid layer structures on the graphene surfaces. Moreover, they share common vibration properties close to the surfaces and thus couple with the vibrations of graphene in a similar way. These common properties explain why they yield similar Kapitza conductance values despite their bulk thermal conductivity differing by more than a factor of two.

In order to study the phenomena where interfaces play a dominant role through molecular simulations, the proper representation of the interfacial thermodynamic properties of a given model is of crucial importance. The use of coarse-grained rather than atomistic models makes it possible to simulate interfacial systems with larger time and length scales. In the present work, we compare the structure and thermodynamic behavior of one atomistic and two single-site coarse-grained models of water on nonpolar surfaces, namely, graphite and the basal plane of molybdenum disulfide. The three models interact with the surfaces through Lennard-Jones potentials parametrized to reproduce recent experimental contact angle measurements. The models form a layered structure close to the surface, which is usually observed on sufficiently attractive nonpolar substrates. However, differences in the structure and thermodynamic behavior are observed between the models. These differences are explained by certain features of the water models, such as short range tetrahedral order and liquid density fluctuations. Besides these results, the approach employed in the present study may be used to assess the ability of coarse-grained models for solid-liquid systems to represent consistent interfacial thermodynamics.

We study the role of solid–liquid interface thermal resistance (Kapitza resistance) on the evaporation rate of droplets on a heated surface by using a multiscale combination of molecular dynamics (MD) simulations and analytical continuum theory. We parametrize the nonbonded interaction potential between perfluorohexane (C₆F₁₄) and a face-centered-cubic solid surface to reproduce the experimental wetting behavior of C₆F₁₄ on black chromium through the solid–liquid work of adhesion (quantity directly related to the wetting angle). The thermal conductances between C₆F₁₄ and (100) and (111) solid substrates are evaluated by a nonequilibrium molecular dynamics approach for a liquid pressure lower than 2 MPa. Finally, we examine the influence of the Kapitza resistance on evaporation of droplets in the vicinity of a three-phase contact line with continuum theory, where the thermal resistance of liquid layer is comparable with the Kapitza resistance. We determine the thermodynamic conditions under which the Kapitza resistance plays an important role in correctly predicting the evaporation heat flux.

When droplets of nanoparticle suspension evaporate from surfaces, they leave behind a deposit of nanoparticles. The mechanism of evaporation-induced pattern formation in the deposit is studied by molecular dynamics simulations for sessile nanodroplets. The influence of the interaction between nanoparticles and liquid molecules and the influence of the evaporation rate on the final deposition pattern are addressed. When the nanoparticle-liquid interaction is weaker than the liquid-liquid interaction, an interaction-driven or evaporation-induced layer of nanoparticles appears at the liquid-vapor interface and eventually collapses onto the solid surface to form a uniform deposit independently of the evaporation rate. When the nanoparticle-liquid and liquid-liquid interactions are comparable, the nanoparticles are dispersed inside the droplet and evaporation takes place with the contact line pinned at a surface defect. In such a case, a pattern with an approximate ring-like shape is found with fast evaporation, while a more uniform distribution is observed with slower evaporation. When the liquid-nanoparticle interaction is stronger than the liquid-liquid interaction, evaporation always occurs with receding contact line. The final deposition pattern changes from volcano-like to pancake-like with decreasing evaporation rate. These findings might help to design nanoscale structures like nanopatterns or nanowires on surface through controlled solvent evaporation.

Owing to its peculiar electronic properties, molybdenum disulfide (MoS₂) has been the subject of a growing number of studies in the recent years. In applications, this material and other transition metal dichalcogenides (TMDs) may have to interact with a liquid or polymer phase as well as solutions of biomolecules. It is therefore of primary importance to understand the wetting and adhesion properties of TMDs. Starting from existing models, we derive Lennard-Jones parameters for the interaction between water and the basal plane of MoS₂ that are consistent with recent wetting experiments. Molecular dynamics simulations indicate that a stack of only two MoS₂ monolayers is necessary to capture the wetting behavior of bulk MoS₂. It is found that the Coulomb interaction between water and monolayer and bilayer MoS₂ plays no role in the related interfacial thermodynamics. Calculations with the optimized parameters show that the depth of the well of the interaction potential between water and bulk MoS₂ is of the order of 8.2 kJ/mol. Such a value is comparable with what was found for graphite and consistent with the fact that the wetting angles of water on graphite and MoS₂ are almost equal. The derivation of the force-field parameters is performed using a methodology which, contrary to previous studies, makes a consistent use of droplet calculations. The results of our work should find application in further simulation studies on the wetting behavior of TMDs and other dispersive materials.

More than two decades ago, the elusiveness of a liquid-vapor equilibrium and a corresponding critical point in simulations of the supposedly simple model of dipolar hard spheres came as a surprise to many liquid matter theorists. van Leeuwen and Smit [Phys. Rev. Lett. 71, 3991 (1993)] showed that a minimum of attractive dispersion interactions among the dipolar particles may be needed to observe regular fluid behavior. Here, we adopt their approach and use an only slightly modified model, in which the original point dipole is replaced by a dipole moment produced by charges that are separated in space, to study the influence of dispersion interactions of variable strength on the coexistence and interfacial properties of a polar fluid. The thermophysical properties are discussed in terms of Guggenheim’s corresponding-states approach. In this way, the coexistence curve, the critical compressibility factor, the surface tension, Guggenheim’s ratio, and modifications of Guldberg’s and Trouton’s rules (related to the vapor pressure and the enthalpy of vaporization) are analyzed. As the importance of dispersion is decreased, a crossover from simple-fluid behavior to that characteristic of strongly dipolar systems takes place; for some properties, this transition is monotonic, but for others it occurs non-monotonically. For strongly dipolar systems, the reduced surface tension is very low, whereas Guggenheim’s ratio and Guldberg’s ratio are found to be high. The critical compressibility factor is smaller, and the coexistence curve is wider and more skewed than for simple fluids. For very weak dispersion, liquid-vapor equilibrium is still observable, but the interfacial tension is extremely low and may, eventually, vanish marking the end of the existence of a liquid phase. We discuss the implications of our findings for real fluids, in particular, for hydrogen fluoride.

The electrowetting behavior of 1-butyl-3-methylimidazolium tetrafluoroborate ([bmim][BF4]) confined between two oppositely charged graphene layers is investigated using molecular dynamics simulations. By introducing charges on the surface counterions are attracted to the surface and co-ions are repelled from it leading to the reduction of the solid-liquid interfacial free energy and consequently the contact angle. Recently, we have shown that changes in the contact angle upon charging the surface are asymmetric with respect to surface polarity, and opposite to the changes in the solid-liquid interfacial free energy. In this work, the asymmetry of the solid-liquid interfacial free energy is shown to originate from differences in structural organization of the ions at the interface, with positively polarized surfaces inducing a more favorable electrostatic arrangement of the ions. Analysis of the liquid structure in the vicinity of the three phase contact line however shows that the ion size-asymmetry, toge...