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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 关于用lammps的compute 命令计算原子位移
51/1返回列表
查看: 4980  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

我爱bigbang

金虫 (小有名气)

[求助] 关于用lammps的compute 命令计算原子位移 已有2人参与

我想要计算原子的位移，使用的是compute dis all displace/atom，但是计算出现错误。ERROR: Dump custom compute does not calculate per-atom vector
部分输入文件如下：
fix 1 all npt temp 300.0 300.1 1.0 x 0.0 0.0  3000 y 0.0 0.0 3000  z 0.0 0.0 3000
fix 2 all temp/rescale 1 300 300.1 0.05 1.0

compute         s all  stress/atom NULL
compute         strea all  reduce sum c_s[1]
variable        stressa equal c_strea/vol
compute         streb all  reduce sum c_s[2]
variable        stressb equal c_streb/vol
compute         strec all  reduce sum c_s[3]
variable        stressc equal c_strec/vol
compute         epot all  pe/atom
compute         csym all centro/atom 4
compute         coord all coord/atom 2.0
compute 1 all pe/atom
compute 2 all ke/atom
compute dis all displace/atom
thermo_style    custom step temp etotal press v_stressa  v_stressb v_stressc vol lx ly lz
dump 2 all custom 1000 dump.2 id type xu yu zu c_1 c_2 c_dis                                         

run  200000
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-05-03 15:25:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

引用回帖:
7楼: Originally posted by 我爱bigbang at 2016-05-06 08:55:08
有没有命令是输出矢量的呢？我用dis dis dis试了一下，但是用ovito看的时候是如下图这样的。但是我做的是压缩晶体，我感觉不应该出现这种情况吧。...

LAMMPS没有那么智能，自己把dump的数据处理一下，用tecplot显示矢量
赞一下 回复此楼
9楼2016-05-06 10:00:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

abinitio

版主 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
没用过此命令。
我想说的是，为什么不通过两个时间的位置坐标来自己 算呢？
赞一下 回复此楼
2楼2016-05-03 15:50:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我爱bigbang

金虫 (小有名气)
引用回帖:
2楼: Originally posted by abinitio at 2016-05-03 15:50:51
没用过此命令。
我想说的是，为什么不通过两个时间的位置坐标来自己 算呢？

我想dump出来的话，可以直接从OVITO中看到。
赞一下 回复此楼
3楼2016-05-03 15:54:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qw8276617

银虫 (小有名气)
你的thermo-style里含有矢量，而dump里不能含有矢量

发自小木虫IOS客户端
赞一下(1人) 回复此楼
4楼2016-05-03 20:55:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 352求调剂 +3 大米饭！ 2026-03-22 3/150 2026-03-22 23:28 by king123！
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 070300化学求调剂 +3 苑豆豆 2026-03-20 3/150 2026-03-22 14:48 by ColorlessPI
[考研] 291 求调剂 +3 化工2026届毕业� 2026-03-21 3/150 2026-03-22 14:26 by ColorlessPI
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[考研] 一志愿华中科技大学071000，求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 261求B区调剂，科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 0856调剂，是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
信息提示
关闭
请填处理意见
关闭