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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】有人用有限差分方法计算过G点声子谱吗？
111/2返回列表12下一页
查看: 2032  |  回复: 10
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhang668

木虫 (正式写手)

[交流] 【求助】有人用有限差分方法计算过G点声子谱吗？

如题，我设置IBRION=5，计算了立方BaTiO3的G点声子谱，结果发现了些问题：
INCAR：SYSTEM = BaTiO3
ISTART = 0
ICHARG = 2
ENCUT = 500
LDA = CA
NELMIN = 4
NSW    = 80
IBRION = 5
POTIM = 0.01
NFREE = 2
PREC   = High
ALGO   = Fast
ISIF   = 2
ISMEAR = -5
SIGMA = 0.01
LREAL  = Auto
EDIFF = 1E-5
EDIFFG =-1E-3
AMIX = 0.2
BMIX = 0.0001
AMIX_MAG = 0.8
BMIX_MAG = 0.0001

就是得到的15个特征向量都是沿一个方向的，这个明显不是简谐振动模式，比如氧八面体扭曲模式O的移动就不可能是在一个方向上的，而且没有发现平移的声子模式。所以请问这个跟声子有什么区别？下面是结果：
Eigenvectors and eigenvalues of the dynamical matrix
----------------------------------------------------


   1 f  =   15.439755 THz    97.010840 2PiTHz  515.014765 cm-1    63.853716 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.000000    0.000000    0.410501  
      2.003000  0.000000  2.003000     0.000000    0.000000   -0.272451  
      0.000000  2.003000  2.003000     0.000000    0.000000   -0.790776  
      2.003000  2.003000  2.003000     0.000000    0.000000    0.363062  
      0.000000  0.000000  0.000000     0.000000    0.000000    0.010895  

   2 f  =   14.403101 THz    90.497354 2PiTHz  480.435728 cm-1    59.566460 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.000000    0.128400    0.000000  
      2.003000  0.000000  2.003000     0.000000   -0.484223    0.000000  
      0.000000  2.003000  2.003000     0.000000    0.834506    0.000000  
      2.003000  2.003000  2.003000     0.000000   -0.228456    0.000000  
      0.000000  0.000000  0.000000     0.000000    0.021191    0.000000  

   3 f  =    9.980407 THz    62.708746 2PiTHz  332.910530 cm-1    41.275660 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000    -0.265607    0.000000    0.000000  
      2.003000  0.000000  2.003000    -0.563172    0.000000    0.000000  
      0.000000  2.003000  2.003000     0.766038    0.000000    0.000000  
      2.003000  2.003000  2.003000     0.143725    0.000000    0.000000  
      0.000000  0.000000  0.000000    -0.069418    0.000000    0.000000  

   4 f  =    9.424248 THz    59.214298 2PiTHz  314.359075 cm-1    38.975572 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.000000    0.000000   -0.339070  
      2.003000  0.000000  2.003000     0.000000    0.000000    0.776425  
      0.000000  2.003000  2.003000     0.000000    0.000000   -0.510194  
      2.003000  2.003000  2.003000     0.000000    0.000000   -0.144225  
      0.000000  0.000000  0.000000     0.000000    0.000000   -0.033118  

   5 f  =    6.907624 THz    43.401880 2PiTHz  230.413517 cm-1    28.567645 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.000000    0.307521    0.000000  
      2.003000  0.000000  2.003000     0.000000   -0.767570    0.000000  
      0.000000  2.003000  2.003000     0.000000   -0.532701    0.000000  
      2.003000  2.003000  2.003000     0.000000   -0.135043    0.000000  
      0.000000  0.000000  0.000000     0.000000    0.119416    0.000000  

   6 f  =    5.812542 THz    36.521276 2PiTHz  193.885510 cm-1    24.038748 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.000000    0.000000    0.273711  
      2.003000  0.000000  2.003000     0.000000    0.000000    0.330617  
      0.000000  2.003000  2.003000     0.000000    0.000000    0.235149  
      2.003000  2.003000  2.003000     0.000000    0.000000    0.472789  
      0.000000  0.000000  0.000000     0.000000    0.000000   -0.732769  

   7 f  =    5.010456 THz    31.481622 2PiTHz  167.130808 cm-1    20.721587 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.000000    0.370569    0.000000  
      2.003000  0.000000  2.003000     0.000000   -0.094861    0.000000  
      0.000000  2.003000  2.003000     0.000000    0.048415    0.000000  
      2.003000  2.003000  2.003000     0.000000    0.515184    0.000000  
      0.000000  0.000000  0.000000     0.000000   -0.765455    0.000000  

   8 f  =    4.076207 THz    25.611565 2PiTHz  135.967632 cm-1    16.857844 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.231970    0.000000    0.000000  
      2.003000  0.000000  2.003000     0.230327    0.000000    0.000000  
      0.000000  2.003000  2.003000     0.069857    0.000000    0.000000  
      2.003000  2.003000  2.003000     0.613178    0.000000    0.000000  
      0.000000  0.000000  0.000000    -0.715732    0.000000    0.000000  

   9 f  =    1.169851 THz     7.350394 2PiTHz   39.022044 cm-1     4.838119 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000    -0.179440    0.000000    0.000000  
      2.003000  0.000000  2.003000    -0.398238    0.000000    0.000000  
      0.000000  2.003000  2.003000    -0.317558    0.000000    0.000000  
      2.003000  2.003000  2.003000    -0.519386    0.000000    0.000000  
      0.000000  0.000000  0.000000    -0.662271    0.000000    0.000000  

  10 f/i=    1.213997 THz     7.627769 2PiTHz   40.494585 cm-1     5.020690 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.000000    0.000000    0.232919  
      2.003000  0.000000  2.003000     0.000000    0.000000    0.386143  
      0.000000  2.003000  2.003000     0.000000    0.000000    0.239573  
      2.003000  2.003000  2.003000     0.000000    0.000000    0.527859  
      0.000000  0.000000  0.000000     0.000000    0.000000    0.678684  

  11 f/i=    1.589322 THz     9.986006 2PiTHz   53.014082 cm-1     6.572911 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.000000    0.365880    0.000000  
      2.003000  0.000000  2.003000     0.000000    0.035667    0.000000  
      0.000000  2.003000  2.003000     0.000000    0.131682    0.000000  
      2.003000  2.003000  2.003000     0.000000    0.669647    0.000000  
      0.000000  0.000000  0.000000     0.000000    0.631738    0.000000  

  12 f/i=    5.231244 THz    32.868872 2PiTHz  174.495495 cm-1    21.634692 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.321368    0.000000    0.000000  
      2.003000  0.000000  2.003000    -0.683471    0.000000    0.000000  
      0.000000  2.003000  2.003000    -0.452513    0.000000    0.000000  
      2.003000  2.003000  2.003000     0.427094    0.000000    0.000000  
      0.000000  0.000000  0.000000     0.205943    0.000000    0.000000  

  13 f/i=   11.525150 THz    72.414655 2PiTHz  384.437620 cm-1    47.664207 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.000000    0.000000    0.766389  
      2.003000  0.000000  2.003000     0.000000    0.000000    0.254010  
      0.000000  2.003000  2.003000     0.000000    0.000000    0.041048  
      2.003000  2.003000  2.003000     0.000000    0.000000   -0.587555  
      0.000000  0.000000  0.000000     0.000000    0.000000    0.034953  

  14 f/i=   12.506382 THz    78.579915 2PiTHz  417.167986 cm-1    51.722257 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.860329    0.000000    0.000000  
      2.003000  0.000000  2.003000    -0.063727    0.000000    0.000000  
      0.000000  2.003000  2.003000     0.320460    0.000000    0.000000  
      2.003000  2.003000  2.003000    -0.388825    0.000000    0.000000  
      0.000000  0.000000  0.000000    -0.043508    0.000000    0.000000  

  15 f/i=   14.030104 THz    88.153744 2PiTHz  467.993886 cm-1    58.023868 meV
             X         Y         Z           dx          dy          dz
      2.003000  2.003000  0.000000     0.000000    0.785974    0.000000  
      2.003000  0.000000  2.003000     0.000000    0.407546    0.000000  
      0.000000  2.003000  2.003000     0.000000   -0.012030    0.000000  
      2.003000  2.003000  2.003000     0.000000   -0.464468    0.000000  
      0.000000  0.000000  0.000000     0.000000    0.016629    0.000000
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-10-10 17:03:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)
顶顶！希望有人关注
赞一下 回复此楼
2楼2010-10-11 13:14:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)
★
zhang668(金币+1):谢谢参与
楼主说的“差分方法”是什么方法？能否先说个大概？
楼主所列的输出数据中，“X Y Z dx dy dz”  具体是什么含义啊？
赞一下(1人) 回复此楼
先进使役材料计算（https://aidme.nimte.ac.cn）
3楼2010-10-11 13:42:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★
zhang668(金币+1):谢谢参与
呵呵，我做个G点的频率计算，但只关注频率值，没注意本征矢。看不懂这个结果，呵呵。

