储能科学与技术, 2023, 12(11): 3287-3298 doi: 10.19799/j.cnki.2095-4239.2023.0488

储能材料与器件

从材料结构特征到功能特征——

一种新型锂离子导体设计方法

赵悠曼,1, 黄旸1,2, 熊立1, 林海军1

1.东莞市创明电池技术有限公司,广东 东莞 523808

2.中国科学院重庆绿色智能技术研究院,重庆 400714

From material structural feature to the functional feature

—A new designing rule for lithium-ion conductors

ZHAO Youman,1, HUANG Yang1,2, XIONG Li1, LIN Haijun1

1.Chuangming Battery Technology Co. Ltd. , Dongguan 523808, Guangdong, China

2.Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences, Chongqing 400714, China

通讯作者: 赵悠曼(1984—),女,硕士,工程师,主要从事锂离子电池的开发及应用,E-mail:youymanm@163.com

收稿日期: 2023-07-17   修回日期: 2023-08-02  

基金资助: 国家自然科学基金青年项目.  双阳离子氧化还原中心的富锂正极材料的理论设计及多谱学实验研究12004382

Received: 2023-07-17   Revised: 2023-08-02  

作者简介 About authors

赵悠曼(1984—),女,硕士,工程师,主要从事锂离子电池的开发及应用,E-mail:youymanm@163.com。 E-mail:youymanm@163.com

摘要

近年来,有理论研究认为硫系锂离子导体材料的硫阴离子框架都具有类bcc (body-centered-cubic)‍的体心立方堆叠结构。因此,在探索新的硫系锂离子导体时,应从材料库中优选那些具有类bcc硫阴离子框架结构的材料。然而对于设计新型锂离子导体而言,这种设计准则过于简略。本团队基于前人工作基础,对硫系化合物中锂离子的传输和扩散机制进行了研究,并讨论了影响锂离子局域跳跃的结构环境和影响锂离子宏观扩散的化学环境。研究发现,在硫系化合物中,不与过渡金属配位的0-TM通道对于锂离子在晶位之间的跳跃具有决定性的作用;而要形成宏观扩散网络,晶体结构中必须含有足够多的0-TM通道。基于这一认识,本团队探索并设计了几种新的硫系锂离子导体,验证了硫系化合物结构特征与其功能特征之间的关联性。这一工作深化了bcc结构设计准则,为发展新型锂离子导体提供了更为明确的理论方向。

关键词: 锂离子导体 ; 离子框架 ; 锂离子扩散 ; 扩散系数

Abstract

In recent years, body-centered-cubic (bcc)‍-type anion framework has been considered as a deft structural design descriptor for identifying unexplored potential sulfide lithium-ion conductors (SICs). However, the bcc guidance must be completed, and other factors impacting lithium diffusion must be considered. Therefore, they are not always useful in the design of new lithium-ion conductors. In this study, based on previous works, we investigated the lithium-ion hopping and diffusion mechanisms in sulfides. We discussed how the local structural environment impacts lithium hopping from site to site in a micro-manner and how the chemical composition impacts diffusion in a macro-manner. In conclusion, we proved that the zero-transition-metal (0-TM) coordinated channel is essential for the facile hopping of lithium ions, and macro-diffusion can be expected in the materials only if the contents of 0-TM channels exceed the penetrating thresholds. Our study provides a more thorough interpretation of the sulfides' lithium diffusion mechanism and directly connects the structural and functional features of the material.

Keywords: lithium ionconductors ; ion framework ; lithium ion diffusion ; diffusion coefficient

