Three-Dimensional Structure of a Simple Liquid at a Face-Centered-Cubic (001) Solid Surface Interface

Abstract

A liquid in the vicinity of a solid-liquid interface (SLI) may exhibit complex structures. In this study, we used molecular dynamics simulations demonstrating for the first time that the liquid adjacent to the SLI can have a two-level structure in some cases: a major structure and a minor structure. Through a time-averaging process of molecular motions, we identified the type of the liquid structure by calculating positions of the maximum liquid density in three spatial dimensions, and these positions were found to distribute in many dispersed zones (called high-density zones (HDZs)). The major structure appears throughout the SLI, while the minor structure only occurs significantly within the third layer. Instead of the previously reported body-centered cubic (BCC) or face-centered-cubic (FCC) types, the major structure was found to show a body-centered tetragonal (BCT) type. The adjacent HDZs are connected by specific junctions, demonstrating that atoms diffuse along some particular high probability paths from one HDZ to another. By considering the three-dimensional liquid density distribution from the continuum point of view, more complete details of the structure and diffusive behavior of liquids in the SLI are also possible to be revealed.

Three-Dimensional Structure of a Simple Liquid at a Face-Centered-

OPEN

Luyao Bao,*, Haibao Hu,*, JunWen, Paavo Sepri & Kai Luo

In order to better understand phenomena such as wetting, segregation, crystal growth, and perhaps applied topics such as Kapitza resistance and slip18, it is important to investigate issues of fundamental importance at the molecular level. Some of these issues involve the identication of the structural order type exhibited by a liquid in contact with a at crystal surface. Will this structure be a simple body-centered cubic (BCC) or face-centered-cubic (FCC)9,10? Will the atoms in the solid-liquid interface (SLI) equally likely to go in all directions when they undergo intermittent hopping11, or are there preferred directions?

The degree of order within a liquid is higher near a solid surface than is exhibited in the short-range order in bulk liquids, which has been observed in both theoretical and experimental investigations6,1217. Recently, main results of ordered structures of a liquid in the SLI include two aspects: one is density uctuations perpendicular to the interface, called the layering order, and the other is in-plane order, which describes the order of liquid occurring parallel to the SLI. The layering order has been well studied by experiments4,1822 and atomistic simulations7,2325. These works have consistently demonstrated that the layering order is universal in the SLI, and its features are sensitive to its environment and properties of both the solid substrate and the bulk liquid.

The in-plane order has been observed experimentally1, and the results are consistent with earlier molecular dynamics (MD) simulations1. In addition to experimental revealing, research could also benet from alternative techniques, atomistic simulations, in obtaining detailed data on the atom positions in the SLI. Atomistic simulations are able to provide atom-resolved insights into phenomena occurring near the interface26. Molecular dynamics simulations15,23,27 indicate that the liquid can form several dierent structures (such as BCC or FCC) under dierent conditions, and that the structural type may dier with dierent solid substrates, or if the liquid condenses to the solid.

Some researchers have described the ordering of liquid structures near the interface being similar to that of a pure solid from the atomistic point of view10. Although this appears to work for the rst few layers, liquid-like features become dominant, as the distance to the interface increases. Other researchers6,7,9,23,27,28 have concentrated

on the liquid structure using a statistical approach. Some have investigated the time-averaged atomistic positions

R. China. *These authors

R

P

SCIENTIFIC REPORTS

1

www.nature.com/scientificreports/

Figure 1. Front view snapshot of the molecular dynamics simulation. The red points are solid atoms in the lower wall and the blue points are liquid atoms. The green line represents a structureless Lennard-Jones (LJ) at upper wall, which constrains the liquid atoms.

in two-dimensional planes parallel to the solid surface. Such two-dimensional representations may miss important features of liquid structures between adjacent planes.

In the present study, we extended our molecular dynamics computations to include the structure of the liquid in three dimensions. Structures and behaviors of liquid argon in contact with a platinum FCC (001) crystal surface are studied. We found that the liquid has a two-level structure, and that liquid atoms near the SLI hop along two paths of high probability. We also discuss the mechanism of the formation of liquid structure in the SLI based on the Boltzmann law.

Methods

