3
大学数学基礎解説
文献あり

作用素で漸化式を解く

163
0

こんにちは,itouです.

導入

微分の世界の作用素と差分の世界の作用素には以下の対応がある.

数列(nの関数)anであって,多項式係数漸化式
L=k=0sfkNkLan=0
(fknの多項式,NNan=an+1なる演算子.なおN1an=an1.)
に対し,指数型母関数

y=n0ann!xn

が満たす微分方程式が対応する.

変数n変数x
ND=ddx
nオイラー作用素ϑ=xddx
nN1x(をかける作用素)

それぞれ見ていく.
Dy=n1an(n1)!xn1=n0an+1n!xn=n0Nann!xn

ϑy=n1an(n1)!xn=n1nann!xn=n0nann!xn

xy=n0ann!xn+1=n0nan1n!xn=n0nN1ann!xn

したがって,以下のような対応が分かる.
L(n,N,nN1)an=0L(ϑ,D,x)y=0

an=(α)n(β)n(γ)n
とすると,
((n+γ)N(n+α)(n+β))an=0
より,
((ϑ+γ)D(ϑ+α)(ϑ+β))y=0
yはガウスの超幾何級数である.上の式がガウスの超幾何微分方程式:
(x(1x)D2+{γ(α+β+1)x}Dαβ)y=0
に一致することは,
Dϑ=(ϑ1)D
に注意して計算するとわかる.

このようにして微分方程式と漸化式の世界をつなぐことができる.調和振動子を生成消滅演算子で解くことを漸化式の世界でも見てみよう.

調和振動子

調和振動子の解き方を軽く復習しておく.

ハミルトニアンは
H^=d2dx2+12x2
である.シュレディンガー方程式:
H^|ψ=E|ψ
を生成消滅演算子を使って解こう.

消滅演算子a^,生成演算子a^を以下のように定義する.
a^:=12(x+ddx)a^:=12(xddx)

交換関係:
[ddx,x]=ddxxxddx=1を用いると,

a^a^=H+12a^a^=H12
と計算できる.よって
H=a^a^+12
と交換関係
[a^,a^]=1

が分かる.

数演算子N^
N^:=a^a^
と定義する.
計算により数演算子がエルミート作用素であることはすぐに分かる.そのため,固有値は実数.固有値をn,固有関数を|nとする.つまり

N^|n=n|n.
N^a^|n,a^|nに作用しよう.

先の交換関係を用いると,
N^a^|n=(n1)a^|nN^a^|n=(n+1)a^|n
がわかる.これはa^|n|n1,a^|n|n+1が定数倍の違いを除いて同じという意味.

内積を自然に定めて計算することで,その定数を見つけることができて,
a^|n=n|n1a^|n=n|n+1

という訳でa^|n1だけ小さい固有値を対応させる.さて,na^|nのノルムの2乗に等しいが,これは0以上.よって「1だけ小さくなる」の繰り返しには限界があり,つまり

a^|0=0
であればよい.これから固有値nが非負整数であったこともわかる.

上の式を解いて,
|0=Aex22

一般に,
(a^)n|0=n!|n

漸化式へ

上でやったことをnの作用素にしよう.各種記号は混同して使う.

H^:=N2+12(nN1)2a^:=12(nN1+N)a^:=12(nN1N)

として上のように論理を展開すればよい.
なお,
DxxD=1

NnN1n=1
に対応する.

結局,
a^|0=0(a^)n|0=n!|n
を得る.

|0=sn

とおけば,a^|0=0
(nN1+N)sn=nsn1+sn+1=0.
これはすぐに解くことができ,
s2n=(1)n(2n1)!!s0s2n+1=(1)n(2n)!!s1
したがって,|nも得られた.

参考文献

投稿日:2024711
更新日:2024718
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

itou
itou
147
13478
数学勉強中. https://twitter.com/G7UOMb0Zd8V7LdP

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. 導入
  2. 調和振動子
  3. 漸化式へ
  4. 参考文献