Radiation of a particle performing helical motion in a multilayer cylindrical waveguide

Armen Grigoryan; Bagrat Grigoryan; Francois Lemery; Hrach Babujyan; Klaus Floettmann; Lusine Aslyan; Mikayel Ivanyan; Suren Arutunian; Vardan Avagyan

REVIEW 4 major objections 6 minor 16 references

Radiation of a particle performing helical motion in a multilayer cylindrical waveguide

T0 review · 4 major / 6 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A new algorithm computes the radiation of a helicoidally moving particle inside an arbitrary multilayer cylindrical waveguide by stitching the free-space helical solution to the wall modes.

desk verdict A useful modal-matching algorithm for helical-particle radiation in multilayer waveguides, but the imported particular solution and a band-edge typo need attention before I'd rely on the numbers. read the letter →

arxiv 2501.12802 v2 pith:UCGXFOYC submitted 2025-01-22 physics.acc-ph

Mikayel Ivanyan , Bagrat Grigoryan , Armen Grigoryan , Lusine Aslyan , Vardan Avagyan , Hrach Babujyan , Suren Arutunian , Klaus Floettmann

show 1 more author

Francois Lemery

This is my paper · ORCID

classification physics.acc-ph PACS 41.60.-m41.20.Jb

keywords helicalmotionradiationcylindricalwaveguidemultilayerwallpartialregionsmethodfree-spaceparticularsolutiondispersionmatrixundulatorresonantfrequencies

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper claims that the radiation field of a point charge spiraling along the axis of an infinite cylindrical waveguide with an arbitrary multilayer side wall can be computed by a direct algorithm. The algorithm combines the recently obtained exact solution for helical motion in free space, used as the particular solution of Maxwell's equations inside the cavity, with the standard partial-regions matching of fields across the cylindrical layers. The result is a four-by-four linear system $\hat{D}\,\hat{X}=\hat{S}$; solving it gives the amplitudes of the waveguide modes and hence the full radiation field. If the algorithm is correct, it turns a previously partial problem into a routine calculation for any number and type of layers, which matters for designing narrow-band radiation sources based on helical undulators.

What carries the argument

The two load-bearing objects are the multilayer dispersion matrix $\hat{D} = \hat{Q} \hat{W}_H - \hat{W}_J$, built from products of layer transfer matrices $\hat{Q}_i$ whose elements are combinations of Bessel and Hankel functions, and the free-space helical solution of reference [7] used as a particular solution. The transfer matrices encode the field matching across each cylindrical layer boundary; the identity $U_{12}U_{21}-U_{11}U_{22} = -4/(\pi^2 \nu_i^2 a_i a_{i+1})$ keeps the arithmetic under control.

What would settle it

Take the single-layer copper waveguide parameters of Section 7 ($v=0.99c$, $v_z=0.98c$, $l=5\,\text{cm}$, $a_1=1\,\text{cm}$) and simulate the radiation field with an independent time-domain Maxwell solver; the predicted resonant frequencies and the amplitudes of the radial electric field component at $r=0.5a_1$ should match the paper's figures. A mismatch in the location or shape of the resonances would indicate the particular solution or the matrix assembly is wrong.

Watch

Extended reading notes

Core claim

The central claim, in the paper's own terms, is that the complete solution of the inhomogeneous Maxwell equations for a point charge in helical motion inside an infinite multilayer cylindrical waveguide is obtained as $\hat{X} = \hat{D}^{-1} \hat{S}$, where $\hat{D}(k, p_m, \nu_{m,0})$ is the four-by-four multilayer dispersion matrix constructed by the partial-regions method and $\hat{S}$ is the vector of tangential field components of the free-space helical radiation solution of reference [7] evaluated at the inner wall radius $a_1$. The radiation field in the vacuum cavity is the sum of the homogeneous waveguide modes with these amplitudes plus the free-space particular solution itself. The resonant frequencies are given by $\det \hat{D}=0$ and are independent of the spiral radius, while the field amplitudes carry the spiral dependence through $\hat{S}$.

