忍者ブログ
20080511~ 13と7と11の倍数の論理積は13と7と11の積の倍数である。 和ァ・・・
[2170] [2169] [2168] [2167] [2166] [2165] [2164] [2163] [2162] [2161] [2160]
前の日記の続きです。
井戸型貞子たん
井戸型ポテンシャルでシュレディンガー方程式を解いた波動関数についてなんですが
深さが無限ではなく有限だったらどうなるのかっちゅうのを
シミュレーションで解いてみてました。


波動関数は数値計算で求めたんですが
固有値問題に関しては解析的にも計算してありまして
2012/12/1114日を参照しますと

[k1*cos(k2*a)-k2*sin(k2*a)]*[k1*sin(k2*a)+k2*cos(k2*a)]=0

を満足する波数k1とk2が固有状態というわけです。
k1とk2は、それぞれ(ħk1)2=2m(V0-E)、(ħk2)2=2mEを満足するものとします。

多項式の積がゼロになればいいという方程式になっているので
どちらかの多項式がゼロを満たせば固有状態を実現します。

この2つの多項式は実はそれぞれ、偶関数の波動関数と奇関数の波動関数が満たすべき固有条件に対応しています。
井戸の中の粒子の持つエネルギーを横軸に取った、上記の固有状態判別の式のグラフは以下のようになるのですが
固有値判別式

分けて考えると、次の2本のうちどちらかがx軸を横切ればおkということになります。
固有値判別式


解析計算と数値計算での固有値
二分法でのシミュレーションによる数値計算でも固有条件を求めてみましたが
照らし合わせてみるとだいたいあってることがわかります。

(ディラック定数ħと粒子の質量mを1、井戸の幅2a=2、ポテンシャルの深さV0=100、距離の刻み幅d=0.01として計算しました)
(精度が悪いのは刻み幅だけの問題じゃないかもしれません。循環参照を利用しているので計算順序の考慮など、その辺の甘さもあるような気がします)

シミュレーションで固有状態を求める際には
井戸の端から十分距離をとった場所で波動関数が収束する条件を利用します。

以下のように、偶関数・奇関数が交互に現れているのがわかると思います。
偶関数と奇関数の波動関数の固有値問題

(波動関数にゲタを履かせて、粒子の持っているエネルギーを表現しました。)

粒子の持つエネルギーがポテンシャルを越えると自由粒子のように振舞うこともわかるかと思います。

ブログランキング・にほんブログ村へ
素敵飲料

拍手[0回]

PR

コメント


コメントフォーム
お名前
タイトル
文字色
メールアドレス
URL
コメント
パスワード
  Vodafone絵文字 i-mode絵文字 Ezweb絵文字


忍者ブログ [PR]
カレンダー
10 2024/11 12
S M T W T F S
1 2
3 4 5 6 7 8 9
10 11 12 13 14 15 16
17 18 19 20 21 22 23
24 25 26 27 28 29 30
ブログランキング
ブログランキング参戦中
にほんブログ村 アニメブログ 深夜アニメへ
にほんブログ村 漫画ブログ SF・ファンタジー漫画へ
にほんブログ村 科学ブログ 自然科学へ
よかったらポチッとお願いします^^
最新CM
[12/30 buy steroids credit card]
[09/26 Rositawok]
[03/24 hydraTep]
[03/18 Thomaniveigo]
[03/17 Robertaverm]
最新TB
プロフィール
HN:
量子きのこ
年齢:
43
性別:
男性
誕生日:
1981/04/04
職業:
WinDOS.N臣T
趣味:
妄想・計算・測定・アニメ
自己紹介:
日記タイトルの頭についてるアルファベットは日記の番号です
26進数を右から読みます
例:H→7番目、XP→15(P)×26+23(X)=413番目。
A=0とする仕様につき一番右の桁はAにできませんのでご了承くださいズコー
バーコード
ブログ内検索
アクセス解析