- ベストアンサー
ハバードハミルトニアンについて
今、一次元の金属を考えて、同じ原子の上に電子がきた時のみ 相互作用を考えるハバード模型について勉強を始めました。 この系のハミルトニアンは H=Σ_k (E_k)(n_k)+Σ_i U(n_i,↑)(n_i,↓) １項目は電子の運動エネルギー、２項目がクーロン相互作用 エネルギーの項です。 ここの第２項のn_i,↑はi番目の原子上にある上向きスピンを持つ 電子の数（0か1）です。 ここでこのn_i,σ (σ=↑,↓)をフーリエ変換 n_i,σ=(1/NA)Σ_q (n_q,σ)exp(iqr_i)を使って書き直すと 第２項は (1/NA)Σ_q (n_q,↑)(n_-q,↓)という風に成るらしいのです。 自分自身でやってみようとすると Σ_i(1/NA)^2Σ_qΣ_q'(n_q,σ)exp(iqr_i)(n_q',σ)exp(iq'r_i) となり、フーリエ成分を代入しただけで詰まってしまいます。 (n_q,↑)(n_-q,↓)の２成分以外は消えてしまっているのですが、 ただ代入するというはじめから間違っているのでしょうか・・・ どうかご教授願います。
- 物理学
- 回答数4
- ありがとう数3
- みんなの回答 （4）
- 専門家の回答
質問者が選んだベストアンサー
その他の回答 (3)
- eatern27
- ベストアンサー率55% (635/1135)
各原子が等間隔に並んでいる場合には、 Σ_j exp(ikr_j) って、単なる等比級数ですからね。 >状況が異なるのかもしれませんが、k=πのときの項と >相殺してしまいそうな感じがしてうまくもっていけずにいます・ k=πの時の項と何が相殺すると思っていて、そう思う理由はなんでしょうか。 #n_qがqにどう依存しているか分からない以上、最初にqに関する和を取る方法は上手く計算できないような気がするのですが。。。
お礼
回答ありがとうございました。 和に関して勘違いしておりました。実際にはiに関しての 和のところを、なぜかずっとkに関する和と思いながら考えてました・・・。
- leo-ultra
- ベストアンサー率45% (228/501)
> q+q'の和を実行した時に、q+q'=0を対称にして正負の項 > でキャンセルされるということでしょうか？ 納得がいかないのであれば、これはいろいろな教科書に載っていると思いますよ。本で探したらいかがですか？ > 結果の(n_q,↑)(n_-q,↓)という項で、qの時に↑、-qの時に↓が > 残る理由が分かりません。（パウリの排他則か何かでしょうか？？？） ｑの時に↓、－ｑの時に↑でもいいんです。ｑについて和をとるわけですから。そもそもハバードが片方の演算子は↑、もう片方は↓に なっているでしょう。
お礼
回答ありがとうございました。 いろいろと勘違いしながら考えていました・・・。
- leo-ultra
- ベストアンサー率45% (228/501)
exp(iqri)とexp(iq'ri)を一緒にして下さい。 exp(i(q+q')ri) これをΣq'する。 ヒント：デルタ関数
補足
q+q'=0の時のみ値が残ってくる理由を考えてみたのですが、 q+q'の範囲は第一ブリュアンゾーン（第２？）で q+q'の和を実行した時に、q+q'=0を対称にして正負の項 でキャンセルされるということでしょうか？ またNo.1の方への補足にも書いたのですが、 結果の(n_q,↑)(n_-q,↓)という項で、qの時に↑、-qの時に↓が 残る理由が分かりません。（パウリの排他則か何かでしょうか？？？） どうかご指摘、ご教授お願い致します。
関連するQ&A
- 準自由電子近似の途中計算の質問
ゼミでテキスト発表があるのですが、シュレーディンガー方程式の途中計算で躓いています。 問題の箇所は シュレーディンガー方程式に３次元のフーリエ級数展開をした波動関数Ψ(r)とポテンシャルV(r)を代入し、この式でexp(k-(2π/a)n・r)の係数について恒等式を作るところです。 Ψ(r)=exp(ikr)ΣAn exp(-i(2π/a)n・r) V(r)=ΣVn exp(-i(2π/a)n・r) シュレーディンガー方程式 -h^2/2m▽^2ψ(r)+V(r)ψ(r)=Eψにこれらを代入して exp(k-(2π/a)n・r)の係数について恒等式を作ると [E-(h(エイチバー)^2/2m)k'^2]An=ΣAn'Vn-n' （ k'=k-(2π/a)n ） となるそうです。この恒等式を作る途中でシュレーディンガー方程式の第２項をどう処理すれば上記のようになるのかわかりません。どうしたらよいでしょうか？教えてください
- ベストアンサー
- 物理学
- 電子が落ち込まない証明
古典論では原子核を回る負電荷の電子は電磁波（エネルギー）を放出して、しだいに運動エネルギーが弱まって、静電相互作用によって原子核に落ち込んでしまうと考えられています。量子論では、量子化エネルギーE＝nh^2/(8mL^2)の式を利用することによって、電子が原子核に落ち込まないことを証明できるらしいです。でも、どうやったらいいのか分かりません。式変形が分かりません。教えてください。 ｎ：エネルギー準位 ｈ：プランク定数 ｍ：電子の質量 L：電子の動ける範囲
- 締切済み
- 化学
- 光電効果
光電効果の過程について質問です。 上司が、 光電効果とは原子が光子を吸収し、そのエネルギーの担い手として軌道電子を放出する相互作用で、K軌道の放出割合が多いのは原子核から近いためである。光子が直接軌道電子に当たるというご認識をよくしやすい。 と話をしていました。 大学では光子と物質の相互作用は、対象が４つで(1)軌道電子(2)原子核の作る場(3)核子(4)中間子であると習いました。 僕の認識では光電効果は(1)との相互作用の結果、光子が吸収される。というものです。 みなさんの解釈を教えて下さい。またこの場合光子のモーメンタムはどのように保存さているのでしょうか？
- ベストアンサー
- 物理学
- 一次元原子鎖のハミルトニアン
一次元原子鎖のハミルトニアン H=Σ_n (p_n)^2 / (2M) + K / 2Σ_n (q_n - q_n+1)^2 を量子化すると、H = Σ_k ħω_k (((b_k)^+)b_k +1/2),[b_k,(b_k')^+] = δ_kk' となることを計算したいのですが、教えてもらえないでしょうか。 参考となるページか何かでもいいです。お願いします。
- ベストアンサー
- 物理学
- 原子の励起
電圧をかけて高エネルギーを持たせた電子が水素原子の電子にぶつかると，水素原子の電子はエネルギーを得て，軌道が外側に遷移しますよね。そして，水素原子が励起状態になって，そのままでは不安定なので，安定な状態に戻る（より内側の軌道に戻る）過程でそのエネルギー差に応じた光子が出ますよね。 この軌道はn=1,n=2,n=3, …ですよね。それと，水素の電子はK殻にありますよね。ということは，n=1, n=2, …, n=∞の軌道はK殻の中にあるいうことですか？？ それと，いっぱい電子を持っている原子が励起状態になると，最外電子殻の電子（たち）のみが外側の軌道に遷移するのですか？？ 教えてください<(__)>
- ベストアンサー
- 物理学
- ハミルトニアンのユニタリー変換
原子核が原点に位置する，電子数 Z 個の原子と，光との相互作用を表すハミルトニアンを考えているのですが，ミニマル結合ハミルトニアンのユニタリー変換の過程が良く分かりません．教えていただけないでしょうか？ ミニマル結合ハミルトニアンを H' = (1/2m) Σ_{j} { p_{j}+ e A (r_{j}) }^{2} + (1/2) ∫σ(r) φ(r) dr + (1/2) (εE^{2} + (1/μ)B^{2}) とします．ややこしい式ですが，右辺第 1 項にだけ注目するので，説明は控えさせて下さい． このハミルトニアンに対し U = exp {- (ie/h~) Σ_{j} ∫_{0}^{1} r_{j}・A (λ r_{j}) dλ} = exp {- (ie/h~) Σ_{j} B (r_{j})} なるユニタリー演算子を用いて H = U^{-1} H' U というユニタリー変換を考えます． ただし，h~ = h/2π であり，また ∫ の積分範囲は 0 から 1 まで．また r_{j} は電子の座標です．A (r) はベクトルポテンシャル演算子です．B は，A の積分表記が複雑なので置き換えただけの演算子です．また p_{j} = p = - ih~ ∇_{j} は運動量演算子です． H' の右辺第一項に対しユニタリー変換を行うと U~{-1} { p + e A (r_{j}) } U = p - e ∇_{j} B + e A (r_{j}) となるらしいのですが，なぜこうなるのかが分かりません．具体的には，U^{-1} e A U = e A となるのは分かるのですが，U^{-1} p U = p - e ∇ B となるのが，なぜだか分かりません． どなたか，ご教授いただけないでしょうか？
- ベストアンサー
- 物理学
- 三角関数の和に関する証明
現在フーリエ解析関係の本を読んでいるのですが、 証明の途中で Σ[k=1~n]exp(-i2πk/n) = 0 iは虚数単位 という式が出てきました。 これを示したいと思います。 オイラーの公式を使って Σ[k=1~n]exp(-i2πk/n) ＝ Σ[k=1~n]cos(2πk/n) - iΣ[k=1~n]sin(2πk/n) 第二項が0であるということは容易に示せるのですが、 第一項が0であるということがどうも示せなくて困っています。 お力を貸していただけたらと思います。よろしくお願い致します。
- ベストアンサー
- 数学・算数
- 多電子原子の電子配置
多電子原子では、電子数が増えるにつれ、電子は順に1s→2s→2p→3s→3p→4s→3d→4p・・・の軌道を占めていくと聞きます(n+l則: Krechkuwsky則)。水素原子における電子状態の計算では、2s軌道と2p軌道は縮退していますが、多電子系では2sのエネルギー準位が2pより低くなるのはなぜでしょうか？ 多電子系では、電子間のクーロン相互作用のために縮退が解けているというのは予想できますが、2s軌道のエネルギー準位が2pのそれよりも低くなる(逆ではない)事を直感的にイメージする方法はあるでしょうか？
- ベストアンサー
- 化学
- 光電効果とコンプトン効果について
光電効果やコンプトン効果は直接電子と相互作用しているのですか？それとも原子と相互作用しエネルギーを電子に与えているのですか？ 以下のような説明が混在しており理解できません。 １．光子が物質中の電子（軌道電子）と衝突し、光子が持っていたエネルギーを電子が獲得して。。。（中略）電子が物質（原子の軌道）から飛び出す現象のこと。 ２．光子が原子により吸収され消滅し、そのエネルギーを受けて軌道電子が放出される。 もし、２．なら光子は原子のどこに吸収されているのでしょうか？軌道電子？原子核？ また、光電効果は内殻の電子のほうが放出されやすいとなっておりますが、なぜかわかりません。光子のエネルギーが束縛エネルギー以上であれば、K殻でもL殻でもさほど変わらないように思います。外殻にいくほど束縛エネルギーが低いのでコンプトン効果がおこりやすいから、でしょうか？ コンプトンについても直接電子と相互作用しているのでしょうか？また、なぜ外殻の電子のほうがはじき飛ばされやすいのかがよくわかりません。 あと内殻の電子が空位になった結果特性X線が発生するのはなんとなくわかるのですが、オージェ電子を放出したときは内殻の空位のところはどうやって埋められるのでしょうか？ 質問ばかりですがよろしくお願いいたします。
- 締切済み
- 物理学
- 逆行列のアルゴリズム
現在逆行列を求めるアルゴリズムを勉強してるんですが、全然わかりません。 とあるサイトさんから逆行列を求めるソースを拝借させていただきました。 実際に２行２列の値を入れて手計算でプログラムを追っていったら23行目で対角成分に1を代入したり、33行目で対角成分以外の成分に0を代入しているのはわかるんですが、わかるのはその程度で、もっと詳しくアルゴリズムを教えていただきたいです。 わかりやすいサイトなどでも結構なんで教えていただけると助かります！ よろしくおねがいします。 #include <stdio.h> 01:int main(void) 02:{ 03: 04: double exp_At[2][2]={{1,3},{2,1}}; //入力用の配列 05: double inv_exp_At[2][2]; //ここに逆行列が入る 06: double buf; //一時的なデータを蓄える 07: int i,j,k; //カウンタ 08: int num=2; //配列の次数 09: //単位行列を作る 10: for(i=0;i<num;i++) 11: { 12: for(j=0;j<num;j++) 13: { 14: inv_exp_At[i][j]=(i==j)?1.0:0.0; 15: } 16: } 17: /* 掃き出し法 */ 18: for(i=0; i<num; i++) /* 逆行列をexp(-At)求める */ 19: { 20: buf=1/exp_At[i][i]; 21: for(j=0; j<num; j++) 22: { 23: exp_At[i][j] *= buf; 24: inv_exp_At[i][j] *= buf; 25: } 26: for(j=0; j<num; j++) 27: { 28: if(i != j) 29: { 30: buf=exp_At[j][i]; 31: for(k=0; k<num; k++) 32: { 33: exp_At[j][k] -= exp_At[i][k] * buf; 34: inv_exp_At[j][k] -= inv_exp_At[i][k] * buf; 35: } 36: } 37: } 38: } 39://逆行列を出力 40: 41: printf("逆行列exp(-At)=\n"); 42: 43: for(i=0;i<num;i++) 44: { 45: for(j=0;j<num;j++) 46: { 47: printf(" %f",inv_exp_At[i][j]); 48: } 49: printf("\n"); 50: } 51:}
- 締切済み
- C・C++・C#
補足
kについての和を実行した時に k=0以外はキャンセルされるということでしょうか・・・ 状況が異なるのかもしれませんが、k=πのときの項と 相殺してしまいそうな感じがしてうまくもっていけずにいます・・・ また結果の(n_q,↑)(n_-q,↓)で、qの時には↑、-qの時には↓が残る 理由が分かりません・・・（パウリの排他則か何かでしょうか？？？） どうかご指摘、ご教授お願い致します。