A new and effective method for thermostatting confined fluids

Electropumping of nanofluidic water by linear and angular momentum coupling: theoretical foundations and molecular dynamics simulations

In this article we review the relatively new phenomenon of electropumping in nanofluidic systems, in which nonzero net flow results when polar molecules are rotated by external electric fields. The flow is a consequence of coupling of the spin angular momentum of molecules with their linear streaming momentum. By devising confining surfaces that are asymmetric - specifically one surface is more hydrophobic compared to the other - unidirectional flow results and so pumping can be achieved without the use of pressure gradients. We first cover the historical background to this phenomenon and follow that with a detailed theoretical description of the governing hydrodynamics. Following that we summarise work that has applied this phenomenon to pump water confined to planar nanochannels, semi-functionalised single carbon nanotubes and concentric carbon nanotubes. We also report on the energy efficiency of this pumping technique by comparisons with traditional flows of planar Couette and Poiseuille flow, with the surprising conclusion that electropumping at the nanoscale is some 4 orders of magnitude more efficient than pumping by Poiseuille flow.

Molecular Insights into the Regulatable Interfacial Property and Flow Behavior of Confined Ionic Liquids in Graphene Nanochannels

The understanding of confined structure and flow property of ionic liquid (IL) in a nanochannel are essential for the efficient application of ILs in the green chemical processes. In this work, the ionic structure and various flow behaviors of ILs inside graphene nanochannels via molecular dynamics simulations are shown. The effect of the nanochannel structure on confined flow is explored, showing that the width mainly heightens the viscosity while the oxidation degree primarily enhances the interfacial friction coefficient. Tuning the width and oxidation degree of nanochannel, three different flow behaviors including Poiseuille, partial plunger and full plunger flow can be achieved, where the second one does not occur in water or other organic solvents. To describe the special flow behavior, an effective influence extent of the nanochannel (w _EIE ) is defined, whose value can distinguish the above flows effectively. Based on w _EIE , the phase diagrams of flow behavior for the nanochannel structure and pressure gradient are obtained, showing that the critical pressure gradient decreases with width and increases with the oxidation degree. Based on the quantitative relations between confined structures, viscosity, friction coefficient, flow behavior, and nanochannel structure, the intrinsic mechanism of regulating the flow behavior and rational design of nanochannel are finally discussed.

Electrokinetic flow of an aqueous electrolyte in amorphous silica nanotubes

We study the pressure-driven flow of aqueous NaCl in amorphous silica nanotubes using nonequilibrium molecular dynamics simulations featuring both polarizable and non-polarizable molecular models.

The confined [Bmim][BF4] ionic liquid flow through graphene oxide nanochannels: a molecular dynamics study

Our work quantitatively shows how hydroxyls influence the flow behavior of ionic liquids in nanochannels.

Kinetic-contact-driven gigantic energy transfer in a two-dimensional Lennard-Jones fluid confined to a rotating pore

A two-dimensional Lennard-Jones system in a circular and rotating container has been studied by means of molecular dynamics technique. A nonequilibrium transition to the rotating stage has been detected in a delayed time since an instant switching of the frame rotation. This transition is attributed to the increase of the density at the wall because of the centrifugal force. At the same time the phase transition occurs, the inner system changes its configuration of the solid-state type into the liquid type. Impact of angular frequency and molecular roughness on the transport properties of the nonrotating and rotating systems is analyzed.

Effects of Confinement on the Dielectric Response of Water Extends up to Mesoscale Dimensions

