Reversible Surface Energy Storage in Molecular-Scale Porous Materials.

Forcible wetting of hydrophobic pores represents a viable method for energy storage in the form of interfacial energy. The energy used to fill the pores can be recovered as pressure-volume work upon decompression. For efficient recovery, the expulsion pressure should not be significantly lower than the pressure required for infiltration. Hysteresis of the wetting/drying cycle associated with the kinetic barrier to liquid expulsion results in energy dissipation and reduced storage efficiency. In the present work, we use open ensemble (Grand Canonical) Monte Carlo simulations to study the improvement of energy recovery with decreasing diameters of planar pores. Near-complete reversibility is achieved at pore widths barely accommodating a monolayer of the liquid, thus minimizing the area of the liquid/gas interface during the cavitation process. At the same time, these conditions lead to a steep increase in the infiltration pressure required to overcome steric wall/water repulsion in a tight confinement and a considerable reduction in the translational entropy of confined molecules. In principle, similar effects can be expected when increasing the size of the liquid particles without altering the absorbent porosity. While the latter approach is easier to follow in laboratory work, we discuss the advantages of reducing the pore diameter, which reduces the cycling hysteresis while simultaneously improving the stored-energy density in the material.

Solvent-Solvent Correlations across Graphene: The Effect of Image Charges.

Wetting experiments show pure graphene to be weakly hydrophilic, but its contact angle (CA) also reflects the character of the supporting material. Measurements and molecular dynamics simulations on suspended and supported graphene often reveal a CA reduction due to the presence of the supporting substrate. A similar reduction is consistently observed when graphene is wetted from both sides. The effect has been attributed to transparency to molecular interactions across the graphene sheet; however, the possibility of substrate-induced graphene polarization has also been considered. Computer simulations of CA on graphene have so far been determined by ignoring the material's conducting properties. We improve the graphene model by incorporating its conductivity according to the constant applied potential molecular dynamics. Using this method, we compare the wettabilities of suspended graphene and graphene supported by water by measuring the CA of cylindrical water drops on the sheets. The inclusion of graphene conductivity and concomitant polarization effects leads to a lower CA on suspended graphene, but the CA reduction is significantly bigger when the sheets are also wetted from the opposite side. The stronger adhesion is accompanied by a profound change in the correlations among water molecules across the sheet. While partial charges on water molecules interacting across an insulator sheet attract charges of the opposite sign, apparent attraction among like charges is manifested across the conducting graphene. The change is associated with graphene polarization, as the image charges inside the conductor attract equally signed partial charges of water molecules on both sides of the sheet. Additionally, using a nonpolar liquid (diiodomethane), we affirm a detectable wetting translucency when liquid-liquid forces are dominated by dispersive interactions. Our findings are important for predictive modeling toward a variety of applications including sensors, fuel cell membranes, water filtration, and graphene-based electrode materials in high-performance supercapacitors.

Anisotropic structure and dynamics of water under static electric fields.

We study the structure and dynamics of water subject to a range of static external electric fields, using molecular dynamics simulations. In particular, we monitor the changes in hydrogen bond kinetics, reorientation dynamics, and translational motions of water molecules. We find that water molecules translate and rotate slower in electric fields because the tendency to reinstate the aligned orientation reduces the probability of finding a new hydrogen bond partner and hence increases the probability of reforming already ruptured bonds. Furthermore, dipolar alignment of water molecules with the field results in structural and dynamic anisotropies even though the angularly averaged metrics indicate only minor structural changes. Through comparison of selected nonpolarizable and polarizable water models, we find that the electric field effects are stronger in polarizable water models, where field-enhanced dipole moments and thus more stable hydrogen bonds lead to slower switching of hydrogen bond partners and reduced translational mobility, compared to a nonpolarizable water model.

Dynamical insights into the mechanism of a droplet detachment from a fiber.

