「ポテンシャルV(x、y)は 

{0<=x、y<=L}のとき0
それ以外の領域は∞

のときのエネルギー固有値と波動関数を求めよ」

という問題なんですがよくわかりません。
周期的境界条件ってこの場合ありますか?
流れだけでもいいですので教えてください。

このQ&Aに関連する最新のQ&A

A 回答 (5件)

結論から言うと、X(L)=0 、Y(0)=0です。

(この境界条件を課さないと解けません)
sotobayasiさんがつまづかれたのは式解釈の単純な誤りによります。

ポテンシャルを与えている式ですが、題意の式は
 0≦x≦L かつ 0≦y≦L
 の領域で0、
 それ以外の領域で無限大
との意味に解釈すべきです。(でないと、2次元量子井戸にならないですよね)
もし出題者がポテンシャルをsotobayasiさんの解釈のように
 0≦xで0 また y≦Lで0
と与えたかったのならばLなんて値を持ち出さずに、ポテンシャルを
 0≦xで0 また y≦0で0
と与えたはずです。(∵単なる座標の平行移動なので、本質的に同じ問題)

そもそもポテンシャル0の部分が半無限に広がっているのであれば波動関数はどこまでも広がってしまい、規格化条件∫|φ|^2 dr=1を適用できません。
またエネルギー準位も好きなものを取れますから固有値はなんでもよくなってしまいます。

以下の図はポテンシャルの井戸を上から覗いた図です。(■が無限にポテンシャルの高いところ、□がポテンシャル0のところだと思って下さい)
【正しい解釈】
y
↑→x
■■■■■■
■■■■■■
■■□□■■L
■■□□■■0
■■■■■■
■■■■■■
  0 L

【誤った解釈】
y
↑→x
■■■■■■
■■■■■■
■■□□□□L
■■□□□□
■■□□□□→無限にポテンシャル0の部分が続く
■■□□□□
  0  ↓こちらも無限にポテンシャル0の部分が続く
    • good
    • 0
この回答へのお礼

どうもありがとうございました。
おっしゃる通り勘違いでした。確かに井戸型ではないですよね。
何度もすいませんでした。

お礼日時:2001/08/05 02:17

No.2のbrogieです。



本質的なことではありまえんが、Umadaさんの回答の式でdと書くべきところを∂とウッカリ書かれてようですので、念のため書いておきます。
(2)式からは∂のところはdです。例えば、(4a)式は
-A^2(1/X)(d^2X/dx^2)=E1(4a)
になると思います。

失礼しました。
    • good
    • 1
この回答へのお礼

ご丁寧にありがとうございます。
変数分離でやればなんとか解けました。

お礼日時:2001/08/05 02:19

周期境界条件をとってクローニッヒ ‐ ペニー のモデル井戸型ポテンシャルを考えて、


井戸の深さ(壁の高さ)を∞にしても解は得られますが、
普通はそうはしないのではないでしょうか?

{0<=x、y<=L}以外の領域は∞ なので、
そこでの波動関数(?)を考えてみれば、
おのずと、x,y=0,Lのときの境界条件がわかると思います
(周期境界条件ではなく、
 剛体境界(??)条件(rigid bounray condition)で解けば良く、
 {0<=x、y<=L}の中だけ考えれば良いことがわかります。)

ちなみに、こういう四角の箱型や、円でない変な形の場合は
きれいな準位に分かれなくなって、量子カオスが起こってきます。
    • good
    • 0
この回答へのお礼

どうもありがとうございました。

お礼日時:2001/08/05 02:15

この程度の問題は量子力学のテキストには必ず記載されています。

図書館などで調べて下さい。従って、皆さんからは、回答してもらえないのではないでしょうか?

ここでは、ヒントだけ書いておきます。
2次元の問題ですから、
(1)ψ(x,y,t)=φ(x,y)T(t)
とおく、
(2)φ(x,y)=X(x)Y(t)
とおく、
X,Y,Tの微分方程式が求まります。
この微分方程式を解いて、
境界条件
(3)X(0)=X(L)=0,Y(0)=Y(L)=0
規格化条件
(4)∫│X(x)│^2dx=1など
以上の条件で解いていきますと、固有関数と固有値が求まります。

では、頑張って解いて見て下さい。分からない時は図書館行きです。自分の手で解いてみることです。
    • good
    • 0
この回答へのお礼

どうもありがとうございました。
どうやら単純に式の解釈を間違っていたみたいです(恥)。
どうもお手数かけました。

お礼日時:2001/08/05 02:12

井戸型ポテンシャルはSchroedinger方程式を解く問題としては、比較的とっつきやすいものですね。



1. 波動関数をφ(x,y)と置き、Schroedinger方程式を書き下す。
-A^2 ∇^2 φ(x,y)=E φ(x,y) (1)
ここに、A^2は (Planck定数を2πで割ったもの)^2/2m なる定数、mは粒子の質量。またこのとき定数をAでなくA^2とおいたのは後の計算上の技巧のため。

