版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

国花1989

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 768.8
红花: 1
帖子: 109
在线: 56.8小时
虫号: 2001124
注册: 2012-09-14
性别: MM
专业: 凝聚态物性I:结构、力学和

[求助] VASP+DFPT+phonopy(graphene) 已有2人参与

问题：计算石墨烯在ZA支有虚频，以下是我的计算过程和结果，望大神们指导···并欢迎探讨···

我的计算采用VASP+DFPT+phonopy（IBRION=8然后用phonopy处理）
超胞用的6*6*1,K点为9*9*1
这是我的INCAR
system=osc
ISTART=0
ICHARG=2
ENCUT=500
PREC = Accurate

ISMEAR=0
SIGMA=0.1
EDIFF=1E-08
IALGO=38
NELM=60

NSW=1
IBRION=8
EDIFFG=-1E-04
#POTIM=0.2

#ISIF=2
ADDGRID = .TRUE.

LREAL=.FALSE.
LWAVE=.FALSE.
LCHARG=.FALSE.

GGA=PE
#LMAXTAU=6
LSUBROT=T #
图片是我的计算结果（声子普）

我的计算输出为
entering main loop
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1    0.272182855342E+04 0.27218E+04 -0.10197E+05  2080 0.106E+03
DAV: 2 -0.158855122404E+03 -0.28807E+04 -0.27773E+04  3040 0.252E+02
DAV: 3 -0.692955056835E+03 -0.53410E+03 -0.53287E+03  2400 0.121E+02
DAV: 4 -0.708964673237E+03 -0.16010E+02 -0.15955E+02  2400 0.276E+01
DAV: 5 -0.709397935827E+03 -0.43326E+00 -0.43292E+00  2816 0.403E+00 0.376E+01
DAV: 6 -0.677015902513E+03 0.32382E+02 -0.30477E+01  2144 0.100E+01 0.221E+01
DAV: 7 -0.664463657658E+03 0.12552E+02 -0.41458E+01  2240 0.115E+01 0.568E+00
DAV: 8 -0.664220895831E+03 0.24276E+00 -0.90711E-01  2304 0.195E+00 0.548E-01
DAV: 9 -0.664172638386E+03 0.48257E-01 -0.72933E-02  2336 0.478E-01 0.436E-01
DAV:  10 -0.664099768725E+03 0.72870E-01 -0.84114E-02  2176 0.397E-01 0.358E-01
DAV:  11 -0.664112393672E+03 -0.12625E-01 -0.22063E-02  2208 0.237E-01 0.218E-01
DAV:  12 -0.664136461507E+03 -0.24068E-01 -0.45457E-03  2176 0.128E-01 0.794E-02
DAV:  13 -0.664148178612E+03 -0.11717E-01 -0.58636E-04  2080 0.385E-02 0.493E-02
DAV:  14 -0.664168062230E+03 -0.19884E-01 -0.29950E-03  2048 0.562E-02 0.168E-02
DAV:  15 -0.664171036448E+03 -0.29742E-02 -0.31765E-04  2144 0.190E-02 0.466E-03
DAV:  16 -0.664171436932E+03 -0.40048E-03 -0.22567E-05  2016 0.597E-03 0.404E-03
DAV:  17 -0.664171668598E+03 -0.23167E-03 -0.11372E-05  1952 0.482E-03 0.158E-03
DAV:  18 -0.664171760753E+03 -0.92155E-04 -0.18372E-07  2016 0.886E-04 0.930E-04
DAV:  19 -0.664171817344E+03 -0.56591E-04 -0.58797E-07  1952 0.968E-04 0.499E-04
DAV:  20 -0.664171832640E+03 -0.15296E-04 0.15921E-08  1600 0.358E-04 0.212E-04
DAV:  21 -0.664171834878E+03 -0.22384E-05 0.10949E-07  1728 0.142E-04 0.871E-05
DAV:  22 -0.664171835375E+03 -0.49700E-06 0.12023E-07  1216 0.818E-05 0.405E-05
DAV:  23 -0.664171835374E+03 0.11933E-08 0.12113E-07  1632 0.309E-05 0.335E-05
DAV:  24 -0.664171835378E+03 -0.46202E-08 0.12142E-07  1376 0.143E-05 0.323E-05
DAV:  25 -0.664171835380E+03 -0.13497E-08 0.12267E-07  1504 0.710E-06 0.323E-05
DAV:  26 -0.664171835379E+03 0.78944E-09 0.12150E-07  1568 0.408E-06 0.323E-05
DAV:  27 -0.664171835379E+03 -0.12005E-09 0.12105E-07  1440 0.260E-06 0.323E-05
DAV:  28 -0.664171835379E+03 0.27285E-09 0.12203E-07  1408 0.177E-06 0.323E-05
DAV:  29 -0.664171835379E+03 -0.18554E-09 0.12074E-07  1440 0.130E-06 0.323E-05
DAV:  30 -0.664171835379E+03 -0.34925E-09 0.12095E-07  1408 0.916E-07 0.323E-05
DAV:  31 -0.664171835379E+03 -0.21828E-10 0.12119E-07  1440 0.694E-07 0.323E-05
DAV:  32 -0.664171835379E+03 0.22919E-09 0.12216E-07  1440 0.529E-07 0.323E-05
DAV:  33 -0.664171835380E+03 -0.47658E-09 0.12042E-07  1376 0.414E-07 0.323E-05
DAV:  34 -0.664171835379E+03 0.57116E-09 0.12199E-07  1440 0.319E-07 0.323E-05
DAV:  35 -0.664171835379E+03 -0.76398E-10 0.12135E-07  1440 0.252E-07 0.323E-05
DAV:  36 -0.664171835379E+03 0.24374E-09 0.12205E-07  1408 0.196E-07 0.323E-05
DAV:  37 -0.664171835380E+03 -0.81127E-09 0.12073E-07  1440 0.156E-07 0.323E-05
DAV:  38 -0.664171835379E+03 0.17099E-09 0.12154E-07  1440 0.123E-07 0.323E-05
DAV:  39 -0.664171835379E+03 0.40382E-09 0.12159E-07  1344 0.985E-08 0.323E-05
DAV:  40 -0.664171835379E+03 -0.98225E-10 0.12241E-07  1280 0.775E-08 0.323E-05
DAV:  41 -0.664171835379E+03 -0.33833E-09 0.12078E-07  1344 0.626E-08 0.323E-05
DAV:  42 -0.664171835379E+03 0.43656E-10 0.12095E-07  1312 0.495E-08 0.323E-05
DAV:  43 -0.664171835379E+03 0.36744E-09 0.12211E-07  1280 0.402E-08 0.323E-05
DAV:  44 -0.664171835379E+03 -0.10186E-09 0.12245E-07  1344 0.319E-08 0.323E-05
DAV:  45 -0.664171835379E+03 0.12369E-09 0.12187E-07  1344 0.263E-08 0.323E-05
DAV:  46 -0.664171835379E+03 0.72760E-10 0.12178E-07  1312 0.209E-08 0.323E-05
DAV:  47 -0.664171835379E+03 0.47294E-10 0.12281E-07  1344 0.171E-08 0.323E-05
DAV:  48 -0.664171835379E+03 -0.36380E-09 0.12229E-07  1280 0.137E-08 0.323E-05
DAV:  49 -0.664171835379E+03 0.50568E-09 0.12234E-07  1248 0.113E-08 0.323E-05
DAV:  50 -0.664171835379E+03 0.18190E-09 0.12249E-07  1248 0.913E-09 0.323E-05
DAV:  51 -0.664171835379E+03 -0.25466E-09 0.12181E-07  1248 0.757E-09 0.323E-05
DAV:  52 -0.664171835379E+03 -0.43292E-09 0.12121E-07  1248 0.617E-09 0.323E-05
DAV:  53 -0.664171835379E+03 0.45475E-09 0.12276E-07  1248 0.521E-09 0.323E-05
DAV:  54 -0.664171835379E+03 -0.29468E-09 0.12108E-07  1248 0.430E-09 0.323E-05
DAV:  55 -0.664171835379E+03 0.28740E-09 0.12326E-07  1248 0.367E-09 0.323E-05
DAV:  56 -0.664171835379E+03 -0.37835E-09 0.12157E-07  1248 0.306E-09 0.323E-05
DAV:  57 -0.664171835380E+03 -0.44020E-09 0.12191E-07  1248 0.267E-09 0.323E-05
DAV:  58 -0.664171835379E+03 0.32014E-09 0.12200E-07  1248 0.229E-09 0.323E-05
DAV:  59 -0.664171835379E+03 0.66939E-09 0.12163E-07  1152 0.201E-09 0.323E-05
DAV:  60 -0.664171835379E+03 -0.56025E-09 0.12214E-07  1152 0.175E-09
1 F= -.66417184E+03 E0= -.66416732E+03  d E =-.902703E-02
Linear response reoptimize wavefunctions to high precision
DAV: 1 -0.664171835379E+03 -0.12733E-09 0.12201E-07  1056 0.157E-09
DAV: 2 -0.664171835379E+03 0.37835E-09 0.12343E-07  1056 0.139E-09
DAV: 3 -0.664171835379E+03 -0.25466E-09 0.12138E-07  1056 0.127E-09
Linear response DOF=          2
Linear response progress:
  Degree of freedom: 1/  2
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.148813587331E+00 -0.14881E+00 0.32617E-02  9250 0.857E-01
RMM: 2 -0.122463619931E+00 0.26350E-01 -0.10851E-02  5855 0.224E-01 0.259E-01
RMM: 3 -0.119917510682E+00 0.25461E-02 -0.39110E-03  5328 0.100E-01 0.130E-01
RMM: 4 -0.118931791076E+00 0.98572E-03 -0.56747E-04  5382 0.418E-02 0.493E-02
RMM: 5 -0.118582905734E+00 0.34889E-03 -0.59806E-05  5363 0.155E-02 0.870E-03
RMM: 6 -0.118603091486E+00 -0.20186E-04 -0.35631E-06  5781 0.362E-03 0.213E-03
RMM: 7 -0.118619378680E+00 -0.16287E-04 -0.11216E-06  5941 0.115E-03 0.713E-04
RMM: 8 -0.118619293187E+00 0.85493E-07 -0.15409E-06  6010 0.383E-04 0.271E-04
RMM: 9 -0.118619719101E+00 -0.42591E-06 -0.15407E-06  6115 0.183E-04 0.121E-04
RMM:  10 -0.118619285641E+00 0.43346E-06 -0.13532E-06  6166 0.146E-04 0.607E-05
RMM:  11 -0.118619171547E+00 0.11409E-06 -0.13236E-06  6445 0.119E-04 0.321E-05
RMM:  12 -0.118619180885E+00 -0.93376E-08 -0.15839E-06  6833 0.100E-04 0.191E-05
RMM:  13 -0.118619210317E+00 -0.29432E-07 -0.14772E-06  7348 0.808E-05 0.109E-05
RMM:  14 -0.118619231149E+00 -0.20832E-07 -0.16626E-06  7741 0.688E-05 0.545E-06
RMM:  15 -0.118619219601E+00 0.11549E-07 -0.14926E-06  7782 0.649E-05 0.361E-06
RMM:  16 -0.118619249265E+00 -0.29665E-07 -0.17996E-06  8223 0.602E-05 0.262E-06
RMM:  17 -0.118619204794E+00 0.44472E-07 -0.13753E-06  8067 0.594E-05 0.173E-06
force on displaced ion    1 direction 1 : -61.436  -0.000 0.000
Linear response progress:
  Degree of freedom: 2/  2
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.110434693204E+00 -0.11043E+00 0.27705E-03  3660 0.622E-01
RMM: 2 -0.107524025501E+00 0.29107E-02 -0.37807E-03  3225 0.109E-01 0.175E-01
RMM: 3 -0.106516338541E+00 0.10077E-02 -0.44343E-04  2624 0.406E-02 0.944E-02
RMM: 4 -0.105721296554E+00 0.79504E-03 -0.45926E-04  2742 0.385E-02 0.145E-02
RMM: 5 -0.105922642210E+00 -0.20135E-03 -0.24097E-05  2998 0.912E-03 0.302E-03
RMM: 6 -0.105895526110E+00 0.27116E-04 -0.10421E-06  3271 0.188E-03 0.114E-03
RMM: 7 -0.105911113437E+00 -0.15587E-04 -0.14701E-07  2970 0.633E-04 0.161E-04
RMM: 8 -0.105909853470E+00 0.12600E-05 0.38011E-08  3426 0.143E-04 0.781E-05
RMM: 9 -0.105910392499E+00 -0.53903E-06 0.47833E-08  3420 0.747E-05 0.415E-05
RMM:  10 -0.105910044845E+00 0.34765E-06 -0.10475E-08  3584 0.532E-05 0.204E-05
RMM:  11 -0.105909834667E+00 0.21018E-06 0.24963E-08  3609 0.517E-05 0.138E-05
RMM:  12 -0.105909960334E+00 -0.12567E-06 -0.16925E-08  3804 0.357E-05 0.668E-06
RMM:  13 -0.105909931785E+00 0.28549E-07 0.67272E-09  3895 0.316E-05 0.351E-06
RMM:  14 -0.105909900178E+00 0.31607E-07 0.18023E-08  4110 0.265E-05 0.207E-06
RMM:  15 -0.105909909342E+00 -0.91648E-08 -0.11498E-08  4031 0.245E-05 0.107E-06
force on displaced ion    1 direction 3 : 0.000 0.000 -14.353
Linear response finished

