凌仕剛,高　健,褚　賡,黃　杰,肖睿娟,歐陽楚英,李　泓,陳立泉

(1.中國科學(xué)院物理研究所，北京 100190)(2.江西師范大學(xué)，江西 南昌 330027)

?

特約專欄

第一作者：凌仕剛，男，1988年生，博士研究生

高通量計(jì)算在鋰電池材料篩選中的應(yīng)用

凌仕剛1,高健1,褚賡1,黃杰1,肖睿娟1,歐陽楚英2,李泓1,陳立泉1

(1.中國科學(xué)院物理研究所，北京 100190)(2.江西師范大學(xué)，江西 南昌 330027)

摘要：先進(jìn)電池技術(shù)是未來十年世界各國前沿技術(shù)競爭的制高點(diǎn)，鋰電池由于具有電壓高、比能量高、充放電壽命長、工作溫度范圍寬等優(yōu)點(diǎn)，已成為許多移動(dòng)電子產(chǎn)品、電動(dòng)汽車以及風(fēng)電和光伏電儲(chǔ)能的首選。因此,更快的開發(fā)出能量密度高、安全性好、充放電速度快的高性能鋰電池材料十分必要。借助高通量計(jì)算方法，通過設(shè)置合理的篩選條件，可加快鋰電池材料的研發(fā)。通過使用自主編寫的計(jì)算軟件及高通量自動(dòng)化計(jì)算流程，以無機(jī)材料晶體結(jié)構(gòu)數(shù)據(jù)庫中現(xiàn)有的結(jié)構(gòu)為基礎(chǔ)，針對(duì)電極材料和固體電解質(zhì)材料分別應(yīng)具有的特點(diǎn)，實(shí)現(xiàn)了對(duì)鋰電池材料的初步篩選。計(jì)算結(jié)果與現(xiàn)有的實(shí)驗(yàn)數(shù)據(jù)匯集到一起，形成了龐大的數(shù)據(jù)庫。進(jìn)一步借助統(tǒng)計(jì)方法進(jìn)行數(shù)據(jù)挖掘，有望揭示材料的結(jié)構(gòu)-性能關(guān)系及內(nèi)在物理規(guī)律，縮短材料從研發(fā)到應(yīng)用的全過程，為開發(fā)全新電池材料乃至電池體系提供有力支持。

關(guān)鍵詞：鋰電池材料；高通量計(jì)算；熱力學(xué)計(jì)算；鍵價(jià)和方法

1前言

能源的清潔高效利用離不開高性能儲(chǔ)能器件。從鉛酸電池、鎘鎳電池到鎳氫電池，再到鋰離子電池，化學(xué)電源技術(shù)得到了迅猛發(fā)展，能量密度也大幅提高。與其他蓄電池比較，鋰電池具有電壓高、比能量高、充放電壽命長、無記憶效應(yīng)、無污染、快速充電、自放電率低、工作溫度范圍寬和安全可靠等優(yōu)點(diǎn)[1]。鋰電池的平均比能量是鎘鎳電池的2.6倍、氫鎳電池的1.75倍，其工作電壓為3.6 V，是鎘鎳、氫鎳電池的3倍[2]。在民用領(lǐng)域，鋰電池已經(jīng)被廣泛應(yīng)用于各種電子產(chǎn)品、電動(dòng)工具[2]，并已成功應(yīng)用于電動(dòng)車、混合動(dòng)力車等[3]，同時(shí)鋰電池或鈉電池體系也是風(fēng)電和光伏電儲(chǔ)能的首選方案[4]。

當(dāng)前，清潔高效的能源儲(chǔ)存技術(shù)已成為工業(yè)可持續(xù)發(fā)展的核心關(guān)鍵技術(shù)之一，因此研發(fā)具有更高能量密度、更高功率密度、更長壽命、更加安全的可充放鋰電池儲(chǔ)能器件十分必要?，F(xiàn)有鋰電池正極材料體系主要包括層狀材料(如LiCoO2)、尖晶石結(jié)構(gòu)材料(如LiMn2O4)以及橄欖石結(jié)構(gòu)材料(如LiFePO4)[5]。這3種材料的嵌鋰電位平臺(tái)分別為3.9 V、4.1 V和3.6 V，其容量分別為130～140 mAh·g-1、120～130 mAh·g-1和150～160 mAh·g-1。在這些材料體系基礎(chǔ)上，還有其他通過改性的正極材料：比如層狀的Li(NixCoyMn1-x-y)O2材料，其比容量能達(dá)到200 mAh·g-1甚至更高[6]；通過對(duì)LiMn2O4材料摻雜，可以提高其電位到4.5 V[7]。然而，由于這些材料體系的結(jié)構(gòu)穩(wěn)定性、可轉(zhuǎn)移電子和離子數(shù)的限制，經(jīng)過近20年的研究，其能量密度已經(jīng)接近其理論極限。要進(jìn)一步提高電池的能量密度等性能，必須探索出全新的鋰電池正極材料體系。

就負(fù)極材料而言，目前的炭負(fù)極材料的應(yīng)用潛力也已經(jīng)接近其理論極限，且應(yīng)用過程中面臨著各種技術(shù)挑戰(zhàn)與難題。比如，石墨負(fù)極材料表面在使用過程中SEI(Solid Electrolyte Interphase)膜生長的控制、首次充放電效率的提高、對(duì)正極提供的Li源的消耗導(dǎo)致可逆容量損失等[8]。這些難題是炭負(fù)極材料本征特性帶來的，要徹底解決這些問題，必須使用全新的負(fù)極材料體系。硅負(fù)極材料雖然理論容量高達(dá)4 200 mAh·g-1，但其嵌Li過程的體積膨脹(～400%)導(dǎo)致其循環(huán)性能的控制非常難，循環(huán)過程中Si和Li的反復(fù)合金化/去合金化過程的機(jī)械分裂和粉末化，將導(dǎo)致材料結(jié)構(gòu)的坍塌,進(jìn)而使得材料的容量急劇下降[9]。而Si材料嵌Li過程的體積膨脹也是材料的本征屬性，難以從根源上徹底解決，因此迫切需要尋找性能更加優(yōu)越的負(fù)極材料。

除了提高電池的能量密度，鋰電池的安全性能也急待提高。電池的安全性能在很大程度上和使用的液態(tài)有機(jī)電解液體系有關(guān)。使用固體電解質(zhì)體系的全固態(tài)鋰離子電池，是從根源上解決電池的安全性問題的方案[10-11]。使用固體電解質(zhì)材料體系的一個(gè)最大問題是鋰離子在固體電解質(zhì)材料中的離子輸運(yùn)很慢。一般來說，鋰離子電導(dǎo)率在固體電解質(zhì)中要比在液體電解液中低1~2個(gè)數(shù)量級(jí)[12]。最近，日本一研究小組在NatureMaterials[13]上報(bào)道了他們開發(fā)的Li10GeP2S12材料, 該材料在室溫下鋰離子的電導(dǎo)率可以達(dá)到10-2S/cm，說明固體電解質(zhì)材料在鋰離子電導(dǎo)率方面已能滿足電池應(yīng)用的要求，但該材料在自身的穩(wěn)定性、和電極材料界面的兼容性等其它方面還存在一些問題[14]，開發(fā)兼具高離子電導(dǎo)率、高穩(wěn)定性、高機(jī)械強(qiáng)度的固體電解質(zhì)材料勢在必行。

