## Wetting behavior of supercooled water droplets impinging on nanostructured graphite surface: A molecular dynamics study

## Amir Afshar and Dong Meng

Mississippi State University, Starkville, MS, 39762, USA

In this work, we studied the wetting behavior of impinged and stagnant supercooled water nanodroplet at the atomistic scale using molecular dynamics simulations. We show that water droplet represents  a  retraction  behavior  from  surface  groove  nanostructure  after  indicating  maximum penetration.  There  is  a  threshold  value  for  the  value  of  impinging  velocity,  v*,  below  which penetration is not observed which depends on size of the surface roughness. This indicates grooves smaller  than  a  certain  size  can  stop  water  from  penetration,  considering  v~0.5M  in  practical situations, which is such a useful information for design of nanostructured ice-phobic substrates. The resistance  presumably  comes  from  interfacial  tension  between  solid  substrate  and  liquid  water droplet. The extent of water penetration significantly depends on droplet impinging velocity and the size of surface roughness, for impinging velocities higher than V*. Our result indicates the dominant role of droplet temperature on time required droplet recede from surface groove structure rather than droplet impinging velocity and the surface groove nanostructure.

## 1. Introduction

Water droplets spreading and impacting on a solid surface represent substantial dynamic wetting phenomena in which there is an intersection between physics, chemistry and mechanics. They are pertinent for engineering applications such as surface coating and inkjet printing . While 1 static wetting is fairly well explained by Young's contact angle equation, the dynamics of water droplet wetting is more intricated, since the involved physical processes in the local region near the  moving  contact  line  become  significant 2,3 .  The  wetting  behavior  of  water  droplet  on  solid surfaces is thoroughly understood on the macroscopic scale 4,5 ; however there are many other new wetting phenomena happened at the nanoscale which is entirely different from the behaviors at the micro- and mesoscales. For numerous nanoscale applications the wetting of fluids on different surfaces, such as graphite  and metals  has been investigated. 6 7

Theoretically,  for an inclusive understanding of the impacting and spreading process of water  droplet  on  solid  surface,  some  significant  physical  property  name  as  droplet  viscosity, inertial effect and surface tension between water droplet and solid surface needs to be considered. In terms of numerical modeling in approaching dynamic wetting some different methods such as volume of fluid (VOF) 8 and finite elements modeling (FEM)  have been used. These frequently 9 used numerical methods provide information of the wetting process at the mesoscale level. In this work, we employed molecular dynamics (MD) simulation method which delivers information of dynamic wetting at the atomistic level. To the best of our knowledge, there have so far been no MD simulation incorporated to investigate the dynamic wetting of impinged water droplet on solid surfaces which has a numerous application in aeronautics, wind turbines and solar panels industries to name a few. At the nanoscale, MD is the best computational method to examine and predict the dynamic wetting properties of nanodroplets on solid surfaces 10 .

On the nanoscale, the wetting property of a solid surface is significantly investigated by the  atomistic  interactions  between  the  liquid  and  the  surface.  Specifically,  the  predication  of wettability  for  a  system  consisting  liquid  droplet  and  solid  surface,  necessarily  requires  some detailed system information but not limited to surface structure, system temperature and chemical composition, which is implicitly contained in the atomistic simulation of materials. Thus, MD is an  appropriate  simulation  tool  to  explore  the  nature  of  the  surface  wetting  resulting  from  the interactions between  liquid droplet and solid surface at the atomistic level. There are some works used MD simulations to study the wetting behavior of nanodroplet on various surfaces. Sedighi et. al. 11 and Park et. el. 12 studied the effect of strength of interfacial energy with respect to the dropletsurface  and found increasing the interaction energy  resulted in significant enhancement on the wettability.

The main objective of the present work is to investigate the dynamic wetting behaviors of impinged supercooled water nanodroplet on textured graphite surfaces using MD simulations. For this, the effect of system physical properties namely as surface roughness, droplet temperature and impinging velocity on dynamic wetting are also examined. The results are discussed in terms of evolution of depth of water penetration, droplet retracting time from surface and the maximum depth of penetration in which droplet experienced during dynamic wetting.

## 2. Model and simulation method

