 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1339 | 回復(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é)果只有最后一個(gè)晶格參數(shù)的計(jì)算。請問腳本那里出了問題？ 還有如果不想得到體模量，只想優(yōu)化，什么時(shí)候才有做狀態(tài)方程擬合晶格參數(shù)必要，一般直接做全優(yōu)化可以嗎？還有做狀態(tài)方程的時(shí)候，沒有優(yōu)化自由度，OUTCAR最后FORCE比較大，這個(gè)能量也不能說明問題吧 |
回復(fù)此樓» 猜你喜歡
0703化學(xué)321分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有19人回復(fù)
-
298求調(diào)劑
已經(jīng)有5人回復(fù)
-
346求調(diào)劑 一志愿070303有機(jī)化學(xué)
已經(jīng)有6人回復(fù)
-
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行
已經(jīng)有3人回復(fù)
-
324求調(diào)劑
已經(jīng)有11人回復(fù)
-
307求調(diào)劑
已經(jīng)有13人回復(fù)
-
071010 323 分求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料調(diào)劑
已經(jīng)有9人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 如何求解描述振動(dòng)的二階微分方程
已經(jīng)有21人回復(fù)
- 橢圓形偏微分方程標(biāo)準(zhǔn)化后該如何解？
已經(jīng)有6人回復(fù)
- 有關(guān)磁流體方程計(jì)算的詢問
已經(jīng)有8人回復(fù)
- 怎么才能做出三個(gè)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ù)
- 求助這個(gè)非線性方程改怎么解�。�？
已經(jīng)有4人回復(fù)
- 求教origin擬合簡單的滲流方程（又稱逾滲方程），做復(fù)合材料的應(yīng)該會(huì)遇到吧
已經(jīng)有20人回復(fù)
- 【求助】如何計(jì)算實(shí)驗(yàn)數(shù)據(jù)與已知方程的相關(guān)系數(shù)？
已經(jīng)有6人回復(fù)
zhang668
木蟲 (正式寫手)
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 11 (小學(xué)生)
- 貴賓: 0.061
- 金幣: 1735
- 散金: 1294
- 紅花: 11
- 帖子: 843
- 在線: 324.9小時(shí)
- 蟲號(hào): 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í)
- 蟲號(hào): 675997
- 注冊: 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
zhang668(金幣+3): 2010-04-24 10:16
|
別人的腳本，不要拿過來就用，最好試一試。 你可以找一個(gè)算完的OUTCAR試一試，確定腳本沒問題了，再進(jìn)行計(jì)算。 EtVo.dat是空的，就說明是E和V的變量沒有賦進(jìn)去值。 我看你用的符號(hào)都是' '，應(yīng)該是一個(gè)`,一個(gè)是' 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
wangzg9292
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3933
- 散金: 56
- 帖子: 348
- 在線: 89.6小時(shí)
- 蟲號(hào): 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í)
- 蟲號(hào): 675997
- 注冊: 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
5樓2010-04-24 11:23:40
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +8 | hanamiko 2026-03-27 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +11 | 芋泥小鈴鐺 2026-03-28 | 11/550 |
|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +10 | 景隨風(fēng) 2026-03-29 | 15/750 |
|
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +4 | 大饞小子 2026-03-29 | 4/200 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
|
[考研] 環(huán)境工程 085701，267求調(diào)劑 +6 | minht 2026-03-29 | 6/300 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +10 | 呼呼？~+123456 2026-03-24 | 10/500 |
|
|
|
[考研]
 085602 289分求調(diào)劑
+8
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑！�。�！ +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|