Non-Continuum Intercalated Water Diffusion Explains Fast Permeation

Oct 25, 2017 - The non-continuum collective diffusion of water intercalation between graphene layers facilitates fast molecular transport due to reduc...
0 downloads 0 Views 3MB Size
www.acsnano.org

Non-Continuum Intercalated Water Diffusion Explains Fast Permeation through Graphene Oxide Membranes Shuping Jiao and Zhiping Xu* Applied Mechanics Laboratory, Department of Engineering Mechanics, and Center for Nano and Micro Mechanics, Tsinghua University, Beijing 100084, China S Supporting Information *

ABSTRACT: Recent experimental studies have revealed unconventional phase and transport behaviors of water confined within lamellar graphene oxide membranes, which hold great promise not only in improving our current understanding of nanoconfined water but also in developing high-performance filtration and separation applications. In this work, we explore molecular structures and diffusive dynamics of water intercalated between graphene or graphene oxide sheets. We identify the monolayer structured water between graphene sheets at temperature T below Tc = ∼315 K and an interlayer distance d = 0.65 nm, which is absent as the sheets are oxidized. The non-continuum collective diffusion of water intercalation between graphene layers facilitates fast molecular transport due to reduced wall friction. This solid-like structural order of intercalated water is disturbed as T or d increases to a critical value, with abnormal declines in the coefficients of collective diffusion. Based on a patched model of graphene oxide sheets consisting of spatially distributed pristine and oxidized regions, we conclude that the non-continuum collective diffusion of intercalated water can explain fast water permeation through graphene oxide membranes as reported in recent experimental studies, in stark contrast to the conventional picture of pressure-driven continuum flow with boundary slip, which has been widely adopted in literature but may apply only at high humidity or in the fully hydrated conditions. KEYWORDS: graphene, graphene oxide, intercalated water, collective diffusion, permeation

T

contrast with those of bulk water and are controlled by a number of parameters such as pressure, temperature, interlayer spacings, and the interfacial commensurability between the water lattices and solid walls.10,11 The solid nature of intercalated water implicates interesting kinetic and dynamical behaviors that differ from the conventional picture of molecular liquids, as demonstrated in recent studies for water on solid surfaces or within nanopores.18−23 Ultrahigh permeation rates of water vapor through GO membranes were reported recently, with an enhancement factor of ∼107 compared to the non-slip viscous-flow prediction.24 This unexpected enhancement was attributed to the slippery nature of water flow between pristine graphene regions in the GO sheets, which offer smooth potential energy surfaces, highly reduced wall friction, and the capillary pressure acting as an additional driven force for the flow. However, either the form of water intercalation or the microstructure of GO membranes has not been resolved with direct experimental

he presence of water intercalation in nanoconfined channels and slits has received much attention because of its relevant applications in filtration, separation, and energy storage. Nanofluidic applications based on lamellar membranes have been proposed for highly selective transport of gas, liquid molecules, and ions as well as efficient heat exchange.1−7 With molecular structures distinctly different from that in the bulk,8−10 water intercalated within lamellae such as graphene oxide (GO) membranes is expected to exhibit unique diffusive, flow behaviors and unconventional mechanisms for selective mass transport. Recent experimental and theoretical studies have demonstrated the existence of two-dimensional (2D) ordered water structures while intercalated within graphitic, MoS2 layers or between graphene and mica/MoS2 substrates, as a result of the nanoconfinement and/or pressure.11−16 Although the nature of “ice” formation under nanoconfinement is still not clear from experimental evidence, and distinction between ice and metastable water needs further clarification,12,17 computational studies have demonstrated the existence of ordered water structures as intercalated between paralleled hydrophobic plates, in forms of square, hexagon, pentagon, or mixed lattices.8−10 More interestingly, the phase behaviors and transitions of nanoconfined water shows high © 2017 American Chemical Society

Received: July 31, 2017 Accepted: October 25, 2017 Published: October 25, 2017 11152

DOI: 10.1021/acsnano.7b05419 ACS Nano 2017, 11, 11152−11161

Article

Cite This: ACS Nano 2017, 11, 11152-11161

Article

ACS Nano

Figure 1. (a) Molecular structures and the hydrogen-bonding network of intercalated water confined between two graphene sheets. (b) 2D structure factors of intercalated water between graphene, calculated at T = 300, 310, 315 K at d = 0.65 nm and NW = 986.

evidence. At high humidity or a fully hydrated condition, the interlayer gallery of the membrane may be filled up by water, and thus a pressure-driven flow mechanism can be manifested in a percolated network. However, for vapor- or liquid-phase water transport in a partially hydrated membrane, the continuum flow picture may fail to capture the underlying physics. In this situation, as the interlayer distance in GO membranes is much shorter than the mean free paths of gas molecules at ambient conditions, the hydrodynamic pressure between water intercalations cannot be transmitted, and their collective diffusion instead of advection may contribute dominantly to the mass transport through the membrane as a percolated flow network cannot be attained. However, this diffusion-controlled mechanism of water transport through GO membranes has not been discussed yet. In this work, we explore the molecular structures and collective diffusion of finite-size water clusters intercalated between two parallel graphene or GO layers, by performing molecular dynamics (MD) simulations. The results are discussed to elucidate the underlying mechanism of efficient

mass transport through GO membranes. We calculate the coefficients of collective diffusion and wall friction, which demonstrate abnormal dependence on temperature and the interlayer distance, conforming to the undergoing solid−liquid structural transition. The contribution of collective water diffusion to mass transport through the GO membrane is estimated, which explains the fast water permeation measured in experiments, offering an alternative picture to the continuum flow with boundary slip.

RESULTS AND DISCUSSION Molecular Structures and Solid−Liquid Phase Transition of Intercalated Water. To model the intercalated water, we place a droplet of NW water molecules between two graphene sheets at a specific interlayer distance d. The top and bottom layers are constrained in plane. Experimental evidence shows that the d in GO membranes increases with the relative humidity (RH), due to the formation of water capillary in hydrophilic regions with oxygen-rich functional groups.4,24 Monolayer structured water with a 2D quasi-square lattice is 11153