Quantifying the detachment behavior of a droplet from a fiber is important in many applications such as fog harvesting, oil-water separation, or water management in fuel cells. When the droplets are forcibly removed from hydrophilic fibers, the ease of detachment strongly depends on droplet volume and the rate of the process controlled by the applied force. Experiments, conducted on a ferrofluid under magnetic force, as well as continuum level calculations from fluid mechanics have so far been unable to resolve the time-dependent dynamics of droplet detachment and, most importantly, to assess the role of the applied force as the key determinant of the volume of the droplet residue remaining on the fiber after detachment. In the present work, we study the mechanism of water droplet detachment and retention of residual water on smooth hydrophilic fibers using nonequilibrium molecular dynamics simulations. We investigate how the applied force affects the breakup of a droplet and how the minimal detaching force per unit mass decreases with droplet size. We extract scaling relations that allow extrapolation of our findings to larger length scales that are not directly accessible by molecular models. We find that the volume of the residue on a fiber varies nonmonotonically with the detaching force, reaching the maximal size at an intermediate force and associated detachment time. The strength of this force decreases with the size of the drop, while the maximal residue increases with the droplet volume, V, sub-linearly, in proportion to the V2/3.

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 solvent-averaged 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.

Extent of Surface Force Additivity on Chemically Heterogeneous Substrates at Varied Orientations.

Surface interactions between chemically mixed surfaces, as well as those among dissolved biomolecules, comprise distinct contributions from polar and hydrophobic moieties. These contributions are often context dependent. Approximate compliance to the Cassie additivity equation for the wetting free energies on mixed surfaces in water is, however, indicative of similarly additive forces between individual surface elements, suggesting a quadratic interpolation model for total force from the forces between pure surfaces. We use molecular dynamics/umbrella sampling simulations of parallel and nonparallel mixed surfaces with demonstrable Cassie-like behavior to verify how well the total surface force between the heterogeneous, molecularly rough surfaces can be approximated as a combination of forces among the homogeneous ones. When accounting for dissimilar distances of approach between functional groups of different types, our results for graphene surfaces with mixed methyl and nitrile coating show such a superposition to provide a reasonable first order approximation of interactions between the platelets. Deviations from additivity are more prominent in parallel-plate configurations, at high content of hydrophobic groups, and small separations. The inclusion of water polarizability does not visibly alter the observed behavior regardless of platelet orientations. The outcome of this study determines the necessary molecular conditions for observing force additivity that emphasize the context dependence of hydrophobic interaction in the presence of polar groups. This notion provides guidelines for the syntheses of new, chemically heterogeneous materials with tailored function-oriented properties in aqueous media.

Dynamic Response in Nanoelectrowetting on a Dielectric.

Droplet spreading at an applied voltage underlies the function of tunable optical devices including adjustable lenses and matrix display elements. Faster response and the enhanced resolution motivate research toward miniaturization of these devices to nanoscale dimensions. The response of an aqueous nanodroplet to an applied field can differ significantly from macroscopic predictions. Understanding these differences requires characterization at the molecular level. We describe the equilibrium and nonequilibrium molecular dynamics simulations of nanosized aqueous droplets on a hydrophobic surface with the embedded concentric electrodes. Constant electrode potential is enforced by a rigorous account of the metal polarization. We demonstrate that the reduction of the equilibrium contact angle is commensurate to, and adjusts reversibly with, the voltage change. For a droplet with O(10) nm diameter, a typical response time to the imposition of the field is of O(10(2)) ps. Drop relaxation is about twice as fast when the field is switched off. The friction coefficient obtained from the rate of the drop relaxation on the nonuniform surface, decreases when the droplet approaches equilibrium from either direction, that is, by spreading or receding. The strong dependence of the friction on the surface hydrophilicity points to the dominance of the liquid-surface friction at the drop's perimeter as described in the molecular kinetic theory. This approach enables correct predictions of trends in dynamic responses associated with varied voltage or substrate material.

Universal Repulsive Contribution to the Solvent-Induced Interaction Between Sizable, Curved Hydrophobes.