2. 変数分離型の解を求める。
すなわちφ(x,y)=X(x) Y(y)とおき、上の式に代入する。
-A^2 (Y(∂^2 X/∂x^2) +X(∂^2 Y/∂y^2)) = E X(x) Y(y) (2)
両辺をXYで除して
-A^2 ((1/X)(∂^2 X/∂x^2) +(1/Y)(∂^2 Y/∂y^2)) = E (3)
第一項はxのみの関数、第二項はyのみの関数であり、(3)が恒等的に成り立つためには
-A^2 (1/X)(∂^2 X/∂x^2) =E1 (4a)
-A^2 (1/Y)(∂^2 Y/∂y^2) =E2 (4b)
が必要十分である。ここにE1, E2は定数であり、E1+E2=Eの関係を満たす。

この方程式は簡単な線形2階微分方程式に帰着しますからすぐに解けますね((4a)の両辺にXをかけるだけ)。普通の2次元膜の振動と本質的に同じ問題です。

3. 境界条件を検討する。
井戸の外ではではVは無限大です。計算してみると分かるのですがここでは波動関数は恒等的に0になります。波動関数が境界x,y=0およびx,y=Lで滑らかにつながるという要請を考えると、
X(0)=X(L)=0 (5a)
Y(0)=Y(L)=0 (5b)
が境界条件として課されます。

4. 固有値を求める
これは(4a)(4b)式に、得られた波動関数X(x), Y(y)を代入することで求まります。
E1, E2をそれぞれ求めて、最後に足せばよいわけです。

周期的境界条件はこの場合関係ありません。(なぜなら、井戸は一個しかないわけですから)
もし井戸が多数個並んでいて、かつ井戸の障壁が有限の高さであれば周期的境界条件が関係してきます。この場合は粒子が取りうるエネルギーに幅が生じるます。(固体のバンド構造の話へとつながっていきます)

*解き方の大筋は上でよいはずですが、細かいところでタイプミス、計算ミスをしているかも知れません。検算しながら読んで頂ければ幸いです。
    • good
    • 0
この回答へのお礼

ありがとうございます。
変数分離のとこまではいけたんです。
>X(0)=X(L)=0
>Y(0)=Y(L)=0
とありますが X(L)=0 、Y(0)=0ですか?
x>=0
y<=L
でポテンシャルは0という条件からX(L)=0 、Y(0)=0と
いう境界がひけるんですか?
何度もすいませんできればお答え願います。

お礼日時:2001/08/04 13:16

このQ&Aに関連する人気のQ&A

お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!

このQ&Aを見た人が検索しているワード

このQ&Aと関連する良く見られている質問

Q物理です x^2+y^2<=1 x>=0 y>=0で与えられる重心を 求める問題で重心のx座標を

物理です
x^2+y^2<=1 x>=0 y>=0で与えられる重心を
求める問題で重心のx座標を
1/S∮(0→1)x√1-x^2となっているのですが
なぜこうなるのかがよく分かりません
解説お願いします

Aベストアンサー

重心は、任意の点の周りのモーメントを考えたときに、「微小部分の重量のモーメントの総和=全重量が重心位置にある場合のモーメント」となる点です。

 与えられたのは、半径 1 の 1/4 円の扇型です。その「微小部分」を、x座標を x ~ x+dx の「縦割り」部分にすると、面積は「高さ」が √(1 - x) 、幅が dx ですから
 ΔS = √(1 - x)*dx
です。
 この部分原点回りのモーメントの「腕の長さ」は x ですから、物理的な「力」を考えるために密度を ρ として、モーメントは
  ρ*xΔS = ρ*x√(1 - x)*dx
です。従って、「微小部分の重量のモーメントの総和」は
  ∫[0~1] ρ*x√(1 - x) dx    (1)
です。

 これに対して、「全重量が重心位置にある場合のモーメント」は、重心の x 座標を x0 とすると
  ρ*S*x0     (2)

(1)と(2)が等しくなるので
  ρ*S*x0 = ∫[0~1] ρ*x√(1 - x) dx

 従って
  x0 = (1/S)∫[0~1] x√(1 - x) dx

 S は 1/4 円なので
   S=(1/4)パイr^2 = パイ/4
ですね。

重心は、任意の点の周りのモーメントを考えたときに、「微小部分の重量のモーメントの総和=全重量が重心位置にある場合のモーメント」となる点です。

 与えられたのは、半径 1 の 1/4 円の扇型です。その「微小部分」を、x座標を x ~ x+dx の「縦割り」部分にすると、面積は「高さ」が √(1 - x) 、幅が dx ですから
 ΔS = √(1 - x)*dx
です。
 この部分原点回りのモーメントの「腕の長さ」は x ですから、物理的な「力」を考えるために密度を ρ として、モーメントは
  ρ*xΔS = ρ*x√(1 - x)*dx
です。従っ...続きを読む

Q運動方程式 m(d^2x/dt^2)+kx=0 の固有角振動数、固有周

