24小时热门版块排行榜    

查看: 3310  |  回复: 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

木虫 (正式写手)

引用回帖:
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的回帖
查看全部 10 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +11 医学老男孩 2026-08-13 29/1450 2026-08-14 00:04 by 大快人心啦
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 28/1400 2026-08-13 17:51 by 医学老男孩
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于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
信息提示
请填处理意见