1次元更新接触過程の改良された生存結果 — 算術的・非原子的到着時間分布への臨界感染パラメータの有限性

§01 接触過程と更新接触過程——マルコフ性の喪失がもたらす豊かさと困難

接触過程（contact process）は 1974 年に Ted Harris が導入した相互作用粒子系の古典的モデルです。整数直線 $\mathbb{Z}$ の各格子点は「感染」（$\eta(x) = 1$）または「健康」（$\eta(x) = 0$）の二状態を取ります。時間の経過に伴い、感染点 $x$ は独立な指数分布 $\mathrm{Exp}(1)$ に従う待機時間の後に回復し、同時に各隣接点 $y \sim x$ に独立なレート $\lambda > 0$ の指数時計で感染を試みます。感染率 $\lambda$ はモデルの唯一のパラメータであり、$\lambda$ の増大に伴って「感染の拡散力」が強まります。

この古典的接触過程の本質的な特徴は「マルコフ性」です。将来の状態は現時点の状態のみに依存し、過去の感染履歴は関係ありません。この性質は指数分布の「無記憶性（memoryless property）」から直接導かれます——過去にどれだけ長く感染していようとも、残余の回復待機時間の分布は常に $\mathrm{Exp}(1)$ に等しい。このマルコフ性こそが古典的接触過程の解析を扱いやすくし、半世紀以上の豊富な理論の蓄積を可能にしてきました。

更新接触過程（renewal contact process）はこのマルコフ的回復を一般化します。各感染点の回復時刻は、独立な**更新過程**（renewal process）によって支配されます。より具体的には、感染点 $x$ が時刻 $s$ に感染した場合、その回復時刻は $s + \tau_1$ です。ここで $\tau_1, \tau_2, \tau_3, \ldots$ は**到着時間分布**（interarrival distribution）$\mu$ に従う独立同分布列であり、$\tau_1$ はその最初の要素です。$\mu = \mathrm{Exp}(1)$ のとき更新接触過程は古典的接触過程に完全に一致しますが、一般の $\mu$ に対してはマルコフ性が失われます。

非マルコフ的なダイナミクスは数学的に著しく困難です。通常の接触過程では「現在の配置 $\eta \in \{0,1\}^{\mathbb{Z}}$」が状態空間を構成しますが、更新接触過程では各感染点が「いつ感染したか」という情報——言い換えると現在の「年齢（age）」——が状態記述に必要となります。状態空間は無限次元的な拡張を要し、通常のマルコフ連鎖理論が直接適用できません。この「マルコフ性の喪失」こそが本論文を含む更新接触過程の研究における中心的な挑戦です。

更新過程の理論では、到着時間分布 $\mu$ の尾部（tail）の挙動が系の長時間挙動を左右します。平均到着時間 $m = \int_0^\infty t \, \mu(dt)$ が有限のとき、更新定理（renewal theorem）は単位時間当たりの更新頻度の収束 $U([t, t+h]) / h \to 1/m$（$t \to \infty$）を保証します。重尾の場合（$m = +\infty$）には更新頻度が低下しますが、感染の持続が長引きやすくなるという直感も働きます。本論文の問いはこの直感の精密化です：$\mu$ がどのような性質を持てば、感染は十分高い感染率のもとで正の確率で永遠に生存できるのか。

§02 臨界感染パラメータ $\lambda_c(\mu)$ と生存問題の定式化

更新接触過程（そして古典的接触過程も含む相互作用粒子系全般）の中心的な問いは「生存（survival）か絶滅（extinction）か」という相転移の解析です。感染率 $\lambda$ を変化させると過程の定性的な振る舞いが質的に変化します。

感染率が十分小さい場合、任意の有限初期配置から出発した過程は有限時間で確率 1 で絶滅します（全格子点が健康状態に達する）。一方、感染率が十分大きい場合には、原点のみを感染点として出発した過程が正の確率で永遠に生存することが期待されます。この相転移の閾値が**臨界感染パラメータ**（critical infection parameter）$\lambda_c(\mu)$ です。

古典的接触過程（$\mu = \mathrm{Exp}(1)$）では $\lambda_c < +\infty$ が長年にわたり証明されており（1 次元では $\lambda_c \approx 1.64$）、$\lambda > \lambda_c$ での生存は確立された事実です。しかし更新接触過程では、$\mu$ の選び方によって $\lambda_c(\mu) = +\infty$ となる場合があります。これは感染率をどれだけ大きくしても常に絶滅してしまうことを意味し、「生存不能（no survival for any λ）」という極端な状況です。

