Determination of the phonon sidebands in the photoluminescence spectrum of semiconductor nanoclusters from ab initio calculations

We propose a theoretical approach based on (constrained) density functional theory and the Franck-Condon approximation for the calculation of the temperature dependent photoluminescence of nanostructures. The method is computationally advantageous and only slightly more demanding than a standard density functional theory calculation and includes transitions into multiphonon final states (higher class transitions). We use the approach for Si and CdSe colloidal nanoclusters (NCs) with up to 693 atoms and obtain very good agreement with experiment which allows us to identify specific peaks and explain their origin. Generally, breathing type modes are shown to dominate the phonon replicas, while optical modes have significant contributions for CdSe NCs and play a lesser role in Si NCs. We obtain significant anti-Stokes peak starting at 140K for Si NC explaining the broadening observed in the corresponding experiment. We also apply the method to small molecular-like carbon structures (diamondoids), where electron-phonon coupling is typically large, and find that multiphonon processes (up to class 4) are very relevant and necessary to compare favorably with experiment. While it is crucial to include these multiphonon states in the small diamondoids with few tens of atoms, neglecting them in only marginally larger Si$_{87}$H$_{76}$ and Cd$_{43}$Se$_{44}$H$^*_{76}$ (and larger) quantum dots represents a good approximation.
Comment: Corrected version of Physical Review B 106, 245404 (2022). An Erratum of this article will be published shortly with PRB as well. This pdf is provided for convenience, please cite the original PRB reference including the erratum reference