24小时热门版块排行榜

返回列表

ultraq

至尊木虫 (著名写手)

应助: 30 (小学生)
金币: 14971.9
散金: 2033
红花: 15
帖子: 2378
在线: 1050.9小时
虫号: 489270
注册: 2008-01-02
性别: GG
专业: 功能陶瓷

[求助] VASP+Phonopy算力常数很慢已有3人参与

第一次用Phonopy算声子谱，发现在算力常数特别慢，48小时了才算了几个displacement，按照这架势下去，一个星期都不够啊

我是2个10核cpu并联，怀疑是INCAR有问题，请高手看看：
# General/output settings:
SYSTEM    = ，，，
LORBIT    = 11    ! write magnetic moments in OUTCAR
ENCUT = 500.0    ! in eV
ISMEAR  = 0       ! 0 - gauss ; 1 - MP ;  -5 - tetrahedron
SIGMA = 0.05    !0.1       ! in eV
ISPIN = 2
MAGMOM  =  32*0 16*5 16*0 80*0
#LSORBIT =    T #spin-orbit coupling
#MAGMOM  =  12*0  0 0 9 0 0 -9 0 0 -9 0 0 9 36*0
#LCALCPOL = .TRUE.
#LEPSILON = .TRUE.
GGA=PS

# LDA+U settings:
LDAU    = .TRUE.
LDAUTYPE = 2 ! Dudarev approach, only U-J is used
LDAUL =  -1  2 -1  -1
LDAUU =  0.0 4.5 0.0 0.0

NELMDL=-5
EDIFF=1E-8
NELM=150
LWAVE=.FALSE.
LCHARG=.FALSE.
VOSKOWN=1

PREC=Accurate
NGX=82
NGY=106
NGZ=78

#AMIX  =  0.2
#BMIX  =  0.0001 ! almost zero, but 0 will crash some versions
#AMIX_MAG  =  0.8
#BMIX_MAG  =  0.0001 ! almost zero, but 0 will crash some versions

# Ionic/cell relaxation:
NSW    = 1
IBRION  = 8    ! 0 - MD
            ! 1 - RMM-DIIS method
            ! 2 - CG method
            ! 3 - damped MD method
ISIF = 3    ! relax ions + calculate stress tensor
EDIFFG  = -0.005 ! forces criteria - 10 meV/Ang
#NPAR = 4
LMAXMIX = 4 ! for s p  LMAXMIX to 2,d: LMAXMIX = 4,  f: LMAXMIX = 6

LREAL= .FALSE
ADDGRID = .TRUE.
LASPH = .TRUE.
#ISYM =  -1 0-nonsym 1-usesym 2-fastsym

谢谢！

回复此楼

» 收录本帖的淘帖专辑推荐

仿真建模与计算

phonopy

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有8人回复
最失望的一年已经有12人回复
存款400万可以在学校里躺平吗已经有29人回复
求推荐英文EI期刊已经有5人回复
请教限项目规定已经有4人回复
国自然申请面上模板最新2026版出了吗？已经有20人回复
26申博已经有3人回复
基金委咋了？2026年的指南还没有出来？已经有10人回复
基金申报已经有6人回复
疑惑？已经有5人回复

1楼 2015-07-18 19:40:23

已阅关注TA 给TA发消息送TA红花 TA的回帖

awdawdawdd

金虫 (小有名气)

应助: 24 (小学生)
金币: 679
散金: 200
红花: 6
帖子: 111
在线: 67小时
虫号: 2237233
注册: 2013-01-12
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
ultraq: 金币+10, ★有帮助, 谢谢！ 2015-07-19 08:40:03

Phonon calculation本来就是很要时间的啊，尤其是supercell建的比较大的情况
这是官方给出的例子：
PREC = Accurate
ENCUT = 500
IBRION = 8
  EDIFF = 1.0e-08
  IALGO = 38
ISMEAR = 0; SIGMA = 0.1
  LREAL = .FALSE.
ADDGRID = .TRUE.
  LWAVE = .FALSE.
LCHARG = .FALSE.

和LZ的INCAR比较的话发现lz开了自旋（ISPIN=2），还有DFT+U。自旋肯定会加大计算量，DFT+U没算过不清楚
而且我还发现你设置了ionic relaxation相关项（NSW ISIF EDIFFG），这些貌似不用设置的吧，可以注释掉

LZ嫌慢的话可以用finite displacement方法来算（就是直接用phonopy生成的POSCAR-*）这个可以多个displacement一起算，速度就会快很多啦当然取决于lz的计算资源咯