Hydrodynamic slip characteristics of shear-driven water flow in nanoscale carbon slits.

This paper reports on the effects of shear rate and interface modeling parameters on the hydrodynamic slip length (LS) for water-graphite interfaces calculated using non-equilibrium molecular dynamics. Five distinct non-bonded solid-liquid interaction parameters were considered to assess their impact on LS. The interfacial force field derivations included sophisticated electronic structure calculation-informed and empirically determined parameters. All interface models exhibited a similar and bimodal LS response when varying the applied shear rate. LS in the low shear rate regime (LSR) is in good agreement with previous calculations obtained through equilibrium molecular dynamics. As the shear rate increases, LS sharply increases and asymptotes to a constant value in the high shear regime (HSR). It is noteworthy that LS in both the LSR and HSR can be characterized by the density depletion length, whereas solid-liquid adhesion metrics failed to do so. For all interface models, LHSR calculations were, on average, ∼28% greater than LLSR, and this slip jump was confirmed using the SPC/E and TIP4P/2005 water models. To address the LS transition from the LSR to the HSR, the viscosity of water and the interfacial friction coefficient were investigated. It was observed that in the LSR, the viscosity and friction coefficient decreased at a similar rate, while in the LSR-to-HSR transition, the friction coefficient decreased at a faster rate than the shear viscosity until they reached a new equilibrium, hence explaining the LS-bimodal behavior. This study provides valuable insights into the interplay between interface modeling parameters, shear rate, and rheological properties in understanding hydrodynamic slip behavior.

Thermal transport across flat and curved gold-water interfaces: Assessing the effects of the interfacial modeling parameters.

The present investigation assesses a variety of parameters available in the literature to model gold-water interfaces using molecular dynamics simulations. The study elucidates the challenges of characterizing the solid-liquid affinity of highly hydrophilic gold-water interfaces via wettability. As an alternative, the local pairwise interaction energy was used to describe the solid-liquid affinity of flat and curved surfaces, where for the latter, the calculation of a contact angle becomes virtually impossible. Regarding the heat transfer properties of different interface models (flat and curved), partly conclusive trends were observed between the total pairwise interaction energy and the thermal boundary conductance. It was observed that the solid surface structure, interfacial force field type, and force field parameters created a characteristic bias in the interfacial water molecules (liquid structuring). Consequently, a study of the liquid depletion layer provided better insight into the interfacial heat transfer among different interfaces. By computing the density depletion length, which describes the deficit or surplus of energy carries (water molecules) near the interface, a proper characterization of the thermal boundary conductance was obtained for the different gold-water interfaces. It was observed that the interfacial heat transfer is favored when the water molecules organize in cluster-like structures near the interface, by a surplus of water molecules at the interface, i.e., lower density depletion length, and by the closeness of water to the solid atoms.

Effects of Moisture and Synthesis-Derived Contaminants on the Mechanical Properties of Graphene Oxide: A Molecular Dynamics Investigation.

This paper reports on the effects of the chemical composition of graphene oxide (GO) sheets on the mechanical properties of bulk GO. Three key factors were analyzed: (i) the oxygenated functional groups' concentration, (ii) the content of intersheet water (moisture), and (iii) the presence of residual contaminants observed from the synthesis of GO. Molecular dynamics simulations using the reactive force field ReaxFF were conducted to model tensile strength, indentation, and shear stress tests. The structural integrity of the carbon basal plane was the primary variable that determined mechanical behavior of GO slabs. Hydrogen-bond networks played an essential role in the tensile fracture mechanism, delaying the onset of fracture whenever strong hydrogen bonds existed in the intersheet space. The presence of interlayer sulfate ion contaminants negatively impacted the tensile strength, stiffness, and toughness of GO. Moreover, it was observed that intersheet sulfate ions improved the resistance to fracture of GO at low sulfur concentrations, while lower fracture strains were observed beyond a critical concentration. Alike the tensile stress findings, the indentation properties were determined by the integrity of the carbon basal plane. Our findings agree with experimental mechanical property measurements and reveal the importance of considering synthesis-derived contaminants in molecular models of GO.

Effects of the Interfacial Modeling Approach on Equilibrium Calculations of Slip Length for Nanoconfined Water in Carbon Slits.

In this investigation, equilibrium molecular dynamics simulations were conducted to assess the influence of the interface modeling approach on the calculation of hydrodynamic slip in carbon nanochannels. A Green-Kubo formalism was implemented for the calculation of the slip length in water confined by graphite layers. The nonbonded interactions between solid and liquid atoms (interface models) were modeled using parameters optimized to represent the wetting behavior and adsorption energy curves from electronic structure calculations. Conventional carbon-oxygen-only interaction models were compared against comprehensive models able to represent the molecular-orientation-dependent energy of interaction. Quasi-universal relationships built under the premise of the slip length dependence on the water-graphite affinity and characterized by macroscopic wettability were critically assessed. It was found that the wetting behavior cannot fully characterize the hydrodynamic slip because interface models that produced the same surface wettability yielded different values of the friction coefficient. Alternatively, the density depletion length, used to characterize the interfacial liquid structuring and the availability of momentum carriers (interfacial water molecules), was able to accurately represent the slip length trends independently of the interface model. These findings reassert the importance of physically sound interface models to study interfacial transport properties and the need of reliable parameters and characterization procedures to support theoretical models that seek to unveil the inconsistencies in hydrodynamic slip calculations.