Please wait a minute...
欢迎访问过程工程学报, 今天是
图/表 详细信息
短链多硫化物在离子液体中聚集行为的分子动力学模拟
胡天媛, 王艳磊, 霍锋, 何宏艳
过程工程学报    2021, 21 (7): 847-856.   DOI:10.12034/j.issn.1009-606X.220221
摘要   （391） HTML （9） PDF（pc） （1273KB）（180）  

离子液体因其优异的物化性质、能抑制多硫化物溶解等特点，近年来被广泛应用于锂硫电池电解液中。在电池充放电产物中，难溶性Li2S和Li2S2易聚集沉积在电极表面，影响电池性能，而目前关于其团聚行为与电解液性质的微观机理研究较少。本工作利用量化计算和分子动力学模拟分析了短链Li2S和Li2S2在离子液体中的微观结构以及形成团簇的情况。通过分析体系的微观结构发现，阳离子中主要与S作用的是侧链甲基，短链多硫化物之间Li-S作用远强于与阴离子的Li-O作用。团簇尺寸分布的结果表明，短链多硫化物在[TFSI]型离子液体中易形成多分子的大团簇，Li2S2体系比Li2S体系中的大团簇比例更高；离子液体阴离子配位能力越强，形成大的Li2S团簇比例越少，但阴离子的构型特点和作用形式也会对团簇的尺寸结构造成影响。