声子普在ZA支有虚频，（1）二维材料用VASP的DFPT（IBRION=8）计算有什么要注意的吗？
                           （2）我的Degree of freedom:  2，我知道二维材料的自由度为2，但是计算声子普是不是三个方向都要震动呢？大神们的计算结果也这样吗？我是刚入门，很纠结··

望大神指导·····

graphene-2fc-DFPT.jpg

graphene-2fc-DFPT-zoom.jpg

回复此楼

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:

DFPT=phonopy计算错误band.conf not found 已经有4人回复
用VASP-DFPT+phonopy计算石墨烯声子谱求助已经有37人回复
关于VASP +phonopy 的两种计算方法 IBRION的设置问题已经有14人回复
vasp+dfpt算声子性质用phonopy软件处理数据时得到band图真的不忍直视，它怎么了？？已经有17人回复
计算声子谱的方法比较 vasp+phonon，espresso，castep, abinit, siesta 已经有0人回复
VASP 设置IBRION=8 没有力常数输出已经有32人回复
vasp+phonopy用dfpt计算声子谱的时候需要建多大的超胞呢已经有2人回复
ZnO声子谱已经有14人回复
【求助】如何用vasp5.2中的DFPT方法算Born effective charge？已经有11人回复

有梦就不要怕痛！！

