文献あり

楕円積分の特殊値を求める

267

はじめに

　この記事では楕円積分 $K (k)$ に対して
$\frac{K^{'} (k_{N})}{K (k_{N})} = \sqrt{N}$
が成り立つときの $K (k_{N})$ の値について解説していきます。
　前回の記事ではLattice Sumの $L$ 関数への分解
$\underset{m, n}{\sum^{'}} \frac{1}{(m^{2} + N n^{2})^{s}} = \sum_{2 ∤ d ∣ N} L_{\pm d} L_{\mp 4 N / d}$
を考えることでラマヌジャンの不変量 $G_{N}, g_{N}$ を求めたのでした。
　実はこの公式からはラマヌジャンの不変量だけでなく $τ = \sqrt{- N}$ における
$K = \frac{π}{2} θ_{3} (τ)^{2}$
の値を求めることもできます。万能かよ。

アイゼンシュタイン級数とクロネッカーの極限公式

$\underset{m, n}{\sum^{'}} \frac{1}{(m^{2} + N n^{2})^{s}} = \sum_{2 ∤ d ∣ N} L_{\pm d} L_{\mp 4 N / d}$
といった分解公式を考えていたZuckerやGlasserらはこの左辺の級数をLattice Sumと呼んでおり、一般に正定値の二次形式 $Q = a x^{2} + b x y + c y^{2}$ に対し
$ζ_{Q} (s) = \underset{m, n}{\sum^{'}} \frac{1}{Q (m, n)^{s}}$
という級数の $L$ 関数への分解を考えていました。この関数はEpsteinのゼータ関数とも呼ばれます。
　 $Q$ は正定値であることから $D = b^{2} - 4 a c < 0$ なのでこれは
$Q = a (m + \frac{b + \sqrt{D}}{2 a} n) (m + \frac{b - \sqrt{D}}{2 a} n) = a | m + n τ |^{2}$
ただし
$τ = \frac{b + \sqrt{D}}{2 a}$
と変形することができ、Epsteinのゼータ関数は
$ζ_{Q} (s) = \underset{m, n}{\sum^{'}} \frac{1}{a^{s} | m τ + n |^{2 s}}$
と表すことができます。
　一般に $τ \in H$ として $a^{s}$ および $y^{s} = Im (τ)^{s}$ を乗じたもの
$E (τ, s) = \underset{m, n}{\sum^{'}} \frac{y^{s}}{| m τ + n |^{2 s}}$
を実解析的アイゼンシュタイン級数と呼びます。流儀によってはこれを $2 ζ (2 s)$ で除したもの
$E^{'} (τ, s) = \frac{1}{2} \sum_{gcd (m, n) = 1} \frac{y^{s}}{| m τ + n |^{2 s}}$
も実解析的アイゼンシュタイン級数と言いますが、この記事では扱いません。
　さてLattice Sumを一般化した実解析的アイゼンシュタイン級数には次のような公式があります。

クロネッカーの極限公式

　 $η$ をデデキントのイータ関数
$η (τ) = q^{- \frac{1}{12}} \prod_{n = 1}^{\infty} (1 - q^{2 n}) (q = e^{i π τ})$
とすると
$lim_{s \to 1} (E (τ, s) - \frac{π}{s - 1}) = 2 π (γ - \log 2 - \log (\sqrt{y} | η (τ) |^{2}))$
が成り立つ。

イータ関数の特殊値

　クロネッカーの極限公式については証明しませんが、この $τ = \sqrt{- N}$ の場合を考えることで以下の公式が得られます。

$lim_{s \to 1} (\frac{\sqrt{N}}{π} \underset{m, n}{\sum^{'}} \frac{1}{(m^{2} + N n^{2})^{s}} - \frac{1}{s - 1}) = 2 (γ - \log 2 - \frac{1}{2} \log N - \log η (\sqrt{- N})^{2})$
が成り立つ。

$\begin{array}{rcl} lim_{s \to 1} (\underset{m, n}{\sum^{'}} \frac{N^{\frac{s}{2}}}{(m^{2} + N n^{2})^{s}} - \frac{π}{s - 1}) \\ = & lim_{s \to 1} ((\underset{m, n}{\sum^{'}} \frac{\sqrt{N}}{(m^{2} + N n^{2})^{s}} - \frac{π}{s - 1}) \frac{N^{\frac{s}{2}}}{N^{\frac{1}{2}}} + \frac{π}{s - 1} (\frac{N^{\frac{s}{2}}}{N^{\frac{1}{2}}} - 1)) \\ = & lim_{s \to 1} (\underset{m, n}{\sum^{'}} \frac{\sqrt{N}}{(m^{2} + N n^{2})^{s}} - \frac{π}{s - 1}) + \frac{π}{2} \log N \end{array}$
に注意するとわかる。

　ここで前回の記事の定理5を思い出しましょう。

　 $(\frac{a}{n})$ をクロネッカー記号とし
$L_{d} (s) = \sum_{n = 1}^{\infty} (\frac{d}{n}) \frac{1}{n^{s}}$
とおく。このとき $N ≢ 3 (\mod 4)$ なる平方因子を持たない便利数 $N \neq 1$ に対してその素因数 $p \geq 3$ の個数を $r$ とおくと
$\underset{m, n}{\sum^{'}} \frac{1}{(m + N n)^{s}} = \frac{1}{2^{r - 1}} \sum_{2 ∤ d ∣ N} L_{\pm d} L_{\mp 4 N / d}$
が成り立つ。ただし符号は $\pm d \equiv 1 (\mod 4)$ となるように取るものとした。

　この公式とクロネッカーの極限公式を合わせることで以下のように $η (\sqrt{- N})$ を計算することができます。

　定理2の条件下で
$\log η (\sqrt{- N})^{2} = \frac{1}{2^{r + 1}} \sum_{n = 1}^{4 N} (\frac{- 4 N}{n}) \log Γ (\frac{n}{4 N}) - \frac{\sqrt{N}}{2^{r} π} \sum_{\begin{matrix} 2 ∤ d ∣ N \\ d \neq 1 \end{matrix}} L_{\pm d} (1) L_{\mp 4 N / d} (1) - \frac{1}{2} \log 8 π N$
が成り立つ。

　この証明は後半に回すとして、まずこの公式の効力を見ていきましょう。
　 $Γ$ 関数の倍数公式
$Γ (n z) = \frac{n^{n z - \frac{1}{2}}}{(2 π)^{\frac{n - 1}{2}}} \prod_{k = 0}^{n - 1} Γ (z + \frac{k}{n})$
から得られる等式
$\sum_{\begin{matrix} n = 1 \\ (n, 4 N) = 1 \end{matrix}}^{4 N} \log Γ (\frac{n}{4 N}) = \frac{φ (4 N)}{2} \log 2 π (N \neq 2)$
(これも後半で証明する)を用いて定理3式を整理すると
$\log η (\sqrt{- N})^{2} = \frac{1}{2^{r}} \sum_{\begin{matrix} n = 1 \\ (n, 4 N) = 1 \end{matrix}}^{4 N} \frac{1 + (- 4 N | n)}{2} \log Γ (\frac{n}{4 N}) - \frac{\sqrt{N}}{2^{r} π} \sum_{\begin{matrix} 2 ∤ d ∣ N \\ d \neq 1 \end{matrix}} L_{\pm d} (1) L_{\mp 4 N / d} (1) - \frac{1}{2} \log 8 π N - \frac{φ (4 N)}{2^{r + 2}} \log 2 π$
となり、さらに
$M_{D} = {\begin{array}{cl} h_{D} \log ε_{D} & D > 0 \\ 2 h_{D} / w_{D} & D < 0 \end{array}$
とおくと二次体の類数公式から次のように整理できます。

$\log η (\sqrt{- N})^{2} = \frac{1}{2^{r}} \sum_{\begin{matrix} n = 1 \\ (n, 4 N) = 1 \end{matrix}}^{4 N} \frac{1 + (- 4 N | n)}{2} \log Γ (\frac{n}{4 N}) - \frac{1}{2^{r}} \sum_{\begin{matrix} 2 ∤ d ∣ N \\ d \neq 1 \end{matrix}} M_{\pm d} M_{\mp 4 N / d} - \frac{1}{2} \log 8 π N - \frac{φ (4 N)}{2^{r + 2}} \log 2 π$

となりまず。
　例えば $N = 5, 13 (r = 1)$ の場合を考えると
$\begin{array}{rcl} \log η (\sqrt{- 5})^{2} & = & \frac{1}{2} \log (Γ (\frac{1}{20}) Γ (\frac{3}{20}) Γ (\frac{7}{20}) Γ (\frac{9}{20})) - \frac{1}{4} \log \frac{\sqrt{5} + 1}{2} - \frac{1}{2} \log 40 π - \log 2 π \\ η (\sqrt{- 5})^{2} & = & \frac{1}{π} {(\frac{\sqrt{5} + 1}{2})}^{- \frac{1}{4}} {(\frac{Γ (\frac{1}{20}) Γ (\frac{3}{20}) Γ (\frac{7}{20}) Γ (\frac{9}{20})}{160 π})}^{\frac{1}{2}} \\ \log η (\sqrt{- 13})^{2} & = & \frac{1}{2} \sum_{a \in A} \log Γ (\frac{a}{52}) - \frac{1}{4} \log \frac{\sqrt{13} + 3}{2} - \frac{1}{2} \log 104 π - 3 \log 2 π \\ η (\sqrt{- 13})^{2} & = & \frac{1}{π} {(\frac{\sqrt{13} + 3}{2})}^{- \frac{1}{4}} {(\frac{\prod_{a \in A} Γ (\frac{a}{52})}{6656 π^{5}})}^{\frac{1}{2}} \\ (A & = & {1, 7, 9, 11, 15, 17, 19, 25, 29, 31, 47, 49}) \end{array}$
と計算できるわけです。やば。

楕円積分の特殊値

　そしてイータ関数の特殊値とラマヌジャンの不変量
$\begin{array}{rcl} f (τ) & = & q^{- \frac{1}{24}} \prod_{n = 1}^{\infty} (1 + q^{2 n - 1}) (q = e^{π i τ}) \\ G_{N} & = & 2^{- \frac{1}{4}} f (\sqrt{- N}) \end{array}$
を用いるとこの記事の定理3の証明から $τ = \sqrt{- N}$ において
$\begin{array}{rcl} K (k_{N}) & = & \frac{π}{2} θ_{3}^{2} \\ = & \frac{π}{2} η (\sqrt{- N})^{2} f (\sqrt{- N})^{4} \\ = & π η (\sqrt{- N})^{2} G_{N}^{4} \end{array}$
と楕円積分の特殊値を求めることができます。
　例えば $N = 5, 13$ の場合を考えると
$G_{5}^{4} = \frac{\sqrt{5} + 1}{2}, G_{13}^{4} = \frac{\sqrt{13} + 3}{2}$
であったので上の結果から
$\begin{array}{rcl} K (k_{5}) & = & (\sqrt{5} + 2)^{\frac{1}{4}} {(\frac{Γ (\frac{1}{20}) Γ (\frac{3}{20}) Γ (\frac{7}{20}) Γ (\frac{9}{20})}{160 π})}^{\frac{1}{2}} \\ K (k_{13}) & = & (5 \sqrt{13} + 18)^{\frac{1}{4}} {(\frac{Γ (\frac{1}{52}) Γ (\frac{7}{52}) Γ (\frac{9}{52}) Γ (\frac{11}{52}) Γ (\frac{15}{52}) Γ (\frac{17}{52}) Γ (\frac{19}{52}) Γ (\frac{25}{52}) Γ (\frac{29}{52}) Γ (\frac{31}{52}) Γ (\frac{47}{52}) Γ (\frac{49}{52})}{6656 π^{5}})}^{\frac{1}{2}} \end{array}$
と計算できます。なんてこった。
　ここでは紹介しませんが他の便利数 $N$ に対しても同様に $η (\sqrt{- N})$ および $K (k_{N})$ を求めることができます(ただ $N$ が偶数のときの $G_{N}$ の値は前回の記事とは別の方法で求めておく必要があります)。またいくつかの $N$ に対する $K (k_{N})$ の値については Wolfram MathWorld に列記されています。

定理3の証明

　クロネッカーの極限公式および定理2より
$\begin{array}{rcl} 4 \log η (\sqrt{- N}) & = & 2 γ - \log 4 N - lim_{s \to 1} (\frac{\sqrt{N}}{2^{r - 1} π} ζ (s) L_{- 4 N} (s) - \frac{1}{s - 1}) - \frac{\sqrt{N}}{2^{r - 1} π} \sum_{\begin{array}{c} 2 ∤ d ∣ N \\ d \neq 1 \end{array}} L_{\pm d} (1) L_{\mp N / d} (1) \end{array}$
が成り立つので、前回の記事での議論から
$L_{- 4 N} (1) = \frac{2^{r - 1} π}{\sqrt{N}}$
であったことに注意すると以下の等式を示せばよい。

$lim_{s \to 1} (ζ (s) \frac{L_{- 4 N} (s)}{L_{- 4 N} (1)} - \frac{1}{s - 1}) = 2 γ + \log 2 π - \frac{1}{2^{r}} \sum_{n = 1}^{4 N} (\frac{- 4 N}{n}) \log Γ (\frac{n}{4 N})$
が成り立つ。

$lim_{s \to 1} (ζ (s) - \frac{1}{s - 1}) = γ$
であることが知られている(後述)ので
$\begin{array}{rcl} lim_{s \to 1} (ζ (s) \frac{L_{- 4 N} (s)}{L_{- 4 N} (1)} - \frac{1}{s - 1}) & = & lim_{s \to 1} ((ζ (s) - \frac{1}{s - 1}) \frac{L_{- 4 N} (s)}{L_{- 4 N} (1)} + \frac{1}{s - 1} (\frac{L_{- 4 N} (s)}{L_{- 4 N} (1)} - 1)) \\ = & γ + \frac{L_{- 4 N}^{'} (1)}{L_{- 4 N} (1)} \end{array}$
が成り立つ。
　また虚二次体の判別式 $D$ に対してクロネッカー記号 $(D | n)$ は法 $D$ の原始的奇ディリクレ指標を定めることが知られており、そのガウス和の符号は
$τ ((\frac{D}{\cdot})) = \sum_{n = 1}^{| D |} (\frac{D}{n}) ζ_{D}^{n} = \sqrt{D}$
と決定されるらしいのでこの記事より
$\begin{array}{rcl} \frac{L_{- 4 N}^{'} (1)}{L_{- 4 N} (1)} & = & γ + \log 2 π - \frac{π i}{τ ((- 4 N | \cdot)) L_{- 4 N} (1)} \sum_{n = 1}^{4 N} (\frac{- 4 N}{n}) \log Γ (\frac{n}{4 N}) \\ = & γ + \log 2 π - \frac{π \sqrt{- 1}}{\sqrt{- 4 N}} \cdot \frac{\sqrt{N}}{2^{r - 1} π} \sum_{n = 1}^{4 N} (\frac{- 4 N}{n}) \log Γ (\frac{n}{4 N}) \\ = & γ + \log 2 π - \frac{1}{2^{r}} \sum_{n = 1}^{4 N} (\frac{- 4 N}{n}) \log Γ (\frac{n}{4 N}) \end{array}$
が成り立ち、主張を得る。

おまけ

$lim_{s \to 1} (ζ (s) - \frac{1}{s - 1}) = γ$
が成り立つ。

　リーマンゼータ関数の関数等式
$ζ (s) = 2^{s} π^{s - 1} \sin (\frac{π s}{2}) Γ (1 - s) ζ (1 - s)$
において $s \to 1$ の挙動を考えると
$\begin{array}{rcl} Γ (1 - s) & = & \frac{Γ (2 - s)}{1 - s} \\ = & \frac{Γ (1) - Γ^{'} (1) (s - 1)}{1 - s} + O (s - 1) \\ = & \frac{1}{1 - s} - γ + O (s - 1) \end{array}$
なので
$f (s) = 2^{s} π^{s - 1} \sin (\frac{π s}{2}) ζ (1 - s)$
とおくと $f (1) = - 1$ より
$\begin{array}{rcl} lim_{s \to 1} (ζ (s) - \frac{1}{s - 1}) & = & lim_{s \to 1} (- (Γ (s) + \frac{1}{s - 1}) \frac{f (s)}{f (1)} + \frac{1}{s - 1} (\frac{f (s)}{f (1)} - 1)) \\ = & γ + \frac{f^{'} (1)}{f (1)} \\ = & γ + \log 2 π - \frac{ζ^{'} (0)}{ζ (0)} = γ \end{array}$
とわかる(最後の等号については昔の記事のコピペを貼っておきます)。

\frac{ζ^{'} (0)}{ζ (0)} = \log 2 π

の証明

$η (x) = (1 - 2^{1 - s}) ζ (s) = \sum_{n = 1}^{\infty} \frac{1}{n^{s}} - 2 \sum_{n = 1}^{\infty} \frac{1}{(2 n)^{s}} = \sum_{n = 1}^{\infty} \frac{(- 1)^{n - 1}}{n^{s}}$
とおくとこれは $Re (s) > 0$ で収束するが、
$2 η (s) = \sum_{n = 1}^{\infty} \frac{(- 1)^{n - 1}}{n^{s}} + (1 - \sum_{n = 1}^{\infty} \frac{(- 1)^{n - 1}}{(n + 1)^{s}}) = 1 + \sum_{n = 1}^{\infty} (- 1)^{n - 1} (\frac{1}{n^{s}} - \frac{1}{(n + 1)^{s}})$
と変形すれば $Re (s) > - 1$ で収束することになる。

収束性について

$\begin{array}{rcl} | \sum_{n = 1}^{\infty} (- 1)^{n - 1} (\frac{1}{n^{s}} - \frac{1}{(n + 1)^{s}}) | \\ = & | \sum_{n = 1}^{\infty} ((\frac{1}{(2 n - 1)^{s}} - \frac{1}{(2 n)^{s}}) - (\frac{1}{(2 n)^{s}} - \frac{1}{(2 n + 1)^{s}})) | \\ = & | \sum_{n = 1}^{\infty} s (\int_{2 n - 1}^{2 n} \frac{d x}{x^{s + 1}} - \int_{2 n - 1}^{2 n} \frac{d x}{(x + 1)^{s - 1}}) | \\ = & | \sum_{n = 1}^{\infty} s (s + 1) \int_{2 n - 1}^{2 n} \int_{x}^{x + 1} \frac{1}{y^{s + 2}} d y d x | \\ \leq & | s (s + 1) | \sum_{n = 1}^{\infty} \frac{1}{(2 n - 1)^{Re (s) + 2}} < \infty \end{array}$
とわかる。

　またその項別微分

2 η^{'} (s) = \sum_{n = 1}^{\infty} (- 1)^{n} (\frac{\log n}{n^{s}} - \frac{\log (n + 1)}{n^{s}})

も

Re (s) > - 1

で収束するので

s = 0

を代入することで

\begin{array}{rcl} 2 η^{'} (0) & = & \sum_{n = 1}^{\infty} (- (\log (2 n - 1) - \log (2 n)) + (\log (2 n) - \log (2 n + 1))) \\ = & \log (\prod_{n = 1}^{\infty} \frac{(2 n)^{2}}{(2 n - 1) (2 n + 1)}) \\ = & \log \frac{π}{2} \end{array}

を得る。ここで

\sin \frac{π}{2} = \frac{π}{2} \prod_{n = 1}^{\infty} (1 - \frac{1}{(2 n)^{2}}) = \frac{π}{2} \prod_{n = 1}^{\infty} \frac{(2 n - 1) (2 n + 1)}{(2 n)^{2}} = 1

であること(ウォリスの公式)を用いた。

　いま

2 η^{'} (s) = 2 ((1 - 2^{1 - s}) ζ (s))^{'} = 2^{1 - s} (\log 2) 2 ζ (s) + 2 (1 - 2^{1 - s}) ζ^{'} (s)

であったので

ζ (0) = - \frac{1}{2}

に注意して

s = 0

を代入することで

2 η^{'} (0) = - 2 \log 2 - 2 ζ^{'} (0) = \log \frac{π}{2}

すなわち

ζ^{'} (0) = - \frac{1}{2} \log 2 π

ひいては

\frac{ζ^{'} (0)}{ζ (0)} = \log 2 π

を得る。

$\sum \log Γ (n / 4 N)$ の計算

　 $n$ が素数べき $p^{e}$ のとき
$\sum_{\begin{matrix} k = 1 \\ (k, n) = 1 \end{matrix}}^{n} \log Γ (\frac{k}{n}) = \frac{φ (n)}{2} \log 2 π - \frac{1}{2} \log p$
　 $n$ が二つ以上の素因数を持つとき
$\sum_{\begin{matrix} k = 1 \\ (k, n) = 1 \end{matrix}}^{n} \log Γ (\frac{k}{n}) = \frac{φ (n)}{2} \log 2 π$
が成り立つ。

　求める式を $a_{n}$ とおくと
$\begin{array}{rcl} \sum_{d | n} a_{d} & = & \sum_{d | n} a_{n / d} \\ = & \sum_{d | n} \sum_{\begin{array}{c} k = 1 \\ (k, n / d) = 1 \end{array}}^{n / d} \log Γ (\frac{d k}{n}) \\ = & \sum_{d | n} \sum_{\begin{array}{c} k = 1 \\ (k, n) = d \end{array}}^{n / d} \log Γ (\frac{k}{n}) \\ = & \sum_{k = 1}^{n} \log Γ (\frac{k}{n}) \end{array}$
が成り立つ。これはガンマ関数の倍数公式
$Γ (n z) = \frac{n^{n z - \frac{1}{2}}}{(2 π)^{\frac{n - 1}{2}}} \prod_{k = 0}^{n - 1} Γ (z + \frac{k}{n})$
において $z = 1 / n$ を代入することで
$\sum_{k = 1}^{n} \log Γ (\frac{k}{n}) = \frac{n - 1}{2} \log 2 π - \frac{1}{2} \log n$
と計算できるのでメビウスの反転公式より
$a_{n} = \sum_{d | n} μ (\frac{n}{d}) (\frac{d - 1}{2} \log 2 π - \frac{1}{2} \log d)$
を得る。また
$\begin{array}{rcl} \sum_{d | n} μ (d) & = & 0 (n > 1) \\ \sum_{d | n} μ (\frac{n}{d}) d & = & φ (n) \end{array}$
より
$a_{n} = \frac{φ (n)}{2} \log 2 π - \frac{1}{2} \sum_{d | n} μ (\frac{n}{d}) \log d$
と計算できる。
　ここで
$\log b_{n} = \sum_{d | n} μ (\frac{n}{d}) \log d$
とおくとメビウスの反転公式より
$n = \prod_{d | n} b_{d}$
が成り立ち、ここから $b_{1} = 1$ および $b_{p^{e}} = p$ がわかる。また
$n = \prod_{p ∣ n} p^{e_{p}}$
と素因数分解したとき、
$\prod_{p | n} \prod_{e = 1}^{e_{p}} b_{p^{e}} = n$
となることに注意すると数学的帰納法より $n$ が二つ以上の素因数を持つときは $b_{n} = 1$ となることがわかる。
　以上より主張を得る。

参考文献

[1]

J. M. Borwein, P. B. Borwein, Pi and the AGM: A Study in Analytic Number Theory and Computational Complexity, Wiley-Interscience, 1987, pp. 296-298

投稿日：2023年2月15日

この記事を高評価した人

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

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

バッジはありません。

バッチを贈って投稿者を応援しよう

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

投稿者

子葉

1105

274415

主に複素解析、代数学、数論を学んでおります。私の経験上、その証明が簡単に探しても見つからない、英語の文献を漁らないと載ってない、なんて定理の解説を主にやっていきます。同じ経験をしている人の助けになれば。最近は自分用のノートになっている節があります。

他の人のコメント

コメントはありません。

読み込み中

子葉

楕円積分の特殊値を求める

はじめに
アイゼンシュタイン級数とクロネッカーの極限公式
イータ関数の特殊値
楕円積分の特殊値
定理3の証明
$\sum\log\G(n/4N)$の計算
参考文献

楕円積分の特殊値を求める

はじめに

アイゼンシュタイン級数とクロネッカーの極限公式

イータ関数の特殊値

楕円積分の特殊値

定理3の証明

おまけ

∑log⁡Γ(n/4N)の計算

参考文献

この記事を高評価した人

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

投稿者

コメント

他の人のコメント

$\sum \log Γ (n / 4 N)$ の計算