In this study, we modeled water molecules using TIP4P/Ice atomistic water model 13 . In comparison  with  other  transferrable  atomistic  water  models  like  as  TIP4P/2005 14 , TIP4P 15 , TIP3P 16 and  SPC 17 ,  this  water  model  appropriately  reproduces  structure  and  phase  transition temperature of ice. It accurately  reproduces transition temperature of 269.8 ±1 K 18 at  1  bar  and density  of  0.906  g/cm   for  hexagonal  ice  which  is  favorably  comparable  to  the  corresponding 3 experimental  values.  Thus,  this  model  is  suitable  for  investigation  of  wetting  behavior  of supercooled liquid water on nanotextured substrate. Interaction parameters of TIP4P/ice force field are listed in Table 1 where dOM is the distance between the oxygen atoms and the massless charged point in the water model, and qH is the electron charge carried by each hydrogen atom in the water molecules.

Table 1. Interaction parameters of the water model 13

| Water     | 𝜺 𝒐-𝒐 [kcal   | 𝝈 𝑶-𝑶 [Å]   | q H (e)   | d OM [Å]   |
|-----------|---------------|-------------|-----------|------------|
| model     | mol -1 ]      |             |           |            |
| TIP4P/Ice | 0.21084       | 3.1668      | 0.5897    | 0.1577     |

Table 2. Lennard-Jones parameters for aromatic carbon atoms from OPLS force field 19

| Graphene carbon   |   𝝈 𝑪-𝑪 [Å] |   𝜺 𝑪-𝑪 [kcal mol -1 ] |
|-------------------|-------------|------------------------|
|                   |        3.55 |                   0.07 |

In our study, we introduced graphite blocks as the substrate. It consists of several graphene pillars with interlayer distance of 0.335nm. Each graphene sheet is made by several carbon atoms in a honey comb molecular structure which are bonded to their nearest neighbors at the equilibrium bond length and angle of 0.1418nm and 120  respectively ◦ 20 . All graphene sheets are electronically neutral during our simulation and are extended along 'z' direction of simulation box. They are modeled using OPLS (Optimized Potential for Liquid Simulation) force field with non-bonded interaction parameters given in Table 2.

Interaction between water molecules and the nanotextured graphite substrate at atomistic level  are  considered  through  the  Lennard-Jones  (LJ)  non-bonded  interaction  potential  existing between oxygen and carbon atoms in our system. Equation (1) shows the LJ potential governs the interactions existing in the system. In our study, we adopted energy well depth (ε) and the interparticle distance (σ) between oxygen  and carbon atoms by  the Lorentz-Berthelot  mixing rule 21 given  by  Equation  (2)  and  (3),  which  we  employed  in  our  fracture  work 22 and  also  the corresponding values reported by Werder et. al (2003) . 6

Table  3  provides  the  values  of  non-bonded  interaction  potentials  between  oxygen  and carbon atoms in our system which are interacted by LJ potential.

$$U _ { c o } = 4 \varepsilon _ { c - o } \left [ \left ( \frac { \sigma _ { c - o } } { r } \right ) ^ { 1 2 } - \left ( \frac { \sigma _ { c - o } } { r } \right ) ^ { 6 } \right ]$$

$$\sigma _ { c - o } \, = \, \frac { 1 } { 2 } \, ( \sigma _ { c - c } \, + \, \sigma _ { o - o } )$$

$$\varepsilon _ { c - o } = \sqrt { \varepsilon _ { c - c } + \varepsilon _ { o - o } }$$

All  simulations  performed  using  the  open-source  molecular  dynamics  (MD)  parallel simulator code, LAMMPS 23 . Long-range electrostatic interactions between water/ice molecules are calculated using particle-particle particle-mesh (PPPM) method 24 . Water molecules are treated as the rigid body by keeping the bond length and bond angle fixed using the SHAKE algorithm 25 . In our study, graphene pillars are kept fixed in their lattice positions throughout the simulations. It has been verified that keeping graphene sheets fixed during the simulation does not influence the wetting behavior of supercooled liquid water on graphite substrate; though it remarkably reduces the computational cost 6,26 .

In this study, equilibrated structure of water droplet is obtained by supercooling a slab of pre-equilibrated ice Ih 27 at different supercooling temperatures of 250K and 265K. The modeled hexagonal Ih ice consists of 8640 number of water molecules has a prismatic face and surface area of 6.6 × 9.1 nm 2 .  This face  was brought in contact with graphite substrates during the wetting simulation. The thickness of the ice cube is 4.9 nm being perpendicular to the ice-substrate contact area.