由此可見，當(dāng)前鋰電池產(chǎn)業(yè)面臨的瓶頸問題是現(xiàn)有的電池材料體系的本征特性,決定了現(xiàn)有電池體系的潛力有限，不能滿足高能量密度、高功率密度、高安全性能、低成本等多方面的需求。而要從根源上解決這些問題，必須開發(fā)出全新的鋰電池材料，進(jìn)而開發(fā)出全新的鋰電池體系，如全固態(tài)鋰電池體系[11]，Li-O2電池體系[15]，Li-S電池體系[16]，Li-CO2電池體系[17]等。然而，傳統(tǒng)的材料研發(fā)方法周期很長，過去20年對(duì)鋰電池材料的開發(fā)沿用“炒菜”式的實(shí)驗(yàn)組合方法，僅開發(fā)了6種能夠?qū)嶋H應(yīng)用的電池材料，其它工作大多是對(duì)這幾種材料的改性或修飾。近年來，隨著研究人員對(duì)鋰電池中熱力學(xué)及動(dòng)力學(xué)過程理解的不斷深入，加之計(jì)算機(jī)運(yùn)算能力的飛速發(fā)展，理論計(jì)算已被有效的應(yīng)用到電池材料的研究中[18-20]。采用理論模擬的方法一方面可以對(duì)已有的材料進(jìn)行模擬和計(jì)算，分析它們的微觀物理機(jī)理，對(duì)相應(yīng)的物理現(xiàn)象加以解釋，為實(shí)驗(yàn)提供理論根據(jù)；另一方面能預(yù)測未知材料的性能，為實(shí)驗(yàn)指導(dǎo)方向。對(duì)鋰電池材料的計(jì)算具體包括：計(jì)算材料的電子結(jié)構(gòu)[21]，分析其電子導(dǎo)電性[18]；模擬離子在材料中擴(kuò)散，分析材料的離子導(dǎo)電性[22]；通過模擬充、放電過程，計(jì)算脫/嵌鋰電位，分析材料的體積和結(jié)構(gòu)變化，以此來說明它的結(jié)構(gòu)穩(wěn)定性和循環(huán)性[23]；研究電極材料表面以及表面和電解液之間相互作用，探索電池材料失效的機(jī)理等[24]。自然界中有著成千上萬的化合物，已經(jīng)研究過的只是少數(shù)。對(duì)已存在物質(zhì)的基本性質(zhì)進(jìn)行計(jì)算，有可能從中篩選出可以用于鋰電池的新材料。對(duì)每種物質(zhì)的計(jì)算無論是從方法選擇、計(jì)算參數(shù)、數(shù)據(jù)輸入等方面都具有一定的相似性，因此采用自動(dòng)化的操作來調(diào)用軟件、輸入數(shù)據(jù)、控制運(yùn)算步驟、分析計(jì)算結(jié)果就十分必要。這種高度自動(dòng)化的運(yùn)算過程稱為高通量計(jì)算。高通量計(jì)算具有高效能、可并行、可擴(kuò)展等優(yōu)點(diǎn)，有助于我們設(shè)計(jì)一些篩選材料體系的方案。借助高通量的計(jì)算，可以有目標(biāo)的去發(fā)現(xiàn)一些新材料體系，同時(shí)研究這些新體系，進(jìn)一步增加對(duì)鋰電池材料的認(rèn)識(shí)，從而形成不斷深入的研究和開發(fā)過程，最終有助于我們更快的開發(fā)出能量密度高、安全性好、充放電速度快的高性能鋰電池材料。

2高通量計(jì)算發(fā)展概述

高通量方法在生物學(xué)、藥學(xué)和化學(xué)領(lǐng)域已得到廣泛的應(yīng)用，例如高通量測序(High-Throughput Sequencing)方法可以并行的合成多達(dá)50萬條DNA分子并進(jìn)行序列測定[25-26]；高通量篩選(High-Throughput Screening)技術(shù)則通過分子水平和細(xì)胞水平的實(shí)驗(yàn)，以自動(dòng)化操作同時(shí)制備、檢測成千上萬的樣品，實(shí)現(xiàn)藥物分子的篩選[27]。在材料科學(xué)中，1995年Xiang等人[28]將薄膜沉積和物理掩膜技術(shù)相結(jié)合，在200 μm×200 μm大小的樣品上制備出一系列不同組分的BiSrCaCuO和YBaCuO超導(dǎo)薄膜，實(shí)現(xiàn)了對(duì)材料的組合制備和高通量篩選，為新材料的發(fā)現(xiàn)提供了一種新的思路和方法。

在鋰電池材料設(shè)計(jì)開發(fā)的過程中，采用高通量計(jì)算的理論方法始于MIT的Ceder研究組，他們從2010年開始開展了稱之為“鋰離子電池材料基因工程”的研究[33]。該方法通過對(duì)含聚陰離子XO4(X=P，S，As，Si)的化合物中元素的替代,來產(chǎn)生新化合物，計(jì)算了其能量密度、電壓、脫Li后的體積變化等參數(shù)，據(jù)此篩選新材料[34]。他們從自然界存在的礦物Sidorenkite結(jié)構(gòu)出發(fā)，對(duì)其進(jìn)行元素替換，構(gòu)造出270多種組分和結(jié)構(gòu)，并計(jì)算其性能，從中篩選出了幾種材料：Li3Mn(CO3)(PO4)、Li2V(CO3)(PO4)、Li3V(CO3)(SiO4)等。應(yīng)用類似辦法，他們還對(duì)含有(PmOn)陰離子團(tuán)的各種過渡金屬與Li的化合物進(jìn)行了篩選[35]。同時(shí)，在上述研究過程中，他們也初步發(fā)展了一套系統(tǒng)的性能計(jì)算、數(shù)據(jù)生成、數(shù)據(jù)分析的辦法[33]。日本的Tanaka研究組將第一性原理分子動(dòng)力學(xué)模擬與機(jī)器學(xué)習(xí)技術(shù)相結(jié)合，研究了LISICON類型材料的離子輸運(yùn)，通過將理論計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)庫結(jié)合,可有效預(yù)測LiO1/2-AOm/2-BOn/2體系不同組分的離子電導(dǎo)率，有助于加快固態(tài)電解質(zhì)的開發(fā)[36]。

上述研究工作中所用到的高通量計(jì)算方法都具有圖1[33]所示的共同特點(diǎn)。首先從外部的結(jié)構(gòu)數(shù)據(jù)庫中選擇數(shù)據(jù)，產(chǎn)生可被計(jì)算軟件調(diào)用的輸入文件,并進(jìn)行計(jì)算得到材料相應(yīng)的性質(zhì)數(shù)據(jù)，運(yùn)算結(jié)果保存到數(shù)據(jù)庫以備進(jìn)一步的分析，獲得的新知識(shí)可擴(kuò)充原先的數(shù)據(jù)庫并有助于更為準(zhǔn)確的數(shù)據(jù)選擇。由此可見，形成并完善一套程序化的高通量計(jì)算流程，把各種計(jì)算軟件包、編寫的單一功能的計(jì)算程序或指令(如輸入文件生成、材料性質(zhì)模擬、結(jié)果分析程序等)和計(jì)算硬件設(shè)備等關(guān)聯(lián)起來，使得整個(gè)過程能夠自動(dòng)完成，是實(shí)現(xiàn)高通量計(jì)算提高整體效率的關(guān)鍵。上面提到的由Curtarolo等人設(shè)計(jì)的AFLOW[30]和由Ceder小組設(shè)計(jì)的計(jì)算流程[33]均實(shí)現(xiàn)了晶體結(jié)構(gòu)數(shù)據(jù)庫與第一性原理計(jì)算程序VASP[37]的自動(dòng)調(diào)用功能，從而通過高通量計(jì)算獲得材料的形成能等熱力學(xué)數(shù)據(jù)以及電子結(jié)構(gòu)的信息。但由于篩選不同類型的材料時(shí)所關(guān)注的性質(zhì)往往是不同的，希望調(diào)用的材料數(shù)據(jù)庫、模擬方法、分析方法都可能有所變化，所以針對(duì)各種不同的材料還需設(shè)計(jì)相應(yīng)軟件以實(shí)現(xiàn)特定的自動(dòng)化運(yùn)算流程。