Load-bearing premise

The free-space helical radiation solution imported from reference [7] is assumed to be the exact, complete particular solution inside the cavity; if it is not, the computed mode amplitudes are wrong even though the resonance frequencies from $\det \hat{D}=0$ would remain unchanged.

Editorial extensions

If this is right

Resonant frequencies of the radiation are determined solely by the waveguide structure, via $\det\hat{D}=0$, and do not depend on the spiral radius; the field amplitudes at those resonances are set by the free-space particular solution.
For single-layer resistive waveguides the spectral lines acquire finite width, and the high-frequency (forward-radiated) branch is more sensitive to wall resistivity and to added layers.
Adding a thin dielectric layer inside a resistive waveguide leaves the number and location of resonances nearly unchanged but can substantially reduce the attenuation decrement of forward radiation.
An NEG coating intended to maintain vacuum can strongly distort forward radiation at higher conductivity, so its thickness must be kept small.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Although the paper does not state it, the same $\hat{X}=\hat{D}^{-1}\hat{S}$ template would work for any source whose free-space multipole expansion is known: only the driving vector $\hat{S}$ changes.
The authors note the phase velocities of both branches are synchronous with $v_z$ but leave the proof out; verifying this equality on the computed eigenvalues would be a direct, low-cost consistency check.
The method's reliance on an imported free-space solution suggests that its accuracy for a given layer stack can be tested against a finite-difference time-domain simulation before committing to a physical design.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Two things to know before you read this. First, it is a real advance over the prior conference paper [6]: it removes the indefinite function in that solution by importing the exact free-space helical solution from the authors' own [7] and using it as the particular solution. The explicit amplitude formulas for one- and two-layer walls (26)-(31) are new and practically useful, and the numerical examples for resistive, copper-dielectric, and three-layer NEG-coated waveguides show the expected narrow resonances. The transfer-matrix matching is standard, and the paper is honest about what it does and does not derive.

Second, the load-bearing particular solution from [7] is not derived or validated anywhere in this paper. That is the main soft spot. In (21)-(23) the source vector S is assembled from [7]'s field components, so any error in the coefficients (20) propagates directly into the computed field amplitudes. The resonance equation D=0 is independent of S, so the resonant frequencies can be right even when the fields are wrong—an easy failure mode to miss. The paper does not check S against Maxwell's equations or against any known limit, and the only external reference is a self-cited arXiv preprint. I would want that fixed before trusting the numbers.

There is also a concrete error in Eq. (32). The allowed band from Re{nu_m,0}>0, Im{nu_m,0}=0 is m*omega0/(1 ± v_z/c), not m*omega0/(1 ± v_z^2/c^2). The printed band is too narrow by a factor (1 ± beta^2)/(1 ± beta). This looks like a typo rather than a deep flaw, but it should be corrected. Minor issues: the figure captions for Figs. 5 and 6 list permittivity as "10 microns", the numerical examples have no quantitative baseline comparison, and the claimed mutual cancellation of TM/TE bursts at the band edges is only demonstrated graphically.

Who should read this: accelerator and THz source designers who want to optimize layer thickness, permittivity, and NEG coatings for narrow-band helical undulator radiation. They will get a workable algorithm and explicit formulas to build on.

My recommendation: send it to peer review. The method is sound and the contribution is meaningful, but the authors need to either derive or independently cross-check the particular solution from [7], fix Eq. (32) and the typos, and add at least one quantitative validation against a known limit, such as an ideally conducting waveguide or a resistive-wall case.

Referee Report

4 major / 6 minor

