应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » ME/Gulp » Gulp的势函数的拟合
171/2返回列表12下一页
查看: 4330  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

流浪猫小琪

银虫 (正式写手)

[求助] Gulp的势函数的拟合 已有1人参与

最近在做分子力场的拟合，主要是Buckingham势函数的拟合，用于熔体计算，希望有经验的前辈可以帮助解答一下几点疑问。

1、关于gulp拟合，是不是应该首先计算熔体中离子所带的电荷，可是castep计算的Mulliken电荷在不同结构中波动较大，那么应该以哪种电荷计算的电荷来描述不同结构的电荷比较合理呢？

2、MS例子中给出的例子中只有使用Bulk Modulus和Elastic Constants和fit Structure三个，那么要得到一个较为靠谱的势函数参数，这几个足够吗？

3、Gulp中observe中的Fit total energy应该怎么输入呢，这个是不是结构的相互作用能，假设需要Fit total energy，那么是不是需要一个其他结构，以其与Fit structure的结构的能量差作为输入数据呢？

4、由于拟合的主要目的是用于高温熔体的模拟，这种方式拟合的参数能不能直接使用？

还有两个金币，加上今天刚领了5个，希望有前辈可以来解答
回复此楼

» 收录本帖的淘帖专辑推荐

动力学模拟

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2017-09-16 23:41:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
4楼: Originally posted by 流浪猫小琪 at 2017-09-18 00:13:39
我现在没有拟合电荷，电荷直接采用的Bader电荷。前辈可否告知一下如何利用介电常数来拟合core-shell电荷的大致步骤？...

通过密度泛函微扰论(DFPT)的方法计算可以直接计算出介电常数的。Bader电荷给出的是价电荷。

在core-shell里，shell的电荷不一定是价电荷。 shell上的电荷，导致体系的极化，需要和实验测（或计算出来的）的极化相关物理量，比如介电常数进行比较。
赞一下(1人) 回复此楼
8楼2017-09-18 20:09:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
你要加入适当的关键词：split .

发自小木虫Android客户端
赞一下(1人) 回复此楼
16楼2017-10-17 11:38:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

valenhou001

至尊木虫 (职业作家)
1.要拟合电荷的话，需基于core-shell模型，拟合介电常数或压电常数来得到core, shell的电荷，这个最不好拟合。如果只是Buckingham 势函数里的几个参数，体弹性模量，弹性常数，结构参数的利用是就可以的。

发自小木虫Android客户端
赞一下 回复此楼
2楼2017-09-17 08:07:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
本帖仅楼主可见
赞一下
3楼2017-09-17 14:33:56
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

流浪猫小琪

银虫 (正式写手)
引用回帖:
2楼: Originally posted by valenhou001 at 2017-09-17 08:07:52
1.要拟合电荷的话，需基于core-shell模型，拟合介电常数或压电常数来得到core, shell的电荷，这个最不好拟合。如果只是Buckingham 势函数里的几个参数，体弹性模量，弹性常数，结构参数的利用是就可以的。
...

我现在没有拟合电荷，电荷直接采用的Bader电荷。前辈可否告知一下如何利用介电常数来拟合core-shell电荷的大致步骤？
赞一下 回复此楼
4楼2017-09-17 23:13:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

流浪猫小琪

银虫 (正式写手)
引用回帖:
3楼: Originally posted by xmc057711 at 2017-09-17 14:33:56
低温下拟合的对于低温肯定可以，对于高温液体你需要验证能不能描述高温的关键性质

请问前辈，应该做怎么样的验证来说明适合？按照我的离解，高温之下熔体的结构有所变化，但是离子之间的距离变化不太大，低温得到的势函数应该还是可以描述这个距离的原子之间的相互作用的。
赞一下 回复此楼
5楼2017-09-17 23:18:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
本帖仅楼主可见
赞一下
6楼2017-09-18 07:49:01
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

匿名

用户注销 (正式写手)
本帖仅楼主可见
赞一下
7楼2017-09-18 07:52:17
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

引用回帖:
5楼: Originally posted by 流浪猫小琪 at 2017-09-18 00:18:37
请问前辈，应该做怎么样的验证来说明适合？按照我的离解，高温之下熔体的结构有所变化，但是离子之间的距离变化不太大，低温得到的势函数应该还是可以描述这个距离的原子之间的相互作用的。...

我也是比较怀疑“buckingham potential"模拟高温熔体性质的可靠性。因为这种势函数实在是太简单了。在高温熔体下，结果与0k的是结构应该是变化很大的。
赞一下 回复此楼
9楼2017-09-18 20:13:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

流浪猫小琪

银虫 (正式写手)
引用回帖:
8楼: Originally posted by valenhou001 at 2017-09-18 20:09:50
通过密度泛函微扰论(DFPT)的方法计算可以直接计算出介电常数的。Bader电荷给出的是价电荷。

在core-shell里，shell的电荷不一定是价电荷。 shell上的电荷，导致体系的极化，需要和实验测（或计算出来的）的极化 ...

对于极化电荷我了解一些。我是想请教一下：假设我们已经有了介电常数，应该怎么样在gulp里面拟合，我还没有做过这个。另外，这个介电常数是不是静介电常数，castep计算出来的那个？
赞一下 回复此楼
10楼2017-09-20 20:47:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 流浪猫小琪 的主题更新
171/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