Dynamic Wetting of Nanodroplets on Smooth and Patterned Graphene

5 days ago - Wettability of graphene-coated surface has gained significant attention due to the practical applications of graphene. In terms of static...
1 downloads 6 Views 1MB Size
Subscriber access provided by UNIV OF SCIENCES PHILADELPHIA

C: Surfaces, Interfaces, Porous Materials, and Catalysis

Dynamic Wetting of Nanodroplets on Smooth and Patterned Graphene-Coated Surface Shih-Wei Hung, and Junichiro Shiomi J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b01344 • Publication Date (Web): 02 Apr 2018 Downloaded from http://pubs.acs.org on April 2, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Dynamic Wetting of Nanodroplets on Smooth and Patterned Graphene-Coated Surface Shih-Wei Hung, Junichiro Shiomi,*

Department of Mechanical Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-8656, Japan

AUTHOR INFORMATION

Corresponding Author: [email protected]

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT Wettability of graphene-coated surface has gained significant attention due to the practical applications of graphene. In terms of static contact angle, the wettability of graphene-coated surfaces has been widely investigated by experiments, simulations, and theories in recent years. However, the studies of dynamic wetting on graphene-coated surfaces are limited. In the present study, molecular dynamics simulation is performed to understand the dynamic wetting of water nanodroplet on graphene-coated surfaces with different underlying substrates, copper, and graphite, from a nanoscopic point of view. The results show that the spreading on smooth graphene-coated surfaces can be characterized by the final equilibrium radius of droplet and inertial time constant. The dynamics of spreading of water nanodroplet on graphene-coated surface with patterned nanostructures on graphene is also carried out. Unlike the smooth graphene-coated surfaces, dynamic wetting on the patterned graphene-coated surfaces depends on the underlying substrate: the inhibition of wetting by the surface nanostructures is stronger for a hydrophobic underlying substrate than the hydrophilic substrate.

2

ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

INTRODUCTION Vast research efforts have targeted on graphene owing to its unique electronic, optical, chemical and mechanical properties.1–4 The advanced coating technique using graphene also benefits from its extraordinary properties, e.g. excellent chemical resistance, impermeability to gases, anti-bacterial properties, and mechanical strength.5 The practical application of graphene coatings requires further understanding of how the single-atom-thick sheet affects the surface properties of the underlying substrate, such as wetting. Recently, there have been extensive researches6–19 in quantifying the wetting behavior of graphene-coated surfaces in terms of static contact angle. An early work of Rafiee et al.11 found the “wetting transparency” of graphene, i.e. the wetting behavior of underlying substrate, where van der Waals forces dominate the wetting, is not significantly altered by the addition of graphene. However, the follow-up experimental studies demonstrate the graphene is not completely but partially transparent to wetting.9,10,12,14,15 Computer simulations6,7,12,13,15 and theoretical calculations6,8,12,15 also support the partial wetting transparency of graphene. Those researches indicate the monolayer graphene acts like a barrier, reducing the water-underlying substrate interaction. The progress in quantifying the static wetting on graphene-coated surface has been substantial, but the reported study of dynamic wetting on graphene-coated surface is 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

still limited. The dynamic wetting is a multiscale and complex problem as the classical hydrodynamic description breaks down at the moving contact line, and the nanoscopic features at the molecular scale needs to be considered.20 Thus, the dynamic wetting shows different behaviors at different length scales due to the size effect. A recent study of capillary filling in nanochannels21 demonstrates that the dynamic contact angle at the nanoscale differs significantly from the static contact angle. While molecular kinetic theory22 has been formulated to understand the advance of contact line in terms of kinetics of molecules, the most direct method to investigate the details of the dynamics is molecular dynamics (MD) simulation, which have been widely used to study dynamic wetting phenomena, particularly at the nanoscale.23–33 The results of MD simulations on flat surfaces are generally in agreement with the molecular kinetic theory.24–26,33 Although there is a limitation in the diameter of the simulated droplet, typically in the order of tens of nanometers (nanodroplet), MD simulations can trace the evolution of the droplet spreading, from which scaling of contact radius R can be directly obtained. When comparing with the macroscopic experiments of initial rapid wetting on flat surface that are typically studied with millimeter droplets, they exhibit some similarities in the scaling despite the large difference in the length scale; R evolves with time, t, as R~t1/2 for completely wetting surfaces23,34,35 in analogy to the

