Evaluation of the AMOEBA force field for simulating metal halide perovskites in the solid state and in solution

In this work, we compare the existing nonpolarizable force fields developed to study the solid or solution phases of hybrid organic-inorganic halide perovskites with the AMOEBA polarizable force field. The aim is to test whether more computationally expensive polarizable force fields like AMOEBA offer better transferability between solution and solid phases, with the ultimate goal being the study of crystal nucleation, growth, and other interfacial phenomena involving these ionic compounds. In the context of hybrid perovskites, AMOEBA force field parameters already exist for several elements in solution, and we decided to leave them unchanged and to only parameterize the missing ones (Pb²⁺ and CH₃NH₃ ⁺ ions) in order to maximize transferability and avoid overfitting to the specific examples studied here. Overall, we find that AMOEBA yields accurate hydration free energies (within 5%) for typical ionic species while showing the correct ordering of stability for the different crystal polymorphs of CsPbI₃ and CH₃NH₃PbI₃. Although the existing parameters do not accurately reproduce all transition temperatures and lattice parameters, AMOEBA offers better transferability between solution and solid states than existing nonpolarizable force fields.