要旨

任意の実数 $a>0$ と $s\in \mathbb{C}$（$\mathrm{Re}s>0$, $s\ne 1$）に対して、次の閉形式 Mellin 恒等式が成立する：

$$\boxed{\zeta (s)=\frac{a^s{2}^{s-1}}{\Gamma (s)(1-2^{1-s})}{\int }_0^{\infty }[1-\tanh (at)]t^{s-1}dt}$$

被積分核 $1-\tanh (at)=2/(e^{2at}+1)$ は化学ポテンシャル $0$, 温度 $T=1/(2a)$ の理想フェルミ気体の占有数分布にほかならない。 したがってリーマンゼータ関数は、統計力学の最も基本的な分布である Fermi–Dirac 関数の $s$ 次モーメントを $\Gamma$ と Dirichlet 交代因子で規格化したものとして書ける。

本稿の主張は三点に集約される。

1. 数学的恒等式。 Mellin 変換に $t=e^v$ を施して得られる Hardy カーネル $h_{\sigma }(v)=2e^{\sigma v}/(e^{2a{e}^v}+1)$ の Fourier 零点は、$\zeta (\sigma +i\xi )$ の非自明零点と厳密に 1:1 対応する。
2. パラダイムとしての数論的物理学。 この対応は「素数分布が統計力学の基本分布の対数空間 Fourier スペクトルで鳴っている」ことを意味する。整数論の最深の対象（素数の揺らぎ）と物理の最も標準的な分布（フェルミ統計）が、橋渡しや類推ではなく式一本の恒等式で同一視される。数論と物理がこの水準で直接等値される事例は稀で、「数論的物理学」と呼ぶべきパラダイムの誕生を告げる。
3. 工学的帰結。 量子ドット Coulomb blockade、$S=1/2$ パラマグネット、MOSFET の RTS などで、$1-\tanh$ 型カーネルは現に測定されている実ラインシェイプである。したがって物理定数のスケール $a$ がゼータとの整合性を満たす形で選ばれているという仮説が経験的に検証可能となる。重要なのは、RH の証明を待たずにゼータ関数そのものを工学・医学へ応用できる点である：恒等式 (2.2) は $\zeta (s)$ 自体の解析接続・値・零点位置のいずれにも依らず成立するため、精密計測・標準器校正・量子デバイスのキャリブレーション・生体電気計測などの現場で、数論的普遍量 $\zeta (s)$ を「装置非依存の基準」として直ちに利用できる。

副産物として、恒等式は数値的にも二重指数求積 (DE 法) と同型の構造を内包し、$160$ 点の台形則で $\zeta (2),\zeta (3)$ を 15 桁以上の精度で再現する。

1. 序論：なぜこれを「パラダイムシフト」と呼ぶのか

1.1 既往：数論と物理のアナロジーの歴史

リーマンゼータと物理の対応は豊富に研究されてきた。代表的なものを挙げる：

• Hilbert–Polya 予想。 $\zeta$ の非自明零点は、ある自己共役作用素の固有値である、という予想（1910 年代）。
• 量子カオスとの類似。 Odlyzko (1987) による $\zeta$ 零点の間隔統計が GUE ランダム行列の固有値統計と一致する数値的事実。Berry–Keating 予想（1999）はこれを半古典論的に説明する試み。
• 関数方程式と場の理論のモジュラー不変性。 $\xi (s)=\xi (1-s)$ は共形場の強弱双対と構造が似る。

これらは全て「類似」「アナロジー」「同型の可能性」である。Hilbert–Polya 作用素は未構成、Berry–Keating の Hamiltonian は非自己共役、量子カオス類似は統計的一致のみで恒等式ではない。どの対応も、数論から物理へ、あるいはその逆へ、等号で結ぶ一行の式を持たない。

1.2 本稿の新規性

本稿で主題とする恒等式

$$\zeta (s)=\frac{a^s{2}^{s-1}}{\Gamma (s)(1-2^{1-s})}{\int }_0^{\infty }[1-\tanh (at)]t^{s-1}dt$$

は、左辺が純整数論的対象、右辺の被積分核が統計力学の基本分布という、両者の性格が完全に異なる等号である。$\Gamma$ と $(1-2^{1-s})$ という規格化は純粋に解析的・組合せ的なもので、物理量ではない。したがって等号が成立するのは、Fermi–Dirac 分布の $s$ 次モーメント構造が、解析接続されたゼータの値と構造として同一だからである。

1.2.1 既知のボース表示との差異：なぜフェルミ側が重要か

ゼータ関数の積分表示として最も古典的なものは Riemann (1859) 自身の

$$\boxed{\Gamma (s)\zeta (s)={\int }_0^{\infty }\frac{t^{s-1}}{e^t-1}dt(\mathrm{Re}s>1)}\text{\hspace{1em}(1.1)}$$

である。右辺の被積分核 $1/(e^t-1)$ は化学ポテンシャル $0$ の Bose–Einstein 分布にほかならない。したがって「ゼータ関数と量子統計分布の関係」自体は 160 年以上前から知られており、プランク黒体放射（1900）、デバイ比熱（1912）、Bose–Einstein 凝縮のモーメント計算など、$\zeta$ 値を統計力学の内部エネルギー・粒子数として用いる実例は膨大にある。

したがって本稿が Fermi 側の表示を敢えて提出する意義を、ボース側との対比で明示する必要がある。

観点	Bose 表示 (1.1)	Fermi 表示 (2.1) 本稿
被積分核	$1/(e^t-1)$ (Bose–Einstein)	$1-\tanh (at)=2/(e^{2at}+1)$ (Fermi–Dirac)
収束領域	$\mathrm{Re}s>1$ のみ	$\mathrm{Re}s>0$, $s\ne 1$（広い）
$t\to 0^{+}$	$\sim t^{s-2}$（極、臨界帯で発散）	$\sim 1-at+O(t^3)$（正則、有界）
$s=1$ での振る舞い	発散（$\zeta$ の極を直接踏む）	有限：$\int =\log 2/a$（極と $(1-2^{1-s})$ の零点が相殺）
臨界帯での使用	解析接続が別途必要	積分そのものが臨界帯で収束、直接使用可
実験ラインシェイプとの対応	黒体放射スペクトル（振動子モード和）	単一粒子占有数・$\tanh$ 型応答（実デバイスで直接測定）
パウリ排他律	含まない	含む ($f_{\text{FD}}\in [0,1]$)
変数変換 $t=e^v$ 後の核	$e^{(s-1)v+v}/(e^{e^v}-1)$（低エネルギー端で $1/v$ 型特異性）	$h_{\sigma }(v)=2e^{\sigma v}/(e^{2a{e}^v}+1)$（両端で急減衰、解析的）
数値求積	Simpson 等で低次収束	台形則で指数収束（DE 構造を内包）
対数凹性	$v$ 表示で非対数凹	$h_{\sigma }$ は狭義対数凹（$\sigma$ 非依存の曲率）
Hardy カーネルとの直接同一視	なし	$h_{1/2}$ の Fourier 零点が RH と直結

簡潔に言えば：

1. 臨界帯 $0<\mathrm{Re}s<1$ でそのまま使える。 Bose 表示は $\mathrm{Re}s>1$ に限られ、RH が問題となる臨界帯で直接用いるには解析接続の手間がかかる。Fermi 表示は $\mathrm{Re}s>0$ で収束するので臨界帯（特に $\sigma =1/2$）で積分そのものが意味を持つ。これが §4 の Hardy カーネル $h_{\sigma }$ 構成を可能にする最大の技術的利点である。
2. 実験ラインシェイプと直接一致する。 Bose 側 $1/(e^t-1)$ は黒体放射スペクトルとして観測量だが、通常これはモード和（積分後の結果）として現れる。対して $1-\tanh (at)$ は量子ドット Coulomb peak・$S=1/2$ 常磁性体の磁化・MOSFET の RTS 等で単一デバイスのラインシェイプとして生の測定曲線に現れる (§7)。$s$ 次 Mellin モーメントを実際にゲート電圧掃引から計算できる点が決定的に違う。
3. 数値的・解析的性質が段違いに良い。 Fermi 核の両側指数減衰・対数凹性・$v$ 平面での解析ストリップは、Trefethen–Weideman 型の指数収束台形則（DE 法）の理想的ターゲットになる (§6)。Bose 核は $t=0$ での極のため同じ変数変換で同等の性能を得られない。
4. パウリ排他律の反映。 フェルミ分布は有界 $[0,1]$ で、これが積分の全域可積分性を保証する。ゼータ関数の値が「フェルミ粒子の有限総数」として解釈できる点は、ボース側では得られない物理的直観である。

