Monte-Carlo simulation combined with density functional theory to investigate the equilibrium thermodynamics of electrode materials: lithium titanates as model compounds. Issue 28 (2nd July 2019)

• Record Type: Journal Article
• Authors: Ozaki, Hiroyuki; Tada, Kohei; Kiyobayashi, Tetsu
• Is Part Of: Physical chemistry chemical physics. Volume 21:Issue 28(2019)
• Page Start: 15551
• Page End: 15559
• Publication Date: 2019-07-02
• Subjects: Chemistry, Physical and theoretical -- Periodicals; 541.3
• DOI: 10.1039/c9cp02107h
• Languages: English
• ISSNs: 1463-9076

Abstract

Abstract : The present study proposes a method to describe the equilibrium thermodynamics of a model electrode material, LTO, based on the Monte-Carlo simulation (MC), for which the energetic parameters are determined by the density functional theory (DFT). Abstract : By using lithium titanate (LTO) as a model electrode material, the present study proposes a method to describe its equilibrium thermodynamics based on the Monte-Carlo simulation (MC), for which the energetic parameters are determined by the density functional theory (DFT). The electrochemical potential profile is simulated by a simple topological model which consists only of three parameters representing the Li site energies; namely, the potential energy of the 8a site ( ε 8a ), the difference in the site energy between the 8a and 16c sites (Δ ε ) and the repulsion between two Li atoms situated at the adjacent 8a and 16c sites ( J ). Parameter physics by the MC revealed that the term Δ ε plays a decisive role, with a collateral effect from J, for characterizing the shape of the potential profile whereas the term ε 8a determines its position along the electrochemical potential. For instance, if Δ ε exceeds the thermal energy at the temperature under consideration, i.e., if Δ ε > 3 kT, the first-order phase transition takes place during which two phases coexist, resulting in a plateau region in the potential profile. On the other hand, if Δ ε < 3 kT, the lithiation of LTO is viewed as a phenomenon above the critical point, above which the material is in a homogeneous uniphasic state. A multiple regression analysis of a set of the total energy calculated by DFT allows us to determine these energetic terms. The MC simulation with the determined parameters well reproduces the shape and position of the experimental potential profile of LTO. Since the determined value, Δ ε /eV ∼ 0.4, far exceeds the thermal energy at ambient temperature, the potential plateau of LTO is explained by the first-order phase transition as long as the equilibrium state is concerned.