24小时热门版块排行榜    

查看: 2907  |  回复: 9

hgfrtuug

银虫 (初入文坛)

[求助] 频率分析计算时虚频过大怎么办 已有1人参与

13 f/i=    0.295923 THz     1.859342 2PiTHz    9.870942 cm-1     1.223841 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000 15.000000    -0.008433    0.022563   -0.307981
      0.935885  1.621000 12.403390    -0.015800    0.032033   -0.401723
      1.871770  0.000000 17.596610    -0.013088    0.028841   -0.438809
      0.935885  1.621000 16.290960    -0.016631    0.043391   -0.534236
      1.871770  0.000000 13.709040    -0.018835    0.046475   -0.508204

  14 f/i=    3.550214 THz    22.306655 2PiTHz  118.422402 cm-1    14.682512 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000 15.000000    -0.549241    0.821441    0.025603
      0.935885  1.621000 12.403390    -0.072545    0.065059    0.009396
      1.871770  0.000000 17.596610    -0.070262    0.064003    0.008797
      0.935885  1.621000 16.290960    -0.039250    0.009443    0.015130
      1.871770  0.000000 13.709040    -0.044662    0.010918    0.015682

  15 f/i=    4.741184 THz    29.789735 2PiTHz  158.148854 cm-1    19.607966 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000 15.000000     0.830121    0.552150    0.013299
      0.935885  1.621000 12.403390     0.034219    0.033258    0.000083
      1.871770  0.000000 17.596610     0.029014    0.034399   -0.002635
      0.935885  1.621000 16.290960    -0.031509   -0.004526    0.002676
      1.871770  0.000000 13.709040    -0.021512   -0.005668    0.005238

该结构虚频过大,已经把第一个原子位置移动过重新进行优化,频率计算,结果虚频还是很大,该怎么办?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

漫天飘雪

版主 (知名作家)

认真做事,踏实做人

可能本身就是不稳定的吧
http://blog.sciencenet.cn/u/pfliu89
2楼2017-11-07 16:50:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

前田庆次

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
不擅长声子计算,但算过几次,解决方法
0. 再次优化
1. 扩胞
2. 好像有个什么冷冻原子法一类的可以消虚频,但我不会
长枪一横花飘零,松风追月伴我行。无双人间世如梦,倾奇万世永留名。
3楼2017-11-07 17:09:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

你这问题说得不清楚,到底是虚频多了还是不该有虚频?
如果你是优化结构的,首先要保证模型是合理的,初始结构不合理优化未必收敛到你要的结构,其次,如果初始结构合理,可以提高收敛精度继续优化试一下,或者按照相应虚频给出的dx,dy,dz,取其一半,添加到POSCAR相应原子坐标上去,再优化。注意坐标形式,如果POSCAR里是分数坐标,则取dx.dy.dz一半后须转换为分数坐标单位再相加。可以从最大虚频尝试一下。
4楼2017-11-07 17:10:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

如果是过渡态消除多余虚频,也可以仿照微调坐标来做
5楼2017-11-07 17:12:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

你的体系里不会只有5个原子吧?如果是,你需要固定其中一个
6楼2017-11-07 17:14:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hgfrtuug

银虫 (初入文坛)

引用回帖:
4楼: Originally posted by hakuna at 2017-11-07 17:10:13
你这问题说得不清楚,到底是虚频多了还是不该有虚频?
如果你是优化结构的,首先要保证模型是合理的,初始结构不合理优化未必收敛到你要的结构,其次,如果初始结构合理,可以提高收敛精度继续优化试一下,或者按照 ...

算过了,虚频还是超级大

发自小木虫Android客户端
7楼2017-11-08 15:54:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaq2800

铁杆木虫 (知名作家)

提高精度再算
8楼2017-11-08 16:29:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinm158

新虫 (正式写手)

引用回帖:
7楼: Originally posted by hgfrtuug at 2017-11-08 15:54:15
算过了,虚频还是超级大
...

请问楼主最后解决了吗

发自小木虫IOS客户端
9楼2019-09-06 20:03:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李佳佳佳

金虫 (正式写手)

Alice

停不下的马蹄,滴答滴答~
10楼2019-09-17 09:30:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hgfrtuug 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jRl3mE6ddZGq 2026-08-03 4/200 2026-08-04 04:54 by k7OM8YghWbkC
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 04:50 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:46 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 5/250 2026-08-04 04:44 by k7OM8YghWbkC
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 9/450 2026-08-04 03:56 by kf2781974
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 03:35 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见