How Do Different Liquid Metal Films Coalesce on ... - ACS Publications

Mar 21, 2018 - Droplet coalescence plays an indispensable role in diverse natural and industrial processes. However, it is still unclear that how the ...
1 downloads 3 Views 5MB Size
Subscriber access provided by UNIV OF SCIENCES PHILADELPHIA

C: Surfaces, Interfaces, Porous Materials, and Catalysis

How Do Different Liquid Metal Films Coalesce on Carbon Substrates? Xingfan Zhang, Tao Li, Yifan Li, Yunrui Duan, and Hui Li J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b12756 • Publication Date (Web): 21 Mar 2018 Downloaded from http://pubs.acs.org on March 27, 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 40 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

How do Different Liquid Metal Films Coalesce on Carbon Substrates? Xingfan Zhang, Tao Li, Yifan Li, Yunrui Duan, and Hui Li* Key Laboratory for Liquid-Solid Structural Evolution and Processing of Materials, Ministry of Education, Shandong University, Jinan 250061, People’s Republic of China * Corresponding author: [email protected]

Abstract

Droplet coalescence plays an indispensable role in diverse natural and industrial processes. However, it is still unclear that how the wettability and microstructure of the substrates influences the coalescence. Molecular dynamics simulations were used to investigate the coalescence behaviors of liquid metal films on carbon substrates. Two different coalescence regimes, climbing-coalescence and coating-coalescence were respectively observed in Cu-Al and Cu-Ag systems on graphene. The balance of metal-substrate and metal-metal interactions was regard as the key factor to determine the coalescence regimes. Additionally, the effects of surface roughness and atomic arrangement of substrates on the movements of droplets were also analyzed on other carbon materials. Our findings offer a thorough understanding of the relationship between coalescence and wettability, and could have important implications in self-assembly, 3D printing, and microfluidic devices.

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

1.

Introduction

As one of the most common phenomena in nature, droplet coalescence has attracted our interests for a rather long time. Coalescence also plays an indispensable role in industrial production, such as petroleum industry1, self-assembly2, 3D printing technology3, and microfluidic devices4. Usually, during the coalescence process, a liquid bridge forms between the two droplets, and then the merged droplet transforms from ellipsoid to sphere5. The change in topology during coalescence is mainly controlled by three forces: viscous force, inertia, and surface tension. At first, coalescence was thought to have just two dynamical regimes: an inertial regime for low-viscosity or inviscid droplets, and a viscous regime for viscous droplets6-9. But during the past few years, more and more regimes of coalescence have been discovered in other fluids. For example, an abnormal crossover from the viscous regime to the inertial regime was observed by Paulsen et al.10 in salt water. They also discovered a third regime called “inertially limited viscous regime”11 in silicone oil, which is different from the two classic regimes. These new-found regimes indicate that the mechanism of droplet coalescence still needs to be explored. Coalescence of liquid metals plays a crucial role in metallurgy, especially with the rapid development of 3D printing technology, liquid metals have been applied to build free standing microstructures12, which requires in-depth research of coalescence to broaden the applications in this field. Nevertheless, experiments mainly focus on the structure and properties of the coalescent products, details of the coalescence mechanism are still limited. Therefore, many scholars selected to use molecular

ACS Paragon Plus Environment

Page 2 of 40

Page 3 of 40 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

dynamics (MD) simulations to solve some problems, which provide an opportunity to thoroughly understand the process from the perspective of atomic motion13-20. Besides, how the substrate affect coalescence is still poorly understood. On one hand, surface microstructure of substrates has long been recognized to have a great influence on contact angles21. People have been committed to control the wettability of droplets by modifying the surface microstructure of substrates to produce superhydrophobic or superhydrophilic materials22-24. On the other hand, coalescence of droplets on substrates25 or confined in two substrates26 were found not to consist with the classic coalescence regimes using sessile drops, which could also be associated with wettability. In fact, wetting and coalescence are closely related. For example, Au, Cu or some other liquid metal films which have large surface tension can convert their surface deformation energy into kinetic energy and jump off from the non-wettable substrates, which is called dewetting27-30. Dewetting has important applications in self-assembly nanotechnology31-34, and how to control the coalescence of droplets becomes the key factor to fabricate desired shapes and patterns2,

35-36

. In our previous works, we

reported that nanopillared graphene substrates can enhanced the dewetting property of liquid Cu films, which can be used to promote the coalescence of Cu droplets37-38. In this paper, we studied the coalescence of Cu films with other liquid metal films on carbon-based substrates, which shows some novel coalescence behaviors. The Cu droplet climbed the Al droplet to achieve coalescence, while the Ag droplet coated the Cu droplet to get coalescence due to the wettability difference on carbon-based

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

substrates. We probed the relationship between coalescence and wettability based on the atomic diffusion and interaction. Additionally, we also analyzed how the microstructure of the substrates affects the movement of droplets during the coalescence.

2.

Methods

MD simulations were performed using the large-scale atomic/molecular massively parallel simulator (LAMMPS) package39 to study the wettability of metal droplets on carbon-based substrates and their coalescence behaviors. At first, pure metal blocks were melted to 1500 K and then relaxed for 100 ps in the isobaric-isothermal (NPT) ensemble to obtain liquid metals. Then, liquid metal films with the same thickness (15 Å) and diameter (100 Å) were extracted from the liquid metal blocks and placed at a distance of 2.0 Å above the carbon substrates to investigate the wetting behaviors in the canonical (NVT) ensemble. Substrates were fixed during the simulations for the convenience of accurately measuring the contact angles and observing the morphological evolution of droplets. After 300 ps of simulations, the systems reached the equilibrium states, in which the contact angles and the diffusion coefficients can be obtained. In a coalescence set-up, two liquid metal films with the same size were placed together above the substrates at a distance of 2 Å at the initial stage. The Cu-Ag and Cu-Al systems were respectively used to study the different coalescence regimes. The simulations of coalescence were also carried out in the canonical ensemble, while the temperature was maintained constant at 1500 K using the

