We report a model describing the molecular orientation disorder in
CH3NH3PbI3, solving a classical Hamiltonian parametrised with electronic
structure calculations, with the nature of the motions informed by ab-initio
molecular dynamics. We investigate the temperature and static electric field
dependence of the equilibrium ferroelectric (molecular) domain structure and
resulting polarisability. A rich domain structure of twinned molecular dipoles
is observed, strongly varying as a function of temperature and applied electric
field. We propose that the internal electrical fields associated with
microscopic polarisation domains contribute to hysteretic anomalies in the
current--voltage response of hybrid organic-inorganic perovskite solar cells
due to variations in electron-hole recombination in the bulk.Comment: 10 pages; 4 figures, 2 SI figure