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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 【讨论】势参数拟合软件
81/1返回列表
查看: 1111  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

楚辰

铜虫 (初入文坛)

[交流] 【讨论】势参数拟合软件

想跟大家讨论下，MD计算时大家都用什么势参数拟合软件？
本人用DLPOLY计算合金界面，用滴是GULP，有没有用别的滴同志？

本帖关键词：合金 GULP 势参数 软件 界面
回复此楼

» 猜你喜欢
1楼 2009-12-02 17:22:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

优秀版主
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zeoliters(金币+2,VIP+0):多谢回帖交流！ 12-8 15:06
有针对特定势形式的调势代码
如针对EAM势的 FORCEMATCH调试代码
改改就能适用到其他形式的势了
不过相对来讲公开的代码不多 留心就能发现
赞一下(11人) 回复此楼
2楼2009-12-02 17:46:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

楚辰

铜虫 (初入文坛)
引用回帖:
Originally posted by lei0736 at 2009-12-2 17:46:
有针对特定势形式的调势代码
如针对EAM势的 FORCEMATCH调试代码
改改就能适用到其他形式的势了
不过相对来讲公开的代码不多 留心就能发现

那针对FS势有没有特定的调试代码呢？这些代码一般在哪能够找到的
赞一下 回复此楼
3楼2009-12-03 10:24:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

优秀版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
应该也有 清华柳百新院士的综述 Physics Reports 455 (2008) 1–134
里面调势的程序就是我师兄用MATLAB写的
赞一下(1人) 回复此楼
4楼2009-12-03 10:34:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

doudouding

木虫 (小有名气)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zeoliters(金币+3,VIP+0):多谢回帖交流！ 12-8 15:06
调出来的势参数的可靠性由什么来验证？
GULP里面通常是用结构参数、介电常数、弹性常数这些物理量来拟合，多数情况下拟合的结果并不是非常好的结果。文献里对同样的体系会给出不同的势参数，都可以和实验值符合的很好，到底那个是更合理的呢？总不能先预设一个我需要的研究结果再看哪个更符合我的预期这样的方法来选取吧。
另外还有一个问题，就是掺杂材料中，掺杂原子和周围原子之间的势参数应该由什么确定，比如example里的MgO中掺杂Li原子，Li-O之间的势参数又该由什么确定呢？
赞一下(11人) 回复此楼
5楼2009-12-06 19:35:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

优秀版主
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zeoliters(金币+5,VIP+0):多谢回帖交流！ 12-8 15:06
调出来的势参数的可靠性由什么来验证？
GULP里面通常是用结构参数、介电常数、弹性常数这些物理量来拟合，多数情况下拟合的结果并不是非常好的结果。
----------------------------------------------
一般也就这些了 如果想拟合一个适用范围更广的势 那就不仅要拟合块体性质 还有表面能 表面驰豫 结构能量差 熔点 缺陷形成能 等等 拟合的结果不够好可能是势形式的限制 也许你的体系并不适合用这个势来描述 牛头不对马嘴是很自然的事。比如你用LJ势去描述金属，不论怎么调整柯西关系总是不正确的。还有你说的非常好是什么程度？既然是经验公式，想精确到小数点后若干位当然是不可能的


文献里对同样的体系会给出不同的势参数，都可以和实验值符合的很好，到底那个是更合理的呢？总不能先预设一个我需要的研究结果再看哪个更符合我的预期这样的方法来选取吧。
---------------------------------------------------------
这个可以接受啊 比如对金属就有EAM，MEAM, FS, TB, TBSMA, SC等等势形式。当然其中大部分本质上是一样的。没有最合理，只有更合理。基本原则就是尽可能多的物理性质符合实验或第一原理值。

另外还有一个问题，就是掺杂材料中，掺杂原子和周围原子之间的势参数应该由什么确定，比如example里的MgO中掺杂Li原子，Li-O之间的势参数又该由什么确定呢？
------------------------------------------
Li-O之间的势应该属于离子势吧 研究不多 BM可以做吧 BOP势也对氧化物有较好的描述。
赞一下(11人) 回复此楼
6楼2009-12-07 10:01:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

doudouding

木虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
GULP里表面能怎么计算？输出文件里只有lattice energy
赞一下(1人) 回复此楼
7楼2009-12-15 21:27:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+2,VIP+0):谢谢 12-16 12:37
表面能的计算不是默认的计算步骤，需要你用Keyword去设置它。详见说明书。
赞一下(11人) 回复此楼
8楼2009-12-16 09:07:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 楚辰 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