J. S. Hansen

How fast does water flow in carbon nanotubes

The purpose of this paper is threefold. First, we review the existing literature on flow rates of water in carbon nanotubes. Data for the slip length which characterizes the flow rate are scattered over 5 orders of magnitude for nanotubes of diameter 0.81-10 nm. Second, we precisely compute the slip length using equilibrium molecular dynamics (EMD) simulations, from which the interfacial friction between water and carbon nanotubes can be found, and also via external field driven non-equilibrium molecular dynamics simulations (NEMD). We discuss some of the issues in simulation studies which may be reasons for the large disagreements reported. By using the EMD method friction coefficient to determine the slip length, we overcome the limitations of NEMD simulations. In NEMD simulations, for each tube we apply a range of external fields to check the linear response of the fluid to the field and reliably extrapolate the results for the slip length to values of the field corresponding to experimentally accessible pressure gradients. Finally, we comment on several issues concerning water flow rates in carbon nanotubes which may lead to some future research directions in this area.

Slip length of water on graphene: Limitations of non-equilibrium molecular dynamics simulations

Data for the flow rate of water in carbon nanopores is widely scattered, both in experiments and simulations. In this work, we aim at precisely quantifying the characteristic large slip length and flow rate of water flowing in a planar graphene nanochannel. First, we quantify the slip length using the intrinsic interfacial friction coefficient between water and graphene, which is found from equilibrium molecular dynamics (EMD) simulations. We then calculate the flow rate and the slip length from the streaming velocity profiles obtained using non-equilibrium molecular dynamics (NEMD) simulations and compare with the predictions from the EMD simulations. The slip length calculated from NEMD simulations is found to be extremely sensitive to the curvature of the velocity profile and it possesses large statistical errors. We therefore pose the question: Can a micrometer range slip length be reliably determined using velocity profiles obtained from NEMD simulations? Our answer is not practical, if not impossible based on the analysis given as the results. In the case of high slip systems such as water in carbon nanochannels, the EMD method results are more reliable, accurate, and computationally more efficient compared to the direct NEMD method for predicting the nanofluidic flow rate and hydrodynamic boundary condition.

Interfacial slip friction at a fluid-solid cylindrical boundary

Recently we proposed a method to calculate the interfacial friction coefficient between fluid and solid at a planar interface. In this work we extend the method to cylindrical systems where the friction coefficient is curvature dependent. We apply the method to methane flow in carbon nanotubes, and find good agreement with non-equilibrium molecular dynamics simulations. The proposed method is robust, general, and can be used to predict the slip for cylindrical nanofluidic systems.

Electropumping of water with rotating electric fields

Pumping of fluids confined to nanometer dimension spaces is a technically challenging yet vitally important technological application with far reaching consequences for lab-on-a-chip devices, biomimetic nanoscale reactors, nanoscale filtration devices and the like. All current pumping mechanisms require some sort of direct intrusion into the nanofluidic system, and involve mechanical or electronic components. In this paper, we present the first nonequilibrium molecular dynamics results to demonstrate that non-intrusive electropumping of liquid water on the nanoscale can be performed by subtly exploiting the coupling of spin angular momentum to linear streaming momentum. A spatially uniform rotating electric field is applied to water molecules, which couples to their permanent electric dipole moments. The resulting molecular rotational momentum is converted into linear streaming momentum of the fluid. By selectively tuning the degree of hydrophobicity of the solid walls one can generate a net unidirectional flow. Our results for the linear streaming and angular velocities of the confined water are in general agreement with the extended hydrodynamical theory for this process, though also suggest refinements to the theory are required. These numerical experiments confirm that this new concept for pumping of polar nanofluids can be employed under laboratory conditions, opening up significant new technological possibilities.

Molecular dynamics study of nanoconfined water flow driven by rotating electric fields under realistic experimental conditions

In our recent work, J. Chem. Phys. 2013, 138, 154712, we demonstrated the feasibility of unidirectional pumping of water, exploiting translational-rotational momentum coupling using nonequilibrium molecular dynamics simulations. Flow can be sustained when the fluid is driven out of equilibrium by an external spatially uniform rotating electric field and confined between two planar surfaces exposing different degrees of hydrophobicity. The permanent dipole moment of water follows the rotating field, thus inducing the molecules to spin, and the torque exerted by the field is continuously injected into the fluid, enabling a steady conversion of spin angular momentum into linear momentum. The translational-rotational coupling is a sensitive function of the rotating electric field parameters. In this work, we have found that there exists a small energy dissipation region attainable when the frequency of the rotating electric field matches the inverse of the dielectric relaxation time of water and when its amplitude lies in a range just before dielectric saturation effects take place. In this region, that is, when the frequency lies in a small window of the microwave region around ∼20 GHz and amplitude ∼0.03 V Å(-1), the translational-rotational coupling is most effective, yielding fluid velocities of magnitudes of ∼2 ms(-1) with only moderate fluid heating. In this work, we also confine water to a realistic nanochannel made of graphene giving a hydrophobic surface on one side and β-cristobalite giving a hydrophilic surface on the other, reproducing slip-and-stick velocity boundary conditions, respectively. This enables us to demonstrate that in a realistic environment, the coupling can be effectively exploited to achieve noncontact pumping of water at the nanoscale. A quantitative comparison between nonequilibrium molecular dynamics and analytical solutions of the extended Navier-Stokes equations, including an external rotating electric field has been performed, showing excellent agreement when the electric field parameters match the aforementioned small energy dissipation region.

Structural and dynamical properties for confined polymers undergoing planar Poiseuille flow

