Electrokinetic flow of 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. Different pressures, electrolyte concentrations and pore sizes are examined. Our results indicate a flow that deviates considerably from the predictions of Poiseuille fluid mechanics. Due to preferential adsorption of the different ionic species by surface SiO! or SiOH groups, we find that a significant electric current is generated, but with opposite polarities using polarizable vs. fixed charge models for water and ions, emphasizing the need for careful parameterization in such complex systems. We also examine the influence of partial deprotonation of the silica surface, and we find that much more current is generated in a dehydrogenated nanopore, even though the overall efficiency remains low. These findings indicate that different methods of nanopore preparation, which can produce a range of surface properties, should be examined more closely in the related experimental methods to generate electrokinetic current

Field exposed water in a nanopore: liquid or vapour?

We study the behavior of ambient temperature water under the combined effects of nanoscale confinement and applied electric field. Using molecular simulations we analyze the thermodynamic causes of field-induced expansion at some, and contraction at other conditions. Repulsion among parallel water dipoles and mild weakening of interactions between partially aligned water molecules prove sufficient to destabilize the aqueous liquid phase in isobaric systems in which all water molecules are permanently exposed to a uniform electric field. At the same time, simulations reveal comparatively weak field-induced perturbations of water structure upheld by flexible hydrogen bonding. In open systems with fixed chemical potential, these perturbations do not suffice to offset attraction of water into the field; additional water is typically driven from unperturbed bulk phase to the field-exposed region. In contrast to recent theoretical predictions in the literature, our analysis and simulations confirm that classical electrostriction characterizes usual electrowetting behavior in nanoscale channels and nanoporous materials.Comment: 20 pages, 6 figures + T.O.C. figure, in press in PCC

Effect of Field Direction on Electrowetting in a Nanopore

We manifest a significant influence of field direction and polarity on surface wetting, when the latter is tuned by application of an external electric field. Thermodynamics of field-induced filling of hydrocarbon-like nanopores with water is studied by open ensemble molecular simulation. Increased field strength consistently results in water-filling and electrostriction in hydrophobic nanopores. A threshold field commensurate with surface charge density of about one elementary charge per 10 nm2 suffices to render prototypical paraffin surfaces hydrophilic. When a field is applied in the direction perpendicular to the confining walls, the competition between orientational polarization and angle preferences of interfacial water molecules relative to the walls results in an asymmetric wettability of opposing surfaces (Janus interface). Reduction of surface free energy observed upon alignment of confinement walls with field direction suggests a novel mechanism whereby the applied electric field can operate selectively on water-filled nanotubes while empty ones remain unaffected

Curvature dependence of the effect of ionic functionalization on the attraction among nanoparticles in dispersion

Solubilization of nanoparticles facilitates nanomaterial processing and enables new applications. An effective method to improve dispersibility in water is provided by ionic functionalization.We explore how the necessary extent of functionalization depends on the particle geometry. Using molecular dynamics/umbrella sampling simulations, we determine the effect of the solute curvature on solventaveraged interactions among ionizing graphitic nanoparticles in aqueous dispersion. We tune the hydrophilicity of molecular-brush coated fullerenes, carbon nanotubes, and graphane platelets by gradually replacing a fraction of the methyl end groups of the alkyl coating by the ionizing –COOK or –NH3Cl groups. To assess the change in nanoparticles’ dispersibility in water, we determine the potential-of-mean-force profiles at varied degrees of ionization. When the coating comprises only propyl groups, the attraction between the hydrophobic particles intensifies from spherical to cylindrical to planar geometry. This is explained by the increasing fraction of surface groups that can be brought into contact and the reduced access to water molecules, both following the above sequence. When ionic groups are added, however, the dispersibility increases in the opposite order, with the biggest effect in the planar geometry and the smallest in the spherical geometry. These results highlight the important role of geometry in nanoparticle solubilization by ionic functionalities, with about twice higher threshold surface charge necessary to stabilize a dispersion of spherical than planar particles. At 25%–50% ionization, the potential of mean force reaches a plateau because of the counterion condensation and saturated brush hydration. Moreover, the increase in the fraction of ionic groups can weaken the repulsion through counterion correlations between adjacent nanoparticles. High degrees of ionization and concomitant ionic screening gradually reduce the differences among surface interactions in distinct geometries until an essentially curvature-independent dispersion environment is created. Insights into tuning nanoparticle interactions can guide the synthesis of a broad class of nonpolar nanoparticles, where solubility is achieved by ionic functionalization

Nanoconfined water under electric field at constant chemical potential undergoes electrostriction

Electric control of nanopore permeation by water and solutions enables gating in membrane ion channels and can be exploited for transient surface tuning of rugged substrates, to regulate capillary permeability in nanofluidics, and to facilitate energy absorption in porous hydrophobic media. Studies of capillary effects, enhanced by miniaturization, present experimental challenges in the nanoscale regime thus making molecular simulations an important complement to direct measurement. In a molecular dynamics (MD) simulation, exchange ofwater between the pores and environment requires modeling of coexisting confined and bulk phases, with confined water under the field maintaining equilibrium with the unperturbed environment. In the present article, we discuss viable methodologies for MD sampling in the above class of systems, subject to size-constraints and uncertainties of the barostat function under confinement and nonuniform-field effects. Smooth electric field variation is shown to avoid the inconsistencies of MD integration under abruptly varied field and related ambiguities of conventional barostatting in a strongly nonuniform interfacial system. When using a proper representation of the field at the border region of the confined water, we demonstrate a consistent increase in electrostriction as a function of the field strength inside the pore open to a field-free aqueous environment

Interplay of local hydrogen-bonding and long-ranged dipolar forces in simulations of confined water

