seminar:al結晶
差分
このページの2つのバージョン間の差分を表示します。
両方とも前のリビジョン前のリビジョン次のリビジョン | 前のリビジョン次のリビジョン両方とも次のリビジョン | ||
seminar:al結晶 [2021/06/21 14:56] – [実習2] kimi | seminar:al結晶 [2021/06/21 14:59] – [状態密度] kimi | ||
---|---|---|---|
行 59: | 行 59: | ||
====== 実習2 ====== | ====== 実習2 ====== | ||
* エネルギーカットオフの値と全エネルギーの関係をグラフにせよ。 | * エネルギーカットオフの値と全エネルギーの関係をグラフにせよ。 | ||
- | ==== エネルギーカットオフ ==== | + | ===== エネルギーカットオフ |
箱(単位格子)を周期的に並べる事により電荷密度をフーリエ展開し、同時に波動関数もフーリエ展開してハミルトニアンの行列要素を計算するが、フーリエ展開は無限級数であるので計算機で取り扱うためにはこれをどこかで打ち切らなければならない。打ち切る波数をエネルギーに換算したものをカットオフエネルギーと呼び、次のようにJacapoのコンストラクタに電荷密度についてカットオフエネルギーは'' | 箱(単位格子)を周期的に並べる事により電荷密度をフーリエ展開し、同時に波動関数もフーリエ展開してハミルトニアンの行列要素を計算するが、フーリエ展開は無限級数であるので計算機で取り扱うためにはこれをどこかで打ち切らなければならない。打ち切る波数をエネルギーに換算したものをカットオフエネルギーと呼び、次のようにJacapoのコンストラクタに電荷密度についてカットオフエネルギーは'' | ||
<code python> | <code python> | ||
行 66: | 行 66: | ||
数字の単位はeVである。 | 数字の単位はeVである。 | ||
- | ===== 実習3 ===== | + | ====== 実習3 |
* BZKPointsを(1, | * BZKPointsを(1, | ||
- | ==== K点サンプリング ==== | + | ===== K点サンプリング |
一酸化炭素のような分子の計算の場合、周期的に並べられた箱(Atoms)と箱の間には電子の行き来は無いと考えられるし、また行き来は無いと考えられる程度に箱を大きくしなければならないが、結晶の場合はそもそも箱は結晶の単位格子でしかないから、電子は隣り合う箱の間を行き来する。すなわち「運動量」を持っており、いろいろな運動量についてエネルギーを計算する必要がある。量子力学では運動量は波数kで表されるのでK-Pointなどとも呼ばれる。正確にはブリリュアンゾーン(BZ)の範囲内の連続的な波数について計算する必要があるが、これを有限個のK点で近似する。K点の個数は | 一酸化炭素のような分子の計算の場合、周期的に並べられた箱(Atoms)と箱の間には電子の行き来は無いと考えられるし、また行き来は無いと考えられる程度に箱を大きくしなければならないが、結晶の場合はそもそも箱は結晶の単位格子でしかないから、電子は隣り合う箱の間を行き来する。すなわち「運動量」を持っており、いろいろな運動量についてエネルギーを計算する必要がある。量子力学では運動量は波数kで表されるのでK-Pointなどとも呼ばれる。正確にはブリリュアンゾーン(BZ)の範囲内の連続的な波数について計算する必要があるが、これを有限個のK点で近似する。K点の個数は | ||
<code python> | <code python> | ||
行 84: | 行 84: | ||
$\vec{p}=\hbar\vec{k}=(0, | $\vec{p}=\hbar\vec{k}=(0, | ||
- | ==== 状態密度 ==== | + | ===== 状態密度 |
* 状態密度の定義$$ | * 状態密度の定義$$ | ||
行 90: | 行 90: | ||
* 状態密度データファイルの作成\\ 正常終了したプログラムの出力結果のncファイルが存在するディレクトリで次のコマンドを実行する。< | * 状態密度データファイルの作成\\ 正常終了したプログラムの出力結果のncファイルが存在するディレクトリで次のコマンドを実行する。< | ||
$ getdos input_nc_file.nc output_text_file.text | $ getdos input_nc_file.nc output_text_file.text | ||
- | </ | + | </ |
<note warning > | <note warning > | ||
出力ファイルに指定したファイル名のファイルが既にあるときは上書きされてしまう!\\ | 出力ファイルに指定したファイル名のファイルが既にあるときは上書きされてしまう!\\ | ||
詳細は[[: | 詳細は[[: | ||
</ | </ | ||
- | - 正常終了したプログラムの出力結果のncファイルが存在するディレクトリで次のコマンドを実行する。< | ||
- | $ getdos input_nc_file.nc output_text_file.text | ||
- | </ | ||
- | - CyberDuck等でoutput_text_file.textをMac側に転送する。\\ output_text_file.textの中にはエネルギーと状態密度がx-yグラフを描く形式で収められている。これは単なるテキストデータなので適当に加工することも可能。 | ||
- | - Abscissa等でoutput_text_file.textのグラフをつくる。 | ||
- | |||
seminar/al結晶.txt · 最終更新: 2022/08/23 13:34 by 127.0.0.1