固体物理
そもそもの「フーリエ変換」と「フーリエ級数」の違いについては以下のサイトが詳しい。 フーリエ変換とDFTの関係 離散化と周期化ざっとまとめると、 フーリエ変換:x(連続、非周期)→ k(連続、非周期) フーリエ級数:x(連続、周期)→ k(離散、非周期)…
ASEを使ってQuantum ESPRESSOを動かすチュートリアルをやってみた。 Espresso — ASE documentation 第一原理計算高速チュートリアル · 物性実験家のための無料でできる第一原理計算入門すると、計算したフェルミエネルギーをバンド図のプロットで設定しよう…
二体演算子は、(Fermion)粒子のペアが重要であり、順列には依らない。 N粒子系に作用する二体演算子は、の全ての粒子ペアを取ったもので表せる。 の位置表示について対角的だとすると、 N粒子系のSlater行列式は、 したがって、 最初は粒子のペアだったの…
前回、繰り込み因子と自己エネルギーの関係についてまとめた。 koideforest.hatenadiary.com 自己エネルギーの虚部がゼロ()であるところでは、スペクトル関数がデルタ関数で表せるため、(準)粒子のエネルギーが上手く定義出来ている。 その(準)粒子エ…
ずっと自分の中でよく分かってなかった。 Mahanをパラパラ読んでたらちゃんと書いてあった。係数には流儀があると思うが、自由状態において係数がとなるように(振動数での逆Fourier的な積分が1になるように)スペクトル関数を定義すると、 具体例として自由…
www.jstage.jst.go.jpCurve fitting (CF)だと最小のR因子を与えるパラメータ(の組)しか最終結果が得られないが、Thorough search (TS)だと(大雑把に言って)他の極小値もわかる、というもの。 具体的には、R因子に上限(例えばRヒストグラム(ある幅を持…
単一三角格子において、各サイトへのホッピングをで表すと、ハミルトニアンは と書ける。 このハミルトニアンの固有値と固有ベクトルを求める。三次方程式の解き方を忘れたので、以下のサイトを参照。 三次方程式の解き方と例題3問 | 高校数学の美しい物語 …
あるセル内の散乱によるセルT行列を と置くと、全てのサイトによる散乱Green関数 は次のように書ける。 ここで、 は自由散乱状態におけるGreen関数である。また、セルに属するサイトをという風に表した。セルT行列 はサイトT行列 を用いて次のように書ける。…
T行列は元々サイトポテンシャルの展開で以下の様に定義されている。 しかし、「違うサイトへの伝搬」と「同じサイト内での多重散乱」が混在しており、実際の計算には非常に不便である。 そこで、同じサイトのポテンシャル散乱を繰り込んだサイトT行列を定義…
Poisson方程式の解は、 と書ける(原子単位系)。多重極子展開および角運動量展開を利用すると、 ここで、元々 の時、 が半径まで値を持つ場合において、その内側で が定数になる の形を求めることを試みる。 定数であるから、微分するとポテンシャルはゼロ…
前回、Sommerfeld展開についてまとめた。 koideforest.hatenadiary.com 今回は、(状態密度)として適用し、化学ポテンシャルの温度依存性について調べる。 ここで、 とすると、Tayler展開を用いて次のように書ける。 ここで、自由電子の状態密度が であるた…
Sommerfeld展開は、Fermi分布関数 が掛かった関数の積分値に関する展開式で、以下のように表される。 室温程度であれば、(が十分滑らかな関数である限り) 第二項までの打ち切りが良い近似として成り立つことが知られている。 この展開を導出してみることにす…
そもそもの大前提として、系が周期性を持ち、ハミルトニアンが(周期性を保つ)並進操作に対して不変であるとする。 平面波で基底を展開すると、強く局在したcore状態との直交性を記述しようとするために高い波数が必要になってしまい、実際に計算するにはと…
演算子、もしくは確率変数をと置き、それらの期待値(平均値)をそれぞれすると、 という式が成り立つ。 これはまだ厳密である。ここで、「は平均値にそれぞれ近いから、「揺らぎ」である, は凄く小さい」と考えるのが平均場近似である。 ポイントは、という…
ハバードモデルにあるオンサイトクーロン項の波数表示 を証明する。証明には、以下の関係式が必要である。 これらは以下の記事で扱ったことがある。 koideforest.hatenadiary.com koideforest.hatenadiary.comよって、
Hohenberg-Kohnの第一定理は以下の様なものである。 「外場と基底状態電子密度は一対一対応する。」 これは、背理法で証明されることがほとんどだろう。しかし、背理法そのものについて学ぶことは、あまりない気がしたので、ここでまとめる。「命題」を、「…
学生の時に教授から言われた「DFTでは強配位子場しか計算出来ない」というフレーズをふと思い出した。 その時は何を言われているのかよくわからなかったが、今は「何も工夫しなければまぁそうだろう」と思う。いくつか鍵となる概念がある。 DFTは、Kohn-Sham…
とあるスクールで、講師が事前に用意したOriginのマクロを使ってXMCDのSum ruleから磁気モーメントを求めるという講習があったが、Originを使いたいと全く思わなかったため、一人でpythonスクリプトを黙々と作っていた。ローレンツ関数やら誤差関数やらで適…
「EXAFS振動をフーリエ変換する」というフレーズはX線吸収分光をやっていると腐る程見かけるが、実際にはArtemisやらLarchやらのソフトで流し込んで、実際に気にするところはパラメータのフィッティングだけということがほとんどな気がする。 ここでは「そも…
四つ足で御馴染みのメタンの分子軌道が、固体におけるLiebフラットバンドが現れる構造と同じということで、テンションが上がった。以下の簡単な飛び移りだけを考えたハミルトニアン行列の固有値と固有ベクトルを求める。 対角成分(各原子上でのエネルギー)…
を原点から測った距離とし、原点から見て球対称な関数をとする。 原点から見て位置にあるサイトがあるとする。 サイトを原点に取り直した任意の位置ベクトルをと定義する。 やりたいことは、をサイトから見ると球対称ではないので、サイトから見たときの球対…
半径Rの構造クラスターをASEで作る。 例として単体アルミニウムを使う。 import math from ase.spacegroup import crystal from ase.visualize import view # parameters LC = 4.05 R = 7. NSL = math.ceil( 2. * R / LC ) # make bulk structure aluminium …
今更ながらIPythonデータサイエンスクックブックを購入した。もちろん私費である。せっかくなので一次元のシュレーディンガー方程式でも解こうかと思った。 一階微分方程式に対しては、かくあきさんのサイトで紹介されているローレンツアトラクターのソース…
xy平面内でのみ回転している、大きさ1の回転を考える。 つまり、 ここから元のベクトルを復元すると、ただし、微分でしか定義されていないので、一意には決まらずに任意性が残る。 例えば以下のようなものも可能である。今回は上のものを考える。これをgnupl…
とが一緒に入っている積分: ガンマ関数に持って行けないか、疑ってみましょう。ちなみにガンマ関数の定義は、 の指数がだったかだったかややこしいし、の肩がプラスだったかマイナスだったか忘れるけど、とにかくとの組み合わせ => と思っておけば見通しは…