 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1188 | 回復(fù): 5 | |||
[交流]
【求助】關(guān)于復(fù)雜的結(jié)構(gòu)參數(shù)的優(yōu)化 已有3人參與
|
|
在侯柱峰老師的文章中看到有這么一句話（《VASP入門指南》-page25）：得到E-V曲線后，通過狀態(tài)方程擬合得到平衡狀態(tài)下的體積，計算出腳本中變量i的值。這個怎么計算啊？？？ 并改變i的循環(huán)值，這個是什么意思，具體是怎么操作的？？？ 再次運行run_cell計算一次，得到其他結(jié)構(gòu)參數(shù)c和位置u。。。這個從哪里可以看出？或者怎樣能夠得到？？ （腳本為書中例子）： #!/bin/sh rm WAVECAR for i in 2.81 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 do cat > POSCAR < Mg-hex $i 0.0 -1.0 0.0 0.8660254037844 0.5 0.0 0.0 0.0 1.6230529595 2 Direct 0.6666666666666667 0.3333333333333333 0.750 0.3333333333333333 0.6666666666666667 0.250 ! cp INCAR.relax INCAR echo "a = $i angstrom " ; time vasp cp CONTCAR POSCAR cp INCAR.static INCAR echo "a = $i angstrom " ; time vasp E=‘grep "TOTEN" OUTCAR | tail -1 | awk ’{printf "%12.6f \n", $5 }’‘ V=‘grep "volume" OUTCAR | tail -1 | awk ’{printf "%12.4f \n", $5 }’‘ echo $V $E >>EtVo.dat done [ Last edited by qiqi2926 on 2010-6-26 at 14:02 ] |
回復(fù)此樓» 猜你喜歡
085600 295分求調(diào)劑
已經(jīng)有14人回復(fù)
-
環(huán)境工程297分求調(diào)劑一志愿杭高院
已經(jīng)有3人回復(fù)
-
理學(xué)07化學(xué) 303求調(diào)劑
已經(jīng)有7人回復(fù)
-
290分調(diào)劑求助
已經(jīng)有6人回復(fù)
-
362求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑推薦 材料 304
已經(jīng)有17人回復(fù)
-
304求調(diào)劑
已經(jīng)有4人回復(fù)
-
張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298
已經(jīng)有8人回復(fù)
-
070300化學(xué)354求調(diào)劑
已經(jīng)有12人回復(fù)
-
08工科，295，接受跨專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助，高手幫幫忙，關(guān)于勢參數(shù)的問題
已經(jīng)有5人回復(fù)
- 結(jié)構(gòu)優(yōu)化參數(shù)設(shè)置
已經(jīng)有6人回復(fù)
- 關(guān)于鋁和金屬氧化物表面的理論計算
已經(jīng)有7人回復(fù)
- 最復(fù)雜情況下的參數(shù)優(yōu)化
已經(jīng)有11人回復(fù)
- 關(guān)于發(fā)布國家重點基礎(chǔ)研究發(fā)展計劃2011年項目申報指南
已經(jīng)有7人回復(fù)
- 【求助】VASP怎么做晶格參數(shù)的結(jié)構(gòu)優(yōu)化
已經(jīng)有17人回復(fù)
- 【原創(chuàng)】關(guān)于結(jié)構(gòu)弛豫設(shè)置及收斂判據(jù)選擇的一點體會
已經(jīng)有140人回復(fù)
hedaors
禁蟲 (正式寫手)
隱清堂主
- 應(yīng)助: 11 (小學(xué)生)
- 貴賓: 0.225
- 金幣: 3948.8
- 散金: 201
- 紅花: 8
- 帖子: 855
- 在線: 271.1小時
- 蟲號: 449622
- 注冊: 2007-11-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2010-06-26 13:04:56
3樓2010-06-26 14:03:48
hedaors
禁蟲 (正式寫手)
隱清堂主
- 應(yīng)助: 11 (小學(xué)生)
- 貴賓: 0.225
- 金幣: 3948.8
- 散金: 201
- 紅花: 8
- 帖子: 855
- 在線: 271.1小時
- 蟲號: 449622
- 注冊: 2007-11-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
gavinliu7390(金幣+2):謝謝交流！ 2010-06-26 15:22:42
cenwanglai(金幣+1):恩,費心了! 2010-06-26 19:33:51
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
gavinliu7390(金幣+2):謝謝交流！ 2010-06-26 15:22:42
cenwanglai(金幣+1):恩,費心了! 2010-06-26 19:33:51
|
for i in 2.81 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 這個語句是一個循環(huán)，每改變一個 i 值，做一個計算。 你沒用過 Linux ? 你需要下功夫補補 Linux 的 Shell 腳本知識了，用 VASP 必須要會的。 這有個不錯的教程: http://tldp.org/LDP/abs/html/ 也可一看看中文版 http://wenku.baidu.com/view/ab2cfe21dd36a32d73758196.html |
4樓2010-06-26 14:20:28
★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
cenwanglai(金幣+1):呵呵,直接用二次曲線擬合得到i就可以了. 2010-06-26 19:34:27
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
cenwanglai(金幣+1):呵呵,直接用二次曲線擬合得到i就可以了. 2010-06-26 19:34:27
|
在小木蟲上搜murn.f程序，用來擬合晶格常數(shù)，具體的用法就在侯老師那本書的后面。 擬合出來晶格常數(shù)的值（即i的值）后，在做一次結(jié)構(gòu)弛豫即INCAR.relax便可得到其他結(jié)構(gòu)參數(shù)的值了，具體怎么得到要依據(jù)晶胞的結(jié)構(gòu)特征。不知我這樣理解是否正確，還請高手指點，僅供樓主參考，共同學(xué)習(xí)！ |
5樓2010-06-26 17:15:10
6樓2010-06-28 09:33:59
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 qiqi2926 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 295分求調(diào)劑 +12 | W55j 2026-03-30 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 290分調(diào)劑求助 +6 | 吉祥止止陳 2026-03-25 | 6/300 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +17 | 荷包蛋hyj 2026-03-26 | 17/850 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +8 | 小黃鴨寶 2026-03-30 | 8/400 |
|
|
|
[考研] 機械學(xué)碩總分317求調(diào)劑�。。�！ +5 | Acaciad 2026-03-25 | 5/250 |
|
|
[考研] 274求調(diào)劑 +6 | xiao愛同學(xué) 2026-03-30 | 6/300 |
|
|
|
[考研] 看看我！301求調(diào)劑 081700 +3 | 一個小明同學(xué) 2026-03-24 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3 | 小圓圈圈ooo 2026-03-30 | 3/150 |
|
|
|
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +10 | Koxui 2026-03-28 | 10/500 |
|
|
|
[考研] 286求調(diào)劑 +5 | Faune 2026-03-30 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +10 | lin0039 2026-03-26 | 10/500 |
|
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +7 | 小蠟新筆 2026-03-29 | 7/350 |
|
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +7 | 爭取九點睡 2026-03-28 | 8/400 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|