| 查看: 492 | 回复: 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
» 猜你喜欢
310求调剂
已经有14人回复
-
290调剂生物0860
已经有28人回复
-
山东双非院校考核超级无底线，领导幸灾乐祸，教师遭殃恐
已经有4人回复
-
279学硕食品专业求调剂院校
已经有14人回复
-
0831生医工第一轮调剂失败求助
已经有16人回复
-
电子信息270求调剂
已经有18人回复
-
食品与营养（0955）271求调剂
已经有5人回复
-
293求调剂
已经有11人回复
-
322求调剂
已经有13人回复
-
272分材料子求调剂
已经有50人回复
