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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (4133)
>导师招生 (770)
>虫友互识 (627)
>文献求助 (411)
>考研 (375)
>招聘信息布告栏 (211)
>硕博家园 (208)
>休闲灌水 (176)
>考博 (132)
>论文投稿 (117)
>博后之家 (88)
>教师之家 (49)
>公派出国 (45)
>SciFinder/Reaxys (41)
>基金申请 (36)
>找工作 (36)

返回列表

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

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

回复此楼

» 猜你喜欢

求助调剂，跨调已经有20人回复
291 求调剂已经有40人回复
300求调剂已经有5人回复
急需调剂已经有3人回复
22408 312求调剂已经有21人回复
恳请有学校收留已经有4人回复
297，工科调剂? 已经有9人回复
291求调剂已经有5人回复
085404 22408 309分求调剂已经有10人回复
296求调剂已经有13人回复

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

求助用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:结构、力学和

引用回帖:

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

查看全部 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 个回答

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

[考研] 材料工程281还有调剂机会吗 +43	xaw. 2026-04-11	44/2200	2026-04-15 12:46 by 西北望—风沙
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +26	慕绝cc 2026-04-09	30/1500	2026-04-14 18:50 by 蔡苏阳
[考研] 药学305求调剂 +10	玛卡巴卡boom 2026-04-10	10/500	2026-04-14 15:55 by zs92450
[考研] 290求调剂 +21	luoziheng 2026-04-10	23/1150	2026-04-14 15:49 by zs92450
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11	不会吃肉 2026-04-12	11/550	2026-04-13 21:59 by bljnqdcc
[考研] 一志愿西交机械专硕求调剂 +9	求上岸的小王 2026-04-10	9/450	2026-04-13 16:08 by jiangguiquan11
[考研] 一志愿双非085400电子信息344 求调剂，对材料和化学方向也感兴趣 +12	无情的小羊 2026-04-09	13/650	2026-04-13 14:17 by 张zhihao
[考研] 344 材料专业求调剂211 无地域要求 +8	hualkop 2026-04-11	8/400	2026-04-12 22:24 by fqwang
[考研] 本科南方医科大学一志愿985 药学学硕284分求调剂 +5	弱水听文 2026-04-09	5/250	2026-04-12 13:16 by shengxi123
[考研] 一志愿西北工业大学289 085602 +33	yang婷 2026-04-10	34/1700	2026-04-12 08:11 by Art1977
[考研] 267求调剂 +8	再忙也要吃饭啊 2026-04-09	8/400	2026-04-11 21:42 by cfdbai
[考研] 284求调剂 +11	archer.. 2026-04-09	12/600	2026-04-11 20:23 by 蓝云思雨
[考研] 本人女孩 +7	吼吼， 2026-04-10	9/450	2026-04-11 14:45 by ACS Nano——
[考研] 机械专硕270求调剂，接受跨专业 +12	老师看看我吧aba 2026-04-09	14/700	2026-04-11 10:21 by laoshidan
[考研] 297求调剂 +9	Kwgyz 2026-04-09	9/450	2026-04-11 10:09 by zhq0425
[考研] 22408 327分求调剂 +4	韵风kon 2026-04-10	4/200	2026-04-11 09:51 by 猪会飞
[考研] 337求调剂 +4	研s. 2026-04-10	4/200	2026-04-11 08:57 by zhq0425
[考研] 342电子信息专硕求调剂 +9	你让我怎么荔枝 2026-04-10	10/500	2026-04-11 08:33 by zhq0425
[考研] 265求调剂 +12	风说她早忘了 2026-04-10	13/650	2026-04-10 18:56 by chemisry
[考研] 青岛科技大学材料学院，环境学院调剂补录4月10日以前都可以 +3	1青科大。 2026-04-09	5/250	2026-04-10 09:58 by 翩翩一书生