| 查看: 2059 | 回復(fù): 11 | ||||
gemucai木蟲 (正式寫手)
|
[求助]
用自己生成的贗勢算能帶效果很好,算聲子卻一團糟
|
| 我用的是PWSCF,能帶計算和聲子計算都比較熟悉了。但最近我自己生成贗勢之后,算出的能帶與別人文獻(xiàn)比較幾乎無差別,但gamma點的聲子頻率卻都是負(fù)的幾千1/cm。我做過了優(yōu)化,力和壓力都幾乎為0了,但還是有如此大的虛頻。為何電子和聲子計算會差別如此之大?有可能是我的贗勢哪里沒做好才導(dǎo)致了DFPT過程的異常? |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗: +174 |
| 啥的贗勢?就是你在郵件列表里說的那個C的全電子贗勢么?一般能帶計算不過是total energy的計算,對于體系的弛豫與否不是那么敏感~所以測試贗勢時一般都是至少測試晶格常數(shù)、彈性常數(shù)的~至于聲子,那是相當(dāng)?shù)拿舾衺 |

木蟲 (正式寫手)
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗: +174 |
|
Mo的贗勢很不好生成,尤其考慮到半芯態(tài)時。 這里有兩個從CPMD轉(zhuǎn)來的LDA模守恒贗勢,06的是沒有半芯態(tài),14的是有半芯態(tài)的。我近來在做其他的計算,還沒時間測試這兩個贗勢的效果和差別,如果有興趣,可以測試一下~ 順便,如果測試了,還希望早點發(fā)上測試結(jié)果圖來,也有利于其他使用QE的人~免得贗勢不好讓太多人浪費了時間,或者免得大家都得自己跑測試~ |

金蟲 (著名寫手)
木蟲 (正式寫手)
木蟲 (正式寫手)
|
那兩個贗勢好像有問題。我測試4p5 5s1組態(tài)的時候,ld1.x報錯 WARNING! Expected number of nodes: 0= 3- 2- 1, number of nodes found: 3. Setting wfc to zero for this iteration. (This warning will only be printed once per wavefunction) WARNING! Expected number of nodes: 0= 1- 0- 1, number of nodes found: 3. Setting wfc to zero for this iteration. (This warning will only be printed once per wavefunction) Message from routine run_pseudo: Errors in PS-KS equations WARNING! Expected number of nodes: 0= 2- 1- 1, number of nodes found: 3. Setting wfc to zero for this iteration. (This warning will only be printed once per wavefunction) 相反,對于Mo的陽離子組態(tài)卻誤差非常。 n l nl e AE (Ry) e PS (Ry) De AE-PS (Ry) 3 2 4D 1( 5.00) -0.78364 -0.78177 -0.00187 1 0 5S 1( 0.00) -0.74605 -0.74501 -0.00105 2 1 5P 1( 0.00) -0.44623 -0.44648 0.00024 ----------------------------------------------------------------------------------------- n l nl e AE (Ry) e PS (Ry) De AE-PS (Ry) 3 2 4D 1( 4.00) -1.54772 -1.55039 0.00267 1 0 5S 1( 0.00) -1.32339 -1.31887 -0.00452 2 1 5P 1( 0.00) -0.94261 -0.94129 -0.00132 至于放在晶體中的情形還沒有測試。 |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗: +174 |

木蟲 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +4 | 阿斯蒂芬2004 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
[考研] 求調(diào)劑 +5 | 蘆lty 2026-03-25 | 6/300 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績338 +4 | YYYYYNNNNN 2026-03-26 | 4/200 |
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
[考研] 材料277求調(diào)劑 +5 | min3 2026-03-24 | 5/250 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
[考研] 315分,誠求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|