$\lambda_c(\mu) = +\infty$ となる状況の直感的な理解は次のとおりです。$\mu$ の裾が非常に重い場合（たとえば $\mu([t, \infty)) \sim t^{-\alpha}$、$\alpha \leq 1$）、ある感染点が極めて長時間にわたって感染し続ける確率が非無視的に大きくなります。しかし同時に、周囲への感染が成功しても隣接点がすぐに回復してしまうことで、感染の「橋頭堡」を確立することが難しくなります。この両効果のバランスが $\lambda_c(\mu)$ の有無を決めます。

Fontes、Marchetti、Mountford、Vares による先行研究（2023 年）は、この問題に対して重要な前進をもたらしました。彼らの研究では、更新測度が絶対連続成分を持つ場合や、到着時間分布が特定の正則性条件を満たす場合に $\lambda_c(\mu) < +\infty$ が証明されています。しかし、算術的分布（格子上に台を持つ分布）と特異連続分布（ルベーグ連続でも原子もない分布、例えばカントール分布）への対応は不完全なまま残されていました。本論文はこの「隙間」を埋めることを目標とします。

生存問題の数学的な定式化では、原点のみを初期感染点とした過程 $\{\xi_t\}_{t \geq 0}$（$\xi_0 = \{0\}$）に対して、生存事象 $\{\xi_t \neq \emptyset, \, \forall t \geq 0\}$ の確率の正値性を問います。この定式化は「弱生存（weak survival）」と呼ばれ、過程が永遠に消えないことを意味します。強生存（strong survival）——感染域が空間的に広がり続けること——はより強い条件ですが、1 次元では弱生存と強生存は等価であることが知られています。この等価性が 1 次元の設定を特別にする理由の一つです。

§03 主定理の二本柱——算術的分布と更新測度の原子成分による統一条件

本論文の中心的な成果は二つの定理として表現されます。これらは到達範囲の異なる独立した結果でありながら、証明の戦略を共有する相互補完的な構造を持ちます。

**第一の主結果（算術的分布への有限性）**: $\mu$ が非退化の算術的到着時間分布であるならば、$\lambda_c(\mu) < +\infty$。

「算術的（arithmetic）」とは、$\mu$ の台が $h\mathbb{Z}_{\geq 0} = \{0, h, 2h, 3h, \ldots\}$ の形の格子に含まれることを意味します（ここで $h > 0$ は**格子スパン**）。言い換えると、更新の到着時刻が常に $h$ の整数倍であることです。整数値確率分布（$h = 1$）は算術的分布の典型例であり、コイン投げ回数、離散的なイベント時刻など、実際の応用でよく現れます。「非退化」とは $\mu$ が単一の原子 $\delta_{kh}$（$k \geq 1$）に集中しないことを意味し、すべての到着間隔が定数である退化した場合を除外します。

算術的分布の難しさは逆説的です。格子構造によって時刻が離散化され、更新の到着が予測可能なタイミングに集中します。これは一方では組合せ的な議論を助けますが、他方では「到着が集中する時刻」に更新測度が大きな原子を持つことを意味し、確率解析上の特殊な扱いが必要となります。

**第二の主結果（一般条件と非原子的分布への応用）**: 更新測度 $U$ の原子成分 $U^{\mathrm{at}}$ が、十分短い区間上で一様に小さい——すなわち適切な $\varepsilon, \delta > 0$ が存在して $\sup_{t \geq 0} U^{\mathrm{at}}(\{t\}) < \varepsilon$（またはその変形）——ならば $\lambda_c(\mu) < +\infty$。

ここで更新測度 $U$ は $U(A) = \sum_{n=0}^\infty P(S_n \in A)$ であり、その原子成分 $U^{\mathrm{at}}$ は $U$ の孤立点への集中を測ります。直感的には「特定の時刻 $t$ に確率的に到着しやすさの上界」を与える量です。$\mu$ が原子を持たない（非原子的な）場合、更新測度 $U$ の原子成分は消える——すなわち $U^{\mathrm{at}} = 0$——ことが古典的な更新論の結果から導かれます。したがって、第二の定理の条件は非原子的なすべての $\mu$ に対して自動的に成立します。

この「非原子的分布」のクラスは絶対連続分布（密度関数を持つ分布）だけでなく、**特異連続分布（singular continuous distribution）**をも含みます。特異連続分布の典型例はカントール分布であり、ルベーグ測度に対して絶対連続でも原子も持ちません。従来の証明手法はしばしば絶対連続性を仮定しており、特異連続分布は「灰色地帯」として残されていました。本論文の第二の結果がこの灰色地帯を完全に照らします。

二つの結果をつなぐ概念上の橋渡しは「原子成分の制御」です。算術的分布の場合には格子構造を利用して原子の集中を別の角度から制御し、非原子的分布の場合には $U^{\mathrm{at}} = 0$ という直接的な事実を用います。この統一的な視点が本論文の理論的な凝集力を生み出しています。

§04 証明の三本柱——更新測度の局所推定・再生有向浸透・パイアルス論法