DOI: 10.1021/acsnano.7b05419 ACS Nano 2017, 11, 11152−11161

Article

ACS Nano

Figure 2. Dependence of collective diffusion coefficients D of intercalated water on (a) the number of water molecules NW at d = 0.65 nm and T = 300 K and (b) the interlayer distance d between graphene sheets with NW = 221−486 and T = 300 K. (c) The dependence of D on NW at d = 0.65 nm, d = 0.9 nm, and T = 300 K and (d) on d for GO sheets with NW = 986 and T = 300 K. Here d for GO is defined as the distance between two basal carbon sheets.

Figure 3. (a) The dependence of collective diffusion coefficients D of intercalated water on the temperature T at d = 0.65 nm and NW = 986. (b) Friction coefficients calculated for the intercalated water at d = 0.65 nm and NW = 986, which demonstrate distinct dependence on T; λ is friction coefficient, and A is area of the water intercalation.

identified at d = 0.65 nm and T = 300 K (Figure 1a), which aligns with recent simulation and experimental results.9,12 This value of d corresponds to a local energy minimum of graphene layers containing the intercalation with NW = 58−986 water molecules (Figure S1). However, this regular 2D lattice structure no longer exists as T or d increases, which can be quantified from the changes in the structure factors (Figures 1b and S2). Specifically, the order-to-disorder structural transition is first-order,25 occurring as T increases to ∼310−325 K for intercalated water with d = 0.65 nm and NW = 58−986. The transition temperature Tc slightly decreases with NW. At T = 300 K the structural transition also occurs by tuning d, at dc = 0.77 nm for NW = 221, 486, and 986 (Figure S2). These results indicate that the ordered 2D quasi-square structure is sensitive

to temperature and the geometric confinement. Our MD simulation results show that as the interlayer distance d increases beyond that of the bilayer, the distinct layering order is lost. Comparative simulations are carried out for water intercalated between GO layers, where no ordered lattice structure forms because of the surface roughness due to the presence of hydroxyl functionalization as well as the hydrogenbond (H-bond) interaction between the hydroxyl groups and water molecules that perturbs the H-bond network within intercalation (Figure 1a). Diffusive Behaviors of Intercalated Water. We measure the coefficients of collective diffusion D for intercalated water droplets from the MD-simulated trajectories, by using the Einstein’s relation. The results show distinct dependence of D 11154

DOI: 10.1021/acsnano.7b05419 ACS Nano 2017, 11, 11152−11161

Article

ACS Nano

Figure 4. (a) Non-Gaussian parameter α2 calculated for d = 0.65 nm and NW = 986 at temperature close to Tc, the transition temperature. (b) Distribution of discretized molecular displacements within a time step of td = 30 ps, dx and dy are the in-plane displacements along the x and y directions. For T < Tc = 315 K, subplots show heavy-tailed distributions that deviate significantly from Gaussian diffusion, while for T ≥ Tc, the Gaussian behavior is recovered. (c) The distribution and (d) time evolution of the angle between the water lattice and graphene lattice (see the inset in panel c) at d = 0.65 nm and NW = 986.

on the size of droplet at d = 0.65 nm and T = 300 K, which decreases rapidly with NW (Figure 2a). With NW changes from 58 to 986, D decreases from 30.87 × 10−5 cm2/s to 1.49 × 10−5 cm2/s (Figure 2a). This size effect was previously reported for water droplets diffusing on the top of graphene,18 which indicates that the collective diffusion is driven by thermal fluctuation, and the value of D would then scale inversely with the mass. For comparison, the molecular self-diffusion coefficient of water in the bulk phase is D0 = 2.5 × 10−5 cm2/s at ambient conditions.18 On the other hand, increasing d from 0.65 to 0.78 nm enhances D by 1−3 times, to D = ∼10−4 cm2/s for NW = 221−486 (Figure 2b), resulting from the weakened confinement. These results suggest that the intercalated water could diffuse in a collective way that is much faster than the self-diffusion of water molecules in the bulk phase, which may lead to rapid water transport in graphitic channels. This enhanced transport could be attributed to the fact that the collective diffusion of intercalated water, which is absent in the bulk, is critically modulated by the friction at the solid−liquid interface, while the pristine graphene provides an atomistically smooth surface with ultralow friction that facilitates fast collective mass transport over it. The coefficient

of self-diffusion in the intercalated water, however, is much lower (by two orders) than the value in the bulk within the temperature range where the solid-like structure of intercalated water is preserved. To gain more insight into the nature of collective diffusion, we then explore the T-dependence of D for water intercalation with NW = 986. As shown in Figure 3a, D increases with T in general, but exhibits a sudden drop in its value at T = 305−315 K, which corresponds to the critical temperature for the orderto-disorder transition in intercalated water (Figure 1b). This finding reveals that the collective diffusion of intercalated water is closely associated with its molecular structures, and the change in phases of water leads to a shift in the diffusion mechanisms. To further illustrate this fact, we carry out MD simulations for droplet sizes with NW = 221 and 486 (Figure S3). The simulation results demonstrate abnormal reduction of D as well at T = 320−325 and 315−320 K, respectively, which is consistent with the critical temperature Tc predicted for the solid−liquid structural transition (Figure S4). Moreover, we find that the statistical error in D near Tc (Tc = 300−315 K for NW = 986) is significant, and the values of D along the x (armchair) and y (zigzag) directions (Dxx, Dyy) are anisotropic 11155

DOI: 10.1021/acsnano.7b05419 ACS Nano 2017, 11, 11152−11161

Article

ACS Nano

