锂电池的计算计算｜常用锂电参数计算公式|上海羊羽卓进出口贸易有限公司

计算｜常用锂电参数计算公式

来源：江子才

（1）电极材料的理论容量

电极材料理论容量，即假定材料中锂离子全部参与电化学反应所能够提供的容量，其值通过下式计算：

其中，法拉第常数(F)代表每摩尔电子所携带的电荷，单位C/mol，它是阿伏伽德罗数NA=6.02214 ×1023mol-1与元电荷e=1.602176 × 10-19 C的积，其值为96485.3383±0.0083 C/mol

故而，主流的材料理论容量计算公式如下:

LiFePO4摩尔质量157.756 g/mol，其理论容量为：

同理可得：三元材料NCM(1:1:1)(LiNi1/3Co1/3Mn1/3O2 ) 摩尔质量为96.461g/mol，其理论容量为278 mAh/g,LiCoO2摩尔质量97.8698 g/mol，如果锂离子全部脱出，其理论克容量274 mAh/g.

石墨负极中，锂嵌入量最大时，形成锂碳层间化合物，化学式LiC6，即6个碳原子结合一个Li。6个C摩尔质量为72.066 g/mol，石墨的最大理论容量为：

对于硅负极，由5Si+22Li++22e- ↔ Li22Si5 可知， 5个硅的摩尔质量为140.430 g/mol，5个硅原子结合22个Li，则硅负极的理论容量为：

这些计算值是理论的克容量，为保证材料结构可逆，实际锂离子脱嵌系数小于1，实际的材料的克容量为：材料实际克容量=锂离子脱嵌系数 × 理论容量

（2）电池设计容量

电池设计容量=涂层面密度×活物质比例×活物质克容量×极片涂层面积

其中，面密度是一个关键的设计参数，主要在涂布和辊压工序控制。压实密度不变时，涂层面密度增加意味着极片厚度增加，电子传输距离增大，电子电阻增加，但是增加程度有限。厚极片中，锂离子在电解液中的迁移阻抗增加是影响倍率特性的主要原因，考虑到孔隙率和孔隙的曲折连同，离子在孔隙内的迁移距离比极片厚度多出很多倍。

（3）N/P比

负极活性物质克容量×负极面密度×负极活性物含量比÷（正极活性物质克容量×正极面密度×正极活性物含量比）

石墨负极类电池N/P要大于1.0，一般1.04~1.20，这主要是出于安全设计，主要为了防止负极析锂，设计时要考虑工序能力，如涂布偏差。但是，N/P过大时，电池不可逆容量损失，导致电池容量偏低，电池能量密度也会降低。

而对于钛酸锂负极，采用正极过量设计，电池容量由钛酸锂负极的容量确定。正极过量设计有利于提升电池的高温性能：高温气体主要来源于负极，在正极过量设计时，负极电位较低，更易于在钛酸锂表面形成SEI膜。

（4）涂层的压实密度及孔隙率

在生产过程中，电池极片的涂层压实密度计算公式：

而考虑到极片辊压时，金属箔材存在延展，辊压后涂层的面密度通过下式计算：

涂层由活物质相、碳胶相和孔隙组成，孔隙率计算公式：

其中，涂层的平均密度为：

（5）首效

首效=首次放电容量/首次充电容量

日常生产中，一般是先化成再进行分容，化成充入一部分电，分容补充电后再放电，故而：

首效=（化成充入容量+分容补充电容量）/分容第一次放电容量

（6）能量密度

体积能量密度(Wh/L)=电池容量(mAh)×3.6(V)/（厚度(cm)*宽度(cm)*长度(cm)）

质量能量密度(Wh/KG)=电池容量(mAh)×3.6(V)/电池重量

来源：江子才

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

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

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

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

计算材料的充放电平台

锂离子电池又称作“摇摇椅(rocking chair)”。这是因为在锂离子电池中，锂离子通过电解质在正负极间穿梭。充电时，锂离子从正极材料(如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.3 eV。

锂离子的传输

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

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

另外一种常用的计算方法是AIMD。由于其基于DFT，所以力场函数(Force filed)不用外加设定，这相比于经典的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.3 eV。在vasp中，考虑范德华修正的相关参数为IVDW。

材料稳定性的计算

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

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

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

结论

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

参考文献

[1]M. Armand and J.-M. Tarascon, nature, 2008, 451, 652.

[2]A. Urban, D.-H. Seo and G. Ceder, npj Computational Materials, 2016, 2.

[3] A. Jain, G. Hautier, S. P. Ong, C. J. Moore, C. C. Fischer, K. A. Persson and G. Ceder, Physical Review B, 2011, 84, 045115.

[4] N. Kamaya, K. Homma, Y. Yamakawa, M. Hirayama, R. Kanno, M. Yonemura, T. Kamiyama, Y. Kato, S. Hama and K. Kawamoto, Nature materials, 2011, 10, 682.

[5] G. A. Tritsaris, E. Kaxiras, S. Meng and E. Wang, Nano letters, 2013, 13, 2258-2263.

[6]X. He, Y. Zhu, A. Epstein and Y. Mo, npj Computational Materials, 2018, 4.

[7] L.-C. Yin, J. Liang, G.-M. Zhou, F. Li, R. Saito and H.-M. Cheng, Nano Energy, 2016, 25, 203-210.

[8]Q. Zhang, Y. Wang, Z. W. Seh, Z. Fu, R. Zhang and Y. Cui, Nano letters, 2015, 15, 3780-3786.

[9] S. Ping Ong, L. Wang, B. Kang and G. Ceder, Chemistry of Materials, 2008, 20, 1798-1807

文章来源：材料牛

锂电池的计算计算｜常用锂电参数计算公式

计算｜常用锂电参数计算公式

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