[论文解读] 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的闭式表达式，提供了显式的四阶展开示例，并通过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)及其空间导数表示，从而在粒子加速器中实现高效仿真。
• 通过显式的四阶展开示例展示该方法，并验证其与麦克斯韦方程组的一致性。
• 支持未来设施（如中微子工厂和μ子对撞机）中先进的μ子离子化冷却与束流动力学研究。

提出的方法

• 通过在z=0处对Bz(r,θ,z)进行泰勒级数展开，假设中平面具有反对称性（即在z=0处Br=Bθ=0）。
• 应用真空麦克斯韦方程组∇·B=0和∇×B=0，建立不同z方向上场分量之间的关系，从而递归确定高阶项。
• 将Br、Bθ和Bz表示为Bz(r,θ,0)的径向、角向和z向导数的无穷级数，其中微分算子( ∂²/∂r² + (1/r)∂/∂r + (1/r²)∂²/∂θ² )的n次幂被使用。
• 利用算子展开方法生成至四阶的项，其中Br和Bθ的项包含z的奇次幂，而Bz的项包含z的偶次幂。
• 以确保展开至偶数阶时∇·B=0，展开至奇数阶时∇×B=0的方式实现该展开。
• 使用Mathematica中的符号计算验证该方法，并将其应用于ICOOL中模拟的逆旋回场，以实现μ子追踪。

实验结果

研究问题

• RQ1 如何仅基于z=0中平面上的磁场信息，准确地在圆柱坐标系中将其展开至任意阶次？
• RQ2 在仅已知Bz(r,θ,0)的条件下，Br、Bθ和Bz在z≠0处的解析形式是什么，同时满足麦克斯韦方程组？
• RQ3 展开阶次如何影响模拟中磁场的散度与旋度？
• RQ4 高阶场展开是否能提高逆旋回μ子冷却系统中粒子轨道追踪的精度？
• RQ5 何种场展开算子结构可确保在给定阶次内与∇·B=0和∇×B=0保持一致？

主要发现

• 本文推导出Br、Bθ和Bz的闭式表达式，其以Bz(r,θ,0)及其空间导数表示，适用于z方向的任意阶次。
• 该展开在截断至偶数阶时保证∇·B=0，在截断至奇数阶时保证∇×B=0，从而保持了基本场约束的正确性。
• 提供了显式的四阶展开，表明Br和Bθ是z的奇函数，而Bz是偶函数，与中平面的对称性一致。
• 该方法可在不依赖数值微分或场图的情况下，实现逆旋回场中高精度的粒子追踪。
• 该方法在模拟的ICOOL反旋回场中成功应用，显示μ子轨道在约200 ns内收敛，验证了场模型的物理一致性。
• 推导出的场展开算子( ∂²/∂r² + (1/r)∂/∂r + (1/r²)∂²/∂θ² )^n 允许使用符号计算系统地生成更高阶项。