原子間距離と面間距離アルミニウム結晶は面心立方格子でありその格子常数(立方体の単位格子をとったときの一辺)は0.40496nm(4.0496Å)である。立方体の形にとると、1の単位格子に4つの原子が含まれている。1の単位格子に1つの原子が含まれるように単位格子をとるには、のように単位格子ベクトルをとればよい。0.0))box1.append(atom1)のように定義される。電子配置アルミニウムの原子番号は13であるから、当然電子の個数は13である。この13個の電子がのような電子配置をとっている。3p状態が中途半端につまっているということは、主量子数が3の状態すなわち3sと3pに入っている3個の電子が価電子、主量子数が3未満すなわち1sと2sと2pに入っている10個の電子が内殻電子ということになる。原子番号の大きな原子、すなわちd状態やf状態が関係してくる場合にはこのように単純には言えないが、Alの場合は単純にこの3個の価電子に対する電子状態だけを計算しておけばよい。スピンの重複度を考えると一つの状態に二つの電子が入れるので、価電子数3の半分を切り上げた2個の状態を最低限計算する必要がある。ところで、結晶の価電子の状態は当然原子の価電子の状態である3sと三つの3p状態、合わせて4個の原子状態が混成した状態であろうとも考えられる。どのような混成が起こっているかわからない以上、一応4個の電子状態を計算しておくのがよいと思われる。4)のように定義する。後はbox1とsolver1を結びつけて、全エネルギーを計算する。一酸化炭素のような分子の計算の場合、周期的に並べられた箱(Atoms)と箱の間には電子の行き来は無いと考えられるし、また行き来は無いと考えられる程度に箱を大きくしなければならないが、結晶の場合はそもそも箱は結晶の単位格子でしかないから、電子は隣り合う箱の間を行き来する。「運動量」を持っており、いろいろな運動量についてエネルギーを計算する必要がある。量子力学では運動量は波数kで表されるのでK-Pointなどとも呼ばれる。正確にはブリリュアンゾーン(BZ)の範囲内の連続的な波数について計算する必要があるが、これを有限個のK点で近似する。Monkhorst-PackのようにJacapoのコンストラクタにkptsキーワードで指定する。この例ではBZ内に2x2x2=8個のK点をとってその運動量を持った電子のエネルギーを計算する。実際にどのようなK点について計算しているかについては実行結果のファイルの内容を確かめよ。1))即ちブリリュアンゾーン内のΓ点についての計算を行う。カットオフ箱(単位格子)を周期的に並べる事により電荷密度をフーリエ展開し、同時に波動関数もフーリエ展開してハミルトニアンの行列要素を計算するが、フーリエ展開は無限級数であるので計算機で取り扱うためにはこれをどこかで打ち切らなければならない。打ち切る波数をエネルギーに換算したものをカットオフエネルギーと呼び、次のようにJacapoのコンストラクタに電荷密度についてカットオフエネルギーはdwと波動関数についてカットオフエネルギーはpwというキーワードで設定する。出力ファイルに指定したファイル名のファイルが既にあるときは上書きされてしまう!詳細は別項を参照。正常終了したプログラムの出力結果のncファイルが存在するディレクトリで次のコマンドを実行する。output_text_file.textCyberDuck等でoutput_text_file.textをMac側に転送する。output_text_file.textの中にはエネルギーと状態密度がx-yグラフを描く形式で収められている。これは単なるテキストデータなので適当に加工することも可能。Abscissa等でoutput_text_file.textのグラフをつくる。演習問題PlaneWaveCutoffの値と全エネルギーの関係をグラフにせよ。(3,3,3)...と変えながら全エネルギーと状態密度の変化を調べよ。全エネルギーの変化(縦軸)のグラフを作成せよ。kimiソースの表示以前のリビジョン文書の先頭へ。