ACS Paragon Plus Environment

Page 4 of 40

Page 5 of 40 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

Nosé-Hoover thermostat40. The Newton's equation of motion was calculated using the velocity-Verlet algorithm with a time step of 1.0 fs throughout the simulations. AIREBO potential41 was used to describe the C-C interaction of carbon nanomaterials. The Cu-Ag interaction was described by embedded atom method (EAM) potential42. An angular-dependent (ADP) potential was used to calculate the Cu-Al interaction43. Additionally, we also chose the 6-12 Lennard-Jones (L-J) potentials with a well depth ε=0.010 eV and size parameter σ=3.225 Å to describe the Cu-C interaction28; a well depth ε=0.015 eV and size parameter σ=2.900 Å to describe the Ag-C interaction44; a well depth ε=0.0309 eV and size parameter σ=3.422 Å to describe the Al-C interaction45. With these parameters, the simulated contact angles and diffusion coefficients are in good agreements with the experimental ones, which guarantee the accuracy of the results. Mean square displacements (MSD) and the diffusion coefficients (D) of metals were calculated according to the Einstein diffusion law46: MSD =< |r t − r 0 | > D = lim



→  

(1) (2)

Where ri(t) is the position of atom i at time t and denotes an average overall atoms, qi is numerical constant that depends on dimensionality, qi = 6 for 3-dimensional diffusion in this work. Interatomic pair potential interaction energy used to explain the origin of different coalescence regimes is given as:

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

! $

E = 4ε  " #



! % "

# & , ( < ()

Page 6 of 40

(3)

Where ε is the well depth of the potential wall, σ is the finite distance at which the inter-particle potential is zero, and rc is the cutoff distance which is selected as 10.0 Å. The calculated metal-metal interaction energy values only include the pair potential portion of the EAM interactions, which may cause a certain margin of numerical error, but have no influence on the results.

3.

Results & discussion

The contact angles (CAs) of Al, Ag, and Cu droplets on the graphene substrate at 1500 K are shown in Fig. 1 (a) to illustrate the wettability difference of the metal droplets. At the nanoscale, the CAs of the droplets keep fluctuating so that we measured five values from the last five picoseconds and took the average value as the result. The simulated CAs of Al, Ag and Cu on graphene are 82.1°, 130.3° and 138.4°, respectively, which are in good agreement with the experimental values of 85°47, 129°48, and 140°49, demonstrating the accuracy of the simulations. At the equilibrium states, the Cu and Ag droplets with poor wettability form a sphere and slightly attach to graphene, while the Al droplet with good wettability keeps hemispherical and closely contacts with the substrate. The Al droplet has the best wettability due to the largest metal-carbon interaction energy as shown in Fig. 1 (b). Besides, we used the equilibrium states as the initial structures to calculate the time-evolution mean square displacement (MSD) and the diffusion coefficients (D) of metals on graphene. DAl, DAg, and DCu are computed as 0.854 Å2/ps, 0.433 Å2/ps and 0.445 Å2/ps, respectively.

ACS Paragon Plus Environment

Page 7 of 40 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 other works which investigated pure metal systems, DAl was calculated as 1.245 Å2/ps (1193 K)50, DAg was calculated as 0.44 Å2/ps (1450 K)51 , and DCu was measured as 0.60 Å2/ps (1533 K)52. By contrast, the calculated diffusion coefficients in this work are smaller than those in other pure metal systems, because the existence of metal-C interactions limits the diffusion of metal atoms on the substrate. Additionally, we compared the time-evolution MSD and metal-substrate interaction energy (Emetal-C) between the droplets and liquid films in Fig. 1 (b) and Fig. 1 (c). It can be seen that the atoms in the Cu and Ag films diffuse much faster than in the droplets because the large surface tension drives the films to rapidly turn into droplets. But the Al film does not show such huge difference due to the largest attraction from the substrate. This result illustrates that the substrate will hinder the atomic diffusion of droplets in different degrees depending on the magnitude of the metal-substrate interaction energy. Coalescence behaviors of two identical Al-Al, Ag-Ag, and Cu-Cu films on graphene are shown in Fig. 1 (d) to compare with the subsequent different metal systems. Different colors of atoms are used to distinguish the atoms in the two droplets. It can be seen that all of the coalescence processes are similar which have three stages: contraction, adjoining, and horizontal interdiffusion. Among these three systems, the coalescence rate of two Al droplets is the fastest, which can be seen from the growth rates of the liquid bridges and the final mixture uniformity of atoms, indicating that the droplets which have better wettability coalesce faster on substrates. These results indicate that the changes in metal types in the coalescence of two identical metal films

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

will cause the differences in the atomic diffusion rate, coalescence speed and the contact angle on the substrate, but will not change the coalescence regimes. However, the coalescence behaviors of different metal films are quite different. Figure 2 shows the coalescence of Cu-Al films and Cu-Ag films on graphene. From Fig. 2 (a), it can be seen that the coalescence of Cu and Al films follows the “climbing-coalescence” regime with three stages: contraction, climbing, and vertical interdiffusion. At first, both kinds of films begin to contract under the surface tension. At 30 ps, the Cu film turns into a sphere and slightly attaches to the substrate, while the Al film changes into a hemisphere and closely contacts with the substrate due to their wettability difference. Next, the Cu droplet gradually moves to the top of the Al droplet along its surface, like mountain climbing. Finally, the Cu and Al atoms interdiffuse along the vertical direction to complete the coalescence. Fig. 2 (b) shows the coalescence process of Cu and Ag films. Different from the Cu-Al system, the Cu and Ag droplets display a coating-coalescence regime which includes four stages: contraction, adjoining, coating, and radial interdiffusion. It is noteworthy that the coating and radial interdiffusion stages do not exist in the coalescence of two identical droplets. After the adjoining stage, the identical droplets just interdiffuse horizontally as shown in Fig. 1 (d). But in the Cu-Ag system, after 70 ps, like the phagocytosis process53 in biology, the Ag droplet devours the Cu droplet step by step. At about 200 ps, an Ag spherical shell has formed to fully wrap the Cu droplet. Finally, the two types of atoms interdiffuse along the radial direction and complete their coalescence.

