 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1337 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
zhang668木蟲 (正式寫手)
|
[交流]
【求助】做狀態(tài)方程計(jì)算時(shí)候得到的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ù)的計(jì)算。請問腳本那里出了問題？ 還有如果不想得到體模量，只想優(yōu)化，什么時(shí)候才有做狀態(tài)方程擬合晶格參數(shù)必要，一般直接做全優(yōu)化可以嗎？還有做狀態(tài)方程的時(shí)候，沒有優(yōu)化自由度，OUTCAR最后FORCE比較大，這個能量也不能說明問題吧 |
回復(fù)此樓» 猜你喜歡
343求調(diào)劑085601
已經(jīng)有3人回復(fù)
-
08工科求調(diào)劑286
已經(jīng)有3人回復(fù)
-
337求調(diào)劑
已經(jīng)有3人回復(fù)
-
329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601）
已經(jīng)有10人回復(fù)
-
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有16人回復(fù)
-
291求調(diào)劑
已經(jīng)有6人回復(fù)
-
299求調(diào)劑
已經(jīng)有10人回復(fù)
-
面上5B能上會嗎？
已經(jīng)有4人回復(fù)
-
322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級。
已經(jīng)有4人回復(fù)
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有14人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 如何求解描述振動的二階微分方程
已經(jīng)有21人回復(fù)
- 橢圓形偏微分方程標(biāo)準(zhǔn)化后該如何解？
已經(jīng)有6人回復(fù)
- 有關(guān)磁流體方程計(jì)算的詢問
已經(jīng)有8人回復(fù)
- 怎么才能做出三個9的線性方程�。�？��！愁死我了
已經(jīng)有11人回復(fù)
- 求助Ta2O5的色散方程
已經(jīng)有4人回復(fù)
- PRSV方程中的k
已經(jīng)有4人回復(fù)
- 大家看看我的excell擬合方程，怎么能合理的解釋相關(guān)性呢？
已經(jīng)有8人回復(fù)
- matlab擬合方程參數(shù)時(shí)初值的選擇
已經(jīng)有15人回復(fù)
- R語言圖書-求解差分方程和并行計(jì)算！
已經(jīng)有43人回復(fù)
- 求助這個非線性方程改怎么解��？？
已經(jīng)有4人回復(fù)
- 求教origin擬合簡單的滲流方程（又稱逾滲方程），做復(fù)合材料的應(yīng)該會遇到吧
已經(jīng)有20人回復(fù)
- 【求助】如何計(jì)算實(shí)驗(yàn)數(shù)據(jù)與已知方程的相關(guān)系數(shù)？
已經(jīng)有6人回復(fù)
wangzg9292
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3933
- 散金: 56
- 帖子: 348
- 在線: 89.6小時(shí)
- 蟲號: 723844
- 注冊: 2009-03-16
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2010-04-24 09:05:33
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時(shí)
- 蟲號: 675997
- 注冊: 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
zhang668(金幣+3): 2010-04-24 10:16
|
別人的腳本，不要拿過來就用，最好試一試。 你可以找一個算完的OUTCAR試一試，確定腳本沒問題了，再進(jìn)行計(jì)算。 EtVo.dat是空的，就說明是E和V的變量沒有賦進(jìn)去值。 我看你用的符號都是' '，應(yīng)該是一個`,一個是' E='grep ''TOTEN'' OUTCAR | tail -1 | awk '{printf ''%12.6f \n'' , $5}' ' 應(yīng)該是E=`grep ''TOTEN'' OUTCAR | tail -1 | awk '{printf ''%12.6f \n'' , $5}' ` 擬合物態(tài)方程絕對不是為了獲得B0， 大部分的目的還是認(rèn)為壓力算的不準(zhǔn)，不如擬合物態(tài)方程算出的壓力準(zhǔn)。至于擬合物態(tài)方程也緊緊限于立方結(jié)構(gòu)，非立方結(jié)構(gòu)由于加壓之后，很可能晶格常數(shù)的比例變了，必須全優(yōu)化后，才能擬合物態(tài)方程，也有不擬合的。 至于你說的FORCE很大， 那最大的原因就是你沒有relax離子，否則不可能很大，或者你的收斂參數(shù)不夠高。對應(yīng)的就是isif=2， nsw=60 立方體系也可以全優(yōu)化，但是那樣比較浪費(fèi)時(shí)間，沒有這樣的快！ |
2樓2010-04-23 23:35:28
zhang668
木蟲 (正式寫手)
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 11 (小學(xué)生)
- 貴賓: 0.061
- 金幣: 1735
- 散金: 1294
- 紅花: 11
- 帖子: 843
- 在線: 324.9小時(shí)
- 蟲號: 733242
- 注冊: 2009-03-28
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2010-04-24 11:08:41
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時(shí)
- 蟲號: 675997
- 注冊: 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
5樓2010-04-24 11:23:40
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +6 | 26要上岸 2026-03-28 | 6/300 |
|
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3 | llss0711 2026-03-28 | 6/300 |
|
|
|
[考研] 0856材料化工調(diào)劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|
|
[考研] 調(diào)劑求院校招收 +6 | 鶴鯨鴿 2026-03-28 | 6/300 |
|
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
|
[考研] 085600，材料與化工321分求調(diào)劑 +9 | 大饞小子 2026-03-28 | 9/450 |
|
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研]
材料學(xué)碩，求調(diào)劑
 6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|