24小时热门版块排行榜    

查看: 2338  |  回复: 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.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 07:01 by YXM1FDCSE8Bn
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 07:01 by YXM1FDCSE8Bn
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 06:53 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:43 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 04:04 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 25/1250 2026-11-01 03:30 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[基金申请] 论文的挂名怎么评判贡献?基金申请中,排名第一的挂名文章 +3 biily 2026-10-25 3/150 2026-10-30 18:33 by 旁观2020
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见