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.

Effect of an external electric field on the dynamics and intramolecular structures of ions in an ionic liquid.

Simulations of the ionic liquid 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide in an external electric field have been performed using a Drude particle polarizable force field. The structure of the ions has been analyzed, with close attention paid to the configurations of the ions. The "charge arm" concept is used to explain certain changes of these configurations that would be difficult to rationalize otherwise, e.g., trans → cis isomerization of the bis(trifluoromethylsulfonyl)imide anion and extension of the alkyl chain of the cation. It has also been shown that the ions orient themselves so that their charge arms align with and stretch out along the field, and these effects occur at lower external electric field strengths than cause a change in the inherent diffusion of the ions. The dynamics of the system parallel and perpendicular to the field were analyzed, and it was found that the applied field affected the diffusion normal to the field. This is explained as a secondary effect of a change in the ion cage formed by the surrounding counterions of a given ion in the ionic liquid. The breakdown of the ion cages was rationalized by correlating changes to the inherent diffusion of the ions with other changes to the diffusion and bulk structure of the liquid, as well as considering the average forces on the ions compared to the force the ions would be expected to experience in an electric field. Parallel to the field, a drift was observed at every electric field studied. In electric fields with no changes to the ion cage structure, the relationship between the drift and electric field was found to be linear, becoming nonlinear as the ion cage structure breaks down.

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.

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

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.

Structure and lifetimes in ionic liquids and their mixtures.

With the aid of molecular dynamics simulations, we study the structure and dynamics of different ionic liquid systems, with focus on hydrogen bond, ion pair and ion cage formation. To do so, we report radial distribution functions, their number integrals, and various time-correlation functions, from which we extract well-defined lifetimes by means of the reactive flux formalism. We explore the influence of polarizable force fields vs. non-polarizable ones with downscaled charges (±0.8) for the example of 1-butyl-3-methylimidazolium bromide. Furthermore, we use 1-butyl-3-methylimidazolium trifluoromethanesulfonate to investigate the impact of temperature and mixing with water as well as with the chloride ionic liquid. Smaller coordination numbers, larger distances, and tremendously accelerated dynamics are observed when the polarizable force field is applied. The same trends are found with increasing temperature. Adding water decreases the ion-ion coordination numbers whereas the water-ion and water-water coordination is enhanced. A domain analysis reveals that the nonpolar parts of the ions are dispersed and when more water is added the water clusters increase in size. The dynamics accelerate in general upon addition of water. In the ionic liquid mixture, the coordination number around the cation changes between the two anions, but the number integrals of the cation around the anions remain constant and the dynamics slow down with increasing content of the chloride ionic liquid.

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.

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 of water 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.

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.