波の表し方

波の表し方

ひよこ さんの書込 (2008/09/12(Fri) 20:52)

一辺が -l/2<x,y,z<l/2 で囲まれる立方体(三次元)の中で,ポテンシャルが \frac{l}{2}<|x|,|y|,|z| のとき V=\infty だとすると,シュレディンガー方程式の解として,複素数( e^{ir\cdot k} の形)を選ぶのは不適なのでしょうか? 周期的条件 u(x,y,z)=u(x+l,y+l,z+l) の下での解ならば複素数を選ぶのが都合がいいと教科書に書いてあるのですが,複素数を使用するのは単純に計算が簡単になるからなのでしょうか? 計算の簡略化ならば最初の問題でも複素数形の解を求められるのでしょうか? 電磁波などで複素数を使う理由がまだうまく理解できてません.どなたか教えていただけませんか?

Re: 波の表し方

toorisugari no Hiro さんのレス (2008/09/13(Sat) 11:45)

少し混乱されているようですね.

質問は (1) 無限に深い立方体の井戸型ポテンシャルにおいてシュレディンガー方程式の解として,複素数を選ぶのは何故不適切なのか?

  1. 一般に偏微分方程式を解くのに何故複素数を使うのか?

の二つですね.

(1)は後回しにして,(2)を答えます. 境界が無視できる(無限に広いor周期境界条件)場合,偏微分方程式を解くのに複素数を使うのは,解法が簡単になるからです.

\mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}} は微分演算子の固有関数

\nabla \mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}}=\mathrm{i}\bm{k} \mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}}

になるので,偏微分方程式,たとえば,

\frac{\partial^2 \phi}{\partial t^2} &= \nabla^2 \phi\\\mathrm{i} \frac{\partial \phi}{\partial t} &= -\nabla^2 \phi\\\frac{\partial \phi}{\partial t} + c \frac{\partial \phi}{\partial x} &= 0

にモード展開

\phi(\bm{r},t) = \sum_{\bm{k}} \hat \phi_{\bm{k}}(t) \mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}}

を代入すれば,(対角化された)常微分方程式系

\frac{\mathrm{d}^2 \hat \phi_{\bm{k}}}{\mathrm{d} t^2} &= -k^2 \hat \phi_{\bm{k}}\\\mathrm{i} \frac{\mathrm{d} \hat \phi_{\bm{k}}}{\mathrm{d} t} &= k^2 \hat \phi_{\bm{k}}\\\frac{\mathrm{d} \hat \phi_{\bm{k}}}{\mathrm{d} t} &= -\mathrm{i}ck_x \hat \phi_{\bm{k}}

になってしまい,簡単に解けるからです.

Re: 波の表し方

ひよこ さんのレス (2008/09/13(Sat) 23:39)

返信ありがとうございます.お時間とらせてすいません. > \mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}} は微分演算子の固有ベクトル

のところは \nabla \mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}}=\bm{k} \mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}} ということでしょうか? 申し訳ないのですが 数学については初心者レベルですので,モード展開についてあまりよくわかりません. 編微分方程式の解が一般的に,級数和で表されることは昔やったことがあるので,モード展開もそのプロセスで解釈してよろしいのですか?

結論についてですが,複素数を使うというのはやはり計算が簡単になるからであって,波(一般の編微分方程式の解として波を選ぶのは適切ですか?)の性質は,オイラーの公式により線形性が保たれているので,解析をする上で物理的意味が崩れないので,一般的に複素数を使ったほうが良いということでしょうか?

Re: 波の表し方

mNeji さんのレス (2008/09/14(Sun) 08:44)

ひよこさん,初めまして.横から失礼します.

量子力学は苦手なのですが,ちょっとコメントします.

[1] 無限大ポテンシャルについて

古典的イメージでは,「運動エネルギ+ポテンシャル・エネルギ=全エネルギ」です.立方体の内部では,各エネルギが有限値で在り得ると思われます.粒子は立方体の壁まで来ると,ポテンシャル障壁に跳ね返されるだけですよね.

量子力学とて「ポテンシャル障壁に跳ね返される」のは同じでも,その表現が違うのだと考えられないでしょうか? もしもその考え方が納得出来なければ,立方体の外側のポテンシャルを, V_{\mathrm{external}} と置いて計算し,その後に V_{\mathrm{external}} \to \infty とすれば,解り易いかもしれませんね.

[2] \nabla \mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}}=\bm{k} \mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}}

これは, \nabla \mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}}=\mathrm{i}\bm{k} \mathrm{e}^{\mathrm{i}\bm{k}\cdot\bm{r}} の事ですよね.

数学的な詳しい論議は置いといたとして,今必要なのは,定常状態のシュレディンガー(波動)方程式での表現と解の関係が大切だとおもいます.

Re: 波の表し方

yama さんのレス (2008/09/14(Sun) 17:08)

