AkaiKKRのspc計算の入出力フォーマット

AkaiKKR(machikaneyama)は26 August 2015からブロッホスペクトル関数の入出力形式が変更されました。
簡単にバンド図が描けるようになり便利なのですが、一方で、これまでのほうが便利だったケースもありえます。
今回は、これを22 May 2015のバージョンと同じフォーマットに戻す方法について書きます。

入力フォーマットは source/spmain.f でサブルーチン readk を呼び出す際の引数 iread を1に指定します。
出力フォーマットは source/spmain.f でサブルーチン wrtspc を呼び出している行をコメントアウトしてしまい、代わりに22 May 2015で出力をしている部分のソースコードをコピペします。


AkaiKKRのブロッホスペクトル関数


AkaiKKR(machikaneyama)は通常のバンド図を出力する代わりにブロッホスペクトル関数を出力します。cpa2002v009cの22 May 2015までは、計算するk点をすべて手動で入力しなければなりませんでした。そのためAkaiKKRでSrTiO3ペロフスカイトなどではgnuplotでバンド構造をプロットするための補助的なScilabスクリプトを使っています。

これに対して26 Augst 2015では、他の一般的なDFTパッケージで見られるように、特徴的なk点の座標とその分割数を指定することで、バンド図を描くことができるようになりました(参考: AkaiKKRでバンド構造(分散関係))。

便利になった反面、特徴的なk点をつなぐパスに沿った計算以外を行う場合、例えばAkaiKKRで銅と銅亜鉛合金のフェルミ面のような計算はやりにくくなりました。また、それ以外にもこれまで書いたスクリプトを再利用するのにも不便です。
そんな折、26 August 2015のソースコードを編集して、以前の入出力形式に戻す方法を教えていただいたので、今回は、この方法についてまとめます。

入力ファイルの形式


spc計算のためのk点指定方法は source/spmain.f の382行目(ぐらい)にある call readk(vkp,coa,boa,nk1,nk3,nk3x,nkmx,kcrt,25,kblst,korder,3) の最後の引数の3を変更することで変えることができます。引数は1, 2, 3の三種類から選ぶことができ、その内容はsource/readk.fのコメントに書かれています。
c-----------------------------------------------------------------------
c This program read-in the k-points according either format (a),
c (b) or (c).
c
c (a) simply input all the k-point to be calculated.
c
c (b) input k-point corresponding to the Gamma, K, M, etc., with
c the number of devision of the k-mesh between that point and
c the previous point. For example, input data such as
c 0.0 0.0 0.0 0
c 1.0 0.0 0.0 50
c
c means that the line connecting (0, 0, 0) and (1, 0, 0)
c is divided into 50 pieces, and on each grid points
c including both endpoints the Bloch spectrum function
c will be calculated.
c
c (c) input k-point corresponding to the Gamma, K, M, etc., with
c the total number of k-point to be used as the first data.
c For example, input data such as
c 101
c 0.0 0.0 0.0
c 1.0 0.0 0.0
c 1.0 0.5 0.0
c
c means that the line connecting (0, 0, 0), (1, 0, 0), and
c (1, 0.5, 0) is divided into more or less equidistant 100
c pieces, and on each grid points including both endpoints
c the Bloch spectrum function will be calculated.
c
c In all formats, data will be read-in untill any data that
c do not follow the above format be met. Which format should be
c adopted is controled by a parameter iread. If iread=1,
c format (a), iread=2 for (b), and iread=3 for (c)
c All the data shoulbe be given in the unit of 2pi*(1/a, 1/b, 1/c),
c i.e. actual k vector is 2pi*(kx/a, ky/b, kz/c), where a, b, c
c are the lattice constant along x, y, and z direction.
c
c Coded by H. Akai, 4 June 2015, Tokyo
c-----------------------------------------------------------------------

今回は、すべてのk点を手動で指定する(a)に戻すので、最後の引数である iread の3を1に変更します。
c     call readk(vkp,coa,boa,nk1,nk3,nk3x,nkmx,kcrt,25,kblst,korder,3)
call readk(vkp,coa,boa,nk1,nk3,nk3x,nkmx,kcrt,25,kblst,korder,1)


出力ファイルの形式


