Water flow in graphene nanochannels driven by imposed thermal gradients: the role of flexural phonons.

Accurate control of fluid transport in nanoscale structures is key to enable the design of foreseeable nanofluidic devices with applications in many fields such as chip cooling, energy conversion, drug delivery and medical diagnosis. Here, inspired by the experimental observation of intrinsic thermal ripples in graphene and by recent advances in the manipulation of 2D nanomaterials, we introduce a graphene-based thermal nanopump which produces controlled and continuous liquid flow in nanoslit channels. We investigate the performance of this thermal nanopump employing large scale molecular dynamics simulations. Upon systematically imposing thermal gradients, a net water flow towards the low-temperature zone is observed, achieving flow velocities up to 4 m s-1. We observe that water flow rates increase monotonically due to larger ripple fluctuations on the graphene layers as higher thermal gradients are applied. Moreover, we find that the out-of-plane flexural phonons in graphene are responsible for flow generation wherein lower frequency phonon branches are activated with higher imposed thermal gradients. Furthermore, by modifying the wettability of the channel walls, an increase of 50% in the water flow rates is observed, showing that the efficiency of the proposed thermal pump can be enhanced by tuning the channel wall hydrophobicity. Our results indicate that thermal gradients can be employed to drive continuous water flow in graphene nanoslit channels with potential applications in nanofluidic devices.

Effect of charge inversion on nanoconfined flow of multivalent ionic solutions.

A comprehensive understanding of fluid dynamics of dilute electrolyte solutions in nanoconfinement is essential to develop more efficient nanofluidic devices. In nanoconduits, the electrical double layer can occupy a considerable part of the channel cross-section, therefore the transport properties of a nanoconfined electrolyte solution can be altered by interfacial phenomena such as the charge inversion (CI). CI is an electrokinetic effect that has been associated with the presence of hydrated multivalent cations in nanoconfinement. Here, all-atom molecular dynamics simulations are employed to study the structure and dynamics of aqueous multivalent electrolyte solutions within slit-shaped silica channels. All simulations are conducted for more than 100 ns to capture the equilibrium ion distribution, the interfacial hydrodynamic properties, and to reveal the influence of CI on nanoconfined fluid transport. The electrolyte solutions consist of water as solvent, chloride as co-ion and different counter-ions, i.e., sodium, magnesium and aluminum. We find that the interfacial viscosity is related to the concentration and valence of the counter-ions in the solution. Our results suggest that higher CI is correlated to the presence of a layer of fluid with augmented viscosity adjacent to the channel wall. As the thickness of this interfacial high-viscosity fluid increases, lower flow rates are measured whereas higher interfacial viscosities and friction coefficients are computed.

Effect of an external electric field on capillary filling of water in hydrophilic silica nanochannels.

Development of functional nanofluidic devices requires understanding the fundamentals of capillary driven flow in nanochannels. In this context, we conduct molecular dynamics simulations of water capillary imbibition in silica nanoslits under externally applied electric (E) fields with strengths between 0 and 1 V nm-1. For increasing E-fields, we observe a systematic lowering in the meniscus contact angle and a decrease in the corresponding water filling rates. These results contrast markedly the classical Washburn-Bosanquet's equation which predicts an increase in filling rates for lower water contact angles. Our study provides evidence that the observed decrease in water filling rates can be attributed to the interplay between two underlying mechanisms, a reduced fluidity of interfacial water and a systematic alignment of the water molecules in the bulk as a response to the particular strength of the applied E-field. Therefore, during water capillary filling a constant E-field applied in the direction parallel to the water imbibition leads to a lower than expected filling rate caused by a viscosity increase in the bulk and an altered solid-liquid friction on the channel walls. These coupled mechanisms governing capillarity under the action of applied E-fields could be manipulated for controlling imbibition of polar liquid solutions in nanofluidic devices.

Water thermophoresis in carbon nanotubes: the interplay between thermophoretic and friction forces.

Thermophoresis is the phenomenon wherein particles experience a net drift induced by a thermal gradient. In this work, molecular dynamics simulations are conducted to study with atomistic detail the thermophoresis of water nanodroplets inside carbon nanotubes (CNTs) and its interplay with the retarding liquid-solid friction. Different applied temperatures, thermal gradients, and droplet sizes are used to reveal the dynamics of the two kinetic regimes of the thermophoretic motion in CNTs. The results indicate that during the droplet motion, the thermophoretic force is independent of the velocity of the droplet, whereas the magnitude of the retarding friction force exhibits a linear dependence. In fact, in the initial regime the magnitude of the friction force increases linearly with the droplet velocity, until the thermophoretic force is balanced by the friction force as the droplet reaches its terminal velocity in the final regime. In addition, an increase in the magnitude of the thermophoretic force is found for longer water droplets. These findings provide a deeper understanding of liquid transport driven by temperature gradients in nanoconfined geometries where liquid-solid interfaces govern fluidics.

Carbon Nanotubes as Thermally Induced Water Pumps.

Thermal Brownian motors (TBMs) are nanoscale machines that exploit thermal fluctuations to provide useful work. We introduce a TBM-based nanopump which enables continuous water flow through a carbon nanotube (CNT) by imposing an axial thermal gradient along its surface. We impose spatial asymmetry along the CNT by immobilizing certain points on its surface. We study the performance of this molecular motor using molecular dynamics (MD) simulations. From the MD trajectories, we compute the net water flow and the induced velocity profiles for various imposed thermal gradients. We find that spatial asymmetry modifies the vibrational modes of the CNT induced by the thermal gradient, resulting in a net water flow against the thermal gradient. Moreover, the kinetic energy associated with the thermal oscillations rectifies the Brownian motion of the water molecules, driving the flow in a preferred direction. For imposed thermal gradients of 0.5-3.3 K/nm, we observe continuous net flow with average velocities up to 5 m/s inside CNTs with diameters of 0.94, 1.4, and 2.0 nm. The results indicate that the CNT-based asymmetric thermal motor can provide a controllable and robust system for delivery of continuous water flow with potential applications in integrated nanofluidic devices.