ACS Paragon Plus Environment

Page 8 of 40

Page 9 of 40 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 interatomic interaction energy graphs in Fig. 3 can help us further understand the origin of the coalescence regimes. The calculated interatomic interaction energy describes the sum of all the pair potential interaction energy between the two types of atoms. In these graphs, a bigger negative value of interaction energy indicates the existence of stronger attraction between the two components. Higher Emetal-C means the attraction from the substrates is higher, and lower Emetal-metal means the interatomic bonding force is lower, which both leads to smaller contact angle of droplet on the substrates. In the Cu-Al system, comparing EAl-C and ECu-C in Fig. 3 (a), it can be seen that the substrate has much larger attraction to the Al atoms than the Cu atoms. As a result, the Al droplet keeps attaching to the substrate with small contact angle, while the Cu droplet keeps spherical and slightly attach to the substrate with large contact angle. During their coalescence, the Al atoms easily intrude into the interface of Cu-C and replace the Cu atoms to contact with the substrate due to the relatively higher Al-C interaction energy. After the contraction stage, the contact area between the Cu droplet and the substrate gets smaller and the Cu-Al interaction gets larger, the attraction from the Al droplet exceeds that from the substrate soon, causing the “mountain climbing” behavior of the Cu droplet. Accordingly, we can conclude that the difference in metal-substrate interaction energy determines whether the climbing phenomenon will occur. In the Cu-Ag system, ECu-C and EAg-C are both low as shown in Fig. 3 (b), so the two droplets both have large contact angles on graphene and the climbing phenomenon does not occur. However, ECu-Cu almost doubles EAg-Ag as shown in Fig. 3 (c), which is

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

the key factor of the coating-coalescence. The relatively higher ECu-Cu indicates the bonding force between the Cu atoms is higher. As a result, the Cu atoms tightly assemble as a spherical droplet, which makes it difficult for Ag atoms to embed into the Cu droplet but easy to move around the surface of the Cu droplet to get coalescence. The movement of metal atoms during the coating stage can be observed in Fig. 3 (d)-(f) in which slices of atomic velocity vectors are shown. These 5-angstrom-thick slices were cut along the directions perpendicular to the Cu-Ag interface. At 70 ps, the Cu atoms at the center of the droplet rush toward Ag under the inertial force, facilitating the coating behavior of the Ag atoms. From 70 ps to 200 ps, the Ag atoms gradually move around the surface of the Cu droplet and form a liquid shell to coat the Cu droplet. Comparing the length of the arrows, it can be seen that the atoms at the surface of the droplet and near the Cu-Ag interface diffuse faster than other atoms. It is concluded that the difference in metal-metal interaction energy determines whether the coating phenomenon will occur. When the two droplets both have low Emetal-C and the effect of the substrate can be ignored, according to the difference in metal-metal interaction energy, there are three kinds of metal droplet coalescence. The Cu-Ag system describes the first case of large difference in metal-metal interaction energy, in which the Ag droplet completely coats the Cu droplet to achieve coalescence. On the contrary, if the two metals have almost the same metal-metal interaction energy, the coating phenomenon will not occur, which can refer to the coalescence of two identical metal films, such as the coalescence of Cu-Cu, Ag-Ag,

ACS Paragon Plus Environment

Page 10 of 40

Page 11 of 40 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

and Al-Al films. In this condition, the coalescence process only has three stages: contraction, adjoining, and horizontal interdiffusion. Finally, if the difference between EA-A and EB-B is small, it is expected to be an intermediate state between the non-coating coalescence and the complete-coating coalescence. For example, if EA-A is slightly higher than EB-B, droplet B will partly coat droplet A during their coalescence. This case may need further verification in other metal systems. However, these conclusions are only suitable for the systems which have negligible droplet-substrate interactions. When the substrate has larger attraction to the droplets, the coalescence behaviors can be different. For example, in the Cu-Al system, the high EAl-C causes the coating phenomenon to be not obvious. We selected another EAM potential54 to compare ECu-Cu and EAl-Al at the initial stage because the used many-body ADP potential could not calculate the pair potential interaction energy. ECu-Cu is calculated as about -8000 eV and EAl-Al is about -2400 eV, indicating that the interatomic bonding force of the Cu atoms is much larger than Al. According to the above conclusions, the Al droplet is supposed to fully coat the Cu droplet. However, the coating phenomenon is covered up by the climbing-coalescence, mainly because of the high EAl-C. As shown in Fig. 2 (a), during the climbing-coalescence, the Cu droplet keeps spherical before the 150 ps, and the Al droplet only partly coats the Cu droplet from the bottom. Due to the large attraction from the substrate, the Al droplet closely attaches to the substrate, preventing it from further wrapping the Cu droplet, so the coating phenomenon is not obvious. As a result, under the balance of Emetal-metal and Emetal-C, the Cu and Al droplet shows the climbing-coalescence from an overall

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

Page 12 of 40

perspective. This case illustrates that although the difference between EA-A and EB-B is the key factor of the coating-coalescence, the droplet-substrate interactions may greatly affect the coalescence regimes, which guides us to further investigated the effect of the substrate microstructure in the Cu-Al system. After coating-coalescence, the stabilized EAg-C and ECu-C are -37.32 eV and -21.45 eV in the Cu-Ag system as shown in Fig. 3 (b), which vary from the results in Fig. 1 (b) calculated from the single droplet systems (EAg-C and ECu-C are -22.81 eV and -37.67 eV). The relative magnitude of EAg-C and ECu-C have changed because the Ag atoms partially

replace

the

Cu

atoms

to

contact

with

the

substrate

during

coating-coalescence, resulting in the increased EAg-C and decreased ECu-C. Additionally, the atomic diffusion was also found to be affected by the coalescence behaviors. The Cu atoms diffuse faster than Al atoms during coalescence as shown in Fig. 4 (a), which is similar to the single-film systems in Fig. 1 (c). But in the Cu-Ag system in Fig. 4 (b), the relatively larger MSD of Ag than Cu indicates that the D of Ag atoms is higher than that of Cu, which is contrary to the results of single-film systems. This difference can also be ascribe to the coating behavior of the Ag atoms, which significantly increases their diffusion rate during the coating-coalescence. Next, we investigated the coalescence on rough surfaces: horizontally-placed carbon nanotubes (HCNT) substrate, vertically-placed carbon nanotubes (VCNT) substrate, in which the carbon nanotubes are placed in parallel and perpendicular to the coalescence direction, and capped-carbon-nanotube-pillared graphene (PG) substrate, which are shown in Fig. 5 (a). Firstly, we selected the Cu-Al system to illustrate how

ACS Paragon Plus Environment

Page 13 of 40 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 substrates affect the coalescence of droplets with different wettability. Snapshots in Fig. 5 (b) show that the climbing-coalescence regime still exists under these three conditions, but the speed of coalescence and the motion of droplets in X-axis are different. It takes the shortest time to coalesce on the HCNT substrate because the grooves on the substrate are parallel to the coalescence direction, which have almost no barriers for the droplet movement. The arrangement of the carbon nanotubes side by side forms some periodic peaks and valleys along the X-axis of the VCNT substrate, which make the droplet move difficultly. This difference illustrates that on anisotropic surfaces, the location of films will affect the coalescence speed. On the PG substrate, some Al atoms penetrate into the gaps of the rough surface and hardly move just like being pinned. Such pinning effect caused by rough surfaces is remarkable for the droplet motion at the nanoscale. From the center-of-mass displacement (CMD) in Fig. 5 (b), the Al droplet is pinned so tightly that it cannot move freely on the PG and VCNT substrates. But the Cu droplet is not that case as shown in Fig. 5 (c). On the contrary, the rougher the surface is, the faster the Cu droplet moves. Due to the large surface tension and low metal-substrate interaction energy, the Cu droplet keeps spherical instead of penetrating into the gaps of the rough surfaces. The pillared surfaces decrease the contact area between Cu and C and further reduce the attraction from the substrate, thus the Cu droplet can move more freely without the restriction from substrates. By comparing the CMD curves in Fig. 5 (c)-(d), the CMD in X-axis of the Cu droplet on the three substrates follows the order PG > VCNT > HCNT, while for the Al droplet follows the converse order HCNT >

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

VCNT > PG. Therefore, it is concluded that rough surfaces can induce the droplets with good wettability to penetrate into the gaps of the surfaces due to the strong metal-substrate interaction, which will pin the droplets and hinder the coalescence. On the contrary, for the droplets with poor wettability, the metal-substrate interaction is too weak to overcome the surface tension of the droplet. Rough surfaces can reduce the contact area with the droplets instead, and further facilitate the coalescence. We also explored how the atomic arrangement affects the coalescence on smooth substrates by simulating the coalescence behaviors of the Cu and Al films on the (1 0 0), (1 1 0), and (1 1 1) surfaces of diamond (DM). As shown in Fig. 6 (a), the coalescence on DM (1 0 0) is distinctly different from other conditions that the coalescent droplet is closer to the right side of the substrate. CMD graphs in Fig. 6 (b)-(c) show that on DM (1 1 0) and DM (1 1 1), the two droplets move toward each other to achieve coalescence. But on DM (1 0 0), the Al droplet hardly moves along the X-axis and the coalescence is accomplished solely by the movement of the Cu droplet. We found the reason from the atomic arrangement of the substrates. The surface packing densities of DM (1 0 0), DM (1 1 0) and DM (1 1 1) are 15.72, 22.37, and 18.27 atoms/nm2, respectively, demonstrating that there are fewest carbon atoms to contact with the droplets on DM (1 0 0). It can also be proved in Fig. 6 (d) that the interaction energy of Al-C on DM (1 0 0) is the lowest (the interaction energy on DM (1 1 1) is the highest because the top two atomic layers are very close, which evidently increases the interaction). During the contraction stage of coalescence, the Al droplet moves upward for a larger distance on DM (1 0 0) than other conditions as

ACS Paragon Plus Environment

Page 14 of 40

Page 15 of 40 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

shown in Fig. 6 (e). Additionally, it was also found that the single Al droplet on DM (1 0 0) shows the worst wettability that it has the largest contact angle and the highest altitude as shown in Fig. 6 (f). As a result, it is difficult for the Cu droplet to climb over the higher Al droplet, and for the higher Al droplet to intrude into the bottom of the Cu droplet. Therefore, the Cu droplet not only climbs, but also pushes the higher Al droplet to move along the opposite direction. Finally, the two droplets move together along the X-axis and the compound droplet stays near the right side of the DM (1 0 0) substrate. But on DM (1 1 0) and DM (1 1 1), due to the better wettability caused by the lower Al-C interaction energy, the Al droplet is shorter in height, so that the pushing force becomes evidently smaller. As a consequence, the two droplets move toward each other to achieve coalescence. The result illustrates that the atomic arrangement of the smooth substrates will affect the coalescence behaviors of droplets by affecting their wettability. In the Cu-Ag system, the substrate microstructure has much less influence on the droplet coalescence behaviors than in the Cu-Al system. Fig. 7 (a) and Fig. 7 (b) shows the CMD of the Cu and Ag droplets in X-axis during coalescence, respectively, in which solid lines represent the cases on rough surfaces (HCNT, VCNT and PG), and dash lines represent the cases on smooth surfaces (DM (1 0 0), DM (1 1 0), and DM (1 1 1)). It can be seen that the CMD curves do not show evident difference whether on the rough or the smooth substrates, because the low ECu-C and EAg-C significantly reduce the contact area between the droplets and the substrates. These results indicate that the changes in the microstructure of the substrates have little

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

