24小时热门版块排行榜    

查看: 3268  |  回复: 9

chufengj

新虫 (正式写手)

[求助] dump命令无法导出文件 已有1人参与

-------------------------------------------------------以下是疑惑---------------------------------------------------------------------------

我想做的是导出不同温度弛豫时,铜原子的势能,用compute计算然后用dump导出。

可是现在dump无法导出,但是它又没报错,所以在想会不会是导出到其他文件夹里面了。

还有,我用thermo输出热力学信息,里面的E_Pair、E_mol分别代表什么?

请各位指导,谢谢!

dump命令无法导出文件
in 1.png


dump命令无法导出文件-1
in2.png


dump命令无法导出文件-2
out1.png


dump命令无法导出文件-3
out2.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luftf

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
chufengj: 金币+50, ★★★★★最佳答案, 解决问题了 2015-11-23 22:27:43
lammps是从上到下逐行执行的,所以得注意命令的顺序
你的run在dump前面,已经run了,所以就不dump了
把dump那一行剪切贴到run的前面就行了
2楼2015-11-21 18:24:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chufengj

新虫 (正式写手)

引用回帖:
2楼: Originally posted by luftf at 2015-11-21 18:24:20
lammps是从上到下逐行执行的,所以得注意命令的顺序
你的run在dump前面,已经run了,所以就不dump了
把dump那一行剪切贴到run的前面就行了

你好,我还有两个问题想请教您。

用thermo输出的热力学信息里面的E_pair和E_mol是啥意思?

还有我用dump导出的*.atom 文件打不开。
3楼2015-11-23 22:31:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chufengj

新虫 (正式写手)

引用回帖:
2楼: Originally posted by luftf at 2015-11-21 18:24:20
lammps是从上到下逐行执行的,所以得注意命令的顺序
你的run在dump前面,已经run了,所以就不dump了
把dump那一行剪切贴到run的前面就行了

你好,我还有两个问题想请教您。

用thermo输出的热力学信息里面的E_pair和E_mol是啥意思?

还有我用dump导出的*.atom 文件打不开,一点击就是火狐浏览器打开了~
4楼2015-11-23 22:31:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chufengj

新虫 (正式写手)

引用回帖:
2楼: Originally posted by luftf at 2015-11-21 18:24:20
lammps是从上到下逐行执行的,所以得注意命令的顺序
你的run在dump前面,已经run了,所以就不dump了
把dump那一行剪切贴到run的前面就行了

我想得到的是compute计算出的势能。

因为我想做一个能量随时间或者温度的变化曲线。
5楼2015-11-23 22:37:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luftf

木虫 (正式写手)

引用回帖:
5楼: Originally posted by chufengj at 2015-11-23 22:37:45
我想得到的是compute计算出的势能。

因为我想做一个能量随时间或者温度的变化曲线。...

1,关于能量
etotal是系统总能量=系统势能+动能
E_pair是对势的能量。比如你只用LJ势函数,就统计的全部是对势
E_mol是分子能,包括键长、键角、二面角变化需要的能量。
如果要模拟碳管、石墨之类的或者碳氢化合物、烃类等,就可能需要统计键长、键角、二面角能量
而就你上传的算例而言,铜是金属,用的是EAM势函数,所以无所谓键长、键角、二面角能量

如果要做一个能量随时间或者温度的变化曲线,直接用etotal好了

2,关于输出文件,
打不开是啥意思?难道你是双击*.atom文件?如果是双击的话,肯定打不开。。。
如果你想要查看系统的形貌随时间的演化,去下载vmd或atomeye打开
如果你想要统计数据,下载个文本编辑器,例如ultraedit之类的


3,另外帮你检查了输入文件,其中可能的问题:
(1)timestep, 0.01太大了,系统容易崩溃,一般用0.001
(2)势函数,铜的势函数不要用系统自带的Cu_u3.eam,当然如果只是拿来练手就没事,如果用来发paper可能会被审稿人揪着不放的
(3)系统的温控并不好,3000步的时候389度,怀疑系统是否已经平衡,建议增加弛豫时间,或者调整nvt那一行的参数。
6楼2015-11-24 00:22:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kenybin

木虫 (小有名气)

认认真真做每一件事
7楼2015-11-24 07:09:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老夫子成

金虫 (正式写手)

8楼2015-11-24 07:09:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chufengj

新虫 (正式写手)

引用回帖:
6楼: Originally posted by luftf at 2015-11-24 00:22:09
1,关于能量
etotal是系统总能量=系统势能+动能
E_pair是对势的能量。比如你只用LJ势函数,就统计的全部是对势
E_mol是分子能,包括键长、键角、二面角变化需要的能量。
如果要模拟碳管、石墨之类的或者碳氢化 ...

你好,先谢谢你细致的回答。

还遇到两个问题不太明白:

1.我用atomeye打开,*.atom文件出现关于PATH的提示。

2.我给变量赋予不同的温度,然后在不同温度下弛豫,想得到不同温度下的能量,结果算出的能量全部相同,不知道是为什么
dump命令无法导出文件-4
atom文件打开失败.png


dump命令无法导出文件-5
compute计算结果.png

9楼2015-11-27 23:02:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chufengj

新虫 (正式写手)

引用回帖:
6楼: Originally posted by luftf at 2015-11-24 00:22:09
1,关于能量
etotal是系统总能量=系统势能+动能
E_pair是对势的能量。比如你只用LJ势函数,就统计的全部是对势
E_mol是分子能,包括键长、键角、二面角变化需要的能量。
如果要模拟碳管、石墨之类的或者碳氢化 ...

这是in文件
dump命令无法导出文件-6
in1.png


dump命令无法导出文件-7
in2.png

10楼2015-11-27 23:06:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chufengj 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 04:07 by 6vVgjDL4CnGu
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见