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

登录注册

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

回复此楼

» 猜你喜欢

药学求调剂已经有14人回复
327求调剂已经有27人回复
0854求调剂已经有11人回复
急需调剂已经有5人回复
300求调剂已经有7人回复
271求调剂已经有33人回复
273求调剂已经有8人回复
材料相关专业344求调剂双非工科学校或课题组已经有23人回复
297，工科调剂? 已经有3人回复
求调剂学校已经有16人回复

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

求助用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的回帖

sunyang1988

至尊木虫 (著名写手)

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

★ ★
zzy870720z(金币+2):谢谢指教 2010-10-12 14:19:46

你关掉对称性再试试，ISYM=0。如果不关掉对称性，这种方法移动原子位置后，所用的仍然是原来由对称性决定的的不可约K点
另外计算振动频率是不是完全弛豫之后再做这一步计算效率高些

赞一下(1人)

回复此楼

心之所向无惧无悔

6楼2010-10-12 09:40:47

已阅回复此楼关注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 (小学生)
金币: 3585.1
散金: 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 个回答

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

[考研] 273求调剂 +6	白居不易. 2026-04-09	8/400	2026-04-15 22:02 by wooluyong
[考研] 297，工科调剂? +3	河南农业大学-能 2026-04-14	3/150	2026-04-15 22:02 by noqvsozv
[考研] 284求调剂 +21	让我上岸吧阿西 2026-04-09	21/1050	2026-04-15 21:53 by noqvsozv
[考研] 22专硕求调剂 +9	haoyun上岸 2026-04-11	11/550	2026-04-15 14:47 by elainzaizai
[考研] 085500求调剂材料 +11	易11122 2026-04-09	11/550	2026-04-14 17:59 by lhj2009
[考研] 各位老师好，求调剂，本科211，一志愿天津大学生物与医药学硕，差两名录取。 +11	路六六jjj 2026-04-13	11/550	2026-04-14 16:01 by zs92450
[考研] 药学305求调剂 +10	玛卡巴卡boom 2026-04-10	10/500	2026-04-14 15:55 by zs92450
[考研] 366求调剂 +11	不知名的小卅 2026-04-11	11/550	2026-04-14 15:50 by zs92450
[基金申请] RY：中国产出的科学垃圾论文，绝对数量和比例都世界第一 +6	zju2000 2026-04-14	17/850	2026-04-14 14:34 by jurkat.1640
[考研] 材料085601调剂 +32	何润采123 2026-04-10	34/1700	2026-04-14 08:47 by 木木mumu～
[考研] 本科西工大 324求调剂 +5	wysyjs25 2026-04-10	5/250	2026-04-13 23:08 by pies112
[考研] 346分，工科0854求调剂，专硕 +6	moser233 2026-04-12	7/350	2026-04-12 22:11 by fqwang
[考研] 一志愿厦大0856，306求调剂 +15	Bblinging 2026-04-11	15/750	2026-04-11 22:53 by 314126402
[考研] 352 求调剂 +6	yzion 2026-04-11	8/400	2026-04-11 16:24 by 明月此时有
[考研] 22408 327分求调剂 +4	韵风kon 2026-04-10	4/200	2026-04-11 09:51 by 猪会飞
[考研] 311求调剂 +13	xyp想读书 2026-04-10	14/700	2026-04-11 09:41 by 猪会飞
[考研] 281求调剂 +11	觉得好的吧 2026-04-10	11/550	2026-04-11 09:35 by 逆水乘风
[考研] 22408 366分，本科211，一志愿西工大 +4	Rubt 2026-04-09	4/200	2026-04-10 19:51 by chemisry
[考研] 314求调剂 +14	weltZeng 2026-04-09	14/700	2026-04-09 23:14 by wolf97
[考研] 已调剂 +18	柴郡猫_ 2026-04-09	19/950	2026-04-09 22:10 by 柴郡猫_