influence on the coalescence rates of two droplets with low metal-substrate interaction energy. Contact angles (CAs) of the coalescent droplets on different substrates were measured from 200 ps to 1000 ps. In the Cu-Al system, the CAs rise gradually before 700 ps because the Cu atoms diffuse downward as shown in Fig. 8 (a). The CAs show very little change after 700 ps, and the average CAs of the Cu-Al droplets from 700 ps to 1000 ps on G, HCNT, VCNT, and PG are 88.6°, 93.3°, 91.4°, and 93.6°, respectively. But in the Cu-Ag system, as shown in Fig. 8 (b), the CAs become stable after 500 ps, earlier than those in the Cu-Al system, indicating that the Cu and Ag droplets coalesce faster. The average CAs of the Cu-Ag droplets from 700 ps to 1000 ps on G, HCNT, VCNT, and PG are 128.7°, 132.1°, 133.3°, and 134.6°, respectively. In each system, the CA of droplet on the G substrate is the smallest and on the PG substrate is the largest, indicating rough substrates can increase the CAs of metal droplets. Furthermore, compared to the CAs of the pure metal droplets (The CAs of Al, Ag and Cu are 82.1°, 130.3°, and 138.4°, respectively), the CA of the alloy droplet is closer to the lower value of the two components.

4.

Conclusion

In summary, molecular dynamics simulations were performed to study the coalescence behaviors of liquid metal films on carbon substrates. Two different coalescence regimes, climbing-coalescence and coating-coalescence were observed on graphene. The Cu droplet climbs to the top of the Al droplet to achieve coalescence

ACS Paragon Plus Environment

Page 16 of 40

Page 17 of 40 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 Cu-Al system, while the Ag atoms coats the Cu droplet to get coalescence in the Cu-Ag system. The balance of metal-substrate and metal-metal interactions was regard as the key factor to determine the coalescence regimes. It is concluded that the difference in Emetal-C between the two metals determines whether the climbing phenomenon will occur, and the difference in Emetal-metal determines whether the coating phenomenon will occur. Additionally, we further studied the effect of surface microstructure of the substrates on coalescence. By comparing the coalescence behaviors of Cu and Al films on rough substrates, it is concluded that rough surfaces will pin the good-wettability droplets and restrict their movements, but make the poor-wettability droplets move more freely during coalescence. The atomic arrangement of the smooth substrates was also found to affect the coalescence behaviors of droplets by affecting their wettability. These results offer an applicable method from the perspective of interatomic interaction to investigate droplet coalescence, which could be extended to other systems which have non-bonded liquid-solid interactions, and could have certain guiding significance for new applications in self-assembly, 3D printing, and microfluidic devices.

Acknowledgements

The authors would like to acknowledge the support from the National Natural Science Foundation of China (Grant No.51671114). This work is also supported by the Special Funding in the Project of the Taishan Scholar Construction Engineering and National Key Research Program of China (Grant No. 2016YFB0300501).

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.

He, L.; Lin, F.; Li, X.; Sui, H.; Xu, Z. Interfacial Sciences in Unconventional

Petroleum Production: From Fundamentals to Applications. Chem. Soc. Rev. 2015, 44, 5446-5494. 2.

Timonen, J. V.; Latikka, M.; Leibler, L.; Ras, R. H.; Ikkala, O. Switchable Static

and Dynamic Self-Assembly of Magnetic Droplets on Superhydrophobic Surfaces. Science 2013, 341, 253-257. 3.

Martin, J. H.; Yahata, B. D.; Hundley, J. M.; Mayer, J. A.; Schaedler, T. A.;

Pollock, T. M. 3d Printing of High-Strength Aluminium Alloys. Nature 2017, 549, 365-369. 4.

Jahn, A.; Stavis, S. M.; Hong, J. S.; Vreeland, W. N.; DeVoe, D. L.; Gaitan, M.

Microfluidic Mixing and the Formation of Nanoscale Lipid Vesicles. ACS Nano 2010, 4, 2077-2087. 5.

Ristenpart, W.; McCalla, P.; Roy, R.; Stone, H. Coalescence of Spreading

Droplets on a Wettable Substrate. Phys. Rev. Lett. 2006, 97, 064501. 6.

Eggers, J.; Lister, J. R.; Stone, H. A. Coalescence of Liquid Drops. J. Fluid Mech.

1999, 401, 293-310. 7.

Duchemin, L.; Eggers, J.; Josserand, C. Inviscid Coalescence of Drops. J. Fluid

Mech. 2003, 487, 167-178. 8.

Yao, W.; Maris, H.; Pennington, P.; Seidel, G. Coalescence of Viscous Liquid

Drops. Phys. Rev. E 2005, 71, 016309. 9.

Case, S. C.; Nagel, S. R. Coalescence in Low-Viscosity Liquids. Phys. Rev. Lett.

ACS Paragon Plus Environment

Page 18 of 40

Page 19 of 40 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

