🔢 mathematics

On the series expansion of the prime zeta function about $s=1$ and its coefficients

この論文は、リーマンゼータ関数のスティルチェス係数に類似する一般公式を導出し、メルテンスの定理を高次へ一般化するとともに、数値的に検証した素数ゼータ関数の $s=1$ における対数特異点まわりの級数展開を扱っています。

原著者： Artur Kawalec

公開日 2026-03-24

📖 1 分で読めます🧠 じっくり読む

原著者： Artur Kawalec

原論文は CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) でライセンスされています。 ✨ これは以下の論文のAI生成解説です。著者が執筆または承認したものではありません。技術的な正確性については原論文を参照してください。免責事項の全文を読む

この論文は、数学の奥深くにある「素数（2, 3, 5, 7...）」の謎を解き明かすための、新しい「地図の描き方」を提案するものです。

専門用語を避け、日常の風景や料理に例えて、この研究が何をしているのかを説明します。

1. 舞台：素数の「お宝地図」

まず、素数（2, 3, 5, 7...）は、数学という巨大な森に散らばっている「お宝」だと想像してください。
数学者たちは、これらの素数がどのように分布しているかを理解するために、**「素数ゼータ関数（P(s)）」**という特別なコンパスを使います。

このコンパスは、ある特定の場所（ $s=1$ という点）で非常に乱暴な動きをします。まるで、コンパスの針がその地点で激しく揺れ動き、針が折れてしまうかのように、**「特異点（カオスな場所）」**になっているのです。

2. 問題：カオスな場所をどう扱うか？

この論文の著者（アルトゥル・カワレック氏）は、このカオスな場所（ $s=1$ ）を「無視」するのではなく、**「そのカオスを分解して、規則性を見つけ出す」**ことに成功しました。

従来の考え方： 「ここはカオスだから、近づくと計算が破綻する」として避けていた。
この論文のアプローチ： 「このカオスには、実は『対数（log）』という特定の形をした『嵐』と、その後に続く『静かな余韻（係数）』がある」と見抜いた。

3. 解決策：嵐を「分解」する

著者は、このカオスな場所を以下のように 2 つに分けて考えました。

嵐そのもの（対数項）：
$s=1$ に近づくと、値が無限大に跳ね上がる部分。これは「 $\log(\frac{1}{s-1})$ "という形をしており、予測可能な「嵐の強さ」です。
嵐の後に残る「静かな余韻」（係数 $\alpha_n$ ）：
嵐の強さを差し引いた後に残る、非常に細かい数字の並び。これがこの論文の真骨頂です。

著者は、この「静かな余韻」を構成する数字（ $\alpha_0, \alpha_1, \alpha_2...$ ）を、「メルトンの定理」という古いレシピを、より高度なレベルに拡張した新しい公式を使って計算できることを証明しました。

4. 具体的な発見：新しい「味付け」のレシピ

この「静かな余韻」の数字たちは、これまで誰も正確に計算したことがないものばかりでした。

$\alpha_0$ （最初の余韻）：
素数の逆数を足し合わせた値から、ある定数（オイラー定数）を引いたもの。これは「メッセル＝メルトンの定数」として知られていますが、著者はこれを新しい視点で再定義しました。
$\alpha_1, \alpha_2...$ （その後の余韻）：
ここが今回の新発見です。これらは、素数の対数（ $\log p$ ）を何乗かして足し合わせた値から、ある特定の形を引いた極限として表せます。

著者は、これらを**「リーマン・ゼータ関数（素数ゼータの親分のような存在）」**の性質を使って、非常に高速に計算できる新しい公式も導き出しました。
まるで、手作業で 1 粒ずつ数えていたお米を、最新の機械で一瞬で計量できるようになったようなものです。

5. 結果：高精度な「数値のリスト」

論文の最後には、この新しい公式を使って計算した、**「最初の 10 個の余韻の数字」**が、30 桁もの高精度でリストアップされています。

これらは、単なる数字の羅列ではなく、素数の分布という「森の地形」をより深く理解するための、新しい**「等高線」**のようなものです。

まとめ：この論文は何をしたのか？

一言で言えば、**「素数が集まる場所の『カオス』を、嵐と余韻に分け、その『余韻』の正体を解き明かす新しい計算式を発見した」**という研究です。

昔：「ここはカオスでわからない」と言っていた。
今：「カオスの正体は『嵐』と『余韻』だった。そして『余韻』の値は、この新しいレシピ（公式）を使えば、超高精度で計算できる！」

これにより、数学者たちは素数の分布について、これまで以上に詳細で精密な地図を描くことができるようになりました。まるで、荒れ狂う海を航海する際に、波の揺れ（対数項）を予測し、その後に残る静かな波の動き（係数）まで正確に把握できるようになったようなものです。

Artur Kawalec による論文「On the series expansion of the prime zeta function about s = 1 and its coefficients（素数ゼータ関数の s=1 における級数展開とその係数）」の技術的サマリーを以下に示します。

1. 研究の背景と問題設定

対象関数: 素数ゼータ関数 $P(s) = \sum_{p} p^{-s}$ （ $p$ は素数）。この関数は $\Re(s) > 1$ で絶対収束し、メビウス反転公式を用いることで $\Re(s) > 0$ へ解析接続可能ですが、 $s=1/n$ やゼータ関数の非自明な零点 $\rho_n$ において対数特異点（分岐点）を持ち、多価関数となります。
問題点: $s=1$ における $P(s)$ の振る舞いは、メルテンスの定理（Mertens's Theorems）などの既知の結果と関連していますが、 $s=1$ における対数特異点を取り除いた後の正則部分の級数展開係数（ $\alpha_n$ ）に対する一般的な公式や、その高次項の精密な計算値は十分に確立されていませんでした。
目的: $s=1$ における $P(s)$ の対数特異点を取り除いた形式で級数展開を行い、その展開係数 $\alpha_n$ に対する一般公式を導出すること、およびその係数を高精度で計算すること。

