| 查看: 486 | 回复: 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
» 猜你喜欢
323求调剂
已经有6人回复
-
一志愿北京化工大学 070300 学硕 336分 求调剂
已经有4人回复
-
352求调剂
已经有3人回复
-
一志愿东华大学化学070300，求调剂
已经有8人回复
-
277材料科学与工程080500求调剂
已经有7人回复
-
317求调剂
已经有18人回复
-
293求调剂
已经有5人回复
-
280分求调剂 一志愿085802
已经有7人回复
-
0854电子信息求调剂
已经有3人回复
-
263求调剂
已经有4人回复
