ラメ関数

数学において、ラメ関数(Lame function)または楕円調和関数(ellipsoidal harmonic function )は、二階常微分方程式であるラメ方程式の解である。これはガブリエル・ラメ( Gabriel Lamé )の論文(1837年)で導入された。ラメ方程式は、楕円座標におけるラプラス方程式に適用される変数分離法に現れる。特殊な場合には、解はラメ多項式と呼ばれる多項式で表すことができる。  

ラメ方程式

ラメの式は

d2yd×2++B×y0{\displaystyle {\frac {d^{2}y}{dx^{2}}}+(A+B\wp (x))y=0,}

ここで、ABは定数であり、はワイエルシュトラスの楕円関数である。最も重要なケースは 、楕円正弦関数であり、整数nと楕円係数に対してであるときであり、この場合、解は複素平面全体で定義される有理型関数に拡張される。Bの他の値では、分岐点を持つ。 {\displaystyle \wp}B×κ2スン2×{\displaystyle B\wp (x)=-\kappa ^{2}\operatorname {sn} ^{2}x}スン{\displaystyle \operatorname {sn} }κ2nn+12{\displaystyle \kappa ^{2}=n(n+1)k^{2}}{\displaystyle k}

独立変数を に置き換えると、ラメ方程式は代数形式で次のように書き直すこともできる。 t{\displaystyle t}tスン×{\displaystyle t=\operatorname {sn} x}

d2ydt2+121te1+1te2+1te3dydt+Bt4te1te2te3y0{\displaystyle {\frac {d^{2}y}{dt^{2}}}+{\frac {1}{2}}\left({\frac {1}{t-e_{1}}}+{\frac {1}{t-e_{2}}}+{\frac {1}{t-e_{3}}}\right){\frac {dy}{dt}}-{\frac {A+Bt}{4(t-e_{1})(t-e_{2})(t-e_{3})}}y=0,}

変数変換後、これはホイン方程式の特殊なケースになります。

ラメ方程式のより一般的な形は、楕円体方程式または楕円波動方程式であり、次のように書くことができる(ここでは、上に示したようにではなく、と書くことに注意)。 Λ{\displaystyle \Lambda}{\displaystyle A}

d2yd×2+Λκ2スン2×Ω24スン4×y0{\displaystyle {\frac {d^{2}y}{dx^{2}}}+(\Lambda -\kappa ^{2}\operatorname {sn} ^{2}x-\Omega ^{2}k^{4}\operatorname {sn} ^{4}x)y=0,}

ここで、 はヤコビ楕円関数の楕円係数であり、と は定数である。 の場合、方程式は を伴うラメ方程式となる。 の場合、方程式はマシュー方程式に簡約される。{\displaystyle k}κ{\displaystyle \kappa }Ω{\displaystyle \オメガ}Ω0{\displaystyle \Omega =0}Λ{\displaystyle \Lambda =A}Ω00κ2hΛ2h2λ×z±π2{\displaystyle \Omega =0,k=0,\kappa =2h,\Lambda -2h^{2}=\lambda ,x=z\pm {\frac {\pi }{2}}}

d2ydz2+(λ2h2cos2z)y=0.{\displaystyle {\frac {d^{2}y}{dz^{2}}}+(\lambda -2h^{2}\cos 2z)y=0.}

ラメ方程式のワイエルシュトラス形式は計算には全く不向きである(Arscottも191ページで述べている)。最も適した形式は、上述のヤコビ形式である。代数形式や三角形式も扱いにくい。ラメ方程式は、量子力学において、様々な周期的および非調和ポテンシャルに対するシュレーディンガー方程式の古典解(周期インスタントン、バウンス、バブルと呼ばれる)に関する微小変動の方程式として現れる。[ 1 ] [ 2 ]

漸近展開

