1
|
Mattoni A, Argiolas S, Cozzolino G, Dell'Angelo D, Filippetti A, Caddeo C. Many-Body MYP2 Force-Field: Toward the Crystal Growth Modeling of Hybrid Perovskites. J Chem Theory Comput 2024. [PMID: 39066691 DOI: 10.1021/acs.jctc.4c00391] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
Hybrid perovskites are well-known for their optoelectronic and photovoltaic properties. Molecular dynamics simulations allow the study of these soft and ionic crystals by including dynamical effects (e.g., molecular rotations, octahedra tilting, ionic diffusion and hysteresis), yet the high computational cost restricts the use of accurate ab initio forces for bulk or small atomic systems. Hence, great interest exists in the development of classical force-fields for hybrid perovskites of low and linear scaling computational cost, via both empirical methods and machine-learning. This work aims at extending the transferability of our MYP0 model, which has been successfully tailored to methylammonium lead iodide (MAPI) and applied to the study of molecular rotations, vibrations, diffusion of defects, and many other properties. The extended model, named MYP2, improves the description of inorganic or hybrid fragments and the processes of crystal formation while preserving a good description of bulk properties. More importantly, it allows for the direct simulation of the crystal growth of cubic MAPI from deposition of PbI and MAI precursors on the surfaces. Our findings pave the way toward classical force-fields able to model the microstructure evolution of hybrid perovskites and the crystalline synthesis from deposition in vacuo.
Collapse
Affiliation(s)
- Alessandro Mattoni
- CNR - Istituto Officina dei Materiali (IOM), Cagliari, Cittadella Universitaria, Monserrato (CA) 09042, Italy
| | - Simone Argiolas
- CNR - Istituto Officina dei Materiali (IOM), Cagliari, Cittadella Universitaria, Monserrato (CA) 09042, Italy
- Dipartimento di Fisica, Università degli Studi di Cagliari, Cittadella Universitaria, Monserrato (CA) 09042, Italy
| | - Giacomo Cozzolino
- CNR - Istituto Officina dei Materiali (IOM), Cagliari, Cittadella Universitaria, Monserrato (CA) 09042, Italy
| | - David Dell'Angelo
- CNR - Istituto Officina dei Materiali (IOM), Cagliari, Cittadella Universitaria, Monserrato (CA) 09042, Italy
- Dipartimento di Fisica, Università degli Studi di Cagliari, Cittadella Universitaria, Monserrato (CA) 09042, Italy
| | - Alessio Filippetti
- CNR - Istituto Officina dei Materiali (IOM), Cagliari, Cittadella Universitaria, Monserrato (CA) 09042, Italy
- Dipartimento di Fisica, Università degli Studi di Cagliari, Cittadella Universitaria, Monserrato (CA) 09042, Italy
| | - Claudia Caddeo
- CNR - Istituto Officina dei Materiali (IOM), Cagliari, Cittadella Universitaria, Monserrato (CA) 09042, Italy
| |
Collapse
|
2
|
Baldwin WJ, Liang X, Klarbring J, Dubajic M, Dell'Angelo D, Sutton C, Caddeo C, Stranks SD, Mattoni A, Walsh A, Csányi G. Dynamic Local Structure in Caesium Lead Iodide: Spatial Correlation and Transient Domains. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2024; 20:e2303565. [PMID: 37736694 DOI: 10.1002/smll.202303565] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/28/2023] [Revised: 07/31/2023] [Indexed: 09/23/2023]
Abstract
Metal halide perovskites are multifunctional semiconductors with tunable structures and properties. They are highly dynamic crystals with complex octahedral tilting patterns and strongly anharmonic atomic behavior. In the higher temperature, higher symmetry phases of these materials, several complex structural features are observed. The local structure can differ greatly from the average structure and there is evidence that dynamic 2D structures of correlated octahedral motion form. An understanding of the underlying complex atomistic dynamics is, however, still lacking. In this work, the local structure of the inorganic perovskite CsPbI3 is investigated using a new machine learning force field based on the atomic cluster expansion framework. Through analysis of the temporal and spatial correlation observed during large-scale simulations, it is revealed that the low frequency motion of octahedral tilts implies a double-well effective potential landscape, even well into the cubic phase. Moreover, dynamic local regions of lower symmetry are present within both higher symmetry phases. These regions are planar and the length and timescales of the motion are reported. Finally, the spatial arrangement of these features and their interactions are investigated and visualized, providing a comprehensive picture of local structure in the higher symmetry phases.
Collapse
Affiliation(s)
- William J Baldwin
- Department of Engineering, University of Cambridge, Cambridge, CB2 1PZ, UK
| | - Xia Liang
- Department of Materials, Imperial College London, London, SW7 2AZ, UK
| | - Johan Klarbring
- Department of Materials, Imperial College London, London, SW7 2AZ, UK
- Department of Physics, Chemistry and Biology (IFM), Linköping University, Linköping, SE-581 83, Sweden
| | - Milos Dubajic
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Cambridge, CB3 0AS, UK
| | | | - Christopher Sutton
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, SC, 29208, USA
| | - Claudia Caddeo
- CNR-IOM, Unitá di Cagliari, Monserrato, Caligari, 09042, Italy
| | - Samuel D Stranks
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Cambridge, CB3 0AS, UK
| | | | - Aron Walsh
- Department of Materials, Imperial College London, London, SW7 2AZ, UK
| | - Gábor Csányi
- Department of Engineering, University of Cambridge, Cambridge, CB2 1PZ, UK
| |
Collapse
|
3
|
Ahlawat P. Crystallization of FAPbI3: Polytypes and stacking faults. J Chem Phys 2023; 159:151102. [PMID: 37846954 DOI: 10.1063/5.0165285] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Accepted: 09/15/2023] [Indexed: 10/18/2023] Open
Abstract
Molecular dynamics simulations are performed to study the crystallization of formamidinium lead iodide. From all-atom simulations of the crystal growth process and the δ-α-phase transitions, we try to reveal the formation of various stack-faulted intermediate defected structures and report various polytypes of formamidinium lead iodide that are observed from simulations.
Collapse
Affiliation(s)
- Paramvir Ahlawat
- SNSF Post-doc Mobility Fellow, Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom and Institute of Chemical Sciences and Engineering, Ecole Polytechnique Federale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
4
|
Wang L, Wang Y, Wang CS, Hao Q. An Application of Dipole–Dipole Interaction Model in Stacking Dimers Including Protonated Arginine Residue. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2022. [DOI: 10.1134/s0036024422090205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
5
|
Eatmon Y, Romiluyi O, Ganley C, Ni R, Pelczer I, Clancy P, Rand BP, Schwartz J. Untying the Cesium "Not": Cesium-Iodoplumbate Complexation in Perovskite Solution-Processing Inks Has Implications for Crystallization. J Phys Chem Lett 2022; 13:6130-6137. [PMID: 35759533 DOI: 10.1021/acs.jpclett.2c01344] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
We illustrate the critical importance of the energetics of cation-solvent versus cation-iodoplumbate interactions in determining the stability of ABX3 perovskite precursors in a dimethylformamide (DMF) solvent medium. We have shown, through a complementary suite of nuclear magnetic resonance (NMR) and computational studies, that Cs+ exhibits significantly different solvent vs iodoplumbate interactions compared to organic A+-site cations such as CH3NH3+ (MA+). Two NMR studies were conducted: 133Cs NMR analysis shows that Cs+ and MA+ compete for coordination with PbI3- in DMF. 207Pb NMR studies of PbI2 with cationic iodides show that perovskite-forming Cs+ (and, somewhat, Rb+) do not comport with the 207Pb chemical shift trend found for Li+, Na+, and K+. Three independent computational approaches (density functional theory (DFT), ab initio Molecular Dynamics (AIMD), and a polarizable force field within Molecular Dynamics) yielded strikingly similar results: Cs+ interacts more strongly with the PbI3- iodoplumbate than does MA+ in a polar solvent environment like DMF. The stronger energy preference for PbI3- coordination of Cs+ vs MA+ in DMF demonstrates that Cs+ is not simply a postcrystallization cation "fit" for the perovskite A+-site. Instead, it may facilitate preorganization of the framework precursor that eventually transforms into the crystalline perovskite structure.
Collapse
Affiliation(s)
- Yannick Eatmon
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Oluwaseun Romiluyi
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Connor Ganley
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Ruihao Ni
- Department of Materials Science and Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - István Pelczer
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Paulette Clancy
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Barry P Rand
- Department of Electrical Engineering and Andlinger Center for Energy and the Environment, Princeton University, Princeton, New Jersey 08544, United States
| | - Jeffrey Schwartz
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
6
|
Gale JD, LeBlanc LM, Spackman PR, Silvestri A, Raiteri P. A Universal Force Field for Materials, Periodic GFN-FF: Implementation and Examination. J Chem Theory Comput 2021; 17:7827-7849. [PMID: 34735764 DOI: 10.1021/acs.jctc.1c00832] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
In this study, the adaption of the recently published molecular GFN-FF for periodic boundary conditions (pGFN-FF) is described through the use of neighbor lists combined with appropriate charge sums to handle any dimensionality from 1D polymers to 2D surfaces and 3D solids. Numerical integration over the Brillouin zone for the calculation of π bond orders of periodic fragments is also included. Aside from adapting the GFN-FF method to handle periodicity, improvements to the method are proposed in regard to the calculation of topological charges through the inclusion of a screened Coulomb term that leads to more physical charges and avoids a number of pathological cases. Short-range damping of three-body dispersion is also included to avoid collapse of some structures. Analytic second derivatives are also formulated with respect to both Cartesian and strain variables, including prescreening of terms to accelerate the dispersion/coordination number contribution to the Hessian. The modified pGFN-FF scheme is then applied to a wide range of different materials in order to examine how well this universal model performs.
Collapse
Affiliation(s)
- Julian D Gale
- Curtin Institute for Computation, School of Molecular and Life Sciences, Curtin University, PO Box U1987, Perth, Western Australia 6845, Australia
| | - Luc M LeBlanc
- Curtin Institute for Computation, School of Molecular and Life Sciences, Curtin University, PO Box U1987, Perth, Western Australia 6845, Australia
| | - Peter R Spackman
- Curtin Institute for Computation, School of Molecular and Life Sciences, Curtin University, PO Box U1987, Perth, Western Australia 6845, Australia
| | - Alessandro Silvestri
- Curtin Institute for Computation, School of Molecular and Life Sciences, Curtin University, PO Box U1987, Perth, Western Australia 6845, Australia
| | - Paolo Raiteri
- Curtin Institute for Computation, School of Molecular and Life Sciences, Curtin University, PO Box U1987, Perth, Western Australia 6845, Australia
| |
Collapse
|
7
|
Jing Z, Liu C, Ren P. Advanced Electrostatic Model for Monovalent Ions Based on Ab Initio Energy Decomposition. J Chem Inf Model 2021; 61:2806-2817. [PMID: 34096706 PMCID: PMC8323402 DOI: 10.1021/acs.jcim.1c00426] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Ions play important roles in the structures and functions of biomolecules. In biomolecular simulations, ions either directly interact with biomolecules or provide an ionic environment that influences electrostatic interactions of solutes. The AMOEBA+ water model has demonstrated significant advancement of the classical force field for describing molecular interactions due to its improvements on the functional forms to account for essential physics. This work expands the applicability of the AMOEBA+ model toward alkali metal (Li, Na, K, Rb, and Cs) and halogen (F, Cl, Br, and I) ions. Various quantum chemical data on ion-ion and ion-water interactions, experimental ion hydration free energies, and lattice energies of salt crystals are used in the parametrization. The final parameters are verified with other properties outside of the parametrization data, including lattice energies of additional salt crystals and ionic activity coefficients in solution. The new model captures a wide range of ion properties from the gas phase to solution phase and crystals. More importantly, AMOEBA+ provides energy components that are consistent with ab initio energy decomposition. Thus, we expect AMOEBA+ to be more general, transferable, and valuable for the interpretation of intermolecular forces in efficient classical simulations.
Collapse
Affiliation(s)
- Zhifeng Jing
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712
| |
Collapse
|
8
|
Pols M, Vicent-Luna JM, Filot I, van Duin ACT, Tao S. Atomistic Insights Into the Degradation of Inorganic Halide Perovskite CsPbI 3: A Reactive Force Field Molecular Dynamics Study. J Phys Chem Lett 2021; 12:5519-5525. [PMID: 34096726 PMCID: PMC8215647 DOI: 10.1021/acs.jpclett.1c01192] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Accepted: 06/02/2021] [Indexed: 06/12/2023]
Abstract
Halide perovskites make efficient solar cells but suffer from several stability issues. The characterization of these degradation processes is challenging because of the limited spatiotemporal resolution in experiments and the absence of efficient computational methods to study these reactive processes. Here, we present the first reactive force field for molecular dynamics simulations of the phase instability and the defect-induced degradation in CsPbI3. We find that the phase transitions are driven by the anharmonic fluctuations of the atoms in the perovskite lattice. At low temperatures, the Cs cations tend to move away from their preferential positions, resulting in worse contacts with the surrounding metal halide framework which initiates the conversion to a nonperovskite phase. Moreover, our simulations of defective structures reveal that, although both iodine vacancies and interstitials are mobile in the perovskite lattice, the vacancies have a detrimental effect on the stability, leading to the decomposition of perovskites to PbI2.
Collapse
Affiliation(s)
- Mike Pols
- Materials
Simulation & Modelling, Department of Applied Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
- Laboratory
of Inorganic Materials Chemistry, Schuit Institute of Catalysis, Department
of Chemical Engineering and Chemistry, Eindhoven
University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
- Center
for Computational Energy Research, Department of Applied Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
| | - José Manuel Vicent-Luna
- Materials
Simulation & Modelling, Department of Applied Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
- Center
for Computational Energy Research, Department of Applied Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
| | - Ivo Filot
- Laboratory
of Inorganic Materials Chemistry, Schuit Institute of Catalysis, Department
of Chemical Engineering and Chemistry, Eindhoven
University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
- Center
for Computational Energy Research, Department of Applied Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
| | - Adri C. T. van Duin
- Department
of Mechanical Engineering, Pennsylvania
State University, University
Park, Pennsylvania 16802, United States
| | - Shuxia Tao
- Materials
Simulation & Modelling, Department of Applied Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
- Center
for Computational Energy Research, Department of Applied Physics, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands
| |
Collapse
|
9
|
Schriber JB, Nascimento DR, Koutsoukas A, Spronk SA, Cheney DL, Sherrill CD. CLIFF: A component-based, machine-learned, intermolecular force field. J Chem Phys 2021; 154:184110. [PMID: 34241025 DOI: 10.1063/5.0042989] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Computation of intermolecular interactions is a challenge in drug discovery because accurate ab initio techniques are too computationally expensive to be routinely applied to drug-protein models. Classical force fields are more computationally feasible, and force fields designed to match symmetry adapted perturbation theory (SAPT) interaction energies can remain accurate in this context. Unfortunately, the application of such force fields is complicated by the laborious parameterization required for computations on new molecules. Here, we introduce the component-based machine-learned intermolecular force field (CLIFF), which combines accurate, physics-based equations for intermolecular interaction energies with machine-learning models to enable automatic parameterization. The CLIFF uses functional forms corresponding to electrostatic, exchange-repulsion, induction/polarization, and London dispersion components in SAPT. Molecule-independent parameters are fit with respect to SAPT2+(3)δMP2/aug-cc-pVTZ, and molecule-dependent atomic parameters (atomic widths, atomic multipoles, and Hirshfeld ratios) are obtained from machine learning models developed for C, N, O, H, S, F, Cl, and Br. The CLIFF achieves mean absolute errors (MAEs) no worse than 0.70 kcal mol-1 in both total and component energies across a diverse dimer test set. For the side chain-side chain interaction database derived from protein fragments, the CLIFF produces total interaction energies with an MAE of 0.27 kcal mol-1 with respect to reference data, outperforming similar and even more expensive methods. In applications to a set of model drug-protein interactions, the CLIFF is able to accurately rank-order ligand binding strengths and achieves less than 10% error with respect to SAPT reference values for most complexes.
Collapse
Affiliation(s)
- Jeffrey B Schriber
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| | - Daniel R Nascimento
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| | - Alexios Koutsoukas
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Steven A Spronk
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| |
Collapse
|
10
|
Hong SN, Yu CJ, Jong UG, Choe SH, Kye YH. Effect of vacancy concentration on the lattice thermal conductivity of CH 3NH 3PbI 3: a molecular dynamics study. RSC Adv 2021; 11:34015-34023. [PMID: 35497285 PMCID: PMC9042382 DOI: 10.1039/d1ra05393k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Accepted: 08/19/2021] [Indexed: 12/04/2022] Open
Abstract
Hybrid halide perovskites are drawing great interest for photovoltaic and thermoelectric applications, but the relationship of thermal conductivities with vacancy defects remains unresolved. Here, we present a systematic investigation of the thermal conductivity of perfect and defective CH3NH3PbI3, performed using classical molecular dynamics with an ab initio-derived force field. We calculate the lattice thermal conductivity of perfect CH3NH3PbI3 as the temperature increases from 300 K to 420 K, confirming a good agreement with the previous theoretical and experimental data. Our calculations reveal that the thermal conductivities of defective systems at 330 K, containing vacancy defects such as VMA, VPb and VI, decrease overall with some slight rises, as the vacancy concentration increases from 0 to 1%. We show that such vacancies act as phonon scattering centers, thereby reducing the thermal conductivity. Moreover, we determine the elastic moduli and sound velocities of the defective systems, revealing that their slower sound speed is responsible for the lower thermal conductivity. These results could be useful for developing hybrid halide perovskite-based solar cells and thermoelectric devices with high performance. Molecular dynamics simulations with the MYP force field were performed to determine the thermal conductivity of perfect and defective bulk MAPbI3. Thermal conductivity was found to decrease overall as the vacancy concentration increased.![]()
Collapse
Affiliation(s)
- Song-Nam Hong
- Computational Materials Design, Faculty of Materials Science, Kim Il Sung University, Ryongnam-Dong, Taesong District, Pyongyang, Democratic People’s Republic of Korea
| | - Chol-Jun Yu
- Computational Materials Design, Faculty of Materials Science, Kim Il Sung University, Ryongnam-Dong, Taesong District, Pyongyang, Democratic People’s Republic of Korea
| | - Un-Gi Jong
- Computational Materials Design, Faculty of Materials Science, Kim Il Sung University, Ryongnam-Dong, Taesong District, Pyongyang, Democratic People’s Republic of Korea
| | - Song-Hyok Choe
- Computational Materials Design, Faculty of Materials Science, Kim Il Sung University, Ryongnam-Dong, Taesong District, Pyongyang, Democratic People’s Republic of Korea
| | - Yun-Hyok Kye
- Computational Materials Design, Faculty of Materials Science, Kim Il Sung University, Ryongnam-Dong, Taesong District, Pyongyang, Democratic People’s Republic of Korea
| |
Collapse
|
11
|
Affiliation(s)
- Xiaoyang Zhu
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|