Summary. The paper presents an analytical algorithm for computing the radiation field of a point charge moving along a helical trajectory on the axis of an infinite cylindrical waveguide with an arbitrary number of concentric material layers. The homogeneous solution is obtained by standard radial mode matching, leading to the dispersion matrix D, and the particular solution is taken from the authors' previous free-space result [7]. The final amplitudes are given as X = D^{-1} S, with explicit formulas for single- and double-layer walls and numerical examples for resistive, metal-dielectric, and NEG-coated waveguides. The paper also gives a criterion (32) for the allowed frequency band of each multipole harmonic.

Significance. If the derivation is correct, the algorithm is a useful and fairly explicit tool for designing multilayer cylindrical waveguides for helical-undulator radiation, including resonant-frequency equations that are independent of the source particular solution. The manuscript is clearly structured, and the matrix elements for the multilayer transfer are given in closed form, which is a strength. However, the central field-amplitude result depends entirely on the imported free-space solution of [7], which is not derived or independently validated in this manuscript. The significance is therefore conditional: the method is promising and the resonance condition is robust, but the computed amplitudes cannot yet be considered verified.

major comments (4)

[§4, Eqs. (16)–(20)] The particular solution is taken from reference [7] without derivation or independent validation. The source vector S in Eq. (22) is assembled from these free-space amplitudes, so any error in the coefficients of Eq. (20) propagates directly into the computed waveguide amplitudes in Eq. (23). The resonance equation (15) is, as the authors note, independent of S and would remain valid, which makes the potential error easy to miss. Please validate the imported solution by direct substitution into Maxwell's equations, by checking a known limiting case (for example v_phi going to zero with fixed v_z, which should reduce to the linear-motion solution), or by comparison with numerical integration.
[§3, Eq. (8)] The elements beta_13, beta_23, beta_24, and beta_43 are written with the argument a1 of the Hankel functions, but the matrix W_H is defined on the boundary of the outer infinite layer at r = a_{N+1}. This is inconsistent with the text immediately above Eq. (7), and if correct as printed it would make the homogeneous matrix D wrong for multilayer walls, thereby affecting the resonance frequencies obtained from Eq. (15). Please correct these arguments to a_{N+1} or explicitly justify why a1 appears.
[§7, Eq. (32)] The allowed-band inequality is incorrect as printed. From nu_{m,0} = 0 with nu_{m,0} = sqrt(omega^2/c^2 - (omega - m omega_0)^2/v_z^2), the endpoints are omega = m omega_0 / (1 ± v_z/c). The printed denominator 1 ± v_z^2/c^2 gives a band that is too narrow by the factor (1 ± v_z^2/c^2)/(1 ± v_z/c). This error shifts the band limits that are marked with crosses in Figures 2–6 and affects the discussion of the allowed frequency region throughout Section 7.
[§6 and §7] No quantitative validation against any previously known limit is provided. The special-case formulas in Eqs. (26)–(31) could be checked by setting the wall conductivity to infinity and comparing the resonant frequencies and amplitudes with the ideal-waveguide solution of [4], and the free-space limit could be checked by removing the waveguide wall. Please add at least one such test to demonstrate that the algorithm reproduces known results before claiming that the numerical examples are reliable.

minor comments (6)

[Figure 5 and Figure 6 captions] The captions state 'permittivity epsilon_2 = 10 μm', which mixes a dielectric constant with a length unit; the relative permittivity should be dimensionless (probably epsilon_2 = 10). Please correct the captions and check the corresponding text.
[§3, Eq. (7)] The matrices W_J and W_H are displayed in an unusual block format with eight rows. For readability, please present them as explicit 4x4 matrices with column and row labels identifying the tangential field components.
[§6, Eq. (28)] The notations X and X1 are introduced inside a brace expression that is difficult to parse. Please define X and X1 separately or write the expressions as full vectors rather than using the compressed brace notation.
[§3, Eq. (8)] The element beta_23 contains both a1 and a_{N+1}; after correcting the radial argument, please verify that all Bessel and Hankel function arguments in W_H use the same outer radius.
[§7] The statement that every low-frequency mode has a partner at high frequency with a close transverse eigenvalue is presented without proof or quantitative criterion. Please state whether this is a general fact following from Eq. (19) or an observation from the plotted examples.
[Abstract and §1] The term 'modal frequency distributions' is used in the abstract but the paper mainly computes resonant frequencies and amplitude spectra. Consider clarifying the terminology to distinguish the discrete mode spectrum from the continuous spectral distribution of the radiation.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the algorithm composes an explicitly constructed dispersion matrix with an imported free-space helical solution; nothing is fitted and then renamed as a prediction.