To  study  the  wetting  behavior  of  supercooled  liquid  water  droplet  on  nano-textured graphite substrate, we initially tuned the hydrophobicity  of graphite substrate by measuring the equilibrium  contact  angle  of  liquid  water  droplet  on  graphene  sheets.  To  investigate  the equilibrium  contact  angle  of  liquid  water  droplet  on  graphene  sheets,  we  placed  the  preequilibrated  h ice in 2.0nm perpendicular to five horizontal graphene sheets. Following that, ice I cube was melted under NVT constraint at 300K for 2.0ns. Figure  1 shows a snapshot from a molecular structure of the constructed system having equilibrated Ih ice and water droplet on top of graphene sheets.

Fig.1. structure of water molecules on top of five graphene sheets; (a) Ih ice, (b) liquid water droplet at 300K.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00560v1-figure-36bb3a0d-d6b5-438a-9455-1914c3c83866.jpg)

In follows, to measure the equilibrium contact angle, the equilibrated liquid water droplet initially placed on the closest possible distance with graphene sheet where there is no atomistic overlap, through energy  minimization process. Then the water droplet is equilibrated for 1.0ns under NVT simulation at 300K. In our study, interactions between water-substrate is calculated by setting force field parameters between oxygen (O) and carbon (C) atoms of graphene sheets. To model the van der Waals forces between water molecules and carbon atoms of graphene layers, we set the interaction parameters using the mixing rule we applied in our fracture study 22 and the

values  which  are  reported  by  Werder  et.  al.     Table  3  shows  the  force  field  interaction  LJ 6 parameters  between  ice  and  graphite  substrate  which  governs  the  interaction  potential  of  the system.

Table 3. Non-bonded LJ parameters between oxygen and carbon atoms

|               |   𝜺 𝑶-𝑪 [kj mol-1] 𝝈 𝑶-𝑪 |    [Å] |
|---------------|--------------------------|--------|
| Mixing rule   |                   0.5084 | 3.3584 |
| Werder et. al |                   0.392  | 3.19   |

Figure  2  illustrates  the  equilibrium  states  of  liquid  water  droplet  on  five  horizontal graphene sheets using different force field parameters, as reported in Table 3. Since it illustrates, the  water  equilibrium  contact  angle  is  between  50 -60   using  the  values  from  mixing  rules. ◦ ◦ Experiments by Wang et al 28 , and Shin et al. 29 indicate a water droplet equilibrium contact angle of 90-95  on graphite substrate, which is fairly well reproduced in our study by applying the force ◦ field parameters from the work of Werder et. al . It should be noted that the equilibrium water 6 droplet contact angle on graphene sheets strongly depends on the number of graphene sheets, size of the water droplet, atomistic water model and the system temperature 30

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00560v1-figure-3b3589b4-519d-4245-a5f6-30eb87cc8deb.jpg)

Figure 2. contact angle of liquid water droplet on graphene sheets using the force field parameters; (a) based on mixing rule, (b) work of Werder et al. 6

Based on the obtained equilibrium contact angle of liquid water droplet on graphene sheets through using the developed force field parameters, we designed our simulation setup to study the wetting  behavior  of  impinging  supercooled  liquid water  droplet  on  nanostructure  hydrophobic graphite  substrate.  To  investigate  the  wetting  profile  of  supercooled  liquid  water  droplet,  we initially placed prismatic facet of pre-equilibrated ice I h in a distance 2nm perpendicular to the graphene pillars. To ensure that the atoms do not interact with their periodic images at any time during the simulation, x-dimension of simulation box is increased. The ice cube was supercooled at  two  supercooling  temperatures  of 250K and 265K. In  our study,  graphite  substrate  built  by several zig-zag edge type of graphene pillars each having height and width of 5.0nm and 6.6nm respectively along 'x' and 'y' directions which are extended along 'z' direction by 3.35Å 31 .