2008, 100, 084503. 10. Paulsen, J. D.; Burton, J. C.; Nagel, S. R. Viscous to Inertial Crossover in Liquid Drop Coalescence. Phys. Rev. Lett. 2011, 106, 114501. 11. Paulsen, J. D.; Burton, J. C.; Nagel, S. R.; Appathurai, S.; Harris, M. T.; Basaran, O. A. The Inexorable Resistance of Inertia Determines the Initial Regime of Drop Coalescence. Proc. Natl. Acad. Sci. U. S. A. 2012, 109, 6857-6861. 12. Ladd, C.; So, J. H.; Muth, J.; Dickey, M. D. 3d Printing of Free Standing Liquid Metal Microstructures. Adv. Mater. 2013, 25, 5081-5085. 13. Ding, F.; Rosén, A.; Bolton, K. Size Dependence of the Coalescence and Melting of Iron Clusters: A Molecular-Dynamics Study. Phys. Rev. B 2004, 70, 075416. 14. Liao, M.-L.; Ju, S.-P.; Yang, S.-H. Coalescence Behavior of Water Nanoclusters: Temperature and Size Effects. J. Phys. Chem. C 2007, 111, 6927-6932. 15. Zhang, X.-Q.; Iype, E.; Nedea, S. V.; Jansen, A. P.; Szyja, B. M.; Hensen, E. J.; van Santen, R. A. Site Stability on Cobalt Nanoparticles: A Molecular Dynamics Reaxff Reactive Force Field Study. J. Phys. Chem. C 2014, 118, 6882-6886. 16. Zhang, X.-Q.; van Santen, R. A.; Hensen, E. J. Carbon-Induced Surface Transformations of Cobalt. ACS Catal. 2014, 5, 596-601. 17. Wang, B.-B.; Wang, X.-D.; Yan, W.-M.; Wang, T.-H. Molecular Dynamics Simulations on Coalescence and Non-Coalescence of Conducting Droplets. Langmuir 2015, 31, 7457-7462. 18. Wang, J.; Chen, S.; Cui, K.; Li, D.; Chen, D. Approach and Coalescence of Gold Nanoparticles Driven by Surface Thermodynamic Fluctuations and Atomic

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

Interaction Forces. ACS Nano 2016, 10, 2893-2902. 19. Li, T.; Wu, W.; Li, H. Coalescence Behavior of Liquid Immiscible Metal Drops in Two-Wall Confinement. Phys. Chem. Chem. Phys. 2016, 18, 27500-27506. 20. Faccin, G. M.; San-Miguel, M. A.; Andrés, J.; Longo, E.; da Silva, E. Z. Computational Modeling for the Ag Nanoparticle Coalescence Process: A Case of Surface Plasmon Resonance. J. Phys. Chem. C 2017, 121, 7030-7036. 21. Wenzel, R. N. Surface Roughness and Contact Angle. J. Phys. Chem. 1949, 53, 1466-1467. 22. Lee, Y.; Park, S. H.; Kim, K. B.; Lee, J. K. Fabrication of Hierarchical Structures on a Polymer Surface to Mimic Natural Superhydrophobic Surfaces. Adv. Mater. 2007, 19, 2330-2335. 23. Rafiee, J.; Rafiee, M. A.; Yu, Z. Z.; Koratkar, N. Superhydrophobic to Superhydrophilic Wetting Control in Graphene Films. Adv. Mater. 2010, 22, 2151-2154. 24. Mulroe, M. D.; Srijanto, B. R.; Ahmadi, S. F.; Collier, C. P.; Boreyko, J. B. Tuning Superhydrophobic Nanostructures to Enhance Jumping-Droplet Condensation. ACS Nano 2017, 11, 8499-8510. 25. Eddi, A.; Winkels, K. G.; Snoeijer, J. H. Influence of Droplet Geometry on the Coalescence of Low Viscosity Drops. Phys. Rev. Lett. 2013, 111, 144502. 26. Yokota, M.; Okumura, K. Dimensional Crossover in the Coalescence Dynamics of Viscous Drops Confined in between Two Plates. Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 6395-6398.

ACS Paragon Plus Environment

Page 20 of 40

Page 21 of 40 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

27. Habenicht, A.; Olapinski, M.; Burmeister, F.; Leiderer, P.; Boneberg, J. Jumping Nanodroplets. Science 2005, 309, 2043-2045. 28. Fuentes-Cabrera, M.; Rhodes, B. H.; Baskes, M. I.; Terrones, H.; Fowlkes, J. D.; Simpson, M. L.; Rack, P. D. Controlling the Velocity of Jumping Nanodroplets Via Their Initial Shape and Temperature. ACS Nano 2011, 5, 7130-7136. 29. Schebarchov, D.; Lefevre, B.; Somerville, W. R.; Hendy, S. C. Filling a Nanoporous Substrate by Dewetting of Thin Films. Nanoscale 2013, 5, 1949-1954. 30. Afkhami, S.; Kondic, L. Numerical Simulation of Ejected Molten Metal Nanoparticles Liquified by Laser Irradiation: Interplay of Geometry and Dewetting. Phys. Rev. Lett. 2013, 111, 034501. 31. Vanmaekelbergh, D. Nanoelectronics: From Droplets to Devices. Nat. Nanotechnol. 2009, 4, 475-476. 32. Fowlkes, J. D.; Kondic, L.; Diez, J.; Wu, Y.; Rack, P. D. Self-Assembly Versus Directed Assembly of Nanoparticles Via Pulsed Laser Induced Dewetting of Patterned Metal Films. Nano Lett. 2011, 11, 2478-2485. 33. Tesler, A. B.; Maoz, B. M.; Feldman, Y.; Vaskevich, A.; Rubinstein, I. Solid-State Thermal Dewetting of Just-Percolated Gold Films Evaporated on Glass: Development of the Morphology and Optical Properties. J. Phys. Chem. C 2013, 117, 11337-11346. 34. Dieter't Mannetje, S. G.; Lagraauw, R.; Otten, S.; Pit, A.; Berendsen, C.; Zeegers, J.; Van Den Ende, D.; Mugele, F. Trapping of Drops by Wetting Defects. Nat. Commun. 2014, 5, 3559. 35. Huang, J.; Kim, F.; Tao, A. R.; Connor, S.; Yang, P. Spontaneous Formation of

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

