Cubic Perovskite Structure of Black Formamidinium Lead Iodide, α-[HC(NH2)2]PbI3, at 298 K

The structure of black formamidinium lead halide, α-[HC(NH2)2]PbI3, at 298 K has been refined from high resolution neutron powder diffraction data and found to adopt a cubic perovskite unit cell, a = 6.3620(8) Å. The trigonal planar [HC(NH2)2]+ cations lie in the central mirror plane of the unit cell with the formamidinium cations disordered over 12 possible sites arranged so that the C–H bond is directed into a cube face, whereas the −NH2 groups hydrogen bond (NH···I = 2.75–3.00 Å) with the iodide atoms of the [PbI3]− framework. High atomic displacement parameters for the formamidinium cation are consistent with rapid molecular rotations at room temperature as evidenced in ab initio molecular dynamic simulations.

H ybrid halide perovskite solar cells have captured significant research attention due to the demonstration of their use in highly efficient photovoltaic devices produced via low-cost deposition methods. 1−7 The complete structure and dynamics of the systems, covering the inorganic framework, the molecular cation and their mutual interactions, all impact upon the photovoltaic performance of these materials; understanding these fundamental structural and motional behaviors are central to developing this class of materials toward devices. 8,9 The bulk of research has focused on methylammonium lead iodide (MAPI); however, the formamidinium cation HC(NH 2 ) 2 + (FA) also templates a lead iodide perovskite-type structure for HC(NH 2 ) 2 PbI 3 (FAPI) with fully vertex connected PbI 6 octahedra. Solar cells based on FAPI with a maximum power conversion efficiency of 20.1% have been recently demonstrated. 10 The increase of effective cation radius on switching from the methylammonium to the formamidinium cation decreases the optical band gap from 1.54 to 1.47 eV, extending absorption into the near-infrared. 11 The synthesis, structure, and phase behavior of FAPbI 3 has been reported by Stoumpos et al. 12 FAPI crystallizes at room temperature as a nonperovskite, hexagonal yellow phase (δphase), in the space group P6 3 mc, and formed from face-sharing PbI 6 octahedra linked into chains. The structure of the black perovskite phase of FAPI (α-phase), which is formed at higher temperatures, was refined in this work from single crystal X-ray diffraction data using trigonal space groups and 3-fold disordered FA cation. FASnI 3 has been reported as adopting a cubic perovskite structure. 13 The low stability of hybrid perovskite materials toward heat and moisture has been discussed widely and results from the facile phase transitions at or near room temperature and formation of hydrated phases. These stability issues have hindered attempts to scale technology using these materials and serious question marks remain over stability even for sealed devices. 14 The formamidinium cation is believed to exhibit even lower stability than methylammonium. 12 MAPbI 3 and MA x FA 1−x PbI 3 have recently been shown to be efficient halide ion conductors 15 and undergo structural changes under polarization 16 or illumination. 17 Solid solutions, for example, FAPI/MAPbBr 3 have been employed in high efficiency perovskite cells, reportedly displaying greater stability than the parent compounds. 15,18,19 Binek et al. 20 have studied the stabilization of thin films of α-FAPI by inclusion of up to 15% MA + .
We have previously demonstrated the power of neutron powder diffraction in determining the atomic structure in hybrid perovskite compounds as it circumvents many of the factors limiting the usefulness of single crystal X-ray diffraction for analysis of these hybrid systems. As well the identification of light atom positions, including hydrogen, and distinguishing carbon from nitrogen, it also avoids difficulties in structure refinement stemming from multiple twinning or domain structure of crystals. 9 Neutron powder diffraction data were obtained on the HRPD instrument at ISIS as summarized in the Supporting Information. Initial inspection of the data showed far fewer reflections than predicted by the hexagonal unit cell model for FAPI proposed by Stoumpos (see Supporting Information Figures S1 and S2) indicating that the true unit cell was probably simpler than used in that work. Peak positions were extracted for 20 reflections and indexed manually from their dspacings to give a cubic unit cell of dimension ∼6.36 Å. This cell indexed all the reflections observed except for a few minor features associated with diffraction from the vanadium slab can (at, for example, 2.14 Å) and off-center aluminum cryoloop tail (these were readily identified from their d-spacings and varying intensities and positions in the different HRPD data banks). This cubic unit cell is similar to that found for the cubic phase adopted by MAPI above 327 K though it is slightly larger, which is indicative of the larger size of the formamidinium cation compared to methylammonium. In order to check the validity of this new choice of unit cell, powder X-ray diffraction (PXD) data were also collected from a FAPI phase, Supporting Information Figure S3, and could also be fully indexed of a simple cubic unit cell of dimensions 6.36 Å. Note that the sensitivity of neutron diffraction data to light elements, such as C, N, and H, shows that these atom positions can also be defined using the simple cubic description−and any deviations from cubic to hexagonal symmetry would direct diffraction intensity into additional reflections. Le Bail fits to the PXD data were undertaken using both the new cubic unit cell and the original hexagonal unit cell and showed an excellent fit (χ 2 = 2.18) using the simple cubic cell, Supporting Information, Figure S4.
A structure model for cubic FAPI was developed using that of cubic MAPI with the lead iodide framework described using nontilted PbI 6 octahedra in the space group Pm-3m. The position of a planar formamidinium cation within this cubic framework was energy minimized within density functional theory (DFT, PBEsol functional) within the regular 6.36 Å cube and found to adopt an orientation with the N−CH−N molecular ion fragment lying in the unit cell midplane, the (200) plane, with the C−H bond pointing directly toward the cube face (thereby minimizing CH···I interactions), and the C− N bonds directed toward adjacent cube faces. This allows the −NH 2 groups to be orientated toward the unit cell edges with the potential formation of NH···I hydrogen bonds. Overall, this produces 12 possible orientations for the triangular-planar formamidinium cation in the cubic cage: six orientations of the C−H bond, identified in this work as the y direction, and two orientations for the N−C−N plane, in xy and yz. Reorientation of the organic cation around these 12 orientations would be expected and has been studied by molecular dynamics, vide infra. The crystallographic description of this model in the space group Pm-3m is very simple with all the atoms on a single mirror plane; the carbon, C, and the directly bonded hydrogen, H1, both located on the 6-fold (1/2, y, 1/2), 6f site, with y ∼ 0.55 and 0.64 for C and H1, respectively, both with site occupancy 0.1667 representing the six possible formamidinium orientations. A single nitrogen position, (x, y, 1/2), 24l, with occupancy 0.0833, and two amine hydrogen positions also in the molecular plane on 24l, H2 near (0.81, 0.54, 0.5) and H3 near (0.70, 0.29, 1/2), site occupancies 0.0833. Thus, only 5 atom positions and 8 refinable coordinates are required to define the 12 different orientations of the complete planar [H 2 N−CH−NH 2 ] + cation. Pb was placed on 1a (0,0,0) and I on 3d (1/2, 0,0) and the refinable coordinates for the C, N, and H atoms were initially set at values derived from the energy minimization. Isotropic atomic displacement parameters (ADPs) for all atoms were set initially at 0.05 Å 2 . This model immediately produced a good overall fit to the data and refinement of profile parameters, unit cell parameter, background functions etcetera, led to a χ 2 value of only ∼4. Refinement of lead and iodine ADPs, anisotropic for iodine, reduced χ 2 to 3.5 and showed thermal motion for the iodine atom to be strongly perpendicular to the Pb−I−Pb direction as found previously in MAPI due to octahedral tilting. ADPs were constrained to a single value for all formamidinium atom sites and this value was refined to yield a value of approximately 0.14 Å 2 with χ 2 = 3.05 indicating significant thermal motion and reorientation of this structural element. Final stages of the refinement included an investigation of the formamidinium cation position which was permitted to translate as a rigid body in the y direction; a translation of the whole molecule, with the C−H bond moving toward the cube face by approximately 0.25 Å, yielded the final structural model given in Table 1 and χ 2 = 2.62. Final profile fits to the neutron powder diffraction data are given in Figure 1.  Figure 2. The FAPI lattice parameter, 6.3620(8) Å, and cell volume, 257.51(5) Å 3 , are significantly larger than those of MAPI (at 298 K), the reduced lattice parameter is 6.304 Å and cell volume 250.5 Å 3 , reflecting the larger organic cation size. The Pb−I bond lengths are increased to 3.181(4) Å, implying the lead iodide framework is under some expansive stress with weak Pb−I bonds. This stress may be the origin of the facile phase transition to the yellow phase of FAPI where the Pb−I bonds are significantly shorter. Also of interest are the extracted ADPs for the iodine atom site which show the expected strong anisotropic motion perpendicular to the Pb−I−Pb direction, similar values to those found for iodine in MAPI and representing a soft lead iodide framework.
The refined position of the FA + cation within the cubic cage is close to that predicted from the equilibrium athermal structure obtained from DFT-PBEsol calculations. Such calculations tend to underestimate the importance of weak hydrogen bonds such as the NH···I interactions found in FAPI and the cation adopts a site, which allows both terminal amine hydrogen atoms to form weak NH−I hydrogen bonds with the inorganic framework N(H2)···I at 2.992 Å and N(H3)···I at 2.770 Å. This position also maintains the C−H1 to iodide distances at, nonbonded, 3.572 Å. Information on the temporal behavior of the cation positions was obtained from molecular dynamics simulations at 300 K (on the DFT-PBEsol energy landscape; 200 ps production runs using a time-step of 0.5 fs). The molecule is rotationally active around room temperature with a preference for alignment of the C−H bond along ⟨100⟩, in agreement with the structure refinement from powder neutron diffraction, with a time constant of ∼2 ps associated with the rollover process. A video of the rotational dynamics   Figures S6 and S7).
Regarding the low temperature phase behavior, preliminary results indicate that cubic FAPI does undergo a transformation on cooling below room temperature−and by 200 K. Neutron powder data collected at 100 K show a more complex diffraction pattern, which is currently under analysis to determine the unit cell size and any ordering of the FA cation that takes place at low temperatures.
Previous work in the literature describes α-FAPI as being hexagonal, 12 but no evidence of the need to invoke this larger, more complex unit cell was found in this study, which is much more sensitive to the light atom structural elements. The [PbI 3 ] − framework is almost identical in the two studies, with Pb−I bonds near 3.18 Å; PXD patterns generated from the two models are also identical within experimental error. Of note is that faulted and twinned phases at the nanoscale can present themselves as hexagonal to single crystal diffraction methods, such as X-ray and electron diffraction. Recent work on the socalled lonsdaleite phase of diamond found it to consist of solely of cubic phase material that has extensive ⟨113⟩ twins and ⟨111⟩ stacking faults. 21 Given that α-FAPI is normally obtained by a thermally driven phase transition of hexagonal crystals of δ-FAPI it is likely that black FAPI forms as microdomains of twinned or faulted cubic phase material. Indeed our SXD studies of FAPI presented similar unit cells to those reported by Stoumpos but also complex diffraction patterns which could be resolved into multiple cubic twins. It has also been proposed that the conversion of δ-FAPI to α-FAPI might be associated with partial decomposition of the FA cations. 12 Our TGA results from δ-FAPI isolated in this study, Supporting Information Figure S5, show almost no weight loss associated with this process. Indeed, samples of cubic α-FAPI synthesized by the route described in this work (see experimental section) seem relatively stable in comparison with materials described in previous reports. We observe partial, minor conversion back to yellow δ-FAPI after 20 days (see Supporting Information Figure S3).
In conclusion, formamidinium lead iodide adopts a cubic unit cell at room temperature with the FA cations formally disordered over 24 orientations with the perovskitic [PbI 3 ] − framework. The molecules are highly rotationally mobile with a turnover frequency of 0.5 THz around room temperature. The structural similarity of the methylammonium and formamidinium lead perovskites has important implications on the understanding, characterization, and stability of their cubicphase solid-solution, as employed in the current champion efficiency hybrid perovskite solar cells. 10,20