VASP Co3O4晶面优化EDDAV耗时过长

Lei1 · 发表于 Post on 2024-10-19 18:01:59

请教一下关于VASP优化Co3O4 220晶面，EDDAV耗时过长问题，计算所用服务器为128核，是磁矩、加U还是哪里的设置出现了问题？下面是相关文件

3540.out (4 KB, 下载次数 Times of downloads: 3)

INCAR (2.12 KB, 下载次数 Times of downloads: 11)

KPOINTS (77 Bytes, 下载次数 Times of downloads: 9)

OUTCAR (74.94 KB, 下载次数 Times of downloads: 5)

POSCAR (14.7 KB, 下载次数 Times of downloads: 11)

XIEH · 发表于 Post on 2026-2-13 18:25:18

卡开发发发表于 2026-2-5 16:14
你或许可以试试用nupdown约束总磁矩跑几步看看。

感谢发发老师，设置这个参数成功收敛到了mag=168

卡开发发 · 发表于 Post on 2026-2-5 16:14:44

XIEH 发表于 2026-2-3 13:37
发发老师，今天遇到个新问题，mag收敛不到设定的值，比如Fe36O48，Fe设置12*4 24*5，但是结构优化完又算 ...

你或许可以试试用nupdown约束总磁矩跑几步看看。

XIEH · 发表于 Post on 2026-2-3 13:37:23

卡开发发发表于 2026-2-2 04:47
你可以对照下IBZKPT里面k点数目谁大谁小。

发发老师，今天遇到个新问题，mag收敛不到设定的值，比如Fe36O48，Fe设置12*4 24*5，但是结构优化完又算了SCF，还是收敛到170，这个slab还有必要往下算吗？

XIEH · 发表于 Post on 2026-2-3 13:34:37

卡开发发发表于 2026-2-2 04:47
你可以对照下IBZKPT里面k点数目谁大谁小。

抱歉，忙忘了。之前看的时候G方案要多一个gamma点

卡开发发 · 发表于 Post on 2026-2-2 04:47:54

XIEH 发表于 2026-1-31 22:38
发发老师，最近发现个规律，也不知道对不对，跟您探讨一下。使用3x2x1的gamma-centred网格比3x3x1的一般M ...

你可以对照下IBZKPT里面k点数目谁大谁小。

XIEH · 发表于 Post on 2026-1-31 22:38:19

卡开发发发表于 2026-1-31 07:38
EDDAV应该指的是Davidson对角，但为什么这么慢我这里信息太少判断不了。

发发老师，最近发现个规律，也不知道对不对，跟您探讨一下。使用3x2x1的gamma-centred网格比3x3x1的一般MP网格还慢，是不是过gamma点的计算就会很难算啊？

XIEH · 发表于 Post on 2026-1-31 09:31:50

卡开发发发表于 2026-1-31 07:38
EDDAV应该指的是Davidson对角，但为什么这么慢我这里信息太少判断不了。

谢谢您的回复，我刚去检查了K点，应该是用的太大了，因为之前想一步算DOS，发现算不动，改成先算单点了，但是K点没改小

卡开发发 · 发表于 Post on 2026-1-31 07:38:11

XIEH 发表于 2026-1-30 22:01
又来请教您了，发发老师。这个 “EDDAV: cpu time 5157.7598: real time 6363.9024“ 的EDDAV啥意思？ ...

EDDAV应该指的是Davidson对角，但为什么这么慢我这里信息太少判断不了。

XIEH · 发表于 Post on 2026-1-30 22:01:09

卡开发发发表于 2024-10-20 00:49
1、本身原子数很多，大概有220+个原子，并且金属原子要占接近一半，这种体系VASP速度就会比较慢。
2、你的 ...

又来请教您了，发发老师。这个 “EDDAV: cpu time 5157.7598: real time 6363.9024“ 的EDDAV啥意思？我输入文件都和之前一样，不知为什么算一个电子步都要几个小时了

