24小时热门版块排行榜    

查看: 2315  |  回复: 4

amynihao

金虫 (正式写手)

寂寞大虫

[求助] lammps的单位问题 已有1人参与

我做的是金属体系,我用的是LJ势函数,那么我在设置单位的时候,units应该设置为metal呢还是lj呢?
回复此楼
不是要做一个单纯优秀的人,而是要做一个不可替代的人。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

19861226ya

木虫 (著名写手)

姜死了,我想种蒜了!
2楼2018-08-02 12:29:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
amynihao: 金币+15, ★★★很有帮助, 非常感谢 2018-08-02 15:40:04
本帖仅楼主可见
3楼2018-08-02 15:07:55
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

Nater1ver

木虫 (正式写手)

你应该仔细阅读手册上的units章节,单位只是为了压缩数字存储的大小,不至于动不动-10次方-20次方几十位小数,在c++中会这用到长双精度类型甚至位数不够以至于运算容易出错。在我的印象中除了lj以外,就是metal和real单位最为实用,理论上,正如楼上所说的,只要单位统一,你可以用任何一种单位,但是做真实系统,转换成lj是很不方便的,常常用metal或者real。最重要的一点是,你用的势函数很有可能是带有单位的,你应该查看势函数每个参数的说明以及常用例子使用的单位,比如eam,可能有十几到二十多个参数,手册上会标明它们各自的默认单位,你如果对一个metal的势函数用了real单位,那里面的所有energy单位(因为metal和real最大的差别只是energy单位不同),软件是不会转换的,lammps会把一个kcal/mol的能量当成eV计算。所以在使用不同的单位之前,你要自己改势函数文件,把单位换算过来。不要问我怎么知道的

发自小木虫Android客户端
4楼2018-08-02 17:38:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Nater1ver

木虫 (正式写手)

引用回帖:
4楼: Originally posted by Nater1ver at 2018-08-02 17:38:17
你应该仔细阅读手册上的units章节,单位只是为了压缩数字存储的大小,不至于动不动-10次方-20次方几十位小数,在c++中会这用到长双精度类型甚至位数不够以至于运算容易出错。在我的印象中除了lj以外,就是metal和re ...

总而言之,以计算中用unit命令定义的单位为主,所有用到的参数包括势函数都要转换成定义的单位才行

发自小木虫Android客户端
5楼2018-08-02 17:42:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 amynihao 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
信息提示
请填处理意见