Simulating Neutron Protein Crystallography Experiments: Applications to the Development of the NMX Instrument at ESS

ESS の NMX 装置向けに McStas を用いたモンテカルロ中性子線追跡シミュレーションを行い、中性子線分割によるイベント形成の改善や新しいサンプリング手法の導入、さらに空気やビームストップ散乱の影響評価を通じて、タンパク質結晶構造解析データの処理可能性を実証しました。

要約

この論文は、**「実際に実験をする前に、コンピューターの中で『未来の超高性能実験』をシミュレーション（模擬実験）して、その実験室を完成させる」**という、とても面白い取り組みについて書かれています。

専門用語を避け、身近な例え話を使って説明しますね。

🧪 1. 何をやっているの？「 neutron（中性子）なし」での実験

通常、タンパク質の構造を調べるには、巨大な加速器を使って「中性子」という目に見えない粒子をタンパク質にぶつけ、その跳ね返り（回折）をカメラで撮ります。でも、この実験は非常に難しく、まだ実験室（ESS の NMX 装置）が完成していない状態です。

そこで著者たちは、**「実際に中性子を飛ばさなくても、コンピューターの中で完璧な物理法則を再現すれば、実験結果を予測できる」と考えました。
まるで、「実際に飛行機を飛ばさなくても、スーパーコンピュータで風や重力を計算して、飛行機の設計図が完璧かどうかをチェックする」**ようなものです。

🎮 2. 使ったツール：「McStas」という超リアルなゲームエンジン

彼らが使ったのは「McStas（マックスタス）」というソフトです。これは中性子の動きを追跡するシミュレーターです。

• 仕組み: 中性子を「光の粒（レイ）」のように考え、それが実験装置の鏡やスリットをどう通り抜け、タンパク質にどうぶつかるかを計算します。
• すごいところ: 普通の計算だと、タンパク質にぶつかる確率は極めて低く、何十億回も計算しないと結果が出ません。そこで彼らは**「分裂（SPLIT）」**という魔法のような技を使いました。
  • 例え話: 1 人の探偵が事件現場に行くのに 1 時間かかるなら、その探偵を 1 万人にコピーして同時に派遣すれば、1 時間で 1 万倍の情報を集められます。これが「SPLIT」です。計算時間を劇的に短縮し、膨大なデータを瞬時に生成しています。

📸 3. データの加工：「確率」を「写真」に変える

シミュレーションの結果は、最初は「この場所に中性子が当たる確率は 0.0001% です」という数字の羅列（確率データ）です。でも、実際の実験では「ピカッ！」という光（イベント）として記録されます。

• 問題: 膨大な確率データを全部メモリに載せて処理するのは、**「図書館の本を全部一度に机の上に広げようとして、机が崩壊する」**ようなものです。
• 解決策: 彼らは新しい方法（ウェイト付きリザーバーサンプリング）を開発しました。
  • 例え話: 川から水を汲み取る際、川全体を一度に汲み上げるのではなく、**「流れてくる水をリアルタイムで、必要な量だけバケツにすくい取る」**技術です。これにより、何テラバイトもの巨大データでも、普通のコンピューターで処理できるようになりました。

🏗️ 4. 実験室の設計図を完成させる

このシミュレーションを使って、彼らは ESS（欧州中性子源）に建設中の「NMX」という実験装置の設計をテストしています。

• 空気の邪魔: 実験室の空気が中性子を散乱させてノイズになるかどうか、シミュレーションで確認しました。
• カメラの位置: 検出器（カメラ）をどこに置けば、タンパク質の構造が最もくっきり写るか、11 通りの配置をシミュレーションで比較しました。
  • 例え話: 暗闇で写真を撮る際、フラッシュの位置やカメラの角度をシミュレーションで試して、「ここなら一番きれいに写る！」と決めてから、実際にカメラを据え付けます。

🎯 5. 結果：本物そっくりのデータができた

シミュレーションで得られたデータを、実際の解析ソフト（DIALS）にかけてみたところ、**「本物の実験で得られたデータと見分けがつかない」ほど良い結果が出ました。
これにより、実験室が完成する前に、「この装置なら、どんなタンパク質の構造も解明できる！」**という自信が持てたのです。

まとめ

この論文は、**「まだできていない実験装置を、コンピューターの中で完璧に作り上げ、その性能を事前に証明した」**という画期的な成果です。

• 現実: 実験室が完成するまで待たないとわからない。
• この研究: コンピューターの中で実験を何万回も繰り返し、完成した後の未来を先取りして、装置の設計を完璧にした。

これによって、将来、実際に実験が始まった瞬間から、世界中の科学者がすぐに新しいタンパク質の構造を解明できるようになるのです。まるで、**「未来の地図を事前に描き終えて、旅路のトラブルをすべて回避した」**ようなものです。

技術概要

以下は、提示された論文「Simulating Neutron Protein Crystallography: Experiments: Applications to the Development of the NMX Instrument at ESS」の技術的な要約です。

論文の概要

本論文は、スウェーデンの欧州中性子源（ESS）で建設中のマクロ分子回折装置「NMX」の commissioning（試運転・調整）および開発を支援するために、モンテカルロ法を用いた中性子タンパク質結晶回折実験のシミュレーション手法を確立したことを報告しています。実在の中性子源を使用せずに、高度な物理モデルと計算機シミュレーションによって実験データを生成・解析する手法を確立し、NMX 装置の性能評価やデータ収集戦略の策定に貢献しています。

---

1. 背景と課題 (Problem)

