版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

043114076

木蟲 (正式寫手)

應(yīng)助: 49 (小學(xué)生)
金幣: 3364.9
散金: 52
紅花: 23
帖子: 788
在線: 369.5小時(shí)
蟲號(hào): 570592
注冊(cè): 2008-06-07
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

[交流] Vasp NEB 計(jì)算前原子排序腳本已有6人參與

使用vasp NEB計(jì)算反應(yīng)路徑一般會(huì)使用nebmake.pl 進(jìn)行插點(diǎn)生成一系列的POSCAR，常常會(huì)遇到初始結(jié)構(gòu)和最終結(jié)構(gòu)POSCAR的原子順序并不一一對(duì)應(yīng)的情況，這樣插點(diǎn)就會(huì)出現(xiàn)結(jié)構(gòu)混亂。我自己寫了個(gè)排序的腳本用于自動(dòng)排序，供大家參考使用：
使用方法：將下面的腳本拷貝到一個(gè)文件，比如pos_order.sh。根據(jù)自己的情況修改begin和end兩個(gè)變量。運(yùn)行：sh pos_order.sh POSCAR.start POSCAR.end
生成一個(gè)POSCAR.end.new。其中POSCAR.start 和POSCAR.end 是初始結(jié)構(gòu)和最終結(jié)構(gòu)，可以使用其他名字。

CODE:

#!/bin/bash

begin=9 #輸入結(jié)構(gòu)開始行號(hào)

end=133 #輸入結(jié)構(gòu)結(jié)束行號(hào)

headline=`awk "BEGIN{print $begin-1}"`

head -n $headline $2 > $2.new

awk -v begin=$begin -v end=$end '

    (FILENAME==ARGV[1])&&(FNR>=begin)&&(FNR<=end){x0[FNR]=$1;y0[FNR]=$2;z0[FNR]=$3}

    (FILENAME==ARGV[2])&&(FNR>=begin)&&(FNR<=end){x1[FNR]=$1;y1[FNR]=$2;z1[FNR]=$3}

     END{for(i=begin;i<=end;i=i+1){

          rmin=1;

          for(j=begin;j<=end;j=j+1){

            r=(x0[i]-x1[j])^2+(y0[i]-y1[j])^2+(z0[i]-z1[j])^2;

            if(r<rmin){rmin=r;nr=j}

            };

          x2[i]=x1[nr];y2[i]=y1[nr];z2[i]=z1[nr];

          printf("%9.6f%9.6f%9.6f\n",x2[i],y2[i],z2[i])

          };

        }' $1 $2 >> $2.new

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

量化軟件	VASP軟件從入門到高深

» 本帖已獲得的紅花（最新10朵）

afan0012

wangscott

» 猜你喜歡

288求調(diào)劑，一志愿華南理工大學(xué)071005 已經(jīng)有5人回復(fù)
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有12人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有5人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有13人回復(fù)
梁成偉老師課題組歡迎你的加入已經(jīng)有11人回復(fù)
復(fù)試調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷已經(jīng)有11人回復(fù)
321求調(diào)劑已經(jīng)有10人回復(fù)
東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) 已經(jīng)有8人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

VASP CI-NEB找過渡態(tài)求助已經(jīng)有4人回復(fù)
關(guān)于Vasp TST Tools的NEB計(jì)算中的SS NEB計(jì)算問題已經(jīng)有7人回復(fù)
vasp結(jié)構(gòu)優(yōu)化出現(xiàn)問題：算完一個(gè)離子步后就停止已經(jīng)有12人回復(fù)
VASP NEB 不收斂問題已經(jīng)有25人回復(fù)
vasp中做利用NEB或者 Cl-NEB做動(dòng)力學(xué) 需要怎樣的計(jì)算機(jī)配置呢？已經(jīng)有14人回復(fù)
vasp 計(jì)算單分子/原子時(shí)晶胞大小的選取已經(jīng)有6人回復(fù)
用VASP計(jì)算NEB都十幾天了還不收斂，急�。。。� 已經(jīng)有8人回復(fù)
vasp計(jì)算NEB時(shí)，算到第27步，自動(dòng)停了，高人幫我看看怎么回事？已經(jīng)有7人回復(fù)
VASP計(jì)算neb的步驟已經(jīng)有14人回復(fù)
vasp neb 加壓，如何修改INCAR ? 已經(jīng)有12人回復(fù)
請(qǐng)教VASP如何在NEB計(jì)算之后做指前因子計(jì)算？已經(jīng)有7人回復(fù)
vasp計(jì)算neb忽然停掉了！如果想繼續(xù)算我應(yīng)該怎么辦吶？已經(jīng)有8人回復(fù)
畫VASP計(jì)算的原子結(jié)構(gòu)圖的一個(gè)方法（自己瞎想的，如和他人方法雷同請(qǐng)見諒）已經(jīng)有4人回復(fù)
vasp計(jì)算界面結(jié)構(gòu)時(shí)原子弛豫的問題已經(jīng)有12人回復(fù)
vasp計(jì)算過渡態(tài)問題已經(jīng)有13人回復(fù)
vasp用ci-neb算出來的過渡態(tài)有多個(gè)虛頻，怎么辦？已經(jīng)有15人回復(fù)
新手請(qǐng)教用NEB方法計(jì)算energy barrier. 已經(jīng)有15人回復(fù)
vasp用neb計(jì)算遷移勢壘已經(jīng)有11人回復(fù)
vasp NEB 如何選取原胞大小？已經(jīng)有10人回復(fù)
【求助】新手求助vasp計(jì)算氧分子結(jié)合能的問題已經(jīng)有12人回復(fù)
【求助】關(guān)于vasp NEB 已經(jīng)有8人回復(fù)
【其他】vasp neb 計(jì)算討論已經(jīng)有20人回復(fù)
【求助成功】NEB計(jì)算的一個(gè)stdout輸出怎么解讀? 已經(jīng)有6人回復(fù)
【求助成功】怎么樣提高NEB計(jì)算速度? 已經(jīng)有19人回復(fù)

