The development of high efficiency perovskite solar cells has sparked a multitude of measurements on the optical properties of these materials. For the most studied methylammonium(MA)PbI3 perovskite, a large range (6–55 meV) of exciton binding energies has been reported by various experiments. The existence of excitons at room temperature is unclear. For the MAPbX3 perovskites we report on relativistic Bethe-Salpeter Equation calculations (GW-BSE). This method is capable to directly calculate excitonic properties from first-principles. At low temperatures it predicts exciton binding energies in agreement with the reported ‘large’ values. For MAPbI3, phonon modes present in this frequency range have a negligible contribution to the ionic screening. By calculating the polarization in time from finite temperature molecular dynamics, we show that at room temperature this does not change. We therefore exclude ionic screening as an explanation for the experimentally observed reduction of the exciton binding energy at room temperature and argue in favor of the formation of polarons.

No abstract

No abstract

Linear optical properties can be calculated by solving the time-dependent density functional theory equations. Linearization of the equation of motion around the ground state orbitals results in the so-called Casida equation, which is formally very similar to the Bethe-Salpeter equation. Alternatively one can determine the spectral functions by applying an infinitely short electric field in time and then following the evolution of the electron orbitals and the evolution of the dipole moments. The long wavelength response function is then given by the Fourier transformation of the evolution of the dipole moments in time. In this work, we compare the results and performance of these two approaches for the projector augmented wave method. To allow for large time steps and still rely on a simple difference scheme to solve the differential equation, we correct for the errors in the frequency domain, using a simple analytic equation. In general, we find that both approaches yield virtually indistinguishable results. For standard density functionals, the time evolution approach is, with respect to the computational performance, clearly superior compared to the solution of the Casida equation. However, for functionals including nonlocal exchange, the direct solution of the Casida equation is usually much more efficient, even though it scales less beneficial with the system size. We relate this to the large computational prefactors in evaluating the nonlocal exchange, which renders the time evolution algorithm fairly inefficient.