Anomalous ionic transport in tunable angstrom-size water films on silica.

Liquid and ionic transport through nanometric structures is central to many phenomena, ranging from cellular exchanges to water resource management or green energy conversion. While pushing down toward molecular scales progressively unveils novel transport behaviors, reaching ultimate confinement in controlled systems remains challenging and has often involved 2D Van der Waals materials. Here, we propose an alternative route which circumvents demanding nanofabrication steps, partially releases material constraints, and offers continuously tunable molecular confinement. This soft-matter-inspired approach is based on the spontaneous formation of a molecularly thin liquid film onto fully wettable substrates in contact with the vapor phase of the liquid. Using silicon dioxide substrates, water films ranging from angstrom to nanometric thicknesses are formed in this manner, and ionic transport within the film can then be measured. Performing conductance measurements as a function of confinement in these ultimate regimes reveals a one-molecule thick layer of fully hindered transport nearby the silica, above which continuum, bulk-like approaches account for experimental results. Overall, this work paves the way for future investigations of molecular scale nanofluidics and provides insights into ionic transport nearby high surface energy materials such as natural rocks and clays, building concretes, or nanoscale silica membranes used for separation and filtering.

Connection between water's dynamical and structural properties: Insights from ab initio simulations.

SignificanceFirst-principles calculations, which explicitly account for the electronic structure of matter, can shed light on the molecular structure and dynamics of water in its supercooled state. In this work, we use density functional theory, which relies on a functional to describe electronic exchange and correlations, to evaluate which functional best describes the temperature evolution of bulk water transport coefficients. We also assess the validity of the Stokes-Einstein relation for all the functionals in the temperature range studied, and explore the link between structure and dynamics. Based on these results, we show how transport coefficients can be computed from structural descriptors, which require shorter simulation times to converge, and we point toward strategies to develop better functionals.

Enhanced Heat Flow between Charged Nanoparticles and an Aqueous Electrolyte.

Heat transfer through the interface between a metallic nanoparticle and an electrolyte solution has great importance in a number of applications, ranging from nanoparticle-based cancer treatments to nanofluids and solar energy conversion devices. However, the impact of the surface charge and dissolved ions on heat transfer has been scarcely explored so far. In this study, we compute the interface thermal conductance between hydrophilic and hydrophobic charged gold nanoparticles immersed in an electrolyte using equilibrium molecular dynamics simulations. Compared with an uncharged nanoparticle, we report a 3-fold increase of the Kapitza conductance for a nanoparticle surface charge of +320 mC/m2. This enhancement is shown to be approximately independent of the surface wettability, charge spatial distribution, and salt concentration. This allows us to express the Kapitza conductance enhancement in terms of the surface charge density on a master curve. Finally, we interpret the increase of the Kapitza conductance as a combined result of the shift of the water density distribution toward the charged nanoparticle and an accumulation of the counterions around the nanoparticle surface which increase the Coulombic interaction between the liquid and the charged nanoparticle. These considerations help us to apprehend the role of ions in heat transfer close to electrified surfaces.

Complex coupling between surface charge and thermo-osmotic phenomena.

Thermo-osmotic flows, generated at liquid-solid interfaces by thermal gradients, can be used to produce electric currents from waste heat on charged surfaces. The two key parameters controlling the thermo-osmotic current are the surface charge and the interfacial enthalpy excess due to liquid-solid interactions. While it has been shown that the contribution from water to the enthalpy excess can be crucial, how this contribution is affected by surface charge remained to be understood. Here, we start by discussing how thermo-osmotic flows and induced electric currents are related to the interfacial enthalpy excess. We then use molecular dynamics simulations to investigate the impact of surface charge on the interfacial enthalpy excess, for different distributions of the surface charge, and two different wetting conditions. We observe that surface charge has a strong impact on enthalpy excess, and that the dependence of enthalpy excess on surface charge depends largely on its spatial distribution. In contrast, wetting has a very small impact on the charge-enthalpy coupling. We rationalize the results with simple analytical models, and explore their consequences for thermo-osmotic phenomena. Overall, this work provides guidelines to search for systems providing optimal waste heat recovery performance.

Equilibrium molecular dynamics evaluation of the solid-liquid friction coefficient: Role of timescales.