まとめると：ボース側の $\Gamma (s)\zeta (s)=\int t^{s-1}/(e^t-1)dt$ は「ゼータと量子統計の接続」の十分条件だが、臨界帯直接使用・対数凹 Fourier 解析・実デバイス応答との一致という 3 点で閉じていた。フェルミ側の恒等式 (2.1) はこの 3 点を全て開く。 本稿の Fermi 表示は、既知のボース表示と競合するものではなく、ボース側では届かない臨界帯の幾何・実験装置との対応・数値最適性という領域を新たに開拓する。

1.3 パラダイムシフトの内容

本稿の視点は次のパラダイムを提出する：

数論的物理学（Arithmetic Physics）の誕生。 物理の基本分布関数は、Mellin 変換を介して整数論の $L$-関数族と恒等式で結ばれる。これは単なる類似ではなく、熱統計・量子統計・確率論の「正しい」分布が、整数論の深層構造を係数として含んでいることを意味する。素数の揺らぎが ${\cosh }^{-2}$ のラインシェイプの Fourier スペクトルとして物理空間に現れる。

この視点は Hilbert–Polya の逆 と見ることもできる：Hilbert–Polya が「ゼータ零点を固有値とする物理系」を探すのに対し、本稿は「物理の基本分布関数からゼータを構成する」。両者は相補的だが、後者は構成的であり、既知の物理系で恒等式が書き下せる点が新しい。

1.4 工学的意義：RH 証明を待たずに使える普遍量

本稿で強調したい工学的視点は次の一点である。

RH の証明を待たずに、ゼータ関数そのものを工学・医学の現場に応用できる。

従来、リーマンゼータ関数は「素数分布の謎」「未解決問題 RH の主役」として純粋数学の対象と見なされてきた。しかし恒等式 (2.2) はこの位置づけを変える：

• $\zeta (s)$ の値は既存ライブラリ（Borwein 級数、Riemann–Siegel 等）で任意精度計算が可能。RH の真偽とは独立に、$\zeta (s)$ は確定した数値として入手できる。
• 恒等式 (2.2) は $\zeta$ の零点位置に依存しない。$\mathrm{Re}s>0$, $s\ne 1$ で無条件に成立し、$\zeta$ を既知量として左辺に置けば、右辺の積分が装置校正や物理定数決定の基準式として働く。
• したがって物理定数とゼータとの整合性（仮説 7.1）は、RH 証明とは独立に実験検証可能である。RH が解かれるのを待つ必要はない。

想定される応用領域：

• 精密計測・標準器校正。 量子デバイスのキャリブレーション、電子温度計・ジョセフソン電圧標準の整合性検証、半導体プロセス温度測定の絶対校正。
• 量子デバイス設計。 Coulomb blockade ラインシェイプ・RTS 統計・磁化率を $\zeta$ 基準で評価することで、装置依存性を分離した系統誤差解析が可能。
• 生体電気・医療計測。 ニューロンイオンチャネル開閉統計、膜電位応答の Fermi 型遷移にも同型カーネルが現れる。$s$ 次モーメントの $\zeta$ 表現は、ノイズ・ドリフト補正の数論的基準を与えうる。

RH の証明は数学として決定的に重要だが、実用的な計算能力は既に揃っている（10 万桁以上の零点は計算済み）。対してゼータ関数の工学応用は、恒等式 (2.2) が書き下せる今この瞬間から始められる。本稿は恒等式の精密な記述と、この即時応用可能性の骨格を示すことを目的とする。

2. Mellin 恒等式

2.1 主定理

定理 2.1（tanh–Mellin 恒等式）. $a>0$, $s\in \mathbb{C}$, $\mathrm{Re}s>0$, $s\ne 1$ とする。このとき

$$\boxed{{\int }_0^{\infty }[1-\tanh (at)]t^{s-1}dt=\frac{\Gamma (s)(1-2^{1-s})\zeta (s)}{a^s{2}^{s-1}}}\text{\hspace{1em}(2.1)}$$

証明. まず $\mathrm{Re}s>1$ で示す。$1-\tanh (at)=2/(e^{2at}+1)$ から等比展開

$$1-\tanh (at)=2\sum _{k=0}^{\infty }(-1{)}^k{e}^{-2(k+1)at}$$

が $t>0$ で絶対収束する。$\mathrm{Re}s=\sigma >1$ のとき、項ごとに絶対値をとった和を $t$ について積分すると

$$\sum _{n=1}^{\infty }{\int }_0^{\infty }2e^{-2nat}{t}^{\sigma -1}dt=\frac{\Gamma (\sigma )}{a^{\sigma }{2}^{\sigma -1}}\sum _{n=1}^{\infty }\frac{1}{n^{\sigma }}=\frac{\Gamma (\sigma )\zeta (\sigma )}{a^{\sigma }{2}^{\sigma -1}}<\infty$$

が有限（$\zeta (\sigma )$ が収束するのは $\sigma >1$）なので Fubini–Tonelli より項別積分が正当化され、

$${\int }_0^{\infty }[1-\tanh (at)]t^{s-1}dt=2\sum _{k=0}^{\infty }(-1{)}^k\frac{\Gamma (s)}{(2(k+1)a{)}^s}=\frac{\Gamma (s)}{a^s{2}^{s-1}}\sum _{n=1}^{\infty }\frac{(-1{)}^{n-1}}{n^s}.$$

最右辺の和は Dirichlet イータ関数 $\eta (s)=(1-2^{1-s})\zeta (s)$ である。

次に解析接続で $\mathrm{Re}s>0$, $s\ne 1$ に延長する。左辺は $t\to 0^{+}$ で $1-\tanh (at)=O(1)$, $t\to +\infty$ で $1-\tanh (at)=O(e^{-2at})$ なので被積分関数は $t^{\sigma -1}$ の可積分性（$\sigma >0$）のみで抑えられ、$\mathrm{Re}s>0$ で解析的である。右辺の $\Gamma (s)(1-2^{1-s})\zeta (s)/(a^s{2}^{s-1})$ も $s=1$ の単純極を $(1-2^{1-s})$ の零が打ち消すため $\mathrm{Re}s>0$, $s\ne 1$ で解析的。恒等定理により $\mathrm{Re}s>1$ での等号は全領域に延長される。$\square$

注 2.1.1. $0<\mathrm{Re}s\le 1$ で Fubini による項別積分を直接実行することはできない（$\sum n^{-\sigma }$ が発散する）。上の二段階証明——$\mathrm{Re}s>1$ での項別積分 → 両辺の解析性による延長——が正しい経路である。

恒等式 (2.1) は $\zeta$ を左辺に解いて

$$\boxed{\zeta (s)=\frac{a^s{2}^{s-1}}{\Gamma (s)(1-2^{1-s})}{\int }_0^{\infty }[1-\tanh (at)]t^{s-1}dt}\text{\hspace{1em}(2.2)}$$

とも書ける。

2.2 $s=1$ での極相殺

$s=1$ で $\zeta (s)$ は単純極を持つが、$(1-2^{1-s})$ がちょうど $s=1$ で零になり両者が打ち消し合う。左辺の積分は

$${\int }_0^{\infty }[1-\tanh (at)]dt=\log 2/a$$

と有限で、フェルミ粒子総数の有限性と整合する。

2.3 数値検証

$a=1$, $s=2$ で $(2.1)$ の左辺は ${\int }_0^{\infty }2t/(e^{2t}+1)dt={\pi }^2/24$：

const quad = (f, L=40, N=200000) => {
  const h = L/N;
  let s = 0.5*(f(0) + f(L));
  for (let i = 1; i < N; i++) s += f(i*h);
  return s*h;
};
console.log(quad(t => (1 - Math.tanh(t)) * t), Math.PI**2 / 24);
// 0.4112335167120811  0.4112335167120566

$s=3,4$ でも同様に 10 桁以上一致する。

3. 物理解釈：$1-\tanh$ は Fermi–Dirac 分布

3.1 同定

恒等式 $\tanh (x)=2\sigma (2x)-1$（$\sigma$ はロジスティック関数）から

