| 查看: 821 | 回复: 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 |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有8人回复
-
26申博
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有22人回复
-
最失望的一年
已经有4人回复
-
国自然申请面上模板最新2026版出了吗？
已经有19人回复
-
请教限项目规定
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
推荐一本书
已经有13人回复
-
疑惑？
已经有5人回复
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
4楼2016-03-30 18:51:54
2楼2016-03-30 17:02:12
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
3楼2016-03-30 18:50:31
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.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