| 查看: 1403 | 回复: 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人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
【答案】应助回帖
感谢参与，应助指数 +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
20