Figure 5. (a, b) The length l of collective diffusive paths in lamellar GO membranes, estimated by eq 1, for water intercalated between (a) pristine (d = 0.65 nm, T = 300 K) and (b) oxidized (d = 0.9 nm, T = 300 K) graphene sheets. (c, d) Illustration of the diffusive transport processes for water vapor permeation through GO membranes. The arrows highlight the shortest paths of transport. Note that for illustration, a single water molecule is plotted instead of the intercalated clusters, and only one side of the graphene sheet is oxidized in GO, with hydroxyl groups only.

below Tc, but isotropic above Tc. This contrast can also be seen from the distribution of diffusive distance ld (extracted for each 30 ps interval), which demonstrates this anisotropy−isotropy transition at Tc (Figure 4b). This is a direct consequence from the structural change with T, where water intercalation features a crystalline lattice below Tc. The relative orientation of ice lattice aligned to the graphene sheet shows discrete jumps of π/ 3 along the path of diffusion, which matches the lattice symmetry of graphene (Figures 4c,d). Along with these findings, the T-dependence of D can also be attributed to the interfacial friction between intercalated water and the graphene wall. We calculate the interfacial friction coefficient λ for water intercalation with NW = 986 and d = 0.65 nm and find a rapid increase in λ near Tc (T = 310−315 K) that indicates a solid− solid interface is less resistant to the collective diffusion than the solid−liquid interface (Figure 3b). This rapid jump of λ conforms to the abnormal decline of D as we have seen in Figure 3a, along with the structural transition. The order−disorder structural transition could also be triggered by tuning the interlayer distance d. From the MD simulation results, we find that D for the collective diffusion

increases with d in general as the nanoconfinement enforced by the graphene or GO sheets is weakened, and there is a sudden drop at d = ∼0.77 nm where the monolayer solid-like order is lost, as shown in Figure 2b specifically for NW = 221 and 486. This critical value of d aligns with the solid−liquid transition as indicated by the structure factors of intercalated water (Figure S2). These results further confirm our conclusion made on the correlation between collective diffusion of intercalated water and its molecular structures. As there is no in-plane structural order in the intercalated water between oxidized graphene sheets, the efficient mass transport by collective diffusion could be much reduced. We consider oxidized graphene sheets functionalized with hydroxyl groups (O:C = 20%), where the interlayer distance is d = 0.65 and 0.9 nm, respectively. Our MD simulation results show that within the narrow channel (d = 0.65 nm), no cluster of the intercalated water is formed, and thus there is no collective motion. While at a larger interlayer distance d = 0.9 nm, the collective diffusion of intercalated water is observed, with D decreases with NW and increases with d (Figure 2c,d). These results indicate that the intercalated water in hydrated GO 11156

DOI: 10.1021/acsnano.7b05419 ACS Nano 2017, 11, 11152−11161

Article

ACS Nano

structural characteristics and permeability based on our simulation results for the diffusivity. By assuming that the water intercalation is randomly distributed in the interlayer gallery, diffusing through the membrane through the shortest path in the complicate network (Figure 5), the characteristic distance for an intercalation to diffusive across one layer of the stacks is l = 0.25L. In this analysis, we neglect the cross-layer barrier for the collective diffusion of intercalated water droplets, which is reasonable for nanopores or nanoslits with size larger than or comparable with the droplets, or for GO sheets with large lateral span and low density of nanopores, where lateral diffusion dominates the mass transport process. One can then extract the value of l from the permeability Pr measured from experiments, that is

membranes could diffuse 3−4 orders lower than that between graphene sheets at the same interlayer distance and 2−3 orders lower than that for molecular diffusion in bulk water. It should be noticed here that the structure of water intercalation between oxidized graphene sheets is highly sensitive to the spatial distribution of chemical groups, and its effect awaits further investigation. In this work, to exclude the effects of static and thermal corrugations in the highly defective GO sheets, we consider both full and partial planar constraints in the MD simulations, which shows that this effect can be neglected. With the remarkable dependence of D on the molecular structures of intercalated water identified from MD simulations, we then probe the dynamical behaviors of collective diffusion. We track the center-of-mass trajectory of the intercalated droplet, with NW = 986, for example, and find a distinct change in the pattern of diffusive motion near Tc (Figure S5). The distributions of stepwise displacements dx and dy have longer tails at temperature below Tc than that above Tc (Figure 4b), demonstrating deviation from the normal Gaussian behavior before the structural phase transition takes place, which is recovered as T increases beyond Tc. The deviation from Gaussian behavior can be captured by the non-Gaussian parameter α2(t) = 3⟨|r(t) − r(0)|4⟩/(5⟨|r(t) − r(0)|2⟩2) − 1.26,27 Our MD simulation results suggest that the parameter α2 reaches a peak and then decreases at a characteristic time scale of td = ∼20 ps for T > Tc (liquid-like) and ∼100 ps for T < Tc (solid-like) (Figure 4a). Consequently, the heavy-tail nature of diffusion observed here could also be responsible for the faster diffusion of ordered intercalation compared to the normal diffusion, and our total simulation time for the calculation of D, 10 ns, is sufficient to extract the coefficient of collective diffusion using the Einstein’s relation. Mass Transport through Lamellar GO Membranes. Ultrafast mass transport of water vapor through GO membranes was reported recently, and the measured water evaporation rate with the presence of a GO film is almost the same as that in the absence of the film.24 This much-enhanced mass transport was attributed to the highly slippery nature of capillary force driven viscous flow within the percolated 2D channels. However, even considering the finite slip length of water flow between graphene sheets (∼60 nm), the factor of flow enhancement estimated theoretically for the slip flow is ∼102, which is much lower than the factor reported from experiments.20 Moreover, the formation of a continuous percolation network across the whole GO membrane has not yet been validated by any experimental evidence, which could be difficult considering the irregular patterns of oxidized and pristine graphene regions in the GO sheets and the complicate brick-and-motar microstructures, especially at low humidity.28 An alternative view of the advective mass transport process is that water transported may exist in forms of 2D intercalated droplets, which could diffuse in both the oxidized and pristine graphene regions but with contrastive diffusion constants. Recently, an experimental study of nanoscale graphitic channels shows that liquid transport inside the channels is the ratelimiting process, compared to the diffusion and evaporation outside.29 We could thus assume the rates for water vapor to condense into or evaporate from the membrane are high enough so that the diffusive mass transport within the membrane is the rate-limiting process. Specifically, for diffusion between graphene sheets with an average lateral size of L, quantitative relation could be derived between the micro-

(mW /A)/(l 2/4D) = Pr or l 2 = 4DmW /(APr )

(1)

where mW/A is the areal mass density of water intercalation in the interlayer space available for the collective diffusion. The area of the intercalation we studied in this work is ∼10−100 nm2, which is compared to the values reported from experimental studies.28,30−32 Our MD simulation results predict that the value of D is in the range from 1.49 to 30.87 × 10−5 cm2/s for water intercalation with NW = 986−58 and d = 0.65 nm at 300 K (Figure 2a). The measured vapor permeation through GO membranes shows that the values of d and Pr increase with the RH, which reach ∼1.1 nm, 0.48 mgh−1mm−2 at 100% RH, respectively.24 The permeation rate is Pr = ∼0.024 mgh−1mm−2 at the RH corresponding to the condition d = 0.65 nm.24 Although at high humidity, continuous percolated networks of water in the GO membrane could form and the pressure-driven flow may dominate the mass transport process,30 collective diffusion is expected to control at low RH conditions,28 where intercalated water cannot fill the whole transport pathway in the GO membrane. As there are no experimental measurements for the amount of water intercalated in GO, we consider a range of areal number density N/A of water molecules from 1 μm−2 to 103 μm−2 in estimating l, which corresponds to the mass density mW/Ad = 0.088−88 g/m3. The values of l estimated from eq 1 for specific size of the intercalation NW are summarized in Figure 5a. The results show that l decreases with NW and decreases with mW/A for a certain permeation rate Pr measured from experiments. Specifically, the length of diffusing path l ranges from a few nanometers to ∼0.7 μm for NW = 58− 986, and for the diffusion of a single water molecule (NW = 1), l is ∼0.17−5.4 μm. These estimated values of l are in a reasonable range considering the typical GO sheet size in the membrane is ∼1 μm,24 and the smaller value of l than L could be attributed to the presence of trans-sheet nanopores that cut the diffusive paths within the 2D interlayer gallery and the irregular shape of GO sheet. It should be remarked here that this estimation neglects the concentration gradient along the thickness direction, with the assumption of diffusion-limited regime of mass transport across the GO membrane. Intercalated water could also diffuse over the hydrophilic regions of GO sheets that are functionalized with oxygen-rich groups. We extend our estimation of l by considering the collective diffusion of intercalated water droplets within this pathway to the overall mass transport across GO membranes. With the experimentally measured value Pr = 0.024 mgh−1mm−2 at d = 0.65 nm, and using the values of D calculated from our MD simulations (Figure 2c), we estimate that l = ∼0.062−4.55 nm for NW = 986−58, N/A = 1−1000 11157

DOI: 10.1021/acsnano.7b05419 ACS Nano 2017, 11, 11152−11161

Article

ACS Nano μm−2; l could be interpreted here as the averaged interpore distance in the oxidized region considering its highly defective nature.33,34 We also make estimations for the measured permeability Pr = 0.21 mgh−1mm−2 at d = 0.9 nm, which yield l = ∼ 0.060−7.97 nm for NW = 986−58, N/A = 1−1000 μm−2 (Figure 5b). The largest value of l estimated here is comparable to the value of ∼1−10 nm characterized from GO samples, although one should note that these values highly depend on the preparing processes of GO membranes.35−38 This agreement leads to the conclusion that both fast diffusion in the interlayer gallery between pristine graphene regions and slow diffusion a cross the GO sheets contribute to the high permeability rate of water vapor measured in experiments.24 Consequently, our results can be used to explain the ultrahigh rate of water mass transport, through an alternative mechanism manifested by efficient collective diffusion of intercalated water droplets through percolated graphene channels, or nanopores in the oxidized region, as illustrated in Figure 5c,d. This picture differs from the slip-flow view that has been widely adopted in recent studies.20,24,39 A clarification of the underlying molecular processes calls for more detailed experimental characterization of the membrane microstructures as well as the distribution of defects and intercalated water content and size. More detailed investigation of RH-dependent mass transport could also offer some clues that may identify a crossover from the diffusion-dominated mechanism to slip flow. Additional Remarks. It should be noted that our discussion on diffusion-controlled mass transport in GO membranes is limited by excluding a few factors that are related to their complex microstructures. First, both pristine and oxidized regions in the GO sheets have a limited areal span. Consequently, in a pristine channel with size comparable to the intercalated droplet, the collective diffusion could be retarded by the side-pinning effect resulting from the capillary force, and the droplet could even attract the droplet into the oxidized region.39 Second, wide channels could form in the wrinkled regions of GO sheets, and the nanodroplets would then need to diffuse over the surface of GO sheets.18 Third, the layer-bylayer microstructure of GO membranes indicates that the crosslayer transport through nanopores or nanoslits may become the rate-limiting process due to the presence of a cross-layer energy barrier, especially in the situation of very thin membranes with high densities of nanopores or nanoslits.40,41 Although the intercalated droplets could diffuse collectively over pores or slits with size larger than or comparable with them, evaporation− condensation process may be critical for diffusion through small pores. Conclusive arguments on the consideration on these complexities requires further clarification. The understanding of collective diffusion explored in this work could also shed light on the selective gas and ion transport processes through GO membranes that are not fully wetted with water continuum inside. The water intercalation could enhance the contrast in permeability of gas molecules by blocking non-solvable molecules, while carrying solvable ones, although the diffusive coefficient of solvable gases such as CO2 is reduced once immersed in water.42 Ions that are immobile in dry GO membranes could be encapsulated in the droplet and diffuse with the water shell.43,44 The mechanism of collective diffusion identified for the intercalated water and the calculated coefficients of diffusion enables quantitative estimation of the contributions for mass transport through both the pristine and oxidized regions, provided with their areal fraction and spatial patterns identified for the GO membranes.