圖1　高通量計(jì)算篩選材料的數(shù)據(jù)流程圖[33]Fig.1　Data flow diagram of high-throughput computational screening[33]

目前對(duì)于鋰電池材料的高通量篩選，國內(nèi)外的研究工作都處于起步階段，尚存在以下問題：①偏重于材料熱力學(xué)性質(zhì)，如形成能、熱穩(wěn)定性等的計(jì)算，缺乏對(duì)鋰離子輸運(yùn)這一重要的動(dòng)力學(xué)性能的考慮；②篩選所用的約束條件缺乏針對(duì)性和層次性，造成大量重復(fù)及不必要的計(jì)算；③計(jì)算采用的多是基于密度泛函的第一性原理方法，尚未結(jié)合其他更多的計(jì)算方法或軟件；④對(duì)計(jì)算得到的大量數(shù)據(jù)尚未進(jìn)行有效的挖掘。事實(shí)上，材料的性能完全由材料的結(jié)構(gòu)和組分兩個(gè)因素決定，通過大規(guī)模的計(jì)算，挖掘出材料結(jié)構(gòu)、組分和材料性能之間的聯(lián)系，構(gòu)建出材料的“構(gòu)效關(guān)系”，也是高通量計(jì)算的重要目標(biāo)之一。通過對(duì)篩選過程設(shè)立有針對(duì)性的外部約束條件，包括應(yīng)用要求約束、基本物理規(guī)律約束、溫度壓強(qiáng)等外部條件，構(gòu)建材料組分和結(jié)構(gòu)模型。通過高通量計(jì)算技術(shù)來計(jì)算這些結(jié)構(gòu)模型對(duì)應(yīng)的化合物的各種性能，并反饋到結(jié)構(gòu)性能數(shù)據(jù)庫，用于研究鋰電池材料中的物理問題，同時(shí)篩選出滿足應(yīng)用約束條件的材料。并通過各種實(shí)驗(yàn)方法制備這些材料，然后測試和表征這些材料的物理化學(xué)性能，最終達(dá)到理解鋰電池材料物理問題，并發(fā)現(xiàn)新材料的目的。

可以預(yù)見在未來幾年內(nèi)，應(yīng)用高通量計(jì)算方法，結(jié)合實(shí)驗(yàn)研究，將大大加速我們對(duì)鋰電池材料的研發(fā)，加深對(duì)鋰電池中各種問題的理解，揭示材料的結(jié)構(gòu)-性能關(guān)系及內(nèi)在物理規(guī)律，從而縮短材料從研發(fā)到應(yīng)用的全過程。

3鋰電池材料性能的計(jì)算方法

在高通量計(jì)算中，可以根據(jù)所篩選材料的特點(diǎn)綜合使用計(jì)算材料學(xué)中的各種技術(shù)方法和手段，包括熱力學(xué)計(jì)算、第一性原理計(jì)算、分子動(dòng)力學(xué)模擬、蒙特卡羅模擬等。對(duì)于鋰電池而言，電極材料的儲(chǔ)鋰容量、嵌鋰電位、電子導(dǎo)電性、鋰離子擴(kuò)散率、脫嵌鋰過程中的體積變化、結(jié)構(gòu)相變，電解質(zhì)材料的離子電導(dǎo)率、電化學(xué)窗口寬度、結(jié)構(gòu)穩(wěn)定性，電解質(zhì)與電極之間界面的化學(xué)穩(wěn)定性等因素都會(huì)影響到電池的性能。上述的材料性質(zhì)都有相應(yīng)的計(jì)算方法，根據(jù)各種方法的原理，會(huì)得到不同計(jì)算精度的結(jié)果，各種方法適當(dāng)組合，可以構(gòu)建出合理的材料篩選流程。

3.1理論儲(chǔ)鋰容量

電極材料的儲(chǔ)鋰容量與可轉(zhuǎn)移電子數(shù)、可轉(zhuǎn)移離子數(shù)有關(guān)。對(duì)于給定的電極材料，其理論儲(chǔ)鋰容量可由式(1)計(jì)算得到[38]：

Capacity=nF/3.6M(mAh·g-1)

(1)

其中，n為每摩爾所能轉(zhuǎn)移的電子數(shù)目，M為摩爾質(zhì)量(g·mol-1)，F(xiàn)為法拉第常數(shù)。因此，根據(jù)電極反應(yīng)中可以參與反應(yīng)的Li+及電子的數(shù)目，即可得到電極材料的理論儲(chǔ)鋰容量，可用作電極材料初步篩選的參數(shù)之一。

3.2嵌鋰電位

電極材料在脫嵌鋰過程中的平均電位與電極反應(yīng)的吉布斯自由能ΔrGΘ有關(guān)，如式(2)所示：

ΔrGΘ=-nEF

(2)

其中，E即為電極反應(yīng)的熱力學(xué)平衡電壓。一方面，吉布斯自由能的數(shù)據(jù)可以從熱力學(xué)實(shí)驗(yàn)數(shù)據(jù)手冊[39]中獲得，對(duì)于缺少實(shí)驗(yàn)數(shù)據(jù)的反應(yīng)，可以通過第一性原理計(jì)算得到[40]。

3.3離子輸運(yùn)

鋰電池中核心的動(dòng)力學(xué)過程是鋰離子的輸運(yùn)，離子傳輸?shù)目炻吐窂街苯佑绊戨姵爻浞烹姷男阅?。在固體和液體中，離子的傳輸是通過離子向近鄰位置隨機(jī)跳躍的擴(kuò)散過程實(shí)現(xiàn)的。離子在運(yùn)動(dòng)過程中會(huì)受到由周圍離子形成的勢場的阻礙，這種阻礙的大小可以用離子擴(kuò)散活化能表示。對(duì)材料中離子輸運(yùn)的計(jì)算方法可分為基于幾何結(jié)構(gòu)[41-43]、基于成鍵配位[44]和基于能量變化[19,22]3類。基于幾何結(jié)構(gòu)的方法是分析晶格中存在的可供Li+運(yùn)動(dòng)的幾何空間，能定性理解Li+在晶格中運(yùn)動(dòng)的難易程度?；诔涉I配位的鍵價(jià)理論[44]考慮了與Li+成鍵的陰離子對(duì)Li+運(yùn)動(dòng)的制約，通過引入基于Morse勢的勢能函數(shù)，并更進(jìn)一步考慮Li+與其他陽離子的庫侖排斥勢，可以得出Li+在晶格中形成連通路徑時(shí)所需要越過的勢壘大小[45]，由于該方法中使用的Morse勢為經(jīng)驗(yàn)勢，因此勢壘的計(jì)算屬于半定量性質(zhì)。更為精確的基于能量的計(jì)算，可借助過渡態(tài)理論，通過DFT計(jì)算來得到，例如采用爬坡彈性能帶(Climbing Image Nudged Elastic Band)方法優(yōu)化路徑和計(jì)算鋰離子遷移勢壘[46]。需要指出的是，對(duì)于上述幾種方法，隨著計(jì)算精度的提高，所需要的運(yùn)算時(shí)間會(huì)大幅增長，如對(duì)于單一結(jié)構(gòu)的計(jì)算，半定量的鍵價(jià)和方法的計(jì)算量在分鐘量級(jí)，而基于DFT的第一性原理計(jì)算則往往耗時(shí)數(shù)小時(shí)甚至數(shù)天。因此，將各種不同精度的方法合理組合，用于材料篩選的不同階段，是有助于提高篩選效率的辦法。