4

ACS Paragon Plus Environment

Page 4 of 28

Page 5 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

scaling of an inertia-dominated coalescence,36,37 and the exponent decreases with decreasing wettability for partial wetting surfaces.30,31,38,39 The dominant physics that counteracts capillarity in the initial rapid wetting phase is expected to be mostly inertia and, for droplets at the microscale, the time histories of R for different droplet sizes coincide with each other by using the inertial scaling even for partial wetting surfaces.38 However, the inertial scaling does not capture the dependence on the wettability. As for the nanodroplets in MD simulations, no scaling for neither diameter nor wettability has been identified. On graphene-coated surfaces, Andrews et al.23 studied the nanodroplet spreading using MD simulation. Unlike the partial wetting-transparency in static wetting, the dynamic wetting of nanodroplet was found to be independent of the underlying substrates; the contact radius followed R~t1/2 regardless of the number of graphene layers or the wetting property of the underlying substrate. The study shed light on the inertia-dominated dynamic wetting of nanodroplets on graphene-coated surfaces. However, the scaling that incorporates both dependences on diameter and wettability is still unknown, and thus we aim to identify such scaling in this work. Furthermore, we investigate the sensitivity of the dynamic wetting on the surface structure with the scope of controlling the wetting phenomena through the geometry of the graphene

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

layer, following series of works on macroscale droplets where large controllability by surface structures has been demonstarated.40,41 Here, the surface structure is altered by forming patterned defects in the graphene, which can be achieved using the ion beam technique.42,43

SIMULATION DETAILS MD simulations were performed to study the dynamic wetting of water droplet on a solid surface down to the nanoscale. Three sizes of water droplet with initial radii, Ri = 6.4, 9.0, and 12.8 nm, corresponding to 25,000, 50,000, and 100,000 water molecules, were considered to study the size effect. The model surfaces included graphene, graphite, monolayer graphene coated copper (MGC), bilayer graphene coated copper (BGC), and copper. For droplet at the nanoscale, the measured contact angle is affected by the line tension at the three-phase line.44 To eliminate the line tension effect,45 the quasi-two-dimensional (2D) system, droplet in cylindrical geometry, was implemented as shown in Fig. 1. The periodic boundary conditions were applied in all directions. The system size in the y direction was set to 5.92 nm, leading to an infinitely long cylindrical droplet. The system size in the x direction was set to 77.00 nm to ensure that the droplet does not interact with the periodic images. The cylindrical geometry is

6

ACS Paragon Plus Environment

Page 6 of 28

Page 7 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

widely implemented in the studies of dynamic wetting of nanodroplet23,27,31,35 to achieve sufficiently large droplet sizes in a hydrodynamic regime, which is computationally expensive in spherical geometry. Besides, the contact problems such as dynamic wetting are essentially 2D phenomena.35 The previous study31 also demonstrates that even though the evolution of contact radius is different for droplets in the spherical and cylindrical geometries, the local processes at the contact line are independent of the geometry of droplet.

Figure 1. Snapshot of the initial state of dynamic wetting simulation of water droplet (Ri=12.8nm) on monolayer graphene coated copper (MGC) surface.

The Cu(111) surface was chosen and the graphite surface was represented by six graphene

layers.

The

separation

distance

between

copper–graphene

and

graphene-graphene was 0.326 nm46 and 0.34 nm,44 respectively. The extended simple 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