[ Last edited by stractor on 2010-10-11 at 21:42 ]
赞一下(1人) 回复此楼
4楼2010-10-11 21:39:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★
zhang668(金币+1):谢谢参与
引用回帖:
Originally posted by lfhuang at 2010-10-11 13:42:20:
楼主说的“差分方法”是什么方法？能否先说个大概？
楼主所列的输出数据中，“X Y Z dx dy dz”  具体是什么含义啊？

楼主所言有限差分法可能是因为手册中说"To calculate the Hessian matrix, finite differences are used”，不知道是不是该翻成“有限差分法”。
详见http://cms.mpi.univie.ac.at/vasp/vasp/IBRION_5_IBRION_6.html
赞一下(1人) 回复此楼
心之所向无惧无悔
5楼2010-10-12 09:36:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★ ★
zzy870720z(金币+2):谢谢指教 2010-10-12 14:19:46
你关掉对称性再试试，ISYM=0。如果不关掉对称性，这种方法移动原子位置后，所用的仍然是原来由对称性决定的的不可约K点
另外计算振动频率是不是完全弛豫之后再做这一步计算效率高些
赞一下(1人) 回复此楼
心之所向无惧无悔
6楼2010-10-12 09:40:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)
★
zzy870720z(金币+1):谢谢解答 2010-10-12 14:20:04
dx，dy，dy就是振动模式。
对称性和驰豫只影响频率的精确值和振幅，不会影响振动模式啊。所以很奇怪怎么振动模式的位移都是一个方向的。
赞一下(1人) 回复此楼
7楼2010-10-12 13:30:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
引用回帖:
Originally posted by zhang668 at 2010-10-12 13:30:40:
dx，dy，dy就是振动模式。
对称性和驰豫只影响频率的精确值和振幅，不会影响振动模式啊。所以很奇怪怎么振动模式的位移都是一个方向的。

