Abstract:Cu2S in its high-temperature phases is a promising candidate for thermoelectric materials with the combination of a solid S lattice and liquidlike Cu atoms in the context of the so-called phonon-liquid electron-crystal (PLEC) mechanism. The atomic dynamics associated with the high Cu mobility is a critical component in our understanding of the low thermal conductivity in these materials. Among various possible PLEC compounds that have been studied experimentally, hexagonal β-Cu2S is a unique case that exhibits anisotropic diffusion channels for Cu atoms in the presence of crystalline S layers. To unveil the diffusion mechanism for such liquidlike Cu atoms, we present first-principles molecular dynamics simulations over 50 ps for hexagonal Cu2S at 450 K. Quantitative analysis of the atomic radial distributions, mean-square displacements, velocity autocorrelations, and atomic trajectories is reported and confirms the liquid-solid hybrid features. Our simulations reveal the preference of threefold triangular sites by Cu atoms in the S lattice. The triangular sites in the interlayer region do not bind Cu atoms strongly, allowing them to diffuse in the horizontal direction between S layers like a mobile liquid with a calculated diffusion coefficient ∼5×10−6cm2/s. The corresponding atomic trajectories have a wide spread and cannot be described by the Chudley-Elliott jump diffusion model. In contrast, Cu atoms are more strongly confined at the triangular sites within the S planes, and Cu diffusion takes place only when the atom hops out of the S layer and enters the interlayer region. This yields a 50% smaller diffusion coefficient in the vertical direction. The anisotropic diffusion channels for liquidlike Cu atoms in hexagonal β-Cu2S may provide an additional degree of freedom in designing promising systems for future energy applications.

Misfit-dislocation-mediated heteroepitaxial island diffusion

Dynamical Behaviour of Ag and Cu Double-Layer Islands on Fcc (111), Surfaces

Diffusion of small two-dimensional Cu Islands on Cu(111)

Diffusion rates of Cu adatoms on Cu(111) in the presence of an adisland nucleated at FCC or HCP sites

Misfit Dislocations and Adatom Domain Competitions in Cu/Ni (11¯1¯) Heteroepitaxial Growth

Misfit dislocations and adatom domain competitions in Cu/Ni (1 1¯ 1¯) heteroepitaxial growth

Misfit Dislocations And Adatom Domain Competitions In Cu/Ni (1(1)Over-Bar(1)Over-Bar) Heteroepitaxial Growth

Decay mechanism of double-layer islands on close-packed surfaces: Silver on Ag(111) and copper on Cu(111)

Diffusions of small clusters on Pt(111) and Cu(111) surfaces

Self-assembled island formation in heteroepitaxial growth

Strain in heteroepitaxial growth

Dislocation-free 3D islands in highly mismatched epitaxy: An equilibrium study with anharmonic interactions

Ab initio study of subsurface diffusion of Cu on the H-passivated Si(001) surface

Surface Diffusion of Hydrogen Atoms on Cu(111) Studied by Optical Diffraction from Hydrogen Density Patterns Formed on Removable Templates of Xenon Monolayers

Epitaxial growth of Cu on Cu(001): experiments and simulations

Cu/Ag EAM Potential Optimized for Heteroepitaxial Diffusion from ab initio Data

Atomic Superlattice Formation Mechanism Revealed by Scanning Tunneling Microscopy and Kinetic Monte Carlo Simulations

Ab initio study of step formation and self-diffusion on Ag(100)

Diffusion of Cu adatoms and dimers on Cu(111) and Ag(111) surfaces

Simulation of Dimer Diffusion on Metal Fcc (001) Surfaces by Molecular Dynamics

Mechanism for anisotropic diffusion of liquid-like Cu atoms in hexagonal β − Cu 2 S