point charge (SPC/E) model47 was applied for water droplets. The Lennard-Jones (LJ) potential was implemented for the van der Waals interaction between water and solid surfaces. The LJ parameters of the interaction between oxygen atom of water molecule and carbon atom were ε=0.392 kJ/mol and σ=0.319 nm,44 while the LJ parameters were ε=5.067 kJ/mol and σ=0.272 nm48 for the interaction between oxygen atom of water molecule and copper atom. The difference between the water contact angles on rigid and flexible substrates has been shown to be negligible,44 and thus all atoms of solid surfaces were kept fixed at their equilibrium positions to reduce the computational cost. All simulations were performed with the GROMACS 4.6.749 and visualized with the software PyMOL.50 The simulations were carried out in a canonical ensemble with integrating time steps of 1.0 fs. The temperature was kept constant at 300 K using the Nosé−Hoover thermostat51,52 with a time constant of 0.2 ps. A cutoff distance of 2.0 nm was used for the nonbonded van der Waals (vdW) interactions. The electrostatic interaction was treated using particle-mesh Ewald summation method53 with a real-space cutoff distance of 2.0 nm. The water droplet was first equilibrated in the system without solid surfaces to achieve the circular shape. The equilibrated droplet was then placed at a location with

8

ACS Paragon Plus Environment

Page 8 of 28

Page 9 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

a certain distance from the solid surface so that they start to interact. The droplet spreads on the solid surface spontaneously due to the solid-liquid interaction. To measure the contact radius of droplet, a square binning of local number of water molecules along the x and z directions with a bin size of 0.2 nm was performed. The liquid-gas interface was defined to be the contour of density = 500 kg/m3, which approximately corresponds to the middle value of bulk liquid and gas densities of water. The atom positions were recorded every picosecond to evaluate time histories of the contact radius evolution.

RESULTS AND DISCUSSION The static contact angle of droplet with Ri=6.4, 9.0, and 12.8 nm on graphene was calculated to be 89.6±3.5°, 89.9±3.5°, and 89.9±3.4°, respectively. The obtained contact angle on graphene is consistent with the previous simulation with the same model; the contact angle on graphene is around 90°.7,17 The static contact angles of droplet with Ri=12.8 nm are 64.5±3.5°, 78.4±2.5°, and 80.2±2.5° on MGC, BGC, and graphite, respectively. The copper surface is completely wetting, i.e. contact angle is zero. The values of static contact angles are in good agreement with the previous MD study12 using the same force field. The contact angle is increased significantly by the addition of a monolayer graphene. By adding two layers of graphene, the value of 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

contact angle approaches that of graphite, which indicates that the effect from underlying substrate can be shielded by bilayer graphene. Figure 2(a) shows how the contact radius of droplets with various initial radii evolves in time on surfaces with various wettabilities. The late-time spreading shows a dependence on droplet sizes that larger droplets spread faster than smaller ones. For surfaces with different wettabilities, the droplets spread faster on more hydrophilic surfaces. The observation is consistent with the microscale experimental results.38 The previous studies suggest that the spreading at the nanoscale operates in the inertial-dominated regime,23,35 thus the inertial scaling is adapted in Figure 2(b). The characteristic inertia time constant (ρRi3/γ)1/2, where ρ is the liquid density and γ is the liquid-vapor surface tension, comes from the balance of the inertial pressure inside the droplet and the capillary pressure associated with the droplet curvature.35 The values ρ=1000 kg/m3 and γ=72 mN/m are used in this work. The curves of droplets with different sizes on the same surface collapse onto one curve. However, the curves for surfaces with different wettabilities scatter, which was also observed in the previous macroscopic experiments.38 Now, since the final equilibrium radius, Rf, of droplet is related to the wettability of surface, we scale the radius by Rf instead of Ri, shown in Fig. 2(c). Note that the final equilibrium radius cannot be obtained for the copper

10

ACS Paragon Plus Environment

Page 10 of 28

Page 11 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