We  introduced  groove  structures  to  the  graphite  substrate  by  making  some  distances between graphite blocks called as 'd' value. Based on the value of 'd' parameter of 2nm and 4.7nm in this study, number of graphene pillars in each graphite block varied from 22, at d=2nm, to  14  sheets  at  d=4.7nm.  Figure  3  illustrates  a  molecular  schematic  of  the  constructed  system consists of pre-equilibrated ice and supercooled liquid water droplet at 250K on top of graphite substrate having 'd' equals to 2nm and 4.7nm

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00560v1-figure-0204d7dd-7501-4330-bde7-9746bea88ab7.jpg)

Figure 3. Molecular representation of  the  constructed model  system  for wetting simulation; (a) initial system with introducing Ih ice on top of graphite substrate, (b) supercooled liquid water at 250K on top of graphite blocks with 'd' at 2.0nm, (c) supercooled liquid water at 250K on top of graphite blocks with 'd' at 4.7nm

To investigate the effect of impinging velocity on wetting profile of supercooled liquid water droplet, we studied two different systems of impinging state and reference state, where there is no impact velocity. In reference state, equilibrated supercooled liquid water droplet brings into contact with the graphite substrate through energy minimization simulations. In impinging state, compared to our previous study 32 , we introduced collective impinging velocity in terms of Much numbers (M) to the thermal velocities of equilibrated supercooled liquid water in a range of 0.5M, 1.0M, 1.5M, 2.0M, and 2.5M. Simulations are all conducted under NVT constraint with periodic boundary conditions in all directions of simulation box for 10ns with timestep of 0.001ps.

## 3. Results and Discussions

## 3.1. wetting profiles of water droplet

Figure 4 shows a molecular illustration of impinged supercooled liquid water with impact velocity of  0.5M and 2.5M at droplet temperature of 250K on graphite substrates having different surface roughness of 2.0nm and 4.7nm.

To calculate the depth of supercooled water penetration (H) into surface groove nanostructure, we used the developed formula as:

$$H \equiv \frac { 2 \int _ { x = 0 } ^ { x = l g } x \rho ( x ) d x } { \int _ { x = 0 } ^ { x = l g } \rho ( x ) \, d x }$$

In this expression, lg is the length of single graphene sheet, and 𝜌(𝑥) is the density of penetrated water molecules (specifically oxygen atoms) confined by graphene sheets which is varying during the simulation. This is calculated by the given equation 5, in which boundaries along 'z' direction is equal to the size of the surface groove structure in each substrate and 𝜌𝑂 stands for density of oxygen atoms of water molecules.

$$\rho ( x ) \, = \, \frac { 1 } { L _ { y } L _ { z } } \left [ \iiint _ { y l o } ^ { y h i } \rho _ { o } ( x, y, z ) \, d y d z \right ]$$

Figure 4 . schematic of wetting behavior of supercooled liquid water on top of graphite substrate at different 'groove' sizes and impinging velocities at 0.2ns; (a) impinging velocity of 0.5M, d= 2.0nm, (b) impinging velocity of 2.5M, d= 2.0nm, (c) impinging velocity of 0.5M, d = 4.7nm, (d) impinging velocity of 2.5M, d = 4.7nm.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00560v1-figure-4df1b78e-64fb-41d7-bd6b-cdb05d47e604.jpg)

Figure 5 indicates the evolution of depth of penetration for supercooled liquid water during wetting simulation at different droplet temperatures of 250K and 265K with d=2.0nm and 4.7nm. It shows supercooled water droplet demonstrates a gradual retreating behavior from surface texture

structure  once  it  penetrates  through  the  substrate.  This  reflects  the  inherent  hydrophobicity  of graphite substrate as we modeled by calibrating the water-graphene interaction strength (εco). Since it demonstrates, retreating behavior is substantially depend on droplet impinging velocity and the size of surface roughness parameter. Also results shows no significant different in wetting profile of supercooled water droplet having impinging velocity equal and less than 1.0M. In this study, we define impinging velocity of 1.5M as the critical impinging velocity that results in fundamental change in wetting  behavior of supercooled  droplet on graphite  substrate. Having  this in mind, water droplet having higher impinging velocity shows higher depth penetration on surface having bigger  texture  size,  which  is  true  at  different  supercooling  temperature.  To  compare  wetting dynamics of impinged supercooled water with static wetting, where there is no impact droplet velocity,  we  conducted  simulations  for  static  supercooled  liquid  water  droplets  on  graphite substrate at different temperatures of 250K and 265K and provided the wetting profile by the black 'triangle' data points in this figure. Comparing the dynamic with static wetting profiles, it indicates similar wetting profile for impinged water droplet after some simulation times called as 'retracting time' which is strongly depend on the droplet and surface physical properties.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00560v1-figure-fa7ed69a-1497-4be8-924c-106c42715678.jpg)