1楼 2017-03-24 10:07:28

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

回帖支持 ( 显示支持度最高的前 50 名 )

小虫虫LG

木虫 (正式写手)

应助: 114 (高中生)
金币: 1775.9
散金: 1871
红花: 25
帖子: 975
在线: 678.5小时
虫号: 2128430
注册: 2012-11-15
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

感谢参与，应助指数 +1

虚频大约几个cm？小的话当作看不见。

赞一下(2人)

回复此楼

» 本帖已获得的红花（最新10朵）

国花1989

9楼2017-03-24 20:07:36

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

前田庆次

木虫 (著名写手)

应助: 202 (大学生)
金币: 395.2
散金: 1140
红花: 88
帖子: 1344
在线: 455.2小时
虫号: 1985685
注册: 2012-09-08
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
漫天飘雪: 金币+2, 谢谢交流 2017-03-25 17:36:53
国花1989: 金币+10, ★★★★★最佳答案 2017-04-05 11:19:46

不太会算这东西，三个建议
1. 优化了吗
2. 扩胞
3. SIGMA是不是有点大

赞一下(1人)

回复此楼

» 本帖已获得的红花（最新10朵）

国花1989

长枪一横花飘零，松风追月伴我行。无双人间世如梦，倾奇万世永留名。

8楼2017-03-24 19:33:02

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

