NEB的POSCAR处理

冰海里的一条鱼

                好主意，先支持一个
官方ftp上的interpolatePOSCAR
script功能简介：将两个POSCAR线性插值得到n个POSCAR用于NEB计算
script:
#
# interpolate two POSCAR files
# POSCAR files must be concated before passing them to  
# this awk script
# first line of the first POSCAR file should be
# rep number_of_replicase
# if there is any line containing the term center
# differences in the center of mass will be removed
#  
file=$1
if [ ! -x $file ]
then
  usage: interpolatePOS POSCAR1_POSCAR2
fi
awk = 8 )  
          {  
              x2[line-7] = $1; y2[line-7] = $2 ; z2[line-7] = $3  }
             if (line==lines) {
                   print "second set read"
             }  
    }
}
END  {
   lines=lines-7
   for ( line=1; linefile
       for ( line=1; linefile
       }
   }
}'
script用法：./interpolatePOSCAR [file]
interpolatePOSCAR 为script的文件名，file为要处理的POSCAR
第一行写成rep n的格式，把要插值的两个POSCAR文件写到一个file里面
注意要建一些目录00 01 02 ...0n，否则会出错
PS：我的fedroa5.0运行不成功，不过我的ubuntu9。04运行成功了。开心阿。
               
               

本文来自ChinaUnix博客，如果查看原文请点：http://blog.chinaunix.net/u3/104783/showart_2089125.html