運動方程式 m(d^2x/dt^2)+kx=0 の固有角振動数、固有周期、固有振動数ってどうやって求めるんですか?

Aベストアンサー

これは簡単な二次の線形微分方程式です。これを解くにはx=exp(λt)と置けばいいのですが、その解は

x=Asinωt+Bcosωt  ω=√k/m となります。これは加法定理を使って

x=C・sin(ωt+δ) C^2=A^2+B^2 δ=arctan(B/A) と変えることができます。

この式の周期が固有振動周期となり、その値はT=2π/ω ですよね。 固有振動数は1/T 固有角振動数はωになります。

ここで気をつけなければいけないのは、t=πのときも t=0 の値に戻って来ますが、その微分値、つまり速度は符号が逆になっているので元に戻ったことにはならず、t=2πになったとき初めて初期値に戻って来たことになることです。だからωT=2πになります。

Q周期ポテンシャルのフーリエ展開はなぜ、V(r)=Σ_{G} V(G)exp(iG.r)とあらわされるのですか。

固体物理の教科書で逆格子ベクトルのところを読んでいると、逆格子ベクトルをGとして、「周期ポテンシャルのフーリエ展開がV(r)=Σ_{G} V(G)exp(iG・r)となる」と書かれていました。なぜ、このように書けるのかが導出できません。
exp(iG・r)=cos(G・r)+isin(G・r)として、Gはゼロから無限大までそれぞれ代入したものを足し合わせると言う考え方でいいでしょうか。

Aベストアンサー

こんばんは。
まず一変数の場合のフーリエ級数は大丈夫ですか。
(指数関数表示のままの方が便利です)
たとえば周期Lの関数のフーリエ級数展開です。
それができれば、3変数に拡張するだけです。
(y、zを定数とおもってxの関数として級数展開、
つぎはzを定数とおもってyの関数を級数展開、最後zの
関数を級数展開)この時、空間の基底として格子の基本並進ベクトルをとって議論する必要があります。
そして格子に対して逆格子の定義が分れば指数関数の
肩がなぜG・rになるかわかるはずです。
以上がヒントです。

Q波動のy-x図とy-tグラフの作図

今波動を授業でやっているんですが、y-x図とy-t図の2つの図を見ても、どう見ればいいのか分からず、理解できません。
なので、y-x図とy-t図の違いを教えていただきたいです。
また、もしy-x図からy-t図へと作図をしなさいという問題が出た場合、どう書けばいいのかも教えてくだされば嬉しいです。

ぜひ分かる方教えてください。お願いします。

Aベストアンサー

まずは正弦波ですので、(現在の教育課程では習いませんが)
y=A sin(2πx/λ)
y=A sin(2πt/T)
の形で書かれています。
なので、A,λ,Tなどの値を読み取る作業をします。
それに伴い、fやvも計算していきます。

y-xグラフは 要は写真です。波の一瞬の様子を切り取っています。
y-tグラフは 波の見る位置を一点に固定しています。

y-xではtがいつの時か、y-tグラフではxがどこのものか
明記されているはずです。


y-xからy-tに変換するときは、
xがどこのy-tグラフを書くのかおさえて、vはあらかじめ与えられていると思いますので、
λ=vTでTを求めます。
振幅はy-xとy-tでは変わりません。
y-xグラフを少し時間を進めたグラフをイメージします。(連続写真のイメージ)
で、波がsinなのか-sinなのかを判断します。

これで書けるかな?

図が書けないのがつらいですが、
こんなところでしょうか。

QV=V0+at → X=V0t+1/2at^2 ?

タイトルの前者の単位は〔m/s〕ですよね
で、後者の単位は〔m〕ですよね

僕は、〔m/s〕を〔m〕に直したいなら〔s〕をかければいいと思ったので
t(V0+at)をしました
けれどそれだと、後者の式の"1/2"が抜けてしまいます
一体この"1/2"がどこから出てきたのかが疑問です

学校の先生に質問しても、積分がどうとやらといっていてよくわかりませんでした

v-tグラフの面積を利用して出すときは、加速度が斜めで出てくるから
三角形の公式を利用したときに"1/2"を使うということは分かりました
けれど、こうして式で考えようとすると、なぜ1/2が出てくるのかよくわかりません
単純にtをかけるだけではダメなのでしょうか

どなたか分かる方いたら解説お願いします

Aベストアンサー

時刻"0"から"t"までのt秒間に進んだ距離を考える場合、その中間時刻"t/2"の時の速度で"t"秒間進んだ、と考えましょう。
時刻"0"の時の速度でt秒とか、時刻"t"の時の速度でt秒よりも正しそうな気がしませんか?

時刻"t/2"の時の速度はV0+a(t/2)=V0+(1/2)atです。この速度でt秒なら
{V0+(1/2)at}t=V0t+(1/2)at^2
となります。

これは等加速度運動の場合だから成り立つのであり、常に成り立つわけではありません。
ただ、t秒間での移動距離=t秒間での平均速度 × t は必ず成り立ちます。


人気Q&Aランキング

おすすめ情報