| 查看: 1017 | 回复: 5 | ||
zhangguangping木虫 (著名写手)
|
[求助]
用收敛过的WAVECAR初始频率计算不收敛 已有1人参与
|
|
请教大家一个问题：我用10个核跑了一个电偶极矩矫正的slab体系的电子结构，收敛了，得到了WAVECAR，然后用它来初始了相同体系的做频率的计算（也做了偶极矫正），计算的时候用的是8个核，为什么不收敛了呢？ 预期的情况是用收敛过的WAVECAR初始频率计算，会很快收敛开始的电子自洽，然后进入频率计算。现在是开始的电子自洽都不收敛了。不理解 附： 1.偶极矫正的INCAR SYSTEM = M2-Au-E0 IALGO = 48 NELM = 500 EDIFF = 1E-4 EDIFFG = -0.02 NSW = 0 ISMEAR = 1 SIGMA = 0.1 PREC = Accurate ENCUT = 400 IBRION = 2 ISPIN = 1 LCHARG = .FALSE. LREAL = Auto LWAVE = .TRUE. LDIPOL = .TRUE. IDIPOL = 3 EFIELD = 0 2.偶极矫正的频率计算INCAR SYSTEM = M2-Au-freq IALGO = 48 NELM = 500 EDIFF = 1E-4 EDIFFG = -0.02 NSW = 2 ISMEAR = 1 SIGMA = 0.1 PREC = Accurate ENCUT = 400 IBRION = 5 POTIM = 0.01 NFREE = 2 ISPIN = 1 LCHARG = .FALSE. LREAL = Auto LWAVE = .FALSE. LDIPOL = .TRUE. IDIPOL = 3 NWRITE = 3 |
回复此楼» 猜你喜欢
这样的filecode谁见过
已经有13人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
2楼2016-03-30 17:02:12
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
3楼2016-03-30 18:50:31
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
4楼2016-03-30 18:51:54
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
|
现在发现问题的原因了，但是不明白为什么。 原因是采用偶极矫正的WAVECAR再次初始化相同的偶极矫正任务的时候，收敛的WAVECAR不是一个好的初始猜测。这个时候需要用没有偶极矫正的WAVECAR才行。 我用没有偶极矫正收敛的WAVECAR初始偶极矫正的计算或者频率计算（做了偶极矫正） 这个时候很快收敛： entering main loop N E dE d eps ncg rms rms(c) RMM: 1 -0.617285334113E+03 -0.61729E+03 -0.28948E-05 958 0.329E-03 0.356E-01 RMM: 2 -0.633045700125E+03 -0.15760E+02 -0.22436E-01 2091 0.263E-01 0.131E+01 RMM: 3 -0.617613909416E+03 0.15432E+02 -0.14867E-01 2087 0.233E-01 0.325E+00 RMM: 4 -0.617282326516E+03 0.33158E+00 -0.10908E-02 2085 0.729E-02 0.156E-01 RMM: 5 -0.617282748456E+03 -0.42194E-03 -0.12005E-03 1997 0.282E-02 0.204E-01 RMM: 6 -0.617281688730E+03 0.10597E-02 -0.10445E-04 1208 0.850E-03 0.374E-02 RMM: 7 -0.617281678568E+03 0.10162E-04 -0.16863E-05 1042 0.361E-03 1 F= -.61728168E+03 E0= -.61723961E+03 d E =-.126213E+00 Finite differences POTIM= 0.01000 DOF= 78 bond charge predicted |
5楼2016-03-31 06:42:21
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
6楼2016-03-31 06:43:13
20