surface, thus the profiles of copper surface are not included in this scaling analysis. Interestingly, all simulation profiles collapse onto one single curve, which suggests that the appropriate length scale and time scale to describe the dynamic wetting at the nanoscale are Rf and (ρRi3/γ)1/2, respectively. We here call this scaling the “modified inertial scaling”.

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2. Time histories of contact radius of droplets with varying sizes on surfaces with varying wettability (a) dimensional units, (b) inertial scaling, and (c) modified inertial scaling. The numbers in parentheses in (a) denote the static contact angle of Ri=12.8 nm droplet. 12

ACS Paragon Plus Environment

Page 12 of 28

Page 13 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In mind of enhancing the controllability of dynamic wetting, one possible strategy is to nanostructure the surface by restructuring the coated monolayer graphene. To this end, patterned holes with various sizes and shapes are formed to understand the effect of geometries of the surface nanostructures. Note that intentional introduction of such defects are now possible in experiments using the ion beam technique.42,43 Each structure was denoted with R(δx, δy) and T(δx, δy), as shown in Fig. 3, where R and T denote rectangular and triangle and δx and δy are the length of the hole in the x and y directions integer-multiples of 0.43 nm and 0.37 nm, respectively. For example, R(3, 4) indicates that a 1.29 nm × 1.48 nm rectangular hole is periodically repeated 20 times in the x direction and 2 times in the y direction. Two underlying substrate, copper and graphite, are chosen to study the influence of underlying substrate with different wettabilities.

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Top view of the patterned graphene-coated copper surface of (a) rectangular R(3,4) and (b) triangle T(4,4). δx and δy denote the pattern sizes in the x and y directions, respectively.

Figure 4 shows the contact radius profiles of droplet with Ri = 12.8 nm spreading on various patterned graphene-coated surfaces. The results demonstrate the droplets spread more slowly on the patterned-graphene coated surface than that on the smooth graphene-coated surface. Comparison between Fig. 4(a) and Fig. 4(b) indicates that the extent of this wetting inhibition depends on the underlying substrate, where the inhibition is stronger for the case of graphite underlying substrate (Fig. 4(a)) than the copper underlying substrate (Fig. 4(b)). The validity of modified inertial scaling is examined to the patterned graphene-coated surface in Fig. 5. For the graphite underlying substrate case (Fig. 5(a)), the profiles of different sizes of droplet are also analyzed. The results show that the profiles do not collapse as well as in the case of 14

ACS Paragon Plus Environment

Page 14 of 28

Page 15 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

smooth graphene-coated surfaces. For the copper underlying substrate case (Fig. 5(b)), the profiles of patterned graphene-coated surfaces do not collapse and clearly differ from the profile of smooth graphene-coated surface. It indicates that the effect of the underlying substrate are enhanced because of the surface nanostructures. The importance of other effects, for example the line friction force, may be amplified on the patterned graphene-coated surfaces. Thus, although the modified inertial scaling works for droplets on the smooth graphene-coated surfaces, it no longer holds on the patterned graphene-coated surfaces.

15

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. Time histories of contact radius of droplets on various patterned graphene-coated surfaces with (a) graphite and (b) copper as underlying substrates.

16

ACS Paragon Plus Environment

Page 16 of 28

Page 17 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5. Time histories of contact radius of droplets on various patterned graphene-coated surfaces using modified inertial scaling with (a) graphite and (b) copper as underlying substrates.

In order to quantify the effect of pattern size on the dynamic wetting, the spreading speeds are calculated. The spreading speed is defined as the average contact line speed over the initial 500 ps (dash line in Fig. 4), i.e., dividing the contact radius at 500 ps by 500 ps. The multiple regression analysis is applied to relate the spreading speed to the

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

