Several metals and many alloys are in liquid form around room temperature, e.g., mercury (Hg, −38.8 °C), francium (Fr, 8.0 °C), cesium (Sc, 28.5 °C), gallium (Ga, 29.8 °C), the eutectic mercury-based alloys, and the eutectic gallium-based alloys. If eutectic, liquid metal alloys can be in liquid form that has been used in practical applications, replacing mercury. Liquid metals have high thermal and electric conductivity and have been used to conduct heat and electricity between non-metallic and metallic surfaces. They have also been used as thermal interface materials between coolers and processors. Concerning these metals, our understanding of the liquid-vapor interface is critical for proper applications. This entry summarizes the basic features of the density distribution of liquid metal-vapor interface, which are advanced based on pseudo-potential representation and numerical simulation at the University of Chicago.
Advances in the understandings of the structures of liquid metal-vapor and liquid metal-solid interfaces have been made in the past decades . The progress has been made possible by a combination of theoretical computer simulation and experimental studies of grazing incidence x-ray diffraction and reflection [26-34], including simple liquid metals, binary alloys , and ternary alloys .
The characteristic factors of the liquid metal-vapor interface are represented by the density distribution along the normal and the pair correlation function parallel to the liquid-vapor interface. These factors can be used to predict the magnitude of segregation of solute in the liquid-vapor interfaces of alloys and the occurrence of two-dimensional crystallization in the segregated layer in the liquid-vapor interface in dilute alloy systems.
Due to the complexity of the inhomogeneous systems, there is no ideal solution to solve these systems. So far, the most successful theoretical approach is the pseudo-potential representation that is established through a multi-component, disorderly distribution of ion cores and valence electrons. Then, based on the constructed pseudo-potential system Hamiltonian, a self-consistent Monte Carlo simulation is carried out to obtain the longitudinal density distribution.
In many practical applications, local density approximations have been employed to approximate an inhomogeneous system. This is achieved by invoking the properties of a homogeneous fluid and using a local density to approximate the properties of an inhomogeneous system. Within many local density approaches, it appears that the local density with a pseudopotential approximation for ion interaction is quite satisfactory in the descriptions of the behavior of the transversal pair correlation function.
The reported liquid metal-vapor interface systems include both simple liquid metals and their alloys, e.g., alkaline metals, gallium (Ga), aluminum (Al), indium (In), thallium (Tl), mercury (Hg), tin (Sn), lead (Pb), InGa binary alloys, BiGa binary alloys, GaSn binary alloy, dilute Pb in Ga alloy, dilute Tl in Ga Alloy, dilute ternary alloy of Pb and Sn in Ga, and more.
The pseudo-potential system Hamiltonian is given by
where is the momentum of the atom-i with mass m; r is the distance for an electron to the center of the atom. R is the distance between atom-i and atom-j, is the effective pair interaction potential, The pseudo-potential is a structure independent contribution to the energy, including the electron-ion and the ion-ion interactions. The pseudo-potential depends only on the electron density and a reference core density.
2.2. Electron Density Profile
Based on a given ion-ion interaction from the pseudo-potential Hamiltonian, several electron density profiles may be prepared along the z-axis. This can be done before ion-core density simulation to achieve calculational efficiency. These electron density profiles can be prepared around the bulk density of the liquid metal. Based on a model “jellium distribution” which is a rigid electron-profile, the electron density profiles can be obtained by solving the Kohn-Sham equation. A typical representation of the electron density profiles is presented in Fig.1.
Figure 1. A normalized longitudinal electronic density profile of liquid Ga (oscillation) based on a “jellium distribution” (solid) [Phys. Rev. E 1997, 56, 7033].
2.3. Monte Carlo Simulation
Usually, a common model of system simulation consists of a slab that contains N ions and nN electrons, where n is the valency of the metal. The particles are randomly placed on the slab that is parallel to the x-y plane. The simulation box has the 3-D boundaries L0 x L0 x 2L0 in the x, y, and z directions. The side of the box (L0) is selected such that the average density of the ions in the slab matches the density of the liquid metal. A schematic representation of a simulation slab at the initial condition is presented in Fig. 2. It shows a liquid BiGa binary alloy with a dilute Bi (4%) in Ga (96%) at the liquid-vapor interface.
Figure 2. A simulation slab at the outmost layer of the BiGa liquid metal-vapor interface where the XY-plane is parallel to the interface and the coordinates are in atomic units (au) [J. Chem. Phys. 1998, 108, 5055].
The total depth of the simulation slab may be arranged from 10-17 layers which are usually slightly more than the characteristic layering structure of the interface. Each layer is about an atomic diameter in thickness. Each side of the slab at the interface can be about 5-7 layers. The initial configuration eliminates ion-core to ion-core overlaps by a force-biased Monte Carlo simulation with periodic boundary conditions.
The most important feature of the physical structure in the liquid metal-vapor interface may be understood from equilibrium correlations of the particles both in the transverse and the longitudinal direction of the interface.
3.1. The Transverse Density Profile
The pair correlation function of a system describes how the density changes concerning the separation between two particles in the system. For liquid-metal-vapor interfaces, the transverse pair correlation function that describes how the density varies within a density layer parallel to the interface can be calculated from paired particles in the interfacial region
where VT is the volume of all the particles in the layer. Assume that all atoms are in the form of a spherical shell, then Vs is an average volume of the spherical shell between r and (r + Δr), NT is the number of particles in a transverse density layer, and N(r, Δr) is the average number of pairs of particles between r and (r + Δr).
where NT is the total number of particles in the slice, N(r, Δr) is the average number of pairs of particles between r and (r + Δr), VT is the total volume of all the particles in the slice, and Vs is the average volume of the intersection of the spherical shell between r and (r + Δr). A typical air-correlation function of liquid metal (or alloy) is shown in Fig. 3.
Figure 3. Pair-correlation functions of bulk liquid Ga at three different temperatures [Phys. Rev. E 1997, 56, 7033].
3.2. The Longitudinal Density Profile
For a given system, the longitudinal density profile is a key characteristic of the system structure at the interface. It provides a sensitive test for theory in comparison to the experimental observations. Fig. 4 is a typical normalized longitudinal density profile along with its corresponding electron density profile.
Figure 4. The longitudinal density distribution of liquid Ga with the corresponding electron density (dotted line) at the liquid-vapor interface [Phys. Rev. E 1997, 56, 7033].
For most of the reported systems of the liquid metal-vapor interface, either simple liquid metals or liquid metal alloys, pseudo-potential representation provides a qualitative agreement between theoretical predictions and the experimental observations. Fig. 5 shows a representative comparison between the experimental observation and the theoretical prediction from simulation based on pseudo-potential Hamiltonian. It shows that the simulated longitudinal density profile agrees qualitatively well with the experimental observation.
It has been noted that the density distributions are usually not sensitive to temperature. There is hardly any obvious difference between the density profiles at different temperatures. This is both for pure metals and for alloys.
Figure 5. A representative comparison of the density profile between the experiment (solid) and the simulation, using the BiGa binary alloy as an example [J. Chem. Phys. 1998, 108, 5055].
The theoretical description of pseudo-potential Hamiltonian with self-consistent Monte Carlo simulation has been successful in comparison to experimental observations. Studies on the liquid-vapor interface have provided us an improved understanding of these systems. In turn, it leads to successful applications in daily life. Contrarily, our understanding of the structures of liquid metal-solid interfaces has been very limited. This is mainly because of the limited experimental and theoretical studies reported so far.
There have been only a few studies of liquid metal-solid interfaces, as reported in the literature [35-39]. One of the first studies was an experimental observation and a self-consistent Monte Carlo simulation concerning the liquid Hg−Sapphire interface [35-36]. In that system, liquid Hg had an interface on the Sapphire (001) surface. We need focused studies on the liquid metal-solid interface to have a better understanding.
We are interested to know, based on our knowledge of the liquid metal-vapor interface, how to improve our understanding of liquid metal-solid interfaces. Specifically, we would like to know the differences between liquid metal-solid interfaces and liquid metal-vapor interfaces and what are the key elements that lead to these differences.