Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access

Symplectic tracking through curved three-dimensional fields by a method of generating functions

Jie Li1, Kedong Wang1,2,*, Jinloong Lee1, Kai Wang1, Xu Zhang1, Xueqing Yan1,2, and Kun Zhu1,2,†

  • *Contact author: wangkd@pku.edu.cn
  • †Contact author: zhukun@pku.edu.cn

Phys. Rev. Accel. Beams 28, 094001 – Published 2 September, 2025

DOI: https://doi.org/10.1103/sgl1-17y8

Abstract

Symplectic integrators play a pivotal role in the long-term tracking of charged particles within accelerators. To obtain symplectic maps in accurate simulations of single-particle trajectories, two key components are addressed: precise analytical expressions for arbitrary electromagnetic fields and a robust treatment of the equations of motion. In a source-free region, electromagnetic fields can be decomposed into harmonic functions, applicable to both scalar and vector potentials, encompassing both straight and curved reference trajectories. These harmonics are constructed according to boundary surface field data due to the uniqueness theorem. By finding generating functions to satisfy the Hamilton-Jacobi equation via a perturbative ansatz, we derive maps that are symplectic at any order of the expansion. This method yields a direct mapping from initial to final coordinates, bypassing the tracking of intermediate coordinates. Our developed particle-tracking algorithm translates the field harmonics into beam trajectories, offering a high-resolution integration method in accelerator simulations.

View figure in article

Physics Subject Headings (PhySH)

Article Text

