用VASP-DFPT+phonopy计算石墨烯声子谱求助 - 第一性原理 - Vasp&MedeA

magicmonk

至尊木虫 (著名写手)

应助: 172 (高中生)
金币: 11451.3
红花: 24
帖子: 1124
在线: 3108.8小时
虫号: 1191315
注册: 2011-01-17
专业: 理论和计算化学

引用回帖:

3楼: Originally posted by chuanghua304 at 2015-02-01 11:27:03
首先看文献用的是LDA还是PBE？SIGMA = 0.02太小了，用0.1就行了。还有就是看人家都到底是扩的多大的包。

请问为何sigma影响这么大？

我们最深的恐惧不是我们能力不够，我们最深的恐惧是能力超越了极限。

12楼2015-02-02 09:16:59

liqizuiyang

木虫 (著名写手)

应助: 59 (初中生)
金币: 6859.5
散金: 261
红花: 94
帖子: 2093
在线: 708.9小时
虫号: 2366686
注册: 2013-03-21
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

12楼: Originally posted by magicmonk at 2015-02-02 09:16:59
请问为何sigma影响这么大？...

算声子需要先算Hessian矩阵，也就是总能对原子坐标的二阶偏导。

对于金属，其总能与展宽类型关系很大。VASP手册上明确建议金属用MP型展宽。

13楼2015-02-04 11:59:16

chuanghua304

禁虫 (职业作家)

本帖内容被屏蔽

14楼2015-02-06 18:34:53

conperint

金虫 (正式写手)

应助: 6 (幼儿园)
贵宾: 0.001
金币: 480.3
散金: 128
红花: 6
帖子: 472
在线: 385.5小时
虫号: 1935542
注册: 2012-08-12
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

14楼: Originally posted by chuanghua304 at 2015-02-06 18:34:53
IBRION=8和IBRION=6是一样的，都是用SPOSCAR这个文件，两个计算的时间是差不多的，但是6更省内存些。

了解了，谢谢

15楼2015-02-08 15:07:30

小虫虫LG

木虫 (正式写手)

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

引用回帖:

15楼: Originally posted by conperint at 2015-02-08 15:07:30
了解了，谢谢...

你这问题解决了吗？能说一下吗？

16楼2015-02-23 16:58:49

conperint

金虫 (正式写手)

应助: 6 (幼儿园)
贵宾: 0.001
金币: 480.3
散金: 128
红花: 6
帖子: 472
在线: 385.5小时
虫号: 1935542
注册: 2012-08-12
性别: GG
专业: 凝聚态物性 II ：电子结构

我把我算graphene的POSCAR贴出来给大家看看
POSCAR-unitcell如下：

graphene
1.00
   2.4673063800000000 0.0000000000000000 0.0000000000000000
-1.2336531900000000 2.1367500000000000 0.0000000000000000
   0.0000000000000000 0.0000000000000000 20.0000000000000000
C
2
Direct
0.33333333 0.66666667 0.50000000
0.66666667 0.33333333 0.50000000

用phonopy生成的SPOSCAR如下（我用的4*4超胞）

C
1.0
   9.8692255200000005 0.0000000000000000 0.0000000000000000
-4.9346127600000003 8.5470000000000006 0.0000000000000000
   0.0000000000000000 0.0000000000000000 20.0000000000000000
  32
