Double Chooz、残留反ニュートリノ束の初の定量測定を発表

停止した原子炉は静かではない。炉心とプールに蓄積した核分裂生成物は数か月、あるいは数年にわたって崩壊を続け、そのベータ崩壊系列から反ニュートリノが放出される。Double Chooz共同研究グループは、Chooz原子力発電所の両原子炉が同時に停止した機会を利用して得た、この残留束の初の定量測定を発表した。

出典：journals.aps.org

わかりやすく言うと

薪の火を消しても、残り火は長いあいだ赤く燃え続ける。原子炉でも似たことが起きるが、その輝きは目には見えない。炉心に残る照射済み燃料や、プールに保管された使用済み燃料集合体に含まれる核分裂生成物のベータ崩壊系列が、幽霊のような粒子である反ニュートリノの流れを生み出す。

アルデンヌ地方にあるChoozの両原子炉が停止していた有効測定時間17.2日間に、炉心から約400メートル離れ、115メートル水換算の遮蔽で守られた近距離検出器は、バックグラウンド差し引き後の残留候補事象を106 ± 18件分離した。これは一日あたり約六件に相当する。詳細モデルでは信号の44%が使用済み燃料を冷却するプールから来ると予測されるが、検出器がこの内訳を直接分離したわけではない。

想定される実用上の関心は、核物質の監視である。反ニュートリノは相互作用する確率が極めて低いまま物質を通過するため、炉心に立ち入らずに原子炉を観測できる。ただし注意が必要だ。この結果は数十立方メートルのシンチレータと、115メートル水換算の地下遮蔽に依存している。残留束の計算を検証するものではあるが、可搬型装置が発電所の地表で機能することを示したものではない。

技術的な説明

1. 反ニュートリノをどう捉えるのか、そしてなぜ多数が必要なのか。 用いる反応は逆ベータ崩壊（IBD）である：$\bar{\nu}_e + p \rightarrow e^+ + n$νˉe​+p→e++n。電子反ニュートリノが自由陽子、すなわちシンチレータ中の水素原子核と相互作用し、中性子と陽電子を生成する。署名は二重であり、一日あたり数事象という低い率でも測定を可能にするのはこの特徴である。陽電子は直ちに運動エネルギーを落としてから対消滅し、即発信号を作る。中性子は拡散した後に捕獲される。Double Choozの「total neutron capture」法は、内部体積のガドリニウム、水素、炭素への捕獲を利用して遅延信号を作る。二つの信号が0.5から800 µsの時間窓に入り、互いに1.2 m未満であることを要求すると、環境放射能の大部分を除去できる。反ニュートリノのエネルギーは即発信号から$E_{\bar{\nu}_e} \simeq E_{\text{prompt}} + 0.78\ \text{MeV}$Eνˉe​​≃Eprompt​+0.78 MeVとして再構成され、この反応には約1.8 MeVの運動学的しきい値がある。それより下では、このチャンネルは反ニュートリノを検出しない。

2. 残留する輝きはどこから来るのか：三つの系列、三つの時計。 稼働中の原子炉の反ニュートリノ束は、主に短寿命の核分裂生成物に由来し、その寄与は停止後数時間で急減する。検出可能な残りは主に三組の親娘核種から来る。照射中とその後の冷却中、長寿命の親核種が、はるかに短寿命の娘核種を絶えず供給する。娘核種は生成のたびに速やかに消えるため、その放射能は主として親核種の放射能に追随する。停止によって失われるのは、こうした持続的な系列を覆い隠していた支配的な束である。

   核種対	親核種の半減期	娘核種：半減期と最大エネルギー$Q_\beta$Qβ​	役割
   $^{144}\mathrm{Ce} \rightarrow {}^{144}\mathrm{Pr}$144Ce→144Pr	285日	17.3 min · 3.00 MeV	モデルが予測する1–3 MeV信号の約54%
   $^{106}\mathrm{Ru} \rightarrow {}^{106}\mathrm{Rh}$106Ru→106Rh	372日	30.1 s · 3.54 MeV	モデルが予測する信号の約38%
   $^{90}\mathrm{Sr} \rightarrow {}^{90}\mathrm{Y}$90Sr→90Y	28.9年	3.19 h · 2.28 MeV	長期的なモデル予測で支配的になる（十年後に> 90%）

   これらのエネルギーが以後のすべてを決める。$Q_\beta$Qβ​が2から3.5 MeVで、IBDしきい値が1.8 MeVであるため、各ベータスペクトルの高エネルギー側の裾だけが検出可能である。そのため著者らは、検出可能な残留束の98.7%が3 MeV未満にあると予測する。これが解析を即発エネルギー1–3 MeVの窓に限定する根拠であり、同時に解析を難しくする理由でもある。バックグラウンドが最も密集する領域だからだ。

