 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 1318 | 回復: 9 | |||
[交流]
【求助】做狀態(tài)方程計算時候得到的EtVo.dat是空的 已有4人參與
|
|
下面是我按照侯博士教程寫的腳本： #!/bin/sh for i in 7.846 7.946 8.046 8.146 8.246 8.346 8.446 8.546 8.646 8.746 8.846 8.946 9.046 do cat > POSCAR < aaa $i 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 2 4 8 Direct 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.2500000000000000 0.2500000000000000 0.2500000000000000 0.6205000281333946 0.6205000281333943 0.6205000281333943 0.6205000281333944 0.6205000281333944 -1.8615000844001834 0.6205000281333947 -1.8615000844001834 0.6205000281333946 -1.8615000844001834 0.6205000281333947 0.6205000281333944 0.3849999904632571 0.3849999904632571 0.3849999904632571 0.8499000072479258 0.8499000072479258 0.8499000072479256 0.3849999904632571 0.3849999904632571 -1.1549999713897712 0.8499000072479258 0.8499000072479258 -2.5497000217437775 0.3849999904632572 -1.1549999713897714 0.3849999904632571 0.8499000072479259 -2.5497000217437775 0.8499000072479259 -1.1549999713897714 0.3849999904632572 0.3849999904632571 -2.5497000217437775 0.8499000072479260 0.8499000072479258 ! 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 結(jié)果輸出文件EtVo.dat是空的，看了OUTCAR，最終結(jié)果只有最后一個晶格參數(shù)的計算。請問腳本那里出了問題？ 還有如果不想得到體模量，只想優(yōu)化，什么時候才有做狀態(tài)方程擬合晶格參數(shù)必要，一般直接做全優(yōu)化可以嗎？還有做狀態(tài)方程的時候，沒有優(yōu)化自由度，OUTCAR最后FORCE比較大，這個能量也不能說明問題吧 |
回復此樓» 猜你喜歡
280求調(diào)劑
已經(jīng)有7人回復
-
環(huán)境領域全國重點實驗室招收博士1-2名
已經(jīng)有5人回復
-
303求調(diào)劑
已經(jīng)有6人回復
-
290求調(diào)劑
已經(jīng)有10人回復
-
0703化學求調(diào)劑 總分331
已經(jīng)有3人回復
-
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有8人回復
-
0703化學調(diào)劑
已經(jīng)有3人回復
-
材料，紡織，生物（0856、0710），化學招生啦
已經(jīng)有9人回復
-
一志愿中國海洋大學，生物學，301分，求調(diào)劑
已經(jīng)有3人回復
-
0703化學調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有9人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 如何求解描述振動的二階微分方程
已經(jīng)有21人回復
- 橢圓形偏微分方程標準化后該如何解？
已經(jīng)有6人回復
- 有關磁流體方程計算的詢問
已經(jīng)有8人回復
- 怎么才能做出三個9的線性方程�。�？��！愁死我了
已經(jīng)有11人回復
- 求助Ta2O5的色散方程
已經(jīng)有4人回復
- PRSV方程中的k
已經(jīng)有4人回復
- 大家看看我的excell擬合方程，怎么能合理的解釋相關性呢？
已經(jīng)有8人回復
- matlab擬合方程參數(shù)時初值的選擇
已經(jīng)有15人回復
- R語言圖書-求解差分方程和并行計算！
已經(jīng)有43人回復
- 求助這個非線性方程改怎么解��？？
已經(jīng)有4人回復
- 求教origin擬合簡單的滲流方程（又稱逾滲方程），做復合材料的應該會遇到吧
已經(jīng)有20人回復
- 【求助】如何計算實驗數(shù)據(jù)與已知方程的相關系數(shù)？
已經(jīng)有6人回復
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強帖: 3
- 應助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時
- 蟲號: 675997
- 注冊: 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和
zhang668(金幣+3): 2010-04-24 10:16
|
別人的腳本，不要拿過來就用，最好試一試。 你可以找一個算完的OUTCAR試一試，確定腳本沒問題了，再進行計算。 EtVo.dat是空的，就說明是E和V的變量沒有賦進去值。 我看你用的符號都是' '，應該是一個`,一個是' E='grep ''TOTEN'' OUTCAR | tail -1 | awk '{printf ''%12.6f \n'' , $5}' ' 應該是E=`grep ''TOTEN'' OUTCAR | tail -1 | awk '{printf ''%12.6f \n'' , $5}' ` 擬合物態(tài)方程絕對不是為了獲得B0， 大部分的目的還是認為壓力算的不準，不如擬合物態(tài)方程算出的壓力準。至于擬合物態(tài)方程也緊緊限于立方結(jié)構(gòu)，非立方結(jié)構(gòu)由于加壓之后，很可能晶格常數(shù)的比例變了，必須全優(yōu)化后，才能擬合物態(tài)方程，也有不擬合的。 至于你說的FORCE很大， 那最大的原因就是你沒有relax離子，否則不可能很大，或者你的收斂參數(shù)不夠高。對應的就是isif=2， nsw=60 立方體系也可以全優(yōu)化，但是那樣比較浪費時間，沒有這樣的快！ |
2樓2010-04-23 23:35:28
wangzg9292
木蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 3933
- 散金: 56
- 帖子: 348
- 在線: 89.6小時
- 蟲號: 723844
- 注冊: 2009-03-16
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2010-04-24 09:05:33
4樓2010-04-24 11:08:41
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強帖: 3
- 應助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時
- 蟲號: 675997
- 注冊: 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和
5樓2010-04-24 11:23:40
經(jīng)世濟民
禁蟲 (正式寫手)
zhang668(金幣+1): 2010-04-24 11:27
|
本帖內(nèi)容被屏蔽 |
6樓2010-04-24 11:24:40
7樓2010-04-24 11:27:42
8樓2010-04-24 11:29:34
9樓2010-04-24 11:45:33
經(jīng)世濟民
禁蟲 (正式寫手)
★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
|
本帖內(nèi)容被屏蔽 |
10樓2010-04-24 13:11:28
相關版塊跳轉(zhuǎn)
我要訂閱樓主 zhang668 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 290求調(diào)劑 +10 | ADT 2026-03-11 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|
|
[考研] 【0856】化學工程（085602）313 分，本科學科評估A類院�；瘜W工程與工藝，誠求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
|
[考研] 085600材料與化工 +4 | 安全上岸！ 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿南京大學，080500材料科學與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 材料工程327求調(diào)劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 【0703化學調(diào)劑】-一志愿華中師范大學-六級475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
|
[考研] 331求調(diào)劑（0703有機化學 +5 | ZY-05 2026-03-13 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[論文投稿]
投稿問題
 5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|