Spherical truncations of Coulomb interactions in standard models for water permit efficient molecular simulations and can give remarkably accurate results for the structure of the uniform liquid. However truncations are known to produce significant errors in nonuniform systems, particularly for electrostatic properties. Local molecular field (LMF) theory corrects such truncations by use of an effective or restructured electrostatic potential that accounts for effects of the remaining long-ranged interactions through a density-weighted mean field average and satisfies a modified Poisson's equation defined with a Gaussian-smoothed charge density. We apply LMF theory to three simple molecular systems that exhibit different aspects of the failure of a naive application of spherical truncations -- water confined between hydrophobic walls, water confined between atomically-corrugated hydrophilic walls, and water confined between hydrophobic walls with an applied electric field. Spherical truncations of 1/r fail spectacularly for the final system in particular, and LMF theory corrects the failings for all three. Further, LMF theory provides a more intuitive way to understand the balance between local hydrogen bonding and longer-ranged electrostatics in molecular simulations involving water.Comment: Submitted to PNA

Persistence length of a polyelectrolyte in salty water: a Monte-Carlo study

We address the long standing problem of the dependence of the electrostatic persistence length $l_e$ of a flexible polyelectrolyte (PE) on the screening length $r_s$ of the solution within the linear Debye-Huckel theory. The standard Odijk, Skolnick and Fixman (OSF) theory suggests $l_e \propto r_s^2$, while some variational theories and computer simulations suggest $l_e \propto r_s$. In this paper, we use Monte-Carlo simulations to study the conformation of a simple polyelectrolyte. Using four times longer PEs than in previous simulations and refined methods for the treatment of the simulation data, we show that the results are consistent with the OSF dependence $l_e \propto r_s^2$. The linear charge density of the PE which enters in the coefficient of this dependence is properly renormalized to take into account local fluctuations.Comment: 7 pages, 6 figures. Various corrections in text and reference

Dynamic Control of Nanopore Wetting in Water and Saline Solutions under an Electric Field

Field-induced nanopore wetting by aqueous solutions, including electrolytes, provides opportunities for a variety of applications. Con!icting porosity requirements have so far precluded direct implementations of a two-way control: the pores have to be su ciently wide to allow water in#ltration at experimentally relevant voltages but should not exceed the kinetic threshold for spontaneous expulsion in the absence of the #eld. Applicable widths are restricted below a few nanometers. Only a narrow window of #elds and pore geometries can simultaneously satisfy both of the above requirements. Accurate accounts of wetting equilibria and dynamics at nanoscale porosity require molecular level descriptions. Here we use molecular dynamics simulations to study dynamic, #eld-controlled transitions between nanocon#ned liquid and vapor phases in contact with an unperturbed aqueous or electrolyte environment. In nanopores wetted by electrolyte solutions, we observe depletion of salt compared to the bulk phase. The application of a local electric #eld enhances the uptake of water and ions in the con#nement. In systems prone to capillary evaporation, the process can be reversed at su cient strength of the electric #eld. For alternating displacement #eld, we identify the conditions where O (ns) responses of the reversible in#ltration/ expulsion cycle can be secured for experimentally realizable #eld strengths, porosity, and salinity of the solution

Reversible electrowetting transitions on superhydrophobic surfaces

Electric field applied across the interface has been shown to enable transitions from Cassie to Wenzel state on superhydrophobic surfaces with miniature corrugations. Molecular Dynamics (MD) simulations manifest the possibility of reversible cycling between the two states when narrow surface wells support spontaneous expulsion of water in the absence of the field. With approximately 1 nm sized wells between the surface asperities, response times to changes of electric field are of O(0.1) ns, allowing up to GHz frequency of the cycle. Because of orientation preferences of interfacial water in contact with the solid, the phenomenon depends on the polarity of the field normal to the interface. The threshold field strength for the Cassie-to-Wenzel transition is significantly lower for the field pointing from the aqueous phase to the surface, however, once in the Wenzel state, the opposite field direction secures tighter filling of the wells. Considerable hysteresis revealed by the delayed water retraction at decreasing field strength indicates the presence of moderate kinetic barriers to expulsion. Known to scale approximately with the square of the length scale of the corrugations, these barriers preclude the use of increased corrugation sizes while the reduction of the well diameter necessitates stronger electric fields. Field-controlled Cassie-to-Wenzel transitions are therefore optimized by using superhydrophobic surfaces with nanosized corrugations. Abrupt changes indicate a high degree of cooperativity reflecting the correlations between wetting states of interconnected wells on the textured surface

Salt and Water Uptake in Nanocon!nement under Applied Electric Field: An Open Ensemble Monte Carlo Study

Permeation of electrolytes in nanoporous materials underlies many applications in energy and materials technologies. Wetting of apolar nanopores can be enhanced by electric !eld, attracting water and ions from unperturbed electrolyte bath. We study absorption of water and NaCl in the pores by Expanded Ensemble Grand Canonical Monte Carlo simulation, which implements particle insertions and deletions through incremental changes in particles’ coupling with the system. We determine the uptake of water and ions in the pores, and concomitant changes in pore thermodynamics, as functions of !eld strength in the pore and salinity in the external bath. Pressure increase and reduction of wetting free energy, !, in the pore intensify near-quadratically with the !eld. Surprisingly, the in uence of bulk salinity on ! can change qualitatively with pore width and !eld strength. Conforming to Gibbs adsorption isotherm, narrow pores with salt molality below that of the bath experience an increase in ! with rising bulk salinity. The !eld can change salt depletion to excess and consequently reverse the salinity dependence of wetting free energy from increasing to declining function of bulk molality. Field polarity continues to play a role, leading to asymmetric wettability at opposing walls as we previously observed in the absence of ions