#!/bin/csh -f ## *** プロジェクト名 *** set PROJECT="HalfMetal" ## ポテンシャルファイル名 set POTENTIAL=${PROJECT} ## *** 格子定数のリスト (Angstrom) *** set ALATT_LIST=( 4.98 5.45 5.65 5.87 6.06 6.10 6.48 ) set AATOM_LIST=( Cr Mn ) set AANUM_LIST=( 24 25 ) set BATOM_LIST=( P As Sb ) set BANUM_LIST=( 15 33 51 ) ## *** 第一原理計算 **** set i=0 foreach AATOM ( ${AATOM_LIST} ) set i=`echo "${i}+1" | bc -l` set AANUM=${AANUM_LIST[$i]} set j=0 foreach BATOM ( ${BATOM_LIST} ) set j=`echo "${j}+1" | bc -l` set BANUM=${BANUM_LIST[$j]} foreach ALATT ( ${ALATT_LIST} ) set ABOHR=`echo "scale=5; ${ALATT}/0.52917721092" | bc -l` ## 強磁性用入力ファイルの作成 sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}.in | sed 's/'ALATT'/'${ALATT}'/g' | sed 's/'AATOM'/'${AATOM}'/g' | sed 's/'AANUM'/'${AANUM}'/g' | sed 's/'BATOM'/'${BATOM}'/g' | sed 's/'BANUM'/'${BANUM}'/g' > in/${AATOM}${BATOM}${ALATT}.in ## LMD用入力ファイルの作成 sed 's/'ABOHR'/'${ABOHR}'/g' template/${PROJECT}-lmd.in | sed 's/'ALATT'/'${ALATT}'/g' | sed 's/'AATOM'/'${AATOM}'/g' | sed 's/'AANUM'/'${AANUM}'/g' | sed 's/'BATOM'/'${BATOM}'/g' | sed 's/'BANUM'/'${BANUM}'/g' > in/${AATOM}${BATOM}${ALATT}-lmd.in ## 強磁性状態の計算 specx out/${AATOM}${BATOM}${ALATT}.out ## LMD初期ポテンシャルの作成 cp data/${AATOM}${BATOM}${ALATT} data/HalfMetal fmg < HalfMetal.fmg cp data/HalfMetal_lmd data/${AATOM}${BATOM}lmd${ALATT} ## LMD状態の計算 specx out/${AATOM}${BATOM}${ALATT}-lmd.out end end end