1樓 2014-07-28 11:07:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dujunhui

金蟲 (正式寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 1769.3
帖子: 403
在線: 136.3小時(shí)
蟲號(hào): 739582
注冊(cè): 2009-04-03
性別: GG
專業(yè): 無機(jī)材料化學(xué)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

這么好的帖子，怎么沒人頂

贊一下

回復(fù)此樓

2樓2014-09-15 11:10:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zb3ynf

銀蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 390.9
帖子: 39
在線: 57.6小時(shí)
蟲號(hào): 1270709
注冊(cè): 2011-04-19
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

2樓: Originally posted by dujunhui at 2014-09-15 11:10:04
這么好的帖子，怎么沒人頂

無比感謝�。。。。。。。。。。。。�

回復(fù)此樓

3樓2014-09-20 16:19:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

呂小強(qiáng)強(qiáng)

新蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 6498.9
紅花: 1
帖子: 866
在線: 237.7小時(shí)
蟲號(hào): 2978560
注冊(cè): 2014-02-19
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

頂了不過動(dòng)力學(xué)的 FFF TTT 都變沒了有影響嗎

贊一下

回復(fù)此樓

4樓2015-02-07 09:22:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

043114076

木蟲 (正式寫手)

應(yīng)助: 49 (小學(xué)生)
金幣: 3364.9
散金: 52
紅花: 23
帖子: 788
在線: 369.5小時(shí)
蟲號(hào): 570592
注冊(cè): 2008-06-07
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

4樓: Originally posted by 呂小強(qiáng)強(qiáng) at 2015-02-07 09:22:05
頂了不過動(dòng)力學(xué)的 FFF TTT 都變沒了有影響嗎

如果設(shè)定了Selective Dyanamic, 是有影響的. 不過稍微修改一下腳本就可以實(shí)現(xiàn).

贊一下

回復(fù)此樓

5樓2015-02-09 08:57:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

呂小強(qiáng)強(qiáng)

新蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 6498.9
紅花: 1
帖子: 866
在線: 237.7小時(shí)
蟲號(hào): 2978560
注冊(cè): 2014-02-19
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

引用回帖:

5樓: Originally posted by 043114076 at 2015-02-09 08:57:37
如果設(shè)定了Selective Dyanamic, 是有影響的. 不過稍微修改一下腳本就可以實(shí)現(xiàn)....

挺好使，只是有些體系為什么產(chǎn)生的并不一一對(duì)應(yīng)而且POSCAR.NEW轉(zhuǎn)化成GIF后有些原子不見了，potcar估計(jì)也難對(duì)上那順序，是不是那種路徑不存在才會(huì)這樣。

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

6樓2015-02-10 00:00:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

抖抖豆豆

7樓2016-01-24 16:22:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

afan0012

木蟲 (小有名氣)

應(yīng)助: 7 (幼兒園)
金幣: 4328
散金: 7
紅花: 6
帖子: 228
在線: 65.2小時(shí)
蟲號(hào): 3028351
注冊(cè): 2014-03-08
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

送紅花一朵

多謝分享

回復(fù)此樓

8樓2016-01-25 19:44:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangscott

送紅花一朵

謝謝分享

9樓2022-05-27 22:55:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 043114076 的主題更新

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	6/300	2026-03-19 17:18 by fei626-918
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +8	zlingli 2026-03-13	8/400	2026-03-19 16:32 by 輕松不少隨
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 材料專碩英一數(shù)二306 +5	z1z2z3879 2026-03-18	5/250	2026-03-19 07:43 by BruceLiu320
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu

24小時(shí)熱門版塊排行榜

[交流] Vasp NEB 計(jì)算前原子排序腳本 已有6人參與

» 收錄本帖的淘帖專輯推薦

» 本帖已獲得的紅花（最新10朵）

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

[交流] Vasp NEB 計(jì)算前原子排序腳本已有6人參與