3. 1–3 MeVの窓が正確には何を測り、それをどの対照が検証するのか。 信号が予想される窓に解析を限定すると、「探した場所で探したものを見つけただけではないか」という典型的な反論が生じる。著者らは負の対照でこれに答える。3 MeVを超える領域では、バックグラウンド差し引き後の率はゼロと整合する。モデルは小さな残留テールを予測するものの、一日あたり0.1事象未満であり、現在の感度を下回る。この領域で測定可能な超過がないことはバックグラウンドモデルを支持するが、物理的な束が厳密にゼロであることを証明するものではない。

4. バックグラウンドと、地下遮蔽が重要な理由。 一日あたり六事象では、すべてが重要になる。三種類を別々に扱う。偶発事象は、無関係な二つのエネルギー沈着の一方を即発、もう一方を遅延と誤認したもので、除去能力が400を超えるニューラルネットワークで排除し、残りは孤立した沈着率からin situで測定する。相関事象は宇宙線ミューオンに由来する。岩盤内の核破砕で生じ、即発信号を模倣してから捕獲される高速中性子、または$^{9}\mathrm{Li}$9Liのようにシンチレータ中の$^{12}\mathrm{C}$12Cに対するミューオン核破砕で生成されるβ-n同位体である。後者の崩壊だけで電子と中性子の対が生じる。これらは親ミューオンとの時間的・空間的相関でタグ付けされる。近距離検出器を守る115 m水換算の遮蔽はミューオン束を減らし、どの veto でも完全には除去できないこれらの同位体の生成を抑える。

5. 予測：燃料集合体ごとのインベントリ計算。 期待数は概算ではなく、詳細なシミュレーションである。予測IBDスペクトルは
   $\frac{dN_{\bar{\nu}_e}}{dE_{\text{vis}}} = \sum_i \frac{N_p}{4\pi L_i^2}\iint \varphi_i(E,t)\,\sigma_{\text{IBD}}(E)\,P_{ee}(E,L_i)\,\mathcal{D}(E_{\text{vis}}|E)\,dE\,dt,$dEvis​dNνˉe​​​=i∑​4πLi2​Np​​∬φi​(E,t)σIBD​(E)Pee​(E,Li​)D(Evis​∣E)dEdt,
   と書ける。ここで$N_p$Np​は標的陽子数、$L_i$Li​は線源$i$iまでの距離、$\varphi_i$φi​は時間依存の放出スペクトル、$\sigma_{\text{IBD}}$σIBD​は断面積、$P_{ee}$Pee​は振動による生存確率、$\mathcal{D}$Dは検出器応答である。核分裂生成物の放射能は核データライブラリJEFF-3.1.1を用いるAPOLLO-2.8.4とDARWIN-3のコードから得られる。ベータスペクトルにはBESTIOLEライブラリを使い、フェルミ理論に電磁補正、有限核サイズ、原子遮蔽、許容・禁制遷移の形状因子を適用する。各燃料集合体について、全寿命にわたる熱出力を追跡し、各停止期間の開始まで照射履歴を個別にシミュレーションした。