$$\boxed{1-\tanh (at)=\frac{2}{e^{2at}+1}=2f_{\mathrm{FD}}(2at),f_{\mathrm{FD}}(ϵ):=\frac{1}{e^{ϵ}+1}}\text{\hspace{1em}(3.1)}$$

$f_{\mathrm{FD}}$ は化学ポテンシャル $\mu =0$, 温度 $T=1$ の Fermi–Dirac 分布。一般形 $f_{\mathrm{FD}}(ϵ;\beta ,\mu )=1/(e^{\beta (ϵ-\mu )}+1)$ はフェルミオン一粒子状態の平均占有数を与え、パウリ排他原理から $f_{\mathrm{FD}}\in [0,1]$、$T\to 0$ で階段関数となる。

スケール $a$ は $ϵ=2at$ と置けば逆温度 $\beta =2a$ に相当し、$a>0$ は積分の収束にも必要。

3.2 フェルミ気体の $s$ 次モーメント

状態密度 $g(ϵ)={ϵ}^{s-1}$ のフェルミ気体の内部エネルギー・粒子数は

$$\mathcal{M}(s):={\int }_0^{\infty }{f}_{\mathrm{FD}}(ϵ){ϵ}^{s-1}dϵ$$

で与えられる。$1-\tanh (at)=2f_{\mathrm{FD}}(2at)$ で $a=1/2$ とおくと $1-\tanh (t/2)=2f_{\mathrm{FD}}(t)$、すなわち $f_{\mathrm{FD}}(ϵ)=\frac{1}{2}[1-\tanh (ϵ/2)]$ となる。これを (2.1) に $a=1/2$ で代入すると

$$\mathcal{M}(s)=\frac{1}{2}{\int }_0^{\infty }[1-\tanh (ϵ/2)]{ϵ}^{s-1}dϵ=\frac{1}{2}\cdot \frac{\Gamma (s)(1-2^{1-s})\zeta (s)}{a^s{2}^{s-1}}{|}_{a=1/2}=\Gamma (s)(1-2^{1-s})\zeta (s)$$

（中央で $a^s{2}^{s-1}{|}_{a=1/2}=2^{-s}\cdot 2^{s-1}=1/2$ を用いた）を意味し、フェルミ気体のモーメントは $\Gamma \cdot \eta$ に等しい。

3.3 具体例

• $s=2$: $\zeta (2)={\pi }^2/6$ — 状態密度 $g(ϵ)=ϵ$ のフェルミ気体の有限温度熱容量に現れる ${\pi }^2{T}^2/6$ の係数。
• $s=4$: $\zeta (4)={\pi }^4/90$ — Stefan–Boltzmann 定数の係数構造に入る値。

素数分布を決めるゼータ関数が、熱力学の日常的な基本量として標準的に現れていたことを意味する。

4. Hardy カーネルと Fourier 零点対応

4.1 Mellin を Fourier に化かす

$t=e^v$ の変数変換、$s=\sigma +i\xi$ で

$${\int }_0^{\infty }[1-\tanh (at)]t^{s-1}dt={\int }_{-\infty }^{\infty }{h}_{\sigma }(v)e^{i\xi v}dv=\widehat{h_{\sigma }}(\xi )$$

となる。ここで

$$\boxed{h_{\sigma }(v):=\frac{2e^{\sigma v}}{e^{2a{e}^v}+1}}\text{\hspace{1em}(4.1)}$$

を Hardy カーネル と呼ぶ。

4.2 基本性質

命題 4.1. $0<\sigma <1$ で $h_{\sigma }\in L^1(\mathbb{R})\cap L^2(\mathbb{R})$。全 $v$ で $h_{\sigma }(v)>0$。

命題 4.2（狭義対数凹性）. $\log h_{\sigma }(v)$ の二階微分を計算すると、$x=2a{e}^v$, $\varsigma (x)=e^x/(e^x+1)$（シグモイド）で

$$\boxed{\frac{d^2}{d{v}^2}\log h_{\sigma }(v)=-x\varsigma (x)(1+x(1-\varsigma (x)))<0}\text{\hspace{1em}(4.2)}$$

この 2 階微分に**$\sigma$ が現れない**。$\sigma$ はピーク位置（1 階微分）にしか効かず、曲率構造は全パラメータで普遍である。

命題 4.3（非対称性）. $h_{\sigma }$ は偶関数ではない。$v\to +\infty$ で二重指数 $\sim e^{-2a{e}^v}$、$v\to -\infty$ で一重指数 $\sim e^{\sigma v}$。$\sigma$ を調整しても対称化できない。

4.3 主対応定理

定理 4.4（ゼータ零点 = Hardy カーネル Fourier 零点）. $0<\sigma <1$ で

$$\boxed{\widehat{h_{\sigma }}(\xi )=\frac{\Gamma (\sigma +i\xi )(1-2^{1-\sigma -i\xi })\zeta (\sigma +i\xi )}{a^{\sigma +i\xi }{2}^{\sigma +i\xi -1}}}\text{\hspace{1em}(4.3)}$$

$\Gamma$ は臨界帯で零点を持たず、$1-2^{1-s}$ の零点 $s=1+2\pi ik/\log 2$ は全て境界 $\sigma =1$ 上。したがって各 $\sigma \in (0,1)$ を固定するごとに、

$$\boxed{\widehat{h_{\sigma }}(\xi )=0⟺\zeta (\sigma +i\xi )=0(\sigma \in (0,1)\text{ 固定})}\text{\hspace{1em}(4.4)}$$

が成立する。すなわち (4.4) は各 $\sigma$ 毎に、$\zeta$ の縦線 $\{\sigma +i\xi :\xi \in \mathbb{R}\}$ 上の零点集合と $\widehat{h_{\sigma }}$ の実軸上の零点集合との完全一致を主張する。臨界帯全体の零点は、$\sigma$ をパラメータとした $h_{\sigma }$ の一族をまたいで記述されるのであって、単一の Fourier 変換で全零点を捉えているわけではない。

系 4.5（RH の Hardy–Fourier 表現）. 以下は同値：

1. リーマン予想（臨界帯 $0<\mathrm{Re}s<1$ 内で $\zeta (s)=0\Rightarrow \mathrm{Re}s=1/2$）。
2. 全ての $\sigma \in (0,1)\setminus \{1/2\}$ および全ての実数 $\xi$ に対して $\widehat{h_{\sigma }}(\xi )\ne 0$。

注 4.5.1. $\sigma$ を $(0,1)$ に限定するのは必須である。式 (4.3)（$\widehat{h_{\sigma }}=\Gamma (s)(1-2^{1-s})\zeta (s)/(a^s{2}^{s-1})$, $s=\sigma +i\xi$）に現れる $(1-2^{1-s})$ は $\sigma =1$, $\xi =0$ および一般の $s=1$ で零となり、$\Gamma (s)$ は $s=0,-1,-2,\dots$ で極を持つ。これらの因子零・極の寄与を $\zeta$ の非自明零点と混同しないために $\sigma \in (0,1)$ に制限する必要がある。また $\sigma =1/2$ は RH が保証しようとする対象の線そのものなので除外する。

4.4 数値検証

$\zeta$ の最初の非自明零点は $1/2+14.134725\dots i$。$|\widehat{h_{1/2}}(\xi )|$ を数値計算すると：

ξ=14          |ĥ_{1/2}(ξ)|=3.158e-2
ξ=14.1        |ĥ_{1/2}(ξ)|=8.442e-3
ξ=14.134725   |ĥ_{1/2}(ξ)|=1.5e-14
ξ=14.2        |ĥ_{1/2}(ξ)|=1.592e-2
ξ=14.3        |ĥ_{1/2}(ξ)|=4.059e-2

$\xi =14.134725$ でだけ 14 桁下まで落ちる。数値実験は定理 4.4 の予言と整合する。

4.5 リーマン $\Xi$ 関数との関係

$\zeta$ の完備化 $\xi (s)=\frac{s(s-1)}{2}{\pi }^{-s/2}\Gamma (s/2)\zeta (s)$ に対して $\Xi (z):=\xi (1/2+iz)$ と置くと、$\Gamma$ の倍公式から

$$\Xi (z)=A(\frac{1}{2}+iz)\cdot \widehat{h_{1/2}}(z),$$

