Disorder to order: how halide mixing in MAPbI3-xBrx perovskites restricts MA dynamics.

Mixed-halide lead perovskites are of particular interest for the design of tandem solar cells currently reaching record efficiencies. While halide phase segregation upon illumination of mixed perovskites is extensively studied, the effect of halide disorder on A cation dynamics is not well understood, despite its importance for charge carrier diffusion and lifetime. Here, we study the methylammonium (MA) reorientational dynamics in mixed halide MAPbI3-xBrx perovskites by a combined approach of experimental solid-state NMR spectroscopy and molecular dynamics (MD) simulations based on machine-learning force-fields (MLFF). 207Pb NMR spectra indicate the halides are randomly distributed over their lattice positions, whereas PXRD measurements show that all mixed MAPbI3-xBrx samples are cubic. The experimental 14N spectra and 1H double-quantum (DQ) NMR data reveal anisotropic MA reorientations depending on the halide composition and thus associated disorder in the inorganic sublattice. MD calculations allow us to correlate these experimental results to restrictions of MA dynamics due to preferred MA orientations in their local Pb8I12-nBrn "cages". Based on the experimental and simulated results, we develop a phenomenological model that correlates the 1H dipolar coupling and thus the MA dynamics with the local composition and reproduces the experimental data over the whole composition range. We show that the dominant interaction between the MA cations and the Pb-X lattice that influences the cation dynamics is the local electrostatic potential being inhomogeneous in mixed halide systems. As such, we generate a fundamental understanding of the predominant interaction between the MA cations and the inorganic sublattice, as well as MA dynamics in asymmetric halide coordinations.

Exploring Librational Pathways with on-the-Fly Machine-Learning Force Fields: Methylammonium Molecules in MAPbX3 (X = I, Br, Cl) Perovskites.

Two seemingly similar crystal structures of the low-temperature (∼100 K) MAPbX3 (X = I, Br, Cl) perovskites, but with different relative methylammonium (MA) ordering, have appeared as representatives of this orthorhombic phase. Distinguishing them by X-ray diffraction experiments is difficult, and conventional first-principles-based molecular dynamics approaches are often too computationally intensive to be feasible. Therefore, to determine the thermodynamically stable structure, we use a recently introduced on-the-fly machine-learning force field method, which reduces the computation time from years to days. The molecules exhibit a large degree of anharmonic motion depending on temperature: that is, rattling, twisting, and tumbling. We observe the crystal's "librational pathways" while slowly heating it in isothermal-isobaric simulations. Marked differences in the thermal evolution of structural parameters allow us to determine the real structure of the system via a comparison with experimentally determined crystal structures.

Phase Transitions of Hybrid Perovskites Simulated by Machine-Learning Force Fields Trained on the Fly with Bayesian Inference.

Realistic finite temperature simulations of matter are a formidable challenge for first principles methods. Long simulation times and large length scales are required, demanding years of computing time. Here we present an on-the-fly machine learning scheme that generates force fields automatically during molecular dynamics simulations. This opens up the required time and length scales, while retaining the distinctive chemical precision of first principles methods and minimizing the need for human intervention. The method is widely applicable to multielement complex systems. We demonstrate its predictive power on the entropy driven phase transitions of hybrid perovskites, which have never been accurately described in simulations. Using machine learned potentials, isothermal-isobaric simulations give direct insight into the underlying microscopic mechanisms. Finally, we relate the phase transition temperatures of different perovskites to the radii of the involved species, and we determine the order of the transitions in Landau theory.

Assessing Density Functionals Using Many Body Theory for Hybrid Perovskites.

Which density functional is the "best" for structure simulations of a particular material? A concise, first principles, approach to answer this question is presented. The random phase approximation (RPA)-an accurate many body theory-is used to evaluate various density functionals. To demonstrate and verify the method, we apply it to the hybrid perovskite MAPbI_{3}, a promising new solar cell material. The evaluation is done by first creating finite temperature ensembles for small supercells using RPA molecular dynamics, and then evaluating the variance between the RPA and various approximate density functionals for these ensembles. We find that, contrary to recent suggestions, van der Waals functionals do not improve the description of the material, whereas hybrid functionals and the strongly constrained appropriately normed (SCAN) density functional yield very good agreement with the RPA. Finally, our study shows that in the room temperature tetragonal phase of MAPbI_{3}, the molecules are preferentially parallel to the shorter lattice vectors but reorientation on ps time scales is still possible.

Behavior of Methylammonium Dipoles in MAPbX3 (X = Br and I).

Dielectric constants of MAPbX3 (X = Br, I) in the 1 kHz-1 MHz range show strong temperature dependence near room temperature, in contrast to the nearly temperature-independent dielectric constant of CsPbBr3. This strong temperature dependence for MAPbX3 in the tetragonal phase is attributed to the MA+ dipoles rotating freely within the probing time scale. This interpretation is supported by ab initio molecular dynamics simulations on MAPbI3 that establish these dipoles as randomly oriented with a rotational relaxation time scale of ∼7 ps at 300 K. Further, we probe the intriguing possibility of transient polarization of these dipoles following a photoexcitation process with important consequences on the photovoltaic efficiency, using a photoexcitation pump and second harmonic generation efficiency as a probe with delay times spanning 100 fs-1.8 ns. The absence of a second harmonic signal at any delay time rules out the possibility of any transient ferroelectric state under photoexcitation.