3.4電子結(jié)構(gòu)

電池材料的能帶結(jié)構(gòu)也可以為材料篩選提供十分有價(jià)值的信息。電極材料中，過渡金屬與氧的分波態(tài)密度的相對(duì)位置，決定了材料中電荷補(bǔ)償?shù)膩碓春脱醯姆€(wěn)定性[1]，從能態(tài)密度中還能定性預(yù)測電極與電解質(zhì)界面的化學(xué)穩(wěn)定性[47]，此外電池的倍率性能也與電極的電子導(dǎo)電性有關(guān)[18]。電池中的固態(tài)電解質(zhì)材料需要具有電子絕緣的特性，其電化學(xué)窗口的寬度與材料的能隙寬度有關(guān)[11]。因此電子結(jié)構(gòu)的計(jì)算無論對(duì)電極材料還是電解質(zhì)材料的篩選都十分必要。

電子結(jié)構(gòu)的計(jì)算方法也可以分為半經(jīng)驗(yàn)和第一性原理計(jì)算兩類[48]。半經(jīng)驗(yàn)方法常需要用到一些擬和參數(shù)，而基于密度泛函理論的第一性原理計(jì)算僅需要用到基本物理常數(shù)，在知道所需模擬體系的結(jié)構(gòu)和組成的情況下，不需要其他經(jīng)驗(yàn)的或擬和的可調(diào)參數(shù)，因此無論在計(jì)算過程還是在結(jié)果的精確度上都具有半經(jīng)驗(yàn)方法不可比擬的優(yōu)勢。第一性原理從量子力學(xué)的基本方程出發(fā)，求解描述體系的波函數(shù)及相應(yīng)的本征能量，并以此推導(dǎo)體系的各種性質(zhì)，電子態(tài)的求解過程是一個(gè)自洽計(jì)算的過程，因此比半經(jīng)驗(yàn)方法需要更多的計(jì)算時(shí)間[49]。因此有必要在材料篩選的不同階段合理選擇計(jì)算方法。

根據(jù)上述鋰電池材料中各種性質(zhì)的計(jì)算方法，我們針對(duì)電極材料和電解質(zhì)材料分別需要滿足的特性，構(gòu)筑出篩選規(guī)則，對(duì)無機(jī)材料晶體結(jié)構(gòu)數(shù)據(jù)庫[50]里已有的結(jié)構(gòu)進(jìn)行高通量計(jì)算,并研究其性質(zhì)與結(jié)構(gòu)之間的關(guān)系。

4鋰電池材料篩選的基本思路

在鋰電池中，電極材料與電解質(zhì)材料在性能上需要滿足不同的要求。對(duì)于電極材料的研發(fā)，側(cè)重于提高電極材料的比容量和提高電池的工作電壓。對(duì)于固體電解質(zhì)材料的研究，則希望找到離子電導(dǎo)率高，電化學(xué)窗口寬，并能與正負(fù)極穩(wěn)定匹配的材料。此外所有的電池材料都需要具有較低的成本、無毒無放射性，減少對(duì)環(huán)境的污染，以期持續(xù)利用和發(fā)展。

在上述約束條件和篩選目標(biāo)的指導(dǎo)下，我們首先考慮從無機(jī)材料晶體結(jié)構(gòu)數(shù)據(jù)庫(ICSD)[50]中，根據(jù)化合物的組成元素進(jìn)行初步篩選，如圖2所示，再進(jìn)一步分別考慮它們作為電極或電解質(zhì)材料的可能性。具體方案可以用圖3表示。對(duì)于從ICSD中初步篩選出的化合物及其衍生物，首先考慮材料的儲(chǔ)鋰極限容量(優(yōu)先考慮Li含量，其次考慮可轉(zhuǎn)移電子數(shù))。將計(jì)算出的理論容量按照140 mAh·g-1(LiCoO2材料的容量)為參考，如果大于這個(gè)參考值，則首先考慮其作為電極材料應(yīng)用。接下來計(jì)算其脫/嵌Li電位，若電位合適，則繼續(xù)計(jì)算脫Li后體積變化以及結(jié)構(gòu)穩(wěn)定性。如果脫Li后結(jié)構(gòu)不穩(wěn)定，我們則轉(zhuǎn)而考慮其作為固體電解質(zhì)。當(dāng)電位、結(jié)構(gòu)穩(wěn)定性兩個(gè)都符合要求，則繼續(xù)計(jì)算離子擴(kuò)散勢壘等參數(shù)，并判斷離子輸運(yùn)是否符合設(shè)置的約束條件。如果符合離子遷移勢壘要求，則計(jì)算電子結(jié)構(gòu)與晶格動(dòng)力學(xué)參數(shù)，對(duì)前面的結(jié)果進(jìn)行進(jìn)一步修正，得到更為精確的數(shù)據(jù)。

圖2　從ICSD數(shù)據(jù)庫中初步篩選鋰電池材料的方案Fig.2　Scheme of preliminary screening of lithium battery materials from ICSD database

圖3　高通量計(jì)算篩選鋰電池電極材料和固體電解質(zhì)材料的流程Fig.3　Procedures of screening electrode materials and solid electrolytes using high-throughput calculation method

當(dāng)開始計(jì)算的理論容量小于140 mAh·g-1時(shí)，我們直接考慮其為固體電解質(zhì)材料的應(yīng)用。我們首先判斷其中是否含有可變價(jià)元素，如果有，則不再做其它考慮。如果不含可變價(jià)元素，我們則計(jì)算其電子結(jié)構(gòu)，判斷其電子絕緣特性，如果是電子導(dǎo)體或帶隙很窄的半導(dǎo)體，那么這種材料也不能用做固體電解質(zhì)，可以不再考慮。如果符合電子絕緣的條件，則考慮其分解電壓。滿足條件后，則進(jìn)一步考慮Li離子在材料中的遷移勢壘。如果滿足遷移勢壘條件，則進(jìn)行晶格動(dòng)力學(xué)參數(shù)修正。

