| 5 | 1/1 | 返回列表 |
| 查看: 2762 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
zhang668木蟲 (正式寫手)
|
[交流]
【求助】有人用有限差分方法計(jì)算過G點(diǎn)聲子譜嗎?
|
||
|
如題,我設(shè)置IBRION=5,計(jì)算了立方BaTiO3的G點(diǎn)聲子譜,結(jié)果發(fā)現(xià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è)方向上的,而且沒有發(fā)現(xiàn)平移的聲子模式。所以請(qǐng)問這個(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 |
木蟲 (著名寫手)

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +5 | 材料學(xué)257求調(diào)劑 2026-03-23 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +5 |  ̄^ ̄゜汗 2026-03-19 | 7/350 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|