CONCLUSION In this work, we explored diffusion of intercalated water between pristine and oxidized graphene sheets in order to understand water vapor permeation through GO membranes. Our MD simulation results show that the diffusive behaviors of droplets are defined by the structural ordering in it, which demonstrate anomalous deviation from normal Gaussian diffusion and reduced interfacial friction below the critical temperature for solid−fluid structural transition. The calculated coefficients of diffusion are implemented in a patchy structural model of GO sheets to estimate the characteristics of diffusive pathway within the lamellar GO membrane, which explains the high permeation rate identified in recent experimental studies, indicating that the non-continuum diffusive processes play a dominant role in mass transport at low humidity and offers an alternative picture to the continuum slip flow of water across GO membranes at high humidity or in a fully hydrated condition. MODELS AND METHODS Molecular Structures. Models with NW = 58−986 water molecules are constructed, which are intercalated between two 20 × 20 nm graphene or GO sheets (Figure 1a). The number of water molecules, or the lateral span of water droplets considered here, covers the typical range of pristine/oxidized region sizes in the graphene oxide sheets.33,34 Consequently, the structural and dynamical properties of water droplets investigated here could capture the characteristic mass transport processes of water intercalation in graphene oxide membranes. An open boundary is used in the z direction, and periodic boundary conditions (PBCs) are applied in the in-plane x and y directions. The proposed molecular models of GO sheets in the literature consist of hydroxyl and epoxy groups on the basal plane and carbonyl groups at defective sites and open edges.45 For the surface functional groups on graphene, hydroxyl groups are reported to be able to stay enriched in the long-living quasi-equilibrium state,45 and a typical fraction of hydroxyl species relative to the amount of carbon atoms in GO is ∼20%.46 Considering these experimental evidence, we constructed in this work hydroxyl-functionalized graphene with c = nOH/nC = 20%. Here nOH and nC are the numbers of hydroxyl groups and carbon atoms, respectively. The distribution of hydroxyl groups was sampled randomly on both sides of the sheet. Molecule Dynamics Simulations. We perform MD simulations using the large-scale atomic/molecular massively parallel simulator (LAMMPS). 47 The all-atom optimized potentials for liquid simulations (OPLS-AA)48 are used for graphene and GO sheets. The SPC/E model49 of water is used in our study, which is widely adopted for MD simulations of water transport as it predicts reasonable density, diffusivity, and viscosity compared to other models.8,49−52 The effect of water models on the intercalated water phases was discussed in our previous work.25 The conclusions in this work are cross-checked by comparing with results simulated using the TIP4P model. We find that, although the critical temperature and interlayer distance for structural order transition are different in these two different models, the characteristics of observations agree well, demonstrating the robustness of the conclusion we made. Long-range Coulomb interactions are computed by using the particle−particle particle-mesh algorithm (PPPM).53 The interaction between carbon atoms in graphene and oxygen atoms in water is modeled with parameters εC−O = 4.063 meV and σC−O = 0.319 nm, which predict a water contact angle (WCA) of θc,G = 98.4° for graphene, consistent with the value measured experimentally.20,54 The WCA θc,GO for GO is lower than θc,G and decreases with the concentration c of oxygen-rich functional groups. For a typical value of c = 20% for GO, our simulation result is θc,GO = 26.8°, which is also close to recent experimental reports.39 The time step to integrate the equation-ofmotion is 1 fs, with the SHAKE algorithm applied for the stretching terms between oxygen and hydrogen atoms of water to avoid the need 11158

DOI: 10.1021/acsnano.7b05419 ACS Nano 2017, 11, 11152−11161

Article

ACS Nano

Calculation of the Interfacial Friction Coefficients. As a parameter that quantifies energy dissipation at the water/wall interface, the friction coefficient λ can be calculated via the Green−Kubo (GK) relation through the autocorrelation function (ACF), ⟨F(t)F(0)⟩, of fluctuating pairwise forces at equilibrium:

of integration with very short time steps for hydrogen-atom-related high-frequency vibrations. Product simulations are carried out after 1 ns equilibration at specific temperature using the Nosé−Hoover thermostat with a damping time constant of 100 fs. In the simulations, we consider both full and partial (by fixing a few carbon atoms randomly) planar constraints of carbon atoms in the GO sheets that do not lead to notable difference for the conclusions we draw in this work (Figure S6). Key parameters used in our model are listed in Table 1.55

λ=

models

sites

σ (nm)

q (e)

graphene

C

0.06831

0.33997

0

GO

C

0.06831

0.33997

+0.1966

O

0.15535

0.3166

−0.5260

H O H

0 0.15535 0

0 0.3166 0

+0.3294 −0.8472 +0.4236

water

bond length (nm)

ref

lC−C = 0.142 lC−O = 0.141 lO−H = 0.0945

48

lO−H = 0.1

49

48

⟨F(t )F(0)⟩dt

(3)

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsnano.7b05419. Figures S1−S7 offer more details on the simulation results and analysis (PDF)

The L-J parameters for atomic pairs are evaluated using the LorentzBerthelot mixing rules.

AUTHOR INFORMATION

Structure Factor Analysis. To assess the structural order of intercalated water between graphene and GO layers, the in-plane structure factor is calculated as8,56 S(q) =



ASSOCIATED CONTENT

a

⎡ N ⎞2 ⎤ ⎞2 ⎛ NW 1 ⎢⎛⎜ W ⎟⎥ ⎜ ⎟ cos( q · r ) + sin( q · r ) ∑ ∑ i ⎟ i ⎟ ⎥ ⎜ NW ⎢⎣⎜⎝ i = 1 ⎠⎦ ⎝ i=1 ⎠

∫0