出力形式の指定は source/wrtspc.f の中で行われますが source/readk.f ほどには、整備されていないようです。このsource/readk.fが追加されたのは26 August 2015からで22 May 2015の段階では、まだsource/wrtspc.fは作成されておらずファイルの出力も source/spmain.f の中で行われています。26 August 2015では wrtspcのサブルーチンは source/spmain.f の735行目(ぐらい)の call wrtspc(spctrl,is,e,mse,kcrt,kblst,nk3,ef) で呼び出されています。そこでこの呼び出しの部分をコメントアウトしてしまい、代わりに22 May 2015の source/spmain.f の else if(ids .eq. 4) then (748行目ぐらい)と endif (784行目ぐらい)の間をコピー&ペーストします。

こうすることによって26 August 2015を使いながら、ブロッホスペクトル関数の入出力を22 May 2015以前の形式にすることができます。

関連エントリ




参考URL




参考文献/使用機器




フィードバック



にほんブログ村 その他趣味ブログ 電子工作へ

 ↑ 電子工作ブログランキング参加中です。1クリックお願いします。


コメント・トラックバックも歓迎です。 ↓      


 ↓ この記事が面白かった方は「拍手」をお願いします。


tag: AkaiKKR machikaneyama 分散関係 

comment

Secret

FC2カウンター
カテゴリ
ユーザータグ

LTspiceAkaiKKRmachikaneyamaScilabKKRPSoCOPアンプCPA強磁性PICモンテカルロ解析常微分方程式odeトランジスタecalj状態密度DOSインターフェース定電流スイッチング回路PDS5022半導体シェルスクリプト乱数レベルシフトHP6632A温度解析分散関係I2Cトランジスタ技術R6452A可変抵抗ブレッドボードセミナーバンドギャップ数値積分確率論反強磁性偏微分方程式バンド構造絶縁熱設計非線形方程式ソルバフォトカプラシュミットトリガLEDLM358カオスISO-I2C三端子レギュレータGW近似A/Dコンバータカレントミラーアナログスイッチ数値微分マフィンティン半径TL431発振回路サーボPC817CUSB直流動作点解析74HC4053補間FFTBSch開発環境パラメトリック解析2ちゃんねるチョッパアンプ量子力学bzqlty電子負荷イジング模型LDA標準ロジックアセンブラ基本並進ベクトルブラべ格子単振り子熱伝導位相図TLP621キュリー温度繰り返し状態方程式MaximaVESTAスイッチト・キャパシタ相対論FETランダムウォークスピン軌道相互作用SMP六方最密充填構造抵抗不規則合金ewidthスレーターポーリング曲線GGAラプラス方程式cygwingfortranQSGW失敗談コバルト条件分岐TLP521テスタLM555Writer509TLP552格子比熱マントルデータロガー自動計測詰め回路ガイガー管ダイヤモンドQNAPMCUFXA-7020ZR過渡解析三角波UPSNE555固有値問題熱力学ブラウン運動フェルミ面awk起電力第一原理計算OpenMPfsolveubuntu最大値xcrysden最小値最適化仮想結晶近似VCA差し込みグラフスーパーセル井戸型ポテンシャル平均場近似シュレディンガー方程式FSMフラクタルOPA2277固定スピンモーメント2SC1815全エネルギー合金multiplotgnuplotc/aTeX結晶磁気異方性interp1ウィグナーザイツ胞初期値マンデルブロ集合疎行列面心立方構造fcc不純物問題非線型方程式ソルバフィルタL10構造PGA半金属二相共存SICZnOウルツ鉱構造BaO重積分クーロン散乱磁気モーメント電荷密度三次元CIF岩塩構造CapSenseノコギリ波デバイ模型ハーフメタル正規分布フォノンquantumESPRESSOルチル構造スワップ領域リジッドバンド模型edelt縮退キーボード軸ラベルグラフの分割凡例トラックボールPC不規則局所モーメント片対数グラフトランス両対数グラフCK1026MAS830L直流解析Excel円周率パラメータ・モデルヒストグラム日本語最小二乗法等価回路モデルGimp線種シンボルTS-110TS-112PIC16F785LMC662化学反応文字列specx.f入出力ifortマテリアルデザインヒストグラム確率論Realforce等高線ジバニャン方程式P-10Ubuntuナイキスト線図Crank-Nicolson法陰解法熱拡散方程式HiLAPWAACircuit連立一次方程式負帰還安定性境界条件EAGLEMBE関数フィッティング

最新コメント
リンク

にほんブログ村 その他趣味ブログ 電子工作へ