size of pattern, δx and δy, as shown in Fig. 6. Note that the results of triangular patterned cases are not shown in Fig. 6 because the effect of triangular pattern on dynamic wetting is difficult to separate in the x and y directions, respectively. The results reveal that the effect of surface nanostructures on dynamic wetting is more obvious on the graphite underlying substrate case than that on the copper one. The difference in spreading speed is 18.35% for the graphite underlying substrate case and it is only 7.77% for the copper one. It implies that the dynamic wetting can be manipulated more effectively by the design of pattern size for the graphite underlying substrate than the copper one. Also, it is shown that the spreading speed correlates to δy, which is perpendicular to the spreading direction, more than to δx, which is parallel to the spreading direction. The larger values in δy (larger pattern size in the y direction) induce more resistance for droplet spreading on both underlying substrates. The behavior can be attributed to the contact line pinning effect of spreading. The pinning effect has been shown to be important for the spreading on the textured surfaces with micro-structures.41,54 On nanostructured surfaces, the contact angle hysteresis shows the dependence on the area density of defects, indicating that the defects pin the contact line of droplet.55 It was also shown that the macroscopic pinning behavior is preserved for nanostructures with dimensions down to ~200 nm.56 On the other hand,

18

ACS Paragon Plus Environment

Page 18 of 28

Page 19 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

the contributions of δx to the dynamic wetting are opposite for these two underlying substrates, i.e. positive coefficient for the graphite underlying substrate case and negative coefficient for the copper one, as shown in Fig. 6. In order to depict the phenomena clearly, the top views of spreading of water adjacent to the surface, within 0.5 nm, are displayed in Fig. 7 for these two underlying substrates. For graphite underlying substrate case, the contact line is pinned on the boundary of hole in both the x and y directions, as shown from t = 500 ps to 560 ps. The pinning effect in the both directions can also be observed for copper underlying substrate case from t = 430 ps to 490 ps. When water molecules contact the copper surface (t = 520ps), the spreading is enhanced because copper is more hydrophilic than graphene. For graphite underlying substrate case, the larger values in both δx and δy form more resistance for droplet spreading due to the pinning effect. The slower spreading and higher hydrophobicity was also observed on the defective graphene surface.8 For copper underlying substrate case, the larger values in δy slow down the spreading while the larger values in δx enhance the spreading. Due to the competition between these two effects, the spreading processes are insensitive to the pattern sizes in the early regime.

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. Spreading speed of droplets as a function of pattern size, δx and δy, on various patterned graphene-coated surfaces with (a) graphite and (b) copper as underlying substrates.

20

ACS Paragon Plus Environment

Page 20 of 28

Page 21 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 7. Top view of spreading of water molecules adjacent to R(7,8) patterned graphene-coated surface with (a) graphite and (b) copper as underlying substrates. Coated graphene in green, graphite in gray, and copper in yellow.

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

CONCLUSION MD simulations are implemented to study the dynamics wetting of water droplet on graphene-coated surface with different underlying substrates, copper, and graphite. The influence from underlying substrate on the both static contact angle and dynamic spreading can be weakened significantly by the addition of one layer of graphene. The dimensional analysis demonstrates that the dynamic wetting can be characterized by the inertial time constant and final equilibrium radius for the nanodroplets. For patterned graphene-coated surfaces, the modified inertial scaling does not hold, indicating other effects, such as line friction force, become important due to the surface nanostructures. Unlike the smooth graphene-coated surfaces, the dynamic wetting shows strong dependence on the underlying substrate. For the graphite underlying substrate case, the spreading speed decreases and hydrophobicity increases with the increasing size of holes. The dependence of early-time spreading on pattern size is weaker for the copper underlying substrate case, because of the competition between the strong attraction from copper surface and pinning effect. The results indicate that the manipulation of spreading can be achieved by modifying the nanostructures of graphene on hydrophobic underlying substrate.

22

ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Notes The authors declare no competing financial interests.

AKNOWLEDGEMENTS This work was financially support by Grant-in-Aid for JSPS Fellows (Grant No. 26-04364). S.-W. H. was financially supported by JSPS Postdoctoral Fellowship for Overseas Researchers Program (26-04364).

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