Nanoparticle Stripe Patterns through Dewetting. Nat. Mater. 2005, 4, 896. 36. Fuentes-Cabrera, M.; Rhodes, B. H.; Fowlkes, J. D.; López-Benzanilla, A.; Terrones, H.; Simpson, M. L.; Rack, P. D. Molecular Dynamics Study of the Dewetting of Copper on Graphite and Graphene: Implications for Nanoscale Self-Assembly. Phys. Rev. E 2011, 83, 041603. 37. Li, X.; He, Y.; Wang, Y.; Dong, J.; Li, H. Dewetting Properties of Metallic Liquid Film on Nanopillared Graphene. Sci. Rep. 2014, 4, 3938. 38. Li, X.; Ren, H.; Wu, W.; Li, H.; Wang, L.; He, Y.; Wang, J.; Zhou, Y. Wettability and Coalescence of Cu Droplets Subjected to Two-Wall Confinement. Sci. Rep. 2015, 5, 15190. 39. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1-19. 40. Hoover, W. G. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A 1985, 31, 1695. 41. Stuart, S. J.; Tutein, A. B.; Harrison, J. A. A Reactive Potential for Hydrocarbons with Intermolecular Interactions. J. Chem. Phys. 2000, 112, 6472-6486. 42. Williams, P. L.; Mishin, Y.; Hamilton, J. C. An Embedded-Atom Potential for the Cu–Ag System. Modell. Simul. Mater. Sci. Eng. 2006, 14, 817-833. 43. Apostol, F.; Mishin, Y. Interatomic Potential for the Al-Cu System. Phys. Rev. B 2011, 83, 054116. 44. Schebarchov, D.; Hendy, S. C. Uptake and Withdrawal of Droplets from Carbon Nanotubes. Nanoscale 2011, 3, 134-141.

ACS Paragon Plus Environment

Page 22 of 40

Page 23 of 40 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

45. Fang, R.; He, Y.; Zhang, K.; Li, H. Melting Behavior of Aluminum Nanowires in Carbon Nanotubes. J. Phys. Chem. C 2014, 118, 7622-7629. 46. Ould-Kaddour, F.; Levesque, D. Molecular-Dynamics Investigation of Tracer Diffusion in a Simple Liquid: Test of the Stokes-Einstein Law. Phys. Rev. E 2000, 63, 011205. 47. Léger, A.; Weber, L.; Mortensen, A. Influence of the Wetting Angle on Capillary Forces in Pressure Infiltration. Acta. Mater. 2015, 91, 57-69. 48. Weltsch, Z.; Lovas, A.; Takács, J.; Cziráki, Á.; Toth, A.; Kaptay, G. Measurement and Modelling of the Wettability of Graphite by a Silver–Tin (Ag–Sn) Liquid Alloy. Appl. Surf. Sci. 2013, 268, 52-60. 49. Naidich, Y. V.; Kolesnichenko, G. A Study of Wetting of Diamond and Graphite by Fused Metals and Alloys Vii. The Effect of Vanadium, Niobium, Manganese, Molybdenum, and Tungsten on Wetting of Graphite by Copper-Based Alloys. Powder Metall. Met. Ceram. 1968, 7, 563-565. 50. Demmel, F.; Szubrin, D.; Pilgrim, W.-C.; Morkel, C. Diffusion in Liquid Aluminium Probed by Quasielastic Neutron Scattering. Phys. Rev. B 2011, 84, 014307. 51. Hoyt, J.; Asta, M. Atomistic Computation of Liquid Diffusivity, Solid-Liquid Interfacial Free Energy, and Kinetic Coefficient in Au and Ag. Phys. Rev. B 2002, 65, 214106. 52. Pasquarello, A.; Laasonen, K.; Car, R.; Lee, C.; Vanderbilt, D. Ab Initio Molecular Dynamics for D-Electron Systems: Liquid Copper at 1500 K. Phys. Rev.

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

Lett. 1992, 69, 1982. 53. Champion, J. A.; Mitragotri, S. Role of Target Geometry in Phagocytosis. Proc. Natl. Acad. Sci. U. S. A. 2006, 103, 4930-4934. 54. Cai, J.; Ye, Y. Simple Analytical Embedded-Atom-Potential Model Including a Long-Range Force for Fcc Metals and Their Alloys. Phys. Rev. B 1996, 54, 8398.

Figure

ACS Paragon Plus Environment

Page 24 of 40

Page 25 of 40 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 1: (a) Wetting behaviors of liquid metal droplets on graphene. The contact angles of Al, Ag, and Cu are measured as 82.1°, 130.3°, and 138.4°, respectively. (b) Interaction energy between metal and the substrate. Solid lines represent the initial structures are metal droplets at the equilibrium states on graphene, and dash lines represent the initial structures are liquid metal films. (c) Time-evolution mean square displacement (MSD) of Al, Ag and Cu on graphene. (d) Coalescence of two identical

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

Al, Ag or Cu films on graphene.

Figure 2: Two different coalescence regimes on graphene. (a) “Climbing-coalescence” of Cu and Al films. The Al droplet climbs to the top of the Cu droplet and they achieve their coalescence by vertical interdiffusion. (b) “Coating-coalescence” of Cu and Ag films. The Ag droplet coats the Cu droplet step by step and they finally

ACS Paragon Plus Environment

Page 26 of 40

Page 27 of 40 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

interdiffuse along the radial direction to complete their coalescence.

