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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】有人用有限差分方法计算过G点声子谱吗？
111/2返回列表12下一页
查看: 2764  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿武理085500机械专业总分300求调剂 +3 an10101 2026-03-24 7/350 2026-03-25 00:00 by 山鬼0-
[考研] 08工学调剂 +12 用户573181 2026-03-20 17/850 2026-03-24 21:45 by 番茄炖了西红柿
[考研] 07化学280分求调剂 +6 722865 2026-03-23 6/300 2026-03-24 20:58 by allen-yin
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 0703化学求调剂 +4 奶油草莓. 2026-03-22 5/250 2026-03-23 19:37 by pswait
[考研] 考研化学308分求调剂 +7 你好明天你好 2026-03-23 8/400 2026-03-23 18:39 by macy2011
[考研] 上海电力大学材料防护与新材料重点实验室招收调剂研究生（材料、化学、电化学，环境） +3 我爱学电池 2026-03-23 3/150 2026-03-23 17:16 by AZMK
[考研] 0854电子信息求调剂 324 +3 Promise-jyl 2026-03-23 3/150 2026-03-23 13:43 by wangkm
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿华中科技大学，080502，354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 295材料求调剂，一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 0856调剂，是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
信息提示
关闭
请填处理意见
关闭