Figure 5. Wetting profiles of supercooled liquid water as a function of simulation time for (a) d=20Å, T = 250K, (b) d = 20Å, T=265K, (c) d=47Å, T = 250K, (d) d=47Å, T = 265K.

## 3.2. Initial penetration and retracting time profiles of water droplet

Figure  6  illustrates  the  initial  penetration  of  supercooled  liquid  water  as  a  function  of impinging velocity and retracting time (ns) versus initial penetration. In these plots, we define the maximum H value represented in figure 5 as the initial penetration and retracting time as the time required dynamic wetting indicates similar static wetting.

Figure 6. Initial penetration vs. imping velocity and retracting time vs. initial penetration.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00560v1-figure-5c5e02e5-ce38-4bc0-8709-ce4dd04aeefd.jpg)

These results demonstrate maximum penetration substantially depends on size of surface roughness and retracting time is significantly  influenced by  the droplet temperature. Results of maximum penetration vs. impinging velocity verifies the results shown in figure 5. It indicates

droplet temperature and surface roughness initiates to indicate their influence on maximum depth of water penetration (MDP) starting at 1.5M in which MDP is noticeably increased by mounting impinging velocity and surface roughness. Result of retracting time represents supercooled water droplets requires higher time to be retracted from groove structure at lower temperature of 250K. We can address this to the potential of the TIP4P/ice atomistic water model we used to indicate the effect of temperature on viscosity of supercooled water droplet. Also, it illustrates  retracting time is a strongly depend on droplet's temperature and the surface texture size. Figure 7 shows droplet retracting time as a function of impinging velocity. It demonstrates results shown in figure 7. Considering 1.5M as the critical impinging velocity, it verifies the dominant effect of droplet temperature  on  retracting  profile  of  supercooled  water  from  surface  nanotexture  structure.    It indicates water impinging velocity has a considerable effect on retracting profile of water droplets having lower temperature. We can use all this useful information at atomistic scale in designing anti-icing surface at macroscale level.

Figure 7. Retracting time vs. supercooled impact velocity for different surface 'groove' sizes and droplet's temperature.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00560v1-figure-f6088042-0061-42eb-bcd1-e53c71149d90.jpg)

Figure 8. Retracting time as a function of initial penetration of impinged water droplet relative to the nonimpinged droplet.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00560v1-figure-30554af0-b9c3-4305-9d02-1bca2221971f.jpg)

Figure 8 illustrates the retracting as a function of 'pure' depth of penetration, where we subtracted the MDP values from the non-impacted water droplet. Result demonstrates the strong effect of supercooled water droplet's temperature on retracting time, which is remarkably higher at 250K compared to 265K at two 'd' values. This supports the results illustrated in figure 8. In detail, it demonstrates at a specific 'd' value, supercooled water droplet with lower temperature requires more time to retract from surface groove structure.

## References

