[論文レビュー] Magnetic Field Expansion Out of a Plane: Application to Inverse Cyclotron Muon Cooling

T. Hart, Summers, D. J.+1 more|arXiv (Cornell University)|May 13, 2011

Particle accelerators and beam dynamics参考文献 18被引用数 3

ひとこと要約

本稿では、円筒座標系における中間面からの磁場の高次解析的展開を提示し、逆サイクロトロン・ミューオン冷却における荷電粒子軌道の正確なモデル化を可能にする。真空中のマクスウェル方程式と任意の次数のテイラー展開を活用することで、指定された Bz(r,θ,0) から Br, Bθ, Bz の閉形式表現を導出する。4次までの明示的例と Mathematica を用いた検証を実施。この手法により、選択された次数まで一貫した発散なし・渦なしの磁場が保証され、将来のミューオン衝突加速器やニュートリノファクトリーの研究開発における正確なビームダイナミクスシミュレーションに不可欠である。

ABSTRACT

In studies of the dynamics of charged particles in a cyclotron magnetic field, the specified field is generally B_z in the z = 0 midplane where B_r and B_theta are zero. B_r(r, theta, z) and B_theta(r, theta, z) are simple to determine through a linear expansion which assumes that B_z is independent of z. But, an expansion to only first order may not be sufficient for orbit simulations at large z. This paper reviews the expansion of a specified B_z(r, theta, z = 0) out of the z = 0 midplane to arbitrary order, and shows simple examples worked out to 4th order.

研究の動機と目的

• 円筒座標系において中間面 z=0 からの磁場を任意次数にまで系統的に展開する手法を開発すること。
• 展開された磁場が選択された次数まで ∇·B=0 および ∇×B=0 を満たすように保証すること。これは正確な粒子追跡に不可欠である。
• Br, Bθ, Bz に対して、Bz(r,θ,0) 及びその空間微分を用いた閉形式表現を提供することにより、粒子加速器内での効率的なシミュレーションを可能にすること。
• 明示的な4次展開を提示し、マクスウェル方程式との整合性を検証すること。
• 将来のニュートリノファクトリーおよびミューオン衝突加速器のための高度なミューオン電離冷却およびビームダイナミクス研究を支援すること。

提案手法

• Bz(r,θ,z) を z=0 の周りにテイラー級数展開し、中間面における反対称性（z=0 で Br=Bθ=0）を仮定する。
• 真空中のマクスウェル方程式 ∇·B=0 および ∇×B=0 を適用し、z 方向における場の成分を関連づけることで、高次項を再帰的に決定する。
• Br, Bθ, Bz を、Bz(r,θ,0) の径方向、角度方向、z 方向微分を含む無限級数として表現する。この際、微分作用素 ( ∂²/∂r² + (1/r)∂/∂r + (1/r²)∂²/∂θ² ) を n 乗した形を用いる。
• 作用素展開を用いて、4次までの項を生成し、Br および Bθ には z の奇数次、Bz には偶数次を割り当てる。
• 展開が偶数次で切られた場合に ∇·B=0、奇数次で切られた場合に ∇×B=0 となるように実装する。
• シンボリック計算ソフトウェア Mathematica を用いて手法を検証し、ICOOL における逆サイクロトロン磁場のシミュレーションに応用してミューオン追跡を実施した。

実験結果

リサーチクエスチョン

• RQ1 z=0 の中間面上でのみ指定された磁場を、円筒座標系において任意次数に正確に展開する方法は何か？
• RQ2 Bz(r,θ,0) のみが与えられた場合、z≠0 における Br, Bθ, Bz の解析的形は何か？ ただしマクスウェル方程式を満たすものとする。
• RQ3 展開次数がシミュレーションにおける磁場の発散および渦の性質に与える影響は何か？
• RQ4 高次磁場展開は、逆サイクロトロン・ミューオン冷却系における粒子軌道追跡の精度を向上させ得るか？
• RQ5 指定された次数まで ∇·B=0 および ∇×B=0 を満たすように保証する磁場展開作用素の構造は何か？

主な発見

• 本稿では、Bz(r,θ,0) 及びその空間微分を用いた Br, Bθ, Bz の閉形式表現を、z に関して任意次数まで導出している。
• 展開を偶数次で切り捨てた場合に ∇·B=0、奇数次で切り捨てた場合に ∇×B=0 が保証され、基本的な場の制約が維持される。
• 明示的な4次展開を提示し、Br および Bθ が z に関して奇関数であり、Bz が偶関数であることが確認された。これは中間面対称性と整合的である。
• 本手法により、数値微分や場のマップに依存せずに、逆サイクロトロン磁場内での高精度な粒子追跡が可能となった。
• シミュレーテッド・ICOOL反サイクロトロンに応用した結果、ミューオン軌道が約200 ns 以内に収束し、磁場モデルの物理的整合性が検証された。
• 導出された磁場展開作用素 ( ∂²/∂r² + (1/r)∂/∂r + (1/r²)∂²/∂θ² )^n を用いることで、シンボリック計算を用いて高次項を体系的に生成できるようになった。