A computational method for angle-resolved photoemission spectra from repeated-slab band structure calculations

Misa Nozaki,Peter Krüger
2024-02-27
Abstract:A versatile method for angle-resolved photoemission spectra (ARPES) calculations is reported within the one-step model of photoemission. The initial states are obtained from a repeated-slab calculation using the projector-augmented wave (PAW) method. ARPES final states are constructed by matching the repeated-slab eigenstates of positive energy with free electron states that satisfy the time-reversed low-energy electron diffraction boundary conditions. Nonphysical solutions of the matching equations, which do not respect the flux conservation, are discarded. The method is applied to surface-normal photoemission from graphene as a function of photon energy from threshold up to 100 eV. The results are compared with independently performed multiple scattering calculations and very good agreement is obtained, provided that the photoemission matrix elements are computed with all-electron waves reconstructed from the PAW pseudo-waves. However, if the pseudo-waves are used directly, the relative intensity between $\sigma$- and $\pi$-band emission is wrong by an order of magnitude. The graphene ARPES intensity has a strong photon energy dependence including resonances. The normal emission spectrum from the $\pi$-band shows a hitherto unreported, sharp resonance at a photon energy of 31 eV. The resonance is due to a 2$D$ interband transitions and highlights the importance of matrix element effects beyond the final state plane-wave approximation.
Materials Science
What problem does this paper attempt to address?
This paper proposes a new method for calculating angle-resolved photoemission spectroscopy (ARPES) that is specifically designed for electronic structure calculations of layered structures. In the photoemission model, the initial state is obtained from calculations of layered structures using the projector-augmented wave (PAW) method. The final state of ARPES is constructed by matching the eigenstates of the periodic layered structures with the free electron states satisfying the time-reversal low-energy electron diffraction boundary conditions. Non-physical solutions that do not satisfy the conservation of flow matching equations are discarded. This method is applied to the photon energy dependence of the photoemission spectra in the normal direction of graphene surface, ranging from the threshold to 100 eV. The results are compared with independent multiple scattering calculations and show good agreement, but only when the photoemission matrix elements are calculated using the pseudowave reconstructed from the full-electron wave in PAW. If the pseudowave is used directly, the relative intensities of σ- and π-band emissions are incorrect by an order of magnitude. It is found that the ARPES intensity of graphene exhibits strong photon energy dependence, including resonance phenomena. The normal emission spectrum of the π-band shows a sharp resonance at 31 eV photon energy, which is attributed to the importance of interband transitions and matrix element effects beyond the final state plane wave approximation. In conclusion, this paper aims to address the accurate and flexible calculation of ARPES for two-dimensional materials, emphasizing the importance of full-electron wave functions, and demonstrating the effectiveness of their method and new discoveries through examples.