RILsDensity/(g/cm3)D+/×10-11 m2/sD-/×10-11 m2/s
[PP13][TFSI]1.381.843±0.0191.339±0.013
[PP13][OTf]1.251.177±0.0120.734±0.010
View table in article
表2 298 K下MD模拟得到的离子液体密度与自扩散系数
正文中引用本图/表的段落
离子液体因其优异的物化性质、能抑制多硫化物溶解等特点，近年来被广泛应用于锂硫电池电解液中。在电池充放电产物中，难溶性Li2S和Li2S2易聚集沉积在电极表面，影响电池性能，而目前关于其团聚行为与电解液性质的微观机理研究较少。本工作利用量化计算和分子动力学模拟分析了短链Li2S和Li2S2在离子液体中的微观结构以及形成团簇的情况。通过分析体系的微观结构发现，阳离子中主要与S作用的是侧链甲基，短链多硫化物之间Li-S作用远强于与阴离子的Li-O作用。团簇尺寸分布的结果表明，短链多硫化物在[TFSI]型离子液体中易形成多分子的大团簇，Li2S2体系比Li2S体系中的大团簇比例更高；离子液体阴离子配位能力越强，形成大的Li2S团簇比例越少，但阴离子的构型特点和作用形式也会对团簇的尺寸结构造成影响。
最近的研究表明短链Li 2S和Li 2S 2在高供体数(Donor Number, DN)或高介电常数的溶剂或锂盐阴离子体系中的溶解度会增加，部分以离子形式存在。Pan等[2]研究发现高介电常数的溶剂如二甲亚砜(DMSO)与高配位能力的阴离子如三氟磺酸根离子(Tf-)能够有效地提高Li 2S 2的溶解度，放电后的电极电镜照片并未观察到明显的短链多硫化物的沉积。该课题组还研究了溶剂供体数对Li 2S成核过程的影响，发现当使用碳纤维作为正极材料、高供体数溶剂为电解液时，Li 2S最终会呈现出“3D花”的形状[3]。Vijayakumar等[4]和Barchasz等[5]研究发现在DMSO溶剂中唯一稳定存在的多硫化物自由基只有S 3?―。Cuisinier等[6]受到启发测定了高介电常数以及高供体数的电子对给体溶剂中S 3?-的含量，结果表明这些溶剂中S 3?―的含量高于S 42-及S 62-等二价多硫离子的含量，他们认为S 3?-的存在促使Li 2S溶解增加，不会形成电极-Li 2S-电解液的三相界面。Zhang等[7]采用了一种高介电常数的质子型溶剂四甲基脲(TMU)作为锂硫电池电解液，发现其能够诱导S 3?-自由基的产生，在DOL/TMU溶剂体系中，S2-在已有的Li 2S晶核上3D沉积，形成薄而致密但又多孔的析出层。Chu等[8]研究了含有不同供体数锂盐阴离子的电解液体系，发现高供体数锂盐阴离子如Tf-和Br-体系可以有效提高活性物质S的利用率，并通过简单的分子动力学模拟解释了这一现象，认为这是由于高供体数的阴离子与Li+的结合能力强，使Li 2S呈现3D生长，而不是沉积形成致密的保护膜阻挡正极发生电化学反应。
对于离子液体的阴阳离子采用的是联合原子力场，将阳离子侧链上的CH 2, CH 3及[OTf]-和[TFSI]-上的CF 3基团按联合原子处理。[TFSI]-和[OTf]-的力场参数来源于Zhong等[19]建立的联合原子力场，其他三种阴离子采用的是OPLS力场[20]。阳离子中平衡键长、键角由优化的结构得到[B3LYP/6-311+g(d,p)水平]，环上各原子之间键结作用的力参数及LJ参数则参考了OPLS-AA力场中的胺类结构[21]，而侧链则参考了Zhong等[19]开发的咪唑类阳离子联合原子力场的参数，原子电荷由限制性拟合静电势(RESP)得到。基于上述的力场参数，对纯离子液体体系进行MD模拟验证，得到的密度与阴阳离子的自扩散系数见表2。[PP13][TFSI]模拟得到的密度与已有的实验数据的相对偏差小于3%，说明基于所建力场的分子动力学模拟能很好地重现该离子液体的实验密度值。根据与[PP13][TFSI]具有类似结构的[PP14][TFSI]实验中测得298 K时[PP14]+和[TFSI]-的自扩散系数[22]分别为0.98×10-11和0.85×10-11 m2/s，考虑到侧链长的阳离子均比侧链短的阳离子扩散慢，如实验得到1-乙基-3-甲基咪唑双三氟甲基磺酰亚胺比1-丁基-3-甲基咪唑双三氟甲基磺酰亚胺阳离子和阴离子自扩散系数分别高约80%和40%[23]，因此可认为模拟结果在合理范围内。
阴离子与短链多硫化物作用的Li-O RDFs以及短链多硫化物之间的Li-S RDFs如图4(c)和4(d)所示。[OTf]型离子液体中O-Li的峰高总是高于[TFSI]型离子液体，这与量化计算的相互作用能结果一致。所有体系中Li-S RDFs的第一个峰的高度均远超Li-O RDFs，说明在所研究的离子液体和短链多硫化物体系中，Li-S作用占主导，短链多硫化物容易发生自身聚集。对比图4(c)和4(d)可知，ILs-Li 2S 2体系中Li-O RDFs第一个峰普遍比ILs-Li 2S体系高，但出峰位置相近，这是因为Li 2S 2中两个Li原子的距离(DFT优化的结构显示为0.3377 nm)比Li 2S中两个Li原子之间的距离(DFT计算的最优结构显示为0.3927 nm)更近，造成电荷局域密度更高，表现出与O的静电作用更强。Li 2S体系中Li-S的RDFs在约0.25 nm出现非常尖锐的峰，峰高约22，而Li 2S 2体系中Li-S的RDFs在稍远处也出现峰值，峰高在7~10之间，说明Li 2S体系中Li-S作用更强，这从侧面解释了体系中Li-O RDFs峰更低的原因。对于Li 2S体系，不同离子液体对Li-S RDFs的改变很小，曲线几乎重合在一起。而Li 2S 2在[TFSI]型离子液体中Li-S RDFs比在[OTf]型离子液体中的Li-S RDFs稍高，这也是由阴离子-短链多硫化物、短链多硫化物之间这两种分子间作用的竞争关系所决定。
量化计算相互作用能通过Gaussian09软件[ 31 ]完成，所有离子或分子7684结构优化均采用密度泛函理论(Density Functional Theory, DFT)方法，基组方法为B3LYP/6-311+g(d,p)，并进行频率分析确保无虚频.相互作用能在M06-2X/6-311+g(d,p)水平上计算得到. ...
OPLS all-atom model for amines:? resolution of the amine hydration problem
1
1999
... 对于离子液体的阴阳离子采用的是联合原子力场，将阳离子侧链上的CH 2, CH 3及[OTf]-和[TFSI]-上的CF 3基团按联合原子处理.[TFSI]-和[OTf]-的力场参数来源于Zhong等[ 19 ]建立的联合原子力场，其他三种阴离子采用的是OPLS力场[ 20 ].阳离子中平衡键长、键角由优化的结构得到[B3LYP/6-311+g(d,p)水平]，环上各原子之间键结作用的力参数及LJ参数则参考了OPLS-AA力场中的胺类结构[ 21 ]，而侧链则参考了Zhong等[ 19 ]开发的咪唑类阳离子联合原子力场的参数，原子电荷由限制性拟合静电势(RESP)得到.基于上述的力场参数，对纯离子液体体系进行MD模拟验证，得到的密度与阴阳离子的自扩散系数见 表2 .[PP13][TFSI]模拟得到的密度与已有的实验数据的相对偏差小于3%，说明基于所建力场的分子动力学模拟能很好地重现该离子液体的实验密度值.根据与[PP13][TFSI]具有类似结构的[PP14][TFSI]实验中测得298 K时[PP14]+和[TFSI]-的自扩散系数[ 22 ]分别为0.98×10-11和0.85×10-11 m2/s，考虑到侧链长的阳离子均比侧链短的阳离子扩散慢，如实验得到1-乙基-3-甲基咪唑双三氟甲基磺酰亚胺比1-丁基-3-甲基咪唑双三氟甲基磺酰亚胺阳离子和阴离子自扩散系数分别高约80%和40%[ 23 ]，因此可认为模拟结果在合理范围内. ...
Synthesis and characterization of two ionic liquids with emphasis on their chemical stability towards metallic lithium
1
2007
... 对于离子液体的阴阳离子采用的是联合原子力场，将阳离子侧链上的CH 2, CH 3及[OTf]-和[TFSI]-上的CF 3基团按联合原子处理.[TFSI]-和[OTf]-的力场参数来源于Zhong等[ 19 ]建立的联合原子力场，其他三种阴离子采用的是OPLS力场[ 20 ].阳离子中平衡键长、键角由优化的结构得到[B3LYP/6-311+g(d,p)水平]，环上各原子之间键结作用的力参数及LJ参数则参考了OPLS-AA力场中的胺类结构[ 21 ]，而侧链则参考了Zhong等[ 19 ]开发的咪唑类阳离子联合原子力场的参数，原子电荷由限制性拟合静电势(RESP)得到.基于上述的力场参数，对纯离子液体体系进行MD模拟验证，得到的密度与阴阳离子的自扩散系数见 表2 .[PP13][TFSI]模拟得到的密度与已有的实验数据的相对偏差小于3%，说明基于所建力场的分子动力学模拟能很好地重现该离子液体的实验密度值.根据与[PP13][TFSI]具有类似结构的[PP14][TFSI]实验中测得298 K时[PP14]+和[TFSI]-的自扩散系数[ 22 ]分别为0.98×10-11和0.85×10-11 m2/s，考虑到侧链长的阳离子均比侧链短的阳离子扩散慢，如实验得到1-乙基-3-甲基咪唑双三氟甲基磺酰亚胺比1-丁基-3-甲基咪唑双三氟甲基磺酰亚胺阳离子和阴离子自扩散系数分别高约80%和40%[ 23 ]，因此可认为模拟结果在合理范围内. ...
Physicochemical properties and structures of room-temperature ionic liquids. 3. variation of cationic structures
1
2006
... 对于离子液体的阴阳离子采用的是联合原子力场，将阳离子侧链上的CH 2, CH 3及[OTf]-和[TFSI]-上的CF 3基团按联合原子处理.[TFSI]-和[OTf]-的力场参数来源于Zhong等[ 19 ]建立的联合原子力场，其他三种阴离子采用的是OPLS力场[ 20 ].阳离子中平衡键长、键角由优化的结构得到[B3LYP/6-311+g(d,p)水平]，环上各原子之间键结作用的力参数及LJ参数则参考了OPLS-AA力场中的胺类结构[ 21 ]，而侧链则参考了Zhong等[ 19 ]开发的咪唑类阳离子联合原子力场的参数，原子电荷由限制性拟合静电势(RESP)得到.基于上述的力场参数，对纯离子液体体系进行MD模拟验证，得到的密度与阴阳离子的自扩散系数见 表2 .[PP13][TFSI]模拟得到的密度与已有的实验数据的相对偏差小于3%，说明基于所建力场的分子动力学模拟能很好地重现该离子液体的实验密度值.根据与[PP13][TFSI]具有类似结构的[PP14][TFSI]实验中测得298 K时[PP14]+和[TFSI]-的自扩散系数[ 22 ]分别为0.98×10-11和0.85×10-11 m2/s，考虑到侧链长的阳离子均比侧链短的阳离子扩散慢，如实验得到1-乙基-3-甲基咪唑双三氟甲基磺酰亚胺比1-丁基-3-甲基咪唑双三氟甲基磺酰亚胺阳离子和阴离子自扩散系数分别高约80%和40%[ 23 ]，因此可认为模拟结果在合理范围内. ...
Densities and derived thermodynamic properties of imidazolium-, pyridinium-, pyrrolidinium-, and piperidinium-based ionic liquids
1
2008
... Note: the experimental density of [PP13][TFSI] is 1.41 g/cm3 for reference[ 24 ]. ...
Adsorption of insoluble polysulfides Li 2S x ( x=1, 2) on Li 2S surfaces
1
2015
... 对于短链多硫化物，采用全原子力场描述其精细结构.文献[ 25 ]中的晶体结构，Li 2S和Li 2S 2中Li-S平衡键长分别采用0.211和0.224 nm，其对应的键伸缩常数分别为8500和5200 kcal/(mol·nm2)，由量化计算键伸缩振动的频率[ 26 ]得来.平衡键角由优化的结构得到(B3LYP/TZVP水平上)，对应的键角弯曲常数同样由频率计算获得，具体的键角参数如 表3 所示.通过RESP方法得到Li 2S中S和Li的原子电荷分别为-1.4230e和0.7115e，而Li 2S 2中S和Li的原子电荷分别为-0.6139e和0.6139e.考虑到静电作用和范德华作用对性质影响较大，对文献中可得的Li和S原子的LJ参数进行了一系列模拟，各体系的LJ参数及模拟结果见 表4 .将模拟得到的S 8, Li 2S密度与实验值进行对比(S 8和Li 2S的实验密度值分别为2.03和1.66 g/cm3)，最终确定合适的 σ和 ε，即 表4 中System-5所采用的LJ参数. ...
Relating normal vibrational modes to local vibrational modes with the help of an adiabatic connection scheme
1
2012
... 对于短链多硫化物，采用全原子力场描述其精细结构.文献[ 25 ]中的晶体结构，Li 2S和Li 2S 2中Li-S平衡键长分别采用0.211和0.224 nm，其对应的键伸缩常数分别为8500和5200 kcal/(mol·nm2)，由量化计算键伸缩振动的频率[ 26 ]得来.平衡键角由优化的结构得到(B3LYP/TZVP水平上)，对应的键角弯曲常数同样由频率计算获得，具体的键角参数如 表3 所示.通过RESP方法得到Li 2S中S和Li的原子电荷分别为-1.4230e和0.7115e，而Li 2S 2中S和Li的原子电荷分别为-0.6139e和0.6139e.考虑到静电作用和范德华作用对性质影响较大，对文献中可得的Li和S原子的LJ参数进行了一系列模拟，各体系的LJ参数及模拟结果见 表4 .将模拟得到的S 8, Li 2S密度与实验值进行对比(S 8和Li 2S的实验密度值分别为2.03和1.66 g/cm3)，最终确定合适的 σ和 ε，即 表4 中System-5所采用的LJ参数. ...
A second generation force field for the simulation of proteins, nucleic acids, and organic molecules
2
1995
... Force field simulations for LJ parameters (compared with experimental densities of Li 2S and S 8)
本文的其它图/表