| 查看: 2115 | 回复: 5 | |||||
ultraq至尊木虫 (著名写手)
|
[求助]
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 |
» 猜你喜欢
考研调剂
已经有3人回复
-
317分 一志愿南理工材料工程 本科湖工大 求调剂
已经有9人回复
-
297求调剂
已经有10人回复
-
一志愿武汉理工，总分321，英一数二，求老师收留。
已经有10人回复
-
343求调剂
已经有5人回复
-
085601材料工程找调剂
已经有7人回复
-
求调剂，一志愿 南京航空航天大学 ，080500材料科学与工程学硕，总分289分
已经有6人回复
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研经历
已经有5人回复
-
340求调剂
已经有5人回复
-
一志愿北化求调剂
已经有4人回复
awdawdawdd
金虫 (小有名气)
- 应助: 24 (小学生)
- 金币: 679
- 散金: 200
- 红花: 6
- 帖子: 111
- 在线: 67小时
- 虫号: 2237233
- 注册: 2013-01-12
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
ultraq: 金币+10, ★有帮助, 谢谢！ 2015-07-19 08:40:03
感谢参与，应助指数 +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的计算资源咯 |
2楼2015-07-19 03:36:17
magicmonk
至尊木虫 (著名写手)
- 应助: 172 (高中生)
- 金币: 11451.3
- 红花: 24
- 帖子: 1124
- 在线: 3108.8小时
- 虫号: 1191315
- 注册: 2011-01-17
- 专业: 理论和计算化学
3楼2015-07-19 10:14:05
huterx
金虫 (小有名气)
- 应助: 9 (幼儿园)
- 金币: 1485.5
- 散金: 300
- 红花: 2
- 帖子: 193
- 在线: 432.8小时
- 虫号: 1380155
- 注册: 2011-08-26
- 专业: 凝聚态物性 II ：电子结构
4楼2015-07-19 14:30:32
5楼2015-07-23 09:19:57
ultraq
至尊木虫 (著名写手)
- 应助: 30 (小学生)
- 金币: 14971.9
- 散金: 2033
- 红花: 15
- 帖子: 2378
- 在线: 1050.9小时
- 虫号: 489270
- 注册: 2008-01-02
- 性别: GG
- 专业: 功能陶瓷
6楼2015-07-23 09:38:07
10