Solid-liquid friction plays a key role in nanofluidic systems. Following the pioneering work of Bocquet and Barrat, who proposed to extract the friction coefficient (FC) from the plateau of the Green-Kubo (GK) integral of the solid-liquid shear force autocorrelation, the so-called plateau problem has been identified when applying the method to finite-sized molecular dynamics simulations, e.g., with a liquid confined between parallel solid walls. A variety of approaches have been developed to overcome this problem. Here, we propose another method that is easy to implement, makes no assumptions about the time dependence of the friction kernel, does not require the hydrodynamic system width as an input, and is applicable to a wide range of interfaces. In this method, the FC is evaluated by fitting the GK integral for the timescale range where it slowly decays with time. The fitting function was derived based on an analytical solution of the hydrodynamics equations [Oga et al., Phys. Rev. Res. 3, L032019 (2021)], assuming that the timescales related to the friction kernel and the bulk viscous dissipation can be separated. By comparing the results with those of other GK-based methods and non-equilibrium molecular dynamics, we show that the FC is extracted with excellent accuracy by the present method, even in wettability regimes where other GK-based methods suffer from the plateau problem. Finally, the method is also applicable to grooved solid walls, where the GK integral displays complex behavior at short times.

The receding contact line cools down during dynamic wetting.

When a contact line (CL)-where a liquid-vapor interface meets a substrate-is put into motion, it is well known that the contact angle differs between advancing and receding CLs. Using non-equilibrium molecular dynamics simulations, we reveal another intriguing distinction between advancing and receding CLs: while temperature increases at an advancing CL-as expected from viscous dissipation, we show that temperature can drop at a receding CL. Detailed quantitative analysis based on the macroscopic energy balance around the dynamic CL showed that the internal energy change of the fluid due to the change of the potential field along the pathline out of the solid-liquid interface induced a remarkable temperature drop around the receding CL, in a manner similar to latent heat upon phase changes. This result provides new insights for modeling the dynamic CL, and the framework for heat transport analysis introduced here can be applied to a wide range of nanofluidic systems.

Microscopic Origins of the Viscosity of a Lennard-Jones Liquid.

Unlike crystalline solids or ideal gases, transport properties remain difficult to describe from a microscopic point of view in liquids, whose dynamics result from complex energetic and entropic contributions at the atomic scale. Two scenarios are generally proposed: one represents the dynamics in a fluid as a series of energy-barrier crossings, leading to Arrhenius-like laws, while the other assumes that atoms rearrange themselves by collisions, as exemplified by the free volume model. To assess the validity of these two views, we computed, using molecular dynamics simulations, the transport properties of the Lennard-Jones fluid and tested to what extent the Arrhenius equation and the free volume model describe the temperature dependence of the viscosity and of the diffusion coefficient at fixed pressure. Although both models reproduce the simulation results over a wide range of pressure and temperature covering the liquid and supercritical states of the Lennard-Jones fluid, we found that the parameters of the free volume model can be estimated directly from local structural parameters, also obtained in the simulations. This consistency of the results gives more credibility to the free volume description of transport properties in liquids.

Electrokinetic sweeping of colloids at a reactive magnesium oxide interface.

Investigating the electrokinetic (EK) response in the vicinity of interfaces has regained interest due to the development of new membrane based processes for energy harvesting or soil depollution. However, the case of reactive interfaces, ubiquitous in these processes, remains scarcely explored. Here we experimentally investigate the EK response of a model interface between an aqueous electrolyte and a bulk MgO crystal surface (100), for different pH. For that purpose, we use a lab-scale non invasive method to monitor the zeta potential of the interface versus time, by confocal fluorescent particle tracking. An unexpected motion of the particles, repelled and then attracted again by the interface is observed. We attributed this motion to the surface reactivity, inducing ion concentration gradients perpendicular to the interface and subsequent diffusiophoresis of the charged particle. Accordingly, we could describe at a semi-quantitative level the particle dynamics by solving numerically the Poisson-Nernst-Planck equations to establish concentration profile in the system and subsequent diffusiophoretic motion. These experiments open the way to the characterization of both the EK response and the reaction rate in the vicinity of reactive interfaces.

Enhanced local viscosity around colloidal nanoparticles probed by equilibrium molecular dynamics simulations.