周期楕円波動関数、およびそれに伴うラメ関数の漸近展開は、大きな値に対してミュラーによって得られている。[ 3 ] [ 4 ] [ 5 ] 彼が得た固有値に対する漸近展開は、近似的に奇数の整数(境界条件によってより正確に決定される - 以下を参照) で、κ{\displaystyle \kappa }Λ{\displaystyle \Lambda }q{\displaystyle q}

Λ(q)=qκ123(1+k2)(q2+1)q26κ{(1+k2)2(q2+3)4k2(q2+5)}1210κ2{(1+k2)3(5q4+34q2+9)4k2(1+k2)(5q4+34q2+9)384Ω2k4(q2+1)},{\displaystyle {\begin{aligned}\Lambda (q)={}&q\kappa -{\frac {1}{2^{3}}}(1+k^{2})(q^{2}+1)-{\frac {q}{2^{6}\kappa }}\{(1+k^{2})^{2}(q^{2}+3)\\[6pt]&-4k^{2}(q^{2}+5)\}{}-{\frac {1}{2^{10}\kappa ^{2}}}{\Big \{}(1+k^{2})^{3}(5q^{4}+34q^{2}+9)\\&-4k^{2}(1+k^{2})(5q^{4}+34q^{2}+9){}-384\Omega ^{2}k^{4}(q^{2}+1){\Big \}}-\cdots ,\end{aligned}}}

(ここに示されていないもう1つの(5番目の)項はミュラーによって計算されており、最初の3つの項はインス[ 6 ]によっても得られています)。 と では、項が偶数項と奇数項が交互に現れることに注意してください(マシュー関数扁平回転楕円体波動関数長楕円回転楕円体波動関数の対応する計算と同様)。次の境界条件( は完全楕円積分によって与えられる4分の1周期) q{\displaystyle q}κ{\displaystyle \kappa }K(k){\displaystyle K(k)}

Ec(2K)=Ec(0)=0,Es(2K)=Es(0)=0,{\displaystyle \operatorname {Ec} (2K)=\operatorname {Ec} (0)=0,\;\;\operatorname {Es} (2K)=\operatorname {Es} (0)=0,}

主な意味の派生語)

(Ec)2K=(Ec)0=0,(Es)2K=(Es)0=0,{\displaystyle (\operatorname {Ec} )_{2K}^{'}=(\operatorname {Ec} )_{0}^{'}=0,\;\;(\operatorname {Es} )_{2K}^{'}=(\operatorname {Es} )_{0}^{'}=0,}

楕円波動関数をそれぞれ定義する

Ecnq0,Esnq0+1,Ecnq01,Esnq0{\displaystyle \operatorname {Ec} _{n}^{q_{0}},\operatorname {Es} _{n}^{q_{0}+1},\operatorname {Ec} _{n}^{q_{0}-1},\operatorname {Es} _{n}^{q_{0}}}

期間と1つを取得する 4K,2K,2K,4K,{\displaystyle 4K,2K,2K,4K,}q0=1,3,5,{\displaystyle q_{0}=1,3,5,\ldots }

qq0=22π(1+k1k)κ/k(8κ1k2)q0/21[(q01)/2]!××[13(q02+1)(1+k2)25κ+].{\displaystyle {\begin{aligned}&&q-q_{0}=\mp 2{\sqrt {\frac {2}{\pi }}}\left({\frac {1+k}{1-k}}\right)^{-\kappa /k}\left({\frac {8\kappa }{1-k^{2}}}\right)^{q_{0}/2}{\frac {1}{[(q_{0}-1)/2]!}}\times \qquad \qquad \qquad \\[6pt]&&\times \left[1-{\frac {3(q_{0}^{2}+1)(1+k^{2})}{2^{5}\kappa }}+\cdots \right].\qquad \qquad \quad \end{aligned}}}

ここで、上側の符号は解を、下側の符号は解を表す。最終的に展開する と、Ec{\displaystyle \operatorname {Ec} }Es{\displaystyle \operatorname {Es} }Λ(q){\displaystyle \Lambda (q)}q0,{\displaystyle q_{0},}