References (30)

  1. W. Wan, L. Brouwer, S. Caspi, S. Prestemon, A. Gerbershagen, J. M. Schippers, and D. Robin, Alternating-gradient canted cosine theta superconducting magnets for future compact proton gantries, Phys. Rev. ST Accel. Beams 18, 103501 (2015).
  2. Y. Giboudot and A. Wolski, Use of transfer maps for modeling beam dynamics in a nonscaling fixed-field alternating-gradient accelerator, Phys. Rev. ST Accel. Beams 15, 044001 (2012).
  3. A. Wolski, J. Gratus, and R.Tucker, Symplectic integrator for s-dependent static magnetic fields based on mixed-variable generating functions, J. Instrum. 7, P04013 (2012).
  4. B. Erdélyi and M. Berz, Optimal symplectic approximation of Hamiltonian flows, Phys. Rev. Lett. 87, 114302 (2001).
  5. Y. K. Wu, E. Forest, and D. S. Robin, Explicit symplectic integrator for s-dependent static magnetic field, Phys. Rev. E 68, 046502 (2003).
  6. A. Wolski and A. T. Herrod, Explicit symplectic integrator for particle tracking in s-dependent static electric and magnetic fields with curved reference trajectory, Phys. Rev. Accel. Beams 21, 084001 (2018).
  7. K. Skoufaris, Non-linear dynamics modelling in accelerators with the use of symplectic integrators, Ph.D. thesis, University of Crete, 2021, 10.12681/eadd/50152.
  8. M. Tao, Explicit symplectic approximation of nonseparable Hamiltonians: Algorithm and long time performance, Phys. Rev. E 94, 043303 (2016).
  9. J. Qiang, Symplectic multiparticle tracking model for self-consistent space-charge simulation, Phys. Rev. Accel. Beams 20, 014203 (2017).
  10. P. J. Channell and C. Scovel, Symplectic integration of Hamiltonian systems, Nonlinearity 3, 231 (1990).
  11. H. Lustfeld, G. Wüstefeld, P. Meads, U. Funk, and L. Schepers, A proper canonical algebraic mapping transformation, in Proceedings of the 12th International Linac Conference, LINAC-1984, Seeheim, Germany (1984), pp. 130–131, https://accelconf.web.cern.ch/l84/papers/tup0016.pdf.
  12. W. Gabella, R. Ruth, and R. Warnock, Beam dynamics with the Hamilton-Jacobi equation, in Proceedings of the 1989 IEEE Particle Accelerator Conference: Accelerator Science and Technology (1989), Vol. 2, pp. 1310–1312, 10.1109/PAC.1989.73432.
  13. J. Bahrdt and G. Wüstefeld, Symplectic tracking and compensation of dynamic field integrals in complex undulator structures, Phys. Rev. ST Accel. Beams 14, 040703 (2011).
  14. M. Titze, J. Bahrdt, and G. Wüstefeld, Symplectic tracking through straight three dimensional fields by a method of generating functions, Phys. Rev. Accel. Beams 19, 014001 (2016).
  15. D. Veres, T. Vaszary, E. Benedetto, and D. Barna, A new algorithm for optimizing the field quality of curved CCT magnets, IEEE Trans. Appl. Supercond. 32, 1 (2022).
  16. P. Schnizer, E. Fischer, and B. Schnizer, Cylindrical circular and elliptical, toroidal circular and elliptical multipoles fields, potentials and their measurement for accelerator magnets, arXiv:1410.8090.
  17. B. van Milligen and A. Lopez Fraguas, Expansion of vacuum magnetic fields in toroidal harmonics, Comput. Phys. Commun. 81, 74 (1994).
  18. C. E. Mitchell and A. J. Dragt, Accurate transfer maps for realistic beam-line elements: Straight elements, Phys. Rev. ST Accel. Beams 13, 064001 (2010).
  19. V. Punzo, S. Besio, S. Pittaluga, and A. Trequattrini, Solution of Laplace equation on non axially symmetrical volumes, IEEE Trans. Appl. Supercond. 16, 1815 (1990).
  20. A. Kildishev, Application of spheroidal functions in magnetostatics, IEEE Trans. Magn. 40, 846 (2004).
  21. A. A. Wolski, Beam Dynamics in High Energy Particle Accelerators (World Scientific Publishing Co. Pte. Ltd., Singapore, 1984), https://worldscientific.com/doi/pdf/10.1142/p899.
  22. J. Li, K. Wang, M. Easton, W. Huang, K. Wang, S. Cai, T. Zhu, C. Zhang, X. Jiang, Y. Liang, W. Wu, C. Lin, X. Yan, and K. Zhu, Beamline analysis for a laser-driven proton therapy accelerator using superconducting bends, Nucl. Eng. Technol. 57, 103267 (2025).
  23. L. Gambini, Expansion of the magnetic flux density field in toroidal harmonics, Master’s thesis, Università di Bologna, 2018, https://cds.cern.ch/record/2316227.
  24. J. Li, K. Wang, K. Wang, X. Yan, and K. Zhu, Evaluation of field quality for curved magnets, IEEE Trans. Appl. Supercond. 34, 1 (2024).
  25. L. Bojtàr, Efficient evaluation of arbitrary static electromagnetic fields with applications for symplectic particle tracking, Nucl. Instrum. Methods Phys. Res., Sect. A 948, 162841 (2019).
  26. L. Bojtár, Frequency analysis and dynamic aperture studies in a low energy antiproton ring with realistic 3d magnetic fields, Phys. Rev. Accel. Beams 23, 104002 (2020).
  27. J. Li, K. Wang, X. Zhang, K. Zhang, X. Yan, and K. Zhu, Using magnetic moment for field reconstruction in accelerator magnets and particle tracking, IEEE Trans. Magn. 61, 1 (2025).
  28. A. Simona, Numerical methods for the simulation of particle motion in electromagnetic fields, Ph.D. dissertation, Politecnico di Milano, Technische Universität Darmstadt, 2020, 10.25534/tuprints-00011687.
  29. S. Manikonda and M. Berz, An accurate high-order method to solve the Helmholtz boundary value problem for the 3D Laplace equation, Int. J. Pure Appl. Math. 23, 365 (2005).
  30. R. Lindberg and M. Borland, Fringe field maps for symplectic models of general Cartesian dipoles, Phys. Rev. Accel. Beams 26, 114001 (2023).

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation