求助如何修正Buckingham势让其在分子与表面近距离时能体现排斥作用

Jack

最近在用第一性原理软件构造分子与金属表面的相互作用势能函数，如图所示发现Buckingham势能很好的拟合所得的数据点（也试过描述分子与表面的Morse势，9-3，10-4等势能函数，但效果不好）。

但是如这篇博客所说https://blog.sciencenet.cn/blog-3437453-1281254.html，"因为指数项收敛为常数r→0，而r ^ -6项发散，所以当r变小时，白金汉势变得有吸引力。当处理原子间距离非常短的结构时，这可能会成问题，因为任何穿过某个阈值的原子核都将在零距离处牢固地（非物理地）相互束缚。"，所以Buckingham势在距离表面很近时会有一个很强吸引作用，明显与现实的巨大排斥作用不符。

为了处理这个问题，也尝试过以下2种办法，都不太理想，所以请教各位大佬，有什么好的方法或者好的势能函数来拟合分子与金属表面的势能函数。
1. 采用cutoff在表面附近进行截断。对于有些分子，在排斥作用的最高点数值依然很小，进行截断的话势能排斥作用依然不强。后续想尝试计算更多离表面更近的点来拟合参数。
2. 按照文献https://doi.org/10.1080/00268976.2019.1669836，在Buckingham势后面加上一个24-6 LJ potential。拟合后发现，在表面附近依然存在很强的吸引势。

sobereva

设计一个切换函数，随距离切换Buckingham和LJ势的权重

Jack

sobereva 发表于 2022-7-6 07:20
设计一个切换函数，随距离切换Buckingham和LJ势的权重

好方法，谢谢Sob老师！之前一直想着和现存或者文献的势函数拟合，主要是担心审稿人会问自行设计的势函数的合理性。

sobereva

Jack 发表于 2022-7-6 07:43
好方法，谢谢Sob老师！之前一直想着和现存或者文献的势函数拟合，主要是担心审稿人会问自行设计的势函数 ...

能有什么不合理的，全程r范围拟合效果都理想，还能怎么质疑？不仅没什么可质疑的，还是个很有意义的创新。瞎质疑就直接怼回去

Jack

sobereva 发表于 2022-7-7 03:16
能有什么不合理的，全程r范围拟合效果都理想，还能怎么质疑？不仅没什么可质疑的，还是个很有意义的创新 ...

是的，谢谢Sob老师！

莫凡凡

大佬您好，我看您问到buckingham势的问题，虽然我不是很懂（LAMMPS萌新，最近在学）但感觉你好厉害，我最近也在学LAMMPS，要用到buckingham势或morse势，不知大佬可不可以留个联系方式，想具体请教您一下🙏🏻🙏🏻