 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (5234)
- >考研 (919)
- >導(dǎo)師招生 (852)
- >文獻(xiàn)求助 (332)
- >蟲(chóng)友互識(shí) (231)
- >考博 (138)
- >碩博家園 (106)
- >基金申請(qǐng) (90)
- >博后之家 (84)
- >論文投稿 (74)
- >招聘信息布告欄 (63)
- >第一性原理 (61)
- >教師之家 (56)
- >休閑灌水 (49)
- >公派出國(guó) (32)
- >外文書(shū)籍求助 (26)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1329 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
zhang668木蟲(chóng) (正式寫(xiě)手)
|
[交流]
【求助】做狀態(tài)方程計(jì)算時(shí)候得到的EtVo.dat是空的 已有4人參與
|
||
|
下面是我按照侯博士教程寫(xiě)的腳本： #!/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ì)算。請(qǐng)問(wèn)腳本那里出了問(wèn)題？ 還有如果不想得到體模量，只想優(yōu)化，什么時(shí)候才有做狀態(tài)方程擬合晶格參數(shù)必要，一般直接做全優(yōu)化可以嗎？還有做狀態(tài)方程的時(shí)候，沒(méi)有優(yōu)化自由度，OUTCAR最后FORCE比較大，這個(gè)能量也不能說(shuō)明問(wèn)題吧 |
回復(fù)此樓» 猜你喜歡
299求調(diào)劑
已經(jīng)有7人回復(fù)
-
0703化學(xué)調(diào)劑，求導(dǎo)師收
已經(jīng)有4人回復(fù)
-
300求調(diào)劑，材料科學(xué)英一數(shù)二
已經(jīng)有3人回復(fù)
-
材料學(xué)碩，求調(diào)劑
已經(jīng)有6人回復(fù)
-
307求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩331求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）
已經(jīng)有7人回復(fù)
-
一志愿C9材料與化工專業(yè)總分300求調(diào)劑
已經(jīng)有4人回復(fù)
-
321求調(diào)劑
已經(jīng)有4人回復(fù)
-
307求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 如何求解描述振動(dòng)的二階微分方程
已經(jīng)有21人回復(fù)
- 橢圓形偏微分方程標(biāo)準(zhǔn)化后該如何解？
已經(jīng)有6人回復(fù)
- 有關(guān)磁流體方程計(jì)算的詢問(wèn)
已經(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語(yǔ)言圖書(shū)-求解差分方程和并行計(jì)算！
已經(jīng)有43人回復(fù)
- 求助這個(gè)非線性方程改怎么解��？？
已經(jīng)有4人回復(fù)
- 求教origin擬合簡(jiǎn)單的滲流方程（又稱逾滲方程），做復(fù)合材料的應(yīng)該會(huì)遇到吧
已經(jīng)有20人回復(fù)
- 【求助】如何計(jì)算實(shí)驗(yàn)數(shù)據(jù)與已知方程的相關(guān)系數(shù)？
已經(jīng)有6人回復(fù)
zhang668
木蟲(chóng) (正式寫(xiě)手)
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 11 (小學(xué)生)
- 貴賓: 0.061
- 金幣: 1735
- 散金: 1294
- 紅花: 11
- 帖子: 843
- 在線: 324.9小時(shí)
- 蟲(chóng)號(hào): 733242
- 注冊(cè): 2009-03-28
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
9樓2010-04-24 11:45:33
gavinliu7390
木蟲(chóng) (著名寫(xiě)手)
葉落鷹飛
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時(shí)
- 蟲(chóng)號(hào): 675997
- 注冊(cè): 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
zhang668(金幣+3): 2010-04-24 10:16
|
別人的腳本，不要拿過(guò)來(lái)就用，最好試一試。 你可以找一個(gè)算完的OUTCAR試一試，確定腳本沒(méi)問(wèn)題了，再進(jìn)行計(jì)算。 EtVo.dat是空的，就說(shuō)明是E和V的變量沒(méi)有賦進(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)方程絕對(duì)不是為了獲得B0， 大部分的目的還是認(rèn)為壓力算的不準(zhǔn)，不如擬合物態(tài)方程算出的壓力準(zhǔn)。至于擬合物態(tài)方程也緊緊限于立方結(jié)構(gòu)，非立方結(jié)構(gòu)由于加壓之后，很可能晶格常數(shù)的比例變了，必須全優(yōu)化后，才能擬合物態(tài)方程，也有不擬合的。 至于你說(shuō)的FORCE很大， 那最大的原因就是你沒(méi)有relax離子，否則不可能很大，或者你的收斂參數(shù)不夠高。對(duì)應(yīng)的就是isif=2， nsw=60 立方體系也可以全優(yōu)化，但是那樣比較浪費(fèi)時(shí)間，沒(méi)有這樣的快！ |
2樓2010-04-23 23:35:28
wangzg9292
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3933
- 散金: 56
- 帖子: 348
- 在線: 89.6小時(shí)
- 蟲(chóng)號(hào): 723844
- 注冊(cè): 2009-03-16
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2010-04-24 09:05:33
zhang668
木蟲(chóng) (正式寫(xiě)手)
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 11 (小學(xué)生)
- 貴賓: 0.061
- 金幣: 1735
- 散金: 1294
- 紅花: 11
- 帖子: 843
- 在線: 324.9小時(shí)
- 蟲(chóng)號(hào): 733242
- 注冊(cè): 2009-03-28
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2010-04-24 11:08:41
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 305分求調(diào)劑（食品工程） +5 | Sxy112 2026-03-21 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 0856材料專碩353求調(diào)劑 +5 | NIFFFfff 2026-03-20 | 5/250 |
|
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +6 | holy J 2026-03-21 | 6/300 |
|
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
|
[考研] 化學(xué)308分求調(diào)劑 +3 | 你好明天你好 2026-03-23 | 3/150 |
|
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +7 | 生物工程調(diào)劑 2026-03-17 | 12/600 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 9/450 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|