バーンズのG関数

数学において、バーンズの G-関数（バーンズのGかんすう、英: Barnes G-function）G(z) は、スーパー階乗を複素数にまで拡張した特殊関数 である。これはガンマ関数、K関数、グレイシャーの定数に関連するものであり、数学者であるエルンスト・ウィリアム・バーンズ（英語版）にちなみ名付けられた。 これは（初等函数を掛ける違いを除いて）二重ガンマ関数の特殊な場合である。 正式には、バーンズの G-関数は以下のワイエルシュトラスの乗積表示 G ( 1 + z ) = ( 2 π ) z / 2 exp ( − z + z 2 ( 1 + γ ) 2 ) ∏ k = 1 ∞ { ( 1 + z k ) k exp ( z 2 2 k − z ) } {\displaystyle G(1+z)=(2\pi )^{z/2}{\text{exp}}\left(-{\frac {z+z^{2}(1+\gamma )}{2}}\right)\,\prod _{k=1}^{\infty }\left\{\left(1+{\frac {z}{k}}\right)^{k}{\text{exp}}\left({\frac {z^{2}}{2k}}-z\right)\right\}} の形で定義される。ここで γ はオイラーの定数であり、exp(x) = ex は指数関数である。また、∏ は総乗の Π-記法である。 バーンズの G-関数は、正規化条件 G(1) = 1 のもと以下の函数等式 G ( z + 1 ) = Γ ( z ) G ( z ) {\displaystyle G(z+1)=\Gamma (z)\,G(z)} を満たす。このバーンズ函数の満たす函数等式とガンマ函数の満たす函数等式 Γ ( z + 1 ) = z Γ ( z ) {\displaystyle \Gamma (z+1)=z\,\Gamma (z)} との類似性に注目せよ。この函数等式を用いることにより、バーンズ G が整数引数に対して以下の通り G ( n ) = { 0 if  n = 0 , − 1 , − 2 , … ∏ i = 0 n − 2 i ! if  n = 1 , 2 , 3 … {\displaystyle G(n)={\begin{cases}0&{\text{if }}n=0,-1,-2,\dots \\\prod _{i=0}^{n-2}i!&{\text{if }}n=1,2,3\dots \end{cases}}} を値とすることが導かれる（特に、G(1) = G(2) = 1 が従う。またこれにより、 G ( n ) = ( Γ ( n ) ) n − 1 K ( n ) {\displaystyle G(n)={\frac {(\Gamma (n))^{n-1}}{K(n)}}} がわかる。ただし Γ(x) はガンマ関数を、K はK関数を表す。上記の函数等式は、凸条件 d3/dx3G(x) ≥ 0 を追加すれば、一意にバーンズ G-函数を定義する。 バーンズの G-関数に対する差分方程式は、ガンマ関数の函数等式と合わせて、バーンズの G-関数の反射公式（英語版）（相反公式） log ⁡ G ( 1 − z ) = log ⁡ G ( 1 + z ) − z log ⁡ 2 π + ∫ 0 z π x cot ⁡ π x d x {\displaystyle \log G(1-z)=\log G(1+z)-z\log 2\pi +\int _{0}^{z}\pi x\cot \pi x\,dx} (1) を得るのに用いることができる（もともとはヘルマン・キンケリン（英語版）によって証明された）。右辺に現れる対数正接積分は（二次の）クラウセン関数（英語版）を用いると、 2 π log ⁡ ( G ( 1 − z ) G ( 1 + z ) ) = 2 π z log ⁡ ( sin ⁡ π z π ) + Cl 2 ⁡ ( 2 π z ) {\displaystyle 2\pi \log \left({\frac {G(1-z)}{G(1+z)}}\right)=2\pi z\log \left({\frac {\sin \pi z}{\pi }}\right)+\operatorname {Cl} _{2}(2\pi z)} と評価することができる。この結果の証明は、対数余接積分 Lc(z) の以下のような評価と dlog(sin πx)⁄dx = π⋅cot πx なる事実によるものである。部分積分により L c ( z ) = ∫ 0 z π x cot ⁡ π x d x = z log ⁡ ( sin ⁡ π z ) − ∫ 0 z log ⁡ ( sin ⁡ π x ) d x = z log ⁡ ( sin ⁡ π z ) − ∫ 0 z [ log ⁡ ( 2 sin ⁡ π x ) − log ⁡ 2 ] d x = z log ⁡ ( 2 sin ⁡ π z ) − ∫ 0 z log ⁡ ( 2 sin ⁡ π x ) d x {\displaystyle {\begin{aligned}Lc(z)&=\int _{0}^{z}\pi x\cot \pi x\,dx\\&=z\log(\sin \pi z)-\int _{0}^{z}\log(\sin \pi x)\,dx\\&=z\log(\sin \pi z)-\int _{0}^{z}\left[\log(2\sin \pi x)-\log 2\right]\,dx\\&=z\log(2\sin \pi z)-\int _{0}^{z}\log(2\sin \pi x)\,dx\end{aligned}}} から、積分変数の置換 y = 2 π x ⟹ d x = d y / ( 2 π ) {\displaystyle \,y=2\pi x\implies dx=dy/(2\pi )\,} により z log ⁡ ( 2 sin ⁡ π z ) − 1 2 π ∫ 0 2 π z log ⁡ ( 2 sin ⁡ π y 2 ) d y {\displaystyle z\log(2\sin \pi z)-{\frac {1}{2\pi }}\int _{0}^{2\pi z}\log \left(2\sin \pi {\frac {y}{2}}\right)\,dy} を得る。二次のクラウセン関数は積分表示 Cl 2 ⁡ ( θ ) = − ∫ 0 θ log ⁡ | 2 sin ⁡ x 2 | d x {\displaystyle \operatorname {Cl} _{2}(\theta )=-\int _{0}^{\theta }\log \left|2\sin {\frac {x}{2}}\right|\,dx} を持つが、0 < θ < 2π なる区間では（積分内の「半正弦函数」の値域は真に正値であるから）被積分函数の絶対値は取り除けて、しかも真に非零である。この定義と、上記の対数正接積分に関する結果とを比較すれば、明らかに L c ( z ) = z log ⁡ ( 2 sin ⁡ π z ) + 1 2 π Cl 2 ( 2 π z ) {\displaystyle Lc(z)=z\log(2\sin \pi z)+{\frac {1}{2\pi }}\,{\text{Cl}}_{2}(2\pi z)} なる関係式が成り立つ。最後に項を並べ替えて、 2 π log ⁡ ( G ( 1 − z ) G ( 1 + z ) ) = 2 π z log ⁡ ( sin ⁡ π z π ) + Cl 2 ( 2 π z ) {\displaystyle 2\pi \log \left({\frac {G(1-z)}{G(1+z)}}\right)=2\pi z\log \left({\frac {\sin \pi z}{\pi }}\right)+{\text{Cl}}_{2}(2\pi z)} とすれば証明は完了する。□ G ( 1 + z ) = Γ ( z ) G ( z ) {\displaystyle \,G(1+z)=\Gamma (z)\,G(z)\,} なる関係を使い、反射公式を 2 π {\displaystyle \,2\pi \,} で割れば、 log ⁡ ( G ( 1 − z ) G ( z ) ) = z log ⁡ ( sin ⁡ π z π ) + log ⁡ Γ ( z ) + 1 2 π Cl 2 ( 2 π z ) {\displaystyle \log \left({\frac {G(1-z)}{G(z)}}\right)=z\log \left({\frac {\sin \pi z}{\pi }}\right)+\log \Gamma (z)+{\frac {1}{2\pi }}{\text{Cl}}_{2}(2\pi z)} もわかる。 反射式 (1) と同等の式に、ベルヌーイ多項式を用いた式 log ⁡ ( G ( 1 2 + z ) G ( 1 2 − z ) ) = log ⁡ Γ ( 1 2 − z ) + B 1 ( z ) log ⁡ 2 π − 1 2 log ⁡ 2 + π ∫ 0 z B 1 ( x ) tan ⁡ π x d x {\displaystyle \log \left({\frac {G\left({\frac {1}{2}}+z\right)}{G\left({\frac {1}{2}}-z\right)}}\right)=\log \Gamma \left({\frac {1}{2}}-z\right)+B_{1}(z)\log 2\pi -{\frac {1}{2}}\log 2+\pi \int _{0}^{z}B_{1}(x)\tan \pi x\,dx} (2) がある。zを(1/2) − z''に置き換えるとこの式は上に等しい。 テイラーの定理とバーンズのG関数の対数微分により以下の級数展開が分かる。 log ⁡ G ( 1 + z ) = z 2 log ⁡ 2 π − ( z + ( 1 + γ ) z 2 2 ) + ∑ k = 2 ∞ ( − 1 ) k ζ ( k ) k + 1 z k + 1 . {\displaystyle \log G(1+z)={\frac {z}{2}}\log 2\pi -\left({\frac {z+(1+\gamma )z^{2}}{2}}\right)+\sum _{k=2}^{\infty }(-1)^{k}{\frac {\zeta (k)}{k+1}}z^{k+1}.} これは 0 < z < 1 {\displaystyle \,0<z<1\,} において有効であり、ここで ζ ( x ) {\displaystyle \,\zeta (x)\,} はリーマンゼータ関数 ζ ( x ) = ∑ k = 1 ∞ 1 k x {\displaystyle \zeta (x)=\sum _{k=1}^{\infty }{\frac {1}{k^{x}}}} である。級数の両辺を指数関数に代入すると、 G ( 1 + z ) = exp ⁡ [ z 2 log ⁡ 2 π − ( z + ( 1 + γ ) z 2 2 ) + ∑ k = 2 ∞ ( − 1 ) k ζ ( k ) k + 1 z k + 1 ] = ( 2 π ) z / 2 exp ⁡ [ − z + ( 1 + γ ) z 2 2 ] exp ⁡ [ ∑ k = 2 ∞ ( − 1 ) k ζ ( k ) k + 1 z k + 1 ] {\displaystyle {\begin{aligned}G(1+z)&=\exp \left[{\frac {z}{2}}\log 2\pi -\left({\frac {z+(1+\gamma )z^{2}}{2}}\right)+\sum _{k=2}^{\infty }(-1)^{k}{\frac {\zeta (k)}{k+1}}z^{k+1}\right]\\&=(2\pi )^{z/2}\exp \left[-{\frac {z+(1+\gamma )z^{2}}{2}}\right]\exp \left[\sum _{k=2}^{\infty }(-1)^{k}{\frac {\zeta (k)}{k+1}}z^{k+1}\right]\end{aligned}}} となる。ここからワイエルシュトラスの乗積表示の形との比較に関し以下が得られる。 exp ⁡ [ ∑ k = 2 ∞ ( − 1 ) k ζ ( k ) k + 1 z k + 1 ] = ∏ k = 1 ∞ { ( 1 + z k ) k exp ( z 2 2 k − z ) } . {\displaystyle \exp \left[\sum _{k=2}^{\infty }(-1)^{k}{\frac {\zeta (k)}{k+1}}z^{k+1}\right]=\prod _{k=1}^{\infty }\left\{\left(1+{\frac {z}{k}}\right)^{k}{\text{exp}}\left({\frac {z^{2}}{2k}}-z\right)\right\}.} ガンマ関数と同様にバーンズのG関数は引数の整数倍に関して以下の公式を有する。 G ( n z ) = K ( n ) n n 2 z 2 / 2 − n z ( 2 π ) − n 2 − n 2 z ∏ i = 0 n − 1 ∏ j = 0 n − 1 G ( z + i + j n ) {\displaystyle G(nz)=K(n)n^{n^{2}z^{2}/2-nz}(2\pi )^{-{\frac {n^{2}-n}{2}}z}\prod _{i=0}^{n-1}\prod _{j=0}^{n-1}G\left(z+{\frac {i+j}{n}}\right)} ここで K ( n ) {\displaystyle K(n)} は以下で与えられる。 K ( n ) = e − ( n 2 − 1 ) ζ ′ ( − 1 ) ⋅ n 5 12 ⋅ ( 2 π ) ( n − 1 ) / 2 = ( A e − 1 12 ) n 2 − 1 ⋅ n 5 12 ⋅ ( 2 π ) ( n − 1 ) / 2 . {\displaystyle K(n)=e^{-(n^{2}-1)\zeta ^{\prime }(-1)}\cdot n^{\frac {5}{12}}\cdot (2\pi )^{(n-1)/2}\,=\,(Ae^{-{\frac {1}{12}}})^{n^{2}-1}\cdot n^{\frac {5}{12}}\cdot (2\pi )^{(n-1)/2}.} ここで ζ ′ {\displaystyle \zeta ^{\prime }} はリーマンゼータ関数の導関数、 A {\displaystyle A} はグレイシャーの定数である。 バーンズの示した通り、G(z + 1)の対数は log ⁡ G ( z + 1 ) = 1 12 − log ⁡ A + z 2 log ⁡ 2 π + ( z 2 2 − 1 12 ) log ⁡ z − 3 z 2 4 + ∑ k = 1 N B 2 k + 2 4 k ( k + 1 ) z 2 k   +   O ( 1 z 2 N + 2 ) {\displaystyle {\begin{aligned}\log G(z+1)&={\frac {1}{12}}-\log A+{\frac {z}{2}}\log 2\pi +\left({\frac {z^{2}}{2}}-{\frac {1}{12}}\right)\log z\\&\quad -{\frac {3z^{2}}{4}}+\sum _{k=1}^{N}{\frac {B_{2k+2}}{4k\left(k+1\right)z^{2k}}}~+~O\left({\frac {1}{z^{2N+2}}}\right)\end{aligned}}} と漸近展開される。 ここで B k {\displaystyle B_{k}} はベルヌーイ数であり、 A {\displaystyle A} はグレイシャーの定数である。（バーンズの時代には少し紛らわしいことにベルヌイ数 B 2 k {\displaystyle B_{2k}} が ( − 1 ) k + 1 B k {\displaystyle (-1)^{k+1}B_{k}} と書かれることもあったが、現在ではそのような規約は用いられないことに注意）。この漸近展開は |z| が大きいとき、負の実軸を含まない任意の扇形に属する z に対して成り立つ。 対数ガンマの媒介変数表示はバーンズ G-函数を用いて ∫ 0 z log ⁡ Γ ( x ) d x = z ( 1 − z ) 2 + z 2 log ⁡ 2 π + z log ⁡ Γ ( z ) − log ⁡ G ( 1 + z ) {\displaystyle \int _{0}^{z}\log \Gamma (x)\,dx={\frac {z(1-z)}{2}}+{\frac {z}{2}}\log 2\pi +z\log \Gamma (z)-\log G(1+z)} と評価することができる。 その証明は少々間接的である。まずはガンマ函数と G-函数との対数差分 z log ⁡ Γ ( z ) − log ⁡ G ( 1 + z ) {\displaystyle z\log \Gamma (z)-\log G(1+z)} を調べる。ここで、 1 Γ ( z ) = z e γ z ∏ k = 1 ∞ { ( 1 + z k ) e − z / k } {\displaystyle {\frac {1}{\Gamma (z)}}=ze^{\gamma z}\prod _{k=1}^{\infty }\left\{\left(1+{\frac {z}{k}}\right)e^{-z/k}\right\}} であり、γ はオイラーの定数である。 バーンズ函数とガンマ函数に関してヴァイヤストラスの乗積形の対数をとることで z log ⁡ Γ ( z ) − log ⁡ G ( 1 + z ) = − z log ⁡ ( 1 Γ ( z ) ) − log ⁡ G ( 1 + z ) = − z [ log ⁡ z + γ z + ∑ k = 1 ∞ { log ⁡ ( 1 + z k ) − z k } ] − [ z 2 log ⁡ 2 π − z 2 − z 2 2 − z 2 γ 2 + ∑ k = 1 ∞ { k log ⁡ ( 1 + z k ) + z 2 2 k − z } ] {\displaystyle {\begin{aligned}&z\log \Gamma (z)-\log G(1+z)=-z\log \left({\frac {1}{\Gamma (z)}}\right)-\log G(1+z)\\=&-z\left[\log z+\gamma z+\sum _{k=1}^{\infty }{\Bigg \{}\log \left(1+{\frac {z}{k}}\right)-{\frac {z}{k}}{\Bigg \}}\right]\\&\qquad -\left[{\frac {z}{2}}\log 2\pi -{\frac {z}{2}}-{\frac {z^{2}}{2}}-{\frac {z^{2}\gamma }{2}}+\sum _{k=1}^{\infty }{\Bigg \{}k\log \left(1+{\frac {z}{k}}\right)+{\frac {z^{2}}{2k}}-z{\Bigg \}}\right]\end{aligned}}} となり、少し整理して項を並べ替えれば級数展開 ∑ k = 1 ∞ { ( k + z ) log ⁡ ( 1 + z k ) − z 2 2 k − z } = − z log ⁡ z − z 2 log ⁡ 2 π + z 2 + z 2 2 − z 2 γ 2 − z log ⁡ Γ ( z ) + log ⁡ G ( 1 + z ) {\displaystyle {\begin{aligned}&\sum _{k=1}^{\infty }{\Bigg \{}(k+z)\log \left(1+{\frac {z}{k}}\right)-{\frac {z^{2}}{2k}}-z{\Bigg \}}\\&\qquad =-z\log z-{\frac {z}{2}}\log 2\pi +{\frac {z}{2}}+{\frac {z^{2}}{2}}-{\frac {z^{2}\gamma }{2}}-z\log \Gamma (z)+\log G(1+z)\end{aligned}}} を得る。最後に、対数ガンマ函数のヴァイヤストラス乗積形をとって区間 [0,z] 上積分すれば ∫ 0 z log ⁡ Γ ( x ) d x = − ∫ 0 z log ⁡ ( 1 Γ ( x ) ) d x = − ( z log ⁡ z − z ) − z 2 γ 2 − ∑ k = 1 ∞ { ( k + z ) log ⁡ ( 1 + z k ) − z 2 2 k − z } {\displaystyle {\begin{aligned}\int _{0}^{z}\log \Gamma (x)\,dx&=-\int _{0}^{z}\log \left({\frac {1}{\Gamma (x)}}\right)\,dx\\&=-(z\log z-z)-{\frac {z^{2}\gamma }{2}}-\sum _{k=1}^{\infty }{\Bigg \{}(k+z)\log \left(1+{\frac {z}{k}}\right)-{\frac {z^{2}}{2k}}-z{\Bigg \}}\end{aligned}}} となる。二つの評価を等しいと置いて ∫ 0 z log ⁡ Γ ( x ) d x = z ( 1 − z ) 2 + z 2 log ⁡ 2 π + z log ⁡ Γ ( z ) − log ⁡ G ( 1 + z ) {\displaystyle \int _{0}^{z}\log \Gamma (x)\,dx={\frac {z(1-z)}{2}}+{\frac {z}{2}}\log 2\pi +z\log \Gamma (z)-\log G(1+z)} の証明は完成する。□