REFERENCES (1) (2) (3) (4)

Geim, A. K. Graphene: Status and Prospects. Science 2009, 324, 1530–1534. Geim, A. K.; Novoselov, K. S. The Rise of Graphene. Nat. Mater. 2007, 6, 183–191. Lee, C.; Wei, X.; Kysar, J. W.; Hone, J. Measurement of the Elastic Properties and Intrinsic Strength of Monolayer Graphene. Science 2008, 321, 385–388. Elias, D. C.; Nair, R. R.; Mohiuddin, T. M. G.; Morozov, S. V; Blake, P.; Halsall, M. P.; Ferrari, A. C.; Boukhvalov, D. W.; Katsnelson, M. I.; Geim, A. K.; et al. Control of Graphene’s Properties by Reversible Hydrogenation: Evidence for

(5)

Graphane. Science 2009, 323, 610–613. Nine, M. J.; Cole, M. A.; Tran, D. N. H.; Losic, D. Graphene: A Multipurpose

(6)

Material for Protective Coatings. J. Mater. Chem. A 2015, 3, 12580–12602. Driskill, J.; Vanzo, D.; Bratko, D.; Luzar, A. Wetting Transparency of Graphene

(7)

(8) (9)

in Water. J. Chem. Phys. 2014, 141, 18C517. Hung, S.-W.; Hsiao, P.-Y.; Chen, C.-P.; Chieng, C.-C. Wettability of Graphene-Coated Surface: Free Energy Investigations Using Molecular Dynamics Simulation. J. Phys. Chem. C 2015, 119, 8103–8111. Li, X.; Li, L.; Wang, Y.; Li, H.; Bian, X. Wetting and Interfacial Properties of Water on the Defective Graphene. J. Phys. Chem. C 2013, 117, 14106–14112. Li, Z.; Wang, Y.; Kozbial, A.; Shenoy, G.; Zhou, F.; McGinley, R.; Ireland, P.; Morganstein, B.; Kunkel, A.; Surwade, S. P.; et al. Effect of Airborne Contaminants on the Wettability of Supported Graphene and Graphite. Nat.

(10)

Mater. 2013, 12, 925–931. Ondarçuhu, T.; Thomas, V.; Nuñez, M.; Dujardin, E.; Rahman, A.; Black, C. T.;

(11)

Checco, A. Wettability of Partially Suspended Graphene. Sci. Rep. 2016, 6, 24237. 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.

(12)

Raj, R.; Maroo, S. C.; Wang, E. N. Wettability of Graphene. Nano Lett. 2013, 13, 1509–1515.

(13)

Ramos-Alvarado, B.; Kumar, S.; Peterson, G. P. On the Wettability Transparency

(14)

of Graphene-Coated Silicon Surfaces. J. Chem. Phys. 2016, 144, 14701. Shih, C.-J.; Strano, M. S.; Blankschtein, D. Wetting Translucency of Graphene.

(15)

Nat. Mater. 2013, 12, 866–869. Shih, C.-J.; Wang, Q. H.; Lin, S.; Park, K.-C.; Jin, Z.; Strano, M. S.; 24

ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Blankschtein, D. Breakdown in the Wetting Transparency of Graphene. Phys. (16)

Rev. Lett. 2012, 109, 176101. Shin, Y. J.; Wang, Y.; Huang, H.; Kalon, G.; Wee, A. T. S.; Shen, Z.; Bhatia, C. S.; Yang, H. Surface-Energy Engineering of Graphene. Langmuir 2010, 26, 3798–3802.

(17)

Taherian, F.; Marcon, V.; van der Vegt, N. F. A.; Leroy, F. What Is the Contact

(18)

Angle of Water on Graphene? Langmuir 2013, 29, 1457–1465. Vanzo, D.; Bratko, D.; Luzar, A. Wettability of Pristine and

(19)