Figure 3: Interatomic interaction energy and atomic motion during coalescence. (a) Metal-C interaction energy in the Cu-Al system. (b) Metal-C interaction energy in the Cu-Ag system. (c) Metal-metal interaction energy in the Cu-Ag system. Slices of temporal evolution of atomic velocity vectors during the coalescence of Cu and Ag droplets at the Cu-Ag interface at (d) 70 ps, (e) 150 ps, and (f) 200 ps.

Figure 4: Time-evolution mean square displacement (MSD) of metals during the

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

coalescence processes on graphene in the (a) Cu-Al system and (b) Cu-Ag system.

Figure 5: Coalescence of Cu and Al films on rough substrates. (a) The horizontally-placed carbon nanotubes (HCNT) substrate, vertically-placed carbon nanotubes (VCNT) substrate, and capped-carbon-nanotube-pillared graphene (PG) substrate are used to describe different situations of rough substrates. (b) Coalescence processes of the Cu and Al films on the three substrates. (c) Center-of-mass displacements (CMD) of the Cu droplet in X-axis during coalescence on the three

ACS Paragon Plus Environment

Page 28 of 40

Page 29 of 40 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

substrates. (d) CMD of the Al droplet in X-axis on the three substrates.

Figure 6: Coalescence of the Cu and Al films on different crystal faces of diamond. (a)

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

Snapshots of the coalescence processes on the (1 0 0), (1 1 0), and (1 1 1) faces of diamond (b) CMD of the Cu droplet in X-axis on the three substrates. (c) CMD of the Al droplet in X-axis. (d) Al-C interaction energy. (e) CMD of the Al droplet in Z-axis. (f) The wettability of the Al droplet on different faces of diamond.

Figure 7: Coalescence of Cu and Ag films on rough (HCNT, VCNT, and PG) and smooth (DM (1 0 0), DM (1 1 0), and DM (1 1 1)) substrates. (a) Center-of-mass displacements (CMD) of the Cu droplet in X-axis during coalescence on the six substrates. (b) CMD of the Ag droplet in X-axis on the six substrates.

ACS Paragon Plus Environment

Page 30 of 40

Page 31 of 40 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 8: Time-evolutional contact angles (CA) of the coalescent droplets on the G, HCNT, VCNT, and PG substrates from 200 ps to 1000 ps. (a) The coalescent Cu-Al droplet. (b) The coalescent Cu-Ag droplet.

TOC Graphic

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 1: (a) Wetting behaviors of liquid metal droplets on graphene. The contact angles of Al, Ag, and Cu are measured as 82.1°, 130.3°, and 138.4°, respectively. (b) Time-evolution mean square displacement (MSD) of Al, Ag and Cu on graphene. Solid lines represent the initial structures are metal droplets at the equilibrium states on graphene, and dash lines represent the initial structures are liquid metal films. (c) Interaction energy between metal and the substrate. (d) Coalescence of two identical Al, Ag or Cu films on graphene. 177x208mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 32 of 40

Page 33 of 40 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 2: Two different coalescence regimes on graphene. (a) “Climbing-coalescence” of Cu and Al films. The Al droplet climbs to the top of the Cu droplet and they achieve their coalescence by vertical interdiffusion. (b) “Coating-coalescence” of Cu and Ag films. The Ag droplet coats the Cu droplet step by step and they finally interdiffuse along the radial direction to complete their coalescence. 177x204mm (300 x 300 DPI)

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: Interatomic interaction energy and atomic motion during coalescence. (a) Metal-C interaction energy in the Cu-Al system. (b) Metal-C interaction energy in the Cu-Ag system. (c) Metal-metal interaction energy in the Cu-Ag system. Slices of temporal evolution of atomic velocity vectors during the coalescence of Cu and Ag droplets at the Cu-Ag interface at (d) 70 ps, (e) 150 ps, and (f) 200 ps. 177x84mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 34 of 40

Page 35 of 40 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 4: Time-evolution mean square displacement (MSD) of metals during the coalescence processes on graphene in the (a) Cu-Al system and (b) Cu-Ag system. 177x67mm (300 x 300 DPI)

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 5: Coalescence of Cu and Al films on rough substrates. (a) The horizontally-placed carbon nanotubes (HCNT) substrate, vertically-placed carbon nanotubes (VCNT) substrate, and capped-carbon-nanotubepillared graphene (PG) substrate are used to describe different situations of rough substrates. (b) Coalescence processes of the Cu and Al films on the three substrates. (c) Center-of-mass displacements (CMD) of the Cu droplet in X-axis during coalescence on the three substrates. (d) CMD of the Al droplet in Xaxis on the three substrates. 177x178mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 36 of 40

Page 37 of 40 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 6: Coalescence of the Cu and Al films on different crystal faces of diamond. (a) Snapshots of the coalescence processes on the (1 0 0), (1 1 0), and (1 1 1) faces of diamond (b) CMD of the Cu droplet in Xaxis on the three substrates. (c) CMD of the Al droplet in X-axis. (d) Al-C interaction energy. (e) CMD of the Al droplet in Z-axis. (f) The wettability of the Al droplet on different faces of diamond. 177x240mm (300 x 300 DPI)

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 7: Coalescence of Cu and Ag films on rough (HCNT, VCNT, and PG) and smooth (DM (1 0 0), DM (1 1 0), and DM (1 1 1)) substrates. (a) Center-of-mass displacements (CMD) of the Cu droplet in X-axis during coalescence on the six substrates. (b) CMD of the Ag droplet in X-axis on the six substrates. 177x75mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 38 of 40

Page 39 of 40 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 8: Time-evolutional contact angles (CA) of the coalescent droplets on the G, HCNT, VCNT, and PG substrates from 200 ps to 1000 ps. (a) The coalescent Cu-Al droplet. (b) The coalescent Cu-Ag droplet. 177x68mm (300 x 300 DPI)

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 85x34mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 40 of 40