AkaiKKRの使い方

AkaiKKR(machikayema)は、東京大学の赤井久純先生が開発された第一原理計算パッケージです。結晶の電子状態を計算することのできる第一原理計算パッケージは、多くの種類が公開されていますが、AkaiKKRはKKR-CPA法を採用することで化学的な不規則(不純物、不規則合金、固溶体)のバンド計算が可能であるというのが最大の特徴です。また、入力ファイルも単純であり、スペックの低い旧式のwindowsノートパソコンでも動作するといったところも優れています。

今回のエントリでは、これまでに書いたAkaiKKR関連のエントリを整理してリンクします。


1. はじめに

この章では、AkaiKKRのインストール方法とテスト計算、及び、使い方の最初の一歩についてまとめます。AkaiKKRは、公式なマニュアルであるKKR-Green関数法によるバンド計算Machikaneyama2000の使用に関するメモ及び、計算機マテリアルデザイン入門 (大阪大学新世紀レクチャー)密度汎関数法の発展 -マテリアルデザインへの応用を用いて独習することも可能ですが、やはり年に二回開催されているCMD®ワークショップに参加するのが最大の近道です。

2. 基本的な使い方

第一原理計算の出力結果として基本的なものは、全エネルギー、状態密度、バンド分散です。第2章では、面心立方構造(fcc)、体心立方構造(bcc)、六方最密充填構造(hcp)、ダイヤモンド構造、そして立方晶ペロフスカイト構造に関して状態密度とブロッホスペクトル関数(バンド分散)を出力する入力ファイルの例を示します。

3. 応用的な計算

第2章が基礎編とするならば、第3章は応用編となります。例えば平衡格子定数を計算したい場合は、格子定数を変化させながら、全エネルギーが最小となる点を探します。

4. 設定パラメータ

第一原理計算は、格子定数と原子番号以外には入力パラメータが必要ないということになっていますが、実際には計算の精度に関するパラメータなどを指定しなければなりません。この章ではパラメータの意味と決定方法について書きます。

5. その他

6. 参考URL

6.1 AkaiKKR関連

6.2 第一原理計算全般

7. 参考書籍


8. データベース

9. あとがきと更新履歴

私は、AkaiKKRの開発者である赤井先生、及び、大阪大学の固体電子論グループとは何の関係も無い一個人であり、自分の固体物理の勉強のためにAkaiKKRを利用させていただこうとしている身です。更に言うなら、バンド計算屋さんではなくただの電子工作趣味人です。従って、このブログの内容の正しさに関しては一切保障できませんので、真似をされる方は自己責任でお願いします。おそらくブログ内の記述には、間違っている点もたくさんあると思います。間違いを見つけられた方は、該当エントリのコメント欄にてお知らせください。


フィードバック

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

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

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

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


tag: AkaiKKR machikaneyama KKR CPA  

comment

Secret

AkaiKKRについて

はじめまして。
趣味で固体物理の勉強をしようと思い、第一原理計算のソフト(フリー)を探しているうちに、gomisaiさんのページにたどり着きました。
gomisaiさんは固体物理の勉強目的で、AkaiKKRを選ばれたそうですが、このソフトを選ばれた理由を教えて頂ければと思います。

Re: AkaiKKRについて

あたさん、こんにちは。

私がAkaiKKRを選んだのは、色々な要素があります。

教科書的に言うならAkaiKKRが、私の興味のある物質の計算に適したコードであったというのが理由のひつとです。

ご存知かもしれませんが、第一原理計算ソフトにはそれぞれに特徴があり、適する用途と適さない用途があります。

このエントリにも書いてありますが、AkaiKKRはコヒーレントポテンシャル近似(CPA)という、化学的な不規則を、単位セルで計算することができる(スーパーセルを必要としない)という珍しい特徴を持っています。私は不規則合金に興味があったので、AkaiKKRを選んだのは結果的に大正解でした。