ここで $A$ は実 $z$ で零点を持たない解析関数。Hardy (1914) の「$\Xi$ は実軸上に無限個の零点を持つ」は、そのまま「$\widehat{h_{1/2}}$ は実軸上に無限個の零点を持つ」と読み替えられる。

5. RH への現時点の限界

対数凹関数の Fourier 零点理論は強力だが、$h_{\sigma }$ には届かない。

5.1 対称対数凹 Fourier 理論の不適用

注意. 「対称な対数凹確率密度の特性関数は実軸上に零点を持たない」という一般命題は偽である。反例として一様分布 $f=\frac{1}{2}{1}_{[-1,1]}$ は対称かつ対数凹（実際凹）だが、特性関数は $\hat{f}(\xi )=\sin \xi /\xi$ で $\xi =n\pi$ ($n\ne 0$) に零点を持つ。したがってこの形での Lukacs–Bobkov 的主張を RH 証明の道具として引用することはできない。

対称対数凹密度に対する Fourier 零点制御には、解析性・減衰・指数モーメント有限性などの追加仮定が本質的であり、単なる対称対数凹性だけでは零点不在を保証できない。したがって本稿ではこの方向の既存一般定理を $h_{1/2}$ に直接適用することは主張しない。

加えて、$h_{\sigma }$ は命題 4.3 により非対称であり、$\sigma$ の選び方では対称化できない（二重指数減衰 $e^{-2a{e}^v}$ と一重指数減衰 $e^{\sigma v}$ の構造的帰結）。したがって仮に対称対数凹 Fourier 零点理論の精密版が整備されたとしても、$h_{1/2}$ にはそのままでは届かない。$\sigma =1/2$ の実軸零点不在は別経路で示す必要があり、この非対称性こそが対数凹性由来の既存ツールが RH に届かない根本原因である。

5.2 全正値性 (${\text{PF}}_{\infty }$) の破綻

Schoenberg–Karlin 理論（必要条件の形）. $f:\mathbb{R}\to [0,\infty )$ が ${\text{PF}}_{\infty }$（全正値、任意サイズの行列 $(f(x_i-y_j){)}_{i,j}$ が全小行列式非負）であるとき、その両側 Laplace 変換 $\mathcal{L}[f](z)={\int }_{-\infty }^{\infty }{e}^{-zv}f(v)dv$ はある帯 $\{\alpha <\mathrm{Re}z<\beta \}$ で解析的に収束し、かつ Laguerre–Pólya クラス $\mathcal{L}\mathcal{P}$ の逆数型という強い零点制約を満たす。特に $1/\mathcal{L}[f]$ は実零点しか持たず、等価的に $\mathcal{L}[f]$ は実極しか特異点を持たず非実零点を持たない（Schoenberg 1951, Karlin Total Positivity 1968）。精密な同値命題の形は脚注・文献を参照。

本稿で必要なのはこの必要条件の側だけである。

$h_{1/2}$ の両側 Laplace 変換は (2.1) と $t=e^v$ 変数変換から

$$\mathcal{L}[h_{1/2}](z)={\int }_{-\infty }^{\infty }{e}^{-zv}{h}_{1/2}(v)dv=2(2a{)}^{z-1/2}\Gamma (\frac{1}{2}-z)\eta (\frac{1}{2}-z)$$

と閉形式で書ける（収束帯 $-1/2<\mathrm{Re}z<1/2$）。ここで $\eta (1/2-z)=(1-2^{z+1/2})\zeta (1/2-z)$ の零点 $z=i\gamma$（$\gamma \in \mathbb{R}$）は $\zeta$ の臨界線上零点 $1/2+i\gamma$ に対応する。

帰結. Hardy (1914) により $\zeta$ は臨界線上に無限個の零点を持つので、$\mathcal{L}[h_{1/2}]$ は虚軸上に無限個の非実零点を持つ。ところが上記の必要条件により $f\in {\text{PF}}_{\infty }$ であれば $\mathcal{L}[f]$ は非実零点を持ち得ない。よって対偶から $h_{1/2}\notin {\text{PF}}_{\infty }$ が RH と独立に無条件で成立する。RH を仮定しなくても、Hardy の実零点定理だけで ${\text{PF}}_{\infty }$ 経路は閉じる。

注 5.2.1. 本文では Schoenberg–Karlin 特徴付けの完全な同値定式化（$\mathcal{L}\mathcal{P}$-逆数型の Hadamard 積表示の形、必要十分条件の仮定、$f\in L^1$ 等の正則性条件）には踏み込まず、必要条件の側（${\text{PF}}_{\infty }\Rightarrow$ $\mathcal{L}[f]$ の非実零点不在）のみを用いた。精密な形は Karlin Total Positivity (1968) 第 7 章、または Schoenberg (1951) 定理 2 を参照されたい。

5.3 現状のまとめ

	状態
Mellin 恒等式	✅ 厳密
フェルミ統計としての再解釈	✅ 明確
対数凹性 ($\sigma$ 非依存)	✅ 厳密
対称 log-concave 理論の適用	❌ 非対称性で不可
${\text{PF}}_{\infty }$ 経路	❌ Hardy 零点で閉じる

RH の幾何学的描像は驚くほどクリアに描ける（素数分布 = フェルミ気体の Fourier スペクトル）にもかかわらず、非対称対数凹関数の Fourier 零点制御という未踏領域で止まっている。

6. 数値的側面：DE 法と整合的な積分表示

この表示は美学にとどまらず、数値計算でも実験的に高い収束性を示す。以下の §6.1–6.2 は数値観察であり、§6.3 は既知の指数収束理論（Trefethen–Weideman 2014 ほか）との整合性の議論である。厳密な誤差評価は離散化・切り詰め・ストリップ幅を明示した上で別途与える必要がある。

6.1 比較：古典 Riemann vs tanh–Mellin

Riemann の古典形式

$$\boxed{\zeta (s)\Gamma (s)={\int }_0^{\infty }\frac{t^{s-1}}{e^t-1}dt}\text{\hspace{1em}(⋆)}$$

に対して、我々の恒等式を同じ Simpson 則で評価する。$s=2$ での誤差：

$N$	($♡$) tanh–Mellin	($\star$) 古典 Riemann
100	$9.4\times 10^{-7}$	$1.1\times 10^{-4}$
200	$6.8\times 10^{-14}$	$5.3\times 10^{-5}$
500	$3.1\times 10^{-15}$（機械精度）	$2.1\times 10^{-5}$

古典表示が $N=1000$ でも $10^{-5}$ 止まりなのに対し、tanh–Mellin は $N=200$ で機械精度に到達する。

6.2 変数変換 $t=e^v$ による台形則

さらに $t=e^v$ で置換し、$\mathbb{R}$ 全体で台形則を適用すると、$s=3$ の誤差は

$N$	誤差
40	$2.7\times 10^{-5}$
80	$2.7\times 10^{-11}$
160	$4.4\times 10^{-16}$（機械 $\epsilon$）

と $N$ 倍増ごとに誤差が 6 桁ずつ 落ちる。これは Simpson の $O(h^4)$ を遥かに超える。

6.3 高速収束の理由（既知理論との整合性）

$h_{\sigma }(v)$ は

• $v\to +\infty$ で二重指数減衰 $e^{-2a{e}^v}$
• $v\to -\infty$ で指数減衰 $e^{\sigma v}$
• 複素ストリップ $|\mathrm{Im}v|<\pi /2$ で解析的（$e^{2a{e}^v}+1$ の零点から決まる最近極までの距離）

の三条件を満たす。Trefethen–Weideman (2014) の「指数収束台形則」の一般論では、実軸上可積分かつ幅 $d$ の水平ストリップで解析的な関数に対し、刻み $h$ の台形則の誤差が $\mathcal{O}(e^{-2\pi d/h})$ で減衰することが示されている。$h_{\sigma }$ はこの仮定（解析性 $d=\pi /2$、急減衰両側）を満たすため、既知定理の守備範囲内で指数収束が期待できる。

この収束構造は Takahashi–Mori (1974) の二重指数求積法 (DE 法) の典型例と同じ形である。通常 DE 法は変数変換を人為的に挿入して被積分核を二重指数化するが、tanh–Mellin 表示は $v=\log t$ 変換のもとでもともと二重指数構造を持っている点が特徴である。

