| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2761 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
zhang668木蟲(chóng) (正式寫(xiě)手)
|
[交流]
【求助】有人用有限差分方法計(jì)算過(guò)G點(diǎn)聲子譜嗎?
|
||
|
如題,我設(shè)置IBRION=5,計(jì)算了立方BaTiO3的G點(diǎn)聲子譜,結(jié)果發(fā)現(xiàn)了些問(wèn)題: 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個(gè)特征向量都是沿一個(gè)方向的,這個(gè)明顯不是簡(jiǎn)諧振動(dòng)模式,比如氧八面體扭曲模式O的移動(dòng)就不可能是在一個(gè)方向上的,而且沒(méi)有發(fā)現(xiàn)平移的聲子模式。所以請(qǐng)問(wèn)這個(gè)跟聲子有什么區(qū)別?下面是結(jié)果: 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 |
金蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)

至尊木蟲(chóng) (著名寫(xiě)手)
|
樓主所言有限差分法可能是因?yàn)槭謨?cè)中說(shuō)"To calculate the Hessian matrix, finite differences are used”,不知道是不是該翻成“有限差分法”。 詳見(jiàn)http://cms.mpi.univie.ac.at/vasp/vasp/IBRION_5_IBRION_6.html |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +6 | 生物工程調(diào)劑 2026-03-17 | 10/500 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專(zhuān)碩324分,求調(diào)劑 +3 | jiajunX 2026-03-22 | 3/150 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過(guò) 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)(081702)315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
[考研] 321求調(diào)劑 +9 | 何潤(rùn)采123 2026-03-18 | 11/550 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|