The extent of confinement effects on water is not clear in the literature. While some properties are affected only within a few nanometers from the wall surface, others are affected over long length scales, but the range is not clear. In this work, we have examined the dielectric response of confined water under the influence of external electric fields along with the dipolar fluctuations at equilibrium. The confinement induces a strong anisotropic effect which is evident up to 100 nm channel width, and may extend to macroscopic dimensions. The root-mean-square fluctuations of the total orientational dipole moment in the direction perpendicular to the surfaces is 1 order of magnitude smaller than the value attained in the parallel direction and is independent of the channel width. Consequently, the isotropic condition is unlikely to be recovered until the channel width reaches macroscopic dimensions. Consistent with dipole moment fluctuations, the effect of confinement on the dielectric response also persists up to channel widths considerably beyond 100 nm. When an electric field is applied in the perpendicular direction, the orientational relaxation is 3 orders of magnitude faster than the dipolar relaxation in the parallel direction and independent of temperature.

Hydrodynamic slip in silicon nanochannels

Equilibrium and nonequilibrium molecular dynamics simulations were performed to better understand the hydrodynamic behavior of water flowing through silicon nanochannels. The water-silicon interaction potential was calibrated by means of size-independent molecular dynamics simulations of silicon wettability. The wettability of silicon was found to be dependent on the strength of the water-silicon interaction and the structure of the underlying surface. As a result, the anisotropy was found to be an important factor in the wettability of these types of crystalline solids. Using this premise as a fundamental starting point, the hydrodynamic slip in nanoconfined water was characterized using both equilibrium and nonequilibrium calculations of the slip length under low shear rate operating conditions. As was the case for the wettability analysis, the hydrodynamic slip was found to be dependent on the wetted solid surface atomic structure. Additionally, the interfacial water liquid structure was the most significant parameter to describe the hydrodynamic boundary condition. The calibration of the water-silicon interaction potential performed by matching the experimental contact angle of silicon led to the verification of the no-slip condition, experimentally reported for silicon nanochannels at low shear rates.

Surface-Heating Algorithm for Water at Nanoscale

A novel surface-heating algorithm for water is developed for molecular dynamics simulations. The validated algorithm can simulate the transient behavior of the evaporation of water when heated from a surface, which has been lacking in the literature. In this work, the algorithm is used to study the evaporation of water droplets on a platinum surface at different temperatures. The resulting contact angles of the droplets are compared to existing theoretical, numerical, and experimental studies. The evaporation profile along the droplet's radius and height is deduced along with the temperature gradient within the drop, and the evaporation behavior conforms to the Kelvin-Clapeyron theory. The algorithm captures the realistic differential thermal gradient in water heated at the surface and is promising for studying various heating/cooling problems, such as thin film evaporation, Leidenfrost effect, and so forth. The simplicity of the algorithm allows it to be easily extended to other surfaces and integrated into various molecular simulation software and user codes.

Boundary-controlled barostats for slab geometries in molecular dynamics simulations

Molecular dynamics simulation barostat schemes are derived for achieving a given normal pressure for a thin liquid or solid layer confined between two parallel walls. This work builds on the boundary-controlled barostat scheme of Lupkowski and van Swol [J. Chem. Phys. 93, 737 (1990)]. Two classes of barostat are explored, one in which the external load is applied to a virtual regular lattice to which the wall atoms are bound using a tethering potential. The other type of barostat applies the external force directly to the wall atoms, which are not tethered. The extent to which the wall separation distribution is Gaussian is shown to be an effective measure of the quality of the barostat. The first class of barostat can suffer from anomalous dynamical signatures, even resonances, which are sensitive to the effective mass of the virtual lattice, whose value lacks any rigorous definition. The second type of barostat performs much better under equilibrium and wall-sliding nonequilibrium conditions and in not being so prone to resonance instabilities in the wall separation and does not require so many largely arbitrary parameters. The results of exploratory simulations which characterize the dynamical response of the model systems for both dry and wet or lubricated systems using the different barostats are presented. The barostats which have an inherent damping mechanism, such as the ones analogous to a damped harmonic oscillator, reduce the occurrence of large fluctuations and resonances in the separation between the two walls, and they also achieve a new target pressure more quickly. Near a nonequilibrium phase boundary the attributes of the barostat can have a marked influence on the observed behavior.

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.