We focused on the liquid argon conned between an FCC (001) solid platinum bottom wall and a numerical top wall, as shown in Fig.1. The x, y, z dimensions of the simulation region are (29.4, +29.4), (58.0, +10.0), (29.4, +29.4), respectively, whose units are angstroms (). The linear extent of the y direction is sufficient to accommodate the existence of some weak liquid layers in the central part of the simulation box, indicating that the liquid is conned. These system dimensions were chosen to facilitate comparison to the cases where the liquid is completely unconned.

The interaction between atoms in the simulations is taken to be the Lennard-Jones 12-6 potential function (LJ/12-6 potential): U(r)=4[(/r)12 (/r)6]. Here, r is the distance between two argon atoms, is the potential well depth, and is the distance of zero potential. The Lennard-Jones parameters for argon are: =0.01042 eV, = 3.405, and the mass is: m= 39.948g/mol. A cuto distance (rc=2.5) for the LJ/12-6 potential is adopted to improve the computational efficiency. The interaction between two atoms is considered to vanish when the distance between them is larger than rc. We chose the LJ/12-6 potential to describe the liquid argon because this potential has been shown to replicate well the thermodynamic properties of argon obtained experimentally.

The equilibrium temperature and bulk number density of the liquid argon are taken to be 85K and 0.01963,

respectively. A noteworthy observation in Fig.1 is that liquid atoms in the rst and second rows are instantaneously positioned within alternating cells of spacing 3.92. The view seems to show clustering of atoms in these cells, but that is merely a superposition of many atoms along the z direction.

At this point, we provide some more information concerning the thermodynamic state of the liquid argon. Johnson et al.29 improved on the Nicolas30 equation of state for LJ liquids by using new molecular dynamics simulations. He arrived at a new set of parameters for the 33-parameter modied Benedict-Webb-Rubin (MBWR) equation. According to the reported state parameters for liquid argon and the MBWR equation of Johnson et al., the thermodynamic state of the liquid argon in our case lies in the phase diagram region between the coexistence and spinodal curves31, but does not lie close to either curve. Due to the connement, the phase diagram of liquid argon will be slightly dierent from that of the bulk argon32. Therefore, it is reasonable to assume that the thermodynamic state of argon here is in the above-mentioned region of the phase diagram.

The bottom wall consists of nine layers of platinum FCC (001) crystal and contains 4050 atoms (the red points in Fig.1). Note that face-centered atoms are also shown without distinction from the corner atoms. The lattice parameter is denoted by a ( 3.92 ), and the solid is described by the potential of the embedded-atom method (EAM)33. The spacing of red atoms in Fig.1 is thus 1.96. The cuto distance, mass of solid atoms, and

SCIENTIFIC REPORTS

2

www.nature.com/scientificreports/

the tabulated values of the functions are taken from S. M. Foiles et al.33. The top wall (the green line in Fig.1) is modelled to be a at LJ/12-6 potential boundary.

Interactions between solid and liquid atoms play an important role in simulating the liquid structure at the SLI. We chose the LJ/12-6 potential to describe these interactions, as other investigators have successfully applied it to the same topic10,23,34, and even other interfacial phenomena, such as wetting3537 and nanouidics3840 using MD simulations. This potential is less satisfactory to model the liquid-solid interaction especially for experimentalists. However, our paper focus on the general trends of local order rather than attempting to extract physically meaningful numbers2. And the general trends from such simulations are important for understanding the general phenomenon2.

The depth of the LJ/12-6 potential, denoted by Solid-Liquid (SL), can be used to represent the wettability of the

liquid to the solid41. For example, SL is set equal to liquid-Liquid (LL), indicating that the solid substrate is very hydrophilic to the liquid40,42. Under this circumstance, the platinum crystal has no particular inuence other than providing an atomistic solid substrate. Dierent lattice structures of the solid substrates should generate dierent liquid structures at the SLI owing to the change in cell size of the potential eld conguration43. However, such variations are not the focus of our present work. In all cases, the distance of zero potential between platinum atoms is 2.935 , which is used to determine the potential parameter () between platinum and argon atoms using the Lorentz-Berthelot mixing rule.

