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

登录注册

应《网络安全法》要求，自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的回帖

zhang668

木虫 (正式写手)

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

★
zzy870720z(金币+1):谢谢解答 2010-10-12 14:20:04

dx，dy，dy就是振动模式。
对称性和驰豫只影响频率的精确值和振幅，不会影响振动模式啊。所以很奇怪怎么振动模式的位移都是一个方向的。

赞一下(1人)

回复此楼

7楼2010-10-12 13:30:40

已阅回复此楼关注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的回帖

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

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

查看全部 11 个回答

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

[考研] 材料学求调剂 +6	Stella_Yao 2026-03-20	6/300	2026-03-25 00:37 by baoball
[考研] 08工学调剂 +12	用户573181 2026-03-20	17/850	2026-03-24 21:45 by 番茄炖了西红柿
[考研] 384求调剂 +3	子系博 2026-03-22	6/300	2026-03-23 21:45 by 子系博
[考研] 一志愿重庆大学085700资源与环境，总分308求调剂 +7	墨墨漠 2026-03-23	8/400	2026-03-23 20:36 by Creta
[考研] 一志愿武理材料工程348求调剂 +6	￣^￣゜汗 2026-03-19	9/450	2026-03-23 19:53 by pswait
[考研] 291 求调剂 +4	化工2026届毕业� 2026-03-21	5/250	2026-03-23 16:46 by 化工2026届毕业�
[考研] 0854电子信息求调剂 324 +3	Promise-jyl 2026-03-23	3/150	2026-03-23 13:43 by wangkm
[考研] 材料与化工085600，总分304，本科有两篇sci参与，求调剂 +4	幸运的酱酱 2026-03-22	5/250	2026-03-22 20:15 by edmund7
[考研] 求调剂 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 求调剂 +7	Auroracx 2026-03-22	7/350	2026-03-22 12:38 by 素颜倾城1988
[考研] 考研调剂 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 无际的草原
[基金申请] 学校已经提交到NSFC，还能修改吗？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 279求调剂 +5	红衣隐官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 330求调剂0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏击518
[考研] 一志愿西北大学，070300化学学硕，总分287，双非一本，求调剂。 +3	晨昏线与星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 一志愿西南交大，求调剂 +5	材化逐梦人 2026-03-18	5/250	2026-03-21 00:26 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 317求调剂 +5	申子申申 2026-03-19	9/450	2026-03-20 22:26 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*