Nanofluids-dispersions of nanometer-sized particles in a liquid medium-have been proposed for a wide variety of thermal management applications. It is known that a solid-like nanolayer of liquid of typical thicknesses of 0.5-1 nm surrounding the colloidal nanoparticles can act as a thermal bridge between the nanoparticle and the bulk liquid. Yet, its effect on the nanofluid viscosity has not been elucidated so far. In this article, we compute the local viscosity of the nanolayer using equilibrium molecular dynamics based on the Green-Kubo formula. We first assess the validity of the method to predict the viscosity locally. We apply this methodology to the calculation of the local viscosity in the immediate vicinity of a metallic nanoparticle for a wide range of solid-liquid interaction strength, where a nanolayer of thickness 1 nm is observed as a result of the interaction with the nanoparticle. The viscosity of the nanolayer, which is found to be higher than its corresponding bulk value, is directly dependent on the solid-liquid interaction strength. We discuss the origin of this viscosity enhancement and show that the liquid density increment alone cannot explain the values of the viscosity observed. Rather, we suggest that the solid-like structure of the distribution of the liquid atoms in the vicinity of the nanoparticle contributes to the nanolayer viscosity enhancement. Finally, we observe a failure of the Stokes-Einstein relation between viscosity and diffusion close to the wall, depending on the liquid-solid interaction strength, which we rationalize in terms of the hydrodynamic slip.

Liquid-Solid Slip on Charged Walls: The Dramatic Impact of Charge Distribution.

Nanofluidic systems show great promise for applications in energy conversion, where their performance can be enhanced by nanoscale liquid-solid slip. However, efficiency is also controlled by surface charge, which is known to reduce slip. Combining molecular dynamics simulations and analytical developments, we show the dramatic impact of surface charge distribution on the slip-charge coupling. Homogeneously charged graphene exhibits a very favorable slip-charge relation (rationalized with a new theoretical model correcting some weaknesses of the existing ones), leading to giant electrokinetic energy conversion. In contrast, slip is strongly affected on heterogeneously charged surfaces, due to the viscous drag induced by counterions trapped on the surface. In that case slip should depend on the detailed physical chemistry of the interface controlling the fraction of bound ions. Our numerical results and theoretical models provide new fundamental insight into the molecular mechanisms of liquid-solid slip, and practical guidelines for searching new functional interfaces with optimal energy conversion properties, e.g., for blue energy or waste heat harvesting.

Molecular modeling of aqueous electrolytes at interfaces: Effects of long-range dispersion forces and of ionic charge rescaling.

Molecular dynamics simulations of aqueous electrolytes generally rely on empirical force fields, combining dispersion interactions-described by a truncated Lennard-Jones (LJ) potential-and electrostatic interactions-described by a Coulomb potential computed with a long-range solver. Recently, force fields using rescaled ionic charges [electronic continuum correction (ECC)], possibly complemented with rescaling of LJ parameters [ECC rescaled (ECCR)], have shown promising results in bulk, but their performance at interfaces has been less explored. Here, we started by exploring the impact of the LJ potential truncation on the surface tension of a sodium chloride aqueous solution. We show a discrepancy between the numerical predictions for truncated LJ interactions with a large cutoff and for untruncated LJ interactions computed with a long-range solver, which can bias comparison of force field predictions with experiments. Using a long-range solver for LJ interactions, we then show that an ionic charge rescaling factor chosen to correct long-range electrostatic interactions in bulk accurately describes image charge repulsion at the liquid-vapor interface, and the rescaling of LJ parameters in ECCR models-aimed at capturing local ion-ion and ion-water interactions in bulk- describes well the formation of an ionic double layer at the liquid-vapor interface. Overall, these results suggest that the molecular modeling of aqueous electrolytes at interfaces would benefit from using long-range solvers for dispersion forces and from using ECCR models, where the charge rescaling factor should be chosen to correct long-range electrostatic interactions.

Giant Thermoelectric Response of Nanofluidic Systems Driven by Water Excess Enthalpy.

Nanofluidic systems could in principle be used to produce electricity from waste heat, but current theoretical descriptions predict a rather poor performance as compared to thermoelectric solid materials. Here we investigate the thermoelectric response of NaCl and NaI solutions confined between charged walls, using molecular dynamics simulations. We compute a giant thermoelectric response, 2 orders of magnitude larger than the predictions of standard models. We show that water excess enthalpy-neglected in the standard picture-plays a dominant role in combination with the electro-osmotic mobility of the liquid-solid interface. Accordingly, the thermoelectric response can be boosted using surfaces with large hydrodynamic slip. Overall, the heat harvesting performance of the model systems considered here is comparable to that of the best thermoelectric materials, and the fundamental insight provided by molecular dynamics suggests guidelines to further optimize the performance, opening the way to recycle waste heat using nanofluidic devices.

