| 查看: 3013 | 回复: 9 | ||
[求助]
dump命令无法导出文件 已有1人参与
|
|
-------------------------------------------------------以下是疑惑--------------------------------------------------------------------------- 我想做的是导出不同温度弛豫时，铜原子的势能，用compute计算然后用dump导出。 可是现在dump无法导出，但是它又没报错，所以在想会不会是导出到其他文件夹里面了。 还有，我用thermo输出热力学信息，里面的E_Pair、E_mol分别代表什么？ 请各位指导，谢谢！  in 1.png  in2.png  out1.png  out2.png |
回复此楼» 猜你喜欢
一志愿东华大学化学070300，求调剂
已经有8人回复
-
277材料科学与工程080500求调剂
已经有7人回复
-
317求调剂
已经有18人回复
-
293求调剂
已经有5人回复
-
280分求调剂 一志愿085802
已经有7人回复
-
0854电子信息求调剂
已经有3人回复
-
263求调剂
已经有4人回复
-
石河子大学（211、双一流）硕博研究生长期招生公告
已经有3人回复
-
一志愿西安交通大学材料工程专业 282分求调剂
已经有13人回复
-
北科281学硕材料求调剂
已经有6人回复
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
10