where F(t) is the total forces exerted on the water intercalation from its interaction with the graphene sheets, kB is the Boltzmann constant, A is the contact area between intercalated water and graphene or GO sheets, and T is the temperature. We calculate the ACF of F(t) in our 5 ns-long equilibrium MD simulations. It should be pointed out that there is a well-documented difficulty to obtain the GK relation via equilibrium MD simulations due to the finite-size effect, which often leads to vanishing of friction coefficients for a very long time simulation.57,58 As a result, the integration of ACF in eq 3 should be done within a reasonable cutoff time tc to resolve this issue. A widely adopted recipe is to use the time corresponding to the onset of a plateau in the integration as tc.57,58 In our analysis, tc is estimated to be in the range of 1−2 ps.

Table 1. 12-6 Lennard-Jones (L-J) Potential Parameters and Atomic Charges Used in the MD Models of Graphene and GOa ε (kcal/mol)

1 AkΒT

Corresponding Author

*E-mail: [email protected]. ORCID

Zhiping Xu: 0000-0002-2833-1966 (2)

Author Contributions

Z. X. designed the project. Both authors performed the study and wrote the manuscript.

Here rj is position of the j-th oxygen atom in the plane of intercalation, NW is the number of water molecules within the intercalation, q is a planar wave vector, and ⟨···⟩ denotes the time average in thermal equilibrium. Analysis of the Diffusive Behaviors. To elucidate the nature of diffusion, we analyze the distribution of diffusion distance ld in a constant interval of simulation time td. We consider the value of td ranging from 20 to 100 ps. We find that the distribution profile of ld or the stepwise displacements dx and dy depend on the choice of td, but the function that can be used to well fit the data is not affected for td < 100 ps, and thus we choose td = 30 ps for our discussion. Characterization of Molecular and Collective Diffusion. In the Gaussian regime, the molecular diffusion constant D is calculated from the trajectories of water molecules by using the Einstein’s definition relating the correlation function of atomic positions ri, or the mean-square distance (MSD), to the value of D = limt→∞⟨|r(t) − r(0)|2⟩/2dit, where di is the dimension of space one explores, t is the simulation time, and ⟨···⟩ is the ensemble average. The collective motion of intercalated water in the 2D interlayer gallery (di = 2) is calculated from its center-of-mass trajectories. In our MD simulations with a time span of a few nanoseconds, the MSD ⟨|r(t) - r(0)|2⟩ of water molecules is calculated based on the time-series of atomic (oxygen) positions. Time averaging and statistical error analysis are performed for 5 independent MD runs, with 3000 or more time-series in each run that are broken down into segments of 100 ps duration that starts at different simulation time. Molecular structures of the whole system are equilibrated using the Nosé−Hoover thermostat at 290−350 K for 10 ns, to ensure the convergence of calculated values of D as a function of the simulation time (Figure S7). We extract MDsimulated trajectories for the calculation of diffusion constant after the first 1 ns record for equilibration.

Notes

The authors declare no competing financial interest.

ACKNOWLEDGMENTS This work was supported by the National Natural Science Foundation of China through grant no. 11472150, the National Science & Technology Major Project (2016ZX05011003), the Tsinghua University Initiative Scientific Research Program through grant no. 2014z22074, and the State Key Laboratory of Mechanics and Control of Mechanical Structures (Nanjing University of Aeronautics and Astronautics) through grant no. MCMS-0414G01. The computation was performed on the Explorer 100 cluster system at Tsinghua National Laboratory for Information Science and Technology. REFERENCES (1) Kim, H. W.; Yoon, H. W.; Yoon, S.-M.; Yoo, B. M.; Ahn, B. K.; Cho, Y. H.; Shin, H. J.; Yang, H.; Paik, U.; Kwon, S.; Choi, J.-Y.; Park, H. B. Selective Gas Transport through Few-Layered Graphene and Graphene Oxide Membranes. Science 2013, 342, 91−95. (2) Li, H.; Song, Z.; Zhang, X.; Huang, Y.; Li, S.; Mao, Y.; Ploehn, H. J.; Bao, Y.; Yu, M. Ultrathin, Molecular-Sieving Graphene Oxide Membranes for Selective Hydrogen Separation. Science 2013, 342, 95− 98. (3) Joshi, R.; Carbone, P.; Wang, F.-C.; Kravets, V. G.; Su, Y.; Grigorieva, I. V.; Wu, H.; Geim, A. K.; Nair, R. R. Precise and Ultrafast 11159

DOI: 10.1021/acsnano.7b05419 ACS Nano 2017, 11, 11152−11161

Article

