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团簇比例越少，但阴离子的构型特点和作用形式也会对团簇的尺寸结构造成影响。


SimulationσLi/×10-1 nmσS/×10-1 nmεLi/(kcal/mol)εS/(kcal/mol)ρLi2S/(g/cm3)ρS8/(g/cm3)
System-1[20]2.87003.55000.00050.25001.831.71
System-2[27]2.02613.56400.01830.25041.911.70
System-3[28]2.13003.55000.07650.25001.731.71
System-4[20,29]3.60003.30320.00050.25001.711.90
System-5[20,29,30]3.60003.39000.00050.40701.631.98
System-6[27,29,30]3.60003.39000.01830.40701.111.98
View table in article
表4 LJ参数的模拟(与Li 2S和S 8的实验密度对比)
正文中引用本图/表的段落
对于短链多硫化物，采用全原子力场描述其精细结构。文献[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参数。
量化计算相互作用能通过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)
本文的其它图/表