注意. 上の記述は「Trefethen–Weideman の一般定理から指数収束が従う」という既知結果の適用であって、本稿で新たに誤差定数を最適化した定理を主張するものではない。具体的な刻み $h$・切り詰め点 $\pm L$・ストリップ幅 $d$ を明示した非漸近評価は別途整備する必要がある。

観察。 フェルミ気体が熱力学的に「最も単純で普遍」な分布であるように、その Mellin 変換も数値求積と整合的な構造を備えている。この整合性が偶然か必然かは、本稿では判定しない。

7. 実在する物理系と工学的意義

7.1 量子ドット Coulomb blockade

「$1-\tanh (at)$ がフェルミ分布」を抽象モデルにとどめないために、実在する物理系を挙げる。

GaAs/AlGaAs 2 次元電子ガス上の側面量子ドット、WS${}_2$ ファンデルワールス量子ドット、カーボンナノチューブ量子ドット等では、ゲート電圧 $V_g$ を掃引するとコンダクタンス $G$ が Coulomb blockade に由来する鋭いピーク列を呈する。単一準位が熱的に広がる領域での測定ラインシェイプは、Beenakker (1991) が導いた

$$\boxed{G(V_g)=G_{\max }{\cosh }^{-2}\left(\frac{\alpha e(V_g-V_{\text{peak}})}{2k_B{T}_e}\right)}\text{\hspace{1em}(7.1)}$$

で精密にフィットできる。$\alpha$ はゲート→ドットのレバーアーム、$T_e$ は電子温度、$e$ は素電荷。${\cosh }^{-2}={\mathrm{sech}}^2$ は

$$-\frac{d{f}_{\mathrm{FD}}}{dϵ}=\frac{1}{4k_{B}T}{\mathrm{sech}}^2\left(\frac{ϵ}{2k_{B}T}\right)$$

に他ならない。すなわち本稿の被積分核 $1-\tanh (at)$ の導関数が、そのまま実験室で測定される量となっている。

ドット上平均電子数 $⟨N⟩(V_g)={\sum }_n{f}_{\mathrm{FD}}({ϵ}_n-\alpha e{V}_g)$ は $V_g$ の関数として tanh 型ステップ（クーロン階段） を描く。

7.2 典型実験値

希釈冷凍機下の典型パラメータ：

量	代表値
電子温度 $T_e$	$\sim 100$ mK
レバーアーム $\alpha$	$0.05$ – $0.3$
チャージングエネルギー $E_C$	$0.1$ – $1$ meV
ピーク間隔	$\sim$ meV オーダー

フィット式 (7.1) のスケーリング定数は $a=\alpha e/(2k_B{T}_e)$ で、$100\text{mK}$, $\alpha =0.1$ なら $a\sim 5/\text{mV}$。$a$ は連続可変な実験制御量である。

7.3 主恒等式との直接接続

(2.2) で主定理を再掲：

$$\zeta (s)=\frac{a^s{2}^{s-1}}{\Gamma (s)(1-2^{1-s})}{\int }_0^{\infty }[1-\tanh (at)]t^{s-1}dt.$$

量子ドットのゲート応答 $1-\tanh (\alpha e{V}_g/(2k_B{T}_e))$ は、この恒等式の被積分核そのもの。したがって

量子ドット実験で観測されるクーロン階段の $s$ 次 Mellin モーメントは、$\Gamma$・$\eta$ で規格化すれば $\zeta (s)$ に一致する。

状態密度 ${ϵ}^{s-1}$ のフェルミ気体の熱統計量としてゼータが出る、という話は、抽象モデルではなく現に測定されている実デバイスで実現している。

副例として、$S=1/2$ 常磁性体の磁化 $M(B,T)=\frac{1}{2}Ng{\mu }_B\tanh (g{\mu }_{B}B/(2k_{B}T))$ も同じカーネル構造を持ち（Brillouin 関数 $B_{1/2}(x)=\tanh x$）、希薄 Cu${}^{2+}$ 錯体・グラフェン欠陥スピン・常磁性塩等で古典的に検証されている。MOSFET の Random Telegraph Signal にも同型の占有数統計が現れる。

7.4 工学的意義：物理定数 ↔ ゼータ整合性の仮説

ここで本稿の中心的な工学的主張を定式化する。

仮説 7.1（物理–ゼータ整合性予想）. 実測可能な物理系のスケーリング定数 $a$ を精密制御し、被積分核 $1-\tanh (at)$ の $s$ 次 Mellin モーメント

$${\mathcal{I}}_{\exp }(s;a):={\int }_0^{\infty }[1-\tanh (at)]t^{s-1}dt$$

を十分高精度で測定すると、$\zeta (s)$ は

$${\zeta }_{\exp }(s;a):=\frac{a^s{2}^{s-1}}{\Gamma (s)(1-2^{1-s})}{\mathcal{I}}_{\exp }(s;a)$$

として系非依存の普遍量に収束し、精密計算ライブラリの $\zeta (s)$ と任意精度で一致する。

この仮説が現在の実験技術で直接検証されていなくても、数学的には (2.1) の帰結として間違いなく成立する。しかし工学的含意は強烈である：

1. 標準器の数論的基礎づけ。 物理定数（$k_B$, $e$, $\alpha$ 等）の SI 定義は、現在も物理現象（セシウム時計、Klitzing 定数、ジョセフソン定数）に基づく。$\zeta$ は人間が構成した数学的普遍量であり、装置や温度計の特性に完全に独立。${\zeta }_{\exp }(s;a)$ を既知の $\zeta (s)$ と比較することで、逆に $a$（したがって $T_e$, $\alpha$, $k_B$ 等）の絶対校正が可能になる。
2. デバイス非依存のキャリブレーション。 量子ドット、磁化率測定、RTS 統計など異なる物理系で同じ $\zeta (s)$ が出ることが仮説 7.1 の内容。異なる実装で得られる ${\zeta }_{\exp }$ の相互整合性が装置の系統誤差の評価を与える。
3. 医療・生体計測への波及。 ニューロンイオンチャネルの Fermi 型開閉確率、パッチクランプ電流のベース応答、膜電位の sigmoid 型遷移などにも同型カーネルが現れる。$\zeta$ 基準の $s$ 次モーメント評価はドリフト・ノイズ補正の普遍指標を与えうる。
4. RH 証明に依存しない即時応用。 恒等式 (2.2) は $\zeta (s)$ の零点位置に依存せず、$\zeta$ の値自体は既存ライブラリで任意精度計算可能。したがって上記の応用はすべて RH の証明を待たずに開始できる。RH の解決は数学的には決定的だが、工学的には $\zeta$ の数値が確定していれば足り、零点の非自明配置問題とは独立に進められる。

7.5 実験ロードマップの試論

仮説 7.1 を検証するための実験構成を素描する。

• ステップ 1: 単一装置での ${\mathcal{I}}_{\exp }(s;a)$ 測定。 量子ドットの $G(V_g)$ データから $V_g$ の数値 Mellin モーメントを計算し、$a$ を独立測定（熱雑音温度計＋電容性測定によるレバーアーム校正）しておいて ${\zeta }_{\exp }(s;a)$ を組み立てる。$s=2,3,4$ で 4 桁以上の精度での $\zeta$ との一致を目標とする。
• ステップ 2: 異種系の相互整合性。 量子ドットと $S=1/2$ 常磁性体で独立に ${\zeta }_{\exp }$ を測定、相対誤差を評価。両者の系統誤差構造は異なるので整合は系統補正の指標になる。
• ステップ 3: 逆問題。 $\zeta$ を既知として、${\zeta }_{\exp }(s;a)$ から $a$ の絶対値を逆算。独立校正との比較で数論ベースの計測標準としての有用性を判定する。

各ステップの実施可能性は現行技術水準で現実的で、$s=2$ に限定すれば $\sim 10^{-6}$ 相対精度はすでに到達可能な射程と考える（要詳細実験設計）。

8. 副産物と周辺話題

8.1 非零化学ポテンシャルと polylogarithm

化学ポテンシャル $\mu \ne 0$ のフェルミ分布 $1/(e^{ϵ-\mu }+1)$ の $s$ 次モーメントは、$\mathrm{Re}s>0$ のとき一般に polylogarithm で

$${\int }_0^{\infty }\frac{{ϵ}^{s-1}}{e^{ϵ-\mu }+1}dϵ=-\Gamma (s){\mathrm{Li}}_s(-e^{\mu })$$

