| 查看: 470 | 回复: 0 | ||
cgzhang_gg铜虫 (著名写手)
|
[求助]
谁能告诉我为什么第一次计算结果和后三次不同，明明是四次完全一样的计算
|
|
下面是我的in.TiO2 文件，第一次计算结果和后三次明显不同，但是我导出不同循环次的 t.1.0.cfg t.2.0.cfg t.3.0.cfg t.4.0.cfg 构型文件，发现都是一样的，而且通过pair_write 导出的不同循环次的pair 也是一样的，就是通过print 输出的初构型优化前的能量 不同。而且我还试了试，不同模拟原子个数的情况，有时四次的print 输出能量四次都相同，有时不同，但大部分是不同的。非常纠结，测试了很长时间了，真心找不出原因。望大神指导和交流。附件中是我用的TiO2.pair.table **************************************************************************************************************************************** variable a loop 4 units metal atom_style charge boundary p p p lattice custom 4.493 a1 1.0 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 0.669& basis 0.0 0.0 0.0 basis 0.5 0.5 0.5 & basis 0.303 0.303 0.0 basis 0.697 0.697 0.0 & basis 0.803 0.197 0.5 basis 0.197 0.803 0.5 region box block -3 3 -3 3 -3 3 create_box 2 box create_atoms 2 box basis 1 1 basis 2 1 & basis 3 2 basis 4 2 basis 5 2 basis 6 2 mass 1 47.867 # Ti mass 2 15.9994 # O set type 1 charge 2.196 # Ti set type 2 charge -1.098 # O # potential pair_style hybrid/overlay table linear 150000 coul/long 11.0 pair_coeff 1 1 table TiO2.pair.table Ti_Ti 11.0 pair_coeff 1 2 table TiO2.pair.table Ti_O 11.0 pair_coeff 2 2 table TiO2.pair.table O_O 11.0 pair_coeff * * coul/long kspace_style pppm 1.0e-5 run 0 dump 1 all cfg 1 t.$a.*.cfg id type xs ys zs q dump_modify 1 element Ti O variable E equal pe print "@@@@ a=$a before minimum! E=${E}" run 4 neighbor 2.0 bin neigh_modify delay 2 every 1 check yes next a clear jump in.TiO2 ******************************************************************************************************************************* 输出结果： @@@@ a=1 before minimum! E=-17183.269383955812373 @@@@ a=2 before minimum! E=-144792.3913198087248 @@@@ a=3 before minimum! E=-144792.3913198087248 @@@@ a=4 before minimum! E=-144792.3913198087248 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : TiO2.pair.table
2013-06-23 20:59:10, 21.89 M
» 猜你喜欢
职称评审没过，求安慰
已经有49人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
三无产品还有机会吗
已经有6人回复