ACS Nano Molecular Sieving through Graphene Oxide Membranes. Science 2014, 343, 752−754. (4) Mi, B. Graphene Oxide Membranes for Ionic and Molecular Sieving. Science 2014, 343, 740−742. (5) Yumura, T.; Yamasaki, A. Roles of Water Molecules in Trapping Carbon Dioxide Molecules inside the Interlayer Space of Graphene Oxides. Phys. Chem. Chem. Phys. 2014, 16, 9656−9666. (6) Olson, E. J.; Ma, R.; Sun, T.; Ebrish, M. A.; Haratipour, N.; Min, K.; Aluru, N. R.; Koester, S. J. Capacitive Sensing of Intercalated H2O Molecules Using Graphene. ACS Appl. Mater. Interfaces 2015, 7, 25804−25812. (7) Tong, W. L.; Ong, W.-J.; Chai, S.-P.; Tan, M. K.; Huang, Y. M. Enhanced Evaporation Strength through Fast Water Permeation in Graphene-Oxide Deposition. Sci. Rep. 2015, 5, 11896. (8) Giovambattista, N.; Rossky, P. J.; Debenedetti, P. G. Phase Transitions Induced by Nanoconfinement in Liquid Water. Phys. Rev. Lett. 2009, 102, 050603. (9) Zhao, W.-H.; Wang, L.; Bai, J.; Yuan, L.-F.; Yang, J.; Zeng, X. C. Highly Confined Water: Two-Dimensional Ice, Amorphous Ice, and Clathrate Hydrates. Acc. Chem. Res. 2014, 47, 2505−2513. (10) Zhao, W.-H.; Bai, J.; Yuan, L.-F.; Yang, J.; Zeng, X. C. Ferroelectric Hexagonal and Rhombic Monolayer Ice Phases. Chem. Sci. 2014, 5, 1757−1764. (11) Chen, J.; Schusteritsch, G.; Pickard, C. J.; Salzmann, C. G.; Michaelides, A. Two Dimensional Ice from First Principles: Structures and Phase Transitions. Phys. Rev. Lett. 2016, 116, 025501. (12) Algara-Siller, G.; Lehtinen, O.; Wang, F.; Nair, R.; Kaiser, U.; Wu, H.; Geim, A.; Grigorieva, I. Square Ice in Graphene Nanocapillaries. Nature 2015, 519, 443−445. (13) Deng, M.; Kwac, K.; Li, M.; Jung, Y.; Park, H. G. Stability, Molecular Sieving, and Ion Diffusion Selectivity of a Lamellar Membrane from Two-Dimensional Molybdenum Disulfide. Nano Lett. 2017, 17, 2342−2348. (14) Kwac, K.; Kim, I.; Pascal, T. A.; Goddard, W. A.; Park, H. G.; Jung, Y. Multilayer Two-Dimensional Water Structure Confined in MoS2. J. Phys. Chem. C 2017, 121, 16021−16028. (15) Bampoulis, P.; Siekman, M. H.; Kooij, E. S.; Lohse, D.; Zandvliet, H. J. W.; Poelsema, B. Latent Heat Induced Rotation Limited Aggregation in 2D Ice Nanocrystals. J. Chem. Phys. 2015, 143, 034702. (16) Bampoulis, P.; Teernstra, V. J.; Lohse, D.; Zandvliet, H. J. W.; Poelsema, B. Hydrophobic Ice Confined between Graphene and MoS2. J. Phys. Chem. C 2016, 120, 27079−27084. (17) Debenedetti, P. G. Metastable Liquids; Princeton University Press: Princeton, NJ, 1996. (18) Ma, M.; Tocci, G.; Michaelides, A.; Aeppli, G. Fast Diffusion of Water Nanodroplets on Graphene. Nat. Mater. 2016, 15, 66−71. (19) Diallo, S.; Mamontov, E.; Nobuo, W.; Inagaki, S.; Fukushima, Y. Enhanced Translational Diffusion of Confined Water under Electric Field. Phys. Rev. E 2012, 86, 021506. (20) Wei, N.; Peng, X.; Xu, Z. Breakdown of Fast Water Transport in Graphene Oxides. Phys. Rev. E 2014, 89, 012113. (21) Park, J. H.; Aluru, N. Ordering-Induced Fast Diffusion of Nanoscale Water Film on Graphene. J. Phys. Chem. C 2010, 114, 2595−2599. (22) Smit, B.; Maesen, T. L. Molecular Simulations of Zeolites: Adsorption, Diffusion, and Shape Selectivity. Chem. Rev. 2008, 108, 4125−4184. (23) Carrasco, J.; Hodgson, A.; Michaelides, A. A Molecular Perspective of Water at Metal Interfaces. Nat. Mater. 2012, 11, 667−674. (24) Nair, R. R.; Wu, H. A.; Jayaram, P. N.; Grigorieva, I. V.; Geim, A. K. Unimpeded Permeation of Water through Helium-Leak-Tight Graphene-Based Membranes. Science 2012, 335, 442−444. (25) Jiao, S.; Duan, C.; Xu, Z. Structures and Thermodynamics of Water Encapsulated by Graphene. Sci. Rep. 2017, 7, 2646. (26) Rahman, A. Correlations in the Motion of Atoms in Liquid Argon. Phys. Rev. 1964, 136, A405−A411.