full rationale

The paper's derivation chain is not circular. Section 3 explicitly constructs the homogeneous dispersion matrix D, including full matrix elements (Eqs. 5-11), and the resonance condition D=0 is independent of the particle's source. The particular solution for the helical source is taken from the authors' prior work [7], but the paper reproduces the needed amplitudes in Eqs. (16)-(20), and that free-space solution is a parameter-free external result whose stated assumptions (an infinite helical trajectory in free space) do not include the target problem of a multilayer cylindrical waveguide. The central system (Eqs. 21-23), X = D^{-1} S, is a linear stitching of these two independently constructed pieces, not a reduction of the output to the input by definition. No parameter is fitted to a subset of the computed results and then called a prediction; material and geometric parameters are all specified a priori. The heavy self-citation reflects reliance on the authors' earlier Green's-function and transfer-matrix work, but that earlier content is either restated in the present paper or is an independent mathematical input, so it does not constitute circularity. The algebraic band-endpoint formula (32) appears inconsistent with the dispersion relation (19), but that is a correctness concern, not a circularity concern.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

The central claim depends on established electromagnetic theory, the method of partial regions, and the previously derived free-space helical solution. No free parameters are fitted; all numerical examples use physical inputs. The paper does not introduce new physical entities.

assumptions (3)

domain assumption The free-space helical radiation solution of reference [7] is exact and complete.
Imported as the particular solution in Section 4, equations (16)-(20), without re-derivation.
standard math Modal expansions in each partial region converge and can be matched at cylindrical boundaries.
Method of partial regions assumes completeness of Bessel/Hankel expansions; no convergence proof is given.
domain assumption Each wall layer is isotropic with frequency-dependent relative permittivity and permeability.
Stated in Section 2; the derivation does not cover anisotropic or nonlinear materials.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Radiation of a particle performing helical motion in a multilayer cylindrical waveguide." pith.science (2026). https://pith.science/paper/UCGXFOYC

@misc{pith2026250112802,
  author       = {Pith},
  title        = {Pith review of: Radiation of a particle performing helical motion in a multilayer cylindrical waveguide},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UCGXFOYC}},
  note         = {Machine review of arXiv:2501.12802}
}

read the original abstract

- An algorithm for calculating the radiation field of a charged point particle performing a spiral motion in an infinite cylindrical waveguide with a multilayer side wall is found. The number of layers and their filling is arbitrary. The axis of the spiral is aligned with the axis of the waveguide, so that the geometry of the problem has cylindrical symmetry. Explicit expressions for modal frequency distributions and equations for resonant frequencies for single-layer and double-layer waveguides are given. Examples of graphical constructions of modal frequency distributions of modes for single-layer (resistive), double-layer (metal-dielectric) and triple-layer (metal-dielectric with internal NEG coating) waveguides are presented.

Figures

Figures reproduced from arXiv: 2501.12802 by the authors.

