| 查看: 2314 | 回復(fù): 24 | ||||
hanyanli0475銅蟲 (小有名氣)
|
[求助]
聲子計(jì)算又負(fù)頻的出現(xiàn)!!
|
|
今天計(jì)算氮化銦的聲子,輸出后低頻有-26cm-1的負(fù)頻,而且其他低頻的走勢(shì)和文獻(xiàn)也不太一樣,我把輸入文件粘出來(lái),請(qǐng)大家?guī)臀艺艺义e(cuò)誤唄,看看有沒有需要修改的地方,謝謝了!! 自洽計(jì)算: &CONTROL calculation = 'scf' , restart_mode = 'from_scratch' , outdir = '/home/zhao/10.21/tmp/' , pseudo_dir = '/opt/espresso-4.2.1/pseudo/' , prefix = 'InN' , tprnfor= .true. , tstress= .true. , / &SYSTEM ibrav = 4,celldm(1) = 6.893,celldm(3) = 1.605,nat = 4,ntyp = 2, ecutwfc = 45.0 ,ecutrho = 180 , / &ELECTRONS conv_thr = 1.0d-8 , mixing_beta = 0.7 , / ATOMIC_SPECIES In 114.820 In.pbe-d-rrkjus.UPF N 14.007 N.pbe-rrkjus.UPF ATOMIC_POSITIONS crystal In 0.333333487 0.666666513 0.003058286 1 1 1 In 0.666666513 0.333333487 0.503058286 1 1 1 N 0.333333439 0.666666561 0.381941714 1 1 1 N 0.666666561 0.333333439 0.881941714 1 1 1 K_POINTS automatic 4 4 4 1 1 1 聲子計(jì)算 phonon for InN &inputph outdir ='/home/zhao/10.21/tmp/', tr2_ph =1.0d-12, prefix ='InN', fildvscf ='InNdv', amass(1) =114.8200, amass(2) =14.00700, fildyn ='InN.dyn', elph =.false., trans =.true., ldisp =.true., nq1=4,nq2=4,nq3=4 / |

專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |

銅蟲 (小有名氣)

金蟲 (著名寫手)
銅蟲 (小有名氣)

金蟲 (著名寫手)
|
測(cè)試很簡(jiǎn)單,就是把截?cái)鄰男〉竭_(dá)取不同數(shù)值進(jìn)行自洽場(chǎng)迭代計(jì)算,看不同截?cái)嗟臅r(shí)候的能量收斂精度。如果基本前后兩步的數(shù)值基本一致,則可認(rèn)為收斂。能量或者壓強(qiáng)的收斂與截?cái)、K點(diǎn)選擇等有關(guān)系,你應(yīng)該不同的參數(shù)分別測(cè)試。 PS:很多時(shí)候我們計(jì)算的時(shí)候都是不管收不收斂,直接給定數(shù)值就進(jìn)行計(jì)算,這是完全不對(duì)的,不是科學(xué)的態(tài)度。我們應(yīng)該一步一步的來(lái),多看看別人給你的步驟和例子,里面都是有的。 自己想當(dāng)然的算,然后靠別人幫你指出錯(cuò)誤,尤其是很簡(jiǎn)單又討論過無(wú)數(shù)次的錯(cuò)誤,這對(duì)自己而言雖然輕松了,但少了一個(gè)自己探索的過程。而且簡(jiǎn)單重復(fù)的回答,也浪費(fèi)了別人的時(shí)間。歸根結(jié)底,是自身沒有足夠的學(xué)習(xí)和調(diào)研造成許多根本上的卻很無(wú)謂的錯(cuò)誤。 說(shuō)多了。 |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |
|
你這么說(shuō)也不管用啊,畢竟測(cè)試也是個(gè)比較耗時(shí)間精力的過程,并且還可能證明了之前看起來(lái)跟實(shí)驗(yàn)非!胺稀钡慕Y(jié)果有問題了~ 事實(shí)上去年我在實(shí)驗(yàn)室里做組會(huì)報(bào)告時(shí),做的就是如何測(cè)試收斂的問題,被組里幾個(gè)也做計(jì)算的人說(shuō)是很無(wú)聊和沒有特殊的東西的報(bào)告~ 結(jié)果,有人的計(jì)算就是低截?cái)嗟蚄點(diǎn)得到了跟實(shí)驗(yàn)很相近的結(jié)果,而我在測(cè)試他的模型時(shí),加大了K點(diǎn)和截?cái),于是就跟?shí)驗(yàn)值差別非常大了~ 不過,由于老板不懂計(jì)算,也沒有這樣的科學(xué)思想,所以,覺得我加大截?cái)嗪蚄點(diǎn)的做法是沒選對(duì)參數(shù),而那個(gè)哥們做的結(jié)果符合實(shí)驗(yàn),是好結(jié)果…… 總的來(lái)說(shuō),工科的急功近利的思想啊~ |

銅蟲 (小有名氣)

鐵桿木蟲 (正式寫手)
鐵桿木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +13 | rare12345 2026-03-18 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +8 | 何潤(rùn)采123 2026-03-18 | 10/500 |
|
|
[考研] 288求調(diào)劑 +11 | 于海海海海 2026-03-19 | 11/550 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +8 | 妮妮ninicgb 2026-03-15 | 12/600 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|