锂电池领域常用计算方法及计算技巧

随着化石能源的日渐枯竭以及由于燃烧化石能源而带来的种种环境问题，发展下一代新型清洁能源已经迫在眉睫。太阳能，风能虽然由于起环境友好型而受到人们追捧，但其间歇性也限制了它们的应用。另一方面，可充电电池或者超级电容器，由于其高能量密度，便携性，相对安全性，已经逐步发展为下一代新型能源中的佼佼者。作为可充电电池的一种，锂离子电池在近几十年来得到飞速发展，并且已经广泛地应用于汽车，手机，电子，工业等等方面。近些年来，锂离子电池的主要研究方向则是发展高能量密度，高电压型正极材料，以及实现锂金属负极的锂硫电池(Li-S)或锂空电池(Li-O2)[1]。

第一性原理(First-Principlemethod)作为一种从头(abinitio)算法，被广泛应用于材料设计，材料预测，以及解释实验等方面。从薛定谔方程开始，经过JohnPope，WalterKohn等人的发展逐步到可以实现计算机编程的DFT(DensityFunctionalTheory)。只需要几个经典参数，DFT便可以计算出材料的大部分性能。进一步结合Born–Oppenheimer假设，DFT也可以用来模拟分子动力学，即Abinitiomoleculardynamics(AIMD)。在锂电等相关领域，DFT得到广泛的应用。由于其准确性，可以有效的地指导实验，从而节约成本。

本文就DFT计算在锂电相关领域的应用作一探讨，为读者梳理出几中比较常见的计算方法，并稍加提及一些计算技巧。由于笔者经常利用VASP(ViennaAbinitioSimulationPackage)作材料相关计算，以下的讨论基于该软件，当然也可以推广到其他的商业软件或者开源代码。

计算材料的充放电平台

锂离子电池又称作“摇摇椅(rockingchair)”。这是因为在锂离子电池中，锂离子通过电解质在正负极间穿梭。充电时，锂离子从正极材料(如LiCoO2)晶格中脱出，嵌入到负极材料(如石墨)中，使得负极富锂，正极贫锂;放电时，由于负极侧锂的化学势较高，锂离子从负极脱出，插入到正极材料晶格中，从而完成一个循环。其中电解质只是充当传导锂离子的作用。好的电解质需要有高的锂离子电导率，低的电子电导率，还要有较宽的电化学窗口以防止氧化还原反应。

在平衡状态下，锂离子电池的嵌入电压定义为正负极两端的Li化学势之差，即

其中z为转移电荷，F为法拉第常数。积分上式则有

其中为嵌入反应过程体系的总吉布斯自由能变化。在不考虑温度及体积变化的情况下，可以近似为。以锂金属为负极，正极材料为过渡金属氧化物LiMO2的平衡电压则可以计算为

其充电时化学反应式为

其中χ1>χ2，且相对于立方相金属锂。

因此，只要计算出脱出锂离子前后体系的能量，就可以计算出充电工作平台。当然，上式也可以用来计算钠离子电池的电压平台。如Ceder课题组[2]就针对NaxMnO2做了计算并和实验数据比较。可以看到计算结果与很好地与实验放电平台匹配。

此外，考虑到3d轨道的过渡金属，如Co，Fe，Mn等，一般的基于GGA的交换关联泛函已经不能够满足计算要求。这时需要考虑对过渡金属元素加U，[3]即GGA+U。Ceder课题组比较了GGA，GGA+U和更高级的HSE对锂离子嵌入电压平台的计算，结果表明GGA+U能有效提高计算精度。GGA+U的结果和HSE的结果十分接近，表明GGA+U能有效的考虑到3d轨道的影响。在实际应用中，U的数值需要拟合，也可以直接从文献中获取。如对于Fe，一般取U为4.0~4.3eV。

锂离子的传输

锂离子的传输对于电解质和电极材料都有着十分重要的意义，尤其是固态电解质。作为固态电解质中的一种，Li10GeP2S12正是由于其较高的锂离子电导率(12mS/cm)而受到追捧。[4]利用DFT可以有效地计算锂离子在固体中的传输扩散，其中NEB(Nudgedelasticband)和AIMD是两种比较常用的计算方法。

NEB可以用来计算过渡态(Transitionstate)，也可以用来计算锂离子在体系中的扩散。以层状Si为例，[5]锂离子在其表面的扩散可以有以下几种方式，分别对应单层(SL)Si和双层(DL)Si结构。扩散能磊越小，意味着锂离子越容易扩散。可以看出不同的扩散方式以及不同浓度下锂离子的扩散能磊各有不同。在实际计算中，通常先优化初始和末态的结构，然后利用vtst的脚本对中间结构进行差值。提醒读者需要在插完值后对中间结构进行检查。

