ZHANG Yilong, CUI Man'ai, LIU Yanhui. Crystal Structure and Stability of LiAlH4 from First Principles[J]. Chinese Journal of High Pressure Physics, 2018, 32(2): 021103. doi: 10.11858/gywlxb.20170561
Citation: ZHANG Yilong, CUI Man'ai, LIU Yanhui. Crystal Structure and Stability of LiAlH4 from First Principles[J]. Chinese Journal of High Pressure Physics, 2018, 32(2): 021103. doi: 10.11858/gywlxb.20170561

Crystal Structure and Stability of LiAlH4 from First Principles

doi: 10.11858/gywlxb.20170561
Funds:

National Natural Science Foundation of China 11764043

More Information
  • Author Bio:

    ZHANG Yilong(1998—), male, undergraduate, major in computational physics under high pressure.E-mail:122696148@qq.com

  • Corresponding author: CUI Man'ai(1979—), female, master, lecturer, major in computational physics under high pressure.E-mail:macui@ybu.edu.cn; LIU Yanhui(1971—), female, Ph.D, professor, major in computational physics under high pressure.E-mail:yhliu@ybu.edu.cn
  • Received Date: 05 Apr 2017
  • Rev Recd Date: 30 Apr 2017
  • The structural stability of LiAlH4, a promising hydrogen storage material, under high pressure was researched using the ab initio pseudopotential plane wave method.It is found that the phase transition occurs at 1.6 GPa from the α-LiAlH4 phase to the β-LiAlH4 (space group I2/b) phase.This phase transition is identified as first-order in nature with volume contractions of 18%.Moreover, the analysis of the phonon dispersion curves suggests that phase transition is related to the phonon softening.Mulliken population analyses indicated that the ambient phase (α-LiAlH4) is expected to be the most promising candidate for hydrogen storage.

     

  • Light metal complex hydrides are promising materials for solid-state hydrogen due to their gravimetric hydrogen density[1-5].For example, lithium alanate (LiAlH4) and sodium alanate (NaAlH4) have theoretical hydrogen capacities of 10.6% and 7.5% (mass fraction), respectively.Hence, LiAlH4 and NaAlH4 could be viable candidates for practical usage as on-board hydrogen storage.However, a serious problem with these materials is their poor kinetics and lacking reversibility with respect to hydrogen absorption/desorption.Recent experimental evidences show that LiAlH4 and NaAlH4 release 7.9% and 5.6% (mass fraction) of H, respectively.This represents nearly 4 and 5 times more stored hydrogen than LaNi5-based alloys which are currently used in nickel-based hydride batteries.However, it is difficult to accurately identify the positions of the hydrogen atoms by the high-pressure X-ray and neutron diffraction studies.Therefore, exploration of the phase stability and structures has attracted a great deal of attention.

    Under ambient conditions, LiAlH4 crystallizes at a monoclinic structure with a space group P21/c (denoted as α-LiAlH4).The structure consists of AlH4 units separated by Li+ ions, and the hydrogen atoms are arranged around the aluminum atoms in an almost regular tetrahedral configuration.On the theoretical side, the α-LiAlH4 structure transforms into a new tetragonal β-LiAlH4 structure (space group I41/a) at 2.6 GPa with a 17% volume contraction[6].However, the neutron diffraction confirms that the structure of the high pressure β-LiAlD4 has a monoclinic space group I2/b[7].It is determined that there is a reversible phase transition back to the α-LiAlD4 with a slow release of pressure and cooling.Consequently, the Raman measurement confirmed that the ambient structure transform to the β-LiAlH4 structure (space group I2/b) at 3 GPa[8], as well as the previously reported transition pressure between 2.2 and 3.5 GPa[9].

    In this work, we report our results on the geometries and electronic structures of LiAlH4 with pressure using the ab initio calculations.The total-energy calculations for the pressure-induced phase transition in LiAlH4 were performed.The structural transition from α-to β-LiAlH4 happened at 1.6 GPa.The full phonon dispersion curves were calculated, for the first time, to provide the evidence of phonon stability for pressured-induced phase transition.We also discussed the density of states and bond overlap populations about the two phases.

    The calculations were performed within the density functional theory (DFT), using the Perdew-Burke-Ernzerh of generalized gradient approximation functional[10-12].A plane-wave norm-conserving pseudopotential method[13] as implemented in the CASTEP code[14] was employed.We used the plane-wave kinetic-energy cutoff of 850 eV which was shown to give excellent convergence of the total energies and structural parameters.According to the Monkhorst-Pack method, k-point spacing smaller than 0.3 nm-1 was individually adjusted in reciprocal space to the size of each computational cell[15].The optimization was performed, all forces on atoms were converged to be less than 0.1 eV/nm and all the stress components are less than 0.02 GPa.The tolerance in the self-consistent field (SCF) calculation was set to 10-6 eV/atom.The phonon frequencies were calculated by the direct approach, which is based on the first-principles calculations of the total energy, the Hellman-Feynman forces, and the dynamical matrix as implemented in the phonon packages[16-17].Convergence test gave the use of 2×1×1 and 2×2×1 supercell in the phonon calculation.

    We optimize the crystal structures allowing simultaneously variations of unit cell and atomic positions at selected pressures. The calculated lattice parameters of α-LiAlH4 (space group P21/c) and β-LiAlH4 (space group I2/b) are listed in Table 1, which also includes the experimental data for comparison[18].For the α-LiAlH4 phase, the deviations between our calculated results and the experimental values are less than 2%, which is sufficiently accurate.This strongly supports the choice of the pseudopotentials and GGA approximation for the current study.The theoretical generated pressure-volume curve (see Fig. 1) shows that the α-LiAlH4 transforms to a monoclinic phase of β-LiAlH4 with a volume contraction (ΔV/V) of 18%, which suggests that this is a first-order phase transition.The huge volume contraction is consistent with the Raman scatting measurement and the other theoretical results[6, 8].In order to give a clear picture of the structural transition, the difference of the enthalpy (per formula unit) between the α-LiAlH4 and β-LiAlH4 phases as a function of pressure is shown in the insert of Fig. 1.It can be clearly seen that the β-LiAlH4 phase becomes energetically more favourable than the α-LiAlH4 phase above 1.6 GPa, which is the transition pressure from α-LiAlH4 to β-LiAlH4.

    Table  1.  Optimized structural parameters, atomic position parameters for the α-LiAlH4 and the β-LiAlH4 structures
    Phase Unit-cell dimensions Atom coordinates
    α-LiAlH4
    (P21/c)
    a=0.485 1 nm(0.481 7 nm*)
    b=0.781 4 nm(0.780 2 nm*)
    c=0.773 2 nm(0.782 1 nm*)
    Li:(0.585, 0.459, 0.829), (0.560, 0.466, 0.827)*
    Al:(0.159, 0.204, 0.938), (0.139, 0.203, 0.930)*
    H1:(0.193, 0.102, 0.766), (0.183, 0.096, 0.763)*
    H2:(0.377, 0.371, 0.986), (0.352, 0.371, 0.975)*
    H3:(0.254, 0.082, 0.119), (0.243, 0.081, 0.115)*
    H4:(0.822, 0.268, 0.882), (0.799, 0.247, 0.872)*
    β-LiAlH4
    (I2/b)
    a=0.445 2 nm
    b=0.445 9 nm
    c=1.010 2 nm
    β=89.978°
    Li:(0, 0.250, 0.125)
    Al:(0, 0.250, 0.625)
    H1:(0.259, 0.425, 0.542)
    H2:(0.324, 0.508, 0.792)
    Note:"*" represents experimental data from Ref.[18].
     | Show Table
    DownLoad: CSV
    Figure  1.  Volume vs. pressure curves of α-LiAlH4 and β-LiAlH4 phases (Enthalpy difference (per formula unit) between α-LiAlH4 and β-LiAlH4 as a function of pressure is shown in the insert.)

    To investigate the dynamic stability, we calculated the phonon dispersion curves along some high-symmetry lines in the Brillouin zones (BZ) and the corresponding phonon density of states (DOS) for the LiAlH4 structure.No imaginary frequency is observed throughout the whole BZ, indicating that the two novel phases are dynamically stable in the pressure region from this study.We also calculated the phonon dispersion curves at the phase transition pressure for the α-LiAlH4 structure in Fig. 2.For comparison, the stability of phonon dispersion curves for the α-LiAlH4 at 0 GPa are shown in Fig. 2(a).In Fig. 2(b), it is indicated that the phase transition should not be related to the pressure-induced phonon softening.The total and partial densities of states for the two phases are shown in Fig. 3. Analysis of the calculated electronic density of state reveals that the two phases exhibit a common insulating feature with a finite energy gap.The valence band (VB) in α-LiAlH4 phase is split into two regions.Subsequently, the two peaks combine into one broader peak in the β-LiAlH4 phase producing the broadening of VB under high pressure.This is originated from the shortened interatomic distance upon squeezing, for the VB region is mainly dominated by H 1s, and Al 2s, 2p states.It shows that the H and Al atoms become the directional covalent bonds within the AlH4 tetrahedron or AlH4 octahedron layers. The bottom of the conduction band just above the Fermi energy is composed of the Al 2s, 2p, Li 1s, and H 1s states, which are consistent with the ionic bonding between the Li and the AlH4 unit.We found more mixing of the s and p states for Al atom in the β-LiAlH4 phase.Consequently, the electronic transition from the Al-s to -p state is related for the huge volume collapse during the α-LiAlH4 to β-LiAlH4 phase transition.

    Figure  2.  Phonon dispersion relations and density of state in high-symmetry directions for the α-LiAlH4 structure at 0 and 1.6 GPa
    Figure  3.  Calculated total and partial electronic densities of states for α-LiAlH4 structure at 0 GPa (a) and β-LiAlH4 structure at 1.6 GPa (b) respectively

    In order to better understand the bonding interactions between H, Al and Li atoms, we studied the Mulliken charges and the bond overlap population (BOP) P on the basis of Mulliken Populations in Table 2[19-20].The scaled BOP (Ps) is defined in PAl(Li)─Hs=PAl(Li)─H/LAl(Li)─H with L representing the average bond length of Al─H or Li─H, respectively[21].It is found that the PAl─Hs value for the α-LiAlH4 phase (0.506) is smaller than that of the β-LiAlH4 phase (0.534) as the coordination number of Al (4) stays unchanged at the α-β transition.Previous studies have shown that the smaller the BOP, the lower the hydrogen desorption kinetic energy[22].From this point of view, the current study suggests that the activation energy of the α-LiAlH4 is lower than that of the β-LiAlH4. The ambient phase of α-LiAlH4 is expected to be the most promising candidate for hydrogen storage.

    Table  2.  Average net charges, bond length (L) and scaled bond overlap population (Ps) between H, Al, and Li atoms in the α-LiAlH4 (at 0 GPa) and the β-LiAlH4 (at 2.0 GPa) structures
    Phase Average net charge L/nm Ps
    H Al Li Al─H Li─H Al─H Li─H
    α-LiAlH4 -0.48 0.64 1.29 0.161 5 0.188 6 0.506 0.019
    β-LiAlH4 -0.46 0.55 1.28 0.162 8 0.206 7 0.534 -0.063
     | Show Table
    DownLoad: CSV

    We investigated the pressure-induced phase transformations in LiAlH4 using the first-principles based on the density functional theory with the plane-wave basis.The structural transitions from α-LiAlH4 to β-LiAlH4 occurs at 1.6 GPa, accompanied with about 18% volume collapse, originating from the electronic transition of Al-s to -p states.The phase transition should be related to the pressure-induced phonon softening.The BOP analysis shows that the activation energy of the α-LiAlH4 is smaller than that of the β-LiAlH4. The α-LiAlH4 is expected to be the most promising candidate for hydrogen storage.

    Acknowledgements: This work was supported by the High Performance Computing Center of Yanbian University.
  • [1]
    SCHLAPBACH L, ZÜTTEL A.Hydrogen-storage materials for mobile applications[J]. Nature, 2001, 414(6861):353-358. doi: 10.1038/35104634
    [2]
    SCHÜTH F, BOGDANOVIC' B, FELDERHOFF M.Light metal hydrides and complex hydrides for hydrogen storage[J]. Chemical Communications, 2005, 36(2):2249-2258.
    [3]
    ZHU C Y, LIU Y H, DUAN D F.Structural transitions of NaAlH4 under high pressure by first-principles calculations[J]. Physica B:Condensed Matter, 2011, 406(8):1612-1614. doi: 10.1016/j.physb.2011.02.007
    [4]
    DILTS J A, ASHBY E C.Thermal decomposition of complex metal hydrides[J]. Inorganic Chemistry, 1972, 11(6):1230-1236. doi: 10.1021/ic50112a015
    [5]
    DYMOVA T N, ALEKSANDROV D P, KONOPLEV V N, et al.Spontaneous and thermal-decomposition of Lithium Tetrahydroaluminate LiAlH4-the promoting effect of mechanochemical action on the process[J]. Russian Journal of Coordination Chemistry, 1994, 20(4):279-285. http://www.academia.edu/13646109/Hydrogen_production_from_solid_reactions_between_MAlH4_and_NH4Cl
    [6]
    VAJEESTON P, RAVINDRAN P, VIDYA R, et al.Huge-pressure-induced volume collapse in LiAlH4 and its implications to hydrogen storage[J]. Physical Review B, 2003, 68:212101. doi: 10.1103/PhysRevB.68.212101
    [7]
    PITT M P, BLANCHARD D, HAUBACK B C, et al.Pressure-induced phase transitions of the LiAlD4 system[J]. Physical Review B, 2005, 72:214113. doi: 10.1103/PhysRevB.72.214113
    [8]
    CHELLAPPA R S, CHANDRA D, GRAMSCH S A, et al.Pressure-induced phase transformations in LiAlH4[J]. The Journal of Physical Chemistry B, 2006, 110(23):11088-11097. doi: 10.1021/jp060473d
    [9]
    TALYZIN A V, SUNDQVIST B.Reversible phase transition in LiAlH4 under high-pressure conditions[J]. Physical Review B, 2004, 70:180101. doi: 10.1103/PhysRevB.70.180101
    [10]
    HOHENBERG P, KOHN W.Inhomogeneous electron gas[J]. Physical Review B, 1964, 136(3):864-871.
    [11]
    KOHN W, SHAM L.Self-consistent equations including exchange and correlation effects[J]. Physical Review A, 1965, 140(4):1133-1138. http://tu-freiberg.de/sites/default/files/media/institut-fuer-theoretische-physik-10451/Lehre/Dichtefunktionaltheorie/a9rf1a3.pdf
    [12]
    PERDEW J P, WANG Y.Accurate and simple analytic representation of the electron-gas correlation energy[J]. Physical Review B, 1992, 45(23):13244-13249. doi: 10.1103/PhysRevB.45.13244
    [13]
    TROULLIER N, MARTINS J L.Efficient pseudopotentials for plane-wave calculations[J]. Physical Review B, 1991, 43(3):1993-2006. doi: 10.1103/PhysRevB.43.1993
    [14]
    SEGALL M D, LINDAN P J D, PROBERT M J, et al.First-principles simulation:ideas, illustrations and the CASTEP code[J]. Journal of Physics:Condensed Matter, 2002, 14:2717-2744. doi: 10.1088/0953-8984/14/11/301
    [15]
    MONKHORST H J, PACK J D.Special points for Brillouin-zone integrations[J]. Physical Review B, 1976, 13(12):5188-5192. doi: 10.1103/PhysRevB.13.5188
    [16]
    PARLINSKI K, LI Z Q, KAWAZOE Y.First-principles determination of the soft mode in cubic ZrO2[J]. Physical Review Letters, 1997, 78(21):4063-4066. doi: 10.1103/PhysRevLett.78.4063
    [17]
    TOGO A, OBA F, TANAKA I.First-principles calculations of the ferroelastic transition between rutile-type and CaCl2-type SiO2 at high pressures[J]. Physical Review B, 2008, 78(13):134106. doi: 10.1103/PhysRevB.78.134106
    [18]
    HAUBACK B C, BRINKS H W, FJELLVG H.Accurate structure of LiAlD4 studied by combined powder neutron and X-ray diffraction[J]. Journal of Alloys and Compounds, 2002, 346(1):184-189.
    [19]
    SEGALL M D, PICKARD C J, SHAH R, et al.Population analysis in plane wave electronic structure calculations[J]. Molecular Physics, 1996, 89(2):571-577. doi: 10.1080/002689796173912
    [20]
    MULLIKEN R S.Electronic population analysis on LCAO-MO molecular wave functions[J]. The Journal of Chemical Physics, 1955, 23(10):1833-1840. doi: 10.1063/1.1740588
    [21]
    HU C H, CHEN D M, WANG Y M, et al.First-principles investigations of the pressure-induced structural transitions in Mg(AlH4)2[J]. Journal of Physics:Condensed Matter, 2007, 19:176205. doi: 10.1088/0953-8984/19/17/176205
    [22]
    WANG H, LI Q, WANG Y C, et al.High-pressure polymorphs of Li2BeH4 predicted by first-principles calculations[J]. Journal of Physics:Condensed Matter, 2009, 21:385405. doi: 10.1088/0953-8984/21/38/385405
  • Relative Articles

    [1]ZHAO Jinggeng. Crystal Structure and Physica Properties of Perovskite Oxide BaMO3 (M Being Transition Metal)[J]. Chinese Journal of High Pressure Physics, 2024, 38(5): 050103. doi: 10.11858/gywlxb.20240753
    [2]ZHOU Xiao, SONG Shupeng, LIU Huiqi, LU Ze. First Principles Study on the Electronic Structure and Optical Properties of Graphene/MoS2 Heterojunctions with Different Rotation Angles[J]. Chinese Journal of High Pressure Physics, 2024, 38(5): 052201. doi: 10.11858/gywlxb.20240752
    [3]HE Xin, TIAN Hui, WANG Jian, CHEN Wanlei, WEI Zhaoxuan, LIU Jincheng, QI Dongli, SHEN Longhai. Density Generalized Function Theory Study on New MAX Phase M2SeC (M=Zr, Hf) under High Pressure[J]. Chinese Journal of High Pressure Physics, 2023, 37(4): 041102. doi: 10.11858/gywlxb.20230644
    [4]LI Zuo, LIU Yun, LIAO Dalin, CHENG Lihong. First-Principles Study on Structural, Electronic and Optical Properties of G2ZT Crystal under High Pressure[J]. Chinese Journal of High Pressure Physics, 2022, 36(4): 042202. doi: 10.11858/gywlxb.20220514
    [5]QIN Rui. New Developments of Time-Dependent Density Functional Theory and Its Applications[J]. Chinese Journal of High Pressure Physics, 2019, 33(3): 030101. doi: 10.11858/gywlxb.20190747
    [6]SUN Ying, ZHONG Xin, LÜ Jian, MA Yanming. Pressure-Induced Metallization and Novel Superconductivity of Chalcogen Hydrides[J]. Chinese Journal of High Pressure Physics, 2018, 32(1): 010102. doi: 10.11858/gywlxb.20170629
    [7]YAN Hua, LIU Zhi-Yong, SHI Mao-Sheng, LUO Yong-Feng, YAN Peng, GONG Xue-Dong. Theoretical Study on the Reaction of Chlorine Trifluoride Oxidewith Propylene Oxide by Density Functional Theory[J]. Chinese Journal of High Pressure Physics, 2015, 29(3): 232-240. doi: 10.11858/gywlxb.2015.03.011
    [8]HU Yong-Jin, HE Kai-Hua. Electronic Structure and Photoelectric Properties ofZnTe under High Pressure[J]. Chinese Journal of High Pressure Physics, 2014, 28(6): 641-647. doi: 10.11858/gywlxb.2014.06.001
    [9]HE Lin, YIN Jun. Effects of the Vacancy Point-Defect on Electronic Structure and Optical Properties of Diamond under High Pressure[J]. Chinese Journal of High Pressure Physics, 2013, 27(6): 802-806. doi: 10.11858/gywlxb.2013.06.002
    [10]DING Ying-Chun, LIU Hai-Jun, JIANG Meng-Heng, CHEN Min, CHEN Yong-Ming. First-Principles Investigations on Structural Transformation and Electronic Properties of BeP2N4 under High Pressure[J]. Chinese Journal of High Pressure Physics, 2012, 26(6): 674-680. doi: 10.11858/gywlxb.2012.06.012
    [11]WAN Jie, LIU Cheng-Mei, LI Ti, LIU Wei, . Effect of Dynamic High Pressure Microfluidization on the Crystal Structure of Dietary Fiber[J]. Chinese Journal of High Pressure Physics, 2012, 26(6): 639-644. doi: 10.11858/gywlxb.2012.06.007
    [12]WANG Xi-Cheng, GUO Jian-Yun, ZHENG Guang, HE Kai-Hua, CHEN Qi-Li, WANG Qing-Bo, CHEN Jing-Zhong. Density Functional Theory Studies of the Electronic and Optical Properties of Zinc Blende InN under High Pressure[J]. Chinese Journal of High Pressure Physics, 2012, 26(6): 653-660. doi: 10.11858/gywlxb.2012.06.009
    [13]HE Kai-Hua, ZHENG Guang, CHEN Gang, Lü Tao, WAN Miao, XIANG Dong, JI Guang-Fu. First Principles Study of the Elastic, Electronic and Optical Properties of CdO under Pressure[J]. Chinese Journal of High Pressure Physics, 2007, 21(3): 299-304 . doi: 10.11858/gywlxb.2007.03.014
    [14]Lü Meng-Ya, CHEN Zhou-Wen, WANG Wen-Kui, LIU Ri-Ping. Pressure Effect on Geometric Structure of Ge R8 Phase and Its Electronic Properties[J]. Chinese Journal of High Pressure Physics, 2006, 20(3): 291-295 . doi: 10.11858/gywlxb.2006.03.012
    [15]SONG Gong-Bao, HAN Cang-Qiong, LIANG Jing-Kui, RAO Guang-Hui. The Extent and Crystal Structure of the Solid Solution Based on Ca2Pr2Cu5O10[J]. Chinese Journal of High Pressure Physics, 2004, 18(1): 27-35 . doi: 10.11858/gywlxb.2004.01.006
    [16]SONG Gong-Bao, LIANG Jing-Kui, RAO Guang-Hui. Subsolidus Phase Relation and Crystal Structure in the Pr1-xBa2-yCax+yCu3O7 System[J]. Chinese Journal of High Pressure Physics, 2004, 18(4): 309-314 . doi: 10.11858/gywlxb.2004.04.004
    [17]ZHANG Chun-Bin, LI Shao-Meng, TENG Ya-Gang. Electronic Structures and Equation of State for Aluminium under High Pressure[J]. Chinese Journal of High Pressure Physics, 1994, 8(4): 272-278 . doi: 10.11858/gywlxb.1994.04.005
    [18]XIE Yuan-Nan. Electronic Structure and EOS of -SiC[J]. Chinese Journal of High Pressure Physics, 1994, 8(1): 57-59 . doi: 10.11858/gywlxb.1994.01.009
    [19]WANG Xiao-Bing, TIAN De-Cheng. Electronic Structure of TiB2 under High Pressures[J]. Chinese Journal of High Pressure Physics, 1993, 7(1): 75-80 . doi: 10.11858/gywlxb.1993.01.012
    [20]JIN Qing-Hua, DING Da-Tong, WANG Ding-Sheng, LI Guang-Wei. Equation of State and Electronic Transfer for Al, Cu, Pb and Ta under Ultrahigh Pressures[J]. Chinese Journal of High Pressure Physics, 1990, 4(1): 17-23 . doi: 10.11858/gywlxb.1990.01.003
  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(3)  / Tables(2)

    Article Metrics

    Article views(7042) PDF downloads(152) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return