以下のような1次元の井戸型ポテンシャル
              V(x)=0 (-L<x<L)
                =V  (x<-L, x>L)
      ただし 0<E<V
中の質量mの粒子について・・・・
この問題でグラフの交点を求めることによって、固有値が求められますが単純な計算では出せません。よって、以下のような課題を出されました。

上の問題で適当なVとLについて、固有値をニュートン法などの簡単なプログラムを組んで、数値的に計算して求めよ。また、そのときの固有関数を求めてプロットせよ。

以上の問題なのですが、簡単なヒントなのでよろしいので分かる方がいたら教えてください。お忙しい中ありがとうございました。    

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

A 回答 (1件)

> この問題でグラフの交点を求めることによって、


> 固有値が求められますが単純な計算では出せません。
と書かれているところを見ますと,
領域をx<-L, -L<x<L,L<x,の3つに分けて
それぞれの領域の波動関数を接続して tan あるいは cot を含んだ
連立方程式が出るのはご存知なんですよね.

そうすると,あとはそれこそ問題通りに
> 上の問題で適当なVとLについて、
> 固有値をニュートン法などの簡単なプログラムを組んで、
> 数値的に計算して求めよ
だけだと思うんですが...
連立方程式の数値解がわかれば波動関数のグラフもすぐ描けますし.
    • good
    • 0

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

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

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

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

Q一次元井戸型ポテンシャル、井戸の外でのシュレディンガー方程式は?

井戸の深さが無限の一次元井戸型ポテンシャルで、井戸の外において電子が満たすべきシュレディンガー方程式を求める問題があるのですが、井戸の外では波動関数φ(x)=0なのでシュレディンガー方程式は0になると考えたのですが、合っているでしょうか?

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量子力学 井戸型ポテンシャルでの場合分けについて

有限深さの井戸型ポテンシャルの問題はポテンシャルの有る領域と無い領域とで2つのシュレディンガー方程式を立てますよね。 これは当然なんですが、その粒子のエネルギーEとポテンシャル障壁Vとの大小関係、つまりE>Vなのか0<E<Vなのかで場合分けする必要は無いのでしょうか。階段型ポテンシャルやトンネル効果の問題ではエネルギーにおける場合分けはしているのですが、有限井戸型ポテンシャルの問題では場合分けをしてる参考書があまり無いような気がします。授業でも先生はE>V、0<E<Vとで場合分けはしていませんでした。場合分けしないのは階段型とは違って透過や反射に焦点を当てていないからですか?

どなたか宜しくお願いします。

Aベストアンサー

E<Vなら束縛問題、E>Vなら散乱問題になるので性格が全く異なります。
井戸型ポテンシャルはたいがいは束縛問題で扱う物なので、
普通はE<Vの場合しかやらないと思います。
おそらく、講義の最初にそんなことを言っていると思いますけどね。

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井戸型ポテンシャルのアナロジーについて

V(x)=∞ (x<0), V(x)=0 (0<x<L), V(x)=V0 (x>L)の許容エネルギーは深さが同じV0で、幅が2Lの有限井戸型ポテンシャルの奇関数解に対するエネルギーに等しい理由を教えてください。

Aベストアンサー

シュレーディンガー方程式の解そのものはあまり重要ではありません。
シュレーディンガー方程式そのものを見れば十分です。

以下の点は分かりますか?

1.x>0の領域だけを見るとどちらの場合もシュレーディンガー方程式は同じ。
2.x<0でV(x)=∞の方についてはx=0で境界条件(ψ(0)=0)を課す事になる
3.幅2Lの有限井戸型ポテンシャルの奇関数解は2番の境界条件およびx>0でのシュレーディンガー方程式を満たしている。

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無限に深い井戸型ポテンシャルについて

無限に深い井戸型ポテンシャルについて、GaNの層厚d=0.5nm、伝導帯電子の有効質量me=0.2m0の
とき