In addition to the direct attraction, sizable hydrophobes in water experience an attractive force mediated by interfacial water. Using simple geometric arguments, we identify the conditions at which the water-induced interaction between curved hydrocarbon surfaces becomes repulsive. The repulsive contribution arises from the thermodynamic penalty due to the emergence of the liquid/vapor boundary created as water gets expelled between curved hydrophobes. By augmenting the mean field approach with atomistic simulations of pristine and alkyl-coated graphitic nanoparticles in three distinct geometries, spherical, cylindrical and planar, immersed in water, we show the macroscopic thermodynamics remarkably works down to the molecular scale. The new insights improve the prediction and control of wetting and dispersion properties for a broad class of nonpolar nanoparticles.

Electrolyte pore/solution partitioning by expanded grand canonical ensemble Monte Carlo simulation.

Using a newly developed grand canonical Monte Carlo approach based on fractional exchanges of dissolved ions and water molecules, we studied equilibrium partitioning of both components between laterally extended apolar confinements and surrounding electrolyte solution. Accurate calculations of the Hamiltonian and tensorial pressure components at anisotropic conditions in the pore required the development of a novel algorithm for a self-consistent correction of nonelectrostatic cut-off effects. At pore widths above the kinetic threshold to capillary evaporation, the molality of the salt inside the confinement grows in parallel with that of the bulk phase, but presents a nonuniform width-dependence, being depleted at some and elevated at other separations. The presence of the salt enhances the layered structure in the slit and lengthens the range of inter-wall pressure exerted by the metastable liquid. Solvation pressure becomes increasingly repulsive with growing salt molality in the surrounding bath. Depending on the sign of the excess molality in the pore, the wetting free energy of pore walls is either increased or decreased by the presence of the salt. Because of simultaneous rise in the solution surface tension, which increases the free-energy cost of vapor nucleation, the rise in the apparent hydrophobicity of the walls has not been shown to enhance the volatility of the metastable liquid in the pores.

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. Conflicting porosity requirements have so far precluded direct implementations of a two-way control: the pores have to be sufficiently wide to allow water infiltration at experimentally relevant voltages but should not exceed the kinetic threshold for spontaneous expulsion in the absence of the field. Applicable widths are restricted below a few nanometers. Only a narrow window of fields 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, field-controlled transitions between nanoconfined 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 field enhances the uptake of water and ions in the confinement. In systems prone to capillary evaporation, the process can be reversed at sufficient strength of the electric field. For alternating displacement field, we identify the conditions where O (ns) responses of the reversible infiltration/expulsion cycle can be secured for experimentally realizable field strengths, porosity, and salinity of the solution.

Wetting transparency of graphene in water.

Measurements of contact angle on graphene sheets show a notable dependence on the nature of the underlying substrate, a phenomenon termed wetting transparency. Our molecular modeling studies reveal analogous transparency in case of submerged graphene fragments in water. A combined effect of attractive dispersion forces, angle correlations between aqueous dipoles, and repulsion due to the hydrogen-bond-induced orientation bias in polarized hydration layers acting across graphene sheet, enhances apparent adhesion of water to graphene. We show wetting free energy of a fully wetted graphene platelet to be about 8 mNm(-1) lower than for graphene wetted only on one side, which gives close to 10° reduction in contact angle. This difference has potential implications for predictions of water absorption vs. desorption, phase behavior of water in aqueous nanoconfinements, solvent-induced interactions among graphitic nanoparticle and concomitant stability in aqueous dispersions, and can influence permeability of porous materials such as carbon nanotubes by water and aqueous solutions.

Wettability of pristine and alkyl-functionalized graphane.

Graphane is a hydrogenated form of graphene with high bandgap and planar structure insensitive to a broad range of chemical substitutions. We describe an atomistic simulation approach to predict wetting properties of this new material. We determine the contact angle to be 73°. The lower hydrophobicity compared to graphene is explained by the increased planar density of carbon atoms while we demonstrate that the presence of partial charges on carbon and hydrogen atoms plays only a minor role. We further examine the effects of graphane functionalization by alkyl groups of increasing chain lengths. The gradual increase in contact angle with chain length offers a precise control of surface wettability. A saturated contact angle of 114° is reached in butylated form. We find the saturation of contact angle with respect to the length of the functional groups to coincide with the loss of water's ability to penetrate the n-alkyl molecular brush and interact with carbon atoms of the underlying lattice. Since no experimental data have yet become available, our modeling results provide the first estimate of the wettability of graphane. The results also show how its alkyl functionalization provides the basis for a variety of chemical modifications to tune hydrophilicity while preserving the planar geometry of the substrate.

