| 查看: 543 | 回复: 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
» 猜你喜欢
麻烦专家们看看评委们的意见（F口面上）
已经有5人回复
-
2026年叶企孙基金
已经有6人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然面上复盘~欢迎讨论
已经有13人回复
-
基金系统什么内容也没有
已经有9人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有11人回复
-
怎么查啊
已经有6人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
10