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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 用VASP-DFPT+phonopy计算石墨烯声子谱求助
382/4返回列表上一页1234下一页
查看: 11961  |  回复: 37
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

conperint

金虫 (正式写手)
送红花一朵
引用回帖:
9楼: Originally posted by liqizuiyang at 2015-02-01 18:00:26
石墨烯为半金属，用MP型展宽力常数算得更准确。

11楼2015-02-01 18:44:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)
引用回帖:
12楼: Originally posted by magicmonk at 2015-02-02 09:16:59
请问为何sigma影响这么大？...

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

对于金属，其总能与展宽类型关系很大。VASP手册上明确建议金属用MP型展宽。
赞一下 回复此楼
13楼2015-02-04 11:59:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuanghua304

禁虫 (职业作家)
本帖内容被屏蔽
14楼2015-02-06 18:34:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

conperint

金虫 (正式写手)
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小虫虫LG

木虫 (正式写手)
引用回帖:
15楼: Originally posted by conperint at 2015-02-08 15:07:30
了解了，谢谢...

你这问题解决了吗？能说一下吗？
赞一下 回复此楼
16楼2015-02-23 16:58:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

conperint

金虫 (正式写手)
我把我算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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

conperint

金虫 (正式写手)
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huterx

金虫 (小有名气)
顺便问问你 PRIMITIVE_AXIS 这个参数是怎么设定的？
赞一下 回复此楼
19楼2015-04-27 17:35:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

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

IBRION=6就不是DFPT了吧？
赞一下 回复此楼
未来就是现在
20楼2015-05-27 09:40:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 conperint 的主题更新
382/4返回列表上一页1234下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +7 yufeiwaner 2026-08-09 7/350 2026-08-10 09:08 by zju2000
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于代码变化问题，想知道的进来 +15 且听虎啸 2026-08-07 22/1100 2026-08-09 22:13 by 布布和一二
[基金申请] 这样的filecode谁见过 +10 布布和一二 2026-08-08 21/1050 2026-08-09 17:21 by 天神眷顾
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 纯娱乐，不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
关闭
请填处理意见
关闭