普通回帖

zhangfrank

金虫 (正式写手)

应助: 178 (高中生)
金币: 660.2
散金: 224
红花: 53
帖子: 891
在线: 344.5小时
虫号: 2679423
注册: 2013-09-25
性别: GG
专业: 凝聚态物性 II ：电子结构

不了解这个帮忙顶一下吧

赞一下

回复此楼

计算顺利点吧

2楼2017-03-24 10:40:27

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

国花1989

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 768.8
红花: 1
帖子: 109
在线: 56.8小时
虫号: 2001124
注册: 2012-09-14
性别: MM
专业: 凝聚态物性I:结构、力学和

引用回帖:

2楼: Originally posted by zhangfrank at 2017-03-24 10:40:27
不了解这个帮忙顶一下吧

谢谢，顶起来····

赞一下

回复此楼

有梦就不要怕痛！！

3楼2017-03-24 11:44:12

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

léon_WANG

木虫 (小有名气)

应助: 5 (幼儿园)
金币: 2577.5
散金: 20
红花: 4
帖子: 156
在线: 226.1小时
虫号: 2458360
注册: 2013-05-10
专业: 凝聚态物性I:结构、力学和

一点点虚频，提高精度算就行了

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

国花1989

4楼2017-03-24 15:10:26

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

国花1989

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 768.8
红花: 1
帖子: 109
在线: 56.8小时
虫号: 2001124
注册: 2012-09-14
性别: MM
专业: 凝聚态物性I:结构、力学和