Slip divergence of water flow in graphene nanochannels: the role of chirality.

Graphene has attracted considerable attention due to its characteristics as a 2D material and its fascinating properties, providing a potential building block for nanofabrication. In nanochannels the solid-liquid interface plays a non-negligible role in determining the fluid dynamics. Therefore, for an optimal design of nanofluidic devices, a comprehensive understanding of the slippage in a water flow confined between graphene walls is important. In nanoconfinement, experimental and computational studies have found the slip length to increase nonlinearly when the shear rate is larger than a critical value. Here, by conducting molecular dynamics simulations, we study the influence of the graphene crystallographic orientation on the slip boundary conditions inside a nanoslit channel. The flow in channels with heights of 2.0, 2.4 and 2.8 nm is driven parallel to the zig-zag and arm-chair crystallographic directions. We extract flow rates, velocity profiles, slip velocities and slip lengths. The slip velocity displays a linear relationship to the shear stress up to a critical value, which is not size dependent. Moreover, the slip length is found to be shear stress dependent above a critical shear stress value of 0.4 MPa. Furthermore, our results indicate that after this critical shear stress is reached, the flow rates are significantly influenced (up to 10%) by the particular orientation of the graphene topology.

Effect of the meniscus contact angle during early regimes of spontaneous imbibition in nanochannels.

Nanoscale capillarity has been extensively investigated; nevertheless, many fundamental questions remain open. In spontaneous imbibition, the classical Lucas-Washburn equation predicts a singularity as the fluid enters the channel consisting of an anomalous infinite velocity of the capillary meniscus. Bosanquet's equation overcomes this problem by taking into account fluid inertia predicting an initial imbibition regime with constant velocity. Nevertheless, the initial constant velocity as predicted by Bosanquet's equation is much greater than those observed experimentally. In the present study, large scale atomistic simulations are conducted to investigate capillary imbibition of water in slit silica nanochannels with heights between 4 and 18 nm. We find that the meniscus contact angle remains constant during the inertial regime and its value depends on the height of the channel. We also find that the meniscus velocity computed at the channel entrance is related to the particular value of the meniscus contact angle. Moreover, during the subsequent visco-inertial regime, as the influence of viscosity increases, the meniscus contact angle is found to be time dependent for all the channels under study. Furthermore, we propose an expression for the time evolution of the dynamic contact angle in nanochannels which, when incorporated into Bosanquet's equation, satisfactorily explains the initial capillary rise.

Wall embedded electrodes to modify electroosmotic flow in silica nanoslits.

Electroosmotic flow in a silica slit channel with nonuniform surface charge density is investigated. In nanoconfinement, the electrical double layer occupies a non-negligible fraction of the system. Therefore, modifying the charge density on specific locations on the channel wall surface allows effective manipulation of the electroosmotic flow rates. In the present study, extensive (160 ns) nonequilibrium molecular dynamics simulations are conducted to investigate the ability of controlling the electroosmotic flow control in a nanoslit by patterning the surface potential. The mechanism to modify the surface charge consists of a set of charged electrodes embedded within one of the channel walls. The presence of the embedded electrodes results in the redistribution of ions in the electrolyte solution and in the alteration of the electroosmotic flow throughout the nanochannel. Indeed, the results reveal significant changes in the electroosmotic driving force and velocity profiles including local flow reversal. This study provides physical insight into the direct manipulation of the electrokinetic flow in nanoslits.

Early regimes of water capillary flow in slit silica nanochannels.

Molecular dynamics simulations are conducted to investigate the initial stages of spontaneous imbibition of water in slit silica nanochannels surrounded by air. An analysis is performed for the effects of nanoscopic confinement, initial conditions of liquid uptake and air pressurization on the dynamics of capillary filling. The results indicate that the nanoscale imbibition process is divided into three main flow regimes: an initial regime where the capillary force is balanced only by the inertial drag and characterized by a constant velocity and a plug flow profile. In this regime, the meniscus formation process plays a central role in the imbibition rate. Thereafter, a transitional regime takes place, in which, the force balance has significant contributions from both inertia and viscous friction. Subsequently, a regime wherein viscous forces dominate the capillary force balance is attained. Flow velocity profiles identify the passage from an inviscid flow to a developing Poiseuille flow. Gas density profiles ahead of the capillary front indicate a transient accumulation of air on the advancing meniscus. Furthermore, slower capillary filling rates computed for higher air pressures reveal a significant retarding effect of the gas displaced by the advancing meniscus.

Thermophoretic motion of water nanodroplets confined inside carbon nanotubes.

We study the thermophoretic motion of water nanodroplets confined inside carbon nanotubes using molecular dynamics simulations. We find that the nanodroplets move in the direction opposite the imposed thermal gradient with a terminal velocity that is linearly proportional to the gradient. The translational motion is associated with a solid body rotation of the water nanodroplet coinciding with the helical symmetry of the carbon nanotube. The thermal diffusion displays a weak dependence on the wetting of the water-carbon nanotube interface. We introduce the use of the moment scaling spectrum (MSS) in order to determine the characteristics of the motion of the nanoparticles inside the carbon nanotube. The MSS indicates that affinity of the nanodroplet with the walls of the carbon nanotubes is important for the isothermal diffusion and hence for the Soret coefficient of the system.