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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp+phonopy计算声子谱虚频求助
91/1返回列表
查看: 4297  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ft1071

新虫 (初入文坛)

[求助] vasp+phonopy计算声子谱虚频求助 已有2人参与

正在计算Ta2C（六方结构）的热膨胀系数，首先要计算11组体积不同（改变POSCAR基矢缩放系数）的声子谱
然后在缩放比较小的时候就在G点附近产生了虚频
缩放0.95：
vasp+phonopy计算声子谱虚频求助
0.95.jpg
缩放1.1：
vasp+phonopy计算声子谱虚频求助-1
1.1.jpg
缩放1.05没有虚频

扩胞3 3 3 ，INCAR参照phonopy手册， K点取3 3 3

计算声子谱的INCAR：
PREC = Accurate
IBRION = -1
ENCUT = 500
EDIFF = 1.0e-08
ISMEAR = 0; SIGMA = 0.01
IALGO = 38
LREAL = .FALSE.
ADDGRID = .TRUE.
LWAVE = .FALSE.
LCHARG = .FALSE.

之前优化的INCAR：
SYSTEM=TaC
ENCUT=550
ISMEAR = 1
SIGMA = 0.01
NSW=41
IBRION = 2
ISIF = 3
EDIFF=1.0E-4
EDIFFG=-0.001
PREC =Accurate

师兄说只要有虚频结果就不可靠，请问是哪个参数还是优化的问题呢？小白只能仰仗各位了
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-02-02 12:50:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

人人余yan

金虫 (著名写手)
结构优化的时候精度可以再高一些吧。十的负七或者负八

发自小木虫Android客户端
赞一下(2人) 回复此楼
2楼2016-02-02 22:04:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ft1071

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by 人人余yan at 2016-02-02 22:04:23
结构优化的时候精度可以再高一些吧。十的负七或者负八

请问EDIFFG需要改吗？
赞一下 回复此楼
3楼2016-02-03 02:19:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

人人余yan

金虫 (著名写手)
引用回帖:
3楼: Originally posted by ft1071 at 2016-02-03 02:19:02
请问EDIFFG需要改吗？...

精度高些更准确，结构太大的话算的慢，也可以是一下ibrion=8/6，为了消除虚频的话

发自小木虫IOS客户端
赞一下 回复此楼
4楼2016-02-03 05:27:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ft1071

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by 人人余yan at 2016-02-03 05:27:52
精度高些更准确，结构太大的话算的慢，也可以是一下ibrion=8/6，为了消除虚频的话
...

这周我增加K点到5 5 5算了两个点，还是有这样的虚频。
重新优化了一下原胞，这是原先的优化：
Ta C                                 
   1.00000000000000     
     3.1244339623514183   -0.0000000000000000    0.0000000000000000
    -1.5622169811757092    2.7058391838478966   -0.0000000000000000
    -0.0000000000000000   -0.0000000000000000    4.9634578427587872
     2     1
Direct
  0.3333333429999996  0.6666666870000029  0.2526072086980557
  0.6666666269999979  0.3333333129999971  0.7473927913019447
-0.0000000000000000 -0.0000000000000000 -0.0000000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00

EDIFF=10E-8，EDIFFG=-10E-6得到的原胞：

Ta C                                   
   1.00000000000000     
     3.1237671313434077   -0.0000000000000000   -0.0000000000000000
    -1.5618835656717038    2.7052616912549250   -0.0000000000000000
     0.0000000000000000    0.0000000000000000    4.9634839229430820
   Ta   C
     2     1
Direct
  0.3333333429999996  0.6666666870000029  0.2526028572060771
  0.6666666269999979  0.3333333129999971  0.7473971427939231
  0.0000000000000000 -0.0000000000000000  0.0000000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00

因为计算很费时，想先问您一下这个差别算大吗？别的参数是否需要修改呢？
——————————————————————————————————————————
另外IBRION=8/6 的话，就是参照phonopy手册里的VASP-DFPT & phonopy calculation 部分？这是手册给的INCAR，把IBRION改成6就可以了吗？
PREC = Accurate
ENCUT = 500
61
phonopy manual, Release 1.7.4
IBRION = 8
EDIFF = 1.0e-08
IALGO = 38
ISMEAR = 0; SIGMA = 0.1
LREAL = .FALSE.
ADDGRID = .TRUE.
LWAVE = .FALSE.
LCHARG = .FALSE.
赞一下 回复此楼
5楼2016-02-03 17:12:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小虫虫LG

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
ft1071: 金币+3 2016-02-14 19:28:46
一个办法是缩放系数变小点把。比如0.99
赞一下 回复此楼
6楼2016-02-11 16:28:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

navyzhang

木虫 (正式写手)
引用回帖:
6楼: Originally posted by 小虫虫LG at 2016-02-11 16:28:02
一个办法是缩放系数变小点把。比如0.99

请问缩放系数在哪里设置啊？
我用phonopy除了改变K点路径，都用的是默认设置
赞一下 回复此楼
7楼2016-06-24 15:54:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Katherinexhl

新虫 (初入文坛)
蛋蛋小童鞋: 应助指数-1 2016-08-04 10:38:10
引用回帖:
7楼: Originally posted by navyzhang at 2016-06-24 15:54:52
请问缩放系数在哪里设置啊？
我用phonopy除了改变K点路径，都用的是默认设置...

请问问题解决了吗？这11个体积用什么优化的呢？是3吗？
赞一下 回复此楼
8楼2016-08-01 17:07:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Katherinexhl

新虫 (初入文坛)

【答案】应助回帖

请问您是用ISIF=?优化的？
赞一下 回复此楼
9楼2016-08-05 20:01:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ft1071 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料080500调剂求收留 +6 一颗meteor 2026-03-13 6/300 2026-03-20 10:41 by EBSD
[考研] 一志愿武汉理工材料工程专硕调剂 +7 Doleres 2026-03-19 7/350 2026-03-20 10:39 by EBSD
[考研] 295复试调剂 +5 简木ChuFront 2026-03-19 5/250 2026-03-20 10:38 by 无际的草原
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考博] 东华理工大学化材专业26届硕士博士申请 +8 zlingli 2026-03-13 8/400 2026-03-19 16:32 by 轻松不少随
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 287求调剂 +3 晨昏线与星海 2026-03-19 4/200 2026-03-19 12:32 by peike
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[基金申请] 被我言中：新模板不强调格式了，假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856专硕279求调剂 +5 加油加油！? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
信息提示
关闭
请填处理意见
关闭