[論文レビュー] OPAL a Versatile Tool for Charged Particle Accelerator Simulations

Andreas Adelmann, Pedro Calvo+10 more|arXiv (Cornell University)|May 16, 2019

Particle accelerators and beam dynamics参考文献 78被引用数 23

ひとこと要約

OPAL は、3次元空間電荷、外部場、および粒子-物質相互作用を伴うビームダイナミクスの高精度なモデル化を可能にする、オープンソースで並列処理可能かつオブジェクト指向のシミュレーションフレームワークである。線形加速器（OPAL-t）とサイクロトロン/FFA（OPAL-cycl）の両方をサポートし、ベンチマークと実際の応用（例えば PSI Injector II や IsoDAR施設）を通じて、バーチカル運動やマルチターンビームダイナミクスといった複雑な現象のシミュレーションにおいて高い正確性を達成している。

ABSTRACT

Many sophisticated computer models have been developed to understand the behaviour of particle accelerators. Even these complex models often do not describe the measured data. Interactions of the beam with external fields, other particles in the same beam and the beam walls all present modelling challenges. These can be challenging to model correctly even with modern supercomputers. This paper describes OPAL (Object Oriented Parallel Accelerator Library), a parallel open source tool for charged-particle optics in linear accelerators and rings, including 3D space charge. OPAL is built from the ground up as a parallel application exemplifying the fact that high performance computing is the third leg of science, complementing theory and experiment. Using the MAD language with extensions, OPAL can run on a laptop as well as on the largest high performance computing systems. The OPAL framework makes it easy to add new features in the form of new C++ classes, enabling the modelling of many physics processes and field types. OPAL comes in two flavours: OPAL-cycl: tracks particles with 3D space charge including neighbouring turns in cyclotrons and FFAs with time as the independent variable. OPAL-t: models beam lines, linacs, rf-photo injectors and complete XFELs excluding the undulator. The code is managed through the git distributed version control system. A suite of unit tests have been developed for various parts of OPAL, validating each part of the code independently. System tests validate the overall integration of different elements.

研究の動機と目的

• 既存のモデルに見られる制限を克服し、高精度な荷電粒子加速器設計のためのスケーラブルでオープンソースのシミュレーションツールを開発すること。
• 3次元空間電荷、複数回のクーロン散乱、残留ガス相互作用を含む複雑なビームダイナミクスの正確なモデル化を可能にすること。
• 統一的かつ拡張可能なソフトウェアアーキテクチャにより、線形加速器（OPAL-t）とサイクロトロン/FFA（OPAL-cycl）の両方をサポートすること。
• 高性能コンピューティング（HPC）および新興アーキテクチャ（例：GPU）を加速器シミュレーションワークフローに統合すること。
• 主要な加速器プロジェクトにおけるベンチマークと実際の応用を通じてコードの妥当性を検証すること。

提案手法

• OPAL は、自己場および外部場を解くために静電的粒子-セルフ（PIC）法を用い、3次元空間電荷計算を可能にする。
• 入力設定には MAD 言語を拡張して使用し、ラップトップから大規模 HPC システムまで、さまざまな環境でシミュレーションを実行可能である。
• 複数の空間電荷ソルバ（FFTベース、SAAMG、AMR（適応メッシュリファインメント））をサポートしており、精度と性能のニーズに応じた選択が可能である。
• OPAL-cycl は、複数ターンにわたるビームダイナミクスをモデル化するため、時間の独立変数として用い、隣接するビームバンチのカップリングも含む。
• OPAL-t は、ビームライン、RFフォトインジェクタ、および完全なXFELを含むモデルを提供し、シンクロトロン放射とビーム残渣相互作用も考慮している。
• コードは真のオブジェクト指向 C++ アーキテクチャで構築されており、新しい物理モデルや場の種別をモジュラーに拡張可能である。

実験結果

リサーチクエスチョン

• RQ1 線形加速器およびサイクロトロンの両方において、高精度な3次元空間電荷効果をどのように正確にシミュレートできるか。
• RQ2 OPAL は、サイクロトロンにおけるバーチカル運動のような複雑なビームダイナミクスをどの程度正確に再現できるか。
• RQ3 OPAL は、RFフォトインジェクタから完全なXFELに至る多様な加速器構成において、一貫性があり妥当性が確認された結果を達成できるか。
• RQ4 サイクロトロンにおけるマルチターンおよび隣接バンチ効果の取り込みが、ビーム損失およびエミッタンス成長にどのように影響するか。
• RQ5 OPAL は、デスクトップ環境および大規模 HPC システムの両方で、性能とスケーラビリティに優れているか。

主な発見

• OPAL は、PSI Injector II サイクロトロンにおけるバーチカル運動を正確に再現し、最大 2.7 mA の H₂⁺ ビームを制御された1ターン抽出に成功した。
• シミュレーションでは、初期に不一致があるビームであっても、DIC サイクロトロン内で10ターン以内に安定でほぼ円形の分布が形成されることを示した。
• IsoDAR プロジェクトにおいて、OPAL シミュレーションは、6.65 mA のビーム電流がコリメーションとビームクリーニングを経て 5.1 mA の抽出電流に達することを示した。
• OPAL-cycl シミュレーションにより、DSRC において220ターン以降に垂直ビームサイズが最大30 mmまで増大するが、利用可能なサイクロトロン断面内では管理可能であることが確認された。
• OPAL の3次元空間電荷ソルバ（FFT、SAAMG、AMR）は、解析的モデルおよび実験データと比較され、良好な一致を示した。
• OPAL の粒子-物質相互作用モデルは、マルチスラブ構成におけるエネルギー損失と複数回のクーロン散乱を正確に再現した。