Rotation mechanism of methylammonium molecules in organometal halide perovskite in cubic phase: An ab initio molecular dynamics study
Related Experiment Video
Updated: Mar 10, 2026

Low Pressure Vapor-assisted Solution Process for Tunable Band Gap Pinhole-free Methylammonium Lead Halide Perovskite Films
Published on: September 8, 2017
Kohei Shimamura1, Tomoya Hakamata2, Fuyuki Shimojo2
1Graduate School of System Informatics, Kobe University, Kobe 657-8501, Japan.
Abstract:
Rotation of methylammonium (CH3NH3 or MA) molecules is believed to govern the excellent transport properties of photocarriers in the MA lead iodide (MAPbI3) 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-MAPbI3 at room temperature. There are two types of rotational motions of MA molecules in a crystalline PbI3 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 CH3 group (HC) hardly interacts with I atoms, those in the NH3 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 MAPbI3 as well, considering that the atomistic structure and time constants regarding the rotation of MA molecules statistically agree with bulk experiments.
More Related Videos
08:54Vibrational Spectra of a N719-Chromophore/Titania Interface from Empirical-Potential Molecular-Dynamics Simulation, Solvated by a Room Temperature Ionic Liquid
Published on: January 25, 2020
08:30Monovalent Cation Doping of CH3NH3PbI3 for Efficient Perovskite Solar Cells
Published on: March 19, 2017
Related Concept Videos
Valence Bond Theory
Spin–Spin Coupling Constant: Overview
Qualitatively, any spin plus-half nucleus polarizes the spins of its electrons to the minus-half state. Consequently, the paired electron in the hydrogen–carbon bond must...
Stability of Substituted Cyclohexanes
The two chair conformations of cyclohexanes undergo rapid interconversion at room temperature. Both forms have identical energies and stabilities, each comprising equal amounts of the equilibrium mixture. Replacing a hydrogen atom with a functional group makes the two conformations energetically non-equivalent.
For example, in...
ortho–para-Directing Activators: –CH3, –OH, –⁠NH2, –OCH3
Electrophilic Addition to Alkynes: Halogenation
Halogenation is another class of electrophilic addition reactions where a halogen molecule gets added across a π bond. In alkynes, the presence of two π bonds allows for the addition of two equivalents of halogens (bromine or chlorine). The addition of the first halogen molecule forms a trans-dihaloalkene as the major product and the cis isomer as the minor product. Subsequent addition of the second equivalent yields the tetrahalide.
Crystal Field Theory - Tetrahedral and Square Planar Complexes
Crystal field theory (CFT) is applicable to molecules in geometries other than octahedral. In octahedral complexes, the lobes of the dx2−y2 and dz2 orbitals point directly at the ligands. For tetrahedral complexes, the d orbitals remain in place, but with only four ligands located between the axes. None of the orbitals points directly at the tetrahedral ligands. However, the dx2−y2 and dz2 orbitals (along the Cartesian axes) overlap with the ligands less than the dxy,...