送红花一朵

引用回帖:

4楼: Originally posted by léon_WANG at 2017-03-24 15:10:26
一点点虚频，提高精度算就行了

精度已经很高了，优化能量收敛用10-8,力收敛用10-4,超胞是用6*6*1，我在想是不是技术上的问题，精度已经足够了！

发自小木虫Android客户端

赞一下

回复此楼

有梦就不要怕痛！！

5楼2017-03-24 15:25:37

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

123dingyc

新虫 (知名作家)

应助: 0 (幼儿园)
金币: 16124.1
红花: 44
帖子: 6584
在线: 571小时
虫号: 2516270
注册: 2013-06-21
专业: 凝聚态物性I:结构、力学和

再扩包

发自小木虫IOS客户端

回复此楼

» 本帖已获得的红花（最新10朵）

国花1989

6楼2017-03-24 15:28:02

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

国花1989

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 768.8
红花: 1
帖子: 109
在线: 56.8小时
虫号: 2001124
注册: 2012-09-14
性别: MM
专业: 凝聚态物性I:结构、力学和

引用回帖:

6楼: Originally posted by 123dingyc at 2017-03-24 15:28:02
再扩包

Linear response reoptimize wavefunctions to high precision
DAV: 1 -0.664171835379E+03 -0.12733E-09 0.12201E-07  1056 0.157E-09
DAV: 2 -0.664171835379E+03 0.37835E-09 0.12343E-07  1056 0.139E-09
DAV: 3 -0.664171835379E+03 -0.25466E-09 0.12138E-07  1056 0.127E-09
Linear response DOF=          2
Linear response progress:
  Degree of freedom: 1/  2
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.148813587331E+00 -0.14881E+00 0.32617E-02  9250 0.857E-01
RMM: 2 -0.122463619931E+00 0.26350E-01 -0.10851E-02  5855 0.224E-01 0.259E-01
RMM: 3 -0.119917510682E+00 0.25461E-02 -0.39110E-03  5328 0.100E-01 0.130E-01
RMM: 4 -0.118931791076E+00 0.98572E-03 -0.56747E-04  5382 0.418E-02 0.493E-02
RMM: 5 -0.118582905734E+00 0.34889E-03 -0.59806E-05  5363 0.155E-02 0.870E-03
RMM: 6 -0.118603091486E+00 -0.20186E-04 -0.35631E-06  5781 0.362E-03 0.213E-03
二维材料，是不是这里的自由度都是2啊？我在怀疑这里是不是有问题？

赞一下

回复此楼

有梦就不要怕痛！！

7楼2017-03-24 16:40:04

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