另外一种常用的计算方法是AIMD。由于其基于DFT，所以力场函数(Forcefiled)不用外加设定，这相比于经典的MD要更加方便。然而，由于其计算量要比经典MD大很多，所以一般的AIMD只适用于比较小的体系(几百个原子)。基于AIMD，可以利用Einstein关系式计算锂离子的扩散系数[6]，即

其中，MSD为锂离子的总平均自由移动距离，d为维数。

在实际应用中，对上式进行拟合时，需要特别注意起始段和结尾段。Mo等人针对Li7La3Zr2O12，Li10GeP2S12，Li1.33Ti1.67Al0.33(PO4)3材料的计算，提出以下的判据：MSD<0.5a2及?t>0.7ttot的数据舍去。其中a为锂原子位间距，ttot为总的MD时长。

根据Nernst–Einstein关系式，电导率可以计算为

其中N为锂离子数，q为锂离子带电荷，V为体系体积，kB为Boltzmann常数，T为温度。一般的MD都是在高温下(~1000K)计算，然后经过Arrhenius关系式插值得到室温下的锂离子扩散系数，再进一步利用Nernst–Einstein关系式得到电导率。

注意这里只有扩散系数D对温度存在这样的函数，而则应考虑一起。

锂离子或小分子的吸附

做为正极材料，硫单质具有较高的比容量(1673mAh/g)。因而发展下一代可循环利用的商业锂硫电池已经提上日程。然而在反应中，锂离子和硫单质形成的多硫化锂溶于电解质并在正负极间穿梭，导致活性物质损失以及增加电解液的粘度。此外，溶解的多硫化物还会扩散到负极，与负极反应破坏电解质界面膜(SEI)，增大界面电阻。因此在实际应用中，有效的限制多硫化锂的穿梭，使其集中在正极侧是一种有效的方法。

N掺杂石墨烯可以有效地抑制多硫化物的脱附，从而控制穿梭效应。Cheng[7]等人研究了不同N掺杂结构石墨烯对多硫化物的吸附能，并得出结论认为石墨烯掺杂吡啶型聚类N可以有效地吸附多硫化物，如下图左侧所示。

在实际计算中，吸附能定义如下

其中对应于Li2Sx在基底(石墨烯)上吸附之后的结构的，是单独优化基地得到的能量，是多硫化锂小分子的能量。

在实际计算中，考虑到S原子和C原子之间比较强的范德华作用力，一般需要在DFT的计算基础上引入范德华修正。Zhang等[8]研究了考虑范德华修正和不考虑范德华修正下多硫化锂和石墨烯的作用，发现引入范德华修正后吸附能会提高，而且趋势也有所变化。在不考虑范德华修正的情况下，多硫化合物和石墨烯几乎没有作用，吸附能在0.1~0.3eV。在vasp中，考虑范德华修正的相关参数为IVDW。

材料稳定性的计算

材料稳定性是锂离子电池能长时间运行的得力保证，同样也是各种材料可以在实验室合成的基础。由热力学第二定律可知，相对于组分来说，平均自由能是一个凸包(convexhull)，如下图。其中坐落在convexhull线上的点对应稳定相，如下图中的红点;而落在convexhull线以上的点则不稳定，如下图中的蓝点。对于非稳定相，定义其距离convexhull线的垂直高度为energyabovehull，即Ehull。该数值越大，则说明该相越不稳定。一般认为Ehull小于100meV/atom属于亚稳定，而Ehull大于该阈值则认为不稳定。因此计算不同组分下的自由能，并构建convexhull可以成功预测材料的稳定性。

在实际计算中，可以先构造出不同组分的机构，然后利用DFT进行结构优化，得到吉布斯自由能G。以平均到每个原子的自由能为标准，是各个组分的函数。以Li-Fe-P-O为例，构建好的convexhull投影到二维平面中就是常见的三角相图。

Ong等[9]对LiFePO4的相图进行研究，并在此基础上发展出pymatgen软件包。该软件包基于Python并结合MaterialProject数据库，可以很方便的进行一系列高通量计算。

结论

第一性原理在锂电相关领域里面发挥着很重要的作用。利用DFT计算，可以有效地预测材料的电子性质，如带宽\费米能级，HUMO和LUMO等重要特性。此外，结合DFT和AIMD，可以有效的预测锂离子在材料中的传输及扩散，从而预测其锂离子电导率。最后，结合已有的材料数据库和相关开发包，可以方便的计算材料的稳定性。