と表される。$\mu =0$ のとき $-{\mathrm{Li}}_s(-1)=(1-2^{1-s})\zeta (s)=\eta (s)$ に戻り、本稿の恒等式 (2.1) を再現する。特殊なパラメータ化や Mellin–Barnes 表示を経由すれば、Hurwitz ゼータ $\zeta (s,a)$ や Dirichlet $L$-関数との関係式も現れる（例：Lerch 超越関数 $\Phi (z,s,a)$ を介した変形）。いずれも本稿の基本恒等式の自然な拡張系統として整理可能である。

8.2 応答関数の Mellin 階層

$f(t):=1-\tanh (at)$ の $k$ 次導関数は恒等式 (2.1) を $k$ 回部分積分することで、

$$\boxed{{\int }_0^{\infty }{f}^{(k)}(t)t^{s-1}dt=(-1{)}^k\frac{\Gamma (s)\eta (s-k)}{a^{s-k}{2}^{s-k-1}}(\mathrm{Re}s>k)}\text{\hspace{1em}(8.1)}$$

を与える。すなわち 被積分核の微分階数 $k$ だけ $\eta$ の引数がシフトする。境界項は $t\to 0$ で $f^{(j)}(0)$ が全て有限、$t\to \infty$ で $f^{(j)}(t)$ が指数的に減衰することから消える（$\mathrm{Re}s>k$ で $t^{s-1}$ の特異性を吸収）。

これは応答関数論的に自然な階層を与える。$f$ は占有数、$-f^{\prime }=a{\mathrm{sech}}^2(at)$ はその化学ポテンシャル微分であり、Coulomb ピーク（式 (7.1)）・量子容量 $dn/d\mu$・圧縮率 $\partial n/\partial \mu$ などの一次応答そのものである。$f^{\prime \prime }$ はその勾配・曲率、$f^{\prime \prime \prime }$ は高次非線形応答に対応する。したがって

$k$	物理的意味	Mellin モーメントが与えるもの
0	占有数 $⟨n⟩$	$\eta (s)$
1	一次応答（コンダクタンス・圧縮率・量子容量）	$\eta (s-1)$
2	一次応答の勾配・曲率	$\eta (s-2)$
3	高次非線形応答	$\eta (s-3)$

この階層は一つの物理系（例：単一量子ドット）に対して、異なる次数の応答を同一実験で測定すれば異なる $\zeta$ の値が取り出せることを意味する。実験装置の系統誤差は各応答次数で共通する部分と独立な部分に分離されるので、$\zeta (s),\zeta (s-1),\zeta (s-2),\dots$ の整合性（既知値との比較）は RH とは独立な 装置校正指標となる。

数値検証. $k=0,1,2,3$ に対し閉形式導関数を用い、相対誤差は全て $<3\times 10^{-8}$（高次の $s$ では $<5\times 10^{-12}$）で (8.1) を再現した。$k=4,5,6$ は差分微分のため相対精度が劣化するが符号と桁は予測通り。

シミュレーション実験. Coulomb ピーク形状 $G(V_g)\propto {\mathrm{sech}}^2(a(V_g-V_{\text{peak}}))$、$a_{\text{true}}=1.25$、$V_g\in [-12,12]$、$N=2001$ サンプル、加法ガウス雑音 1% の下で、HWHM から $a$ を推定（相対誤差 $8\times 10^{-3}$）し、全実軸上の対称モーメント

$${\mu }_s:={\int }_{-\infty }^{\infty }{\mathrm{sech}}^2(a\tau )|\tau {|}^{s-1}d\tau =2{\int }_0^{\infty }{\mathrm{sech}}^2(a\tau ){\tau }^{s-1}d\tau =\frac{\Gamma (s)\eta (s-1)}{a^s{2}^{s-3}}$$

（$-f^{\prime }=a{\mathrm{sech}}^2(at)$ を (8.1) の $k=1$ に代入し因子 2 を考慮）から $\eta (s-1)$ を抽出し、最終的に $\zeta (s-1)=\eta (s-1)/(1-2^{2-s})$ を復元した：

$s$	復元 $\zeta (s-1)$	真値	相対誤差
3	1.6858	1.6449 ($\zeta (2)$)	2.48%
4	1.2421	1.2021 ($\zeta (3)$)	3.33%
5	1.1276	1.0823 ($\zeta (4)$)	4.18%

1% のノイズ条件で $\zeta (2),\zeta (3),\zeta (4)$ が数 % で復元されることは、仮説 7.1 の弱形（装置キャリブレーション精度での整合）が現実的射程にあることを支持する。誤差の主要成分は $a$ の HWHM 推定誤差が $a^s$ を通じて増幅される寄与で、$a$ の独立校正（温度計＋レバーアーム）で系統的に低減可能である。

8.3 Jensen–Turán 不等式

$\Xi (z)=\sum a_n{z}^{2n}$ に対する Turán 不等式 $a_n^2\ge a_{n-1}{a}_{n+1}$ は RH の必要条件で、Jensen 多項式の全実根性が RH と同値（Polya）。Griffin–Ono–Rolen–Zagier (2019) は各次数で十分大 $n$ の全実根性を証明。本稿の Mellin 分解 $\Xi =A(1/2+iz)\widehat{h_{1/2}}(z)$ は、Turán 条件を $\widehat{h_{1/2}}$ 単独の係数問題と $A$ による補正に分離する視角を与える。

$h_{1/2}$ は命題 4.3 の通り非対称なので $\widehat{h_{1/2}}(z)={\sum }_n{c}_n{z}^n$ の奇数次係数 $c_{2k+1}$ は一般に純虚数（モーメント $\int v^{2k+1}{h}_{1/2}(v)dv\ne 0$ の $i^{2k+1}/(2k+1)!$ 倍）となり、Turán 不等式は $\widehat{h_{1/2}}$ そのものには直接適用できない。したがって偶関数部分 ${\widehat{h_{1/2}}}^{\mathrm{even}}(z):=\frac{1}{2}(\widehat{h_{1/2}}(z)+\widehat{h_{1/2}}(-z))$ の実係数列

$$b_n=\frac{(-1{)}^n}{(2n)!}{\int }_{-\infty }^{\infty }{v}^{2n}{h}_{1/2}(v)dv$$

に対して Turán 条件を評価する（$b_n$ が偶数次係数の係数自体、${\widehat{h_{1/2}}}^{\mathrm{even}}(z)=\sum b_n{z}^{2n}$）。奇関数部分 ${\widehat{h_{1/2}}}^{\mathrm{odd}}$ は $A$ の対応するパリティ成分と結合し $\Xi$ の偶関数性を作る側に回る。

数値観察（$a=1$, $\sigma =1/2$ で実施）. $\Xi$ の Maclaurin 係数 $a_n$ を、$\zeta$ を Euler–Maclaurin 型級数加速（項数 $M=50$）で評価し、$\Xi (z)$ を円周 $|z|=R$ 上で離散 Fourier 展開（点数 $N=2048$、クロスバリデーション半径 $R\in \{1.5,2,2.5,3\}$、主値は $R=1.5$）により $0\le n\le 7$ まで取得した。この条件下で:

• $\Xi$ 係数 $a_n$ の Turán 比 $a_n^2/(a_{n-1}{a}_{n+1})$: $1\le n\le 6$ で $>1$ を維持、$n=7$ でも数値的に $>1$（Griffin–Ono–Rolen–Zagier の理論予想と整合）。
• 偶関数化モーメント $b_n$ の Turán 比 $b_n^2/(b_{n-1}{b}_{n+1})$: $1\le n\le 6$ で $\ge 1$、$n=7$ で $1$ を割る領域に入る兆候が観察される。
• モーメント型 Turán 閾値 ${\mu }_{2n}^2/({\mu }_{2n-2}{\mu }_{2n+2})\ge 1/[(2n+2)(2n+1)]$（${\mu }_{2n}=(-1{)}^n(2n)!a_n$）: $n\le 7$ で全て十分なマージン（最小でも閾値の $10^2$ 倍以上）で成立。

