如何计算离子在聚合物结构上的迁移能垒

Dreamer1

老师好，我这里有一个阴离子聚合物结构，我们预测了Li离子在这个聚合物链段上的吸附位点和可能的迁移方式（如下图所示）。现在我想计算锂离子从其中一个位点迁移到另一个位点的迁移能垒，我的问题是，这个计算该采用gaussian还是cp2k去计算。之前我参加的第4届 CP2K第一性原理计算培训班里“过渡态搜索和反应路径的产生”部分有一些关于离子迁移的计算实例，不知道是否可以参考那部分实例进行计算？

sobereva

Gaussian和CP2K都完全可以
Gaussian的话通常需要你去猜迁移路径上各个过渡态（可能不止一个）结构然后优化过渡态。CP2K用NEB的话，可以自己给出迁移前后的结构并插点，如果中间有中间体、多个过渡态，在NEB能量变化变化曲线上会有所体现（之后再对每个基元过程用CI-NEB或dimer找轨道），因此适合不好估计中间有哪些过渡态/中间体的情况。也可以二者结合，比如CP2K做NEB确定过渡态粗略位置，作为Gaussian的初猜结构用opt=TS优化过渡态。