MORSEの論文を読み、ポテンシャル
V[x_] := d Exp[-2 a (x - x0)] - 2 d Exp[-a (x - x0)]上の波動関数は、
R[x_, n_] := Exp(-x/2)*x^((k - 2 n - 1)/2)
*LaguerreL[k - n - 1, k - 2 n - 1, x]
と与えられることがわかったのですが、
これの直交性を確かめるためにMathematicaで
int[R[x,n]R[x,m],{x,0,∞}]
を行ったのですが、これがなぜか0になりません。
なんでですかね?
ラゲール多項式の直交性は
int[Exp[-x]*x^k*LaguerreL[n,k,x]
*LaguerreL[m,k,x],{x,0,∞}]
が0になることは確認できているので、Mathematicaの計算過程に問題はないと思うのですが…
なので、元の波動関数に問題があるのかも知れなく思えてきました。
原因がわかる方教えてください。
No.1ベストアンサー
- 回答日時:
■R[x,n]の定義式にある間違いその1
×:Exp(-x/2)
○:Exp[-x/2]
ケアレスミスですね。
■R[x,n]の定義式にある間違いその2
×:LaguerreL[k - n - 1, k - 2 n - 1, x]
○:LaguerreL[n, k - 2 n - 1, x]
モースの論文にある generalized Laguerre polynomials の引数の書き方は、MathematicaのLaguerreLのそれとは違います。MathematicaのLaguerreLの定義は
http://en.wikipedia.org/wiki/Laguerre_polynomial …
にあるものと同じはずですので、確認して下さい。
■直交性を確かめるための積分にある間違い
×:R[x,n]R[x,m]
○:R[x,n]R[x,m]/x
ポテンシャルをV[r_] := d Exp[-2 a (r - r0)] - 2 d Exp[-a (r - r0)]と定義すれば、rと波動関数R[x,n]の独立変数xと間の関係は
x ∝ Exp[-a (r - r0)]
のようになります。これを微分するとdx/dr=-a*xになりますから、
∫R(x(r),n)R(x(r),m)dr =∫R(x,n)R(x,m)/(-a*x) dx ∝∫R(x,n)R(x,m)/x dx
より、直交性を確かめるために積分される関数は、R[x,n]R[x,m] ではなくて R[x,n]R[x,m]/x になります。
rの定義域が(-r0,+∞)ではなくて(-∞,+∞)になっている理由については
http://en.wikipedia.org/wiki/Talk:Morse_potentia …
をご覧下さい。-∞<r<+∞なので、0<x<+∞になります。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 量子力学 球面調和関数 導出 方位角成分 微分方程式の解 2 2022/07/02 13:40
- 工学 制御工学の問題について 1 2022/11/01 23:45
- 工学 制御工学の問題について 1 2022/11/01 09:12
- 物理学 この波動関数の複素共役はなんですか? 2 2022/08/17 00:32
- 物理学 ポテンシャルが有限で不連続の時、右側の波動関数をφ1(x)、左側をφ2(x)とする。境界条件の「波動 2 2023/06/04 13:53
- 数学 離散フーリエ逆変換が周波数分割数をNにできる理由について 4 2022/09/18 12:56
- 数学 参考文献の探し方(数学) 1 2022/07/19 01:09
- 物理学 原子核のまわりの電子の運動 が等速円運動であると仮定する。 ただし、電子は、この円軌道 に沿って振動 1 2022/05/01 12:26
- 物理学 「次式で与えられる1次元の波動関数ψ(x,t)が自由電子のシュレディンガー方程式を満たすことを確かめ 2 2023/03/08 12:33
- 物理学 同軸ケーブル 伝送の仕組み TEMモード Maxwell方程式 円柱座標 ポアソン方程式 3 2022/08/16 20:40
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
至急お願いします!数学2の三角...
-
[溶液散乱] 回転半径Rgの計算(...
-
アボガドロ数の計算
-
高校数2です 0≦θく2πのとき次の...
-
ヤング率が変わる原因
-
許容荷重
-
なぜ、鋼材の許容引張応力は許...
-
τ ←この記号について教えてく...
-
菊ワッシャ(菊座金)の取り付け向き
-
材料力学の問題について
-
GT/sをGB/sへの換算方法
-
ワッシャー(座金)の由来
-
1kcalは何kW?という質問
-
ピースに掛かる応力計算
-
【エクセル・カメラ機能】貼り...
-
tex 図のCaptionを中央揃えにし...
-
長方形板のたわみ
-
材料力学の、『棒の両端が壁で...
-
フランジ部の増し締めとは何で...
-
NBR(ニトリルゴム)について
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報