2. 主要な手法とアプローチ

論文では、以下の数学的アプローチを組み合わせています。

級数展開の定義:
$P(s)$ を $s=1$ 周りで以下のように展開します。
$P(s) = \log\left(\frac{1}{s-1}\right) + \sum_{n=0}^{\infty} \frac{\alpha_n}{n!} (s-1)^n$
ここで、 $\log\left(\frac{1}{s-1}\right)$ が特異部分を表し、残りの級数が正則部分となります。
メビウス反転公式の利用:
$P(s) = \sum_{n=1}^{\infty} \frac{\mu(n)}{n} \log \zeta(ns)$ という関係式を用い、リーマンゼータ関数 $\zeta(s)$ の展開と結びつけることで、収束性の良い表現を導きます。
スティルチェス積分（Stieltjes Integration）と素数分布:
素数分布関数 $\pi(x)$ と対数積分関数 $\text{li}(x)$ の差 $f(x) = \pi(x) - \text{li}(x)$ を用いた積分表現を構築します。リーマン予想（RH）を仮定して $f(x) = O(\sqrt{x}\log x)$ と評価し、積分の収束性を保証します。
係数の導出:
- 定理 1: 素数の和に関する極限式として $\alpha_n$ を定義します。
- 定理 2: メビウス関数とゼータ関数の導関数を用いた、より高速に収束する級数表現を導出します。
- 定理 3: 積分形式による代替公式を提示します。

3. 主要な成果と結果

定理 1: 極限による係数の定義

展開係数 $\alpha_n$ は、素数 $p$ に関する和の極限として以下のように表されます。

$n=0$ の場合（メルテンスの定数に関連）:
$\alpha_0 = \lim_{x\to\infty} \left( \sum_{p\le x} \frac{1}{p} - \log\log x \right) - \gamma$
ここで $\gamma$ はオイラー定数です。
$n \ge 1$ の場合:
$\alpha_n = \lim_{x\to\infty} (-1)^n \left( \sum_{p\le x} \frac{\log^n p}{p} - \frac{\log^n x}{n} \right)$
この極限式はスチルチェス定数（Stieltjes constants）の公式と類似していますが、数値計算においては非常に遅く収束します。

定理 2: 高速収束する級数表現（主要な貢献）

定理 1 の極限式は計算が困難であるため、メビウス反転公式を用いた以下の表現を導出しました。これにより、 $\alpha_n$ をゼータ関数の導関数とメビウス関数の級数として計算できます。
$\alpha_n = g_n + \frac{d^n}{ds^n} \left[ \sum_{k=2}^{\infty} \frac{\mu(k)}{k} \log \zeta(ks) \right]_{s\to 1}$
ここで $g_n$ は $\log \zeta(s)$ の $s=1$ における展開係数であり、スチルチェス定数 $\gamma_k$ を用いて表現されます。

この表現は、極限式に比べてはるかに高速に収束するため、高精度な数値計算を可能にします。
具体的には、 $g_n$ がスチルチェス定数の多項式で表されることを示し、 $\alpha_n$ を同様に計算可能な形にしました。

数値計算結果

論文では、上記の定理 2 の公式を用いて、最初の 10 個の係数 $\alpha_0$ から $\alpha_{10}$ までを 30 桁以上の精度で計算しました（Table 1 参照）。

例: $\alpha_0 \approx -0.315718\dots$ （OEIS A143524）、 $\alpha_1 \approx 1.332582\dots$ など。
これらの値は、既存の文献（メルテンスの定理など）で知られている低次項と一致し、高次項については初めて高精度で計算された値として提供されています。

定理 3: 積分表現

係数 $\alpha_n$ に対して、素数分布の誤差項 $f(t) = \pi(t) - \text{li}(t)$ を用いた積分公式も導出されました。
$\alpha_n = (-1)^n \int_{1}^{\infty} \frac{\log^n t - n \log^{n-1} t}{t^2} [\pi(t) - \text{li}(t)] dt$

4. 意義と貢献

メルテンスの定理の一般化:
従来のメルテンスの定理（ $n=0, 1$ のケース）を、任意の次数 $n$ に対する高次の一般化として定式化しました。
計算手法の革新:
素数ゼータ関数の展開係数を計算する際、従来の遅い収束する極限式ではなく、メビウス関数とゼータ関数の導関数を用いた高速収束する級数表現を確立しました。これにより、高次係数の高精度計算が現実的になりました。
解析的性質の解明:
$s=1$ における $P(s)$ の特異点構造を明確にし、正則部分の係数がゼータ関数の性質（特にスチルチェス定数）と密接に関連していることを示しました。
数値データベースの提供:
高精度で計算された係数 $\alpha_n$ の値を提供し、今後の数論研究やアルゴリズム検証のための基準値（ベンチマーク）として機能します。

結論

本論文は、素数ゼータ関数の $s=1$ における局所的な振る舞いを詳細に解析し、その展開係数に対する理論的な一般公式と、実用的な高精度計算手法を確立した点で重要な貢献をしています。特に、スチルチェス定数との関係性を明確にしつつ、数値計算の効率を劇的に改善した手法は、数論における計算研究の進展に寄与するものです。