6. モデルが44%をプールに帰属させる理由。 二つの貯蔵プールは炉心に近いため、約400 m離れた検出器から見た幾何学因子$1/L^2$1/L2は同程度のままである。各集合体は複数回の照射サイクルを経てプールに移され、そこで数年間冷却される。停止した炉心には急速に崩壊する新しいインベントリがある一方、プールには年代の異なる集合体が蓄積し、一部の同位体は長く残る。長期停止では、予測されるプール寄与の合計が炉心の寄与と同程度になる。積分計算は炉心に56%、プールに44%を与えるが、現在の幾何配置ではこの二つの割合を個別に測定できない。

7. 予測を制限するものと、そこから結論できないこと。 予測事象数の規格化不確かさは7.4%に達し、内訳はスペクトルモデル化6.0%、ベースラインの幾何2.9%、核分裂生成物インベントリ2.1%、プール内の燃料集合体インベントリ2.0%、検出器が0.8%未満である。支配項は$^{144}\mathrm{Pr}$144Prのベータ遷移形状に由来し、その不確かさは詳細な核構造計算と$\xi$ξ近似との差として見積もられる。この7.4%はデータとモデルの一致精度でも、インベントリ推定性能でもない。この研究は計数と予測を比較したもので、燃料の質量や組成を逆推定してはいない。

図の表記：$N_{\rm obs}$Nobs​は測定された残留計数、$N_{\rm pred}$Npred​はモデル予測、$H_0(B)$H0​(B)はバックグラウンドのみの仮説である。$5.9\,\sigma$5.9σは$N_{\rm obs}$Nobs​と$H_0(B)$H0​(B)を比較したもので、$N_{\rm pred}$Npred​との比較ではない。

なぜ成功したのか

この測定は、ある状況と幾何配置によって成り立った。状況とは、2017年にChoozの両原子炉が四回同時停止し、そのうち一回は20.8日間続いたことである。近距離検出器は全出力時には一日あたり約900件のIBD事象を数えるのに対し、ここでの残留事象は一日あたり6.2 ± 1.1件なので、これは不可欠な条件だった。幾何配置とは、すでに建設され、較正され、特性が理解された検出器が400 mの位置にあり、長年の振動測定によって系統不確かさが確立されていたことである。

数値比較は直接的である。2011–2012年のデータを用いた先行解析では候補が約二十件しか集まらず、残留スペクトルを抽出するには少なすぎた。今回はバックグラウンド差し引き後の残留候補事象106 ± 18件を、期待値88 ± 7件と比較する。差は$\Delta = 18 \pm 19$Δ=18±19事象であり、両者は整合する。遠距離検出器では期待値14 ± 1件に対して27 ± 13件で、誤差幅が大きく、結論ははるかに緩い。

限界は明示されている。信号解析は1–3 MeVの窓を対象とする。3 MeVを超える領域では差し引き後の率はゼロと整合する一方、モデルには約4.5 MeVまで、一日あたり0.1事象未満の小さな残留テールが残るが、この標本では分解できない。7.4%の規格化不確かさは予測を特徴づけるものであり、実験がインベントリ推定の精度を実証したわけではない。

実証されたことと、将来構想との隔たりを明確にする必要がある。実証されたのは、残留束がバックグラウンドから$5.9\,\sigma$5.9σで識別され、計数$106\pm18$106±18が予測$88\pm7$88±7と整合することだ。これは$\Delta=18\pm19$Δ=18±19だからである。実証されていないのは、燃料の運用監視である。この実証には115 m水換算の地下にある30 m³超のシンチレータを使い、400 mの距離で得られるのは一日あたり約六事象である。地表の小型装置には、宇宙線バックグラウンドを抑え、インベントリ変化への感度を実証する課題が残る。

因果の連鎖