• 中性子マクロ分子結晶学（n-MX）の課題: 水素原子の位置を高精度に決定できる n-MX は X 線や電子顕微鏡にとって重要な補完技術ですが、実験には巨大な結晶（〜1 mm³）と重水素化（deuteration）が必須であり、実験のハードルが高いです。
• NMX 装置の特殊性: ESS の NMX 装置は、長いパルス幅と明るいモデレーター、可変幾何構造の検出器システム（3 つの移動式 GEM パネル）を特徴としており、既存の装置とは異なる設計です。
• シミュレーションの困難さ:
  • 逆格子空間（reciprocal space）を網羅的にサンプリングするには、10⁶ 以上の反射を扱う必要があり、計算コストが膨大になります。
  • 中性子散乱確率は極めて低いため、統計的な信頼性を得るには莫大な数のシミュレーションが必要ですが、メモリと計算時間の制約が壁となります。
  • 建設中の装置であるため、実測データとの比較ができず、シミュレーションデータが現実のデータ処理フロー（DIALS など）と互換性があるか検証する必要があります。

2. 手法 (Methodology)

本研究では、中性子光学シミュレーションソフトウェア McStas を基盤とした以下のワークフローを開発・適用しました。

• 装置モデルの構築:
  • ESS の「バタフライ」モデレーター、ビームライン、チョッパー、コリメーター、試料（1 mm³ のルードレキシン結晶）、および可動式検出器パネルを詳細にモデル化しました。
  • 空気散乱やビームストップ（6LiF 粉末入りアルミ管）からの散乱など、現実的なノイズ源も Union コンポーネントを用いてシミュレーションに組み込みました。
• 計算効率化の工夫（SPLIT 指令）:
  • 結晶からの散乱確率が低いため、結晶に到達した中性子線を多数の「スプリット（コピー）」に分けて処理する McStas の SPLIT 指令を適用しました。これにより、結晶までの軌道計算を省略し、反射リストのキャッシュを活用することで、計算効率を劇的に向上させました（例：SPLIT 11000 でイベント数が約 1000 倍に増加）。
• イベント生成とサンプリング（Single-pass Weighted Reservoir Sampling）:
  • McStas は確率（event probability）を出力しますが、実際の検出器データ（イベント数）に変換する必要があります。
  • 膨大な確率リストをメモリにすべて読み込むことが不可能なため、単一パス重み付きリザーバーサンプリング（single-pass weighted reservoir sampling） アルゴリズムを採用しました。これにより、ストリーミングデータとしてリアルタイムにイベントを生成し、NeXus 形式（TOFRaw）の実際の検出器データと同等の形式に変換しました。
• データ処理パイプライン:
  • 生成されたシミュレーションデータを、実データ処理で用いられる DIALS ソフトウェア（spot finding, indexing, integration）にそのまま入力し、処理の妥当性を検証しました。

3. 主要な貢献 (Key Contributions)

• 現実的なシミュレーションデータの生成: 確率ベースのシミュレーションを、実際の検出器で得られるような「イベントデータ」に変換する完全なパイプラインを確立しました。
• 計算効率の飛躍的向上: SPLIT 指令と単一パスサンプリング法の組み合わせにより、逆格子空間の網羅的なサンプリングを現実的な計算時間で可能にしました。
• NMX 装置設計への応用: 可変検出器配置（11 種類の設定）を一度のシミュレーションで評価可能とし、検出器配置がデータ品質（特に空気散乱の影響）に与える影響を定量的に評価しました。
• オープンソースツールの開発: 構造因子計算、ストリーミングサンプリング、データ削減（ESSNMX）などのツールを公開し、コミュニティへの貢献を行いました。

4. 結果 (Results)

• SPLIT 手法の有効性: SPLIT 値を 10000 に設定することで、検出器表面のほぼ完全なサンプリングが可能となり、計算時間は 10 倍の増加に対してイベント数は 10 倍増加し、効率的な計算が実現されました（Table 1, 2）。
• DIALS による解析: シミュレーションデータを DIALS で処理した結果、スポット発見、インデキシング、積分が正常に行われ、実験データと同様の統計的性質（$I/\epsilon(I)$ など）を示しました（Table 3）。
• ノイズ源の評価: 空気散乱が背景ノイズの主要因であることが確認されました。特に検出器が試料から遠い配置（Configuration 9, 10）では、空気散乱の影響が支配的になり、データ品質が低下することが示されました。
• ビームストップの影響: 現実的なビームストップ（6LiF）を含むシミュレーションでは、ビームストップ自体からの散乱は比較的小さいものの、空気散乱による拡散散乱が $d \approx 1.5$ Å 付近で顕著であることが分かりました。

5. 意義と将来展望 (Significance)

• 装置 commissioning の事前検証: NMX 装置が稼働する前に、データ収集戦略や検出器配置の最適化、解析パラメータの調整をシミュレーションデータを用いて行えるため、実運用時のリスクを低減できます。
• 実験誤差のモデル化: 将来的には、より複雑なノイズ源（H₂O の存在、電子ノイズなど）をシミュレーションに組み込むことで、実験誤差の発生源を特定し、解析アルゴリズムの改善に役立てることが期待されます。
• 計算科学と実験科学の融合: 本論文は、計算機シミュレーションが単なる設計支援ツールではなく、実験データそのものを生成・解析するプラットフォームとして機能しうることを示しました。

総じて、本研究は ESS の NMX 装置の成功に向けた重要な基盤技術を提供し、中性子結晶学の分野において「計算による実験」の新たな標準を確立したと言えます。