| 查看: 473 | 回复: 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
» 猜你喜欢
溴的反应液脱色
已经有7人回复
-
国自然申请面上模板最新2026版出了吗？
已经有11人回复
-
推荐一本书
已经有12人回复
-
基金申报
已经有4人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有4人回复
-
纳米粒子粒径的测量
已经有7人回复
-
常年博士招收(双一流，工科)
已经有4人回复
-
参与限项
已经有5人回复
-
有没有人能给点建议
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复