PDF (3658KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

赵悠曼, 黄旸, 熊立, 林海军. 从材料结构特征到功能特征——[J]. 储能科学与技术, 2023, 12(11): 3287-3298

ZHAO Youman. From material structural feature to the functional feature[J]. Energy Storage Science and Technology, 2023, 12(11): 3287-3298

锂离子电池因其优良的性能,在储能领域得到了广泛的应用[1]。目前绝大部分商用锂离子电池采用的都是可燃性的有机电解液,存在着一定的燃爆风险[2]。随着锂离子电池能量密度和容量的不断升高,有机电解液体系的安全性问题将会变得更加突出[3]。固态电解质由于其极高的燃点、更宽的电化学稳定窗口以及工作温度区间,受到了学术界和产业界的广泛关注,被认为是解决电池安全性问题和进一步提高其性能的重要方案[4]。近年来,人们花费了大量的精力来寻找具有高室温锂离子电导率(>10-1 mS/cm)的材料以及体系,以研究其作为锂离子电池固态电解质的性能[5-8]。在所有候选材料体系中,硫系锂离子导体材料由于其较高的锂离子电导率,具有很好的应用前景[9-11]。近期报道的一些硫系化合物如Li7-x PS6-x Cl x (LPSC)[12],Li10GeP2S12(LGPS)[13]和其Si/Sn掺杂衍生化合物[14-15],Li7P3S11(LPS)[16]都具有很高的室温锂离子电导率,甚至可以比拟液态电解质。

为了更好地探索锂离子导体,人们提出了一些结构设计准则[17-19]。有研究者认为在具有类bcc排布的硫阴离子框架的硫化物中,锂离子具有更低的扩散势垒,因此这种化合物更有可能成为锂离子导体[17]。一些典型的硫系锂离子导体,如LGPS和LPS的硫阴离子框架都具有类bcc的堆叠特征。根据这一理论,人们认为同样具有硫阴离子bcc堆叠结构的LZPS应该也是锂离子导体[20]。然而,理论计算和实验分析均表明,LZPS中锂离子具有很高的扩散势垒,并不满足锂离子导体的要求。而进一步的理论计算则表明,通过在Zn位上引入Li,形成Li-Zn反位缺陷,可以极大地降低锂离子的扩散势垒,使LZPS满足锂离子导体的要求。在这一理论的指导下,人们合成了一系列具有Li-Zn反位缺陷的LZPS富锂相衍生化合物(Li1+2x Zn1-x PS4),并证明它们具有较高的锂离子电导率[21-22]。LZPS的低锂离子电导性和Li1+2x Zn1-x PS4的高锂离子电导性充分说明材料中锂离子的传输不仅受阴离子框架结构的影响,而且受阳离子框架结构的影响。因此,探索材料中阴阳离子的结构排布对于锂离子扩散的影响,对设计新型锂离子导体具有重要理论指导意义。

因此,本工作以LZPS为对象,研究LZPS的整体结构特征以及影响锂离子在晶位之间跳跃的局域结构环境,并利用密度泛函理论(DFT)计算不同局域结构环境下锂离子的迁移势垒。并以此为基础,分析了LZPS是非锂离子导体而其富锂Li1+2x Zn1-x PS4是锂离子导体的原因,提出一种更为明确的硫系锂离子导体的结构设计准则,并依据该结构设计准则,设计了几种新型硫系锂离子导体材料。

1 实验方法

1.1 密度泛函理论(DFT)计算

本工作用VASP软件进行DFT计算[23-25]。采用投影缀加平面赝势(PAW)描述原子核及内层电子之间的相互作用[26-27]。采用广义梯度近似(GGA)里的PBE泛函描述电子的交换-关联作用[28-29]。平面波的截断能为400 eV。在进行几何结构优化时,离子步的收敛标准为0.01 eV/Å,电子步的收敛标准为10-5 eV。

本工作用climb image nudged elastic band (CI-NEB)方法搜索材料中锂离子的最佳扩散路径[30-31]。搜索路径的起点和终点事先经过了结构优化,以保证结果的准确性。搜索路径上插入的帧数由起点和终点之间的距离决定。

对于理论预测的 ε -Li3PS4,本工作用Phonopy软件包计算其声子谱,以判断其结构在动力学上的稳定性[32]。该软件利用直接法(frozen-phonon)在优化后的平衡结构中引入原子位移,计算作用在原子上的Hellmann-Feynman力,进而由动力学矩阵计算出声子色散曲线。这种方法需要构建材料晶体结构的超胞,通常超胞内的原子数要大于100个。P-T相图描述的是材料在不同温度和压强下的吉布斯自由能变化,通过对其进行计算可以判断材料不同相结构随着压强和温度的变化情况。本工作利用Phonopy软件包里的准简谐近似方法(quasi harmonic approximation,QHA)计算材料在不同温度和压强下的吉布斯自由能[33]

本工作利用第一性原理分子动力学(abinitio molecular dynamics,AIMD)计算材料中锂离子的扩散系数。为了进行分子动力学模拟,通常需要建立包含原子数超过100的超胞。为了提高计算的效率,在K空间中只取Γ一个网格点。采用Verlet算法对牛顿运动方程进行积分求解[34],积分的步长为2 fs。分子动力学模拟的初始结构经过第一性原理结构优化,体系位于势能曲面的局域最低处。体系的初始温度设为100 K,体系内原子的初始速度符合玻尔兹曼分布。接着使用速度标度(scaling velocities)的方法将体系加热到设定的温度,并且保证加热的总步长不小于1000步(2 ps)。最后采用正则系综(NVT)对体系在设定温度附近进行热平衡,通过Nosé-Hoover恒温器保持体系温度维持在设定温度附近[35-36]

1.2 富锂相体系的结构计算

富锂相Li1+2x Zn1-x PS4中存在着Li和Zn共占位的现象,在实际的结构中Li和Zn在晶位上的占据是随机的,只需满足其化学计量比即可。但第一性原理的材料结构计算必须是针对有序结构进行,即结构模型中所有晶位上占据的原子类型和原子数目都是确定且唯一的。本工作利用pymatgen软件包模拟出富锂相下的有序子结构以进行结构计算。其基本原理是利用枚举法对存在共占位的晶位上的原子进行枚举,进而获得大量的有序结构实例,接着利用Ewald方法计算这些结构实例的静电能,并根据静电能的大小对其进行排序,最终得到静电能最低的10个结构。再利用DFT计算这10个结构的精确能量,得到能量最低的结构,并认为这个结构是富锂相体系的稳定结构,用于后续的NEB计算和AIMD模拟。

1.3 体系中0-TM扩散通道的渗透阈值计算

本工作中,材料构成锂离子导体的关键在于材料中0-TM微观扩散通道在材料中形成连续的扩散通路,从而支持锂离子在材料体相内持续传输。针对硫系化合物,这一问题可以简化成一个数学问题,即在无限大小的空间网格点中,这种网格点被0-TM通道占据的概率为x(x的大小由材料中锂离子的含量xLi决定),被非0-TM通道占据的概率为1-x。求解当xLi等于多少时,这些0-TM通道能够在空间网格点上形成一个连续的扩散网络。这一问题可以采用渗透理论的相关算法进行求解[37, 38]。0-TM通道形成渗透扩散网络的概率Pw是一个阶跃函数:

PwxLi=  0    if xLi< xc  1    if xLi > xc

xc为渗透阈值,即当材料中锂离子的含量超过xc时,材料中的0-TM通道一定能构成连续的扩散网络,材料成为锂离子导体。本工作采用Newman和Ziff所提出的算法来求解材料的0-TM扩散通道渗透阈值[39]。在实际计算中,由于模型的尺寸限制,Pw并不是一个标准的阶跃函数,而是S形函数。但随着计算模型尺寸的增加,Pw函数的形状越来越接近于阶跃函数。在实际计算中,我们进行了收敛测试,发现S形曲线的形状随着计算模型的尺寸增加迅速接近于阶跃函数的形状,其拐点可以认为是阶跃函数的变化点xc。最终,通过对计算量大小和计算结果精度的权衡,我们选定在10×10×10的超胞上进行渗透阈值xc的计算。

2 结果讨论

2.1 LZPS的结构特征以及锂离子在其中的局域扩散环境

LZPS所属的空间群为I4¯,其晶体结构如图1(a)、(b)所示。由PS4四面体构成非金属层,而LiS4四面体和ZnS4四面体构成金属层。非金属层和金属层沿着c轴交替排布,形成材料的整体框架。在金属层内,每个LiS4与四个ZnS4以角对角的方式相连。同样地,每个ZnS4也与四个LiS4以角对角的方式相连。在非金属层,PS4四面体之间并不直接相连,留下间隙位作为锂离子的扩散位点。

图1

图1   (a), (b) LZPS的晶体结构示意图 (LiS4ZnS4PS4 以及间隙位四面体分别用深绿色、蓝色、青色以及粉色表示,硫阴离子用黄色圆球表示;(c)为图 (b) 中红色虚线框标识部分结构的放大图 (其表示了锂离子跳跃到邻近间隙位的局域结构环境); (d)NEB计算得到的锂离子在图 (c) 所示的局域结构中的迁移路径

Fig. 1   (a), (b) The unit cell of LZPS viewed from different directions. The LiS4, ZnS4, PS4 and interstitial site are colored by dark green, blue, cyan and pastel orange. The sulfur is labeled with yellow spheres. (c) The local environment involving the lithium hopping in the LZPS structure which is extracted from the red box in (b). (d) The diffusion path in the octahedral cage determined by NEB calculation


图1(c)展示了LZPS中锂离子扩散的局域化学环境。其主体由一个硫八面体孔洞构成。这个硫八面体孔洞则由一个LiS4四面体,一个ZnS4四面体,一个PS4四面体以及一个四面体间隙位构成。直观看,四面体(T)内的锂离子可以通过硫八面体(O)迁移到相邻的四面体(T)间隙位上,构成T-O-T的迁移机制,这种机制常见于层状氧化中(例如LiCoO2、三元NMC材料等)。但是,由于八面体可以进一步地分成4个子四面体。因此锂离子在运动时,也有可能仅仅是通过这些子四面体进行迁移,形成T-T的迁移机制。在下文中,我们将利用NEB计算LZPS中锂离子的具体迁移路径[40]

2.2 LZPS中锂离子的扩散机制

图1(d)展示了用NEB方法计算得到的锂离子在硫八面体孔洞(6个顶点分别为S1、S2、S3、S4、S5、S6)中的具体迁移路径。很明显,锂离子并没有经过硫八面体的中心位置。从其迁移路径可以看到,其先经过硫四面体T1(四个顶点分别为S1、S2、S3、S4)然后再经过硫四面体T2(四个顶点分别为S1、S2、S3、S5),最后达到间隙位。因此,LZPS中锂离子的扩散是T-T机制,而不是T-O-T机制。图2(a)中展示了NEB计算的锂离子沿着图1(d)所示的路径进行迁移时所受到的势垒曲线。从图中可以看到,锂离子的最大迁移势垒Ea为0.43 eV,远大于常见硫系锂离子导体中锂离子的迁移势垒(LPS中Ea为0.18 eV)[16],LGPS及其掺杂衍生化合物中Ea为0.22~0.25 eV[13-15]。根据原子扩散系数σ和迁移势垒Ea之间具有e指数幂的关系[σ∝exp(-Ea/kT)],可以粗略地计算得到,在常温下LZPS中锂离子的扩散系数大约比常见的硫系快锂离子导体中锂离子的扩散系数低3个数量级。因此,LZPS是非锂离子导体。

图2

图2   (a) 蓝色实线表示用NEB计算得到的LZPS中锂离子迁移到邻近间隙位的迁移势垒曲线;绿色实线表示在富锂相Li1+2x Zn1-x PS4 中锂离子迁移到邻近间隙位且配位的锌离子被空位取代时的迁移势垒曲线;橘色实线表示在富锂相Li1+2x Zn1-x PS4 中锂离子直接迁移到邻位锌离子空位的迁移势垒曲线;(b) 在硫阴离子严格遵循bcc排布的结构中,锂离子的迁移势垒曲线

Fig. 2   (a) The calculated energy landscape of diffusion path from the Li site to: the interstitial in the stoichiometry LZPS (blue color). The interstitial in a Li-Zn substitution local environment (green color). The Zn site in the same layer (orange color). (All the migration paths are scaled to the same length.) (b) The calculated energy landscape of diffusion path from the Li site to the interstitial site in a perfect bcc sulfur framework


从以上计算可知,LZPS中虽然具有bcc排布的硫阴离子框架,但锂离子在其中的扩散势垒仍很高,其在常温下是贫离子导体。这与文献中具有bcc硫阴离子框架结构的材料是快离子导体的指导理论不相符。这是因为前述的理论中没有考虑到体系中阳离子对于锂离子扩散的影响。为了定量区分材料中硫阴离子和阳离子对于锂离子扩散的影响,我们设计了一个特殊的结构模型。在这个模型中,只有硫阴离子,没有阳离子。而且硫阴离子均位于一个标准的bcc晶胞的格点上,晶胞的体积为n×37.08 Å3n为晶胞中硫阴离子的数量。37.08 Å3称为自由扩散体积,其为LZPS晶胞中单个硫阴离子所占据的空间体积。通过这一换算,可以保证设计的特殊结构模型中,锂离子受到的硫阴离子的作用等效于LZPS中锂离子受到的硫阴离子的作用。在bcc晶胞中,锂离子的最稳态位置为硫四面体中心,其沿着T-T的跳跃路径,从一个硫四面体中扩散到相邻的硫四面体中心。图2(b)展示了用NEB方法计算得到的bcc结构模型中锂离子的扩散势垒曲线。可以看到,其最大扩散势垒Ea仅为0.14 eV,比LZPS中锂离子实际的最大扩散势垒要低很多。虽然LZPS中硫阴离子所组成的类bcc堆叠结构允许锂离子在结构中快速扩散,但由于阳离子的存在,其与锂离子之间强烈的排斥作用使锂离子在结构中的扩散势垒大大增加(从0.14 eV增加到0.43 eV),严重阻碍了锂离子的运动。

前人的工作[20]显示,通过在LZPS中Zn位进行锂元素掺杂,可以形成具有快离子导体特征的富锂相Li1+2x Zn1-x PS4,且随着掺杂量的增加,锂离子的扩散系数迅速增加,甚至可以超过目前已知的几种硫系快离子导体。可以看到,这一材料改性方法只针对阳离子框架,而硫阴离子框架的成分和结构并没有改变。这也进一步论证了我们的观点,阻碍LZPS中锂离子扩散的主要因素来自于阳离子之间的相互作用。通过对LZPS的阳离子框架进行修饰和改性,降低阳离子对于锂离子的束缚作用,可以让LZPS从贫离子导体变成快离子导体。在此基础上,我们对富锂相Li2.0Zn0.5PS4中锂离子的传输机制做了更深入的分析,通过与LZPS中锂离子的扩散机制进行对比,来探究这种掺杂改性降低锂离子扩散势垒的结构和化学因素。

在LZPS的富锂相中,锂离子掺杂到了锌位以及间隙位,在这两类位点上形成了无序的排布,因此在富锂相中,锂离子的扩散途径及机制较为复杂。在本文中,我们基于枚举法研究了Li2.0Zn0.5PS4中几种典型局域结构中锂离子的扩散机制。例如图1(d)的局域结构中的锌离子在富锂相中,可能被锂离子代替,形成LiZn反位缺陷。在这一局域结构中,锂离子从初始位置跳跃到相邻的间隙位时,对其发生排斥作用的阳离子从原来的锌离子变成了锂离子。由于锂离子之间相对较弱的静电相互作用,因此这一发生跳跃的锂离子的扩散势垒会进一步降低,如图2(a)中绿线所示,其最大扩散势垒降低到0.32 eV。在富锂相中,还存在着另外一种跳跃机制,即锌位的锌离子被锂离子取代,形成LiZn反位缺陷后,LiZn可以跳跃到同一层的锂位,如图2(a)中的橙线所示,这一扩散机制的最大扩散势垒为0.21 eV。值得一提的是在LZPS中,锂离子在Li-Zn层内的扩散是无法直接进行的,因为每个Li位的周围被4个Zn位填满,没有可供跳跃的空位,其只能沿着c轴跳跃到非金属层的间隙位后再跳跃到另一层的锂空位,可以看到在LZPS中锂离子的扩散是一维的。而在其富锂相中,由于部分LiZn反位缺陷的存在,锂离子可以在Li-Zn层内直接扩散,并且这种扩散的势垒极低。由于层内扩散的存在,富锂相内锂离子的扩散从一维变成三维。由于扩散系数和扩散网络的维度之间存在着指数的关系[41],因此富锂相中锂离子的扩散系数大大提升。同时,LiZn反位缺陷连接了相邻晶位上的锂离子,促进了锂离子跳跃的协同作用(相邻晶位上的锂离子在极短的时间内依次向邻近的空位上进行跳跃)。根据研究,在典型的快离子导体,如LGPS、LLZO和LATP中,由于协同作用的存在,可以使锂离子的扩散势垒降低0.2~0.3 eV[42]

2.3 富锂相中扩散网络的渗透

快离子导体必须在结构上满足两个条件:①材料中存在微观扩散通道,锂离子在通道的始态和终态位置上进行跳跃的势垒足够低;②这些通道必须互相连接并且贯穿材料的整个框架,以便锂离子能够在材料中连续跳跃(图3)。因此,尽管我们在2.2节计算得知,LZPS的富锂相中存在着具有低势垒的扩散通道,但如果材料中这种扩散通道的含量不足以形成扩散网络,就不能成为快离子导体。这一现象在阳离子无序岩盐相材料中也存在[43-46]。在这类材料中,锂离子通过氧四面体通道在氧八面体晶位之间跳跃,形成O-T-O的跳跃机制。计算表明,只有当氧四面体通道的周围没有过渡金属(TM)存在时,锂离子在其中跳跃的势垒才足够低,可以形成快速扩散。这种没有过渡金属配位的通道被称作0-TM通道。由于0-TM通道周围配位的都是锂离子,因此要增加0-TM通道含量使得其能在材料内形成互相连接的扩散网络,就必须增加材料中锂离子的含量,形成富锂相。前人的计算[43]显示,要保证阳离子无序岩盐相材料中渗透扩散网络的存在,至少要保证9%的过锂量。

图3

图3   富锂相Li1+2x Zn1-x PS4中锂离子扩散通道的分布示意图 [红色方块表示与锌离子配位的硫八面体通道,锂离子无法从这种通道内通过;黄色方块表示不与锌离子配位的硫八面体通道 (0-TM通道);在左图内,0-TM通道没有形成贯穿材料整体的渗透扩散网络,锂离子无法在材料内形成连续性扩散]

Fig. 3   The skecth map of a non-percolated network (left) and a percolated network (right)


上文对LZPS富锂相中锂离子局域扩散机制的计算表明,只有当硫八面体的周围没有锌离子配位时,锂离子经过硫八面体跳跃到邻近的间隙位,LiZn位的势垒才足够低。类比于阳离子无序岩盐相材料中的工作,我们视这种没有锌离子配位的硫八面体为具有0-TM通道。可以看到,在LZPS中也存在着扩散网络的渗透问题,说明只有当锌位上掺杂的锂离子数量够多形成足够量的0-TM通道时,这些0-TM通道相互连接构成贯穿整个材料体系的扩散网络,才能够保证锂离子的持续、大范围扩散。我们在以前的工作[47]中编写了一个计算阳离子无序岩盐相材料中锂离子扩散网络渗透阈值的蒙特卡洛程序。现在只需对这个程序进行很小部分的修改就可移植到LZPS中的锂离子扩散网络的渗透阈值的计算上。计算的结果如图4所示(为了更好地说明计算结果,我们在图例中给出的富锂相化学式是LixLiZn1.5-xLi/2PS4,而不是前文中一直用的Li1+2x Zn1-x PS4)。计算表明,当x=1.48时,LZPS富锂相中的0-TM通道能够形成渗透的扩散网络。也就是说,我们在考虑对LZPS进行掺杂以提高其对锂离子的导通时,至少要掺杂0.48份化学计量数的锂离子,减少0.24份化学计量数的锌离子。我们的计算结果也得到了实验数据的支持。对合成的LZPS富锂相材料的锂离子扩散系数的测量表明,当x大于1.8时材料中锂离子的扩散系数才会有明显的提升[21]。该计算结果略小于实验数据。但在我们的计算结果中,当材料的组分为Li1.48Zn0.76PS4时,虽然材料中具有渗透的扩散网络,但在总计1.48份的锂离子中只有0.25份锂离子与这些渗透网络相连,具有扩散活性。其余的锂离子只在局域位置进行扩散,对于实验所观测到的材料扩散系数并没有贡献。因此,这些富锂组分较低的锂离子扩散系数测量结果是可以解释的。而在Li1.8Zn0.6PS4中,总计1.8份的锂离子中有1.4份的锂离子是具有扩散活性的,因此其具有较高的锂离子扩散系数。

图4

图4   利用蒙特卡洛算法计算得到的富锂相Li1+2x Zn1-x PS4(图中用对锂含量归一化的LixLiZn1.5-xLi/2PS4表示)中0-TM通道在材料中形成渗透扩散网络的概率与材料中锂离子含量之间的关系 (红色实线);以及与渗透网络相连的可形成全局性扩散的锂离子含量

Fig. 4   Periodic wrapping probabilities (red color line) and percolation thresholds (indicated by vertical bars). Accessible amount of lithium connecting to the percolating network (blue color line)


2.4 LZPS的富锂极限相 ε-Li3PS4

通过前述的讨论可以知道,LZPS中锌离子是阻碍锂离子扩散的主要因素。通过在锌位上进行锂离子掺杂形成富锂相,可以提高材料中锂离子的扩散能力。前人的研究表明,Li1+2x Zn1-x PS4中锂离子的掺杂量越多,材料锂离子电导率越高[20]。因此,本文对Li1+2x Zn1-x PS4的极端情况Li3PS4(x=1)进行研究。文献研究表明,Li3PS4这一材料组分在不同的烧结温度和压强下可以形成不同的结构,分别为 α -Li3PS4β -Li3PS4γ -Li3PS4δ -Li3PS4[48-49]。其中, α -Li3PS4β -Li3PS4γ -Li3PS4在常压下随着烧结温度的提高,可以互相转化( α -Li3PS4具有最低的烧结温度)。 δ -Li3PS4是通过理论计算发现的,其只有在高压(>5 GPa)下才能生成。根据相的命名规则,本文中的Li3PS4被命名为 ε 相。对 ε -Li3PS4声子谱的计算表明,这种结构在动力学上是稳定的(声子谱上没有虚频)。而要判断其是否能够在实验室内被合成,还需要对其进行热力学稳定性的计算。由于在常压的材料合成实验中并没有发现 ε -Li3PS4,因此我们推断其必须要在高压下才能合成。因为 γ -Li3PS4也是一种在高压下合成的相结构,我们计算了 γ -Li3PS4ε -Li3PS4的吉布斯自由能,并在不同压力下对它们进行了对比。如图5所示,可以看到在常压下(约等于0 GPa), ε -Li3PS4相在温度低于255 K (-18 ℃)时,其自由能都小于 γ -Li3PS4相。这就意味着在常压下,当合成温度低于255 K时,组分Li3PS4合成的是 ε 相。但这一合成温度过低(低于室温),在材料合成时需要对前驱体材料进行加热等处理以打破原有的晶体结构,一般硫系化合物的合成温度都在373 K(100 ℃)以上,显然在常压下 ε 相是无法合成的。而随着压强的增加, ε -Li3PS4γ -Li3PS4 的转化温度也随之升高。如图5所示,在5 GPa的合成压强下,这一转化温度为315 K (42 ℃),在10 GPa的合成压强下,这一转化温度为445 K (172 ℃)。这说明 ε -Li3PS4在实验中是可以合成的,但其合成条件比 γ -Li3PS4更为苛刻,需要更高的压强。

图5

图5   (a) ε-Li3PS4 的声子谱 (可以看到其上没有虚频)(b), (c), (d) 在常压、5 GPa以及10 GPa的外界条件下 ε-Li3PS4γ-Li3PS4 的吉布斯自由能变化

Fig. 5   (a) The calculated phonon band structures for ε-Li3PS4; (b), (c), (d) The calculated Gibbs free energy of ε-Li3PS4 and γ-Li3PS4 at different pressures


同时,我们用AIMD方法计算了 ε -Li3PS4的室温锂离子电导率。如图6所示, ε -Li3PS4的室温锂离子电导率可以达到207.41 mS/cm,扩散势垒仅为0.132 eV。前人对LZPS富锂相的常温锂离子电导率的计算结果为 σ (Li1.5Zn0.75PS4)=27.7 mS/cm, σ (Li2.0Zn0.5PS4)=53.8 mS/cm, σ (Li2.5Zn0.25PS4)=114 mS/cm。我们对 ε -Li3PS4的计算结果也符合前人计算结果所表现出来的趋势,锌的含量越低,锂离子的扩散系数越高。并且,与Li3PS4的其他相比起来, ε -Li3PS4的锂离子具有最高的导锂离子能力[48, 50-52]。从前人的研究结果看, ε -Li3PS4的锂离子电导率在所有目前已知的材料中是最高的。因此,尽管我们的计算表明这种材料在实验中合成的难度较大,但其仍然值得人们进一步研究。

图6

图6   AIMD模拟所得到的 ε-Li3PS4,以及LZGSLZSS富锂相的阿伦尼乌斯扩散系数拟合图

Fig. 6   Arrhenius plots for the ε-Li3PS4 and the composition engineered LZSS, LZGS obtained by the AIMD simulation


2.5 基于材料结构特征的硫系快离子导体设计

在前面的章节中,我们对LZPS中锂离子的扩散机理进行了深入的分析和讨论,并解释了为何LZPS是贫离子导体而其富锂相却是快离子导体。这种转变的关键在于阳离子排布结构上的变化。当锌离子被移除后,在锂离子的迁移路径上可以形成有利于其快速通过的0-TM通道。因此,在材料中锂离子的扩散不仅受到阴离子框架的制约,也受到阳离子框架的制约,在考虑快离子导体的结构设计准则时,应该将阴阳离子的排布特征同时考虑进去。因此,我们提出了新的硫系快离子导体的结构设计准则。

(1)首先判断候选材料的硫阴离子框架是否具有bcc的排布特征。在这种结构中,锂离子遵循T-T的跳跃机制,硫阴离子对于锂离子跳跃的阻碍最低。

(2)移除材料中的部分过渡金属阳离子,降低过渡金属阳离子与锂离子之间的静电排斥,以在材料中形成有利于锂离子扩散的0-TM通道。

(3)在满足材料热力学和动力学稳定的前提下,移除多的过渡金属阳离子,以形成足够量的0-TM通道,在材料中形成渗透扩散网络,让锂离子能够在材料中形成持续性大范围扩散。

为了进一步验证我们提出的结构设计准则的可行性,我们将其运用到两种新的材料 β -Li2ZnGeS4 (LZGS)和Li2ZnSiS4 (LZSS)上[53- 54],这两种材料从未被作为快离子导体研究过。这两种材料的空间群分别为Pn和Pna21,锂、锌、锗和硅原子都位于互相顶角连接的硫四面体中心。其晶体结构如图7所示。我们的结构匹配算法显示,这两种材料的硫阴离子都具有类bcc的排布(图8)。因此,它们满足硫系快离子导体设计准则的第一条。接下来对阳离子结构排布进行改造。由于LZGS和LZSS较低的结构对称性,其锂离子迁移的局域结构环境更加复杂。如图7所示,LZGS中与锂离子迁移相关的局域结构可以分为三种,第一种是由两个LiS4四面体,一个PS4四面体以及一个空四面体围合而成的硫八面体孔洞,我们用Oct(2LiS4-PS4)来表示。同样地,Oct(2LiS4-ZnS4)表示两个LiS4四面体,一个ZnS4四面体以及一个空四面体围合而成的硫八面体孔洞,而Tet(2LiS4-ZnS4-PS4)则表示两个LiS4四面体,一个ZnS4四面体以及一个PS4四面体围合而成的硫四面体孔洞。

图7

图7   (a), (b) LZGSLZSS的晶胞结构示意图 (LiS4ZnS4GeSLiSi4 用深绿色、蓝色、天青色以及藏青色四面体表示)(c), (d), (e)为图 (a) 中红色箭头所示的锂离子三种迁移路径上的局域结构环境的放大图

Fig. 7   The unit cell of LZGS (a) and LZSS (b). The LiS4, ZnS4, GeS, and LiSi4 are colored by dark green, blue, pastel cyan and dark cyan. (c), (d), (e) The local environment involving the lithium hopping in the LZGS structure marked as path 1, path 2 and path 3 in (a)


图8

图8   用图像匹配算法得到的与LZGSLZSS的硫阴离子框架最匹配的bcc排布结构 (图中的黄色圆球表示材料中硫阴离子的实际晶体位置,红色圆球表示经过图像匹配算法调整后的最符合bcc排布的硫阴离子位置)

Fig. 8   Ideal body-centered-tetragonal (bct) lattice (red) overlaid on the sulfur framework (yellow) of β-Li2ZnGeS4 (a) and Li2ZnSiS4 (b). The root-mean-square (rms) distances from the S atoms in the β-Li2ZnGeS4 and 0.56 Å in the Li2ZnSiS4 to the idealized bcc position are 0.53 Å and 0.56 Å respectively


NEB计算显示,锂离子在Oct(2LiS4-PS4)内的扩散势垒为0.27 eV,在Oct(2LiS4-ZnS4)内的扩散势垒为0.39 eV,在Tet(2LiS4-ZnS4-PS4)内的扩散势垒为0.53 eV(图9)。

图9

图9   利用NEB计算得到的LZGS中锂离子在Oct(2LiS4-PS4) 孔洞内的迁移势垒曲线 (橙色实线);在Oct(2LiS4-ZnS4) 孔洞中的迁移势垒曲线 (蓝色实线);在Tet(2LiS4-ZnS4-PS4) 中的迁移势垒曲线 (绿色实线)。在LZGS富锂相中,部分锌离子被取代形成空位后,锂离子在Oct(2LiS4-VacS4) 孔洞中的迁移势垒曲线 (黑色实线) 以及在Tet(2LiS4-VacS4-PS4) 孔洞中的迁移势垒曲线 (红色实线)

Fig. 9   The calculated energy landscape of diffusion path in the LZGS from the Li site to the adjacent vacancy through the Oct(2LiS4-PS4) cage (orange color). The Oct(2LiS4-ZnS4) cage (blue color). The Tet(2LiS4-ZnS4-PS4) cage (green color). The Oct(2LiS4-VacS4) cage (black color). The Tet(2LiS4-VacS4-PS4) cage (red color)


显然,只有Oct(2LiS4-PS4)支持锂离子的快速扩散。这一计算结果和LZPS富锂相的计算结果一致,锂离子只有在不与过渡金属阳离子配位的硫八面体孔洞内(即0-TM通道)才具有较低的扩散势垒。尽管LZGS中存在0-TM通道,符合我们设计准则的第二条,但这些0-TM通道并不互相连接,一个Oct(2LiS4-PS4)硫八面体孔洞仅与两个Oct(2LiS4-ZnS4)硫八面体孔洞以及两个Tet(2LiS4-ZnS4-PS4) 硫四面体孔洞相连,锂离子并不能在材料中形成持续性的扩散。因此,其在结构上不符合我们设计准则的第三条。AIMD的模拟也显示,材料中的锂离子的扩散非常缓慢(图10)。根据第三条的设计原理,我们必须移除足够量的锌离子,在材料中制造更多的0-TM通道,并形成渗透的扩散网络。如图9所示,当移除与Oct(2LiS4-ZnS4)八面体孔洞配位的锌离子后,锂离子在其中的扩散势垒降低到0.28 eV。而对于Tet(2LiS4-ZnS4-PS4),即使将锌离子移除后,锂离子在其中的扩散势垒仍然高达0.48 eV,锂离子无法快速通过。如上分析,当LZGS中的锌离子被移除后,锂离子可以沿着Oct(2LiS4-PS4)-Oct(2LiS4-ZnVacS4)-Oct(2LiS4-PS4)的路径形成持续性的扩散。上述对于LZGS的分析同样适用于LZSS。在LZSS中,与锂离子迁移相关的局域结构包括Oct(2LiS4-SiS4)、Oct(2LiS4-ZnS4)和Tet(2LiS4-PS4-ZnS4)。NEB计算显示锂离子在不与过渡金属配位的Oct(2LiS4-SiS4)内的迁移势垒为0.15 eV(图11),其从结构上和性能上都符合0-TM通道的特征。同样地,Oct(2LiS4-SiS4)与两个高势垒的Oct(2LiS4-ZnS4)单元以及两个高势垒的Tet(2LiS4-PS4-ZnS4)单元相连,因此锂离子不能形成连续性扩散。根据设计准则的第三条,将Oct(2LiS4-ZnS4)中的锌离子移除后,锂离子在其内的迁移势垒大幅降低(图11),且由于其连接了两个Oct(2LiS4-SiS4)单元,因此锂离子可以形成连续性迁移。

图10

图10   AIMD模拟得到的在600 KLZGS和其富锂相Li2.5Zn0.5Ge0.5P0.5S4LZSS和其富锂相Li2.5Zn0.5Si0.5P0.5S4 中锂离子的均方位移 (MSD) (可以看到,在进行锌位的掺杂替代后,材料中锂离子的扩散能力大幅提高)

Fig. 10   The mean square displacement (MSD) of lithium ions in the of β-Li2ZnGeS4 (left) and Li2ZnSiS4 (right) calculated from the AIMD simulation at 600 K. After the composition engineering, the diffusion is significantly enhanced, proving the success of our designing strategy


图11

图11   利用NEB计算得到的LZSS中锂离子在Oct(2LiS4-SiS4) 孔洞内的迁移势垒曲线 (橙色实线);在Oct(2LiS4-ZnS4) 孔洞中的迁移势垒曲线(蓝色实线);在Tet(2LiS4-ZnS4-PS4) 中的迁移势垒曲线(绿色实线)。在LZGS富锂相中,部分锌离子被取代形成空位后,锂离子在Oct(2LiS4-VacS4) 孔洞中的迁移势垒曲线 (黑色实线)

Fig. 11   The calculated energy landscape of diffusion path in the LZSS from the Li site to the adjacent vacancy through the Oct(2LiS4-SiS4) cage (orange color). The Oct(2LiS4-ZnS4) cage (blue color). The Tet(2LiS4-ZnS4-PS4) cage (green color). The Oct(2LiS4-VacS4) cage (black color)


上述分析显示在LZGS和LZSS中,通过对过渡金属位的掺杂可以形成锂离子的持续性扩散通路。接下来,我们用AIMD进一步展示了这种掺杂的可行性。我们分别设计了LZGS的富锂相Li2.5Zn0.5Ge0.5/P0.5S4和LZSS的富锂相Li2.5Zn0.5Si0.5/P0.5S4,磷在锗位和硅位的掺杂以及锂在锌位上的掺杂是为了保持材料的电荷平衡。我们这里并没有计算验证0.5份的磷在锗位和硅位上的掺杂以及0.5份的锂在锌位上的掺杂是否在固溶体的互溶范围之内,但类似的掺杂方法已经在其他的材料组分上被实验证实是可行的[21- 22, 55-57]。AIMD模拟材料经过掺杂后,材料内锂离子的扩散能力有了显著的提升(图10)。同样地,阿伦尼乌斯拟合的结果显示,这两种材料经过掺杂后都成为快离子导体(如图6所示,Li2.5Zn0.5Ge0.5/P0.5S4的室温锂离子电导率为5.71 mS/cm,Li2.5Zn0.5Si0.5/P0.5S4的室温锂离子电导率为10.12 mS/cm)。这一结果证明我们的硫系快离子导体设计准则是可行的。

实际上,上文讨论的LZPS、LZGS和LZSS最初都不是作为快离子导体使用的。它们最初都是作为非线性光学(NLO)材料,被用来将固定频率的激光转换到其他频率。其他的NLO材料如LiZnSnS4[58]α‍ -Li2ZnGeS4[59]、Li2CdSiS4[60]、Li2CdSn/GeS4[61]、Li2CoSnS4、Li2MnGeS4[62]等都具有与本文所讨论的三种材料类似的晶体结构,它们都是类bcc结构的硫阴离子排布,按照我们的设计准则,它们都是快离子导体材料。下一步工作将对过渡金属位进行合适的掺杂改性,在材料中形成锂离子的渗透扩散网络。

3 结论

在本文中,对前人的bcc快离子导体设计准则进行了进一步的解析和深化,并提出了一个更明确、更具操作性的设计准则。对于硫系快离子导体来说,在结构上应满足硫阴离子框架具有类bcc的排布,在微观局域结构环境上应存在0-TM通道,且这些0-TM通道还必须互相连接,形成布满整个材料的全局渗透扩散网络。因此在进行快离子导体设计时,优选那些具有类bcc排布的硫阴离子框架材料,对其中的过渡金属进行掺杂和替代,在结构中形成尽可能多的0-TM通道,以满足渗透条件,将材料从贫离子导体转变成快离子导体。运用这一设计准则,本文作者寻找和设计了几种新型的硫系锂离子导体,如 ε -Li3PS4、Li2.5Zn0.5Ge0.5/P0.5S4以及Li2.5Zn0.5Si0.5/P0.5S4。AIMD模拟计算的结果显示,这三种材料具有较高的室温锂离子电导率,满足锂离子导体的性能要求。

参考文献

SUN Y K. Direction for development of next-generation lithium-ion batteries[J]. ACS Energy Letters, 2017, 2(12): 2694-2695.

[本文引用: 1]

GOODENOUGH J B, KIM Y. Challenges for rechargeable Li batteries[J]. Chemistry of Materials, 2010, 22(3): 587-603.

[本文引用: 1]

LI M, WANG C S, CHEN Z W, et al. New concepts in electrolytes[J]. Chemical Reviews, 2020, 120(14): 6783-6819.

[本文引用: 1]

LIANG J N, LUO J, SUN Q, et al. Recent progress on solid-state hybrid electrolytes for solid-state lithium batteries[J]. Energy Storage Materials, 2019, 21: 308-334.

[本文引用: 1]

HE Y M, LU C Y, LIU S, et al. Interfacial incompatibility and internal stresses in all-solid-state lithium ion batteries[J]. Advanced Energy Materials, 2019, 9(36): doi: 10.1002/aenm.201901810.

[本文引用: 1]

LV F, WANG Z Y, SHI L Y, et al. Challenges and development of composite solid-state electrolytes for high-performance lithium ion batteries[J]. Journal of Power Sources, 2019, 441: 227175.

YU C, GANAPATHY S, VAN ECK E R H, et al. Accessing the bottleneck in all-solid state batteries, lithium-ion transport over the solid-electrolyte-electrode interface[J]. Nature Communications, 2017, 8(1): 1086.

WAN J, XIE J, MACKANIC D G, et al. Status, promises, and challenges of nanocomposite solid-state electrolytes for safe and high performance lithium batteries[J]. Materials Today Nano, 2018, 4: 1-16.

[本文引用: 1]

LIU L L, XU J R, WANG S, et al. Practical evaluation of energy densities for sulfide solid-state batteries[J]. eTransportation, 2019, 1(C): 100010.

[本文引用: 1]

HAYASHI A, SAKUDA A, TATSUMISAGO M. Development of sulfide solid electrolytes and interface formation processes for bulk-type all-solid-state Li and Na batteries[J]. Frontiers in Energy Research, 2016, 4: 25.

LIU D, ZHU W, FENG Z, et al. Recent progress in sulfide-based solid electrolytes for Li-ion batteries[J]. Materials Science and Engineering: B, 2016, 213: 169-176.

[本文引用: 1]

YU C A, GANAPATHY S, DE KLERK N J J, et al. Unravelling Li-ion transport from picoseconds to seconds: Bulk versus interfaces in an argyrodite Li6PS5Cl-Li2S all-solid-state Li-ion battery[J]. Journal of the American Chemical Society, 2016, 138(35): 11192-11201.

[本文引用: 1]

KAMAYA N, HOMMA K, YAMAKAWA Y, et al. A lithium superionic conductor[J]. Nature Materials, 2011, 10(9): 682-686.

[本文引用: 2]

KUHN A, GERBIG O, ZHU C B, et al. A new ultrafast superionic Li-conductor: Ion dynamics in Li11Si2PS12 and comparison with other tetragonal LGPS-type electrolytes[J]. Physical Chemistry Chemical Physics, 2014, 16(28): 14669-14674.

[本文引用: 1]

BRON P, JOHANSSON S, ZICK K, et al. Li10SnP2S12: An affordable lithium superionic conductor[J]. Journal of the American Chemical Society, 2013, 135(42): 15694-15697.

[本文引用: 2]

SEINO Y, OTA T, TAKADA K, et al. A sulphide lithium super ion conductor is superior to liquid ion conductors for use in rechargeable batteries[J]. Energy & Environmental Science, 2014, 7(2): 627-631.

[本文引用: 2]

WANG Y, RICHARDS W D, ONG S P, et al. Design principles for solid-state lithium superionic conductors[J]. Nature Materials, 2015, 14(10): 1026-1031.

[本文引用: 2]

Clement R J, Lun Z, Ceder G. Cation-disordered rocksalt transition metal oxides and oxyfluorides for high energy lithium-ion cathodes[J]. Energy & environmental science. 2020, 13(2): 345-373.

HE X F, BAI Q, LIU Y S, et al. Crystal structural framework of lithium super-ionic conductors[J]. Advanced Energy Materials, 2019, 9(43): doi: 10.1002/aenm.201902078.

[本文引用: 1]

RICHARDS W D, WANG Y, MIARA L J, et al. Design of Li1+2 xZn1- xPS4, a new lithium ion conductor[J]. Energy & Environmental Science, 2016, 9(10): 3272-3278.

[本文引用: 3]

SUZUKI N, RICHARDS W D, WANG Y, et al. Synthesis and electrochemical properties of I-4-type Li1+2 xZn1- xPS4 solid electrolyte[J]. Chemistry of Materials, 2018, 30(7): 2236-2244.

[本文引用: 3]

KAUP K, LALÈRE F, HUQ A, et al. Correlation of structure and fast ion conductivity in the solid solution series Li1+2 xZn1- xPS4[J]. Chemistry of Materials, 2018, 30(3): 592-596.

[本文引用: 2]

SHYUE P O, WILLIAM D R, ANUBHAV J, et al. Python materials genomics (pymatgen): A robust, open-source python library for materials analysis[J]. Computational Materials Science, 2013, 68: 314-319.

[本文引用: 1]

KRESSE G, HAFNER J. Ab initio molecular dynamics for liquid metals[J]. Physical Review B, 1993, 47(1): 558-561.

KRESSE G, FURTHMÜLLER J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set[J]. Physical Review B, 1996, 54(16): 11169-11186.

[本文引用: 1]

BLÖCHL P E. Projector augmented-wave method[J]. Physical Review B, Condensed Matter, 1994, 50(24): 17953-17979.

[本文引用: 1]

KRESSE G, JOUBERT D. From ultrasoft pseudopotentials to the projector augmented-wave method[J]. Physical Review B, 1999, 59(3): 1758-1775.

[本文引用: 1]

KRESSE G, HAFNER J. Ab initio molecular dynamics for open-shell transition metals[J]. Physical Review B, 1993, 48(17): 13115-13118.

[本文引用: 1]

PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple[J]. Physical Review Letters, 1996, 77(18): 3865-3868.

[本文引用: 1]

HENKELMAN G, UBERUAGA B P, JóNSSON H. A climbing image nudged elastic band method for finding saddle points and minimum energy paths[J]. The Journal of Chemical Physics, 2000, 113(22): 9901.

[本文引用: 1]

HENKELMAN G, JÓNSSON H. Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points[J]. The Journal of Chemical Physics, 2000, 113(22): 9978-9985.

[本文引用: 1]

TOGO A, TANAKA I. First principles phonon calculations in materials science[J]. Scripta Materialia, 2015, 108: 1-5.

[本文引用: 1]

TOGO A, CHAPUT L, TANAKA I, et al. First-principles phonon calculations of thermal expansion in Ti3SiC2, Ti3AlC2, and Ti3GeC2[J]. Physical Review B, 2010, 81(17): 174301.

[本文引用: 1]

VERLET L. Computer "experiments" on classical fluids. I. thermodynamical properties of lennard-jones molecules[J]. Physical Review, 1967, 159(1): 98-103.

[本文引用: 1]

NOSÉ S. A unified formulation of the constant temperature molecular dynamics methods[J]. The Journal of Chemical Physics, 1984, 81(1): 511-519.

[本文引用: 1]

HOOVER W G. Canonical dynamics: Equilibrium phase-space distributions[J]. Physical Review A, 1985, 31(3): 1695-1697.

[本文引用: 1]

STAUFFER D, AHARONY A. Introduction to percolation theory[M]. London: Taylor & Francis, 2018.

[本文引用: 1]

EWING R, HUNT A. Percolation theory for flow in porous media[M]. Berlin, Heidelberg: Springer Berlin Heidelberg, 2009.

[本文引用: 1]

NEWMAN M E, ZIFF R M. Efficient Monte Carlo algorithm and high-precision results for percolation[J]. Physical Review Letters, 2000, 85(19): 4104-4107.

[本文引用: 1]

CEDER G, ONG S P, WANG Y. Predictive modeling and design rules for solid electrolytes[J]. MRS Bulletin, 2018, 43(10): 746-751.

[本文引用: 1]

MEHRER, H. Diffusion in solids: Fundamentals, methods, materials, diffusion-controlled processes[J]. Springer series in solid state science, 2007, 155.

[本文引用: 1]

HE X F, ZHU Y Z, MO Y F. Origin of fast ion diffusion in super-ionic conductors[J]. Nature Communications, 2017, 8: 15893.

[本文引用: 1]

URBAN A, LEE J, CEDER G. The configurational space of rocksalt-type oxides for high-capacity lithium battery electrodes[J]. Advanced Energy Materials, 2014, 4(13): 1400478.

[本文引用: 2]

LEE J, URBAN A, LI X, et al. Unlocking the potential of cation-disordered oxides for rechargeable lithium batteries[J]. Science, 2014, 343(6170): 519-522.

JI H W, URBAN A, KITCHAEV D A, et al. Hidden structural and chemical order controls lithium transport in cation-disordered oxides for rechargeable batteries[J]. Nature Communications, 2019, 10: 592.

URBAN A, MATTS I, ABDELLAHI A, et al. Computational design and preparation of cation-disordered oxides for high-energy-density Li-ion batteries[J]. Advanced Energy Materials, 2016, 6(15): 1600488.

[本文引用: 1]

HUANG Y, LIU L, ZHU Y Y, et al. A new model on cation distribution in cation-disordered Li1+ xTM1-xO2 cathodes[J]. Solid State Ionics, 2020, 351: 115341.

[本文引用: 1]

HOMMA K, YONEMURA M, KOBAYASHI T, et al. Crystal structure and phase transitions of the lithium ionic conductor Li3PS4[J]. Solid State Ionics, 2010, 182(1): 53-58.

[本文引用: 2]

IIKUBO S, SHIMOYAMA K, KAWANO S, et al. Novel stable structure of Li3PS4 predicted by evolutionary algorithm under high-pressure[J]. AIP Advances, 2018, 8(1): 15008.

[本文引用: 1]

KIM J S, JUNG W D, CHOI S J, et al. Thermally induced S‑sublattice transition of Li3PS4 for fast lithium-ion conduction[J]. J. Phys. Chem. Lett. 2018, 9, 5592-5597.

[本文引用: 1]

YANG J J, TSE J S. First-principles molecular simulations of Li diffusion in solid electrolytes Li3PS4[J]. Computational Materials Science, 2015, 107: 134-138.

LIM M S, JHI S H. First-principles study of lithium-ion diffusion in β-Li3PS4 for solid-state electrolytes[J]. Current Applied Physics, 2018, 18(5): 541-545.

[本文引用: 1]

HUANG Y, WU K, CHENG J N, et al. Li2ZnGeS4: A promising diamond-like infrared nonlinear optical material with high laser damage threshold and outstanding second-harmonic generation response[J]. Dalton Transactions, 2019, 48(14): 4484-4488.

[本文引用: 1]

LI G M, CHU Y, ZHOU Z X. From AgGaS2 to Li2ZnSiS4: Realizing impressive high laser damage threshold together with large second-harmonic generation response[J]. Chemistry of Materials, 2018, 30(3): 602-606.

[本文引用: 1]

ZHOU L D, ASSOUD A, SHYAMSUNDER A, et al. An entropically stabilized fast-ion conductor: Li3.25[Si0.25P0.75]S4[J]. Chemistry of Materials, 2019, 31(19): 7801-7811.

[本文引用: 1]

HARM S, HATZ A K, MOUDRAKOVSKI I, et al. Lesson learned from NMR: Characterization and ionic conductivity of LGPS-like Li7SiPS8[J]. Chemistry of Materials, 2019, 31(4): 1280-1288.

KANNO R, MURAYAMA M. Lithium ionic conductor thio-LISICON the Li2S-GeS2-P2S5 system[J]. Journal of the Electrochemical Society, 2001, 148(7): A742-A746.

[本文引用: 1]

LEKSE J W, LEVERETT B M, LAKE C H, et al. Synthesis, physicochemical characterization and crystallographic twinning of Li2ZnSnS4[J]. Journal of Solid State Chemistry, 2008, 181(12): 3217-3222.

[本文引用: 1]

ZHANG J H, CLARK D J, BRANT J A, et al. α -Li2ZnGeS4: A wide-bandgap diamond-like semiconductor with excellent balance between laser-induced damage threshold and second harmonic generation response[J]. Chemistry of Materials, 2020, 32(20): 8947-8955.

[本文引用: 1]

LI G M, CHU Y, LI J, et al. Li2CdSiS4, a promising IR NLO material with a balanced Eg and SHG response originating from the effect of Cd with d10 configuration[J]. Dalton Transactions, 2020, 49(6): 1975-1980.

[本文引用: 1]

LI Y L, FAN W L, SUN H G, et al. Electronic, optical and lattice dynamic properties of the novel diamond-like semiconductors Li2CdGeS4 and Li2CdSnS4[J]. Journal of Physics: Condensed Matter, 2011, 23(22): 225401.

[本文引用: 1]

BRANT J A, CLARK D J, KIM Y S, et al. Outstanding laser damage threshold in Li2MnGeS4 and tunable optical nonlinearity in diamond-like semiconductors[J]. Inorganic Chemistry, 2015, 54(6): 2809-2819.

[本文引用: 1]

/