Lei1 · 发表于 Post on 2024-10-21 10:17:38

卡开发发发表于 2024-10-20 23:37
LREAL=Auto（或者直接写A，因为实际读取可能只读取首字母，比如你写Amagi有可能也不会出错[doge]），没有 ...

非常感谢您的解答，改成Auto,并更换到gamma版本VASP，单步运算100多秒可以结束，到正常值了，感谢

卡开发发 · 发表于 Post on 2024-10-20 23:37:11

Lei1 发表于 2024-10-20 22:49
Global Parameters
SYSTEM = Co3O4-220
ISTART = 1 (Read existing wavefunction, if ther ...

LREAL=Auto（或者直接写A，因为实际读取可能只读取首字母，比如你写Amagi有可能也不会出错[doge]），没有那两个.，True和False情况特殊，在Fortran里面是logical（或是说bool）类型，当然INCAR当中对True和False你直接写True、False或者T、F程序都能正常解读出来。

Lei1 · 发表于 Post on 2024-10-20 22:51:22

zmjsce 发表于 2024-10-20 22:43
可以考虑用小一点的超胞来做，资源不多的情况下100-150原子内差不多了。

好的，后期试一下看看

Lei1 · 发表于 Post on 2024-10-20 22:49:42

卡开发发发表于 2024-10-20 22:10
INCAR文件上传下。

Global Parameters
SYSTEM = Co3O4-220
ISTART =  1          (Read existing wavefunction, if there)
ISPIN  =  2          (Non-Spin polarised DFT)
ICHARG =  1       (Non-self-consistent: GGA/LDA band structures)
LREAL  = .Auto.    (Projection operators: automatic)
ENCUT  =  500       (Cut-off energy for plane wave basis set, in eV)
PREC =  Normal (Precision level: Normal or Accurate, set Accurate when perform structure lattice relaxation calculation)
LWAVE  = .TRUE.       (Write WAVECAR or not)
LCHARG = .TRUE.       (Write CHGCAR or not)
ADDGRID= .TRUE.       (Increase grid, helps GGA convergence)
LVTOT  = .FALSE.    (Write total electrostatic potential into LOCPOT or not)
LVHAR  = .FALSE.    (Write ionic + Hartree electrostatic potential into LOCPOT or not)
LELF = .FALSE.
# NELECT =          (No. of electrons: charged cells, be careful)
# LPLANE = .TRUE.    (Real space distribution, supercells)
# NWRITE = 2          (Medium-level output)
KPAR = 1          (Divides k-grid into separate groups)
# NGXF = 300       (FFT grid mesh density for nice charge/potential plots)
# NGYF = 300       (FFT grid mesh density for nice charge/potential plots)
# NGZF = 300       (FFT grid mesh density for nice charge/potential plots)
NCORE = 8
IVDW = 12

Electronic Relaxation
ISMEAR =  0          (Gaussian smearing, metals:1)
SIGMA  =  0.2       (Smearing value in eV, metals:0.2)
AIGO = Fast
NELM =  300          (Max electronic SCF steps)
NELMIN =  5          (Min electronic SCF steps)
EDIFF  =  1E-05       (SCF energy convergence, in eV)
GGA  =  PE          (PBEsol exchange-correlation)
LDAU = TRUE
LDAUTYPE = 2
LDAUL = 2 -1
LDAUU = 3.5 0.0
LDAUJ = 0.0 0.0
LMAXMIX = 4

Ionic Relaxation
NSW =  300       (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 = -0.05       (Ionic convergence, eV/AA)
ISYM =  0          (Symmetry: 0=none, 2=GGA, 3=hybrids)\
MAGMOM = 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5 5 -5  0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

zmjsce · 发表于 Post on 2024-10-20 22:43:38

可以考虑用小一点的超胞来做，资源不多的情况下100-150原子内差不多了。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] VASP Co3O4晶面优化EDDAV耗时过长

浏览过的版块