これらは浮動小数点倍精度（有効桁約 15）・積分領域 $|z|\le 3$ 内での観察であり、$n\ge 8$ では $a_n$ の計算精度自体がキャンセレーションで劣化する。したがって本節の主張は「$n\le 7$ の範囲での数値観察」に限定され、全次数での成立や $A$ による補正機構の厳密定式化は別途の研究課題である。特に「$\Xi$ では $A$ との畳み込みで回復する」という構造は Cauchy 積としての観察であって、任意の $b_n$ の挙動から $a_n$ の Turán 性を導く一般定理ではない。

8.4 明示公式の熱力学的読み替え

Riemann–von Mangoldt 明示公式の零点項 ${\sum }_{\rho }{x}^{\rho }/\rho$ は、Hardy カーネルの Fourier スペクトルを $\log x$ 軸でサンプリングしたものに相当する。

素数の分布の揺らぎは、対数エネルギー軸上のフェルミ統計の音響である。

8.5 ボース–フェルミ双対とゲージ階層構造

本稿の主恒等式 (2.1) は Fermi 分布の Mellin 表示だが、Riemann (1859) の古典的ボース表示

$$\Gamma (s)\zeta (s)={\int }_0^{\infty }\frac{t^{s-1}}{e^{at}-1}dt\cdot a^s(\mathrm{Re}s>1)\text{\hspace{1em}(8.2)}$$

も任意の $a>0$ でスケール自由度を持つ。本節では両表示のスケール構造を比較し、Fermi 側が単一核で閉じる明示的なスケール自由度を持つ理由と、ボース–フェルミ等式から派生する非自明な解析的事実を整理する。

8.5.1 ゲージ自由度の非対称性

両表示のゲージ構造を対比する:

	Bose (Riemann 1859)	Fermi (本稿 2.1)
スケール $a>0$	持つ（$a^s$ 因子に吸収）	持つ
化学ポテンシャル $\mu$	$\mu \le 0$ 必須（$\mu >0$ で発散）	$\mu$ 任意（分母常に正）
等式の閉包性	代数恒等式 (8.3) で自己相似階層に分解	単一 $a$ で閉じる

化学ポテンシャル $\mu =0$ に固定したとき、両者とも温度スケール $a=1/(2T)$ の自由度を持つ。しかし代数恒等式

$$\frac{1}{e^{ax}-1}=\frac{1}{e^{ax}+1}+\frac{2}{e^{2ax}-1}\text{\hspace{1em}(8.3)}$$

を再帰適用すると

$$\boxed{\frac{1}{e^{ax}-1}=\sum _{k=0}^{\infty }\frac{2}{e^{2^{k+1}ax}+1}}\text{\hspace{1em}(8.4)}$$

が得られる（数値検証: $a=1,x=3$ で両辺 $5.24\times 10^{-2}$ が 3 桁以上一致、残差は階層の切断による）。すなわちボース分布はフェルミ分布の 2 進自己相似階層 $\{a,2a,4a,8a,\dots \}$ の重ね合わせであり、スケール $a$ を一つ選ぶとその整数倍階層が連動して決まる。この意味でボース側のスケールは階層和に暗に含まれるかたちで現れ、フェルミ側のスケールは単一核で閉じる、より明示的なスケール自由度として振る舞う。

系. フェルミ表示 (2.1) を一つ固定すると、対応するボース表示の階層 (8.4) が一意に決まる。逆向きには、ボース表示は 2 進階層 (8.4) の総和として与えられるため、単一のフェルミ成分（例えば $k=0$ 成分）を自然に選び出すには、階層中のどの項を特権化するかの追加規約が必要となる。

8.5.2 等号から解くスケール比

両表示を同じ $\zeta (s)$ に等号で結び、独立なスケール $a_B$（Bose）と $a_F$（Fermi）を割り当てる:

$$\frac{\Gamma (s)\zeta (s)}{a_B^s}=\frac{2^{1-s}\Gamma (s)(1-2^{1-s})\zeta (s)}{a_F^s}$$

$\zeta (s)\ne 0$ なる領域で整理すると

$$\boxed{{\left(\frac{a_F}{a_B}\right)}^s=2^{1-s}(1-2^{1-s})}\text{\hspace{1em}(8.5)}$$

を得る。右辺は $s$ 依存であり、任意の $s$ で成立する単一比 $a_F/a_B$ は存在しない。しかし各 $s$ で比

$$R(s):=\frac{a_F}{a_B}(s)=[2^{1-s}(1-2^{1-s}){]}^{1/s}$$

が定義される（以下 $s>1$ の実軸上、正の実 $1/s$ 乗 branch を取る）。この関数の解析は以下の非自明な帰結を与える。

事実 1 ($s=2$ での厳密 $1/2$). $R(2)=[2^{-1}\cdot 2^{-1}{]}^{1/2}=1/2$。すなわちモーメント次数 $s=2$（分散・面積モーメント相当）では任意の Bose スケール $a_B$ に対し $a_F=a_B/2$ が厳密に成立する。この関係は $a$ 非依存で、Fermi 分布の特徴的スケールが Bose 分布の半分に対応する現象を普遍的に表す。

事実 2 (UV 固定点). $s\to \infty$ で

$$R(s)=[2^{1-s}(1-2^{1-s}){]}^{1/s}\sim 2^{(1-s)/s}\to \frac{1}{2}.$$

高次モーメント極限で $a_F/a_B\to 1/2$ に収束し、これは階層 (8.4) の最低次項 $k=0$ がボースの半周期に対応する事実の定量的表現である。

事実 3 (IR 退化). $s\to 1^{+}$ で $R(s)\to 0$。$\zeta$ の単純極 $s=1$ では Bose 側 $\int t^{s-1}/(e^{a_{B}t}-1)dt$ が赤外発散する一方 Fermi 側 $(1-2^{1-s})\zeta (s)$ は有限（$\eta (s)$ 側では $\log 2$ に収束）で、両 prefactor の比 $R(s)$ は実数スケール比として退化する。これは Bose/Fermi の解析構造がこの極近傍で質的に異なることを示しており、物理的には理想気体の Bose–Einstein 凝縮における発散的挙動を想起させる。

事実 4 (極大点の存在). $\log R(s)$ は $s\in (1,\infty )$ で唯一の極大点 $s^{\ast }\approx 3.5519$ を持ち、$R(s^{\ast })\approx 0.5766$。数値的には、黄金比・円周率・$e$ などの既知初等定数との単純な代数関係は少なくとも見出せていない（$\phi ,{\phi }^2,\log \phi ,\pi /e$ などとの比較で有効数字 4 桁以上一致する組合せは数値的に検出されない）。$s^{\ast }$ の数論的性格は未解明であり、閉形式の存在有無は今後の課題である。

事実 5 (正の実スケール比の破綻). $s=1/2$（Riemann 臨界線上）で

$$2^{1-s}(1-2^{1-s}){|}_{s=1/2}=\sqrt{2}\cdot (1-\sqrt{2})\approx -0.5858<0$$

となり、$s<1$ の広い領域で内部式が負となる。したがって正の実数比 $R(s)>0$ として実解析的に連続する $a_F/a_B$ のパラメータ化は $s>1$ に限定され、臨界帯 $0<\mathrm{Re}s<1$ では成立しない。

ただしこれは両表示そのものの解析接続が壊れることを意味しない。(8.5) の両辺は複素 $s$ に対し branch を選べば複素値関数として延長でき、$\zeta (s)$ および Mellin 積分の解析接続自体は独立に行える。破綻するのはあくまで**「正の実数温度比を共通に保つ実解析的整合性」**であり、Riemann 臨界線上の情報を得るには両表示をこの実数比の枠組みの外で（複素 branch 込み、または独立な解析接続で）扱う必要がある。

8.5.3 物理的解釈

事実 1–5 を統合すると、$R(s)$ は次のような $s$ 軸上の流れを描く:

領域	$R(s)$ の挙動	物理的対応
$s\to 1^{+}$ (IR)	$0$	Bose 凝縮発散, Fermi 退化
$s^{\ast }\approx 3.55$	$0.577$ (極大)	未解明
$s\to \infty$ (UV)	$1/2$	階層 (8.4) の最低項 $k=0$ が支配

$s$ はモーメント次数であって結合定数ではないが、この構造は繰り込み群の UV/IR 不動点構造と形式的類似を持つ。$s$ 軸を「有効理論の記述次元」と読む枠組みがあれば、ボース–フェルミ双対はその枠組みの内部で自然な RG 流を持つ解析的対象となる。この解釈の物理的根拠の確立は今後の課題である。