Alkyl-Functionalized Graphane. J. Chem. Phys. 2012, 137, 34707. Wang, S.; Zhang, Y.; Abidi, N.; Cabrales, L. Wettability and Surface Free

(20)

Energy of Graphene Films. Langmuir 2009, 25, 11078–11081. Bonn, D.; Eggers, J.; Indekeu, J.; Meunier, J.; Rolley, E. Wetting and Spreading.

(21)

Rev. Mod. Phys. 2009, 81, 739–805. Yang, M.; Cao, B.-Y.; Wang, W.; Yun, H.-M.; Chen, B.-M. Experimental Study on Capillary Filling in Nanochannels. Chem. Phys. Lett. 2016, 662, 137–140

(22) Blake, T. D. The Physics of Moving Wetting Lines. J. Colloid Interface Sci. 2006, 29, 1–13. (23) Andrews, J. E.; Sinha, S.; Chung, P. W.; Das, S. Wetting Dynamics of a Water Nanodrop on Graphene. Phys. Chem. Chem. Phys. 2016, 18, 23482–23493. (24) Bertrand, E.; Blake, T. D.; Coninck, J. De. Influence of Solid–liquid Interactions on Dynamic Wetting: A Molecular Dynamics Study. J. Phys. Condens. Matter 2009, 21, 464124. (25) De Coninck, J.; Blake, T. D. Wetting and Molecular Dynamics Simulations of Simple Liquids. Annu. Rev. Mater. Res. 2008, 38, 1–22. (26) de Ruijter, M. J.; Blake, T. D.; De Coninck, J. Dynamic Wetting Studied by Molecular Modeling Simulations of Droplet Spreading. Langmuir 1999, 15, 7836–7847. (27) He, G.; Hadjiconstantinou, N. G. A Molecular View of Tanner’s Law: Molecular Dynamics Simulations of Droplet Spreading. J. Fluid Mech. 2003, 497, 123–132. (28) Heine, D. R.; Grest, G. S.; Webb, E. B. Spreading Dynamics of Polymer Nanodroplets in Cylindrical Geometries. Phys. Rev. E 2004, 70, 11606. (29) Heine, D. R.; Grest, G. S.; Webb, E. B. Surface Wetting of Liquid Nanodroplets: Droplet-Size Effects. Phys. Rev. Lett. 2005, 95, 107801. (30) Johansson, P.; Carlson, A.; Hess, B. Water-Substrate Physico-Chemistry in Wetting Dynamics. J. Fluid Mech. 2015, 781, 695–711. (31) Nakamura, Y.; Carlson, A.; Amberg, G.; Shiomi, J. Dynamic Wetting at the 25

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nanoscale. Phys. Rev. E 2013, 88, 33010. (32) Samsonov, V. M.; Ratnikov, A. S. Comparative Molecular Dynamics Study of Simple and Polymer Nanodroplet Spreading. Colloids Surfaces A Physicochem. Eng. Asp. 2007, 298, 52–57. (33) Seveno, D.; Dinter, N.; De Coninck, J. Wetting Dynamics of Drop Spreading. New Evidence for the Microscopic Validity of the Molecular-Kinetic Theory. Langmuir 2010, 26, 14642–14647. (34) Biance, A.-L.; Clanet, C.; Quéré, D. First Steps in the Spreading of a Liquid Droplet. Phys. Rev. E 2004, 69, 16301. (35) Winkels, K. G.; Weijs, J. H.; Eddi, A.; Snoeijer, J. H. Initial Spreading of Low-Viscosity Drops on Partially Wetting Surfaces. Phys. Rev. E 2012, 85, 55301. (36) Duchemin, L.; Eggers, J.; Josserand, C. Inviscid Coalescence of Drops. J. Fluid Mech. 2003, 487, 167–178. (37) Eggers, J.; Lister, J. R.; Stone, H. A. Coalescence of Liquid Drops. J. Fluid Mech. 1999, 401, 293–310. (38) Bird, J. C.; Mandre, S.; Stone, H. A. Short-Time Dynamics of Partial Wetting. Phys. Rev. Lett. 2008, 100, 234501. (39) Samsonov, V. M.; Zhukova, N. A.; Dronnikov, V. V. Molecular-Dynamic Simulation of Nanosized Droplet Spreading over a Continual Solid Surface. Colloid J. 2009, 71, 835–845. (40) Nita, S.; Do-Quang, M.; Wang, J.; Chen, Y.; Suzuki, Y.; Amberg, G.; Shiomi, J. Electrostatic Cloaking of Surface Structure for Dynamic Wetting. Sci. Adv. 2017, 3, e1602202. (41) Wang, J.; Do-Quang, M.; Cannon, J. J.; Yue, F.; Suzuki, Y.; Amberg, G.; Shiomi, J. Surface Structure Determines Dynamic Wetting. Sci. Rep. 2015, 5, 8474. (42) Emmrich, D.; Beyer, A.; Nadzeyka, A.; Bauerdick, S.; Meyer, J. C.; Kotakoski, J.; Gölzhäuser, A. Nanopore Fabrication and Characterization by Helium Ion Microscopy. Appl. Phys. Lett. 2016, 108, 163103. (43) Kalhor, N.; Boden, S. A.; Mizuta, H. Sub-10 Nm Patterning by Focused He-Ion Beam Milling for Fabrication of Downscaled Graphene Nano Devices. Microelectron. Eng. 2014, 114, 70–77. (44) Werder, T.; Walther, J. H.; Jaffe, R. L.; Halicioglu, T.; Koumoutsakos, P. On the Water−Carbon Interaction for Use in Molecular Dynamics Simulations of Graphite and Carbon Nanotubes. J. Phys. Chem. B 2003, 107, 1345–1352. (45) Weijs, J. H.; Marchand, A.; Andreotti, B.; Lohse, D.; Snoeijer, J. H. Origin of 26

ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Line Tension for a Lennard-Jones Nanodroplet. Phys. Fluids 2011, 23, 22001. (46) Xu, Z.; Buehler, M. J. Interface Structure and Mechanics between Graphene and Metal Substrates: A First-Principles Study. J. Phys. Condens. Matter 2010, 22, 485301. (47) Berendsen, H. J. C.; Grigera, J. R.; Straatsma, T. P. The Missing Term in Effective Pair Potentials. J. Phys. Chem. 1987, 91, 6269–6271. (48) Zhu, S.; Philpott, M. Interaction of Water with Metal Surfaces. Jounal Chem. Phys. 1994, 100, 6961. (49) Pronk, S.; Páll, S.; Schulz, R.; Larsson, P.; Bjelkmar, P.; Apostolov, R.; Shirts, M. R.; Smith, J. C.; Kasson, P. M.; van der Spoel, D.; et al. GROMACS 4.5: A High-Throughput and Highly Parallel Open Source Molecular Simulation Toolkit. Bioinformatics 2013, 29, 845–854. (50) Schrodinger LLC. The PyMOL Molecular Graphics System, Version 0.99. 2015. (51) Hoover, W. G. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A 1985, 31, 1695–1697. (52) Nose, S. A Unified Formulation of the Constant Temperature Molecular Dynamics Methods. J. Chem. Phys. 1984, 81, 511. (53) Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An N⋅log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089. (54) Paxson, A. T.; Varanasi, K. K. Self-Similarity of Contact Line Depinning from Textured Surfaces. Nat. Commun. 2013, 4, 1492. (55) Ramos, S.; Tanguy, A. Pinning-Depinning of the Contact Line on Nanorough Surfaces. Eur. Phys. J. E 2006, 19, 433–440. (56) Wong, T. S.; Huang, A. P. H.; Ho, C. M. Wetting Behaviors of Individual Nanostructures. Langmuir 2009, 25, 6599–6603.

27

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Graphic

28

ACS Paragon Plus Environment

Page 28 of 28