计算化学公社

标题: 求助，没有势函数参数的原子怎么找 [打印本页]

作者
Author: 无日白昼 时间: 2024-4-11 15:06
标题: 求助，没有势函数参数的原子怎么找
【求助】想问一下，例如Pb、Sn这类原子的势函数参数要去哪里找呢，用过了sobtop，GAff、amber力场文件里面都没有，uff力场的话精度太低了，这种情况要怎么办呢？

作者
Author: 小嘴抹了蜜 时间: 2024-4-11 17:09
我是在openkim的官网里找的势函数，你可以去搜一下看看，可以根据你选择的原子来给你显示相关的势函数

作者
Author: sobereva 时间: 2024-4-12 05:36
看具体这些元素出现在什么体系里、关注的具体是什么参数（LJ还是成键相关的项），没法一概而论
诸如只有一个这种元素出现在某个刚性区域里，而且离分子表面还挺远，那么这个原子的LJ参数准不准就无所谓，跟其它分子作用距离肯定很远，而且此时相关的成键项参数也无所谓，用谐振势维持成刚性就行了

作者
Author: 无日白昼 时间: 2024-4-12 09:39

sobereva 发表于 2024-4-12 05:36
看具体这些元素出现在什么体系里、关注的具体是什么参数（LJ还是成键相关的项），没法一概而论
诸如只有一 ...

你好老师，我想关注的是PbI2、SnCl2的结晶情况，这种情况应该是需要LJ参数准确度，想问问这种情况去哪里找合适的势函数参数呢？

作者
Author: 无日白昼 时间: 2024-4-12 09:39

小嘴抹了蜜发表于 2024-4-11 17:09
我是在openkim的官网里找的势函数，你可以去搜一下看看，可以根据你选择的原子来给你显示相关的势函数

好的感谢，我去看看！

作者
Author: 无日白昼 时间: 2024-4-12 09:44

无日白昼发表于 2024-4-12 09:39
好的感谢，我去看看！

我去找了一下，相关的只有EAM，不太适合我需要的LJ势函数参数，请问还有什么别的来源可以找到LJ参数吗？

作者
Author: 小嘴抹了蜜 时间: 2024-4-12 21:34

无日白昼发表于 2024-4-12 09:44
我去找了一下，相关的只有EAM，不太适合我需要的LJ势函数参数，请问还有什么别的来源可以找到LJ参数吗？

你可以到网上找一下，好像有一个计算lj参数的MATLAB脚本，可以直接算出来lj的参数

作者
Author: 小嘴抹了蜜 时间: 2024-4-12 21:45

无日白昼发表于 2024-4-12 09:44
我去找了一下，相关的只有EAM，不太适合我需要的LJ势函数参数，请问还有什么别的来源可以找到LJ参数吗？

或者有个lammpsstube里面的potential你去找一下看看

作者
Author: sobereva 时间: 2024-4-13 06:01

无日白昼发表于 2024-4-12 09:44
我去找了一下，相关的只有EAM，不太适合我需要的LJ势函数参数，请问还有什么别的来源可以找到LJ参数吗？

google搜相关体系模拟的文献看里面的参数，搜不到就自己根据其实验数据（如密度等）进行拟合或参数优化，可以参考主流的水模型原文了解小分子参数化的思想和方式

作者
Author: Weldingspock 时间: 2024-4-13 09:47
这是要做钙钛矿的模拟吗

作者
Author: 无日白昼 时间: 2024-4-15 09:27

Weldingspock 发表于 2024-4-13 09:47
这是要做钙钛矿的模拟吗

是的！请问您有相关的参数吗？

作者
Author: Weldingspock 时间: 2024-4-15 09:53

无日白昼发表于 2024-4-15 09:27
是的！请问您有相关的参数吗？

经验力场我只听说过一个叫MYP的，剩下的都是自己搞的机器学习势。重原子的势函数太少了

作者
Author: 无日白昼 时间: 2024-4-15 15:45

Weldingspock 发表于 2024-4-15 09:53
经验力场我只听说过一个叫MYP的，剩下的都是自己搞的机器学习势。重原子的势函数太少了

MYP我在查文献的时候看过，但是精度不太够。机器学习势我没有了解过，我去看看，谢谢你！

作者
Author: 小嘴抹了蜜 时间: 2024-4-17 17:05
你好，我找到一个matlab脚本获取LJ参数的链接了，你可以试试 https://zhuanlan.zhihu.com/p/116913952，（ps：我不知道能不能发链接，如果不可以，希望管理员见谅。）

作者
Author: 无日白昼 时间: 2024-4-22 14:54

小嘴抹了蜜发表于 2024-4-17 17:05
你好，我找到一个matlab脚本获取LJ参数的链接了，你可以试试 https://zhuanlan.zhihu.com/p/116913952，（ ...

你好，这个之前看到过，但是这个脚本原理就是从已有的UFF力场里面找出自己需要的原子LJ参数，UFF力场精度有点低，可能不太适用。很感谢你的帮助！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)