版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhang668

木虫 (正式写手)

1ST强帖: 3
应助: 11 (小学生)
贵宾: 0.061
金币: 1735
散金: 1294
红花: 11
帖子: 843
在线: 324.9小时
虫号: 733242
注册: 2009-03-28
性别: GG
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】有人用有限差分方法计算过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

回复此楼

» 猜你喜欢

生物学学硕求调剂已经有10人回复
上海电力大学材料防护与新材料重点实验室招收调剂研究生（材料、化学、电化学，环境）已经有4人回复
材料学求调剂已经有6人回复
303求调剂已经有5人回复
一志愿武理085500机械专业总分300求调剂已经有7人回复
考研调剂已经有4人回复
281求调剂已经有4人回复
0805 316求调剂已经有6人回复
085601求调剂总分293英一数二已经有3人回复
08工学调剂已经有17人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助用MS计算声子谱出错的问题已经有15人回复
质谱硕士毕业找工作难吗已经有32人回复
“声子液体”优化热电性能已经有9人回复
R语言图书-求解差分方程和并行计算！已经有43人回复
有关矢量场的推导计算过程，求高手提点！已经有4人回复
今天用CASTEP计算声子谱，出现下列提示，有能看懂的大侠帮忙解释下已经有11人回复
关于用超胞方法计算声子的一些经验已经有21人回复
用siesta画声子谱的问题已经有12人回复
MS计算晶体声子谱的时候，横坐标有负能量值，为什么？已经有11人回复
castep计算Fe2O3声子谱失败，原因求教~~ 已经有10人回复
如何分析声子散射谱在G点的LO, TO分裂已经有8人回复
【求助】castep计算声子谱，为什么不能进行振动分析已经有9人回复
【求助】请教CASTEP计算声子谱后出现虚频已经有6人回复
【求助】关于PWSCF介电常数和声子谱的计算已经有7人回复

1楼 2010-10-10 17:03:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

应助: 39 (小学生)
金币: 3577.6
散金: 2141
红花: 37
帖子: 1432
在线: 681.4小时
虫号: 505824
注册: 2008-02-17
性别: GG
专业: 金属材料表面科学与工程

★
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的回帖

查看全部 11 个回答

zhang668

木虫 (正式写手)

1ST强帖: 3
应助: 11 (小学生)
贵宾: 0.061
金币: 1735
散金: 1294
红花: 11
帖子: 843
在线: 324.9小时
虫号: 733242
注册: 2009-03-28
性别: GG
专业: 凝聚态物性 II ：电子结构

顶顶！希望有人关注

赞一下

回复此楼

2楼2010-10-11 13:14:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

stractor

金虫 (著名写手)

1ST强帖: 6
应助: 91 (初中生)
金币: 1487.6
散金: 666
红花: 30
沙发: 1
帖子: 2224
在线: 426.8小时
虫号: 394087
注册: 2007-06-06
专业: 半导体材料

★
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

至尊木虫 (著名写手)

1ST强帖: 1
应助: 24 (小学生)
贵宾: 0.984
金币: 9397.7
散金: 1593
红花: 23
沙发: 1
帖子: 1744
在线: 1255.8小时
虫号: 965316
注册: 2010-03-09
性别: GG
专业: 凝聚态物性I:结构、力学和

★
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的回帖

查看全部 11 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 一志愿武理085500机械专业总分300求调剂 +3	an10101 2026-03-24	7/350	2026-03-25 00:00 by 山鬼0-
[考研] 281求调剂 +3	Koxui 2026-03-24	4/200	2026-03-24 22:59 by barlinike
[考研] 340求调剂 +4	话梅糖111 2026-03-24	4/200	2026-03-24 21:41 by yuyu98412
[考研] 一志愿河北工业大学0817化工278分求调剂 +7	jhybd 2026-03-23	12/600	2026-03-24 09:03 by jhybd
[考研] 269求调剂 +4	我想读研11 2026-03-23	4/200	2026-03-23 21:25 by pswait
[考研] 一志愿重庆大学085700资源与环境，总分308求调剂 +7	墨墨漠 2026-03-23	8/400	2026-03-23 20:36 by Creta
[考研] 350求调剂 +6	weudhdk 2026-03-19	6/300	2026-03-23 15:47 by tangyuan0840221
[考研] 311求调剂 +3	26研0 2026-03-20	3/150	2026-03-22 14:46 by ColorlessPI
[考研] 材料学硕301分求调剂 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 一志愿东华大学控制学硕320求调剂 +3	Grand777 2026-03-21	3/150	2026-03-21 19:23 by 简之-
[考研] 296求调剂 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 学员8dgXkO
[考研] 求调剂 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 085601调剂 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 304求调剂 +6	曼殊2266 2026-03-18	6/300	2026-03-21 00:32 by JourneyLucky
[考研] 南京大学化学376求调剂 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 329求调剂 +9	想上学吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] 一志愿西南交通专硕材料355 本科双非求调剂 +5	西南交通专材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 261求B区调剂，科研经历丰富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky
[考研] 材料学硕318求调剂 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad