| 查看: 2944 | 回复: 9 | ||
[求助]
dump命令无法导出文件 已有1人参与
|
luftf
木虫 (正式写手)
- 应助: 51 (初中生)
- 金币: 3257.7
- 散金: 685
- 红花: 19
- 帖子: 541
- 在线: 200.2小时
- 虫号: 821742
- 注册: 2009-08-04
- 性别: GG
- 专业: 金属材料的力学行为
2楼2015-11-21 18:24:20
3楼2015-11-23 22:31:20
4楼2015-11-23 22:31:43
5楼2015-11-23 22:37:45
luftf
木虫 (正式写手)
- 应助: 51 (初中生)
- 金币: 3257.7
- 散金: 685
- 红花: 19
- 帖子: 541
- 在线: 200.2小时
- 虫号: 821742
- 注册: 2009-08-04
- 性别: GG
- 专业: 金属材料的力学行为
|
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
kenybin
木虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 2547.8
- 红花: 7
- 帖子: 105
- 在线: 32.5小时
- 虫号: 3886912
- 注册: 2015-05-23
- 性别: GG
- 专业: 基础物理学

7楼2015-11-24 07:09:15
老夫子成
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2138.4
- 散金: 72
- 红花: 5
- 帖子: 857
- 在线: 96.7小时
- 虫号: 3477523
- 注册: 2014-10-15
- 专业: 零件加工制造
8楼2015-11-24 07:09:38
9楼2015-11-27 23:02:40
10楼2015-11-27 23:06:27











回复此楼