- (1) Bonn, D.; Eggers, J.; Indekeu, J.; Meunier, J.; Rolley, E. Wetting and Spreading. Rev. Mod. Phys. 2009, 81 (2), 739-805. https://doi.org/10.1103/RevModPhys.81.739.
- (2) Sui, Y.; Ding, H.; Spelt, P. D. M. Numerical Simulations of Flows with Moving Contact Lines. Annu. Rev. Fluid Mech. 2014, 46 (1), 97-119. https://doi.org/10.1146/annurev-fluid-010313141338.
- (3) Snoeijer, J. H.; Andreotti, B. Moving Contact Lines: Scales, Regimes, and Dynamical Transitions. Annu. Rev. Fluid Mech. 2013, 45 (1), 269-292. https://doi.org/10.1146/annurev-fluid-011212140734.
- (4) Stillinger, F. H. Water Revisited. Science 1980, 209 (4455), 451-457. https://doi.org/10.1126/science.209.4455.451.
- (5) Kropman, M. F.; Bakker, H. J. Dynamics of Water Molecules in Aqueous Solvation Shells. Science 2001, 291 (5511), 2118-2120. https://doi.org/10.1126/science.1058190.
- (6) 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 (6), 1345-1352. https://doi.org/10.1021/jp0268112.
- (7) 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. Mater. 2013, 12 (10), 925-931. https://doi.org/10.1038/nmat3709.
- (8) Malgarinos, I.; Nikolopoulos, N.; Marengo, M.; Antonini, C.; Gavaises, M. VOF Simulations of the Contact Angle Dynamics during the Drop Spreading: Standard Models and a New Wetting Force Model. Adv. Colloid Interface Sci. 2014, 212, 1-20. https://doi.org/10.1016/j.cis.2014.07.004.
- (9) Sprittles, J. E.; Shikhmurzaev, Y. D. Finite Element Simulation of Dynamic Wetting Flows as an Interface Formation Process. J. Comput. Phys. 2013, 233, 34-65. https://doi.org/10.1016/j.jcp.2012.07.018.
- (10)   Zhang, J.; Borg, M. K.; Reese, J. M. Multiscale Simulation of Dynamic Wetting. Int. J. Heat Mass Transf. 2017, 115, 886-896. https://doi.org/10.1016/j.ijheatmasstransfer.2017.07.034.
- (11)   Wu, C.-D.; Kuo, L.-M.; Lin, S.-J.; Fang, T.-H.; Hsieh, S.-F. Effects of Temperature, Size of Water Droplets, and Surface Roughness on Nanowetting Properties Investigated Using Molecular Dynamics Simulation. Comput. Mater. Sci. 2012, 53 (1), 25-30. https://doi.org/10.1016/j.commatsci.2011.09.024.
- (12)   Park, J.-H.; Inam, E.; Abdullah, M. H.; Agustiyani, D.; Duan, L.; Hoang, T. T.; Kim, K.-W.; Kim, S. D.; Nguyen, M. H.; Pekthong, T.; et al. Implications of Rainfall Variability for Seasonality and Climate-Induced Risks Concerning Surface Water Quality in East Asia. J. Hydrol. 2011, 400 (3), 323-332. https://doi.org/10.1016/j.jhydrol.2011.01.050.
- (13)   Abascal, J. L. F.; Sanz, E.; García Fernández, R.; Vega, C. A Potential Model for the Study of Ices and Amorphous Water: TIP4P/Ice. J. Chem. Phys. 2005, 122 (23), 234511. https://doi.org/10.1063/1.1931662.
- (14)   Abascal, J. L. F.; Vega, C. A General Purpose Model for the Condensed Phases of Water: TIP4P/2005. J. Chem. Phys. 2005, 123 (23), 234505. https://doi.org/10.1063/1.2121687.
- (15)   Huggins, D. J. Correlations in Liquid Water for the TIP3P-Ewald, TIP4P-2005, TIP5P-Ewald, and SWM4-NDP Models. J. Chem. Phys. 2012, 136 (6), 064518. https://doi.org/10.1063/1.3683447.
- (16)   Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79 (2), 926-935. https://doi.org/10.1063/1.445869.
- (17)   Berendsen, H. J. C.; Postma, J. P. M.; Gunsteren, W. F. van; Hermans, J. Interaction Models for Water in Relation to Protein Hydration. In Intermolecular Forces; The Jerusalem Symposia on Quantum Chemistry and Biochemistry; Springer, Dordrecht, 1981; pp 331-342. https://doi.org/10.1007/978-94-015-7658-1\_21.