(1)基底状態、第一励起状態および第二励起状態のエネルギー固有値E1,E2,E3をeVの単位であらわすとどうなるのですか?

(2)GaNの伝導帯の3次元状態密度および、この問題のような2次元の状態密度を計算した場合の状態密度とエネルギーの関係はどうなるのでしょうか?

自分で勉強してみたものの、無限に深い井戸型ポテンシャルだけは良くわかりません。どなたか教えていただけるとさいわいです。

Aベストアンサー

(1)「井戸の壁の高さが無限大」から言えることは、「壁面において波動関数の値はゼロでなくてはならない」です。これが計算の過程において境界条件として課されます。
ご質問の場合の量子井戸内の波動関数は以下のようになります。左側の壁をx=0の原点に取り、また井戸の壁を太線(┃)で示しています。図示の制約から正確には描けませんが雰囲気は分かっていただけると思います。古典的な振動のイメージで言うならば「壁面の位置と振動の節が一致している」ということです。

[基底状態]
・      ・
:↑波動関数 :
┃      ┃
┃  __  ┃
┃ /  \ ┃ 
┃/    \┃←壁の位置では波動関数ψ=0
┃      ┃
x=0     x=d


[第一励起状態]
・      ・
:      :
┃      ┃
┃ _    ┃
┃/ \   ┃←壁の位置では波動関数ψ=0
┃   \ /┃
┃     ̄ ┃
     

[第二励起状態]
・      ・
:      :
┃      ┃
┃      ┃
┃/\  /\┃←これまた、壁の位置では波動関数ψ=0
┃  \/  ┃
┃      ┃