The authors present the results from nonequilibrium molecular dynamics simulations for the structural and dynamical properties of highly confined linear polymer fluids undergoing planar Poiseuille flow. They study systems confined within pores of several atomic diameters in width and investigate the dependence of the density profiles, the mean squared radius of gyration, the mean squared end-to-end distance, streaming velocity, strain rate, shear stress, and streaming angular velocity as functions of average fluid density and chain length. Their simulation results show that, sufficiently far from the walls, the radius of gyration for molecules under shear in the middle of the pore follows the power law Rg=ANbnu, where Nb is the number of bonds and the exponent has a value of 0.5 which resembles the value for a homogeneous equilibrium fluid. Under the conditions simulated, the authors find the onset of flat velocity profiles but with very little wall slippage. These flat profiles are most likely due to the restricted layering of the fluid into just one or two molecular layers for narrow pore widths compared to chain length, rather than typical plug-flow conditions. The angular velocity is shown to be proportional to half the strain rate in the pore interior when the chain length is sufficiently small compared to the pore width, consistent with the behavior for homogeneous fluids in the linear regime.

A new and effective method for thermostatting confined fluids

We present a simple thermostatting method suitable for nanoconfined fluid systems. Two conventional strategies involve thermostatting the fluid directly or employing a thermal wall that couples only the wall atoms with the thermostat. When only a thermal wall is implemented, the temperature control of the fluid is true to the actual experiment and the heat is transferred from the fluid to the walls. However, for large or complex systems it can often be computationally prohibitive to employ thermal walls. To overcome this limitation many researchers choose to freeze wall atoms and instead apply a synthetic thermostat to the fluid directly through the equations of motion. This, however, can have serious consequences for the mechanical, thermodynamic, and dynamical properties of the fluid by introducing unphysical behaviour into the system [Bernardi et al., J. Chem. Phys. 132, 244706 (2010)]. In this paper, we propose a simple scheme which enables working with both frozen walls and naturally thermostatted liquids. This is done by superimposing the walls with oscillating particles, which vibrate on the edge of the fluid control volume. These particles exchange energy with the fluid molecules, but do not interact with wall atoms or each other, thus behaving as virtual particles. Their displacements violate the Lindemann criterion for melting, in such a way that the net effect would not amount to an additional confining surface. One advantage over standard techniques is the reduced computational cost, particularly for large walls, since they can be kept rigid. Another advantage over accepted strategies is the opportunity to freeze complex charged walls such as β-cristobalite. The method furthermore overcomes the problem with polar fluids such as water, as thermalized charged surfaces require higher spring constants to preserve structural stability, due to the effects of strong Coulomb interactions, thus inevitably degrading the thermostatting efficiency.

Local linear viscoelasticity of confined fluids.

In this paper the authors propose a novel method to study the local linear viscoelasticity of fluids confined between two walls. The method is based on the linear constitutive equation and provides details about the real and imaginary parts of the local complex viscosity. They apply the method to a simple atomic fluid undergoing zero mean oscillatory flow using nonequilibrium molecular dynamics simulations. The method shows that the viscoelastic properties of the fluid exhibit dramatic spatial changes near the wall-fluid boundary due to the high density in this region. It is also shown that the real part of the viscosity converges to the frequency dependent local shear viscosity sufficiently far away from the wall. This also provides valuable information about the transport properties in the fluid, in general. The viscosity is compared with predictions from the local average density model. The two methods disagree in that the local average density model predicts larger viscosity variations near the wall-fluid boundary than what is observed through the method presented here.

An extended analysis of the viscosity kernel for monatomic and diatomic fluids

We present an extended analysis of the wavevector dependent shear viscosity of monatomic and diatomic (liquid chlorine) fluids over a wide range of wavevectors and for a variety of state points. The analysis is based on equilibrium molecular dynamics simulations, which involve the evaluation of transverse momentum density and shear stress autocorrelation functions. For liquid chlorine we present the results in both atomic and molecular formalisms. We find that the viscosity kernel of chlorine in the atomic representation is statistically indistinguishable from that in the molecular representation. The results further suggest that the real space viscosity kernels of monatomic and diatomic fluids depend sensitively on the density, the potential energy function and the choice of fitting function in reciprocal space. It is also shown that the reciprocal space shear viscosity data can be fitted to two different simple functional forms over the entire density, temperature and wavevector range: a function composed of n-Gaussian terms and a Lorentzian-type function. Overall, the real space viscosity kernel has a width of 3-6 atomic diameters, which means that the generalized hydrodynamic constitutive relation is required for fluids with strain rates that vary nonlinearly over distances of the order of atomic dimensions.

Rotational viscosity of fluids composed of linear molecules: an equilibrium molecular dynamics study

In this paper, we investigate the rotational viscosity for a chlorine fluid and for a fluid composed of small linear molecules by using equilibrium molecular dynamics simulations. The rotational viscosity is calculated over a large range of state points. It is found that the rotational viscosity is almost independent of temperature in the range studied here but exhibits a power-law dependency on density. The rotational viscosity also shows a power-law relationship with the molecular length, and the ratio between the shear and rotational viscosities approaches 0.5 for the longest molecule studied here. By changing the number of atoms or united atomic units per molecule and by keeping the molecule length fixed, we show that fluids composed of molecules which have a rodlike shape have a lower rotational viscosity. We argue that this phenomenon is due to the reduction in intermolecular connectivity, which leads to larger fluctuations around the values possessed by the fluid on average. The conclusions here can be extended to fluids composed of uniaxial molecules of arbitrary length.