| 查看: 2313 | 回復(fù): 24 | ||||
hanyanli0475銅蟲 (小有名氣)
|
[求助]
聲子計(jì)算又負(fù)頻的出現(xiàn)!!
|
|
今天計(jì)算氮化銦的聲子,輸出后低頻有-26cm-1的負(fù)頻,而且其他低頻的走勢(shì)和文獻(xiàn)也不太一樣,我把輸入文件粘出來,請(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)唷點(diǎn)選擇等有關(guān)系,你應(yīng)該不同的參數(shù)分別測(cè)試。 PS:很多時(shí)候我們計(jì)算的時(shí)候都是不管收不收斂,直接給定數(shù)值就進(jìn)行計(jì)算,這是完全不對(duì)的,不是科學(xué)的態(tài)度。我們應(yīng)該一步一步的來,多看看別人給你的步驟和例子,里面都是有的。 自己想當(dāng)然的算,然后靠別人幫你指出錯(cuò)誤,尤其是很簡(jiǎn)單又討論過無數(shù)次的錯(cuò)誤,這對(duì)自己而言雖然輕松了,但少了一個(gè)自己探索的過程。而且簡(jiǎn)單重復(fù)的回答,也浪費(fèi)了別人的時(shí)間。歸根結(jié)底,是自身沒有足夠的學(xué)習(xí)和調(diào)研造成許多根本上的卻很無謂的錯(cuò)誤。 說多了。 |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |
|
你這么說也不管用啊,畢竟測(cè)試也是個(gè)比較耗時(shí)間精力的過程,并且還可能證明了之前看起來跟實(shí)驗(yàn)非!胺稀钡慕Y(jié)果有問題了~ 事實(shí)上去年我在實(shí)驗(yàn)室里做組會(huì)報(bào)告時(shí),做的就是如何測(cè)試收斂的問題,被組里幾個(gè)也做計(jì)算的人說是很無聊和沒有特殊的東西的報(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é)果…… 總的來說,工科的急功近利的思想啊~ |

銅蟲 (小有名氣)

鐵桿木蟲 (正式寫手)
鐵桿木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +5 | 陌の森林 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +3 | 何潤(rùn)采123 2026-03-18 | 3/150 |
|
|
[考研] 266求調(diào)劑 +5 | 陽(yáng)陽(yáng)哇塞 2026-03-14 | 9/450 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-12 | 5/250 |
|
|
[論文投稿]
投稿問題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|