逆にAkaiKKRに適さない計算もあります。
AkaiKKRの使うKKR法は、全電子法の一種ですが、残念ながらフルポテンシャルではありません。マフィンティン近似か原子球近似(ASA)のどちらかを利用します。これらは、金属のよう最密充填に近い構造ではよい近似ですが、グラファイトのように結晶構造が「すかすか」な結晶の計算には適さないようです。
そして有限温度の物性を見積もるのにも向いていないようです。
第一原理分子動力学(MD)はもちろんのこと、密度汎関数摂動論(DFPT)でフォノンのバンド構造を計算することもできません。
ただ、これらの問題もまったく手が出ないわけではなくAkaiKKRでグラファイトのバンド構造( http://gomisai.blog75.fc2.com/blog-entry-641.html )やAkaiKKRで金属の熱物性( http://gomisai.blog75.fc2.com/blog-entry-646.html )などで挑戦はしています。



教科書的な理由のほかに、現実的な理由として、日本語のマニュアルがあるコードの中で比較的インストールや使用が簡単であるという点があります。
色々なコードは、それぞれ物性値を計算するために便利な機能があったりしますが、やはり基本となるのは「全エネルギー」「バンド構造」と「状態密度」です。これらは、おおよそほとんどどのような第一原理バンド計算ソフトでも計算できるはずなので、基本的な計算をすることが目的なら、機能の豊富さよりもインストールと使用が簡単なものを選ぶのがよい事になります。
なんといっても年2回開催されているCMDワークショップ( http://phoenix.mp.es.osaka-u.ac.jp/CMD/ )で使い方を習うことができるのが大きいです。私もすでに2回ほど参加させていただいています。(第23回CMDワークショップ参加報告 http://gomisai.blog75.fc2.com/blog-entry-565.html )
またCMDワークショップに参加せずとも、インストールや実行のやり方にはほとんど悩まないはずです。コンパイラはgfortranでもifortranでもmakefileを書き換えるだけでコンパイルできます。(つまづくとしたらopenmpの設定のあたりぐらいと思います。)

No title

AkaiKKR選定理由の他、第一原理計算全般についていろいろとご教示頂き、ありがとうございました。

第一原理計算およびソフトについては全くと言っていいほど知らないため、前半の部分はとても参考になりました。とはいっても、書かれている内容を本当に理解するまでには、時間がかかりそうだと思いますが・・・。

使う立場にとっては、そのソフトが実用的なものか(きちんと計算が出来るものなのか)、マニュアルが整備されてあるのかがとても重要なのですが、gomisaiさんの説明でAkaiKKRが安心して使えるものだと納得できました。
第一候補として検討してみたいと思います。

Re: No title

あたさん、こんにちは。

ねがてぃぶろぐの記事が、AkaiKKRのチュートリアルとして役に立てば幸いです。
私も初心者なので内容が怪しい記事もあるかと思いますが、何か疑問点や気づいた点があれば、その記事にコメントをしてください。

また、どの計算コードを選ぶにしても、最初の1種類に慣れれば、他のコードのセットアップの方法も理解しやすくなるかと思います。実際私も最近ではフォノンの計算をすることを考えてQuantum ESPRESSOを触り始めました。
FC2カウンター
カテゴリ
ユーザータグ

LTspiceAkaiKKRmachikaneyamaScilabKKRPSoC強磁性CPAPICOPアンプecalj常微分方程式モンテカルロ解析状態密度トランジスタodeDOSインターフェース定電流スイッチング回路PDS5022半導体シェルスクリプト分散関係レベルシフト乱数HP6632AR6452A可変抵抗トランジスタ技術温度解析ブレッドボードI2C反強磁性確率論数値積分セミナーバンドギャップバンド構造偏微分方程式非線形方程式ソルバ熱設計絶縁ISO-I2Cカオス三端子レギュレータLM358GW近似マフィンティン半径A/DコンバータフォトカプラシュミットトリガLEDPC817C発振回路数値微分直流動作点解析サーボカレントミラーTL431アナログスイッチUSB74HC4053bzqltyVESTA補間電子負荷アセンブライジング模型BSch量子力学単振り子2ちゃんねるチョッパアンプLDA開発環境基本並進ベクトルFFT標準ロジックブラべ格子パラメトリック解析抵抗SMPMaxima失敗談ラプラス方程式繰り返し位相図スイッチト・キャパシタ熱伝導状態方程式キュリー温度gfortranコバルトTLP621不規則合金Quantum_ESPRESSO六方最密充填構造ランダムウォーク相対論ewidthスピン軌道相互作用FETQSGWVCAcygwinスレーターポーリング曲線GGA仮想結晶近似PWscfシュレディンガー方程式LM555ハーフメタル固有値問題NE555最小値ガイガー管QNAPUPS自動計測ダイヤモンドマントルTLP552格子比熱最適化MCU井戸型ポテンシャル最大値xcrysdenCIF条件分岐詰め回路フェルミ面差し込みグラフスーパーセルfsolveブラウン運動awk過渡解析起電力三角波第一原理計算FXA-7020ZRWriter509Ubuntuテスタ熱力学データロガーTLP521OpenMPubuntu平均場近似MAS830LトランスCK1026PIC16F785PGA2SC1815EAGLEノコギリ波負帰還安定性ナイキスト線図MBEOPA2277P-10フィルタCapSenseAACircuitLMC662文字列固定スピンモーメントFSMTeX結晶磁気異方性全エネルギーc/a合金multiplotgnuplot非線型方程式ソルバL10構造正規分布等高線ジバニャン方程式初期値interp1fcc面心立方構造ウィグナーザイツ胞半金属デバイ模型電荷密度重積分SIC二相共存磁気モーメント不純物問題PWgui擬ポテンシャルゼーベック係数ZnOウルツ鉱構造edeltquantumESPRESSOフォノンリジッドバンド模型スワップ領域BaO岩塩構造ルチル構造ヒストグラム確率論マテリアルデザインフラクタルマンデルブロ集合キーボードRealforceクーロン散乱三次元疎行列縮退化学反応関数フィッティング最小二乗法Excel直流解析PCTS-110TS-112日本語パラメータ・モデル等価回路モデルcif2cell入出力陰解法熱拡散方程式HiLAPW両対数グラフCrank-Nicolson法連立一次方程式specx.fifort境界条件片対数グラフグラフの分割円周率ヒストグラム不規則局所モーメントGimpシンボル軸ラベル凡例線種トラックボール

最新コメント
リンク

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