Experiments and atomistic simulations have suggested the existence of a direct correlation between the wetting properties of a surface and heat transfer across it. In this investigation, molecular dynamics simulations of surface wettability and solid-liquid thermal transport were conducted in order to better understand the relationship between the surface chemistry and thermal transport. The wettability transparency of graphene-coated surfaces was considered in order to investigate heat transfer across a complex interface with similar wettability as a bare surface. The results indicate that the relationship between the interfacial heat transfer and wettability is not universal. The density depletion length was found to reconcile the thermal boundary conductance calculations for different bare and graphene-coated silicon surfaces.

The wettability of graphitic carbon and silicon surfaces was numerically and theoretically investigated. A multi-response method has been developed for the analysis of conventional molecular dynamics (MD) simulations of droplets wettability. The contact angle and indicators of the quality of the computations are tracked as a function of the data sets analyzed over time. This method of analysis allows accurate calculations of the contact angle obtained from the MD simulations. Analytical models were also developed for the calculation of the work of adhesion using the mean-field theory, accounting for the interfacial entropy changes. A calibration method is proposed to provide better predictions of the respective contact angles under different solid-liquid interaction potentials. Estimations of the binding energy between a water monomer and graphite match those previously reported. In addition, a breakdown in the relationship between the binding energy and the contact angle was observed. The macroscopic contact angles obtained from the MD simulations were found to match those predicted by the mean-field model for graphite under different wettability conditions, as well as the contact angles of Si(100) and Si(111) surfaces. Finally, an assessment of the effect of the Lennard-Jones cutoff radius was conducted to provide guidelines for future comparisons between numerical simulations and analytical models of wettability.

Equilibrium and nonequilibrium molecular dynamics simulations were conducted in order to evaluate the hypothesis that the hydrodynamic slip length is a surface property. The system under investigation was water confined between two graphite layers to form nanochannels of different sizes (3-8 nm). The water-carbon interaction potential was calibrated by matching wettability experiments of graphitic-carbon surfaces free of airborne hydrocarbon contamination. Three equilibrium theories were used to calculate the hydrodynamic slip length. It was found that one of the recently reported equilibrium theories for the calculation of the slip length featured confinement effects, while the others resulted in calculations significantly hindered by the large margin of error observed between independent simulations. The hydrodynamic slip length was found to be channel-size independent using equilibrium calculations, i.e., suggesting a consistency with the definition of a surface property, for 5-nm channels and larger. The analysis of the individual trajectories of liquid particles revealed that the reason for observing confinement effects in 3-nm nanochannels is the high mobility of the bulk particles. Nonequilibrium calculations were not consistently affected by size but by noisiness in the smallest systems.

No abstract

In order to better understand the behavior and governing characteristics of the wetting transparency phenomenon observed in graphene-coated surfaces, molecular dynamics simulations were coupled with a theoretical model. Graphene-coated silicon was selected for this analysis, due to potential applications of hybrid silicon-graphene materials as detectors in aqueous environments. The results indicate good agreement between the theory and simulations at the macroscopic conditions required to observe wetting transparency. A microscopic analysis was also conducted in order to identify the parameters, such as the interaction potential energy landscape and the interfacial liquid structure that govern the wetting behavior of graphene-coated surfaces. The interfacial liquid structure was found to be different between uncoated Si(100) and the graphene-coated version and very similar between uncoated Si(111) and the graphene-coated version. However, the concentration of liquid particles for both silicon surfaces was found to be very similar under transparent wetting conditions.