炉心での核分裂 → 照射中に、短寿命のベータ崩壊娘核種を供給する長寿命の親核種（$^{144}\mathrm{Ce}$144Ce、$^{106}\mathrm{Ru}$106Ru、$^{90}\mathrm{Sr}$90Sr）が蓄積 → 原子炉停止と支配的な束の急速な消失 → 炉心とプールへ移された集合体ですでに成立していた系列が残存 → 2017年に合計24.4日の同時停止四回 → 近距離検出器の有効測定時間17.2日 → 1から3 MeVで残留候補事象$106\pm18$106±18、バックグラウンドのみとは$5.9\,\sigma$5.9σで不整合 → 3 MeV超に予測される一日あたり0.1事象未満の小さなテールは未分解 → 集合体ごとのモデル（$88\pm7$88±7、$\Delta=18\pm19$Δ=18±19）と整合 → 残留束計算を検証する初の定量的基準 → 炉心とプールの分離、および専用検出器に関する将来研究。

こぼれ話

Double Choozは2011–2012年のデータでもこの測定を試みていた。当時、共同研究グループが集めた残留候補は約二十件にすぎなかった。近距離検出器はまだ利用できず、同時停止時間も短すぎてスペクトルを抽出できなかった。2017年には近距離検出器が稼働し、20.8日間の一回を含む四回の同時停止から17.2日の有効測定時間が得られた。この二つの変化が定量測定を可能にした。

遺産と現在のデータ

この測定は、残留束の総和計算と実測計数を比較する初の公表された定量的基準点を提供する。予測の規格化不確かさは7.4%で、観測された比較は$\Delta=18\pm19$Δ=18±19事象である。この二つを混同してはならない。

応用面で最も直接利用できる新情報は、予測された内訳である。モデルは停止中のサイトからの信号の44%を炉心ではなくプールに帰属させる。炉心に残る燃料だけをモデル化する監視手法は、期待信号を大幅に過小評価することになる。一方、この論文から地表の小型検出器が達成できる感度を定量化することはできない。その種のデータは論文になく、著者らは地下装置の有効体積から外挿していない。

研究者の視点 — 未解決の問い

以下は将来考えられる展開であり、本研究の結果ではない。

• 決定的なスペクトル検証。 統計を増やし、$^{144}\mathrm{Ce}$144Ceと$^{106}\mathrm{Ru}$106Ruの相対的な時間変化を数か月追跡すれば、モデルの時間形状を検証できる。$^{90}\mathrm{Sr}$90Sr–$^{90}\mathrm{Y}$90Y対の支配を観測するには数年を要する。予測が90%を超えるのは十年以上後だからである。
• 幾何配置によってプールと炉心を分ける。 ここでの56%/44%はモデル予測であって、測定による分離ではない。400 mの距離では38 mの差を分解できない。検出器をはるかに近くに置くか、異なる角度に二台置けば、計算で振り分ける代わりに二つの寄与を実験的に分離できる可能性がある。
• 核構造というボトルネック。 支配的な不確かさは検出器ではなく、$^{144}\mathrm{Pr}$144Prのベータ遷移形状に由来する。このベータスペクトルを専用に測定すれば6.0%の項を低減できる可能性があるが、その場合にモデルが到達する最終精度を本研究は定量化していない。

出典

2026年八月16日のファクトチェック監査で確認した参考文献。

• Double Chooz Collaboration, “First Measurement of Neutrino Emissions from Spent Nuclear Fuel by the Double Chooz Experiment”, Physical Review Letters 137, 061803, 2026年八月4日 — 査読済み論文。DOI：10.1103/dr26-j19g
• 同じ研究のプレプリント版：arXiv:2510.04869 — 全文を確認（要旨、手法、系統不確かさの表）。

背景文献

1. P. Vogel、J. F. Beacom — Angular distribution of neutron inverse beta decay — 逆ベータ崩壊の運動学と基準断面積。

信頼性に関する声明

査読済み論文から再現した計数、統計的有意度、エネルギー窓、幾何配置、不確かさ成分には高い信頼を置く。応用範囲への信頼は中程度である。論文は核保障措置に言及するが、地表の小型検出器もインベントリ推定も試験していない。プールに帰属する44%はモデル出力であり、測定による分離ではない。この初の定量解析には、まだ独立した追試がない。一意性確認：登録簿A001–A018に同一のDOI、arXiv識別子、または等価な中心結果はなく、判定はUNIQUE。DOI、日付、数式、リンク、重複がないことを確認済み。