8.5.4 副作用: 応答関数階層の二重化

§8.2 の応答関数階層 (8.1) はフェルミ側 $f(t)=1-\tanh (at)$ の微分系列だが、ボース側 $g(t)=1/(e^{at}-1)$ にも同様の階層が存在する:

$${\int }_0^{\infty }{g}^{(k)}(t)t^{s-1}dt=(-1{)}^k\frac{\Gamma (s)\zeta (s-k)}{a^{s-k}}(\mathrm{Re}s>k+1)\text{\hspace{1em}(8.6)}$$

Fermi 側の $\eta$ に対しボース側は $\zeta$ そのものが出る点が差異で、代数恒等式 (8.3) の微分から両階層の係数関係式が従う。したがって Bose 応答と Fermi 応答を独立に測定できる実験系では、$\zeta (s-k)$ と $\eta (s-k)$ が原理上それぞれ独立に制約され、校正に使える独立な関係式が増える可能性がある。実際に自由度が倍増するかは、測定系でのノイズ相関や $a_B,a_F$ を別個に校正できるかに依存する。これは §7 の工学的意義に対する副次的拡張の候補である。

9. 結論

1. 恒等式 (2.2) は Fermi–Dirac 分布の $s$ 次 Mellin モーメントと $\zeta (s)$ を完全に同一視する。
2. Hardy カーネル $h_{\sigma }$ の Fourier 零点は $\zeta$ の非自明零点と 1:1 対応し、RH は $h_{\sigma }$ の Fourier 非零性問題に書き換わる。
3. 対数凹性は成立するが非対称性により現時点の古典定理だけで RH は解けない。
4. 数値的には DE 構造を内包、$160$ 点台形則で機械精度。
5. 物理的には量子ドット・常磁性体・RTS で実測されている核であり、仮説 7.1 の下で物理定数とゼータの整合性は実験的に検証可能。
6. パラダイムとしての数論的物理学：本恒等式は類似ではなく恒等式の水準で数論と統計力学を等値し、Hilbert–Polya 以降で初の恒等式レベルでの橋渡しを与える。
7. 工学的意義：$\zeta$ の値は RH とは独立に既存ライブラリで任意精度計算可能であり、恒等式 (2.2) は零点位置にも依存しない。したがってRH の証明を待たずに、ゼータ関数そのものを精密計測・標準器校正・量子デバイス設計・生体電気計測へ応用できる。

リーマン予想は、式のレベルでは統計力学の最も基本的な分布関数の Fourier 性質の問題である。そしてこの問題は、もはや純粋数学のものとは言い切れない。

参考文献

数論・ゼータ関数

1. E. C. Titchmarsh, The Theory of the Riemann Zeta-function, 2nd ed., revised by D. R. Heath-Brown, Oxford University Press, Oxford (1986). ISBN 0-19-853369-1.
2. H. M. Edwards, Riemann's Zeta Function, Pure and Applied Mathematics 58, Academic Press, New York (1974); reprinted by Dover, Mineola, NY (2001). ISBN 0-486-41740-9.
3. A. Ivić, The Riemann Zeta-Function: Theory and Applications, John Wiley & Sons, New York (1985); reprinted by Dover, Mineola, NY (2003). ISBN 0-486-42813-3.
4. H. L. Montgomery, R. C. Vaughan, Multiplicative Number Theory I. Classical Theory, Cambridge Studies in Advanced Mathematics 97, Cambridge University Press (2007). ISBN 978-0-521-84903-6.
5. G. H. Hardy, Sur les zéros de la fonction $\zeta (s)$ de Riemann, Comptes Rendus Acad. Sci. Paris 158 (1914), 1012–1014.
6. G. Pólya, Über die algebraisch-funktionentheoretischen Untersuchungen von J. L. W. V. Jensen, Kgl. Danske Vidensk. Selsk. Math.-Fys. Medd. 7 (1927), 3–33.
7. M. Griffin, K. Ono, L. Rolen, D. Zagier, Jensen polynomials for the Riemann zeta function and other sequences, Proc. Natl. Acad. Sci. USA 116 (2019), no. 23, 11103–11110. DOI: 10.1073/pnas.1902572116. arXiv:1902.07321.

零点統計とランダム行列

8. A. M. Odlyzko, On the distribution of spacings between zeros of the zeta function, Math. Comp. 48 (1987), no. 177, 273–308. DOI: 10.1090/S0025-5718-1987-0866115-0.
9. M. V. Berry, J. P. Keating, The Riemann zeros and eigenvalue asymptotics, SIAM Review 41 (1999), no. 2, 236–266. DOI: 10.1137/S0036144598347497.
10. M. L. Mehta, Random Matrices, 3rd ed., Pure and Applied Mathematics 142, Elsevier/Academic Press, Amsterdam (2004). ISBN 978-0-12-088409-4.

対数凹性・全正値性・特性関数

11. I. J. Schoenberg, On Pólya frequency functions I. The totally positive functions and their Laplace transforms, J. Anal. Math. 1 (1951), 331–374. DOI: 10.1007/BF02790092.
12. S. Karlin, Total Positivity, Vol. I, Stanford University Press, Stanford (1968).
13. E. Lukacs, Characteristic Functions, 2nd ed., Griffin, London (1970). ISBN 978-0852640517.
14. S. G. Bobkov, Asymptotic expansions for products of characteristic functions under moment assumptions of non-integer orders, in: E. Carlen, M. Madiman, E. Werner (eds.), Convexity and Concentration, IMA Volumes in Mathematics and its Applications 161, Springer, New York (2017), 297–357. DOI: 10.1007/978-1-4939-7005-6_11.

統計力学

15. L. D. Landau, E. M. Lifshitz, Statistical Physics, Part 1, 3rd ed., Course of Theoretical Physics 5, Pergamon Press, Oxford (1980). ISBN 0-08-023039-3.
16. R. K. Pathria, P. D. Beale, Statistical Mechanics, 3rd ed., Elsevier/Academic Press (2011). ISBN 978-0-12-382188-1.

数値解析（二重指数求積）

17. H. Takahasi, M. Mori, Double exponential formulas for numerical integration, Publ. Res. Inst. Math. Sci. (RIMS, Kyoto Univ.) 9 (1974), no. 3, 721–741. DOI: 10.2977/prims/1195192451.
18. L. N. Trefethen, J. A. C. Weideman, The exponentially convergent trapezoidal rule, SIAM Review 56 (2014), no. 3, 385–458. DOI: 10.1137/130932132.

量子ドット・Coulomb Blockade

19. C. W. J. Beenakker, Theory of Coulomb-blockade oscillations in the conductance of a quantum dot, Phys. Rev. B 44 (1991), no. 4, 1646–1656. DOI: 10.1103/PhysRevB.44.1646.
20. H. van Houten, C. W. J. Beenakker, A. A. M. Staring, Coulomb-blockade oscillations in semiconductor nanostructures, in: H. Grabert, M. H. Devoret (eds.), Single Charge Tunneling, NATO ASI Series B 294, Plenum, New York (1992), 167–216. arXiv:cond-mat/0508454.
21. A. A. M. Staring, Coulomb-Blockade Oscillations in Quantum Dots and Wires, PhD thesis, Technische Universiteit Eindhoven (1992). Available from Philips Research Laboratories / Eindhoven University of Technology institutional repositories.
22. X.-X. Song, Z.-Z. Zhang, J. You, D. Liu, H.-O. Li, G. Cao, M. Xiao, G.-P. Guo, Temperature dependence of Coulomb oscillations in a few-layer two-dimensional WS${}_2$ quantum dot, Sci. Rep. 5 (2015), 16113. DOI: 10.1038/srep16113. arXiv:1510.03583.
23. D. Maradan, L. Casparis, T.-M. Liu, D. E. F. Biesinger, C. P. Scheller, D. M. Zumbühl, J. D. Zimmerman, A. C. Gossard, GaAs quantum dot thermometry using direct transport and charge sensing, J. Low Temp. Phys. 175 (2014), no. 5–6, 784–798. DOI: 10.1007/s10909-014-1169-6. arXiv:1401.2330.
24. I. L. Aleiner, P. W. Brouwer, L. I. Glazman, Quantum effects in Coulomb blockade, Phys. Rep. 358 (2002), no. 5–6, 309–440. DOI: 10.1016/S0370-1573(01)00063-1.