Shear force measurement of the hydrodynamic wall position in molecular dynamics.

Flows in nanofluidic systems are strongly affected by liquid-solid slip, which is quantified by the slip length and by the position where the slip boundary condition applies. Here, we show that the viscosity, slip length, and hydrodynamic wall position (HWP) can be accurately determined from a single molecular dynamics (MD) simulation of a Poiseuille flow, after identifying a relation between the HWP and the wall shear stress in that configuration. From this relation, we deduce that in gravity-driven flows, the HWP identifies with the Gibbs dividing plane of the liquid-vacuum density profile. Simulations of a generic Lennard-Jones liquid confined between parallel frozen walls show that the HWP for a pressure-driven flow is also close to the Gibbs dividing plane (measured at equilibrium), which therefore provides an inexpensive estimate of the HWP, going beyond the common practice of assuming a given position for the hydrodynamic wall. For instance, we show that the HWP depends on the wettability of the surface, an effect usually neglected in MD studies of liquid-solid slip. Overall, the method introduced in this article is simple, fast, and accurate and could be applied to a large variety of systems of interest for nanofluidic applications.

Viscosity and self-diffusion of supercooled and stretched water from molecular dynamics simulations.

Among the numerous anomalies of water, the acceleration of dynamics under pressure is particularly puzzling. Whereas the diffusivity anomaly observed in experiments has been reproduced in several computer studies, the parallel viscosity anomaly has received less attention. Here we simulate viscosity and the self-diffusion coefficient of the TIP4P/2005 water model over a broad temperature and pressure range. We reproduce the experimental behavior and find additional anomalies at negative pressure. The anomalous effect of pressure on dynamic properties becomes more pronounced upon cooling, reaching two orders of magnitude for viscosity at 220 K. We analyze our results with a dynamic extension of a thermodynamic two-state model, an approach which has proved successful in describing experimental data. Water is regarded as a mixture of interconverting species with contrasting dynamic behaviors, one being strong (Arrhenius) and the other fragile (non-Arrhenius). The dynamic parameters of the two-state models are remarkably close between experiment and simulations. The larger pressure range accessible to simulations suggests a modification of the dynamic two-state model, which in turn also improves the agreement with experimental data. Furthermore, our simulations demonstrate the decoupling between viscosity η and self-diffusion coefficient D as a function of temperature T. The Stokes-Einstein relation, which predicts a constant Dη/T, is violated when T is lowered, in connection with the Widom line defined by an equal fraction of the two interconverting species. These results provide a unifying picture of thermodynamics and dynamics in water and call for experiments at negative pressure.

What Controls Thermo-osmosis? Molecular Simulations Show the Critical Role of Interfacial Hydrodynamics.

Thermo-osmotic and related thermophoretic phenomena can be found in many situations from biology to colloid science, but the underlying molecular mechanisms remain largely unexplored. Using molecular dynamics simulations, we measure the thermo-osmosis coefficient by both mechanocaloric and thermo-osmotic routes, for different solid-liquid interfacial energies. The simulations reveal, in particular, the crucial role of nanoscale interfacial hydrodynamics. For nonwetting surfaces, thermo-osmotic transport is largely amplified by hydrodynamic slip at the interface. For wetting surfaces, the position of the hydrodynamic shear plane plays a key role in determining the amplitude and sign of the thermo-osmosis coefficient. Finally, we measure a giant thermo-osmotic response of the water-graphene interface, which we relate to the very low interfacial friction displayed by this system. These results open new perspectives for the design of efficient functional interfaces for, e.g., waste-heat harvesting.

Electro-osmosis at surfactant-laden liquid-gas interfaces: beyond standard models.

