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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

043114076

木蟲 (正式寫手)

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

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

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

CODE:

#!/bin/bash

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

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

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

» 猜你喜歡

求調(diào)劑有機化學考研356分已經(jīng)有4人回復(fù)
0703 化學求調(diào)劑，一志愿山東大學 342 分已經(jīng)有4人回復(fù)
086000生物與醫(yī)藥調(diào)劑已經(jīng)有7人回復(fù)
305求調(diào)劑已經(jīng)有7人回復(fù)
348求調(diào)劑已經(jīng)有3人回復(fù)
332求調(diào)劑已經(jīng)有6人回復(fù)
各位老師好，我的一志愿為北京科技大學085601材料專碩已經(jīng)有6人回復(fù)
085600 材料與化工 329分求調(diào)劑已經(jīng)有10人回復(fù)
356求調(diào)劑已經(jīng)有3人回復(fù)
一志愿武漢理工，總分321，英一數(shù)二，求老師收留。已經(jīng)有7人回復(fù)

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

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

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

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

043114076

木蟲 (正式寫手)

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

引用回帖:

4樓: Originally posted by 呂小強強 at 2015-02-07 09:22:05
頂了不過動力學的 FFF TTT 都變沒了有影響嗎

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

贊一下

回復(fù)此樓

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

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

查看全部 9 個回答

dujunhui

金蟲 (正式寫手)

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

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

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

贊一下

回復(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小時
蟲號: 1270709
注冊: 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的回帖

呂小強強

新蟲 (正式寫手)

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

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

頂了不過動力學的 FFF TTT 都變沒了有影響嗎

贊一下

回復(fù)此樓

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

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

查看全部 9 個回答

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

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

[考研] 調(diào)劑310 +6	溫柔的晚安 2026-03-25	7/350	2026-03-28 16:04 by qq2000
[考研] 0703本科鄭州大學求調(diào)劑 +3	nhj_ 2026-03-25	3/150	2026-03-28 13:24 by Iveryant
[考研] 求調(diào)劑 +6	蘆lty 2026-03-25	7/350	2026-03-28 13:13 by 唐沐兒
[考研] 070300求調(diào)劑306分 +4	26要上岸 2026-03-27	4/200	2026-03-28 13:06 by 唐沐兒
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6	烤麥芽 2026-03-27	8/400	2026-03-28 10:40 by 烤麥芽
[考研] 311求調(diào)劑 +3	希望上岸阿小楊 2026-03-23	3/150	2026-03-28 07:57 by 熱情沙漠
[考研] 287求調(diào)劑 +10	land xuxu 2026-03-26	10/500	2026-03-27 15:33 by 帕爾馬拉特
[考研] 316求調(diào)劑 +5	Pigcasso 2026-03-24	5/250	2026-03-27 12:10 by zhshch
[考研] 考研調(diào)劑 +10	呼呼？~+123456 2026-03-24	10/500	2026-03-27 11:46 by wangjy2002
[考研] 求調(diào)劑 +3	劉柯@ 2026-03-24	4/200	2026-03-27 11:28 by shangxh
[考研] 312求調(diào)劑 +9	上岸吧ZJY 2026-03-22	13/650	2026-03-27 11:24 by sanrepian
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 0703化學338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 325求調(diào)劑 +3	Aoyijiang 2026-03-23	3/150	2026-03-26 20:46 by 不吃魚的貓
[考研] 化學調(diào)劑一志愿上海交通大學336分-本科上海211 +4	小魚愛有機 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 0854電子信息求調(diào)劑 +7	α____ 2026-03-22	9/450	2026-03-25 13:37 by α____
[考研] 292求調(diào)劑 +4	鵝鵝鵝額額額額?/a> 2026-03-24	4/200	2026-03-24 16:41 by peike
[考研] 求調(diào)劑一志愿武漢理工大學材料工程（085601） +5	WW.' 2026-03-23	7/350	2026-03-24 14:50 by sprinining
[考研] 085404電子信息284分求調(diào)劑 +4	13659058978 2026-03-24	4/200	2026-03-24 12:15 by syl20081243
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

043114076

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

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

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

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

043114076

dujunhui

zb3ynf

呂小強強

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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

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