在計(jì)算方法上，根據(jù)化合物中Li的含量以及可變價(jià)元素可能轉(zhuǎn)移的電子數(shù)目來計(jì)算材料的理論儲(chǔ)鋰容量[38]，采用熱力學(xué)基本數(shù)據(jù)來計(jì)算電極反應(yīng)的能量密度，并獲得材料的平均嵌鋰電位[38]，對(duì)于缺少熱力學(xué)數(shù)據(jù)的電極反應(yīng)，采用DFT方法計(jì)算嵌鋰電位[40]，同時(shí)得到脫/嵌鋰前后的體積變化。對(duì)于材料是否電子絕緣，以及對(duì)應(yīng)的分解電位，也采用DFT計(jì)算來進(jìn)行預(yù)測。對(duì)于鋰離子傳導(dǎo)的動(dòng)力學(xué)性質(zhì)，則采用半經(jīng)驗(yàn)與DFT相結(jié)合的方法，先通過運(yùn)算速度較快的鍵價(jià)和方法[44]對(duì)材料的離子輸運(yùn)性質(zhì)進(jìn)行初步篩選，將具有連通的離子通道并且活化能較低的化合物找出來，再選擇更為精確的DFT計(jì)算，采用NEB方法[46]計(jì)算準(zhǔn)確的遷移勢壘，對(duì)于包含無序占位或部分占位的化合物，則考慮采用第一性原理的分子動(dòng)力學(xué)(FPMD)[49]來進(jìn)行計(jì)算。DFT計(jì)算主要采用VASP程序[37]，基于鍵價(jià)和的離子通道計(jì)算則采用自行編寫的BVpath進(jìn)行。同時(shí)，我們根據(jù)圖3的篩選流程編寫了相應(yīng)腳本，實(shí)現(xiàn)了對(duì)ICSD數(shù)據(jù)庫中的晶體結(jié)構(gòu)數(shù)據(jù)自動(dòng)轉(zhuǎn)化、對(duì)計(jì)算軟件的自動(dòng)調(diào)用、以及根據(jù)篩選條件將計(jì)算結(jié)果自動(dòng)分類存儲(chǔ)的過程。

5高通量計(jì)算篩選鋰電池電極材料

提高鋰電池能量密度的主要方法之一是提高電極材料的儲(chǔ)鋰容量。電極材料的儲(chǔ)鋰容量與可轉(zhuǎn)移電子數(shù)、可轉(zhuǎn)移離子數(shù)有關(guān)。目前幾種重要的正極材料(LiCoO2、LiMn2O4、LiFePO4)以及與其類似組成的材料，其電子轉(zhuǎn)移都是通過過渡金屬的氧化還原來實(shí)現(xiàn)。根據(jù)等式(1)和(2)可以從熱力學(xué)數(shù)據(jù)計(jì)算各種化學(xué)反應(yīng)對(duì)應(yīng)的能量密度。

根據(jù)化合物在25 ℃時(shí)的標(biāo)準(zhǔn)吉布斯自由焓數(shù)據(jù)以及密度數(shù)據(jù)[39]，Zu等人[38]計(jì)算了1 172種化學(xué)反應(yīng)所對(duì)應(yīng)的能量密度。圖4顯示了部分鋰電池體系的計(jì)算結(jié)果，當(dāng)反應(yīng)物的形成能較低而產(chǎn)物的形成能較高時(shí)，化學(xué)反應(yīng)對(duì)應(yīng)的能量密度較高。計(jì)算的所有體系中，Li/F2體系具有最高的能量密度，6 294 Wh·kg-1。其次是Li/O2體系，其能量密度為5 217 Wh·kg-1。除此之外，一些轉(zhuǎn)化反應(yīng)也顯示出較高的能量密度。Zu等人[38]進(jìn)一步計(jì)算了金屬鋰與各種物質(zhì)發(fā)生轉(zhuǎn)化反應(yīng)的熱力學(xué)平衡電壓，部分結(jié)果如圖5所示。對(duì)于二元過渡金屬化合物與鋰的轉(zhuǎn)化反應(yīng)，對(duì)于同種過渡金屬，當(dāng)轉(zhuǎn)移電子數(shù)目相同時(shí)，熱力學(xué)平衡電壓值具有如下規(guī)律：氟化物>氧化物>硫化物>氮化物>磷化物,并且高氧化態(tài)的材料通常顯示更高的熱力學(xué)平衡電壓值。

圖4　根據(jù)能斯特方程計(jì)算的各種鋰電池體系的能量密度[38]Fig.4　Calculated energy densities of different lithium battery systems according to the Nernst equation[38]

圖5　計(jì)算的金屬鋰與各種二元過渡金屬化合物發(fā)生轉(zhuǎn)化反應(yīng)的熱力學(xué)平衡電壓[38]Fig.5　Calculated thermodynamic equilibrium voltage of conversion reactions between binary transition metal compounds and lithium[38]

根據(jù)熱力學(xué)數(shù)據(jù)計(jì)算的能量密度和熱力學(xué)平衡電壓，可有效用于電極材料(特別是正極材料)的初步篩選。以目前商業(yè)應(yīng)用的正極材料LiCoO2為參考，將儲(chǔ)鋰能力高于此標(biāo)準(zhǔn)的化合物作為進(jìn)一步篩選的候選材料。由于減少電極材料循環(huán)過程中的晶格體積變化對(duì)提高電池循環(huán)壽命至關(guān)重要，因此在接下來的研究中，將在高通量計(jì)算里引入第一性原理方法，對(duì)材料儲(chǔ)鋰前后的體積變化和熱力學(xué)穩(wěn)定性進(jìn)行計(jì)算，以實(shí)現(xiàn)對(duì)上述候選材料的進(jìn)一步篩選。

6高通量計(jì)算篩選鋰電池固態(tài)電解質(zhì)材料

從長遠(yuǎn)來看，開發(fā)全固態(tài)鋰離子電池是解決目前鋰電池發(fā)展瓶頸的一個(gè)思路。首先，全固態(tài)鋰離子電池的安全性能非常好。當(dāng)前鋰離子電池采用液態(tài)電解液體系，采用的可燃的有機(jī)溶劑是電池體系中最大的安全隱患。而全固態(tài)鋰離子電池采用固態(tài)電解質(zhì)材料，安全問題可以從根源上一次性解決。其次，全固態(tài)鋰離子電池可以采用金屬Li為負(fù)極，其理論容量是鋰電池體系的最高極限，而且其電位也是所有負(fù)極材料中最低的(這樣可以提高電池體系的整體的電壓)。因此，從總體上可以大大提高電池體系的能量密度。再次，由于負(fù)極材料采用了金屬Li，這就給了正極材料更寬的選擇——正極材料可以不作為Li源的提供者。

然而，要開發(fā)全固態(tài)鋰離子電池體系，面臨著許多的挑戰(zhàn)，其中最關(guān)鍵的挑戰(zhàn)在于開發(fā)出可以應(yīng)用的固體電解質(zhì)材料體系。固體電解質(zhì)需要兼具高離子電導(dǎo)率、高穩(wěn)定性、高機(jī)械強(qiáng)度等特點(diǎn)。目前的固體電解質(zhì)體系有3類：①無機(jī)晶態(tài)固體電解質(zhì)材料，包括石榴石(Garnets)結(jié)構(gòu)材料[51]、金紅石(Rutiles)結(jié)構(gòu)材料[52]、沸石(Zeolites)類結(jié)構(gòu)材料[53]等；②玻璃態(tài)的(Glassy)無定形態(tài)材料體系，這些材料通常都是Li的S化物和P化物(Li2S、P2S5、 Li3P7、Li7P3S11、 Li4P2S7等)[12]；③ 聚合物(PEO)固態(tài)電解質(zhì)材料體系[54]。這些材料體系都有其各自的優(yōu)點(diǎn)和缺點(diǎn)，但目前尚沒有獲得可實(shí)際應(yīng)用的固體電解質(zhì)材料，因此發(fā)現(xiàn)新的固態(tài)電解質(zhì)材料體系勢在必行。

