| 6 | 1/1 | 返回列表 |
| 查看: 868 | 回復(fù): 5 | ||
zhangguangping木蟲(chóng) (著名寫(xiě)手)
|
[求助]
用收斂過(guò)的WAVECAR初始頻率計(jì)算不收斂 已有1人參與
|
|
請(qǐng)教大家一個(gè)問(wèn)題:我用10個(gè)核跑了一個(gè)電偶極矩矯正的slab體系的電子結(jié)構(gòu),收斂了,得到了WAVECAR,然后用它來(lái)初始了相同體系的做頻率的計(jì)算(也做了偶極矯正),計(jì)算的時(shí)候用的是8個(gè)核,為什么不收斂了呢? 預(yù)期的情況是用收斂過(guò)的WAVECAR初始頻率計(jì)算,會(huì)很快收斂開(kāi)始的電子自洽,然后進(jìn)入頻率計(jì)算,F(xiàn)在是開(kāi)始的電子自洽都不收斂了。不理解 附: 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.偶極矯正的頻率計(jì)算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 |


木蟲(chóng) (著名寫(xiě)手)

木蟲(chóng) (著名寫(xiě)手)

木蟲(chóng) (著名寫(xiě)手)
|
現(xiàn)在發(fā)現(xiàn)問(wèn)題的原因了,但是不明白為什么。 原因是采用偶極矯正的WAVECAR再次初始化相同的偶極矯正任務(wù)的時(shí)候,收斂的WAVECAR不是一個(gè)好的初始猜測(cè)。這個(gè)時(shí)候需要用沒(méi)有偶極矯正的WAVECAR才行。 我用沒(méi)有偶極矯正收斂的WAVECAR初始偶極矯正的計(jì)算或者頻率計(jì)算(做了偶極矯正) 這個(gè)時(shí)候很快收斂: 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 |

木蟲(chóng) (著名寫(xiě)手)

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|