An initial period of computation is completed rst in order to achieve an equilibrium state as follows. Solid atoms are arranged within the bottom wall in an FCC crystal structure (a = 3.92 ), and the liquid atoms are arranged internally in an FCC structure with a lattice parameter of 1.34a. The temperature of the liquid is then increased to 2000K to produce complete disorder in the liquid10. The positions of solid atoms are kept unchanged in this process. Aer the liquid reaches an initial equilibrium, the liquid is then cooled to 85K, and the temperature of the bottom wall is set at 85K to simulate a exible wall. The thermostat method used in all these processes is that of Nos-Hoover. Aer the system reaches this second equilibrium, the thermostat is removed to make the system into an NVE ensemble. All the simulations were carried out using the LAMMPS code44 with a timestep of 0.002ps. The three spatial dimensions (x, y, z) are dened such that y is the perpendicular distance away from the wall. Periodic boundary conditions are imposed at the extremities of both x and z directions for solid and liquid alike. Whenever an atom crosses such a boundary, it will renter into the simulation box at the opposite boundary. While the atoms move in a random thermal sense, there is no macroscopic cross-ow imposed either in the x or z directions.

The inuence of periodic boundary conditions on the thermodynamic state of system becomes insignicant in a statistical sense as the system exceeds certain spatial dimensions. In order to verify this, we have run four simulations with argon to observe the eect on the pressure at dierent temperatures in the range 85115K. The rst simulation, (a), was with the same dimensions as the simulation reported herein. The other simulations were with dimensions larger than the rst: (b) double the rst simulation in the x direction; (c) double the rst simulation in the z direction; and (d) double the rst simulation in both x and z directions. The resulting average pressures of argon diered insignicantly among these simulations. Therefore, we conclude that the sizes used in the simulation reported here are sufficiently large to model the bulk liquid argon.

The three-dimensional volume is divided into an equal mesh with spacings: x =0.392 , y =0.4 , and z= 0.392. The time-averaged density distribution, (x, y, z), is calculated at corresponding cells. To help with the analysis, (y) is also calculated by time averaging the number density of atoms in a slice of width y=0.01 parallel to the solid surface.

Results

We rst investigated the three-dimensional liquid structure at the SLI for two dierent cases, namely: (1) a high degree of incommensurability in size between solid and liquid atoms40,43, and (2) a very hydrophilic substrate.

The layering order is characterized by the decaying oscillation of the liquid density (denoted as (y) in this paper) perpendicular to wall. In this sense, (y) is considered to be a continuous function, and it represents both the solid-like and liquid-like features of the liquid itself in the SLI. We extended this computational method to three dimensions, so that the time-averaging procedure produced a probabilistic function of the liquid density distribution in space, (x, y, z). Here, the density at each location represents the total mass-centers of atoms per unit volume statistically. We are able to identify both the uid structure and atom mobility using (x, y, z). For the bulk crystal without any lattice faults, the three-dimensional uid density consists of many dispersed spheres whose centers congregate within the potential function lattice cells generated by the solid atoms oscillating around their lattice points. We can identify the structure near the SLI based on the positions of the density spheres, and we can measure the atom mobility quantitatively from the dynamic details of the spheres. For the bulk region of the liquid, (x, y, z) is uniform in space because atoms are able to diuse in all directions without bias.

We dene the boundaries of liquid layers parallel to the wall where the density distribution, (y), achieves successive minima. These planes appear predominantly blue in color as surmised from Fig.2(b,c). The planes shown in Fig.2(a) represent cuts through the central portions of the liquid layers. The rst layer is comprised of the liquid between the wall and rst minimum of (y). We examined the rst six layers and denoted them as Nos. 1~6 successively away from the wall. The liquid is clearly distributed into periodic cells of high density, which we call high-density zones (HDZs). Surprisingly, it appears that the liquid structure consists of two dierent types of HDZs rather than a single one. For the rst type, the liquid density is high and the structure extends throughout the SLI region (Fig.2). Therefore, we call these to comprise the major HDZs. For the second type, the liquid density is much lower than in major HDZs, and thus we call them the minor HDZs. The minor HDZs are shown in Fig.2(c) within the third layer signicantly, and perhaps in the h layer.

SCIENTIFIC REPORTS

3

www.nature.com/scientificreports/