離子輸運(yùn)性質(zhì)是決定固態(tài)電解質(zhì)性能的最主要因素，因此計(jì)算材料的鋰離子通道及活化能大小是進(jìn)行電解質(zhì)材料篩選的首要步驟。初步的篩選從計(jì)算精度較低但運(yùn)算速度很快的鍵價(jià)和方法開始。BVpath程序結(jié)合了改進(jìn)的鍵價(jià)和理論[44]，計(jì)算流程如圖6所示。從晶體結(jié)構(gòu)數(shù)據(jù)庫中讀出化合物中的原子坐標(biāo)、占據(jù)率和各原子的名義電荷，據(jù)此計(jì)算鋰離子與陰離子之前的鍵能以及鋰離子與陽離子之間的庫倫排斥勢，以二者之和來做等值面，找出能使鋰離子通道聯(lián)通的最小能量等值面，從而確定擴(kuò)散路徑和勢壘。

圖6　BVpath軟件計(jì)算鋰離子擴(kuò)散路徑及勢壘的流程圖Fig.6　Flow diagram of BVpath code based on bond valence method

圖7　用鍵價(jià)和理論計(jì)算得到的Li2MnO3材料中的鋰離子擴(kuò)散路徑:(a~c)為1.08 eV對(duì)應(yīng)的等值面，顯示了Li+在鋰層內(nèi)的遷移；(d~f)為1.55 eV對(duì)應(yīng)的等值面，顯示了Li+在過渡金屬層與鋰層之間的遷移Fig.7　The migration pathway of Li+ in Li2MnO3 calculated by Bond Valence method:(a~c)the isosurface of 1.08 eV indicates the migration of Li+ within the lithium layer and (d~f) the isosurface of 1.55 eV indicates the migration of a Li+ between the transition metal layer and lithium layer

圖7是用鍵價(jià)和方法計(jì)算的Li2MnO3材料中的鋰離子輸運(yùn)通道。以不同能量值畫出的等值面顯示了Li+在不同遷移路徑上的能量勢壘。從圖(a)~(c)中可以看到，Li+在鋰層之間的遷移的勢壘為1.08 eV，低于Li+在過渡金屬層和鋰層之間的遷移勢壘1.55 eV(圖(d)~(f ))。鍵價(jià)和方法得到的鋰離子輸運(yùn)路徑與采用NEB方法的DFT計(jì)算結(jié)果[55]一致,二者都揭示出了沿著ab方向的層內(nèi)傳導(dǎo)和沿著c方向的層間傳導(dǎo)兩種路徑；鍵價(jià)和方法計(jì)算得到的兩種路徑的勢壘值由于受到計(jì)算方法的局限，在精度上無法與DFT計(jì)算比擬，但仍顯示出層內(nèi)傳導(dǎo)比層間傳導(dǎo)更易于發(fā)生，得到了與DFT計(jì)算結(jié)果一致的趨勢。

通過對(duì)Li2MnO3計(jì)算結(jié)果的分析，表明了鍵價(jià)和方法在揭示晶體結(jié)構(gòu)中可能存在的鋰離子輸運(yùn)通道和半定量計(jì)算其遷移勢壘的有效性。因此，這種方法可以更廣泛的用于鋰電池材料離子輸運(yùn)性質(zhì)的研究?？紤]到離子輸運(yùn)性質(zhì)是固體電解質(zhì)材料研發(fā)中需首要考慮的性質(zhì)，我們編寫了調(diào)用BVpath程序的高通量計(jì)算流程，對(duì)ICSD數(shù)據(jù)庫中800余種含鋰化合物的鋰離子輸運(yùn)通道和遷移能量勢壘進(jìn)行了計(jì)算[56-57]。圖8給出了7種分別屬于不同晶系的化合物的計(jì)算結(jié)果。圖中淺灰色的包絡(luò)面是能量勢壘值所對(duì)應(yīng)的等值面，其中包裹的區(qū)域?yàn)長i+可能運(yùn)動(dòng)的區(qū)域。等值面所包裹的通道的最窄處為空間中的能量鞍點(diǎn)，是遷移路徑中的能量最高點(diǎn)。根據(jù)鍵價(jià)和方法計(jì)算得到的離子通道和勢壘值可以用于判斷材料中鋰離子運(yùn)動(dòng)趨勢及其難易程度，從而用于從晶體結(jié)構(gòu)數(shù)據(jù)庫中初步篩選固體電解質(zhì)。

圖8　用鍵價(jià)和理論計(jì)算得到的7種材料的鋰離子擴(kuò)散通道與勢壘:(a)Li2CO3,Ea=0.60 eV;(b) Li3B7O12,Ea=1.42 eV;(c) Li6Ge2O7,Ea=0.90 eV;(d) γ-LiAlO2,Ea=0.85 eV;(e) α-LiAlO2,Ea=1.89 eV;(f)LiAlSiO2 ,Ea=0.38 eV;(g) LiF ,Ea=1.07 eV[57]Fig.8　The pathways and energy barriers of Li+ migration for seven compounds:(a)Li2CO3,Ea=0.60 eV;(b) Li3B7O12,Ea=1.42 eV;(c) Li6Ge2O7,Ea=0.90 eV;(d) γ-LiAlO2,Ea=0.85 eV;(e) α-LiAlO2,Ea=1.89 eV;(f)LiAlSiO2 ,Ea=0.38 eV;(g) LiF ,Ea=1.07 eV[57]

在進(jìn)一步的高通量計(jì)算中，引入第一性原理方法，計(jì)算各種備選材料的電子結(jié)構(gòu)，獲得能隙寬度的數(shù)據(jù)，以此估算固體電解質(zhì)材料的電化學(xué)窗口寬度，作為下一步的篩選條件。這些進(jìn)一步的高通量計(jì)算和材料篩選工作正在進(jìn)行中，預(yù)期將為鋰電池中固體電解質(zhì)材料的發(fā)現(xiàn)提供有價(jià)值的參考。

7結(jié)語

應(yīng)用高通量計(jì)算方法探索和設(shè)計(jì)鋰電池新材料，是未來材料科學(xué)研究和開發(fā)的一條重要途徑。借助無機(jī)材料晶體結(jié)構(gòu)數(shù)據(jù)庫中現(xiàn)有的和衍生的結(jié)構(gòu)，設(shè)定合理的限制條件，有望篩選出下一代鋰電池的備選材料。作者團(tuán)隊(duì)自主編寫了可實(shí)現(xiàn)晶體結(jié)構(gòu)選擇、輸入文件生成、VASP軟件調(diào)用、計(jì)算結(jié)果分類存儲(chǔ)的自動(dòng)化計(jì)算流程，并編寫了基于鍵價(jià)和方法的鋰離子輸運(yùn)計(jì)算軟件BVpath，同時(shí)也實(shí)現(xiàn)了高通量計(jì)算過程對(duì)該軟件的自動(dòng)調(diào)用。

從材料研發(fā)的角度，這種方法有可能從現(xiàn)有晶體結(jié)構(gòu)中找到或構(gòu)建出新晶體結(jié)構(gòu)，并評(píng)估新結(jié)構(gòu)存在的可能性以及作為鋰電池材料的性能，更快找到高性能、高安全性的正極、負(fù)極與電解質(zhì)材料，促進(jìn)鋰電池產(chǎn)業(yè)的發(fā)展，加速鋰電池材料的開發(fā)過程。

