 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 958 | 回復(fù): 2 | ||||
空谷新苗木蟲 (小有名氣)
|
[求助]
用vasp做結(jié)構(gòu)參數(shù)的優(yōu)化時(shí)輸出文件中只有一組體積-能量數(shù)據(jù)，求解答。
|
|
我是參照侯老師《入門指南》中的一個(gè)腳本來算的，但輸出的EtVo.dat文件中只有一組體積-能量數(shù)據(jù)，麻煩大家看看是不是我的腳本有問題。 腳本： #!bin/sh rm WAVECAR for i in 4.09 4.19 4.29 4.39 4.49 4.59 4.69 4.79 4.89 4.99 5.09do cat > POSCAR <<! TiO2 1 4.594000000000000 0.000000000000000 0.000000000000000 0.000000000000000 4.594000000000000 0.000000000000000 0.000000000000000 0.000000000000000 2.959000000000000 O Ti 4 2 Direct 0.305000 0.305000 0.000000 -0.305000 -0.305000 0.000000 0.195000 0.805000 0.500000 0.805000 0.195000 0.500000 0.000000 0.000000 0.000000 0.500000 0.500000 0.500000 ! 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 INCAR.relax文件： SYSTEM = bulk rutile TiO2 unit cell ENCUT = 450 ISTART = 0 ICHARG = 2 ISMEAR = 1 SIGMA = 0.05 ISIF=5 IBRION = 2 NSW = 300 POTIM = 0.1 PREC = Accurate EDIFF = 1E-6 EDIFFG = 1E-4 INCAR.static文件： SYSTEM = bulk rutile TiO2 unit cell ENCUT = 450 PREC = Accurate ISTART = 0 ICHARG = 2 ISMEAR = -5 KPOINTS文件： Automatic 0 Monkhorst–Pack 5 5 8 0.0 0.0 0.0 輸出的EtVo.dat: 62.4500 -53.101836 62.4500 -53.101836 62.4500 -53.101836 62.4500 -53.101836 62.4500 -53.101836 ······ |
回復(fù)此樓» 猜你喜歡
299求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑
已經(jīng)有6人回復(fù)
-
300求調(diào)劑，材料科學(xué)英一數(shù)二
已經(jīng)有8人回復(fù)
-
招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
070305高分子化學(xué)與物理 304分求調(diào)劑
已經(jīng)有7人回復(fù)
-
289求調(diào)劑
已經(jīng)有13人回復(fù)
-
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分
已經(jīng)有9人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
調(diào)劑求院校招收
已經(jīng)有5人回復(fù)
-
調(diào)劑310
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 關(guān)于vasp中，結(jié)構(gòu)優(yōu)化時(shí)所使用的參數(shù)POTIM
已經(jīng)有3人回復(fù)
- ELK簡(jiǎn)單計(jì)算失敗，求高手解答！
已經(jīng)有4人回復(fù)
- Vasp相比經(jīng)典分子動(dòng)力學(xué)的優(yōu)缺點(diǎn)
已經(jīng)有11人回復(fù)
- 用vasp優(yōu)化CO結(jié)構(gòu)，機(jī)器不算，不知道哪出錯(cuò)了？高手來幫幫我！
已經(jīng)有29人回復(fù)
- P4VASP做pdos交流，及后續(xù)origin計(jì)算d-band center
已經(jīng)有34人回復(fù)
- in文件運(yùn)行失�。╨ost atoms）求助��！
已經(jīng)有11人回復(fù)
- 結(jié)構(gòu)優(yōu)化完結(jié)果老出現(xiàn)ZBRENT: fatal error in bracketing
已經(jīng)有30人回復(fù)
- 用VASP優(yōu)化結(jié)構(gòu)收斂不了，求指教
已經(jīng)有10人回復(fù)
- VASP結(jié)構(gòu)優(yōu)化時(shí)出現(xiàn)的POTIM問題
已經(jīng)有4人回復(fù)
- vasp優(yōu)化結(jié)構(gòu)中出現(xiàn)的問題
已經(jīng)有10人回復(fù)
- 重金請(qǐng)教采用準(zhǔn)諧近似求熱膨脹系數(shù)的方法
已經(jīng)有16人回復(fù)
- 【求助完結(jié)】vasp優(yōu)化結(jié)構(gòu)運(yùn)行一段時(shí)間后沒有出現(xiàn)錯(cuò)誤，但是停止了！
已經(jīng)有20人回復(fù)
- 【求助】VASP怎么做晶格參數(shù)的結(jié)構(gòu)優(yōu)化
已經(jīng)有17人回復(fù)
- 【求助】vasp中如何做結(jié)構(gòu)優(yōu)化？
已經(jīng)有10人回復(fù)
- 【求助成功】VASP中，如何固定磁矩做結(jié)構(gòu)優(yōu)化
已經(jīng)有14人回復(fù)
- 【討論】交流討論下用VASP計(jì)算不同溫度下材料的熱力學(xué)性能
已經(jīng)有19人回復(fù)
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3514.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
★
ljw4010: 金幣+1, 謝謝交流！ 2015-03-27 20:34:25
ljw4010: 金幣+1, 謝謝交流！ 2015-03-27 20:34:25
|
在POSCAR需要讓晶格常數(shù)成為變量，這里我只好固定一下c/a的值。 。。。。。。 rm WAVECAR for i in 4.09 4.19 4.29 4.39 4.49 4.59 4.69 4.79 4.89 4.99 5.09do cat > POSCAR <<! TiO2 1 $i 0.000000000000000 0.000000000000000 0.000000000000000 $i 0.000000000000000 0.000000000000000 0.000000000000000 $[$i*0.6441] 。。。。。。。 不過 $[$i*0.6441]這樣的表達(dá)我不確定行不行，如果不行的話你可以找找linux腳本相關(guān)的資料進(jìn)行修改。 |
空谷新苗2樓2015-03-23 11:10:50
空谷新苗
木蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 3797.8
- 散金: 182
- 紅花: 5
- 帖子: 238
- 在線: 158.8小時(shí)
- 蟲號(hào): 2218173
- 注冊(cè): 2013-01-02
- 性別: GG
- 專業(yè): 光學(xué)
3樓2015-03-23 15:25:38
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 空谷新苗 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +6 | mikasa的圍巾 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +9 | 辰沐5211314 2026-03-26 | 9/450 |
|
|
|
[考研] 調(diào)劑310 +7 | 溫柔的晚安 2026-03-25 | 8/400 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +8 | 張zz111 2026-03-27 | 9/450 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書 2026-03-23 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|