**Figure 2.** Spectral distribution of amplitudes (a) and radial electric component at [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗

**Figure 3.** Spectral distribution of amplitudes (a) and radial electric component at [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗

Figures from the paper (1 more)

**Figure 5.** Figure 5: Spectral distribution of amplitudes (a) and radial electric component at [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

[6]

Radiation of a particle moving along a helical trajectory in a resistive-wall cylindrical waveguide

Ivanyan, M., B. Grigoryan, , A. Sargsyan, A. Grigoryan, K. Floettmann, and F. Lemery, “Radiation of a particle moving along a helical trajectory in a resistive-wall cylindrical waveguide”, in Proc. 13th Int. Particle Accelerator Conf. (IPAC’22), Bangkok, Thailand, JACOW Publishing, Geneva, Switzerland, June 12-17, 2150-2153, 2022

work page 2022
[7]

A new approach to solving the radiation field problem of an extended helical undulator

Ivanyan M. I., B. Grigoryan, A. Grigoryan, L. Aslyan, V. Avagyan, H. Babujyan, S. Arutunian, K. Floettmann, F. Lemery, “A new approach to solving the radiation field problem of an extended helical undulator”, arxiv.org/abs/2410.05730

work page Pith review arXiv
[4]

Synchrotron radiation in a circular metallic waveguide

Dokuchaev, V. P., “Synchrotron radiation in a circular metallic waveguide”, Izv. VUZov Radiofiz, Vol. 44, No 7, 587- 591, 2001

work page 2001
[1]

Radiation of an electron in an infinitely long waveguide

Motz, H. and M. Nakamura, “Radiation of an electron in an infinitely long waveguide”, Annals of Physics, Vol. 7, Issue 1, 84-131, 1959

work page 1959
[2]

Radiation from a Charge Moving Along a Spiral Line Inside a Cylindrical Waveguide

Karapetian, G. G., “Radiation from a Charge Moving Along a Spiral Line Inside a Cylindrical Waveguide”, Izv.Akad.Nauk Arm.SSR Fiz. No 12, 186-190, 1977

work page 1977
[3]

Synchrotron radiation of wiggled electron beam in rectangular waveguide

Haus, H. A. and N. Islam, “Synchrotron radiation of wiggled electron beam in rectangular waveguide”, J. Appl. Phys. Vol. 54, 4784–4793, 1983

work page 1983
[5]

Mode enhancement of THz helical undulator radiation in waveguide

Ivanyan, M. I., A. H. Grigoryan, B. K. Sargsyan, T. L. Vardanyan, A. V. Tsakanian, and V. M. Tsakanov, “Mode enhancement of THz helical undulator radiation in waveguide”, Laser Phys. Vol. 30, 115002, 2020

work page 2020
[8]

Wake fields and ohmic losses in round vacuum chambers

Piwinski A., “Wake fields and ohmic losses in round vacuum chambers”, DESY, Hamburg, Germany, Rep. DESY HERA 92-11, May 1992

work page 1992

Show all 16 references

[9]

Multi-layer tube impedance and external radiation,

Ivanyan, M., E. Laziev, V. Tsakanov, A. Vardanyan, S. Heifets. and A. Tsakanian, “Multi-layer tube impedance and external radiation," Physical Review Special Topics-Accelerators and Beams, Vol. 11, No. 8, 084001, 2008

2008
[10]

Dispersion relations for a cylindrical waveguide with multilayer walls

Ivanyan, M. I., B. Grigoryan, L Aslyan, A. Grigoryan, K. Floettmann, F. Lemery, “Dispersion relations for a cylindrical waveguide with multilayer walls”, in Proc. 14th International Particle Accelerator Conference (IPAC23), Venice, Italia, JACOW Publishing, Geneva, Switzerland...

2023
[11]

Impedances and wakes in high energy particle accelerators

Zotter, B. W. and S. A. Kheifetz, “Impedances and wakes in high energy particle accelerators”, Singapore, World Scien-tific, 1997

1997
[12]

Wake fields in a three-layer cylindrical waveguide

Ivanyan, M., B. Grigoryan, L. Aslyan, A. Grigoryan, A. Vardanyan, V. Avagyan, “Wake fields in a three-layer cylindrical waveguide”, 2023 PhotonIcs & Electromagnetics Research Symposium (PIERS), Prague, Czech Republic, July3- 6, 71-78, 2023

2023
[13]

Electromagnetic theory

Stratton, J., A., “Electromagnetic theory”, McGraw-Hill book Company, Inc., Cambridge, 615 p., 1941

1941
[14]

The short-range resistive wall wakefields,

Bane, K. L. F. and M. Sands, “The short-range resistive wall wakefields," AIP Conference Proceedings, Vol. 367, No. 1, 131{149, American Institute of Physics, 1996

1996
[15]

Wakefelds in conducting waveguides with a lossy dielectric channel,

Ivanyan, M., L. Aslyan, K. Floettmann, F. Lemery, and V. Tsakanov, “Wakefelds in conducting waveguides with a lossy dielectric channel," Physical Review Accelerators and Beams, Vol. 23, No. 4, 041301, 2020

2020
[16]

Characterization of NEG Coatings for SLS 2.0,

Dehler, M., S. Alberti, A. Citterio, M. Hahn, J.-P. Hogge, X. Liu, and H. Marques, “Characterization of NEG Coatings for SLS 2.0," 10th Int. Particle Accelerator Conf. (IPAC'19),1662{1665, JACOW Publishing, Geneva, Switzerland, Melbourne, Australia, May 19-24, 2019

2019

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.

[1] [6]

Radiation of a particle moving along a helical trajectory in a resistive-wall cylindrical waveguide

Ivanyan, M., B. Grigoryan, , A. Sargsyan, A. Grigoryan, K. Floettmann, and F. Lemery, “Radiation of a particle moving along a helical trajectory in a resistive-wall cylindrical waveguide”, in Proc. 13th Int. Particle Accelerator Conf. (IPAC’22), Bangkok, Thailand, JACOW Publishing, Geneva, Switzerland, June 12-17, 2150-2153, 2022

work page 2022

[2] [7]

A new approach to solving the radiation field problem of an extended helical undulator

Ivanyan M. I., B. Grigoryan, A. Grigoryan, L. Aslyan, V. Avagyan, H. Babujyan, S. Arutunian, K. Floettmann, F. Lemery, “A new approach to solving the radiation field problem of an extended helical undulator”, arxiv.org/abs/2410.05730

work page Pith review arXiv

[3] [4]

Synchrotron radiation in a circular metallic waveguide

Dokuchaev, V. P., “Synchrotron radiation in a circular metallic waveguide”, Izv. VUZov Radiofiz, Vol. 44, No 7, 587- 591, 2001

work page 2001

[4] [1]

Radiation of an electron in an infinitely long waveguide

Motz, H. and M. Nakamura, “Radiation of an electron in an infinitely long waveguide”, Annals of Physics, Vol. 7, Issue 1, 84-131, 1959

work page 1959

[5] [2]

Radiation from a Charge Moving Along a Spiral Line Inside a Cylindrical Waveguide

Karapetian, G. G., “Radiation from a Charge Moving Along a Spiral Line Inside a Cylindrical Waveguide”, Izv.Akad.Nauk Arm.SSR Fiz. No 12, 186-190, 1977

work page 1977

[6] [3]

Synchrotron radiation of wiggled electron beam in rectangular waveguide

Haus, H. A. and N. Islam, “Synchrotron radiation of wiggled electron beam in rectangular waveguide”, J. Appl. Phys. Vol. 54, 4784–4793, 1983

work page 1983

[7] [5]

Mode enhancement of THz helical undulator radiation in waveguide

Ivanyan, M. I., A. H. Grigoryan, B. K. Sargsyan, T. L. Vardanyan, A. V. Tsakanian, and V. M. Tsakanov, “Mode enhancement of THz helical undulator radiation in waveguide”, Laser Phys. Vol. 30, 115002, 2020

work page 2020

[8] [8]

Wake fields and ohmic losses in round vacuum chambers

Piwinski A., “Wake fields and ohmic losses in round vacuum chambers”, DESY, Hamburg, Germany, Rep. DESY HERA 92-11, May 1992

work page 1992