求助：使用lammps计算TaC熔点时Ta-C的势函数如何选择及拟合

圆又圆

各位老师好，我最近想利用lammps计算TaC的熔点，在关于Ta-C的势函数选择上查阅了相关文献，但找不到现有的有关Ta-C的势函数参数

1.     查阅文献中发现有关过渡族金属碳化物的研究使用L-J势、MEAM-2nn、以及tersoff势的都有，请问对于这种碳化物选择什么样的势函数更合适？

2.     现阶段我先自行对Ta-C的L-J势进行了拟合，参考教程为【[6.10]--势函数-Cu-C势曲线的DFT计算和LJ势拟合170907】 https://www.bilibili.com/video/BV1YT4y1P7o7/?share_source=copy_web&vd_source=5bcc7c367cef37d2e1fee7fb02eed56c，（如果不能发链接我会删除）利用MS对Cu-C的L-J势进行拟合。我根据教程在晶胞中分别建立了不同距离的Ta-C模型（图1），通过CASTEP计算了模型能量，利用公式(E(Ta-C)-E(Ta)-E(C))/2计算了形成能。但是与教程不同的是我的结构中不同原子间距其能量变化基本没有差别（图2），在MATLAB中的拟合结果也完全不对，是不是这个体系不能这样计算L-J势？

3.     后来我查阅了一些Sober老师对其他有关势函数拟合求助贴的解答，可以用Gaussian做势能面扫描拟合，但大多数大家的求助都是非晶体系，我想求助一下有没有对这方面比较了解的老师及同学，如果我算碳化钽这种晶体的Ta-C键也可以使用这种方法吗？是建立如上图MS中的Ta-C键就可以吗？

4.     最后我查阅的一些文献中有关ZrC，TiN等过渡族金属碳氮化物的是势函数拟合使用的是MEAM-2nn，拟合过程中涉及了使用ATK软件对于ZrCx晶格常数、形成焓进行了计算，如果我想进行Ta-C的MEAM-2nn势函数的拟合，晶格常数和形成焓能否用MS的CASTEP进行计算？

由于我是分子动力学的初学者，组里也没有懂得师兄师姐，所以问题比较多，以上问题我已经自行查阅了一些相关的文献、视频以及网络搜索，但仍然存在一些疑问，想要求助一下各位老师，希望老师们能够给我一些帮助或者建议，不胜感激。