Figure 2. (a) The top six planes show the two-dimensional liquid density, ((x, y0, z), parallel to the wall atthe specied locations, y0, perpendicular to the wall. These six locations correspond to the rst six maximaof (y). a is the lattice constant of lower wall. The brown solid spheres reside in the rst layer of the solid wall and the bottom plane depicts the wall-potential at the location of the second plane. We only show segments of the horizontal planes to make the density distribution clear. The density distribution is shown in the elevated planes with corresponding periodicities and amplitudes. (b,c) represent vertical cuts of liquid density obtained by intersections through the third plane in (a) at the yellow and pink solid lines, respectively. The black solid arrows in (a,b) show the high-probability hopping paths for liquid atoms within the major structure. The magenta solid arrows in (c) show the high-probability hopping paths within the minor structure. The time span used for averaging is 200ns, which is long enough to obtain a smooth distribution for (x, y, z).

Figure 3. Quantitative evidence of the BCT arrangement of the centers of the major HDZs. Nos 16 are the rst six layers successively away from the wall. The error bars are standard deviation from 210 HDZs for corresponding data.

The centers of HDZs are dened to be the positions of liquid density maxima, analogous to the density of the bulk solid mentioned earlier. We dene the liquid structure as the arrangement of the HDZs centers. The presence of two types of HDZs demonstrates that there is a two-level structure for the liquid in the SLI under the conditions of the present simulation. We named the arrangement of the major and minor HDZs the major and minor structure, respectively. Within the SLI region, the transition from high ordering at the interface to the bulk liquid further away appears to have this compound structure.

Figure2(a) indicates that the major structure of liquid is of BCC or BCT types. We conrm this by further analyzing the simulation results. There are 210 major HDZs for each layer in our simulation. First, we identify the positions (xi, yi, zi) of the centers of the major HDZs and compare the xi and zi components with the corresponding components xp_i and zp_i coming from the rst layer at the lower wall. Three dierences in location are calculated as: xi=xixp_i, zi=zizp_i and di=((xi)2+(zi)2)1/2. Then, we obtained the average and variance of these three quantities for the 210 HDZs in each horizontal layer. The results are shown in Fig.3.

SCIENTIFIC REPORTS

4

www.nature.com/scientificreports/

Figure 4. (a) Illustration of the lowest potential distance L1=21/6LL from the reference liquid atom (solid squares) and the distance L2=a/2 from the lowest potential generated by four nearest solid atoms (solid circles) to the same reference liquid atom. a is the lattice constant of the solid; The in-plane order of the rst (b1,c1) and second (b2,c2) layers: (b) SL=LL and SL= 2.954, unconned; (c) SL=7LL and SL=2.954, unconned. The time span used for averaging is 200ns which is long enough to obtain the smooth (x, y, z).

The value of d (average of the individual di) for each layer is much less than the bottom-wall lattice constant a. The corresponding red solid circles in Fig.3 demonstrate that the centers of the major HDZs line up in a vertical structure. In addition, the values of x and z (average of the xi and the zi, respectively) are nearly zero, as expected (the blue and pink solid triangles in Fig.3). The error bars show the range of variation in these measures. Another interesting observation is about the locations of the HDZ planes, yi, corresponding to the positions of the maxima of (y) as mentioned in the caption of Fig.2(a). The dierence () between the values of yi of adjacent layers increases as the layers get further from the bottom wall (the black solid squares in Fig.3). This indicates that the separation between HDZ layers increases progressively away from the wall. These observations demonstrate that the major structure of the liquid in the SLI is BCT instead of BCC or FCC under conditions of the present simulation. The HDZs of the minor structure are also aligned vertically above and below the centers of their unit cell respectively (Fig.2(c)).

The HDZs within the liquid result from the vibrational behavior of atoms similar to that of the solid, and this might be termed a solid-like feature of SLI. However, liquid atoms do not vibrate in the HDZs all the time, and will occasionally hop from one HDZ to another both horizontally and vertically. Such motion has been called the intermittent-hopping behavior11, and this could be considered a liquid-like feature of the SLI. The density distribution is proportional to the probability that liquid atoms appear at any given position. In Fig.2(ac), there are specic junctions that exist between adjacent HDZs, which demonstrate that hopping atoms follow particular high probability paths within one BCT unit cell. In the horizontal direction, atoms hop along the edges of the BCT unit cell (the black horizontal arrows in Fig.2(a)). In the vertical direction, atoms follow two paths: one is along vertices towards the center of the BCT unit cell (the inclined black arrows in Fig.2(b)); and the other is from a vertex of a minor HDZ and towards the adjoining major HDZs of the unit cell (the magenta arrows in Fig.2(c)). We call the hopping paths in the rst and second level structures the major and minor paths, respectively.

Geysermans6,34 et al. have also discussed the diusion paths and in-plane order of the liquid at the SLI using a MD simulation. They proposed that the diusion behavior of liquid atoms relaxes the mist stresses existing in the liquid layer. In our present study, we have been able to supply additional details in three dimensions concerning two types of structures and hopping paths of the liquid at the SLI.

Discussion

The density distribution in general be a direct result of the Boltzmann law which states that density is proportional to exp(U/kBT), where U and T are the potential energy and temperature at some positions, and kB is the Boltzmann constant. Atoms at the bottom wall generate low-potential wells (WLPWs) which arrange with periodicity parallel to the wall (the rst plane in Fig.2(a)). These WLPWs tend to trap liquid atoms and cause them to vibrate around the lowest points of the WLPWs, thereby resulting in the formation of the rst liquid layer (FLL). The FLL will generate a similar potential eld leading to the formation of the second liquid layer. Thus, the structure of the wall is transmitted into the liquid forming the major structure.

The depth in the sequence of low-potential wells (LPWs) gets progressively shallower away from the wall. This can be attributed to the increasingly diusive behavior of atoms evidenced by the change in connections between adjacent HDZs. In Fig.4, we show the relationship of the HDZ locations relative to two characteristic length scales: L1=21/6LL and L2=a/2. The former represents the distance to the minimum point of the LJ/12-6 potential, and the latter represents the diagonal length in the crystal lattice structure. In Fig.4(a), the solid squares represent the locations of the HDZs in the third layer, and the solid circles represent the locations of solid atoms at the wall. Therefore, at the third layer the depth of the LPWs generated by the atoms in the two nearest major HDZs (the two red solid squares in Fig.4(a)) are shallower but comparable to the depth of LPWs generated by the atoms in all four nearest major HDZs (the four solid squares in Fig.4(a)) in this simulation. Some liquid atoms may appear in the LPWs generated by the two nearest major HDZs at a lower probability than the atoms in the major HDZs. Thereby, this may explain the appearance of the second level structure of liquid in the SLI. The same argument does not apply for the rst two liquid layers shown in Fig.2 according to our simulation results.

The conguration of the FLL determines the entire structure of the SLI according to the previous analysis. Therefore, we only need to discuss the formation of the FLL. If one atom in the major HDZ tries to enter into an empty region of the WLPZ to change the conguration of the FLL, the adjacent HDZs will impose a repulsive force to prevent this from happening since L1 is larger than L2 (Fig.4(a)). The intensity of this prevention depends on the probability of atoms appearing in the four nearest major HDZs. This probability depends on the number

SCIENTIFIC REPORTS

5

www.nature.com/scientificreports/

of atoms occupying the FLL. In our simulation, the time-averaged number of atoms in the FLL is 12.4377, which means that the probability of the WLPZ being occupied by atoms is 0.05923 (divide the above by 210, the total number of WLPZs). It turns out that this probability is large enough to maintain the conguration of the FLL as shown in Fig.2(a). From the point of view of the Boltzmann law, the potential generated by the four nearest HDZs counteracts the depth of the WLPZs located at these HDZs statistically. Therefore, the probability that liquid atoms are trapped into these WLPZs is so low that the liquid atoms prefer to enter into other deeper LPZs generated by the FLL or remaining liquid layers. Thus, the conguration of the FLL is maintained.

We verify this assertion by carrying out another simulation where the top wall is 10 higher than previous cases in order to create a free liquid-gas interface. Under this condition, the time-averaged number of atoms in the FLL reduces to 11.0216 and the probability of a WLPZ occupied by atoms reduces to 0.05248. As expected, all the WLPZs are occupied and the HDZs appearing now are weaker than their four surroundings (Fig.4(b1)) owing to the oset of these WLPZs from atoms in the surrounding four strong HDZs. In addition, the HDZs in the successive layers arrange based on the rst layer (Fig.4(b2)). The structure in Fig.4(b1,b2) indicates that the conguration of the two-level structure of liquid in the SLI is not unique but changes with the conditions of the system.

Another reason for the formation of the structure shown in Fig.2 is the mobility of liquid atoms in the FLL, which also leads to similar HDZs in all layers. The solid-liquid interaction is equal to the liquid-liquid interaction. The liquid atoms in the FLL possess enough mobility to hop horizontally and vertically owing to the complete liquid state in the bulk liquid. Therefore, the LPZs each have an equal probability for occupation by atoms at any given time. Hence, the HDZs in every layer have the same structure. If we decrease the mobility of atoms in the FLL by increasing the strength of the wall potential, the HDZs in the rst layer should become dierent from each other and the in-plane order might even disappear. The results of simulations shown in Fig.4(c1,c2) are consistent with our explanations. Here, the solid-liquid interaction in case (c) is 7 times greater than that in case (b). This strength of interaction is approximately equal to the interaction between Pt and Ar atoms according to Lorentz-Berthelot mixing rule, if the interaction between Pt and Ar atoms is modeled by the LJ/12-6 potential. Under this condition, there is no periodic structure of liquid in the solid-liquid interface according to the third simulation results in our paper. As the high temperature leads to appropriate mobility of liquid atoms in the rst liquid layer, the periodic structure including the minor structure is expected to exist in physical Pt/Ar system. In other word, the minor structure predicting from our simulations will exist in the system under (1) a high degree of incommensurability in size between solid and liquid atoms, and (2) a very hydrophilic substrate (SL is approximately equal to LL). In addition, this result also demonstrates that the conguration of the FLL determines the entire structure of the SLI.

The slip phenomenon sometimes appearing in liquid-solid systems directly relates to the liquid structure at the SLI39,45. In the molecular kinetic theory of slip, slip is mainly caused by liquid atoms hopping from one equilibrium site to another, passing through a location of higher energy46,47. Also, the slip between liquid layers at the SLI makes a contribution to the apparent slip at high rates of shear47. Once we know the three-dimensional density distribution of owing liquid at the SLI, we will be able to obtain the rate of liquid atoms hopping from one site to another directly according to the molecular kinetic theory of slip47. The momentum transport between the liquid layer at the SLI is more complicated due to the unknown but signicant48 inuence of liquid structure order on viscosity in sequential layers. Despite of this, it is of great value to analyze the stress tensor by referring to the three-dimensional density distribution.

The implications for segregation have been well commented by Kaplan and Kaumann2. According to their argument, the concept of a decrease in entropy due to ordering should be taken into account for segregation at solid-liquid interfaces, where segregation and ordering both occur. To this end, the complete three-dimensional liquid structure should be helpful in describing the combined ordering-segregation phenomenon using the framework of the multilayer adsorption2.

References

1. Kauffmann, Y. et al. Quantitative analysis of layering and in-plane structural ordering at an alumina-aluminum solid-liquid interface. Acta Mater. 59, 43784386 (2011).

2. Kaplan, W. D. & Kaumann, Y. Structural order in liquids induced by interfaces with crystals. Ann. Rev. Mater. Res. 36, 148 (2006).3. Oh, S. H., Kaumann, Y., Scheu, C., Kaplan, W. D. & Ruhle, M. Ordered liquid aluminum at the interface with sapphire. Science 310, 661663 (2005).

4. Gandman, M., Kaumann, Y., Koch, C. T. & Kaplan, W. D. Direct Quantication of Ordering at a Solid-Liquid Interface Using Aberration Corrected Transmission Electron Microscopy. Phys. Rev. Lett. 110, 086106 (2013).

5. Kelton, K. F. Glass-forming alloys order at the interface. Nat. Mater. 12, 473474 (2013).6. Geysermans, P., Gorse, D. & Pontikis, V. Molecular dynamics study of the solid-liquid interface. J. Chem. Phys. 113, 63826389 (2000).

7. Men, H. & Fan, Z. Atomic ordering in liquid aluminium induced by substrates with mists. Comput. Mater. Sci. 85, 17 (2014).8. Zhang, H. L. et al. Atomic study on the ordered structure in Al melts induced by liquid/substrate interface with Ti solute. Appl. Phys. Lett. 106, 041606 (2015).

9. Hashibon, A., Adler, J., Finnis, M. W. & Kaplan, W. D. Ordering at solid-liquid interfaces between dissimilar materials. Interface Sci. 9, 175181 (2001).

10. Donnelly, S. E. et al. Ordering in a uid inert gas conned by at surfaces. Science 296, 507510 (2002).11. Skaug, M. J., Mabry, J. & Schwartz, D. K. Intermittent Molecular Hopping at the Solid-Liquid Interface. Phys. Rev. Lett. 110, 256101 (2013).

12. Baram, M., Garofalini, S. H. & Kaplan, W. D. Order in nanometer thick intergranular lms at Au-sapphire interfaces. Acta Mater. 59, 57105715 (2011).

13. Hashibon, A., Adler, J., Finnis, M. W. & Kaplan, W. D. Atomistic study of structural correlations at a liquid-solid interface. Comput. Mater. Sci. 24, 443452 (2002).

14. Broughton, J. Q., Bonissent, A. & Abraham, F. F. The fcc (111) and (100) crystalmelt interfaces: A comparison by molecular dynamics simulation. J. Chem. Phys. 74, 40294039 (1981).

15. Davidchack, R. L. & Laird, B. B. Simulation of the hard-sphere crystal-melt interface. J. Chem. Phys. 108, 94529462 (1998).16. Rull, L. F. & Toxvaerd, S. The structure and thermodynamics of a solid-fulid interface.2. J. Chem. Phys. 78, 32733278 (1983).

SCIENTIFIC REPORTS

6

www.nature.com/scientificreports/

17. Shpyrko, O. G. et al. Crystalline surface phases of the liquid Au-Si eutectic alloy. Phys. Rev. B 76, 245436 (2007).18. Borgardt, N. I., Plikat, B., Schroter, W., Seibt, M. & Wagner, T. Atomic structure of the interface between silicon (111) and amorphous germanium. Phys. Rev. B 70, 195307 (2004).

19. Tamam, L. et al. Atomic-Scale Structure of a Liquid Metal-Insulator Interface. J. Phys. Chem. Lett. 1, 10411045 (2010).20. Reedijk, M. F. et al. Structure of liquid Sn on Ge(111). Phys. Rev. B 64, 033403 (2001).21. Arai, S., Tsukimoto, S., Miyai, H. & Saka, H. High-resolution in-situ transmission electron microscopy observation of a solid-liquid interface in the Al-Si system. J. Electron. Microsc. 48, 317321 (1999).

22. Perret, E. et al. X-ray reectivity theory for determining the density prole of a liquid under nanometre connement. J. Synchrotron. Radiat. 17, 465472 (2010).

23. Geysermans, P. Molecular dynamics of the asymmetric solid-liquid interface. Mol. Phys. 103, 27172724 (2005).24. Mishin, Y., Asta, M. & Li, J. Atomistic modeling of interfaces and their impact on microstructure and properties. Acta Mater. 58, 11171151 (2010).

25. Gao, X. Z., Mueser, M. H., Kong, L. T. & Li, J. F. Atomic structure and energetics of amorphous-crystalline CuZr interfaces: a molecular dynamics study. Model. Simul. Mater. Sc. 22, 065007 (2014).

26. Guerdane, M., Teichler, H. & Nestler, B. Local Atomic Order in the Melt and Solid-Liquid Interface Eect on the Growth Kinetics in a Metallic Alloy Model. Phys. Rev. Lett. 110, 086105 (2013).

27. Palafox-Hernandez, J. P., Laird, B. B. & Asta, M. Atomistic characterization of the Cu-Pb solid-liquid interface. Acta Mater. 59,

31373144 (2011).

28. Yang, Y., Olmsted, D. L., Asta, M. & Laird, B. B. Atomistic characterization of the chemically heterogeneous Al-Pb solid-liquid interface. Acta Mater. 60, 49604971 (2012).

29. Johnson, J. K., Zollweg, J. A. & Gubbins, K. E. The lennard-jones equation of state revisited. Mol. Phys. 78, 591618 (1993).30. Nicolas, J. J., Gubbins, K. E., Streett, W. B. & Tildesley, D. J. Equation of state for the Lennard-Jones uid. Mol. Phys. 37, 14291454 (1979).

31. Cosden, I. A. & Lukes, J. R. A hybrid atomistic-continuum model for uid ow using LAMMPS and OpenFOAM. Comput. Phys. Commun. 184, 19581965 (2013).

32. Page, K. S. & Monson, P. A. Phase equilibrium in a molecular model of a uid conned in a disordered porous material. Phys. Rev. E 54, R29R32 (1996).

33. Baskes, M. I. & Daw, M. S. Embedded-atom-method functions for the fcc metals Cu, Ag, Au, Ni, Pd, Pt, and their alloys. Phys. Rev. B 33, 79837991 (1986).

34. Geysermans, P. & Pontikis, V. The structure of the solid-liquid interface: Atomic size eect. J. Phys. IV 12, 239248 (2002).35. Yuan, Q. Z. & Zhao, Y. P. Precursor Film in Dynamic Wetting, Electrowetting, and Electro-Elasto-Capillarity. Phys. Rev. Lett. 104, 246101 (2010).

36. Heine, D. R., Grest, G. S. & Webb, E. B. Surface wetting of liquid nanodroplets: Droplet-size eects. Phys. Rev. Lett. 95, 107801 (2005).

37. Coninck, J. D. & Blake, T. D. Wetting and Molecular Dynamics of Simple Liquids. Ann. Rev. Mater. Res. 38, 122 (2008).38. Thompson, P. A. & Troian, S. M. A general boundary condition for liquid ow at solid surfaces. Nature 389, 360362 (1997).39. Priezjev, N. V. Rate-dependent slip boundary conditions for simple uids. Phys. Rev. E 75, 051605 (2007).40. Yong, X. & Zhang, L. T. Slip in nanoscale shear ow: mechanisms of interfacial friction. Microuid. Nanouid. 14, 299308 (2013).41. de Gennes, P. G. Wetting-statics and dynamics. Rev. Mod. Phys. 57, 827863 (1985).42. Voronov, R. S., Papavassiliou, D. V. & Lee, L. L. Boundary slip and wetting properties of interfaces: Correlation of the contact angle with the slip length. J. Chem. Phys. 124, 204701 (2006).

43. Thompson, P. A. & Robbins, M. O. Shear-ow near solids-epitaxial order and ow boundary-conditions. Phys. Rev. A 41, 68306837 (1990).

44. Plimpton, S. Fast parallel algorithms for short-range molecular-dynamics. J. Comput. Phys. 117, 119 (1995).45. Barrat, J. L. & Bocquet, L. Inuence of wetting properties on hydrodynamic boundary conditions at a uid/solid interface. Faraday Discuss. 112, 119127 (1999).

46. Lichter, S., Martini, A., Snurr, R. Q. & Wang, Q. Liquid slip in nanoscale channels as a rate process. Phys. Rev. Lett. 98, 226001 (2007).47. Wang, F. C. & Zhao, Y. P. Slip boundary conditions based on molecular kinetic theory: The critical shear stress and the energy dissipation at the liquid-solid interface. So Matter 7, 86288634 (2011).

48. Ortiz-Young, D., Chiu, H. C., Kim, S., Voitchovsky, K. & Riedo, E. The interplay between apparent viscosity and wettability in nanoconned water. Nat. Commun. 4, 24822485 (2013).

This work was supported by the National Natural Science Foundation of China (No. 51335010, 51109178), the Fundamental Research Funds for the Central Universities (No. 3102015ZY017).

Author Contributions

H.H. conceived the project and planed the simualtions. L.B. and J.W. performed the simualtions analyzed the data. H.H., P.S. and K.L. contributed to the theoretical discussion of results. L.B. and H.H. wrote the manuscript. All authors reviewed the manuscript.

Competing nancial interests: The authors declare no competing nancial interests.

How to cite this article: Bao, L. et al. Three-Dimensional Structure of a Simple Liquid at a Face-Centered-Cubic (001) Solid Surface Interface. Sci. Rep. 6, 29786; doi: 10.1038/srep29786 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

SCIENTIFIC REPORTS

7