Eulerの定数(Euler's constant)とは、調和数 $H_n=1+\frac12+\cdots+\frac1n$ と自然対数の差の極限 $\gamma=\lim_{n\to\infty}(H_n-\log n)=0.5772\ldots$ のことで、Euler–Mascheroni 定数ともいう。差 $H_n-\log n-\gamma$ は $\frac1{2(n+1)}$ と $\frac1{2n}$ の間にあり、$\log n$ を $\log(n+\frac12)$ に替えると誤差は $\frac1{24n^2}$ 程度に縮む。$\gamma$ は小数部分の積分やゼータ関数の値の級数 $\sum_{k\ge2}(-1)^k\zeta(k)/k$ で表され、$\Gamma'(1)=-\gamma$ や Mertens の定理にも現れるが、有理数かどうかは知られていない。
調和数 $H_n=1+\frac12+\cdots+\frac1n$ は $n\to\infty$ で発散するが、その増え方は対数 $\log n$ とほとんど同じである。差 $H_n-\log n$ は $n=1,10,100,1000$ で $1,\ 0.6264,\ 0.5822,\ 0.5777$ と減っていき、ある定数 $\gamma=0.5772156649\ldots$ に近づく。この $\gamma$ が Euler の定数 である。$\gamma$ は調和数の漸近式に現れるだけでなく、$\Gamma$ 関数の微分係数 $\Gamma'(1)=-\gamma$、Riemann ゼータ関数の $s=1$ での展開、素数についての Mertens の定理など、解析学と数論のいろいろな場所に顔を出す。
この記事では、極限 $\gamma$ が存在することを示すと同時に、差 $H_n-\log n-\gamma$ が $\frac1{2n}$ とほぼ等しいことを上下の評価で示し、$\log n$ を $\log(n+\frac12)$ に取り替えると誤差が $\frac1{24n^2}$ 程度まで小さくなることを証明する。さらに $\gamma$ の積分表示と、ゼータ関数の値を使った級数表示を導く。
$n\geq1$ に対し $H_n:=\sum_{k=1}^{n}\frac1k$ とおく(第 $n$ 調和数)。$\log$ は自然対数である。
極限
$$
\gamma:=\lim_{n\to\infty}\bigl(H_n-\log n\bigr)
$$
を Euler の定数(Euler's constant)または Euler–Mascheroni の定数(Euler–Mascheroni constant)という。
極限の存在は、次節の thm-euler-constant-rate で収束の速さの評価と同時に示す。別の証明が 調和級数 の記事の定理「Eulerの定数の存在」($1-\log2<\gamma<1$ と、すべての $n\geq1$ で $0< H_n-\log n-\gamma<\frac1n$ も示されている)と、和と積分の差 の系「Euler の定数の存在」にもある。以下 $\varepsilon_n:=H_n-\log n-\gamma$ とおく。同じ定数は Rud76 の第 8 章の演習 9(p. 197)や Leb1 の演習 5.4.6(p. 212)でも、$\lim(H_N-\log N)$ として導入されている。
$H_n-\log n$ は面積の差として見ることができる。区間 $[k,k+1]$ の上に高さ $\frac1k$ の長方形を並べると、その面積の合計は $H_n$ である。一方、$t=1$ から $t=n+1$ までの曲線 $y=\frac1t$ の下の面積は $\log(n+1)$ である。長方形から曲線の下を取り除くと、各区間に「曲がった三角形」が 1 つずつ残る。その面積の合計が $H_n-\log(n+1)$ であり、$n\to\infty$ でこれが $\gamma$ に増加しながら近づく。
区間 $[k,k+1]$ の三角形は、ほぼ直角三角形であり、面積はおよそ $\frac12\bigl(\frac1k-\frac1{k+1}\bigr)$ である。これを $k\geq n$ について足すと $\frac1{2n}$ になる。つまり $\gamma$ から $H_n-\log(n+1)$ を引いた残りも、$H_n-\log n$ から $\gamma$ を引いた残りも、だいたい $\frac1{2n}$ になるはずである。次の節では、この見積もりを、曲線 $y=\frac1t$ が下に凸であることを使って正確な不等式にする。
$x>0$ について
$$
\frac{2x}{2+x}<\log(1+x)<\frac{x(2+x)}{2(1+x)}
$$
が成り立つ。
$\log(1+x)=\int_1^{1+x}\frac{dt}{t}$ である。
右側:$(1,1)$ と $\bigl(1+x,\frac1{1+x}\bigr)$ を結ぶ直線を $\ell(t)=1-\frac{t-1}{1+x}$ とすると
$$
\frac1t-\ell(t)=\frac{(1+x)(1-t)+t(t-1)}{t(1+x)}=\frac{(t-1)(t-1-x)}{t(1+x)}
$$
であり、これは $1< t<1+x$ で負である。よって $\int_1^{1+x}\frac{dt}{t}<\int_1^{1+x}\ell(t)\,dt=x-\frac{x^2}{2(1+x)}=\frac{x(2+x)}{2(1+x)}$ である。
左側:$m:=1+\frac x2$、$h:=\frac x2$ とおく。$0< s< h$ について
$$
\frac1{m+s}+\frac1{m-s}=\frac{2m}{m^2-s^2}>\frac2m
$$
である。積分区間を中点 $m$ で折り返して足すと
$$
\int_1^{1+x}\frac{dt}{t}=\int_0^h\Bigl(\frac1{m+s}+\frac1{m-s}\Bigr)ds>\frac{2h}{m}=\frac{x}{1+x/2}=\frac{2x}{2+x}
$$
を得る。$\square$
左の量は区間の中点での高さを使った長方形の面積(中点公式)、右の量は台形の面積(台形公式)である。下に凸な関数では、積分は中点公式より大きく、台形公式より小さい。
極限 $\gamma=\lim_{n\to\infty}(H_n-\log n)$ は存在し、すべての $n\geq1$ について
$$
\frac1{2(n+1)}< H_n-\log n-\gamma<\frac1{2n}
$$
が成り立つ。特に $n=1$ として $\frac12<\gamma<\frac34$ である。
存在。 $a_k:=H_k-\log k$ とおくと
$$
a_k-a_{k+1}=\log\Bigl(1+\frac1k\Bigr)-\frac1{k+1}
$$
である。下の 2 つの段落の計算(lem-euler-log-midpoint-trapezoid による)から
$$
0<\frac1{(2k+1)(k+1)}< a_k-a_{k+1}<\frac1{2k(k+1)}
$$
である。よって $(a_n)$ は狭義単調減少であり、$a_1=1$ と $\sum_{k\geq1}\frac1{2k(k+1)}=\frac12$ から
$$
a_n=1-\sum_{k=1}^{n-1}(a_k-a_{k+1})>1-\sum_{k=1}^{\infty}\frac1{2k(k+1)}=\frac12
$$
と下に有界である。したがって $(a_n)$ は収束し(単調収束定理)、極限 $\gamma$ が定まる。
$\varepsilon_k=a_k-\gamma$ なので $\varepsilon_k-\varepsilon_{k+1}=a_k-a_{k+1}$ であり、$\varepsilon_N\to0$ だから、$N\to\infty$ として $\varepsilon_n=\sum_{k=n}^{\infty}(\varepsilon_k-\varepsilon_{k+1})$ が成り立つ。
上からの評価:lem-euler-log-midpoint-trapezoid の右側を $x=\frac1k$ で使うと、$\log(1+\frac1k)<\frac{2k+1}{2k(k+1)}$ なので
$$
\varepsilon_k-\varepsilon_{k+1}<\frac{2k+1}{2k(k+1)}-\frac1{k+1}=\frac1{2k(k+1)}=\frac1{2k}-\frac1{2(k+1)}
$$
である。$k=n,n+1,\dots$ について足すと、右辺の和は $\frac1{2n}$ に収束するので $\varepsilon_n<\frac1{2n}$ である(各項で狭義の不等号が成り立つので、和でも狭義である)。
下からの評価:lem-euler-log-midpoint-trapezoid の左側を $x=\frac1k$ で使うと、$\log(1+\frac1k)>\frac2{2k+1}$ なので
$$
\varepsilon_k-\varepsilon_{k+1}>\frac2{2k+1}-\frac1{k+1}=\frac1{(2k+1)(k+1)}>\frac1{2(k+1)(k+2)}=\frac1{2(k+1)}-\frac1{2(k+2)}
$$
である($2k+1<2(k+2)$ による)。$k=n,n+1,\dots$ について足すと $\varepsilon_n>\frac1{2(n+1)}$ を得る。$\square$
thm-euler-constant-rate は $H_n=\log n+\gamma+\frac1{2n}+O\bigl(\frac1{n^2}\bigr)$ を意味する($O$ は Landau記号)。実際、$\frac1{2n}-\varepsilon_n$ は $0$ と $\frac1{2n}-\frac1{2(n+1)}=\frac1{2n(n+1)}$ の間にある。
次の定理は、補正項 $\frac1{2n}$ を対数の中に取り込むと、誤差がもう一段小さくなることを示す。$\log n+\frac1{2n}$ と $\log\bigl(n+\frac12\bigr)=\log n+\frac1{2n}-\frac1{8n^2}+\cdots$ は、$\frac1{n^2}$ の項までしか違わない。
すべての $n\geq1$ について
$$
\frac1{24(n+1)^2}< H_n-\log\Bigl(n+\frac12\Bigr)-\gamma<\frac1{24n^2}
$$
が成り立つ。
$d_k:=H_k-\log(k+\frac12)$ とおく。$d_k=(H_k-\log k)-\log(1+\frac1{2k})\to\gamma$ である。$m:=k+1$、$x:=\frac1{2m}$ とおくと $\frac{m+1/2}{m-1/2}=\frac{1+x}{1-x}$ なので
$$
d_k-d_{k+1}=\log\frac{k+3/2}{k+1/2}-\frac1{k+1}=\log\frac{1+x}{1-x}-2x
$$
である。$|x|<1$ での冪級数展開 $\log\frac{1+x}{1-x}=2\sum_{j\geq0}\frac{x^{2j+1}}{2j+1}$(対数関数)から
$$
d_k-d_{k+1}=2\sum_{j\geq1}\frac{x^{2j+1}}{2j+1}
$$
であり、$0< x\leq\frac14$ なので各項は正である。
下からの評価:最初の項だけを残すと $d_k-d_{k+1}>\frac{2x^3}{3}=\frac1{12m^3}$ である。さらに
$$
\frac1{12(k+1)^3}\geq\frac{2k+3}{24(k+1)^2(k+2)^2}=\frac1{24(k+1)^2}-\frac1{24(k+2)^2}
$$
である。実際、分母を払うと $2(k+2)^2\geq(2k+3)(k+1)$、すなわち $2k^2+8k+8\geq2k^2+5k+3$ となり、これは成り立つ。
上からの評価:$j\geq1$ で $\frac1{2j+1}\leq\frac13$ なので
$$
d_k-d_{k+1}<\frac{2x^3}{3}\sum_{j\geq0}x^{2j}=\frac{2x^3}{3(1-x^2)}=\frac1{3m(4m^2-1)}=\frac1{3(k+1)(2k+1)(2k+3)}
$$
である($j\geq2$ の項で $\frac1{2j+1}<\frac13$ なので狭義)。これが $\frac1{24k^2}-\frac1{24(k+1)^2}=\frac{2k+1}{24k^2(k+1)^2}$ 以下であることは、分母を払った $8k^2(k+1)\leq(2k+1)^2(2k+3)$ から分かる。実際、右辺は $8k^3+20k^2+14k+3$、左辺は $8k^3+8k^2$ である。
$d_n-\gamma=\sum_{k\geq n}(d_k-d_{k+1})$ に 2 つの評価を足し合わせると、それぞれ $\frac1{24(n+1)^2}$ と $\frac1{24n^2}$ に収束する差の和になり、主張を得る。$\square$
$n=10$ では $H_{10}=\frac{7381}{2520}=2.9289682\ldots$、$\log10.5=2.3513752\ldots$ なので $d_{10}=0.5775930\ldots$ である。thm-euler-constant-half-shift から
$$
0.5771763\ldots=d_{10}-\frac1{2400}<\gamma< d_{10}-\frac1{2904}=0.5772486\ldots
$$
となり、10 項の和だけで $\gamma$ の小数第 3 位まで決まる。同じ $n=10$ で thm-euler-constant-rate を使うと $0.5763\ldots<\gamma<0.5809\ldots$ にとどまる。実際の誤差を並べると次のとおりである(直接の計算による)。
| $n$ | $H_n-\log n-\gamma$ | $\frac1{2n}$ | $H_n-\log(n+\frac12)-\gamma$ | $\frac1{24n^2}$ |
|---|---|---|---|---|
| $1$ | $0.42278$ | $0.5$ | $0.017319$ | $0.041667$ |
| $10$ | $0.049167$ | $0.05$ | $0.00037733$ | $0.00041667$ |
| $100$ | $0.0049917$ | $0.005$ | $0.0000041252$ | $0.0000041667$ |
以下、実数 $t$ に対し $\lfloor t\rfloor$ を $t$ 以下の最大の整数、$\{t\}:=t-\lfloor t\rfloor$ を小数部分とする。
次の広義積分は収束し、
$$
\gamma=\int_1^{\infty}\Bigl(\frac1{\lfloor t\rfloor}-\frac1t\Bigr)dt=1-\int_1^{\infty}\frac{\{t\}}{t^2}\,dt
$$
が成り立つ。
どちらの被積分関数も $[1,\infty)$ で $0$ 以上なので、上端 $X$ までの積分は $X$ について単調増加であり、$X$ を整数 $N$ に限った極限が有限であれば、広義積分はその値に収束する(広義積分 の命題「非負の関数の広義積分」)。
$[k,k+1)$ では $\lfloor t\rfloor=k$ なので
$$
\int_1^N\Bigl(\frac1{\lfloor t\rfloor}-\frac1t\Bigr)dt=\sum_{k=1}^{N-1}\frac1k-\log N=H_N-\log N-\frac1N
$$
であり、これは $N\to\infty$ で $\gamma$ に収束する。
また $[k,k+1)$ では $\{t\}=t-k$ なので
$$
\int_1^N\frac{\{t\}}{t^2}\,dt=\log N-\sum_{k=1}^{N-1}k\Bigl(\frac1k-\frac1{k+1}\Bigr)=\log N-\sum_{k=1}^{N-1}\frac1{k+1}=\log N-H_N+1
$$
である。よって $1-\int_1^N\frac{\{t\}}{t^2}\,dt=H_N-\log N\to\gamma$ である。$\square$
2 つ目の式から、$1-\gamma=0.4227843\ldots$ は関数 $\{t\}/t^2$ のグラフの下の面積である。同じ形の積分は、和 $\sum_{n\leq x}f(n)$ を積分 $\int f$ と小数部分の積分に分ける Abel の総和公式から自然に現れる(Abelの級数変形法とその応用)。
次の級数表示では、整数 $k\geq2$ について $\zeta(k)=\sum_{n\geq1}n^{-k}$ とする(Riemannゼータ関数)。
$$
\gamma=1-\log2+\sum_{k=2}^{\infty}\frac{(-1)^k\bigl(\zeta(k)-1\bigr)}{k}
\qquad\text{および}\qquad
\gamma=\sum_{k=2}^{\infty}\frac{(-1)^k\zeta(k)}{k}
$$
が成り立つ。左の級数は絶対収束し、右の級数は収束するが絶対収束はしない。
$H_N-\log(N+1)=(H_N-\log N)-\log(1+\frac1N)\to\gamma$ であり、$\log(N+1)=\sum_{n=1}^{N}\log(1+\frac1n)$ なので
$$
\gamma=\sum_{n=1}^{\infty}\Bigl(\frac1n-\log\Bigl(1+\frac1n\Bigr)\Bigr)
$$
である。$n=1$ の項は $1-\log2$ である。$n\geq2$ では $\frac1n\leq\frac12$ なので、$\log(1+x)=\sum_{k\geq1}\frac{(-1)^{k+1}x^k}{k}$($|x|<1$)を $x=\frac1n$ で使って
$$
\frac1n-\log\Bigl(1+\frac1n\Bigr)=\sum_{k=2}^{\infty}\frac{(-1)^k}{k\,n^k}
$$
を得る。ここで
$$
\sum_{n=2}^{\infty}\sum_{k=2}^{\infty}\frac1{k\,n^k}\leq\sum_{n=2}^{\infty}\sum_{k=2}^{\infty}\frac1{n^k}=\sum_{n=2}^{\infty}\frac1{n(n-1)}=1
$$
なので、二重級数 $\sum_{n\geq2}\sum_{k\geq2}\frac{(-1)^k}{kn^k}$ は絶対収束し、和の順序を入れ替えてよい(Fubiniの定理 の例「二重級数の順序交換」)。$\sum_{n\geq2}n^{-k}=\zeta(k)-1$ なので
$$
\sum_{n=2}^{\infty}\Bigl(\frac1n-\log\Bigl(1+\frac1n\Bigr)\Bigr)=\sum_{k=2}^{\infty}\frac{(-1)^k\bigl(\zeta(k)-1\bigr)}{k}
$$
であり、左の式を得る。同じ評価から $\sum_k\frac{\zeta(k)-1}{k}\leq1$ なので、この級数は絶対収束する。
右の式:調和級数 の記事の系「交代調和級数の和」により $\sum_{k\geq1}\frac{(-1)^{k+1}}k=\log2$ なので、$\sum_{k\geq2}\frac{(-1)^k}{k}=-(\log2-1)=1-\log2$ である。収束する 2 つの級数を項ごとに足すと
$$
\sum_{k=2}^{\infty}\frac{(-1)^k\zeta(k)}{k}=\sum_{k=2}^{\infty}\frac{(-1)^k}{k}+\sum_{k=2}^{\infty}\frac{(-1)^k\bigl(\zeta(k)-1\bigr)}{k}=(1-\log2)+(\gamma-1+\log2)=\gamma
$$
である。$\zeta(k)\geq1$ なので $\sum_k\frac{\zeta(k)}{k}\geq\sum_k\frac1k=\infty$ であり、右の級数は絶対収束しない。$\square$
$0<\zeta(k)-1\leq\sum_{n\geq2}n^{-k}$ はおよそ $2^{-k}$ の大きさなので、左の級数は速く収束する。$k=20$ までの和で $0.57721568\ldots$ となり、誤差は $2\times10^{-8}$ 程度である(直接の計算による)。右の級数は項の大きさが $\frac1k$ 程度で、収束は遅い。
Euler の定数は、$\Gamma$ 関数と Riemann ゼータ関数の展開に定数項として現れる。以下はそれぞれの記事で証明されている。
Euler の定数は、整数の平均的な性質を測る式の 2 番目の項として現れる。
$H_n-\log n=\sum_{k=1}^{n}f(k)-\int_1^nf(t)\,dt$($f(t)=\frac1t$)の収束は、$f$ が減少することと $0$ 以上であることから従う(和と積分の差 の定理「和と積分の差の収束」)。また thm-euler-constant-zeta-series の右の級数は絶対収束しない。それぞれの条件を 1 つだけ外すと、次のように結論が成り立たなくなる。
| 外す条件 | 反例 | 成り立たなくなること |
|---|---|---|
| $f$ が減少する | $f(t)=\sin^2(\pi t)$($0$ 以上だが減少しない) | $f(k)=0$ なので和は $0$、積分は $\frac{n-1}2$ で、差は $-\frac{n-1}2\to-\infty$ と発散する |
| $f$ が $0$ 以上 | $f(t)=-t$(減少するが負) | 差は $-\frac{n(n+1)}2+\frac{n^2-1}2=-\frac{n+1}2$ で発散する |
| 級数の絶対収束 | $\sum_{k\geq2}(-1)^k\zeta(k)/k$ の項を、正の項 2 つと負の項 1 つの順に並べ替える | 並べ替えた級数の和は $\gamma$ と一致しない(条件収束する級数は並べ替えで和が変わる) |
最後の行について:正の項は $k$ が偶数の項 $\frac{\zeta(k)}k$、負の項は $k$ が奇数の項 $-\frac{\zeta(k)}k$ である。$\frac{\zeta(k)}{k}=\frac1k+\frac{\zeta(k)-1}{k}$ と分けると、後半は絶対収束するので並べ替えても和は変わらない。前半は交代調和級数から $1$ を除いたものの並べ替えであり、正の項 2 つと負の項 1 つの順に並べると和が $\frac12\log2$ だけ増える。実際、正の項 $\frac1{2j}$ を 2 つ、負の項 $-\frac1{2j+1}$ を 1 つずつ $m$ 組並べた部分和は
$$
\frac12H_{2m}-\Bigl(H_{2m+1}-\frac12H_m-1\Bigr)\to1-\frac12\log2\qquad(m\to\infty)
$$
であり(thm-euler-constant-rate の $H_n=\log n+\gamma+O(\frac1n)$ を代入すると $\log$ と $\gamma$ の項が $\frac12\log(2m)-\log(2m)+\frac12\log m=-\frac12\log2$ に整理される)、並べ替える前の和 $1-\log2$ より $\frac12\log2$ 大きい。
Euler は 1734 年に発表し 1740 年に刊行された論文でこの定数を導入した(Eul40)。記号 $\gamma$ はほかの定数にも使われるので、文献を読むときは式の形 $\lim(H_n-\log n)$ で確かめるのがよい。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する