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

登录注册

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

回复此楼

» 猜你喜欢

求b区院校调剂已经有3人回复
材料与化工考研调剂已经有8人回复
281求调剂已经有5人回复
0854电子信息求调剂 324 已经有4人回复
材料调剂已经有3人回复
一志愿吉林大学材料与化工303分求调剂已经有4人回复
294求调剂材料与化工专硕已经有17人回复
271求调剂已经有4人回复
求调剂323材料与化工已经有4人回复
材料学硕，求调剂已经有9人回复

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

求助用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 ：电子结构

★
ym23(金币+1):欢迎常来交流 2010-10-12 16:04:57

每个模式位移都是沿dx，dy，dz中的一个，就是三个只有1列有值，就是沿一个方向振动。

赞一下(1人)

回复此楼

9楼2010-10-12 14:55:27

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

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

[考研] 一志愿西北大学，070300化学学硕，总分287，双非一本，求调剂。 +4	晨昏线与星海 2026-03-20	4/200	2026-03-25 10:16 by allen-yin
[考研] 303求调剂 +3	蓝山月 2026-03-25	3/150	2026-03-25 10:04 by 雾散后相遇lc
[考研] 318求调剂 +3	plum李子 2026-03-23	3/150	2026-03-25 09:42 by 雾散后相遇lc
[有机交流] 有机合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 材料专硕找调剂 +5	哈哈哈吼吼吼哈 2026-03-23	5/250	2026-03-24 19:07 by 了了了了。。
[考研] 化工专硕求调剂 +3	question挽风 2026-03-24	3/150	2026-03-24 18:48 by jhhcooi
[考研] 0854 考研调剂招生了！AI 方向 +5	pk3725069 2026-03-19	17/850	2026-03-24 17:30 by zhouxuan..
[考研] 300求调剂，材料科学英一数二 +5	leaflight 2026-03-24	5/250	2026-03-24 16:25 by laoshidan
[考研] 341求调剂(一志愿湖南大学070300) +5	番茄头--- 2026-03-22	6/300	2026-03-23 23:45 by Txy@872106
[考研] 材料专硕英一数二306 +8	z1z2z3879 2026-03-18	8/400	2026-03-23 20:49 by baobaoye
[考研] 263求调剂 +6	yqdszhdap－ 2026-03-22	9/450	2026-03-23 12:57 by yqdszhdap－
[考研] 269专硕求调剂 +6	金恩贝 2026-03-21	6/300	2026-03-22 14:31 by ColorlessPI
[基金申请] 山东省面上项目限额评审 +4	石瑞0426 2026-03-19	4/200	2026-03-22 08:50 by Wei_ren
[考研] 材料学硕301分求调剂 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 0703化学297求调剂 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7	墨墨漠 2026-03-20	7/350	2026-03-21 16:36 by barlinike
[考研] 0805材料320求调剂 +3	深海物语 2026-03-20	3/150	2026-03-21 15:46 by 无际的草原
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 考研调剂求学校推荐 +3	伯乐29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky