24小时热门版块排行榜    

查看: 3282  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +16 Lanmanbaby 2026-08-09 26/1300 2026-08-11 17:16 by 天神眷顾
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 yPv6ndDtiRxQ 2026-08-11 3/150 2026-08-11 17:15 by CLKTSgjM0Uqj
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +6 布布和一二 2026-08-10 7/350 2026-08-11 16:57 by Ai19911123
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +7 Tide man 2026-08-10 8/400 2026-08-11 15:49 by 铁帽子农民
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
请填处理意见