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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【其他】在线答疑:VASP的问题
31520/32返回列表首页上一页181920212223下一页末页
查看: 16424  |  回复: 314
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

星愿9997

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
你好，我想问问我用vasp+fropho计算声子谱，得到的结果有虚频，结构优化精度EDIFF=10E-6,EDIFFG=-0.01。怎样消除虚频呢
赞一下(1人) 回复此楼
191楼2011-12-21 15:49:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞
★
youzhizhe(金币+1): 谢谢交流。 2011-12-31 14:16:05
引用回帖:
191楼: Originally posted by 星愿9997 at 2011-12-21 15:49:54:
你好，我想问问我用vasp+fropho计算声子谱，得到的结果有虚频，结构优化精度EDIFF=10E-6,EDIFFG=-0.01。怎样消除虚频呢

如果是再gamma点附近有虚频的话， 那可能是假虚，可以通过加大k点或需要用更大的supercell。    如果是非gamma点的话，有可能真虚。 那就比较复杂了。 需要找到虚频的振动方向，然后通过移动原子，来找到下一个local minima。
赞一下(1人) 回复此楼
真理是一点点接近的！
192楼2011-12-21 22:34:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星愿9997

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
192楼: Originally posted by gavinliu7390 at 2011-12-21 22:34:52:
如果是再gamma点附近有虚频的话， 那可能是假虚，可以通过加大k点或需要用更大的supercell。    如果是非gamma点的话，有可能真虚。 那就比较复杂了。 需要找到虚频的振动方向，然后通过移动原子，来找到下一个l ...

我试过加大K点没有用，我得到的是比较大的虚频。我传了一张我的声子谱图您看看，它是从-8THz到0THz一直有连续的虚频。怎么找到虚频的振动放向，fropho计算后不是只输出dos这一个数据吗？

声子谱图
赞一下(1人) 回复此楼
193楼2011-12-22 09:46:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞
★ ★
youzhizhe(金币+2): 谢谢交流。 2011-12-31 14:16:11
引用回帖:
193楼: Originally posted by 星愿9997 at 2011-12-22 09:46:32:
我试过加大K点没有用，我得到的是比较大的虚频。我传了一张我的声子谱图您看看，它是从-8THz到0THz一直有连续的虚频。怎么找到虚频的振动放向，fropho计算后不是只输出dos这一个数据吗？
[eimg]56/53/1423352_13 ...

你需要看manual， 里面有个开关是控制本征值的， 需要本征值和本征矢， 找到虚频最大的频率的本征矢量。然后再distortion.
从dos里看不出真虚还是假虚。   需要画个能带， 才能分辨。
赞一下(1人) 回复此楼
真理是一点点接近的！
194楼2011-12-22 11:03:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星愿9997

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
194楼: Originally posted by gavinliu7390 at 2011-12-22 11:03:55:
你需要看manual， 里面有个开关是控制本征值的， 需要本征值和本征矢， 找到虚频最大的频率的本征矢量。然后再distortion.
从dos里看不出真虚还是假虚。   需要画个能带， 才能分辨。

这是我计算得到的能带图，您看看。还有哪有manual啊？
赞一下(1人) 回复此楼
195楼2011-12-22 16:51:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星愿9997

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
194楼: Originally posted by gavinliu7390 at 2011-12-22 11:03:55:
你需要看manual， 里面有个开关是控制本征值的， 需要本征值和本征矢， 找到虚频最大的频率的本征矢量。然后再distortion.
从dos里看不出真虚还是假虚。   需要画个能带， 才能分辨。

这是我的声子谱图，刚刚没传上去

声子谱图
赞一下(1人) 回复此楼
196楼2011-12-22 16:54:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞
★ ★
youzhizhe(金币+2): 谢谢交流。 2011-12-31 14:16:15
引用回帖:
196楼: Originally posted by 星愿9997 at 2011-12-22 16:54:47:
这是我的声子谱图，刚刚没传上去

声子谱图

我估计还是你计算的参数有问题。
你扩完的超胞，a，b，c都是多少？
你可以找一下manual， 把负的位移加上。 可能会好些。
赞一下(1人) 回复此楼
真理是一点点接近的！
197楼2011-12-22 23:08:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞
★ ★
youzhizhe(金币+2): 谢谢交流。 2011-12-31 14:16:23
引用回帖:
196楼: Originally posted by 星愿9997 at 2011-12-22 16:54:47:
这是我的声子谱图，刚刚没传上去

声子谱图

More displacements to obtain precise frequencies

You can specify displacements more than independent freedom to obtain precise eigenvalues (frequencies). This treatment works especially near -point. An effective example is  displacements. If your DISP file is,
          1   1.00   0.00   0.00
          1   0.00   1.00   0.00
          1   0.00   0.00   1.00
You can create displacements by hand (or use sed), e.g.,
          1   1.00   0.00   0.00
          1  -1.00   0.00   0.00
          1   0.00   1.00   0.00
          1   0.00  -1.00   0.00
          1   0.00   0.00   1.00
          1   0.00   0.00  -1.00
This setting may be called 'plus-minus'. The plus-minus setting works to increase stability of phonon band structure even when the initial atomic geometry is not enough converged to the equilibrium position. Especially frequency equivalent to approximately zero is sensitive to this treatment. Sometimes fails when forces are converged too much.

加上负方向的受力计算， 会使计算的声子谱更准确。
赞一下(1人) 回复此楼
真理是一点点接近的！
198楼2011-12-22 23:24:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星愿9997

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
197楼: Originally posted by gavinliu7390 at 2011-12-22 23:08:12:
我估计还是你计算的参数有问题。
你扩完的超胞，a，b，c都是多少？
你可以找一下manual， 把负的位移加上。 可能会好些。

我建的超胞就是1*1*1的，a=b=7.612,c=11.213
我传了我的文件，EIGENVECTORS太大没有传。INPHON是计算力时的输入文件，INPHON1是计算声子谱性质的文件
赞一下(1人) 回复此楼
199楼2011-12-23 11:41:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞
★
youzhizhe(金币+1): 谢谢交流。 2011-12-31 14:15:47
引用回帖:
199楼: Originally posted by 星愿9997 at 2011-12-23 11:41:50:
我建的超胞就是1*1*1的，a=b=7.612,c=11.213
我传了我的文件，EIGENVECTORS太大没有传。INPHON是计算力时的输入文件，INPHON1是计算声子谱性质的文件

最好建立2*2*2的超胞！
赞一下(1人) 回复此楼
真理是一点点接近的！
200楼2011-12-23 20:01:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gavinliu7390 的主题更新
31520/32返回列表首页上一页181920212223下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子 +14 botar 2026-07-14 15/750 2026-07-19 00:15 by beefly
[基金申请] 面上的水很深 +13 布布和一二 2026-07-18 18/900 2026-07-18 20:37 by 日月天地
[基金申请] 国社科通讯评审 +4 潜水冠军123 2026-07-18 5/250 2026-07-18 19:58 by chenmodeyu
[基金申请] 5个%2F，非固定段22个字母 +8 jaasxulu 2026-07-15 9/450 2026-07-18 17:40 by dxcgbm
[基金申请] 生命口C13面上会评时间 +7 luckinging 2026-07-14 11/550 2026-07-17 20:43 by 虫卜1988
[教师之家] 内心匮乏 +16 水冰月月野兔 2026-07-12 16/800 2026-07-17 19:25 by otani
[基金申请] 生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？ +6 碧水00 2026-07-16 9/450 2026-07-17 15:21 by jaasxulu
[论文投稿] 投英文期刊如何查重 10+5 不知道取什么名� 2026-07-15 6/300 2026-07-17 09:07 by bobvan
[基金申请] 不要再数国自然申请书的 filecode 的分隔符个数了 +14 sadpencil 2026-07-12 21/1050 2026-07-17 00:15 by wangze12014
[有机交流] 伯胺上甲基求助 100+4 002tmumu 2026-07-14 7/350 2026-07-16 19:47 by 002tmumu
[论文投稿] 跨出版社商投稿 20+5 倪好520 2026-07-13 5/250 2026-07-16 17:24 by 无视军团
[基金申请] 生生命学部会评 +4 wuke100666 2026-07-15 4/200 2026-07-16 00:12 by 啦啦小妮子
[基金申请] H口，面上，时间戳7月7日。 +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[基金申请] 会评 +5 无名者登山 2026-07-15 5/250 2026-07-15 18:15 by 杠就是你对
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[基金申请] E口青C会评结束了吗？ 40+3 小小书虫c 2026-07-14 9/450 2026-07-15 16:04 by 安静的知了
[论文投稿] 农业机械学报投稿周期 +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[基金申请] 国自然面上D口祈祷 +4 monkeynana 2026-07-14 4/200 2026-07-15 09:14 by 白水煮
[基金申请] 没收到消息，再次陪跑 +3 天空星辰fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
[基金申请] 明天E口面上会评 +8 布布和一二 2026-07-12 9/450 2026-07-13 22:59 by QTSorange
信息提示
关闭
请填处理意见
关闭