本論文の証明は互いに補完し合う三つの技術的要素から構成されます。それぞれが単独でも意義深い手法ですが、それらを有機的に組み合わせることで初めて多様な分布クラスを統一的に扱える証明が完成します。

**第一の柱：更新測度の局所推定（local estimates for renewal measures）**

更新測度 $U([t, t+h])$ の精密な推定が証明の出発点です。古典的な更新定理（renewal theorem）は長時間での漸近挙動 $U([t, t+h]) \to h/m$（$t \to \infty$）を保証しますが、有限時間での精密な上界や、原子成分の制御は別途の議論を要します。

算術的分布の場合、格子スパン $h$ での更新点 $U(\{kh\})$ の推定が鍵です。$S_n = \sum_{i=1}^n \tau_i$ が値 $kh$ を取る確率は $P(S_n = kh) = \mu^{*n}(\{kh\})$（$n$ 重畳み込み）であり、$\mu$ の非退化性から「ある $k$ と $k'$ が互いに素」という条件が成立すれば、$\sum_{n} P(S_n = kh)$ の収束性が導かれます。この格子上の再帰確率の精密な評価が、算術的分布の場合の技術的な核心です。

非原子的分布の場合には、$U$ が原子を持たないことから $U(\{t\}) = 0$ が直接成立します。しかし証明では、短区間 $[t, t+\delta]$ 上での $U$ の集中——「局所的な更新頻度の突発的な増大」——が起きないことの保証が必要です。この局所的推定は更新過程の理論における独立した技術的結果として意義を持ちます。

**第二の柱：再生有向浸透モデルへの比較（comparison to regenerative oriented percolation）**

更新接触過程の生存問題を、より扱いやすいモデルである「再生有向浸透（regenerative oriented percolation）」に比較することが証明の核心的なステップです。有向浸透（oriented percolation）は $\mathbb{Z} \times \mathbb{Z}_{\geq 0}$（空間方向の整数直線と非負整数の時間軸の積空間）上のモデルです。各格子点 $(x, t)$ から $(x-1, t+1)$ および $(x+1, t+1)$ への辺がそれぞれ独立に確率 $p$ で開き（open）、情報が時間の正方向に伝播します。有向浸透には明確な臨界確率 $p_c^{\mathrm{OP}}$ が存在し、$p > p_c^{\mathrm{OP}}$ ならば正の確率での生存が保証されることが知られています。

著者の方々が導入する「再生有向浸透」は、有向浸透に更新過程の到着時刻による時間同期を組み込んだモデルです。各サイトの活性化（感染の可能性）が更新過程のランダムな到着時刻と連動するように設計されており、更新接触過程のダイナミクスをこの再生有向浸透に確率的に支配（stochastic domination）できます。感染率 $\lambda$ が十分大きければ再生有向浸透の辺開確率が $p_c^{\mathrm{OP}}$ を超え、ひいては更新接触過程の生存が導かれます。この「上から押さえる比較」を有効にするためには、更新測度の局所推定（第一の柱）によって更新到着の時間的な「集中」が制御されていることが前提となります。

**第三の柱：パイアルス型輪郭論法（Peierls-type contour argument）**

パイアルス論法は統計力学の古典的な技法です。Rudolf Peierls が 1936 年に 2 次元イジングモデルの相転移を証明した際に導入されて以来、接触過程を含む多くの相互作用粒子系への応用が蓄積されています。

基本的なアイデアは「感染が有限時間内に絶滅する」という事象を、時空間上の「輪郭（contour）」——感染域と非感染域の境界曲線——によって記述することです。輪郭の出現確率を「輪郭の大きさ（長さや複雑さ）」に対するペナルティとして上から評価します。感染率 $\lambda$ が大きいほど各輪郭が現れる確率は低下するため、十分大きな $\lambda$ ではすべての輪郭の出現確率の和が 1 未満となり、感染が絶滅しない確率の正値性が示されます。

更新接触過程のパイアルス論法では、輪郭を構成する時間区間が更新過程の到着時刻と整合するよう取ることが不可欠です。算術的分布の場合、輪郭を格子時刻（$h$ の整数倍）に沿って構成することで離散的な組合せ的数え上げが適用できます。非原子的分布の場合には、更新測度の推定によって「輪郭の境界点での更新イベントの密度」が制御され、連続時間上でのパイアルス論法が成立します。この連続時間版パイアルス論法の精密な構成が、本論文の技術的な貢献の一部です。

三本柱の有機的な接続をまとめると：更新測度の局所推定が再生有向浸透への比較を可能にし、有向浸透の生存がパイアルス輪郭論法によって証明された臨界性の超越により確認される——という推論の連鎖として証明が完成します。この構造は、将来の更新接触過程の研究に対して再利用可能な証明の骨格を提供するものでもあります。