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

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

金蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
|
那兩個(gè)贗勢(shì)好像有問(wèn)題。我測(cè)試4p5 5s1組態(tài)的時(shí)候,ld1.x報(bào)錯(cuò) 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) 相反,對(duì)于Mo的陽(yáng)離子組態(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 至于放在晶體中的情形還沒(méi)有測(cè)試。 |
專家顧問(wèn) (著名寫(xiě)手)
![]() |
專家經(jīng)驗(yàn): +174 |

木蟲(chóng) (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過(guò) 2026-03-27 | 6/300 |
|
|
[考研] 求調(diào)劑 +8 | 張zz111 2026-03-27 | 9/450 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +3 | IceGuo 2026-03-26 | 4/200 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚(yú) 2026-03-26 | 3/150 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚(yú)567 2026-03-23 | 4/200 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[有機(jī)交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|