井戸内部における、定常状態のSchroedinger方程式は
(h~^2/2m_e)(d^2/dx^2)ψ+Eψ=0  (1)
です。ψは波動関数、h~はhバー(Planck定数÷2π)です。また「m×e」との混乱を避けるため有効質量はmeでなくm_eと表記しました。
この方程式は2階線形の斉次微分方程式ですから
ψ=A1 sin(kx) + A2 cos(kx)  (2)
の形の解となることは自明です。ここにkは{√(2m_e E)}/h~です。またx=0でψ=0との境界条件がありますからA2=0もすぐに分かります。
さてもう一つの境界条件、x=dでψ=0を検討します。これは(2)にx=dを代入することで
d×{√(2m_e E)}/h~=π, 2π, 3π,...,nπ  (3)
とすぐに書き換えられます。
(3)をエネルギーEの式に変形すると
E=n^2 {(h~π^2}/2m_e d^2  n=1,2,3,...  (4)
ということになります。エネルギーEはnに応じてとびとびの値しか許されないことが分かります。
題意から、n=1が基底状態、n=2が第一励起状態、n=3が第二励起状態に相当します。数字を代入して計算するのはセルフサービスでいいですよね。m_eをkg単位、h~をJ・s単位で代入したならば最終結果はJ単位で出てきますから、電荷素量で除して電子ボルトの単位に直すのをお忘れなく。

(2)これも上記と同様に考えることができます。
3次元の大きな直方体の箱を考えて下さい。各辺の長さをL_x, L_y, L_zとします。この箱の中に電子が存在するとしましょう。箱の内部ではポテンシャルは一様でゼロとし、箱の外ではポテンシャルは∞とします。従って箱の内壁面において電子の波動関数ψが0という境界条件が課されます。
(1)では電子の運動は1次元(1自由度)で考えましたが、今度は3次元で考えなくてはなりません。3次元空間のSchroedinger方程式(定常状態)は箱の内部においては
(h~^2/2m_e)(∇^2)ψ+Eψ=0  (5)
となります。∇^2はラプラシアンです。この方程式はψを変数分離(ψ=X(x)×Y(y)×Z(z))して解くのが常道です。解き方は量子力学の初歩の教科書には大抵載っていますから、詳細はそちらで読んで下さい。
計算結果だけ示しますと、やはりこの場合もエネルギーEの値はとびとびのものに限られ、その値は
E=(h~^2/2m_e)×{(πn_x/L_x)^2 +(πn_y/L_y)^2 +(πn_z/L_z)^2}  (6)
となります。n_x, n_y, n_zは整数で1,2,3,..の値を取ります。これは(4)と同じことで各壁面でψ=0という境界条件が課されるところから生じていますが、今度はx, y, zの3方向ありますからある一つの状態を表現するには整数の3つ組が必要、ということです。

ではいよいよ状態密度Z(E)、すなわち単位体積・単位エネルギー領域で収容できる電子数について考えてみます。
まず波数という物理量を導入します。一般にある波動があった場合、その波長λに対し2π/λという物理量を考え「波数」と呼びます。長さ方向にどれくらい波が詰まっているかを表す量です。
今考えている箱の中の波動関数で、例えばL_x/n_xは半波長に相当する量ですから、この場合の波数(x方向の波数)はπn_x/L_xということになります。同様にy方向についてπn_y/L_y、z方向についてπn_z/L_zがそれぞれ波数です。
波数はx, y, zの各方向についてk_x, k_y, k_zなどと表され、また(k_x, k_y, k_z)のように3つ組にしたものは「波数ベクトル」と呼ばれます。
波数を用いると(6)は
E=(h~^2/2m_e)×{k_x^2 +k_y^2 +k_z^2}  (7)
と変形できます。この式は波数とエネルギーの関係を表すものです。また波数を変数と考えると(7)は波数空間(3次元)での球面の式を表していることになります。(ただし定常状態となり得る波数がπ/L_xなどの整数倍、すなわち格子点に限られることは変わりありません)

これを知った上でエネルギーEとE+dEの間にある状態数Z(E)dEを考えてみましょう。
さてまず(7)で示したように等エネルギー面は波数空間で球面をなします。半径{√(2m_e E)}/h~の球と、半径[√{2m_e (E+dE)]/h~の球の間の体積差を考えると
[{√(2m_e)}/h~]^3 (4π/3){(E+dE)^(3/2)-E^(3/2)}  (8)
ですが、dEはEに比べて十分に小さいので
[{√(2m_e)}/h~]^3 (4π/3){(3/2)√E dE}
=2π[{√(2m_e)}/h~]^3 √E dE  (9)
となります。
一方、許される波数は波数空間で格子点として存在しています。その密度は波数空間の体積π/L_x×π/L_y×π/L_zあたりに1つですが、1つの準位に電子は2つまで入れますから電子の状態数にするには2を掛けます。
これを(9)にかけると、
2×(L_x×L_y×L_z)÷π^3×2π[{√(2m_e)}/h~]^3 √E dE
=(L_x×L_y×L_z)×4√(2m_e) √E÷(π^2×h~^3) dE  (10)
を得ます。n_x, n_y, n_zは正の値に限られますので、(10)の1/8が実際の状態の数ということになります。また状態密度は単位体積で定義していますから、(10)をさらに実空間の体積L_x×L_y×L_zで除して
(1/8)×{4√(2m_e)^3 /π^2×h~^3} E^(1/2) dE  (11)
を得ます。従って状態密度Z(E)は
Z(E)={√(2m_e)^3 /2 π^2×h~^3} E^(1/2)  (12a)
あるいは
Z(E)=4π{√(2m_e)/h}^3 E^(1/2)  (12b)
ということになり、いずれにしてもZ(E)はエネルギーEの平方根に比例することが分かります。ここにhはPlanck定数です。グラフにするとおよそ以下のようになります。

↑Z(E)
│      *
│   *
│ *
│*
└───────→E

以上は3次元でのお話です。では2次元ではどうなるでしょうか。2次元、すなわち薄膜中での電子伝導は(7)式までは同じですがその先が少し異なります。
薄膜はx方向について極めて薄く、y方向およびz方向には十分に大きく広がっているとしましょう(L_x≪L_y, L_z)。すると波数空間において、許されるk_xの値は極めて限られたものになり、その相互の間隔は許されるk_yやk_zの間隔に比べて非常に大きくなります。格子点はy方向とz方向には稠密ですが、x方向には相当に離散的ということです。
この場合エネルギーEとE+dEの間の状態数の差は、球の体積の差でなく円の面積の差で求めることになります。
まずx方向の基底状態はn_x=1について、
E=(h~^2/2m_e)×{(π/L_x)^2 +(πn_y/L_y)^2 +(πn_z/L_z)^2}
{(2m_e E/h~^2)-(π/L_x)^2}=(πn_y/L_y)^2 +(πn_z/L_z)^2  (13)
となります。(13)は球の方程式でなく、円の方程式と看做すべき式です。
波数空間の平面k_x=π/L_xにおいて、半径√[(2m_e E/h~^2)-(π/L_x)^2]の円と半径√[{2m_e (E+dE)/h~^2}-(π/L_x)^2]の円の面積差は明らかに
π{2m_e/h~^2} dE  (14)
です。
またこの平面内において許されるL_y, L_zの面積密度は、π/L_y×π/L_zあたりに2つです。2つになる理由は上記と同じ(スピンの差異で2つの電子が入りうるから)です。
さらにn_y, n_zは正の値に限られることを考えると、許される状態の数は
2×(1/4)×π{2m_e/h~^2}×L_y×L_z÷π^2 dE  (15)
ということになります。実空間の単位体積当たりに直せば
2×(1/4)×{2m_e/h~^2}/(L_x π) dE
=4π m_e/(L_x h^2)  (16)
ということになります。すなわちEの値によらず状態密度は一定という結果になります。
ただし、ある程度Eの値が大きくなってくると今度はn_x=2の場合が入ってきますから、それも含めて状態数を計算する必要が生じてきます。しかしn_x=2であっても(16)の結果は同じであり、n_x=1についての状態密度とn_x=2についての状態密度を単に足し合わせれば全体の状態密度になるということです。
これは言葉で表現するより図を見てもらった方が早いでしょう。以下のような階段状のグラフになります。なお階段の段の高さはどこでも同じです。

↑Z(E)
│          
│      *******
│   *****
│ ***
│ **
└───────────→E

以上はかなり駆け足で説明したのと計算間違いをしている可能性があるのとから、miake-kiyoshiさんご自身で検算しながら、また教科書などで復習しながら読んでいただければ幸いです。

参考URLのページには精細な波動関数の図や状態密度のグラフがありますからぜひ一読ください。

参考URL:http://www.qed.eedept.kobe-u.ac.jp/japanese/semicon/semicon_basics.htm

(1)「井戸の壁の高さが無限大」から言えることは、「壁面において波動関数の値はゼロでなくてはならない」です。これが計算の過程において境界条件として課されます。
ご質問の場合の量子井戸内の波動関数は以下のようになります。左側の壁をx=0の原点に取り、また井戸の壁を太線(┃)で示しています。図示の制約から正確には描けませんが雰囲気は分かっていただけると思います。古典的な振動のイメージで言うならば「壁面の位置と振動の節が一致している」ということです。

[基底状態]
・      ・
:↑...続きを読む

Q単振動の解x(t)=Asin(ωt+φ)において、次の①②を満たす解を求めよ。 ①初期条件x(0)=

単振動の解x(t)=Asin(ωt+φ)において、次の①②を満たす解を求めよ。
①初期条件x(0)=x0、x°(0)=0 (x(t)=dx(t)/dt)

②初期条件x(0)=0、x°(0)=v0

これの途中式と答えを教えて下さい。
お願いします。

Aベストアンサー

x° は時間微分ということですね? x°(t)≡dx/dt ということ。

だとすると、何が分からないのですか? 「微分」が分からないということ?
x(t)=Asin(ωt+φ)
ならば
dx/dt = Aωcos(ωt+φ)
ということは分かりますよね?

問題は「解を求めよ」ですが、「初期条件の値を求めよ」ということでしょうか。
そうであれば
① x(0) = Asin(φ) = x0

 x°(0) = Aωcos(φ) = 0
 → φ= [ (2n+1)/2 ]パイ (n:任意の整数)
このとき sin(φ) = 1 または -1
よって
 x0 = ±A

② x(0) = Asin(φ) = 0
 → φ=nパイ (n:任意の整数)
 x°(0) = Aωcos(φ) = v0
  φ=nパイ (n:任意の整数)なので cos(φ)=1 または -1
よって
 v0 = ±Aω

Q井戸型ポテンシャル(JJサクライ)

JJサクライ下巻のp397(5.1.15)
摂動論ではないのですが、
さらりと触れられている箇所が納得できずにモヤモヤしています。

非常に弱い一次元の井戸型ポテンシャルを考えます。

V=-V0 (|x| < a)
V=0 (|x| > a)

λ>0 の引力に対して
E = - (2ma^2)/h^2 |λV0|^2
というようなエネルギーの束縛状態がある。
(ここでhはh/2πの意味です)

このエネルギーの式はどのように導いたのでしょうか。
単純な井戸型ポテンシャルでもなさそうですし、
トンネル効果と比較してもよく分かりません。
よろしくお願いします。

Aベストアンサー

λが唐突に出てくるのですが、
|x|<aでV=-λV_0、あるいは同じ事ですが、ハミルトニアンをH=p^2/2m + λVとしているという事でいいですね?

エネルギーE(<0)の束縛状態の波動関数(のうち偶関数のもの)は、
x>aで、ψ=Aexp(-ρx)
|x|<aで、ψ=Bcos(kx)
x<-aで、ψ=Aexp(ρx)
のようになっています。ここで、E=-h^2ρ^2/2m=h^2k^2/2m -λV0です。・・・(◎)
※(λV0が大きければ)奇関数の束縛状態も存在しますが、基底状態のパリティが正なので、簡単のため奇関数の束縛状態は考えない事にします。

x=±aで、ψとその微分が連続という条件から、ρ=k tan(ka)が出てきます。・・・(☆)
※実際には、◎と☆の交点が固有状態を与えます。(もちろん、奇関数解で☆に相当する式も考える必要はあります)

λV0に関する最低次の項のみを考えると、
☆よりρ=k^2a.これと◎から、h^2k^2/2m=λV0となります。
故に、E=-h^2ρ^2/2m=-(2ma^2)/h^2 |λV0|^2となります。

λが唐突に出てくるのですが、
|x|<aでV=-λV_0、あるいは同じ事ですが、ハミルトニアンをH=p^2/2m + λVとしているという事でいいですね?

エネルギーE(<0)の束縛状態の波動関数(のうち偶関数のもの)は、
x>aで、ψ=Aexp(-ρx)
|x|<aで、ψ=Bcos(kx)
x<-aで、ψ=Aexp(ρx)
のようになっています。ここで、E=-h^2ρ^2/2m=h^2k^2/2m -λV0です。・・・(◎)
※(λV0が大きければ)奇関数の束縛状態も存在しますが、基底状態のパリティが正なので、簡単のため奇関数の束縛状態は考えない事にします。

x=±a...続きを読む

Q固有値・固有ベクトルの求め方

固有値・固有ベクトルの求め方
ある行列をA、単位行列をE、Aの固有値をλ、固有ベクトルをuとすると、
(λE-A)u=0
を立てて、(λE-A)が逆行列を持たないことから、λはわかりますよね?そこでλを(λE-A)u=0に代入してuを求めると教科書にあるんですが、0しか出てきません…
どうしたら良いのでしょうか?他に方法があるのでしょうか?

Aベストアンサー

#3です。
余計なことを書いたみたいだから、
添付のように訂正します。
スカラーu2倍はすべて固有ベクトルです。


人気Q&Aランキング

おすすめ情報