- (18)   Conde, M. M.; Rovere, M.; Gallo, P. High Precision Determination of the Melting Points of Water TIP4P/2005 and Water TIP4P/Ice Models by the Direct Coexistence Technique. J. Chem. Phys. 2017, 147 (24), 244506. https://doi.org/10.1063/1.5008478.
- (19)   Pykal, M.; Jurečka, P.; Karlický, F.; Otyepka, M. Modelling of Graphene Functionalization. Phys. Chem. Chem. Phys. 2016, 18 (9), 6351-6372. https://doi.org/10.1039/C5CP03599F.
- (20)   Tsai, J.-L.; Tzeng, S.-H.; Tzou, Y.-J. Characterizing the Fracture Parameters of a Graphene Sheet Using Atomistic Simulation and Continuum Mechanics. Int. J. Solids Struct. 2010, 47 (3), 503-509. https://doi.org/10.1016/j.ijsolstr.2009.10.017.
- (21)   Nikitin, A.; Milchevskiy, Y.; Lyubartsev, A. AMBER-Ii: New Combining Rules and Force Field for Perfluoroalkanes. J. Phys. Chem. B 2015, 119 (46), 14563-14573. https://doi.org/10.1021/acs.jpcb.5b07233.
- (22)   Afshar, A.; Zong, J.; Thompson, D. S.; Meng, D. Fracture of Ice at Interfaces from Molecular Dynamics Simulations. In 2018 Atmospheric and Space Environments Conference; American Institute of Aeronautics and Astronautics. https://doi.org/10.2514/6.2018-3017.
- (23)   Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117 (1), 1-19. https://doi.org/10.1006/jcph.1995.1039.
- (24)   Luty, B. A.; van Gunsteren, W. F. Calculating Electrostatic Interactions Using the Particle-Particle Particle-Mesh Method with Nonperiodic Long-Range Interactions. J. Phys. Chem. 1996, 100 (7), 2581-2587. https://doi.org/10.1021/jp9518623.
- (25)   Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H. J. C. Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes. J. Comput. Phys. 1977, 23 (3), 327-341. https://doi.org/10.1016/0021-9991(77)90098-5.
- (26)   Andrews, J. E.; Sinha, S.; Chung, P. W.; Das, S. Wetting Dynamics of a Water Nanodrop on Graphene. Phys. Chem. Chem. Phys. 2016, 18 (34), 23482-23493. https://doi.org/10.1039/C6CP01936F.
- (27)   Afshar, A.; Zhong, J.; Thompson, D. S.; Meng, D. Nanoscale Ice Fracture by Molecular Dynamics Simulations. ArXiv190709018 Cond-Mat Physicsphysics 2019.
- (28)   Wang, S.; Zhang, Y.; Abidi, N.; Cabrales, L. Wettability and Surface Free Energy of Graphene Films. Langmuir 2009, 25 (18), 11078-11081. https://doi.org/10.1021/la901402f.
- (29)   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 (6), 3798-3802. https://doi.org/10.1021/la100231u.
- (30)   Taherian, F.; Marcon, V.; van der Vegt, N. F. A.; Leroy, F. What Is the Contact Angle of Water on Graphene? Langmuir 2013, 29 (5), 1457-1465. https://doi.org/10.1021/la304645w.
- (31)   Ambrosia, M. S.; Ha, M. Y. A Molecular Dynamics Study of Wenzel State Water Droplets on Anisotropic Surfaces. Comput. Fluids 2018, 163, 1-6. https://doi.org/10.1016/j.compfluid.2017.12.013.
- (32)   Thompson, D. S.; Meng, D.; Afshar, A.; Bassou, R.; Zong, J.; Bonaccurso, E.; Laroche, A.; Vercillo, V. Initial Development of a Model to Predict Impact Ice Adhesion Stress. In 2018 Atmospheric and Space Environments Conference; American Institute of Aeronautics and Astronautics. https://doi.org/10.2514/6.2018-3344.
Amir Afshar Dong Meng
Wetting behavior of supercooled water droplets impinging on nanostructured graphite surface: A molecular dynamics study
In this work, we studied the wetting behavior of impinged and stagnant
supercooled water nanodroplet at the atomistic scale using molecular dynamics
simulations. We show that water droplet represents a retraction behavior from
surface groove nanostructure after indicating maximum penetration. There is a
threshold for the value of impinging velocity, v*, below which penetration is
not observed which depends on size of the surface roughness. This indicates
grooves smaller than a certain size can stop water from penetration,
considering v~0.5M in practical situations, which is such a useful information
for design of nanostructured ice-phobic substrates. The resistance presumably
comes from interfacial tension between solid substrate and liquid water
droplet. The extent of water penetration significantly depends on droplet
impinging velocity and the size of surface roughness, for impinging velocities
higher than V*. Our result indicates the dominant role of droplet temperature
on time required droplet recede from surface groove structure rather than
droplet impinging velocity and the surface groove nanostructure.