Phonon Confinement and Transport in Ultrathin Films

Abstract

Thermal transport by phonons in films with thicknesses of less than 10 nm is investigated in a soft system (Lennard-Jones argon) and a stiff system (Tersoff silicon) using two-dimensional lattice dynamics calculations and the Boltzmann transport equation. This approach uses a unit cell that spans the film thickness, which removes approximations related to the finite cross-plane dimension required in typical three-dimensional-based approaches. Molecular dynamics simulations are performed to obtain finite-temperature structures for the lattice dynamics calculations and to predict thermal conductivity benchmarks. Thermal conductivity decreases with decreasing film thickness until the thickness reaches four unit cells (2.1 nm) for argon and three unit cells (1.6 nm) for silicon. With a further decrease in film thickness, thermal conductivity plateaus in argon while it increases in silicon. This unexpected behavior, which we identify as a signature of phonon confinement, is a result of an increased contribution from low-frequency phonons, whose density of states increases as the film thickness decreases. Phonon mode-level analysis suggests that confinement effects emerge below thicknesses of ten unit cells (5.3 nm) for argon and six unit cells (3.2 nm) for silicon. Thermal conductivity predictions based on the bulk phonon properties combined with a boundary scattering model do not capture the low thickness behavior. To match the two-dimensional lattice dynamics and molecular dynamics predictions for larger thicknesses, the three-dimensional lattice dynamics calculations require a finite specularity parameter that in some cases approaches unity. These findings point to the challenges associated with interpreting experimental thermal conductivity measurements of ultrathin silicon films, where surface roughness and a native oxide layer impact phonon transport.

References in corpus (11)

• Heat transport in silicon from first principles calculations
• Thermal Conductivity and Phonon Transport in Suspended Few-Layer Hexagonal Boron Nitride
• Anisotropic in-plane thermal conductivity observed in few-layer black phosphorus
• Quantum mechanical prediction of four-phonon scattering rates and reduced thermal conductivity of solids
• Force and heat current formulas for many-body potentials in molecular dynamics simulation with applications to thermal conductivity calculations
• Efficient molecular dynamics simulations with many-body potentials on graphics processing units
• Nanophononic metamaterial: Thermal conductivity reduction by local resonance
• Cross-plane heat conduction in thin films with ab-initio phonon dispersions and scattering rates
• Two orders of magnitude reduction in silicon membrane thermal conductivity by resonance hybridizations
• Native surface oxide turns alloyed silicon membranes into nanophononic metamaterials with ultra-low thermal conductivity
• Mapping Phonon Modes from Reduced-Dimensional to Bulk Systems

Cited by in corpus (6)

• GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations
• First-principles based analysis of thermal transport in metallic nanostructures: size effect and Wiedemann-Franz law
• Impact of dimensional crossover on phonon transport in van der Waals materials: a case study of graphite and graphene
• Surface phonons limit heat conduction in thin films
• Strain-Gradient-Driven Decoupling of Thermal Suppression from Anisotropy in \b{eta}-Ga2O3
• Phonon-confinement theory of thermal conductivity in ultrathin silicon films