用張量級數(shù)表示特殊函數(shù)
2017-05-12 14:26:48征夏明
中國科技縱橫 2017年5期
征夏明
摘 要：在數(shù)學(xué)和物理的研究過程中，出現(xiàn)了許多特殊函數(shù)，在各領(lǐng)域有其重要作用。一部分特殊函數(shù)是由多項式級數(shù)表示的，本文中，將運用一種新的思路，用對稱張量的級數(shù)表示出球諧函數(shù)和埃爾米特多項式，同時給出了定義高維球諧函數(shù)的方法，這將對理解其數(shù)學(xué)本質(zhì)上有所幫助。
關(guān)鍵詞：對稱張量；球函數(shù)；埃爾米特多項式
中圖分類號：O411.1 文獻標識碼：A 文章編號：1671-2064（2017）05-0252-05
Abstract：Many special functions exist in mathematics， physics or other applications. Some of them are expressed by series of polynomial. In this paper， with a new method， spherical harmonic and Hermite polynomial will be represented.
Key words：Symmetric Tensors； Spherical Harmonic； Hermite Polynomial.
1 球函數(shù)
我們知道球函數(shù)是在球坐標系下用分離變量法解拉普拉斯方程得出的。球坐標系下，拉普拉斯方程為：
（1.1）
其中關(guān)于角變量的部分可由角動量算符的平方表示：
（1.2）
這里的略去了-，只是為了便于表示，并不影響結(jié)果。
將分離為：，帶入（1.1）整理得：
（1.3）
我們想求出的特征向量，這正是所需要的球函數(shù). 解特征方程 （1.4）
同時有 （1.5）
在球坐標系中，一點的位置可由r=r表示。其中=sinθcosφ+sinθsinφ+cosθ，，，分別為x，y，z軸正方向的單位向量，用分量可寫成
。
從而我們可以將關(guān)于θ，φ的函數(shù)F看作關(guān)于的函數(shù)。
筆者在這里不加證明的做出猜想：任何關(guān)于θ和φ的函數(shù)可以展開成關(guān)于的冪級數(shù)：
（1.6）
式中是三位空間的l階張量，式中第l項是l個向量的l重點積。
事實上這是比較直觀的展開，因為任一l階張量與l個單位張量做點積，其結(jié)果確實是一個l次多項式，結(jié)合關(guān)于函數(shù)用多項式展開的知識，這二者應(yīng)當(dāng)是等價的。
思考一下，若想表示球函數(shù)，從幾何上說必須滿足空間中點的中心調(diào)和性和與長度無關(guān)的性質(zhì)，即其被定義在球面上，并且與空間中一點到原點的長度無關(guān)。
這樣我們就可以對（1.6）中的張量做出限制：
1.是對稱張量.
，是關(guān)于的任意排列（交換群）。
2.的跡為零.
，
即任取兩指標，將其設(shè)為相同的指標，再對這個指標求和（就是對任意二指標收縮）所得的l-2階張量是零張量.由于我們已經(jīng)將張量限制為對稱的，因此張量只要對某兩個指標滿足條件，則對任意兩個指標收縮求得的跡都等于0。
（注：，，…是l個不同的指標，與i，j，k等的意義相同，可在1，2，3中取值，下文同理，表示對指標收縮求跡，.為一l-2階張量。）
解，定義（）= （1.7）
對于一個給定的（），可以構(gòu)造一個新函數(shù)（r），
（1.8）
這里的和（1.7）中的是相同的，只是將其與而不是相乘后再對各指標求和，其中，即.
聯(lián)想到前面的（1.1）式，如果將（r）看作其中的j，那么（r）是否是（1.1）的解（r）=0？答案是肯定的. 那么在（1.3）中，R（r）=，F(xiàn)（θ，φ）=（），
可得，
（1.9）
這樣我們求得了的特征值為-l（l+1），對應(yīng)的特征向量即（），這正是球諧函數(shù).從而我們推出了球函數(shù)與對稱且跡等于零的張量間的內(nèi)在關(guān)系.下面證明（r）是（1.1）的解。
（r）=0 從l=0的情況開始，；
這是顯然的，因為是零階張量，為一常數(shù).同樣，當(dāng)l=1，（r）=0；
一般情形時，，是對第m個指標的自變量求偏導(dǎo)，即。
（r）=//
=（+…）+（…）+…）
=
//運用了δ函數(shù)交換指標的性質(zhì)
=（+...）+（（）+（…）+…）
由于我們限定張量是對稱且即等于0的，對任意二指標收縮等于0，從而最后一個等式中第一個括號是等于0的，而后面的部分都是與第一個括號對稱的，只是調(diào)換了指標次序.從而可以斷定（r）=0，QED.
這樣（1.9）的結(jié)論便得以印證。
下面我們思考一下以這種方式定義的球函數(shù)的代數(shù)結(jié)構(gòu).在量子力學(xué)中，我們利用對易關(guān)系和構(gòu)造了的特征向量的階梯結(jié)構(gòu)，對于每一個l值，存在m=-l，-l+1，…l-1，l共2l+1個和共有的線性無關(guān)的特征向量（因為一個分量上的角動量不可能大于系統(tǒng)的總角動量）。
我們也可以從另一個角度理解這個概念，由（1.9）已知l階對稱且跡為零的張量與的點積是和的特征向量，那么對于一個給定的非負整數(shù)l，最多有多少個線性無關(guān)的對稱且跡為零的張量？（換句話說，是所有l(wèi)階線性無關(guān)且跡為零的張量構(gòu)成的線性空間的維數(shù)是多少？）
先考慮l階線性無關(guān)的對稱張量的個數(shù)，記作。事實上這和求構(gòu)造一個l階對稱張量需要的最多獨立元素的個數(shù)是等價的，在下面的說明里可以看出。
從l=0開始，顯然一個元素就能確定，于是。
當(dāng)l=1，本身就是對稱的，其三個分量需要3個不同的數(shù)確定，則.
當(dāng)l=2，對稱矩陣可表示為：
=
需要6個分量確定.顯然所有的三維對稱矩陣構(gòu)成了一個六維的線性空間。
考慮一般情形，因為我們的張量是處于3維空間中的方形張量，每個指標可以從1取至3，這樣可以抽象成一個問題，l個指標分別如何取值的問題，也就是將l個不同的球放進標上1，2，3的箱子里有多少種放法的問題.可以用隔板法解決：將l個球和兩塊隔板放在一起混合后不考慮順序地將l個球取出（對稱），有：==其中是組合數(shù)，這樣可以得出
.
再來分析張量滿足跡等于0的條件.對稱張量滿足
方程等號左端是l-2階張量，容易驗證其仍是對稱的，因為對任意二指標收縮都等于0，這會得到多個方程.其左端都是l-2階張量，則方程的個數(shù)是.最后可得線性無關(guān)的對稱且跡為零的l階張量的最大個數(shù)為，與其他方法給出的結(jié)論相同。
目前為止，我們已經(jīng)說明了對稱且跡等于零的張量的一些性質(zhì)，下面將構(gòu)造出一組滿足這樣條件的2l+1個l階張量，且與球函數(shù)的分析性質(zhì)一一對應(yīng)，并證明與勒讓德多項式的相關(guān)性，從而達到用張量表示球函數(shù)的目的。
先考慮軸對稱的情況，則連帶勒讓德多項式退化為勒讓德多項式.此時（1.1）中的和j（1.3）中的F都是與φ無關(guān)的，我們?nèi)我獾剡xz軸作為對稱軸，當(dāng)然j和F是關(guān)于z軸對稱的.我們想到一種直接的構(gòu)造方式，即多個z軸的單位向量的并積的形式，z軸正方向單位向量記作，用坐標表示為（0，0，1），用分量表示為=，同時有·=cosθ。
我們看到，當(dāng)然是對稱的，下面將證明通過加減運算，將轉(zhuǎn)化為跡為0的張量。
定義映射T：，可將滿足一定條件的對稱張量映為跡等于零的對稱且跡等于零的張量.是所有3維空間下l階對稱張量組成的集合，是所有3維l階對稱且跡為零張量的集合.先看l=0和1的情形，此時，，顯然此時跡已經(jīng)等于0了，有T[1]=1和.再看l=2時，，可寫成矩陣：。
想將的跡化為0的一種形式較為簡單的方法是減去，I是三維的單位矩陣.寫成分量表示為，從而T[]=-.
當(dāng)l=3，，可以想象出，只有一個頂點上的元素是1，其余都是0。同樣可將其減去另一個跡為（0，0，1）的張量，也就是，形狀類似于爪形.所以T[]=。
我們發(fā)現(xiàn)了規(guī)律：將的跡化為零需要加減的張量都是和δ函數(shù)的乘積的形式，并將一直加減至δ函數(shù)的個數(shù)所可能達到的最大值[]（[n]為不超過n的最大整數(shù)）.從而對稱且跡等于零的張量
（1.10）
其中是待定系數(shù)。
在確定這些系數(shù)之前，我們先數(shù)一下這樣的l階張量的個數(shù)，由m個δ函數(shù)，l-2m個單位向量組合而成的l階張量的個數(shù)。當(dāng)m=0，顯然是1個；m=1時，是從l個指標中任選兩個指標，由于δ函數(shù)關(guān)于兩個指標是對稱的，因而與次序無關(guān)，可用組合數(shù)計算出有個這樣的張量（對應(yīng)系數(shù)為的括號中的）.再看m=2的情況，這時需要兩次連續(xù)選取兩個指標，由于δ函數(shù)之間的乘積也是可交換的，所以要再除以（為避免重復(fù)），得出m=2時有個.這樣通過歸納，可以得出由m個δ函數(shù)，l-2m個單位向量組合成的所有可能的l階張量的個數(shù)為：
（1.11）
我們可以通過解方程來確定各待定系數(shù)，也可從其與勒讓德多項式的等價性計算，最后可得出各系數(shù)。
這樣得到
（1.10）
注：式中最后一項的意義是在l為奇數(shù)時出現(xiàn)，偶數(shù)時則不存在.符號表示對的所有可能排列求和。
這是前幾階構(gòu)造出的對稱且跡為零的張量：
.
（1.7）式中，（）是與的點積，顯然這就是勒讓德多項式。
.
事實上，
，
考慮到每一項的l階張量的個數(shù)和其系數(shù)，可以看出：
（）= （1.12）
相差一個因數(shù)，可由羅德里格斯公式算出，易得：
（1.12）
綜上，可以完整的寫出（1.6）
F（）=F（θ）=1+++…
+ （1.6）
這就是軸對稱情形下的球函數(shù)。
下面分析更為一般的情形，即沒有軸對稱的情況。我們知道此時的球函數(shù)。
這個性質(zhì)說明了在坐標系對z軸作一旋轉(zhuǎn)變換時，球函數(shù)的變化過程為一個相位的變化。
我們可以從這個角度構(gòu)造出一般情況張量，即用關(guān)于z軸的旋轉(zhuǎn)變換的特征向量。
設(shè)想一線性變換，對z軸逆時針（正方向）旋轉(zhuǎn)α。
可用矩陣表示為，X=RX（1.13）
其特征方程為det（R-λI）=0，解得特征值λ=cos=cosα±isinα=（1.14）。
在實數(shù)范圍內(nèi)的情形正是我們之前計算的軸對稱情況，特征向量是z軸方向的單位向量，特征值為1.現(xiàn)在我們需要尋找其他滿足條件的復(fù)特征向量。
將λ=代回方程（1.13），解得=，最后得到，，加上，三者兩兩正交，我們發(fā)現(xiàn)
。
嘗試以這樣的方式構(gòu)造：或，由于， ，顯然這樣的張量是對稱且跡等于零的。
為了使結(jié)果與成正比，可令其包含m個或的并積（m>0時用，m<0時用），，同時其余的部分是一個最高次為l的關(guān)于cosθ的多項式，化簡后包含sinθ的m次冪和cosθ的l-m次冪的乘積.由于，可以嘗試用l-m個和m個的并積構(gòu)造張量，.從而可以取其對稱且跡等于零的部分，來匹配滿足上文中討論的滿足條件的張量的代數(shù)結(jié)構(gòu)。
先將（下文中記作）化為對稱的形式，再用變換得到同時滿足分析性質(zhì)和代數(shù)性質(zhì)的張量.定義，M將一張量化為對稱張量。
，就是對所有指標的任意排列求和，再除以l！.
則，再取
，得到我們需要的對稱且跡為零的張量，與（1.7）類似，這里設(shè)（1.15）。
這是前幾階的張量。
計算后對應(yīng)的前幾項為：
注：書寫過于復(fù)雜，不妨記作表示為取一張量對稱且跡為零的部分，這也包括了軸對稱的情形，因為永遠對稱。
然而，對于更大的l，會變得難以計算，可以利用如下的等式：
（1.16）
即不直接計算，間接計算 ，從而求得（1.15）中的球函數(shù).
我們知道球函數(shù)用多項式表示為：
（1.17）
是連帶勒讓德多項式。
同時，可用羅德里格斯公式推知系數(shù)
（1.18）
可得。
現(xiàn)在觀察中或的個數(shù)，若依次增大或減小的個數(shù)，則有2l+1種可能的（這里限定和不能同時存在），對應(yīng)有2l+1種可能的，聯(lián)想到量子力學(xué)中角動量算符的特征向量的階梯結(jié)構(gòu)，其中的關(guān)聯(lián)是顯然的。
當(dāng)用升降算符作用于一個球函數(shù)時，將會得到或的量子態(tài)，其結(jié)果正是將的的個數(shù)增大或減?。╩≥0時作用會將一個替換為，作用將一個替換為；m<0時相反.）再取的對稱且跡為零的部分.
最后，可以寫出以球函數(shù)為基底的張量級數(shù)展式：
（1.19）
事實上這樣的球函數(shù)張量不僅適用于三維空間，也可以推廣至更高維空間，只需取更高維數(shù)中的，和，用類似的方法計算系數(shù)，便可得出高維空間內(nèi)的“超球函數(shù)”。
2 埃爾米特多項式
在之前討論的基礎(chǔ)上，埃爾米特多項式也可以用同樣的方法以張量級數(shù)表示。
我們知道埃爾米特多項式是埃爾米特方程的解：
y-2xy+2ny=0 （2.1）
設(shè)一n階張量是對稱的，可以將上式改寫為張量方程的形式：
（2.1）
式中是Nabla（Del）算符，為，I是d×d的單位矩陣，表示張量積，ζ是方向上的長度為x的向量（x，0，0…，0），用分量表示為=x.
這個方程也可寫成分量的形式：
（2.1）
其各指標可在1，2…d中取值，意義是我們是在d維空間中討論這個問題的。
定義算符
（2.2）
式中最后一項的意義是n為奇數(shù)時ζ存在，偶數(shù)時則不存在.的形式與（1.10）相似，區(qū)別在于其目的不是將對稱張量的跡化為零，而是將其轉(zhuǎn)化為滿足埃爾米特方程的解.能看出每一項的系數(shù)正負交替且以的比例減小。
這是前幾階的埃爾米特對稱張量：
關(guān)于埃爾米特多項式的遞推關(guān)系和微分性質(zhì)同樣有張量的形式：
（2.3）
證明：由（2.2），
（偶數(shù)情形，奇數(shù)同理.），
而
，
通過計算每個括號中相同張量的個數(shù)，并利用恒等式，可得證.□
（2.3）式可得出關(guān)于相鄰的3階埃爾米特對稱張量的遞推公式：
（2.4）
證明：將（2.3）式，
寫成，同時將（2.1）變換為
，
將這三式整理，
得
，
雖然張量積的運算一般不滿足消去律，但對單位矩陣I成立，從而QED□
最后通過搭配系數(shù)，埃爾米特多項式和埃爾米特對稱張量的關(guān)系可表示為：
（2.5）
是方向的單位向量，坐標表示為（1，0…，0），分量表示為。
3 結(jié)語
利用張量級數(shù)表示特殊函數(shù)，其數(shù)學(xué)規(guī)律更容易記憶，物理意義或幾何意義更為鮮明，這是一種新的思路和方法，缺陷在于實際計算中并不易于使用。
注：文中關(guān)于張量運用了多種不同的符號，在不同的計算情形下選擇了不同的表示方法，例如張量就表示這個n階張量的整體，而是這個張量的一個分量，當(dāng)指標確定了表示哪一個分量也就確定了，比如就是第一纖維取1，第二纖維取2…，第n纖維取n的這個分量（第幾纖維相當(dāng)于矩陣中行與列的推廣）.文中出現(xiàn)的所有張量都是方形張量，即所有纖維上的指標所可能取的個數(shù)都是相等的（相當(dāng)于矩陣中方陣的）.可寫成關(guān)于空間中選定基底的線性組合的形式：，當(dāng)然可用愛因斯坦求和約定省略求和符號。
參考文獻
[1]Guth Alan.The Tressless Symmetric Tensor Explanation and standard Spherical Harmonic（.pdf） [EB/OL]. http：//ocw.mit.edu， 2012-10-17.
[2]陳玉麗.補充材料，張量分析初步（.ppt） [EB/OL]. http：// wenku.baidu.com，2013-7-21.
[3]李政道.物理學(xué)中的數(shù)學(xué)方法[M].江蘇：江蘇科學(xué)技術(shù)出版社，1980-6.
[4]西安理工大學(xué)應(yīng)用數(shù)學(xué)系.第七章：埃爾米特（Hermite）多項式（.ppt） [EB/OL]. http：// wenku.baidu.com，2014-11-8.
[5]Gasiorowicz Stephen. Quantum Physics [M].John Wiley & Sons，2003.
