Phonon Anharmonicity in Few-Layer Black Phosphorus

We report a temperature-dependent Raman spectroscopy study of few-layer black phosphorus (BP) with varied incident polarization and sample thickness. The Raman-active modes Ag₁, B_2g, and Ag₂ exhibit a frequency downshift, while their line width tends to increase with increasing temperature. To understand the details of these phenomena, we perform first-principles density functional theory calculations on freestanding monolayer BP. The effect of thermal expansion is included by constraining the temperature-dependent lattice constant. The study of the temperature-induced shift of the phonon frequencies is carried out using ab initio molecular dynamics simulations. The normal-mode frequencies are calculated by identifying the peak positions from the magnitude of the Fourier transform of the total velocity autocorrelation. Anharmonicity induces a frequency shift for each individual mode, and the three- and four-phonon process coefficients are extracted. These results are compared with those obtained from many-body perturbation theory, giving access to phonon lifetimes and lattice thermal conductivity. We establish that the frequency downshift is primarily due to phonon-phonon scattering while thermal expansion only contributes indirectly by renormalizing the phonon-phonon scattering. Overall, the theoretical results are in excellent agreement with experiment, thus showing that controlling phonon scattering in BP could result in better thermoelectric devices or transistors that dissipate heat more effectively when confined to the nanoscale.