Nanoscale wetting under electric field from molecular simulations.

Applying an electric field is a well-established experimental method to tune surface wettability. As accessible experimental length scales become shorter, the modification of interfacial properties of water using electric field must come to grips with novel effects existing at the nanoscale. We survey recent progress in understanding these effects on water interfacial tension and on water-mediated interactions using molecular simulations. We highlight the key role of external conditions in determining the system's response to applied electric field. We further discuss the role of appropriate boundary conditions in modeling polar fluids subject to collective polarization. The work reviewed here broadens the basic understanding of applied and internal field effects that can operate in condensed phase systems, from modulating local hydrophilicity/hydrophobicity of engineered and biological surfaces, to surface manipulation in nanofluidic devices.

Computational probe of cavitation events in protein systems.

Previous all-atom simulations have identified several classes of proteins where hydrophobic de-wetting (cavitation) is at play. Here we develop and validate a computationally fast method that predicts in which protein systems water spontaneously cavitates. We implement a cubic lattice model, which incorporates the protein shape from crystallographic data and the protein-water interactions from thermodynamic data. Combining it with the previously developed coarse-grained model for water, we determine the extent of occupancy of water at protein-protein interfaces and in protein-ligand cavities. The model captures essential findings from all-atom molecular dynamics studies on the same systems by distinguishing protein cavities that dry from those that remain wet. We also interpret the origin of the cavitation inside the melittin tetramer on simple thermodynamic grounds, and show that part of the mellitin surface is sufficiently hydrophobic to trigger cavitation. Using Glauber/Kawasaki dynamics we obtain the time-scales for de-wetting events that are in agreement with those from all-atom simulations. The method can serve as an intermediate step between the necessary initial screening that identifies proteins with abundance of hydrophobic patches using bioinformatics tools (L. Hua, X. H. Huang, P. Liu, R. H. Zhou and B. J. Berne, J. Phys. Chem. B, 2007, 111, 9069), and computationally extensive studies that need to incorporate molecular details (e.g. single mutation studies of amino acid residues).

Probing surface tension additivity on chemically heterogeneous surfaces by a molecular approach.

Surface free energy of a chemically heterogeneous surface is often treated as an approximately additive quantity through the Cassie equation [Cassie ABD (1948) Discuss Faraday Soc 3:11-16]. However, deviations from additivity are common, and molecular interpretations are still lacking. We use molecular simulations to measure the microscopic analogue of contact angle, Θ(c), of aqueous nanodrops on heterogeneous synthetic and natural surfaces as a function of surface composition. The synthetic surfaces are layers of graphene functionalized with prototypical nonpolar and polar head group: methyl, amino, and nitrile. We demonstrate positive as well as negative deviations from the linear additivity. We show the deviations reflect the uneven exposure of mixture components to the solvent and the linear relation is recovered if fractions of solvent-accessible surface are used as the measure of composition. As the spatial variations in polarity become of larger amplitude, the linear relation can no longer be obtained. Protein surfaces represent such natural patterned surfaces, also characterized by larger patches and roughness. Our calculations reveal strong deviations from linear additivity on a prototypical surface comprising surface fragments of melittin dimer. The deviations reflect the disproportionately strong influence of isolated polar patches, preferential wetting, and changes in the position of the liquid interface above hydrophobic patches. Because solvent-induced contribution to the free energy of surface association grows as cos Θ(c), deviations of cos Θ(c) from the linear relation directly reflect nonadditive adhesive energies of biosurfaces.

The influence of molecular-scale roughness on the surface spreading of an aqueous nanodrop.