Direct
  0.0833333325000000  0.1666666675000000  0.5000000000000000
  0.3333333325000000  0.1666666675000000  0.5000000000000000
  0.5833333325000000  0.1666666675000000  0.5000000000000000
  0.8333333325000000  0.1666666675000000  0.5000000000000000
  0.0833333325000000  0.4166666675000000  0.5000000000000000
  0.3333333325000000  0.4166666675000000  0.5000000000000000
  0.5833333325000000  0.4166666675000000  0.5000000000000000
  0.8333333325000000  0.4166666675000000  0.5000000000000000
  0.0833333325000000  0.6666666675000000  0.5000000000000000
  0.3333333325000000  0.6666666675000000  0.5000000000000000
  0.5833333325000000  0.6666666675000000  0.5000000000000000
  0.8333333325000000  0.6666666675000000  0.5000000000000000
  0.0833333325000000  0.9166666675000000  0.5000000000000000
  0.3333333325000000  0.9166666675000000  0.5000000000000000
  0.5833333325000000  0.9166666675000000  0.5000000000000000
  0.8333333325000000  0.9166666675000000  0.5000000000000000
  0.1666666675000000  0.0833333325000000  0.5000000000000000
  0.4166666675000000  0.0833333325000000  0.5000000000000000
  0.6666666675000000  0.0833333325000000  0.5000000000000000
  0.9166666675000000  0.0833333325000000  0.5000000000000000
  0.1666666675000000  0.3333333325000000  0.5000000000000000
  0.4166666675000000  0.3333333325000000  0.5000000000000000
  0.6666666675000000  0.3333333325000000  0.5000000000000000
  0.9166666675000000  0.3333333325000000  0.5000000000000000
  0.1666666675000000  0.5833333325000000  0.5000000000000000
  0.4166666675000000  0.5833333325000000  0.5000000000000000
  0.6666666675000000  0.5833333325000000  0.5000000000000000
  0.9166666675000000  0.5833333325000000  0.5000000000000000
  0.1666666675000000  0.8333333325000000  0.5000000000000000
  0.4166666675000000  0.8333333325000000  0.5000000000000000
  0.6666666675000000  0.8333333325000000  0.5000000000000000
  0.9166666675000000  0.8333333325000000  0.5000000000000000

17楼2015-03-23 17:04:33

conperint

金虫 (正式写手)

应助: 6 (幼儿园)
贵宾: 0.001
金币: 480.3
散金: 128
红花: 6
帖子: 472
在线: 385.5小时
虫号: 1935542
注册: 2012-08-12
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

17楼: Originally posted by conperint at 2015-03-23 17:04:33
我把我算graphene的POSCAR贴出来给大家看看
POSCAR-unitcell如下：

graphene
1.00
2.4673063800000000 0.0000000000000000 0.0000000000000000
-1.2336531900000000 2.136750000000000 ...

扩胞不同，对应band.conf里的设置也有所不同，例如2*2超胞与4*4超胞对应的band.conf是不一样的。4*4的如下
ATOM_NAME = C
DIM = 4 4 1
PRIMITIVE_AXIS = 0.25 0.0 0.0 0.0 0.25 0.0 0.0 0.0 1.0
BAND = 0.0 0.0 0.0 0.08333333 0.08333333 0.0 0.125 0.0 0.0 0.0 0.0 0.0
BAND_POINTS = 200
FORCE_CONSTANTS = READ

2*2的应该如下
ATOM_NAME = C
DIM = 2 2 1
PRIMITIVE_AXIS = 0.5 0.0 0.0 0.0 0.5 0.0 0.0 0.0 1.0
BAND = 0.0 0.0 0.0 0.16666667 0166666667 0.0 0.25 0.0 0.0 0.0 0.0 0.0
BAND_POINTS = 200
FORCE_CONSTANTS = READ

PRIMITIVE_AXIS与BAND这两个参数都有不同。可以这样看，4*0.25=2*0.5=1倍石墨烯单胞晶格常数，所以PRIMITIVE_AXIS这个参数有区别。另外DIM里的数据乘以BAND的坐标应该得到石墨烯单胞的高对称点坐标。所以4*0.0833333=2*0.1666667=0.333333，如你所知，（0.333333,0.333333,0.0）是graphene单胞K点坐标。同理，4*0.125=2*0.25，（0.5,0.0,0.0）是graphene单胞的M点坐标。

以上解释我没法给出确切的物理解释，起到一个方便大家记忆的作用吧

18楼2015-03-23 17:08:48