| 6 | 1/1 | 返回列表 |
| 查看: 1189 | 回復(fù): 5 | |||
qiqi2926銅蟲 (小有名氣)
|
[交流]
【求助】關(guān)于復(fù)雜的結(jié)構(gòu)參數(shù)的優(yōu)化 已有3人參與
|
|
在侯柱峰老師的文章中看到有這么一句話(《VASP入門指南》-page25):得到E-V曲線后,通過狀態(tài)方程擬合得到平衡狀態(tài)下的體積,計(jì)算出腳本中變量i的值。這個(gè)怎么計(jì)算??? 并改變i的循環(huán)值,這個(gè)是什么意思,具體是怎么操作的??? 再次運(yùn)行run_cell計(jì)算一次,得到其他結(jié)構(gòu)參數(shù)c和位置u。。。這個(gè)從哪里可以看出?或者怎樣能夠得到?? (腳本為書中例子): #!/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 ] |
禁蟲 (正式寫手)
隱清堂主
銅蟲 (小有名氣)
禁蟲 (正式寫手)
隱清堂主
|
for i in 2.81 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 這個(gè)語句是一個(gè)循環(huán),每改變一個(gè) i 值,做一個(gè)計(jì)算。 你沒用過 Linux ? 你需要下功夫補(bǔ)補(bǔ) Linux 的 Shell 腳本知識(shí)了,用 VASP 必須要會(huì)的。 這有個(gè)不錯(cuò)的教程: http://tldp.org/LDP/abs/html/ 也可一看看中文版 http://wenku.baidu.com/view/ab2cfe21dd36a32d73758196.html |
|
在小木蟲上搜murn.f程序,用來擬合晶格常數(shù),具體的用法就在侯老師那本書的后面。 擬合出來晶格常數(shù)的值(即i的值)后,在做一次結(jié)構(gòu)弛豫即INCAR.relax便可得到其他結(jié)構(gòu)參數(shù)的值了,具體怎么得到要依據(jù)晶胞的結(jié)構(gòu)特征。不知我這樣理解是否正確,還請高手指點(diǎn),僅供樓主參考,共同學(xué)習(xí)! |
銅蟲 (小有名氣)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +10 | 墨墨漠 2026-03-25 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 生醫(yī)工0831調(diào)劑求推薦 +4 | 小熊睿睿_s 2026-03-27 | 6/300 |
|
|
[考研] 生物學(xué)296求調(diào)劑 +8 | 湯圓包 2026-03-29 | 12/600 |
|
|
[考研] 322求調(diào)劑:一志愿湖南大學(xué) 材料與化工(085600),已過六級。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +28 | RichLi_ 2026-03-25 | 36/1800 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +9 | 柳溪邊 2026-03-26 | 9/450 |
|
|
[考研] 一志愿中國科學(xué)院大學(xué)265求調(diào)劑 +6 | 恬淡ye 2026-03-31 | 7/350 |
|
|
[考研] 一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-30 | 5/250 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +5 | 小張zxy張 2026-03-26 | 10/500 |
|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工(085600)296求調(diào)劑 +25 | 稻妻小編 2026-03-26 | 25/1250 |
|
|
[考研] 0703化學(xué)321分求調(diào)劑 +10 | 三dd. 2026-03-30 | 11/550 |
|
|
[考研] 求調(diào)劑 +10 | 家佳佳佳佳佳 2026-03-29 | 10/500 |
|
|
[考研] 一志愿北京工業(yè)大學(xué),324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|