(27) Kob, W.; Donati, C.; Plimpton, S. J.; Poole, P. H.; Glotzer, S. C. Dynamical Heterogeneities in a Supercooled Lennard-Jones Liquid. Phys. Rev. Lett. 1997, 79, 2827−2830. (28) Rezania, B.; Severin, N.; Talyzin, A. V.; Rabe, J. P. Hydration of Bilayered Graphene Oxide. Nano Lett. 2014, 14, 3993−3998. (29) Radha, B.; Esfandiar, A.; Wang, F.; Rooney, A.; Gopinadhan, K.; Keerthi, A.; Mishchenko, A.; Janardanan, A.; Blake, P.; Fumagalli, L.; et al. Molecular Transport through Capillaries Made with AtomicScale Precision. Nature 2016, 538, 222−225. (30) Lerf, A.; Buchsteiner, A.; Pieper, J.; Schöttl, S.; Dekany, I.; Szabo, T.; Boehm, H. Hydration Behavior and Dynamics of Water Molecules in Graphite Oxide. J. Phys. Chem. Solids 2006, 67, 1106− 1110. (31) Li, Q.; Song, J.; Besenbacher, F.; Dong, M. Two-Dimensional Material Confined Water. Acc. Chem. Res. 2015, 48, 119−127. (32) Medhekar, N. V.; Ramasubramaniam, A.; Ruoff, R. S.; Shenoy, V. B. Hydrogen Bond Networks in Graphene Oxide Composite Paper: Structure and Mechanical Properties. ACS Nano 2010, 4, 2300−2306. (33) Gao, W.; Alemany, L. B.; Ci, L.; Ajayan, P. M. New Insights into the Structure and Reduction of Graphite Oxide. Nat. Chem. 2009, 1, 403−408. (34) Bagri, A.; Mattevi, C.; Acik, M.; Chabal, Y. J.; Chhowalla, M.; Shenoy, V. B. Structural Evolution During the Reduction of Chemically Derived Graphene Oxide. Nat. Chem. 2010, 2, 581−587. (35) Loh, K. P.; Bao, Q.; Eda, G.; Chhowalla, M. Graphene Oxide as a Chemically Tunable Platform for Optical Applications. Nat. Chem. 2010, 2, 1015−1024. (36) Dikin, D. A.; Stankovich, S.; Zimney, E. J.; Piner, R. D.; Dommett, G. H. B.; Evmenenko, G.; Nguyen, S. T.; Ruoff, R. S. Preparation and Characterization of Graphene Oxide Paper. Nature 2007, 448, 457−460. (37) Erickson, K.; Erni, R.; Lee, Z.; Alem, N.; Gannett, W.; Zettl, A. Determination of the Local Chemical Structure of Graphene Oxide and Reduced Graphene Oxide. Adv. Mater. 2010, 22, 4467−4472. (38) Gómez-Navarro, C.; Meyer, J. C.; Sundaram, R. S.; Chuvilin, A.; Kurasch, S.; Burghard, M.; Kern, K.; Kaiser, U. Atomic Structure of Reduced Graphene Oxide. Nano Lett. 2010, 10, 1144−1148. (39) Huang, H.; Song, Z.; Wei, N.; Shi, L.; Mao, Y.; Ying, Y.; Sun, L.; Xu, Z.; Peng, X. Ultrafast Viscous Water Flow through NanostrandChannelled Graphene Oxide Membranes. Nat. Commun. 2013, 4, 2979. (40) Willcox, J. A.; Kim, H. J. Molecular Dynamics Study of Water Flow across Multiple Layers of Pristine, Oxidized, and Mixed Regions of Graphene Oxide. ACS Nano 2017, 11, 2187−2193. (41) Muscatello, J.; Jaeger, F.; Matar, O. K.; Muller, E. A. Optimizing Water Transport through Graphene-Based Membranes: Insights from Nonequilibrium Molecular Dynamics. ACS Appl. Mater. Interfaces 2016, 8, 12330−12336. (42) Jiao, S.; Xu, Z. Selective Gas Diffusion in Graphene Oxides Membranes: A Molecular Dynamics Simulations Study. ACS Appl. Mater. Interfaces 2015, 7, 9052−9059. (43) Sun, P.; Zhu, M.; Wang, K.; Zhong, M.; Wei, J.; Wu, D.; Xu, Z.; Zhu, H. Selective Ion Penetration of Graphene Oxide Membranes. ACS Nano 2013, 7, 428−437. (44) Sun, P.; Zheng, F.; Zhu, M.; Song, Z.; Wang, K.; Zhong, M.; Wu, D.; Little, R. B.; Xu, Z.; Zhu, H. Selective Trans-Membrane Transport of Alkali and Alkaline Earth Cations through Graphene Oxide Membranes Based on Cation-pi Interactions. ACS Nano 2014, 8, 850−859. (45) Kim, S.; Zhou, S.; Hu, Y.; Acik, M.; Chabal, Y. J.; Berger, C.; de Heer, W.; Bongiorno, A.; Riedo, E. Room-Temperature Metastability of Multilayer Graphene Oxide Films. Nat. Mater. 2012, 11, 544−549. (46) Dreyer, D. R.; Park, S.; Bielawski, C. W.; Ruoff, R. S. The Chemistry of Graphene Oxide. Chem. Soc. Rev. 2010, 39, 228−240. (47) Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1−19. (48) Jorgensen, W. L.; Maxwell, D. S.; Tirado-Rives, J. Development and Testing of the OPLS All-Atom Force Field on Conformational 11160

DOI: 10.1021/acsnano.7b05419 ACS Nano 2017, 11, 11152−11161

Article

ACS Nano Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 1996, 118, 11225−11236. (49) Berendsen, H. J. C.; Grigera, J. R.; Straatsman, T. P. The Missing Term in Effective Pair Potentials. J. Phys. Chem. 1987, 91, 6269−6271. (50) Svishchev, I. M.; Kusalik, P. G.; Wang, J.; Boyd, R. J. Polarizable Point-Charge Model for Water: Results under Normal and Extreme Conditions. J. Chem. Phys. 1996, 105, 4742−4750. (51) Baez, L. A.; Clancy, P. Existence of a Density Maximum in Extended Simple Point Charge Water. J. Chem. Phys. 1994, 101, 9837−9840. (52) Wu, Y.; Tepper, H. L.; Voth, G. A. Flexible Simple Point-Charge Water Model with Improved Liquid-State Properties. J. Chem. Phys. 2006, 124, 024503. (53) Hockney, R. W.; Eastwood, J. W. Computer Simulation Using Particles; Taylor & Francis: New York, 1989. (54) Rafiee, J.; Mi, X.; Gullapalli, H.; Thomas, A. V.; Yavari, F.; Shi, Y.; Ajayan, P. M.; Koratkar, N. A. Wetting Transparency of Graphene. Nat. Mater. 2012, 11, 217−222. (55) Shih, C.-J.; Lin, S.; Sharma, R.; Strano, M. S.; Blankschtein, D. Understanding the pH-Dependent Behavior of Graphene Oxide Aqueous Solutions: A Comparative Experimental and Molecular Dynamics Simulation Study. Langmuir 2012, 28, 235−241. (56) Falk, K.; Sedlmeier, F.; Joly, L.; Netz, R. R.; Bocquet, L. Molecular Origin of Fast Water Transport in Carbon Nanotube Membranes: Superlubricity Versus Curvature Dependent Friction. Nano Lett. 2010, 10, 4067−4073. (57) Xiong, W.; Liu, J. Z.; Ma, M.; Xu, Z.; Sheridan, J.; Zheng, Q. Strain Engineering Water Transport in Graphene Nanochannels. Phys. Rev. E 2011, 84, 056329. (58) Bocquet, L.; Barrat, J.-L. Hydrodynamic Boundary Conditions, Correlation Functions, and Kubo Relations for Confined Fluids. Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top. 1994, 49, 3079−3092.

11161

DOI: 10.1021/acsnano.7b05419 ACS Nano 2017, 11, 11152−11161