We examine the effect of nanoscale roughness on spreading and surface mobility of water nanodroplets. Using molecular dynamics, we consider model surfaces with sub-nanoscale asperities at varied surface coverage and with different distribution patterns. We test materials that are hydrophobic, and those that are hydrophilic in the absence of surface corrugations. Interestingly, on both types of surfaces, the introduction of surface asperities gives rise to a sharp increase in the apparent contact angle. The Cassie-Baxter equation is obeyed approximately on hydrophobic substrates, however, the increase in the contact angle on a hydrophilic surface differs qualitatively from the behavior on macroscopically rough surfaces described by the Wenzel equation. On the hydrophobic substrate, the superhydrophobic state with the maximal contact angle of 180 degrees is reached when the asperity coverage falls below 25%, suggesting that superhydrophobicity can also be achieved by the nanoscale roughness of a macroscopically smooth material. We further examine the effect of surface roughness on droplet mobility on the substrate. The apparent diffusion constant shows a dramatic slow down of the nanodroplet translation even for asperity coverage in the range of 1% for a hydrophilic surface, while droplets on corrugated hydrophobic surfaces retain the ability to flow around the asperities. In contrast, for smooth surfaces we find that the drop mobility on the hydrophilic surface exceeds that on the hydrophobic one.

Water-mediated ordering of nanoparticles in an electric field.

Interfacial polar molecules feature a strongly anisotropic response to applied electric field, favoring dipole orientations parallel to the interface. In water, in particular, this effect combines with generic orientational preferences induced by spatial asymmetry of water hydrogen bonding under confined geometry, which may give rise to a Janus interface. The two effects manifest themselves in considerable dependence of water polarization on both the field direction relative to the interface and the polarity (sign) of the field. Using molecular simulations, we demonstrate strong field-induced orientational forces acting on apolar surfaces through water mediation. At a field strength comparable to electric fields around a DNA polyion, the torques we predict to act on an adjacent nanoparticle are sufficient to overcome thermal fluctuations. These torques can align a particle with surface as small as 1 nm2. The mechanism can support electrically controlled ordering of suspended nanoparticles as a means of tuning their properties and can find application in electro-nanomechanical devices.

Microscopic dynamics of the orientation of a hydrated nanoparticle in an electric field.

We use atomistic simulations to study the orientational dynamics of a nonpolar nanoparticle suspended in water and subject to an electric field. Because of the molecular-level effects we describe, the torque exerted on the nanoparticle exceeds continuum-electrostatics-based estimates by about a factor of 2. The reorientation time of a 16.2 x 16.2 x 3.35 A(3) nanoparticle in a field |E| > 0.015 V/A is an order of magnitude less than the field-free orientational time (approximately 1 ns). Surprisingly, the alignment speed is nearly independent of the nanoparticle size in this regime. These findings are relevant for design of novel nanostructures and sensors and development of nanoengineering methods.

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 conditions, and contraction at others. 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 the 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.

Attractive surface force in the presence of dissolved gas: a molecular approach.

Despite widespread evidence of the influence of dissolved air on hydrophobic interaction, the mechanisms of observed effects are still unknown. Although some experiments indicate that adsorbed gases can modify the structure of water next to hydrophobic surfaces, gas effects on measured forces have been observed only at large surface separations. Gas-specific depletion of water at a hydrophobic surface has been detected but was not reproduced in subsequent measurements. We use computer simulations to study short-ranged hydrophobic attraction in the absence and presence of dissolved gas and monitor gas adsorption at molecular resolution inaccessible in experiments. Although we observe a significant accumulation of dissolved gases at hydrophobic surfaces, even in supersaturated gas solutions surface concentrations remain too low to induce any significant change in the local structure of water and short-range surface forces. We present direct calculations of the hydrophobic force between model hydrocarbon plates at separations between 1.5 and 4 nm. Although stronger, the calculated solvation force has a similar decay rate as deduced from recent surface force apparatus measurements at a somewhat lower contact angle. Within the statistical uncertainty, short-range attraction is not affected by the presence of dissolved nitrogen, even in supersaturated solution with a gas fugacity as high as 30 atm. Comparisons of the adsorption behavior of N2, O2, CO2, and Ar reveal similar features in contrast to the peculiar suppression of water depletion reported for an Ar solution in a neutron reflectivity experiment. Our calculations reveal a notable difference between pathways to the capillary evaporation of pure water and gas-phase nucleation in confined supersaturated gas solutions.