| 24小時(shí)熱門(mén)版塊排行榜 |
| 28 | 1/1 | 返回列表 |
| 查看: 6316 | 回復(fù): 27 | |||||
| 本帖產(chǎn)生 1 個(gè) 模擬EPI ,點(diǎn)擊這里進(jìn)行查看 | |||||
lchunmei001木蟲(chóng) (正式寫(xiě)手)
|
[交流]
【求助】gulp中改變溫度,晶格常數(shù)沒(méi)變?
|
||||
|
1.請(qǐng)教一下大家:我在gulp輸入文件中改變溫度,為什么相應(yīng)的晶格常數(shù)沒(méi)有變化? 輸入文件如下: conp opti prop comp cell 4.000000 4.000000 4.000000 90.000000 90.000000 90.000000 fractional 4 Ni core 0.0000000 0.0000000 0.0000000 0.00000000 1.0000 0.00000 Ni core 0.0000000 0.5000000 0.5000000 0.00000000 1.0000 0.00000 Ni core 0.5000000 0.0000000 0.5000000 0.00000000 1.0000 0.00000 Ni core 0.5000000 0.5000000 0.0000000 0.00000000 1.0000 0.00000 species Ni core 0.000 pressure 0.0 temperature 300.0 manybody Ni core Ni core 0.0 12.0 eam_functional square_root Ni core 1.000 eam_density power 6 Ni core 729.6928296 lennard 9 6 Ni core Ni core 1303.09841 0.0 0.0 12.0 dump example17.grs 發(fā)現(xiàn)temperature隨便改為多少,得到的晶格常數(shù)都沒(méi)變?望高手指點(diǎn)指點(diǎn)。我是新手,所以很多不明白!謝謝大家! 2.跑MD是出現(xiàn)以下錯(cuò)誤,不知調(diào)哪些參數(shù)可以解決,試調(diào)整了MD相關(guān)參數(shù),還是有此錯(cuò)誤! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! ERROR : temperature has exceeded maximum allowed set by mdmaxtemp !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
天地靈氣,日月精華 | gulp |
» 搶金幣啦!回帖就可以得到:
+3/1736
+2/52
+1/49
+1/36
+1/34
+1/28
+1/27
+1/25
+1/19
+1/14
+1/7
+1/7
+1/6
+2/4
+1/4
+1/3
+1/3
+1/3
+1/2
+1/1
木蟲(chóng) (著名寫(xiě)手)
你好。昨晚我用你的輸入文件作了測(cè)試,好像確實(shí)有問(wèn)題。我另?yè)Q了Johnson勢(shì)函數(shù),也改過(guò)幾次控溫參數(shù),均發(fā)現(xiàn)有問(wèn)題,有時(shí)候是溫度,有時(shí)候是壓力。但這樣的問(wèn)題對(duì)于晶體的模擬是不該出現(xiàn)的,我一時(shí)也看不出是程序有Bug還是輸入有問(wèn)題。![]() 你用的Sutton-Chen勢(shì)太長(zhǎng)了,運(yùn)行很慢,我沒(méi)有時(shí)間繼續(xù)測(cè)試下去,我建議在找到問(wèn)題癥結(jié)之前,你先把Gulp程序中的mxmdtemp參數(shù)增大,在initial.f90中找到 rmdmaxtemp = 100.0_dp, 把100改為1000,這樣容錯(cuò)性就好一些。注意改后要重新編譯生成新的可執(zhí)行文件。 但如果這樣還是有問(wèn)題,那就真的是不好辦了。另外,請(qǐng)你打開(kāi)mdprop.f90這個(gè)文件,找到大約127行處,有如下代碼: if (tempratio.gt.rmdmaxtemp) then call outerror('temperature has exceeded maximum allowed set by mdmaxtemp',0_i4) call gulpfinish endif 你還可以在call gulpfinish 這一句之前(即call outerror這一句之后),添加如下代碼: write(ioout,'('' ** Time : '',f18.5,'' ps :'')') nsteps*tstep(ncf) write(ioout,*) 'temperature, targettemperature, tempratio:' write(ioout,*) temperature, targettemperature, tempratio 改后要重新編譯。這個(gè)是在出錯(cuò)的時(shí)候輸出錯(cuò)在什么地方。這一句對(duì)改錯(cuò)是沒(méi)什么用的,只是弄明白錯(cuò)的程度有多大,如果輸出的tempratio不是特別大(比如你前面設(shè)的rmdmaxtemp是1000,這里 輸出的tempratio是1020)你就可以再加大rmdmaxtemp,比如設(shè)成1200之類(lèi),如果tempratio大得很,比如2000,3000之類(lèi),你就沒(méi)有必要再搞了,肯定問(wèn)題很大。 再一個(gè)你的輸入文件太麻煩了。對(duì)FCC金屬,可以這樣寫(xiě): nosy cart md conp cell 4.09 4.09 4.09 90.000000 90.000000 90.000000 frac Ag 0.0 0.0 0.0 Ag 0.5 0.0 0.5 Ag 0.5 0.5 0.0 Ag 0.0 0.5 0.5 supercell 7 7 7 libr suttonchen 下面就是pressure=0這些,是你的MD控制參數(shù)。 這樣你換其它FCC體系也容易,只要改晶格參數(shù)、元素名就可以了。suttonchen勢(shì)是在Gulp的勢(shì)函數(shù)庫(kù)當(dāng)中的。用的時(shí)候把Libraries中的suttonchen.lib拷到當(dāng)前運(yùn)行目錄下。用7*7*7的晶格也太大了,大概256或500個(gè),對(duì)于晶格優(yōu)化也就差不多了,864也可以。時(shí)間步數(shù)大概0.001ps就夠了。你對(duì)金屬Ag,應(yīng)該在熔點(diǎn)以下運(yùn)行,溫度不要太高。 當(dāng)然輸入文件的寫(xiě)入是小事,關(guān)鍵還是溫度控不住的情況比較煩人。而且一般是突然發(fā)生的錯(cuò)誤。實(shí)在不行你可以給程序作者發(fā)Email,說(shuō)說(shuō)這個(gè)問(wèn)題,這個(gè)作者Julian Gale是個(gè)很好、很有耐心的人。要不就換DLPOLY做MD吧,DLPOLY計(jì)算MD專(zhuān)業(yè)一點(diǎn)。呵呵呵。 |
木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
|
我試了用MD,改變溫度有影響,但出現(xiàn)上面第二個(gè)錯(cuò)誤,即: !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! ERROR : temperature has exceeded maximum allowed set by mdmaxtemp !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 要不就是: !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! ERROR : volume has exceeded maximum allowed set by mdmaxvolume !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 動(dòng)力學(xué)跑不完,還請(qǐng)指教! [ Last edited by lchunmei001 on 2011-1-7 at 20:32 ] |
木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
|
再請(qǐng)教老虎大王:我也按你上面說(shuō)的,增加體系的原子數(shù)(1300多個(gè))或減小時(shí)間步長(zhǎng),在低溫下確實(shí)可以解決上面出現(xiàn)的兩個(gè)錯(cuò)誤,但溫度稍微高些(如1500K以上),又出現(xiàn)以下錯(cuò)誤。 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! ERROR : temperature has exceeded maximum allowed set by mdmaxtemp !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 動(dòng)力學(xué)跑不完。還請(qǐng)?jiān)僦更c(diǎn)指點(diǎn),謝謝! |
木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (知名作家)
木蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +129 |
|
用準(zhǔn)簡(jiǎn)諧近似(QHA)的方法考慮溫度效應(yīng),注意溫度不要達(dá)到1/2熔點(diǎn),不然非諧效應(yīng)主導(dǎo),用準(zhǔn)簡(jiǎn)諧近似無(wú)法算準(zhǔn): conp free prop comp cell 4.000000 4.000000 4.000000 90.000000 90.000000 90.000000 fractional 4 Ni core 0.0000000 0.0000000 0.0000000 0.00000000 1.0000 0.00000 Ni core 0.0000000 0.5000000 0.5000000 0.00000000 1.0000 0.00000 Ni core 0.5000000 0.0000000 0.5000000 0.00000000 1.0000 0.00000 Ni core 0.5000000 0.5000000 0.0000000 0.00000000 1.0000 0.00000 species Ni core 0.000 pressure 0.0 temperature 300.0 shrink 10 10 10 manybody Ni core Ni core 0.0 12.0 eam_functional square_root Ni core 1.000 eam_density power 6 Ni core 729.6928296 lennard 9 6 Ni core Ni core 1303.09841 0.0 0.0 12.0 dump example17.grs |
銀蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)
銀蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)
銀蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)
|
你好,請(qǐng)問(wèn)問(wèn)題解決了嗎,我最近也遇到了這個(gè)問(wèn)題 發(fā)自小木蟲(chóng)Android客戶(hù)端 |
| 28 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +7 | 紙魚(yú)ly 2026-03-21 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 材料專(zhuān)碩英一數(shù)二306 +8 | z1z2z3879 2026-03-18 | 8/400 |
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專(zhuān)碩324分,求調(diào)劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-17 | 10/500 |
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專(zhuān)業(yè) 281分,一志愿四川大學(xué),誠(chéng)心求調(diào)劑 +11 | 吃吃吃才有意義 2026-03-19 | 11/550 |
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 南昌大學(xué)材料專(zhuān)碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽(yáng)CF 2026-03-18 | 5/250 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 一志愿西南交通 專(zhuān)碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專(zhuān)材355 2026-03-19 | 5/250 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲(chóng)200408204 2026-03-18 | 3/150 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|