Development of a comprehensive PMT optical model for the JUNO experiment

技術要約：JUNO実験のための包括的なPMT光学モデルの開発

問題提起
江門地下ニュートリノ観測装置（JUNO）は、ニュートリノを検出するために17,612本の20インチ光電子増倍管（PMT）を使用している。精密な検出器シミュレーションは、ニュートリノの質量順序の決定および振動パラメータの測定において極めて重要である。しかし、既存の光学モデルは、同型のPMTは同一の光電陰極（PC）および反射防止膜（ARC）特性を持つという簡略化された仮定に依存していた。本研究は、PMT表面における光子検出効率（PDE）の不均一性と、個々のPMT間におけるPDEおよび反射率の著しい変動を考慮することで、こうしたモデルの限界に対処するものである。これらの要因は、JUNOの動作環境（水浸）および入射光子の入射角が、質量測定時の空気中での条件とは大きく異なるため、簡略化されたモデルでは正確に捉えられないイベント依存のPDE変動を理解する上で極めて重要である。

手法
著者らは、大量生産時のテストからの経験的データおよび新しい反射率測定値を取り入れるために、以前のフレームワーク[6]を拡張し、包括的なPMT光学モデルを開発した。その手法は以下のステップで構成される：

1. データの統合： モデルは、2つのソースからのPDEデータを利用する。「コンテナシステム」（全17,612本のPMTの平均PDEを測定）と「スキャニングステーションシステム」（約3,200本のPMTのサブセットに対し、7つの特定の天頂角で測定されたPDEの均一性）である。さらに、積分球を用いたセットアップを用いて、4つの天頂角において669本のPMT（HPK 276本、NNVT normal-QE 115本、およびNNVT high-QE 174本）の新しい反射率データを収集した。
2. モデルの拡張： モデルは、同型のPMTについては屈折率（$n$）および消衰係数（$k$）のスペクトル応答は同一であるが、PCおよびARCの厚さが異なることで不均一性が生じるという仮定に基づいている。厚さの分布は、経験的な二次式を用いて天頂角（$\theta$）の関数としてモデル化されている。
3. パラメータ抽出： $\chi^2$最小化手順により、測定された反射率およびPDEデータを同時にフィッティングし、各PMTのPCおよびARCの厚さマップを制約する。このプロセスにより、量子効率（QE）と集光効率（CE）を測定されたPDEから分離する。
4. シミュレーションと補正： Geant4ベースの光学シミュレーションを用い、直接的な光子の吸収と、PMT構造内での透過光子の吸収を区別した。これにより、入射位置における真の吸収率を特定するためのPDEデータの補正が可能となった。
5. スペクトルの拡張： 分散関係と文献データを用い、光学特性（$n$, $k$, および脱出因子）を元の390–500 nmの範囲から、JUNOの全スペクトルである300–700 nmへと拡張した。

主な貢献

• 厚さのマッピング： 本研究は、全17,612本の20インチPMTに対するPCおよびARCの厚さマップの導出に成功し、これまで考慮されていなかった不均一性を明らかにした。
• 集光効率（CE）の決定： PDEデータを分解し、内部の光学プロセスを考慮することで、著者らは異なるPMTタイプにおける天頂角の関数としてのCEを決定した。
• 包括的な光学パラメータ： 屈折率、消衰係数、および脱出因子を、3つの異なるPMTタイプ（HPKダイノードPMT、NNVT normal-QE PMT、およびNNVT high-QE PMT）について300–700 nmの範囲で評価した。
• 反射率の特性評価： 新しい高精度な反射率測定を確立し、表面曲率の影響を補正した上で、スキャニングステーションでカバーされていないPMTのためのデータセットを提供した。

結果

• 反射率の変動： 包括的なモデルは、簡略化された一様モデルと比較して、著しく異なる反射率を予測する。HPK PMTについては反射率が高くなると予測される一方で、NNVT PMT（normalおよびhigh-QEの両方）については、反射率が約30%から40%低くなることが予測される。
• PDEへの影響： 平均PDEの変化はより緩やかであり、小さな入射角（AOI）では数パーセント程度にとどまるが、大きなAOIでは5%から10%に達する。
• 厚さ分布： フィットされた厚さマップは、異なるPMTタイプにおいて天頂角に対する明確な依存性を示している。例えば、NNVT normal-QE PMTのARC厚さは反射率および吸収率に与える影響が無視できるほど小さかったため、モデル内では固定値として扱われた。
• CEの挙動： 計算されたCEは、一般に天頂角が増加するにつれて減少する。一部で1を超える値が観察されたが、これは同型のPMT間で一定の脱出係数スペクトルを仮定していることに起因しており、利用可能なデータではCEと脱出係数の変動を完全に分離できないことを示している。

意義
本論文は、この包括的なPMT光学モデルが、JUNO検出器シミュレーションの精度向上およびエネルギー応答モデルの精緻化に不可欠であると主張している。一様な光学特性という仮定を超えて、このモデルは光が検出器構成要素とどのように相互作用するかをより現実的に表現する。著者らは、このモデルがJUNOシミュレーションプログラムに統合され、エネルギー分解能、不均一性、および非線形性の評価に用いられる予定であると述べている。さらに、この包括的なアプローチは、精密な光学モデリングを通じて検出器性能の向上を目指す他のPMTベースのアプリケーションにとっても、価値あるリファレンスを提供するものである。本研究は、PDEの変化は中程度であるものの、反射率の予測における大幅な変化があるため、高精度な物理解析にはこの洗練されたモデルが必要であると結論付けている。