Electro-osmosis (EO) is a powerful tool to manipulate liquids in micro and nanofluidic systems. While EO has been studied extensively at liquid-solid interfaces, the case of liquid-vapor interfaces, found e.g. in foam films and bubbles, remains to be explored. Here we perform molecular dynamics (MD) simulations of EO in a film of aqueous electrolyte covered with fluid layers of ionic surfactants and surrounded by gas. Following the experimental procedure, we compute the zeta potential from the EO velocity, defined as the velocity difference between the middle of the liquid film and the surrounding gas. We show that the zeta potential can be smaller or larger than existing predictions depending on the surfactant coverage. We explain the failure of previous descriptions by the fact that surfactants and bound ions move as rigid bodies and do not transmit the electric driving force to the liquid locally. Considering the reciprocal streaming current effect, we then develop an extended model, which can be used to predict the experimental zeta potential of surfactant-laden liquid-gas interfaces.

Decoupling of viscosity and relaxation processes in supercooled water: a molecular dynamics study with the TIP4P/2005f model.

During the last few decades, many experimental and numerical studies have tried to understand the special dynamics of water at low temperatures by measuring structural relaxation times or shear viscosity, but their conclusions strongly depended on the chosen observable and on the range of temperatures considered. Moreover, recent work [J. Chem. Phys. 2013, 138, 12A526] showed that viscosity and relaxation times could decouple at low temperature in a model binary mixture, raising questions on their equivalence to study supercooled water. Here we used molecular dynamics simulations with the promising TIP4P/2005f water force field to investigate the behavior of both the shear viscosity and the relaxation times of water in a large range of temperatures, in order to get a consistent picture of the dynamics of supercooled water. We show that the TIP4P/2005f model reproduces accurately the experimental values of both the viscosity and the diffusion coefficient over a very large range of temperatures. Focusing first on the structural relaxation dynamics, we observe a decoupling between the so-called α- and ß-relaxation times of water at ca. 350 K, suggesting a supercooled-like dynamics over a very large domain of temperatures. By computing shear viscosity over this domain, we compare the accuracy of several phenomenological laws for low temperature dynamics of water to describe both viscosity and α-relaxation time. Unlike what is usually admitted, our tests suggest those quantities are not coupled at low temperatures, and thus should not be considered equivalent. In particular, deviations from the Stokes-Einstein relation appear at lower temperatures for the viscosity than for the α-relaxation time. These results open new perspectives to understand the dynamics of supercooled water and show the performance of the TIP4P/2005f force field to characterize it.

Molecular dynamics analysis of the friction between a water-methanol liquid mixture and a non-polar solid crystal surface.

We performed molecular dynamics analysis of the momentum transfer at the solid-liquid interface for a water-methanol liquid mixture between parallel non-polar solid walls in order to understand the strong decrease of the friction coefficient (FC) induced by the methanol adsorption at the surface observed in our previous work [S. Nakaoka et al., Phys. Rev. E 92, 022402 (2015)]. In particular, we extracted the individual contributions of water and methanol molecules to the total FC and found that the molecular FC for methanol was larger than that for water. We further showed that the reduction of the total solid-liquid FC upon the increase of the methanol molar fraction in the first adsorption layer occurred as a result of a decrease in the molecular number density as well as a decrease in the molecular FCs of both molecules. Analysis of the molecular orientation revealed that the decrease of the molecular FC of methanol resulted from changes of the contact feature onto the solid surface. Specifically, methanol molecules near the solid surface had their C-O bond parallel to the surface with both CH3 and O sites contacting the solid at low methanol molar fraction, while they had their C-O bond outward from the surface with only the CH3 site contacting the solid at higher methanol molar fraction. The mechanisms discussed in this work could be used to search for alternative water additives to further reduce the solid-liquid friction.

Assessment of elastic models in supercooled water: A molecular dynamics study with the TIP4P/2005f force field.

Glass formers exhibit a viscoelastic behavior: at the laboratory time scale, they behave like (glassy) solids at low temperatures and like liquids at high temperatures. Based on this observation, elastic models relate the long time supercooled dynamics to short time elastic properties of the supercooled liquid. In the present work, we assess the validity of elastic models for the shear viscosity and the α-relaxation time of supercooled water, using molecular dynamics simulations with the TIP4P/2005f force field over a wide range of temperatures. We show that elastic models provide a good description of supercooled water dynamics. For the viscosity, two different regimes are observed and the crossover temperature is found to be close to the one where the Stokes-Einstein relation starts to be violated. Our simulations show that only shear properties are important to characterize the effective flow activation energy. This study calls for experimental determination of the high frequency elastic properties of water at low temperatures.