私も横から失礼します. 量子力学の波動関数の場合と電磁波の場合は事情が違うと思います.波動関数は一般には複素数ですが,電磁波は実数で表されます. 従って電磁波の場合は複素数を用いないで方程式を解くこともでき,その場合は三角関数が用いられます.しかし,最終的に実数の解が得られればいいので,途中の計算には複素数を用いてもよく,そのほうが計算が簡単になることが多いようです. たとえば三角関数の場合, \sin x を微分すれば \cos x になりますが,複素数の指数関数では, e^{ix} を微分しても ie^{ix} となって関数の形が変わりません. また, \sin(x+y)=\sin x\cos y+\cos x\sin ye^{i(x+y)}=e^{ix}e^{iy} を比べても複素数を用いるほうが見通しがよいと思います. つまり,電磁波の場合に複素数を用いるのは計算の便宜上のことで,最終的には実数部分をとるなどして実数の解を得るわけです. これに対して波動関数の場合は,計算の便宜上ではなく,複素数の解になるのが一般的です.しかし境界条件によっては実数になる場合もあります.

ご質問の場合について考えてみましょう.この場合のシュレディンガー方程式は,定数を適当に置き換えると \nabla^2u=-k^2u の形になります. u(x,y,z)=X(x)Y(y)Z(z) と置いて変数分離すると,次の方程式が得られます.

X''=-k_1^2X,\quad Y''=-k_2^2Y,\quad Z''=-k_3^2Z,\quad

ただし k_1^2+k_2^2 +k_3^2=k^2 です. この解として X=e^{ik_1x},\quad Y=e^{ik_2y},\quad Z=e^{ik_3z} が考えられます.周期的境界条件の場合は,x,y,zのそれぞれの方向の周期に合うように (k_1,k_2,k_3) の値を定めることができ,その結果

u=XYZ=e^{ik_1x}e^{ik_2y}e^{ik_3z}=e^{i(k_1x+k_2y+k_3z)}=e^{i\bm{kr}}

が境界条件を満たす解(の1つ)になります.

しかしポテンシャルの井戸の場合は,この解は境界条件を満たすことはできません. 波動関数は壁のところで0にならねばならないのに, e^{i\bm{kr}} は決して0にならないからです. この場合,いろいろな \bm k の値について, e^{i\bm{kr}} を重ね合わせて境界条件を満たすようにすることも原理的にはできます.しかし,次のように考えるほうが見通しがよいでしょう. 境界条件を満たす解を見つけるために X=e^{ik_1x} に注目すると, X=e^{-ik_1x} も解であることが分かるので,これらの線形結合 \frac{1}{2}(e^{ik_1x}+e^{-ik_1x})=\cos k_1x をつくると,これは0になりうるので k_1 の値を適当に定めれば |x|=l/2 における境界条件を満たすことができます. Y, Z についても同様に考えて境界条件を満たすように k_2,k_3 を定めることができます.その結果,立方体の内部では

u=XYZ=\cos k_1x\cos k_2y\cos k_3z

が境界条件を満たす解(の1つ)になります.もちろん,この場合は初めから X=\cos k_1x のような実数解を考えることもできるわけです.

Re: 波の表し方

ひよこ さんのレス (2008/09/14(Sun) 18:48)

>>mNejiさん 返信ありがとうございます.

>量子力学とて「ポテンシャル障壁に跳ね返される」のは同じでも,その表現が違うのだと考えられないでしょうか?

なるほど,表現の仕方ということですね.少し,数学的に考えすぎたかもしれません.目が覚めました.

>>yamaさん 返信ありがとうございます.

>これに対して波動関数の場合は,計算の便宜上ではなく,複素数の解になるのが一般的です.しかし境界条件によっては実数になる場合もあります.

これは,非常に興味深いと思いました. 複素数を電磁波のときのように使うのは,意味が違うわけですね.初めて知りました.

ポテンシャル井戸の場合でも,解は複素数の線形結合として現れ,結果的に実数になるということですか. いや,また数学的になりすぎました.ようするに古典論と量子論の表現の違いといったところですか.

なんとなく理解できたような気がします.みなさんありがとうございました.

Re: 波の表し方

yama さんのレス (2008/09/14(Sun) 20:57)

電場や磁場は測定できる物理量なので実数で表されます. しかし,波動関数は測定できる物理量ではないので,実数である必要はありません.

Re: 波の表し方

mNeji さんのレス (2008/09/14(Sun) 22:26)

>ポテンシャル井戸の場合でも,解は複素数の線形結合として現れ,結果的に実数になるということですか. >いや,また数学的になりすぎました.ようするに古典論と量子論の表現の違いといったところですか.

何か抵抗を感じます.この場合,ポテンシャル障壁が無限大でないけれど,井戸の内部での運動エネルギよりは十分に大きい場合を考えますと,井戸の外側での運動エネルギ相当の項が負になり,言い換えると波数の自乗が負になるので,本来は波導関数が \mathrm{e}^{\pm \mathrm{i}\vec k\cdot \vec r} と進行波の重畳で合ったのが, \mathrm{e}^{\pm \vec \kappa\cdot \vec r} のような単調減少関数になる為だとおもいます.敢えていえば,逆トンネル現象とでも言えるかも知れませんね.

でも,類似現象として,金属に入射した電磁場が金属中に定在波で侵入出来ずに,ジュール熱をだしながら内向きに減衰進行する過程があるとおもいます.電磁場は波動関数と違い,本質的に実関数です.それを除けば,真空中での振る舞いは進行波の重ね合わせで洗わせられているのも,この場合の井戸内部(ポテンシャルV=0)の解が同様に進行波の重ねあわせであるのとも同じです.

すこし極論かも知れませんが,量子力学と電磁気学の違いを確実に認識しておく限りにおいて,むしろ両者の物理的な類似性にも注目したいですね.