| 查看: 1441 | 回复: 1 | |||
| 【悬赏金币】回答本帖问题，作者大力2010将赠送您 10 个金币 | |||
[求助]
请教高手，我在用vasp进行结构优化. 优化1周多了，还没结果，怎样解决？ 已有1人参与
|
|||
|
请教高手，我在用vasp进行结构优化. 优化1周多了，还没结果，怎样解决？ 我用的18核进行。 体系一共132个原子，我固定了中部的一些原子（X Y Z三个方向固定）， INCAR 和POSCAR如下所示： Global Parameters ISTART = 1 (Read existing wavefunction; if there) # ISPIN = 2 (Spin polarised DFT) # ICHARG = 11 (Non-self-consistent: GGA/LDA band structures) LREAL = .FALSE. (Projection operators: automatic) # ENCUT = 400 (Cut-off energy for plane wave basis set, in eV) PREC = Normal (Precision level) LWAVE = .TRUE. (Write WAVECAR or not) LCHARG = .TRUE. (Write CHGCAR or not) ADDGRID= .TRUE. (Increase grid; helps GGA convergence) # LVTOT = .TRUE. (Write total electrostatic potential into LOCPOT or not) # LVHAR = .TRUE. (Write ionic + Hartree electrostatic potential into LOCPOT or not) # NELECT = (No. of electrons: charged cells; be careful) # LPLANE = .TRUE. (Real space distribution; supercells) # NPAR = 4 (Max is no. nodes; don't set for hybrids) # NWRITE = 2 (Medium-level output) # KPAR = 2 (Divides k-grid into separate groups) # NGX = 500 (FFT grid mesh density for nice charge/potential plots) # NGY = 500 (FFT grid mesh density for nice charge/potential plots) # NGZ = 500 (FFT grid mesh density for nice charge/potential plots) Electronic Relaxation ISMEAR = 0 (Gaussian smearing; metals:1) SIGMA = 0.05 (Smearing value in eV; metals:0.2) NELM = 60 (Max electronic SCF steps) NELMIN = 6 (Min electronic SCF steps) EDIFF = 1E-08 (SCF energy convergence; in eV) # GGA = PS (PBEsol exchange-correlation) Ionic Relaxation NSW = 100 (Max ionic steps) IBRION = 2 (Algorithm: 0-MD; 1-Quasi-New; 2-CG) ISIF = 2 (Stress/relaxation: 2-Ions, 3-Shape/Ions/V, 4-Shape/Ions) EDIFFG = -2E-02 (Ionic convergence; eV/AA) # ISYM = 2 (Symmetry: 0=none; 2=GGA; 3=hybrids) |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有10人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有8人回复
-
课题组会帮你通过考核么？
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
【答案】应助回帖
感谢参与，应助指数 +1
|
自洽精度太高了，EDIFF = 1E-08，我一般就是1E-05 18核心算132原子也是压力太大了，一般来说1个核心1个原子效率可以达到最大化，最高也不要超过1个核心2原子 How to choose the number of cores Rule of thumb: 1 atom per core = Good 0.5 atom per core = Could work (but bad efficiency and time wasted) < 0.5 atom per core = Don’t do it https://www.pdc.kth.se/software/ ... la/index_using.html |
2楼2021-08-09 22:03:10