Λ±(q)Λ(q0)+(qq0)(Λq)q0+=Λ(q0)+(qq0)κ[1q0(1+k2)22κ126κ2{3(1+k2)2(q02+1)4k2(q02+2q0+5)}+]Λ(q0)2κ2π(1+k1k)κ/k(8κ1k2)q0/21[(q01)/2]![1125κ(1+k2)(3q02+8q0+3)+13.211κ2{3(1+k2)2(9q04+8q0378q0288q087)+128k2(2q03+9q02+10q0+15)}].{\displaystyle {\begin{aligned}\Lambda _{\pm }(q)\simeq {}&\Lambda (q_{0})+(q-q_{0})\left({\frac {\partial \Lambda }{\partial q}}\right)_{q_{0}}+\cdots \\[6pt]={}&\Lambda (q_{0})+(q-q_{0})\kappa \left[1-{\frac {q_{0}(1+k^{2})}{2^{2}\kappa }}-{\frac {1}{2^{6}\kappa ^{2}}}\{3(1+k^{2})^{2}(q_{0}^{2}+1)-4k^{2}(q_{0}^{2}+2q_{0}+5)\}+\cdots \right]\\[6pt]\simeq {}&\Lambda (q_{0})\mp 2\kappa {\sqrt {\frac {2}{\pi }}}\left({\frac {1+k}{1-k}}\right)^{-\kappa /k}\left({\frac {8\kappa }{1-k^{2}}}\right)^{q_{0}/2}{\frac {1}{[(q_{0}-1)/2]!}}{\Big [}1-{\frac {1}{2^{5}\kappa }}(1+k^{2})(3q_{0}^{2}+8q_{0}+3)\\[6pt]&{}+{\frac {1}{3.2^{11}\kappa ^{2}}}\{3(1+k^{2})^{2}(9q_{0}^{4}+8q_{0}^{3}-78q_{0}^{2}-88q_{0}-87)\\[6pt]&{}+128k^{2}(2q_{0}^{3}+9q_{0}^{2}+10q_{0}+15)\}-\cdots {\Big ]}.\end{aligned}}}

マシュー方程式の限界(ラメ方程式を簡約できる)では、これらの表現はマシューの場合の対応する表現に簡約されます(ミュラーによって示されるように)。

注記

  1. ^ HJW Müller-Kirsten,量子力学入門:シュレーディンガー方程式と経路積分、第2版、World Scientific、2012年、 ISBN 978-981-4397-73-5
  2. ^ Liang, Jiu-Qing; Müller-Kirsten, HJW; Tchrakian, DH (1992). 「円周上のソリトン、バウンス、スファレロン」. Physics Letters B. 282 ( 1– 2 ). Elsevier BV: 105– 110. doi : 10.1016/0370-2693(92)90486-n . ISSN 0370-2693 . 
  3. ^ W. Müller, Harald J. (1966). 「楕円波関数の漸近展開とその特性数」. Mathematische Nachrichten . 31 ( 1–2 ). Wiley: 89–101 . doi : 10.1002/mana.19660310108 . ISSN 0025-584X . 
  4. ^ミュラー、ハラルド JW (1966)。 「エルミート関数による楕円波動関数の漸近展開」。数学的表現32 ( 1-2 )。ワイリー: 49–62土井: 10.1002/mana.19660320106ISSN 0025-584X 
  5. ^ミュラー、ハラルド JW (1966)。 「楕円波動関数の漸近展開について」。数学的表現32 ( 3-4 )。ワイリー: 157–172土井: 10.1002/mana.19660320305ISSN 0025-584X 
  6. ^ Ince, EL (1940). 「VII—周期的ラメ関数のさらなる研究」.エディンバラ王立協会紀要. 60 (1). ケンブリッジ大学出版局 (CUP): 83–99 . doi : 10.1017/s0370164600020071 . ISSN 0370-1646 . 

参考文献