Combining molecular dynamics (MD) and continuum simulations, we study the dynamics of propagation of a peeling front in a system composed of multilayered graphene nanosheets completely immersed in water. Peeling is induced by lifting one of the nanosheet edges with an assigned pulling velocity normal to the flat substrate. Using MD, we compute the pulling force as a function of the pulling velocity, and quantify the viscous resistance to the advancement of the peeling front. We compare the MD results to a 1D continuum model of a sheet loaded with modelled hydrodynamic loads. Our results show that the viscous dependence of the force on the velocity is negligible below a threshold velocity. Above this threshold, the hydrodynamics is mainly controlled by the viscous resistance associated to the flow near the crack opening, while lubrication forces are negligible owing to the large hydrodynamic slip at the liquid-solid boundary. Two dissipative mechanisms are identified: a drag resistance to the upward motion of the edge, and a resistance to the gap opening associated to the curvature of the flow streamlines near the entrance. Surprisingly, the shape of the sheet was found to be approximately independent of the pulling velocity even for the largest velocities considered.

@en

Rigid platelike particles displaying interfacial slip can attain a constant orientation in a shear flow when the slip length is sufficiently large. But actual thin particles such as single-layer graphene are flexible and prone to bending deformations when exposed to shear stress. To study the effect of bending deformation on the dynamics of flexible platelike particles with large interfacial slip in a shear flow, we develop a two-dimensional (2D) fluid-structure interaction model. Our model is based on coupling the Euler-Bernoulli beam equation with a boundary integral method to solve the hydrodynamic stress at the particle surface. Emphasis is placed on resolving accurately the stress distribution at the edges of the particle. We find that (i) a stable alignment occurs even for relatively flexible particles and that (ii) edges effects on the shape of the plate are important for values of the length-to-thickness aspect ratio as large as 100. Our results are particularly relevant in view of recent research on the hydrodynamics of suspended flexible sheets made of 2D nanomaterials.

@en

The classical theory by Jeffery predicts that, in the absence of Brownian fluctuations, a thin rigid platelet rotates continuously in a shear flow, performing periodic orbits. However, a stable orientation is possible if the surface of the platelet displays a hydrodynamic slip length comparable to or larger than the thickness of the platelet. In this article, by solving the Fokker-Plank equation for the orientation distribution function and corroborating the analysis with boundary integral simulations, we quantify a threshold Péclet number, above which such alignment occurs. We found that for smaller than, but larger than a second threshold, a regime emerges where Brownian fluctuations are strong enough to break the platelet's alignment and induce rotations, but with a period of rotation that depends on the value of. For below this second threshold, slip has a negligible effect on the orientational dynamics. We use these thresholds to classify the dynamics of graphene-like nanoplatelets for realistic values of and apply our results to the quantification of the orientational contribution to the effective viscosity of a dilute suspension of nanoplatelets with slip. We find a non-monotonic variation of this term, with a minimum occurring when the slip length is comparable to the thickness of the particle.

@en

The adsorption of graphene-oxide (GO) nanoparticles at the interface between water and vapor was analyzed using all-atom molecular simulations for single and multiple particles. For a single GO particle, our results indicate that the adsorption energy does not scale linearly with the surface coverage of oxygen groups, unlike typically assumed for Janus colloids. Our results also show that the surface activity of the particle depends on the number of surface oxygen groups as well as on their distribution: for a given number of oxygen groups, a GO particle with a patched surface was found to be more surface active than a particle with evenly distributed groups. Then, to understand what sets the thickness of GO layers at interfaces, the adsorption energy of a test GO particle was measured in the presence of multiple GO particles already adsorbed at the interface. Our results indicate that in the case of high degree of oxidation, particle-particle interactions at the water-vapor interface hinder the adsorption of the test particle. In the case of a low degree of oxidation, however, clustering and stacking of GO particles dominate the adsorption behavior, and particle-particle interactions favor the adsorption of the test particle. These results highlight the complexity of multiple particle adsorption and the limitations of single-particle adsorption models when applied to GO at a relatively high surface concentration.

@en

Using molecular dynamics simulations we investigate the shear-induced rotational dynamics of a Brownian nanographene (hexabenzocoronene) freely suspended in a liquid. We demonstrate that, owing to a finite hydrodynamic slip at the molecular surface, these flat molecules tend to align with a constant orientation angle instead of performing the classical periodic orbits predicted by Jeffery's theory. Results are extracted for different Péclet numbers and compared to the predictions by a theory developed for a rigid axisymmetric particle with orientation confined to the flow-gradient plane. The theory is based on the resolution of a one-dimensional Fokker-Planck equation for the angle φ made by one of the particle's diameters with the flow direction. Remarkably, our results show that the essential features of the three-dimensional orientational statistics of the nanographene are captured by the one-dimensional model, given that the hydrodynamic velocity is closed in terms of the slip length λ. Finally, we explore the situation in which multiple nanographenes are suspended in the liquid, and show that slip results in a reduction in specific viscosity.

@en

Liquid-phase exfoliation, the use of a sheared liquid to delaminate graphite into few-layer graphene, is a promising technique for the large-scale production of graphene. However, the microscale and nanoscale fluid-structure processes controlling the exfoliation are not fully understood. Here, we perform non-equilibrium molecular dynamics simulations of a defect-free graphite nanoplatelet suspended in a shear flow and measure the critical shear rate γ̇ c needed for the exfoliation to occur. We compare γ̇ c for different solvents, including water and N-methyl-pyrrolidone, and nanoplatelets of different lengths. Using a theoretical model based on a balance between the work done by viscous shearing forces and the change in interfacial energies upon layer sliding, we are able to predict the critical shear rates γ̇ c measured in simulations. We find that an accurate prediction of the exfoliation of short graphite nanoplatelets is possible only if both hydrodynamic slip and the fluid forces on the graphene edges are considered and if an accurate value of the solid–liquid surface energy is used. The commonly used “geometric-mean” approximation for the solid–liquid energy leads to grossly incorrect predictions.@en

The large-scale processing of nanomaterials such as graphene and MoS₂ relies on understanding the flow behaviour of nanometrically-thin platelets suspended in liquids. Here we show, by combining non-equilibrium molecular dynamics and continuum simulations, that rigid nanoplatelets can attain a stable orientation for sufficiently strong flows. Such a stable orientation is in contradiction with the rotational motion predicted by classical colloidal hydrodynamics. This surprising effect is due to hydrodynamic slip at the liquid-solid interface and occurs when the slip length is larger than the platelet thickness; a slip length of a few nanometers may be sufficient to observe alignment. The predictions we developed by examining pure and surface-modified graphene is applicable to different solvent/2D material combinations. The emergence of a fixed orientation in a direction nearly parallel to the flow implies a slip-dependent change in several macroscopic transport properties, with potential impact on applications ranging from functional inks to nanocomposites.

@en