從基礎(chǔ)研究的角度，高通量計(jì)算得到的大量數(shù)據(jù)，將有助于理解與鋰電池有關(guān)的現(xiàn)象，如結(jié)構(gòu)與離子擴(kuò)散、儲(chǔ)鋰容量之間的本質(zhì)關(guān)系。計(jì)算得到的晶體結(jié)構(gòu)、電子結(jié)構(gòu)和各種性質(zhì)，與現(xiàn)有的實(shí)驗(yàn)數(shù)據(jù)匯集到一起，將形成一個(gè)龐大的數(shù)據(jù)庫，為材料學(xué)家及物理學(xué)家提供全面而寶貴的數(shù)據(jù)參考。進(jìn)一步借助統(tǒng)計(jì)方法進(jìn)行數(shù)據(jù)挖掘，還有可能找到?jīng)Q定材料性質(zhì)的因素，如相關(guān)的結(jié)構(gòu)單元或是通過統(tǒng)計(jì)方法定義的關(guān)聯(lián)函數(shù)等，從而找到背后的物理原因。

可以預(yù)見在未來幾年內(nèi)，應(yīng)用高通量計(jì)算方法，結(jié)合實(shí)驗(yàn)研究，將大大加速我們對(duì)鋰電池材料中各種問題的理解，揭示材料的結(jié)構(gòu)-性能關(guān)系及內(nèi)在物理規(guī)律，縮短材料從研發(fā)到應(yīng)用的全過程，為開發(fā)全新電池材料乃至電池體系提供有力支持。

參考文獻(xiàn)References

[1]Goodenough J B, Kim Y. Challenges for Rechargeable Li Batteries[J].ChemistryofMaterials, 2010 (22): 587-603.

[2]Tarascon J M, Armand M. Issues and Challenges Facing Rechargeable Lithium Batteries[J].Nature, 2001 (414): 359-367.

[3]Karden E, Ploumen S, Fricke B,etal. Energy Storage Devices for Future Hybrid Electric Vehicles[J].JPowerSources, 2007 (168): 2-11

[4]Dunn B, Kamath H, Tarascon J M. Electrical Energy Storage for the Grid: A Battery of Choices[J].Science, 2011 (334): 928-935

[5]Li H, Wang Z X, Chen L Q,etal. Research on Advanced Materials for Li-Ion Batteries[J].AdvancedMaterials, 2009 (21): 4 593-4 607.

[6]Chen Y, Wang G X, Konstantinov K,etal. Synthesis and Characterization of LiCoxMnyNi1-x-yO2as a Cathode Material for Secondary Lithium Batteries[J].JournalofPowerSources, 2003 (119): 184-188.

[7]Palacin M R, Chabre Y, Dupont L,etal. On the Origin of the 3.3 and 4.5 V Steps Observed in LiMn2O4-Based Spinels[J].JournalofTheElectrochemicalSociety, 2000 (147): 845-853.

[8]Hu Y, Kong W, Li H,etal. Experimental and Theoretical Studies on Reduction Mechanism of Vinyl Ethylene Carbonate on Graphite Anode for Lithium Ion Batteries[J].ElectrochemistryCommunications, 2004 (6): 126-131.

[9]He Y, Yu X, Wang Y,etal. Alumina-Coated Patterned Amorphous Silicon as the Anode for a Lithium-Ion Battery with High Coulombic Efficiency[J].AdvancedMaterials, 2011 (23): 4 938-4 941.

[10]Xu K. Nonaqueous Liquid Electrolytes for Lithium-Based Rechargeable Batteries[J].ChemicalReviews, 2004 (104): 4 303-4 417.

[11]Quartarone E, Mustarelli P. Electrolytes for Solid-State Lithium Rechargeable Batteries: Recent Advances and Perspectives[J].ChemicalSocietyReviews, 2011 (40): 2 525-2 540.

[12]Mizuno F, Hayashi A, Tadanaga K,etal. New, Highly Ion-Conductive Crystals Precipitated from Li2S-P2S5Glasses[J].AdvancedMaterials, 2005 (17): 918-921.

[13]Kamaya N, Homma K, Kamakawa Y,etal. A Lithium Superionic Conductor[J].NatureMaterials, 2011 (10): 682-686.

[14]Mo Y, Ong S P, Ceder G,etal. First Principles Study of the Li10GeP2S12Lithium Super Ionic Conductor Material[J].ChemistryofMaterials, 2012 (24): 15-17.

[15]Lu J, Li L, Park J B,etal. Aprotic and Aqueous Li-O2Batteries[J].ChemicalReviews, 2014 (114): 5 611-5 640.

[16]Bruce P G, Freunberger S A, Hardwick L J,etal. Li-O2and Li-S Batteries with High Energy Storage[J].NatureMaterials, 2012 (11): 19-29.

[17]Liu Y L, Wang R, Lyu Y C,etal.Rechargeable Li/CO2-O2(2:1) Battery and Li/CO2Battery[J].EnergyEnvironmentScience, 2014 (7): 677-681.

[18]Shi S Q, Liu L J, Ouyang C Y,etal. Enhancement of Electronic Conductivity of LiFePO4by Cr Doping and Its Identification by First-Principles Calculations[J].PhysicalReviewB, 2003 (68): 195 108.

[19]Ouyang C Y, Shi S Q, Wang Z X,etal. First-Principles Study of Li Ion Diffusion in LiFePO4[J].PhysicalReviewB, 2004 (69): 104 303.

[20]Ouyang C Y, Shi S Q, Wang Z X,etal. Ab Initio Molecular-Dynamics Studies on LixMn2O4as Cathode Material for Lithium Secondary Batteries[J].EurophysicsLetters, 2004 (67): 28-34.

[21]Wei Y J, Xu X G, Wang C Z,etal. Electronic Structure of Cubic Li(Fe0.1Mn1.9)O4Studied with M?ssbauer Spectroscopy and First-Principles Calculation[J].AppliedPhysicsLetters, 2003 (83): 1 791-1 793.

[22]Van der Ven A, Ceder G, Asta M,etal. First-Principles Theory of Ionic Diffusion with Nondilute Carriers[J].PhysicalReviewB, 2001(64): 184 307.

[23]Ouyang C Y, Zhong Z Y, Lei M S,etal. Ab Initio Studies of Structure and Electronic Properties of Li4Ti5O12Spinel[J].ElectrochemistryCommunications, 2007(9): 1 107-1 112.

[24]Ouyang C Y, Zeng X M,ljivancanin ?,etal. Oxidation States of Mn Atoms at Clean and Al2O3-Covered LiMn2O4(001) Surfaces[J].JournalofPhysicalChemistryC,2010(114): 4 756-4 759.

[25]Ten Bosch J R, Grody W W. Keeping Up with the Next Generation[J].TheJournalofMolecularDiagnostics, 2008 (10): 484-492.

[26]Tucker T, Marra M, Friedman J M. Massively Parallel Sequencing: The Next Big Thing in Genetic Medicine[J].TheAmericanJournalofHumanGenetics, 2009 (85): 142-154.

[27]Howe D, Costanzo M, Fey P,etal. Big Data: The Future of Biocuration[J].Nature,2008 (455): 47-50.

[29]Andersson M P, Bligaard T, Kustov A,etal. Toward Computational Screening in Heterogeneous Catalysis: Pareto-Optimal Methanation Catalysts[J].JournalofCatalysis, 2006 (239): 501-506.

[30]Curtarolo S, Setyawan W, Wang S D,etal. AFLOWLIB.ORG: A Distributed Materials Properties Repository from High-Throughput Ab Initio Calculations[J].ComputationalMaterialsScience, 2012 (58): 227-235.

