Rotation mechanism of methylammonium molecules in organometal halide perovskite in cubic phase: An ab initio molecular dynamics study

Kohei Shimamura; Tomoya Hakamata; Fuyuki Shimojo; Rajiv K Kalia; Aiichiro Nakano; Priya Vashishta

doi:10.1063/1.4971791

Rotation mechanism of methylammonium molecules in organometal halide perovskite in cubic phase: An ab initio molecular dynamics study

J Chem Phys. 2016 Dec 14;145(22):224503. doi: 10.1063/1.4971791.

Authors

Kohei Shimamura¹, Tomoya Hakamata², Fuyuki Shimojo², Rajiv K Kalia³, Aiichiro Nakano³, Priya Vashishta³

Affiliations

¹ Graduate School of System Informatics, Kobe University, Kobe 657-8501, Japan.
² Department of Physics, Kumamoto University, Kumamoto 860-8555, Japan.
³ Collaboratory for Advanced Computing and Simulations, Department of Physics and Astronomy, Department of Computer Science, and Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089-0242, USA.

PMID: 27984900
DOI: 10.1063/1.4971791

Abstract

Rotation of methylammonium (CH₃NH₃ or MA) molecules is believed to govern the excellent transport properties of photocarriers in the MA lead iodide (MAPbI₃) perovskite. Of particular interest is its cubic phase, which exists in industrially important films at room temperature. In order to investigate the rotational behaviors of the MA molecules, we have performed ab initio molecular dynamics simulations of cubic-MAPbI₃ at room temperature. There are two types of rotational motions of MA molecules in a crystalline PbI₃ cage: reorientation of a whole molecule and intramolecular rotation around the C-N bond within MA molecules. Using a cubic symmetry-assisted analysis (CSAA), we found that the prominent orientation of the C-N bond is the crystalline ⟨110⟩ directions, rather than the ⟨100⟩ and ⟨111⟩ directions. Rapid rotation around the C-N bond is also observed, which easily occurs when the rotational axis is parallel to the ⟨110⟩ directions according to the CSAA. To explain the atomistic mechanisms underlying these CSAA results, we have focused on the relation between H-I hydrogen bonds and the orientation of an MA molecule. Here, the hydrogen bonds were defined by population analysis, and it has been found that, while H atoms in the CH₃ group (HC) hardly interacts with I atoms, those in the NH₃ group (HN) form at least one hydrogen bond with I atoms and their interatomic distances are in a wide range, 2.2-3.7 Å. Based on these findings, we have given a possible explanation to why the ⟨110⟩ directions are preferred. Namely, the atomic arrangement and interatomic distance between MA and surrounding I atoms are most suitable for the formation of hydrogen bonds. In addition to films, these results are potentially applicable to the rotational behaviors in bulk MAPbI₃ as well, considering that the atomistic structure and time constants regarding the rotation of MA molecules statistically agree with bulk experiments.