OJPC  Vol.5 No.4 , November 2015
Molecular Dynamics Simulations of Perovskites: The Effect of Potential Function Representation on Equilibrium Structural Properties
Abstract: The perovskites with general formula ABX3 have been widely used as for materials with their unique properties (ferroelectric, piezoelectric, dielectric, catalytic and so on). Hybrid organolead halide perovskites are a class of semiconductors with ABX3 (X = Cl, Br, and I) structures consisting of lead cations in 6-fold coordination (B site), surrounded by an octahedron of halide anions (X site, face centered) together with the organic components in 12-fold cub octahedral coordination. These hybrid perovskites have a direct band gap, a large absorption coefficient as well as high charge carrier mobility that represent a very attractive characteristic of cost-effective solar cells. Basically, these crystals are inorganic solids of CaTiO3 type held together by bonds that are either ionic or partially ionic and partially covalent. In spite of the partially covalent character of the Ti-O bond, the system is modeled by a two-body central force interatomic potential (the form of the Vashishta and Rahman interatomic potential), which has been used successfully for many materials with a perovskite structure. In the present work using molecular dynamics (MD) simulation method we investigate the dynamical and structural behavior of CaTiO3 perovskite at normal pressure and temperature conditions. The MD calculations were performed on a system of 16,000 particles (3200Ca + 3200Ti + 96,00O), initially in an orthorhombic-Pbnm structure. The orthorhombic MD box had edges Lx = 53.4 Å, Ly = 53.4 Å and Lz = 61.12 Å, which provided a density matching the experimental value of ρ = 4 g/cm3. Starting with this structure and using proposed interatomic potentials the MD system stabilizes at room temperature in its initial configuration. The aim of the present study to explore the effect of potential function representations on structural equilibrium properties for the perovskite models including hybrid halide ones outlined above. Concerning the perovskite equilibrium state we elucidate the role of potential function modification on the atomic pair correlation and structural re-organization. The details of the interatomic potential representation have to be crucially important for obtaining of correct analysis data in crystallic, liquid and amorphous phases including perovskite systems.
Cite this paper: Kholmurodov, K. , Ibragimova, S. , Gladishev, P. , Vannikov, A. , Tameev, A. and Zelenyak, T. (2015) Molecular Dynamics Simulations of Perovskites: The Effect of Potential Function Representation on Equilibrium Structural Properties. Open Journal of Physical Chemistry, 5, 110-121. doi: 10.4236/ojpc.2015.54012.

[1]   Souza, J.A. and Rino, J.P. (2011) A Molecular Dynamics Study of Structural and Dynamical Correlations of CaTiO3. Acta Materialia, 59, 1409-1423.

[2]   Taibi-Benziada, L., Mezroua, A. and Von Der Mühll, R. (2004) CaTiO3 Related Materials for Resonators. Ceramics-Silikáty, 48, 180-184.

[3]   Calleja, M., Dove, M.T. and Salje, E.K.H. (2003) Trapping of Oxygen Vacancies on Twin Walls of CaTiO3: A Computer Simulation Study. Journal of Phyics: Condensed Matter, 15, 2301-2307.

[4]   Cai, B., Xing, Y.D., Yang, Z., Zhang, W.-H. and Qiu, J.S. (2013) High Performance Hybrid Solar Cells Sensitized by Organolead Halide Perovskitest. Energy & Environmental Science, 6, 1480.

[5]   Haruyama, J., Sodeyama, K., Han, L.Y. and Tateyama, Y. (2014) Termination Dependence of Tetragonal CH3NH3PbI3 Surfaces for Perovskite Solar Cells. Journal of Physical Chemistry Letters, 5, 2903-2909.

[6]   Jishi, R.A., Ta, O.B. and Sharif, A.A. (2014) Modeling of Lead Halide Perovskites for Photovoltaic Applications. Journal of Physical Chemistry C, 118, 28344-28349.

[7]   Mashiyama, H., Kawamura, Y., Kasano, H., Asahi, T., Noda, Y. and Kimura, H. (2007) Disordered Configuration of Methylammonium of CH3NH3PbBr3. Determined by Single Crystal Neutron Diffractometry, Ferroelectrics, 348, 182-186.

[8]   Vashishta, P. and Rahman, A. (1978) Ionic Motion in α-AgI. Physical Review Letters, 40, 1337-1340.

[9]   Voronov, V.N. (2006) Ion Mobility and Properties of Perovskite ABX3 Type Systems. Preprint № 000F. Physics Institute SB RAS, Krasnoyarsk, 64 p.

[10]   Razumov, V.F. and Klyuev, M.V., Eds. (2015) Organic and Hybrid Nanomaterials: Synthesis and Application Aspects: Monograph. Ivanovo State University, Ivanovo, 426-556.

[11]   Smith, W. and Forester, T.R. (1996) DL_POLY_2.0: A General-Purpose Parallel Molecular Dynamics Simulation Package. Journal of Molecular Graphics, 14, 136-141.

[12]   Smith, W., Forester, T.R. and Todorov, I.T. (2008) The DL POLY 2 User Manual. STFC Daresbury Laboratory Daresbury, Warrington, Version 2.19.

[13]   Kholmurodov, K., Ed. (2007) Molecular Simulation Studies in Material and Biological Sciences. International Workshop. Nova Science Publishers Ltd., New York, 196 p.

[14]   Kholmurodov, K., Ed. (2009) Molecular Simulation in Material and Biological Research. Nova Science Publishers Ltd., New York, 155 p.

[15]   Kholmurodov, K., Ed. (2011) Molecular Dynamics of Nanobistructures. Nova Science Publishers Ltd., New York, 210 p.

[16]   Kholmurodov, K., Ed. (2013) Models in Bioscience and Materials Research: Molecular Dynamics and Related Techniques. Nova Science Publishers Ltd., New York, 208 p.