[31]Setyawan W, Gaumé R M, Lam S,etal. High-Throughput Combinatorial Database of Electronic Band Structures for Inorganic Scintillator Materials[J].ACSCombinatorialScience, 2011 (13): 382-390.

[32]Yang K, Setyawan W, Wang S,etal. A Search Model for Topological Insulators with High-Throughput Robustness Descriptors[J].NatureMaterials, 2012 (11): 614-619.

[33]Jain A, Hautier G, Moore C J,etal. A High-Throughput Infrastructure for Density Functional Theory Calculations[J].ComputationalMaterialsScience, 2011 (50): 2 295-2 310.

[34]Hautier G, Jain A, Chen H,etal. Novel Mixed Polyanions Lithium-Ion Battery Cathode Materials Predicted by High-Throughput Ab Initio Computations[J].JournalofMaterialsChemistry, 2011 (21): 17 147-17 153.

[35]Hautier G, Jain A, Ong S P,etal.Phosphates as Lithium-Ion Battery Cathodes: An Evaluation Based on High-Throughput Ab Initio Calculations[J].ChemistryofMaterials, 2011 (23): 3 495-3 508.

[36]Fujimura K, Seko A, Koyama K,etal. Accelerated Materials Design of Lithium Superionic Conductors Based on First-Principles Calculations and Machine Learning Algorithms[J].AdvancedEnergyMaterials, 2013 (3): 980-985.

[37]Kresse G, Furthmuller J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set[J].ComputionalMaterialsScience, 1996 (6): 15-50.

[38]Zu C X, Li H. Thermodynamic Analysis on Energy Densities of Batteries[J].EnergyEnvironmentalScience, 2011 (4): 2 614-2 624.

[39]David R.CRCHandbookofChemistryandPhysics, 90thedn[M].Boca Raton:CRC Press, 2009.

[40]Ceder G, Chiang Y-M, Sadoway D R,etal. Identification of Cathode Materials for Lithium Batteries Guided by First-Principles Calculations[J].Nature, 1998 (392): 694-696.

[41]Anurova N A, Blatov V A. Analysis of Ion-Migration Paths in Inorganic Frameworks by Means of Tilings and Voronoi-Dirichlet Partition: A Comparison[J].ActaCrystallographicaSectionB:StructuralScience, 2009 (65): 426-434.

[42]Nuspl G, Takeuchi T, Weiss A,etal. Lithium Ion Migration Pathways in LiTi2(PO4)3and Related Materials[J].JournalofAppliedPhysics, 1999 (86): 5 484-5 491.

[43]Filso M O, Turner M J, Gibbs G V,etal. Visualizing Lithium-Ion Migration Pathways in Battery Materials[J].Chemistry-AEuropeanJournal, 2013 (19): 15 535-15 544.

[44]Brown I D. Recent Developments in the Methods and Applications of the Bond Valence Mode[J].ChemicalReview, 2009 (109): 6 858-6 866.

[45]Adams S, Rao R P. High Power Lithium Ion Battery Materials by Computational Design[J].PhysicaStatusSolidiA, 2011 (208): 1 746-1 753.

[46]Sheppard D, Xiao P, Chemelewski W,etal. A Generalized Solid-State Nudged Elastic Band Method[J].JournalofChemicalPhysics, 2012,136(7): 074 103.

[47]Lepley N D, Holzwarth N A W, Du Y A. Structure, Li+Mobilities, and Interfacial Properties of Solid Electrolytes Li3PS4and Li3PO4from First Principles[J].PhysicalReviewB, 2013,88(10): 104 103.

[48]Singleton J.BandTheoryandElectronicPropertiesofSolids[M].Oxford University Press, 2001.

[49]Martin R M.ElectronicStructure:BasicTheoryandPracticalMethods[M]. Cambridge University Press, 2004.

[50]Inorganic Crystal Structure Database, ICSD. Karlsruhe: Fachinformationszentrum,2008.

[51]Murugan R, Thangadurai, Weppner W. Fast Lithium Ion Conduction in Garnet-Type Li7La3Zr2O12[J].AngewandteChemieInternationalEdition, 2007 (46): 7 778-7 781.

[52]Yates D E, Healy T W. Titanium Dioxide-Electrolyte Interface[J].JournaloftheChemicalSociety,FaradayTransactions1: Physical Chemistry in Condensed Phases, 1980 (76): 9-18.

[53]Noto V D. Zeolitic Inorganic-Organic Polymer Electrolyte Based on Oligo(Ethylene Glycol) 600 K2PdCl4and K3Co(CN)6[J].JournalofPhyscialChemistryB, 2000 (104): 10 116-10 125.

[54]Jacob M M E, Hackett E, Giannelis E P. From Nanocomposite to Nanogel Polymer Electrolytes[J].JournalofMaterialsChemistry, 2002,1(1): 1-5.

[55]Xiao R, Li H, Chen L Q. Density Functional Investigation on Li2MnO3[J].ChemistryofMaterials, 2012 (24): 4 242-4 251.

[56]Gao J, Chu G, He M,etal. Screening Possible Solid Electrolytes by Calculating the Conduction Pathways Using Bond Valence Method[J].ScienceChina:Physics,MechanicsandAstronomy, 2014,57(8): 1 526-1 535.

[57]Zhang Shu(張舒), Wang Shaofei(王少飛), Ling Shigang(凌仕剛),etal. 鋰離子電池基礎(chǔ)科學(xué)問題(X)—— 全固態(tài)鋰離子電池[J].EnergyStorageScienceandTechnology(儲(chǔ)能科學(xué)與技術(shù)), 2014,3(4): 376-394.

(編輯蓋少飛)

Application of High-Throughput Calculations for ScreeningLithium Battery Materials

LING Shigang1, GAO Jian1, CHU Geng1, HUANG Jie1, XIAO Ruijuan1,

OUYANG Chuying2, LI Hong1, CHEN Liquan1

(1.Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China)

(2.Jiangxi Normal University, Nanchang 330027, China)

Abstract：Advanced cell technology is the world’s leading competition in the next ten years. Lithium batteries are the key to the success of portable electronic products, electrical vehicle and the best choices for the energy storage devices of wind and solar power because of their high voltage, high specific energy density, rapid recharge capability, and wide working temperature range. The development of the next generation lithium battery materials is necessary to achieve higher energy density, better safety, and faster charge and discharge speed. The development of new materials is possible to be accelerated by high-throughput calculations and reasonable screening procedures. In this work, using automated calculation flow designed by ourselves, the initial screening of electrode matreials and solid electrolytes is performed from the Inorganic Crystal Structure Database according to the application requirments. The calculation results and experimental data integrate into a large database, which is hopeful to reveal the realtionship between structure and property, thus shortening the process of materials from discovering to application and providing support for the development of new battery materials and even new battery systems.

Key words：lithium battery material; high-throughput calculation; thermodynamic calculation; bond valence method

中圖分類號(hào)：O646.2+1

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1674-3962(2015)04-0272-10

DOI:10.7502/j.issn.1674-3962.2015.04.03

通訊作者：肖睿娟，女，1980年生，副研究員，Email: rjxiao@aphy.iphy.ac.cn

基金項(xiàng)目：國家自然科學(xué)基金(11234013，51325206，51172274);北京市科委重點(diǎn)項(xiàng)目;科技部973項(xiàng)目(2012CB932900)

收稿日期：2014-09-04