| 2 | 1/1 | 返回列表 |
| 查看: 1126 | 回復(fù): 1 | ||
111hzs銅蟲(chóng) (小有名氣)
|
[求助]
不同方法計(jì)算得到Gamma點(diǎn)聲子不宜致 已有1人參與
|
|
大家好,我正在用PWscf計(jì)算某種材料的Gamma點(diǎn)聲子,通過(guò)單點(diǎn)計(jì)算得到的結(jié)果如下: freq ( 1) =0.00 [cm-1] freq ( 2) =0.00 [cm-1] freq ( 3) =0.00 [cm-1] freq ( 4) =176.91 [cm-1] freq ( 5) =176.91 [cm-1] freq ( 6) =269.40 [cm-1] freq ( 7) =269.40 [cm-1] freq ( 8) =338.76 [cm-1] freq ( 9) =340.88 [cm-1] 得到的結(jié)果和別人的工作一致的。 我進(jìn)一步在4*4*1的q-mesh網(wǎng)格上計(jì)算聲子譜,在完成了q2r.x和matdyn.x計(jì)算之后,"matdyn.modes"文件中得到的聲子譜數(shù)據(jù)顯示Gamma點(diǎn)的聲子頻率如下: freq ( 1) = -0.000007 [cm-1] freq ( 2) =-0.000006 [cm-1] freq ( 3) = -0.000001 [cm-1] freq ( 4) = 176.907415 [cm-1] freq ( 5) =269.403108 [cm-1] freq ( 6) =269.403108 [cm-1] freq ( 7) = 298.948294 [cm-1] freq ( 8) =338.760427 [cm-1] freq ( 9) =340.877820 [cm-1] 和我進(jìn)行單點(diǎn)計(jì)算得到的Gamma點(diǎn)聲子頻率有比較大的出入。這是為什么呢? 輸入文件中的一下參數(shù)的設(shè)置如下:ecutwfc = 90, ecutrho = 900, conv_thr = 1.0d-8, mixing_beta = 0.7, etot_conv_thr 和 forc_conv_thr 都用得默認(rèn)值K-mesh為17 17 2 1 1 1,在聲子計(jì)算過(guò)程中求和規(guī)則控制參數(shù)ASR設(shè)置為'simple'。 在計(jì)算開(kāi)始之前,晶體結(jié)構(gòu)被優(yōu)化到最佳狀態(tài),也完成了K-mesh和cut-off收斂性的測(cè)試。 我以為是關(guān)于q-mesh不收斂,但是我重新計(jì)算了q-mesh為4*4*2和8*8*1情況下的聲子譜,得到的結(jié)果和q-mesh為4*4*1的結(jié)果一樣。 為什么通過(guò)單點(diǎn)計(jì)算得到的Gamma點(diǎn)的聲子頻率和通過(guò)q2r.x和matdyn.x計(jì)算得到的Gamma點(diǎn)的聲子頻率有那么大的出入? |
金蟲(chóng) (小有名氣)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 279分求調(diào)劑 一志愿211 +3 | chaojifeixia 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 3/150 |
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級(jí)已過(guò),有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 材料專(zhuān)碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚(yú) 2026-03-18 | 6/300 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專(zhuān)業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|