怎么看出振动的位移都是沿一个方向的呢？能不能解释下结果中每一个参数的含义
赞一下 回复此楼
心之所向无惧无悔
8楼2010-10-12 14:34:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)
★
ym23(金币+1):欢迎常来交流 2010-10-12 16:04:57
每个模式位移都是沿dx，dy，dz中的一个，就是三个只有1列有值，就是沿一个方向振动。
赞一下(1人) 回复此楼
9楼2010-10-12 14:55:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

losenq

新虫 (小有名气)
★
zhang668(金币+1): 谢谢参与
设置ISIM=0

IBRION=6试试
赞一下 回复此楼
10楼2013-10-29 20:37:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhang668 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 博后网上办公系统维护 +13 取名字烦人 2024-06-25 20/1000 2024-06-26 18:32 by Breezy恋
[基金申请] 国基在研影响申请结果吗 +10 WOWO159357 2024-06-26 10/500 2024-06-26 18:19 by 德尚中行
[考研] 刚当完兵回来想考研 +7 五官在线 2024-06-23 18/900 2024-06-26 17:05 by 1158057902
[论文投稿] 投稿求助 +5 平凡的日子 2024-06-19 6/300 2024-06-26 16:38 by 大大熊呀
[基金申请] 博后面上今天下午会公布吗？大家有无消息？ +9 地理学1995 2024-06-24 13/650 2024-06-26 16:02 by kyukitu
[基金申请] 博后面上和特助今天出吗？ +41 逗您玩 2024-06-21 78/3900 2024-06-26 16:00 by HAPPY_0225
[基金申请] 自然基金变化？？ +14 今晚推荐22 2024-06-21 24/1200 2024-06-26 15:47 by 梦渺岚烟
[育儿交流] 暑假孩子一个人在家可以干什么 +4 唠叨的映 2024-06-25 5/250 2024-06-26 15:35 by 361997417
[基金申请] 青年基金E02口青基去年几个函评专家？ +6 他山攻玉之石 2024-06-25 9/450 2024-06-26 15:09 by 他山攻玉之石
[基金申请] 要持续整治滥发“帽子”、“牌子”之风 +6 babu2015 2024-06-25 6/300 2024-06-26 14:52 by felicity6056
[基金申请] 刚刚收到科研之友邮件 +30 olivermiaoer 2024-06-19 44/2200 2024-06-26 13:11 by 芦苇薇薇～
[公派出国] CSC博士联培对将来就业有用吗 +4 也就这样 2024-06-22 4/200 2024-06-26 10:46 by lxs7405
[基金申请] 专科院校申请国自然 +14 孤单12站 2024-06-24 20/1000 2024-06-26 09:09 by chong2024
[考博] 申博好难 +6 自强不息a?a 2024-06-24 6/300 2024-06-25 23:02 by 考研吒儿
[基金申请] 这样的说辞是上会了吗 +12 学员d3zYCz 2024-06-24 15/750 2024-06-25 22:05 by 3115321
[考博] 申请25博士，可以提前进组做科研助理 +6 逐梦途中w 2024-06-22 8/400 2024-06-24 21:49 by keyaner23
[基金申请] 说博后基金7月出的真打电话了吗？ +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[基金申请] 青年和面上，哪个上会难度更大 +12 今晚推荐22 2024-06-21 18/900 2024-06-24 11:08 by 半简体
[基金申请] F03青年基金函评结果 +5 暨阳一只柴 2024-06-19 6/300 2024-06-23 14:30 by adsqsj
[基金申请] 工材口青年基金大概什么样能上会？ +15 今晚推荐22 2024-06-20 21/1050 2024-06-22 23:04 by qbn0326
信息提示
关闭
请填处理意见
关闭