 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 28 | 1/1 | 返回列表 |
| 查看: 6317 | 回復(fù): 27 | |||||
| 本帖產(chǎn)生 1 個(gè) 模擬EPI ，點(diǎn)擊這里進(jìn)行查看 | |||||
[交流]
【求助】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 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 天地靈氣，日月精華 | gulp |
» 猜你喜歡
276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷
已經(jīng)有10人回復(fù)
-
一志愿南航材料專317分求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)學(xué)碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有12人回復(fù)
-
求材料，環(huán)境專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
335求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 在做CASTP模擬時(shí)，建立晶格需要晶格常數(shù)，如何獲取呢？
已經(jīng)有10人回復(fù)
- 用LDA贗勢(shì)優(yōu)化出來(lái)的晶格常數(shù)一定比實(shí)驗(yàn)值小嗎？謝謝！
已經(jīng)有10人回復(fù)
- 晶格常數(shù)優(yōu)化
已經(jīng)有7人回復(fù)
- 【求助】castep 怎么預(yù)測(cè)晶格常數(shù)
已經(jīng)有8人回復(fù)
- 【求助】擬合晶格常數(shù)
已經(jīng)有14人回復(fù)
- 【求助】原子替換后晶格常數(shù)的改變？？？？
已經(jīng)有22人回復(fù)
- 【求助】si晶格常數(shù)
已經(jīng)有7人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
94年女征男友
+1/271
-
2026年國(guó)基求好運(yùn)！
+5/190
-
上海真誠(chéng)脫單，這是一篇很認(rèn)真的相親貼，找一個(gè)靈魂伴侶
+1/167
-
福建師范大學(xué)葉應(yīng)祥課題組招收2026級(jí)博士研究生
+2/132
-
濟(jì)南大學(xué)林秀娟教授招收材料類博士生
+1/82
-
國(guó)自然面上祈福
+1/79
-
有沒(méi)有在香港的蟲友啊
+1/72
-
散金幣，求基金中標(biāo)！！
+1/60
-
武漢高校接受研究生調(diào)劑，材料、化學(xué)、紡織方向
+2/44
-
中國(guó)科學(xué)院杭州醫(yī)學(xué)研究所覃江江課題組招聘博士后
+1/30
-
湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院 課題組招收2026級(jí)食品/生物方向碩士
+1/28
-
浙江大學(xué)光電學(xué)院極端光學(xué)技術(shù)與儀器國(guó)家重點(diǎn)實(shí)驗(yàn)室誠(chéng)聘極紫外阿秒方向博士后
+1/26
-
南京林業(yè)大學(xué)-國(guó)家級(jí)青年人才團(tuán)隊(duì) 招2026級(jí)博士、碩士（合成化學(xué)、植物化學(xué)方向）
+1/15
-
新加坡國(guó)立大學(xué)藥學(xué)系化學(xué)生物學(xué)課題組招PhD
+1/5
-
深圳信息職業(yè)大學(xué)招聘聯(lián)合培養(yǎng)碩博研究生（急）
+1/5
-
211/雙一流石河子大學(xué)化學(xué)化工學(xué)院--電催化，CO2轉(zhuǎn)化
+1/3
-
歡迎生物與醫(yī)藥、藥學(xué)、化學(xué)等相關(guān)專業(yè)的同學(xué)
+1/3
-
5科大招收材料專業(yè)碩士
+1/3
-
香港中文大學(xué)化學(xué)系吳基培教授招募2026年秋季（2026年8月1日）入學(xué)全獎(jiǎng)博士研究生
+1/2
-
上海高校教授長(zhǎng)年招收碩士和博士研究生-動(dòng)力工程與工程熱物理
+1/1
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
ghcacj(金幣+10, 模擬EPI+1): 謝謝 2011-02-24 13:33:37
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
ghcacj(金幣+10, 模擬EPI+1): 謝謝 2011-02-24 13:33:37
你好。昨晚我用你的輸入文件作了測(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之類，如果tempratio大得很，比如2000，3000之類，你就沒(méi)有必要再搞了，肯定問(wèn)題很大。 再一個(gè)你的輸入文件太麻煩了。對(duì)FCC金屬，可以這樣寫： 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)然輸入文件的寫入是小事，關(guān)鍵還是溫度控不住的情況比較煩人。而且一般是突然發(fā)生的錯(cuò)誤。實(shí)在不行你可以給程序作者發(fā)Email，說(shuō)說(shuō)這個(gè)問(wèn)題，這個(gè)作者Julian Gale是個(gè)很好、很有耐心的人。要不就換DLPOLY做MD吧，DLPOLY計(jì)算MD專業(yè)一點(diǎn)。呵呵呵。 |
18樓2011-02-23 09:18:46
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
2樓2011-01-07 12:16:07
3樓2011-01-07 17:27:10
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
4樓2011-01-07 17:39:48
5樓2011-01-07 18:35:00
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
6樓2011-01-07 19:11:23
|
我試了用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 ] |
7樓2011-01-07 19:18:19
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
8樓2011-01-08 11:00:32
|
再請(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)，謝謝！ |
9樓2011-02-17 10:49:47
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
10樓2011-02-18 21:53:25
11樓2011-02-19 18:29:53
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
12樓2011-02-22 13:13:34
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
13樓2011-02-22 13:14:43
14樓2011-02-22 14:06:03
15樓2011-02-22 16:18:32
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
省那些坐標(biāo)干什么？保密么？16樓2011-02-22 17:06:09
17樓2011-02-22 19:20:33
19樓2011-02-23 09:48:29
20樓2014-02-24 19:37:47
oxox6085
專家顧問(wèn) (正式寫手)
-
專家經(jīng)驗(yàn): +129 - 模擬EPI: 1
- 應(yīng)助: 69 (初中生)
- 金幣: 8061
- 帖子: 473
- 在線: 636.8小時(shí)
- 蟲號(hào): 612537
★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
月只藍(lán): 金幣+1, 鼓勵(lì)交流！ 2014-02-25 11:26:27
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
月只藍(lán): 金幣+1, 鼓勵(lì)交流！ 2014-02-25 11:26:27
|
用準(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 |
21樓2014-02-25 06:02:40
22樓2014-02-25 16:46:47
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
23樓2014-03-03 14:14:08
24樓2014-03-03 16:43:48
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
25樓2014-03-05 19:47:58
26樓2014-03-05 22:51:33
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 0.17
- 金幣: 4774.1
- 帖子: 1361
- 在線: 215.2小時(shí)
- 蟲號(hào): 659094
27樓2014-03-11 09:28:03
28樓2018-12-12 17:56:04
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 lchunmei001 的主題更新
| 28 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +10 | Zhao anqi 2026-03-22 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +7 | 紙魚ly 2026-03-21 | 8/400 |
|
|
[考研] 材料專碩找調(diào)劑 +3 | 哈哈哈吼吼吼哈 2026-03-23 | 3/150 |
|
|
|
[考研] 350求調(diào)劑 +6 | weudhdk 2026-03-19 | 6/300 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +5 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 5/250 |
|
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
|
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3 | 1911236844 2026-03-17 | 3/150 |
|
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫郵件了，沒(méi)回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|