Wien-Filter Hamiltonian and Transfer Matrix

Volker Ziemann

REVIEW 2 major objections 7 minor 10 references

Wien-Filter Hamiltonian and Transfer Matrix

T0 review · 2 major / 7 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read This paper derives, from first principles, the first-order six-dimensional transfer matrix of an ideal Wien filter with vertical magnetic field, and shows that the whole matrix is parametrized by the spin-rotation angle $\phi = L/R$.

desk verdict A clean 6D first-order Wien-filter transfer matrix, marred by a sign typo in the central Hamiltonian equation that must be fixed before the derivation works. read the letter →

arxiv 2608.06052 v1 pith:AEJSBARM submitted 2026-08-06 physics.acc-ph physics.app-phphysics.ins-det

Volker Ziemann This is my paper · ORCID

classification physics.acc-phphysics.app-phphysics.ins-det

keywords WienfiltertransfermatrixHamiltonianbeamopticsspinprecessiondispersionlongitudinalcouplingsymplecticmap

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 derives, from first principles, the first-order six-dimensional transfer matrix of an ideal Wien filter whose magnetic field points vertically. The central claim is that the full 6×6 matrix is fixed by a single inverse length $R = \gamma_0(B\rho)/B_0$, equivalently by the spin-rotation angle $\phi = L/R$, so the filter's focusing, dispersion, and longitudinal coupling are all determined by the same parameter that governs spin precession. A sympathetic reader would care because Wien filters are standard in spin-polarized beams, yet their effect on ordinary beam optics is usually treated only in reduced dimensions; this gives a complete matrix that existing beam-optics codes can use. The paper also points out that reversing the field polarity reverses the sign of the dispersion and timing-coupling elements while leaving the focusing block unchanged, a systematic effect that high-precision experiments may need to account for.

What carries the argument

The load-bearing object is the simplified Hamiltonian $$H = \frac{\$delta^{2}$}{2\$gamma_0^{2}$} + \frac{1}{2}x'^2 + \frac{1}{2}y'^2 + \frac{$x^{2}$}{$2R^{2}$} + \frac{x\delta}{\gamma_0 R},$$ with $R = \gamma_0(B\rho)/B_0$. It turns the ideal Wien filter into a drift plus a horizontal focusing term and an energy-dispersion term, and it generates the transfer matrix Eq. (14) whose entries are $\sin$ and $\cos$ of $\phi = L/R$. The identity that carries the argument is that $\phi$, defined by the geometry and field strengths, is exactly the spin-rotation angle given by the spin-precession equation, which is what makes the transfer matrix already parametrized by spin rotation.

What would settle it

Measure the full 6×6 response of a real Wien filter by scanning entrance values of $x$, $x'$, $y$, $y'$, $z$, and $\delta$, fit the exit coordinates, and compare every entry with Eq. (14). The clearest test is the polarity flip: Eq. (14) predicts that $R_{16}$, $R_{26}$, $R_{51}$, and $R_{52}$ change sign when the field polarity is reversed while the 4×4 focusing block stays fixed; if those elements do not flip, or if the longitudinal element $R_{56} = R\sin(\phi)/\gamma_0^2$ disagrees, the matrix is not the full story.

Watch

Extended reading notes

Core claim

Starting from the general accelerator Hamiltonian, the paper imposes the Wien condition $B_0 = -E/(\beta_0 c)$ so that the reference trajectory is straight, expands to second order in the dynamical variables, and solves Hamilton's equations. The result, Eq. (14), maps entrance coordinates to exit coordinates: $x$ and $x'$ rotate through the angle $\phi = L/R$ with focusing, $y$ drifts, $\delta$ is constant, and $z$ receives coupling terms from $x$, $x'$, and $\delta$. The key claim is that $R = \gamma_0(B\rho)/B_0$ is not the ordinary magnetic bending radius but an effective scale, and that the spin-precession angle of an electron in the filter is exactly $\phi = L/R$; therefore the transfer matrix is already parametrized by the spin-rotation angle. The author notes that the transverse elements agree with an earlier five-dimensional derivation, and that the result extends earlier four- and five-dimensional work to full six-dimensional phase space.

Load-bearing premise

The whole matrix rests on treating the Wien filter as ideal: perfectly uniform electric and magnetic fields that cancel exactly on the reference trajectory, abrupt hard edges with no fringe fields, no electrode misalignment, and small-angle paraxial motion.

Editorial extensions

If this is right

Any first-order beam-optics code can model an ideal Wien filter in full 6D phase space by inserting Eq. (14) as a standard transfer matrix.
For a given spin-rotation angle the optics are fixed: a 90-degree rotator creates a dispersion of order $R/\gamma_0$ and couples a 1 mm horizontal offset into about 0.72 mm of longitudinal path change for the example parameters.
Reversing the Wien-filter polarity to flip the spin leaves the 4×4 focusing block unchanged but changes the sign of the dispersion elements ($R_{16}$, $R_{26}$) and of the longitudinal coupling elements ($R_{51}$, $R_{52}$), so the beam's energy and timing response flips systematically.
A Wien filter with horizontal magnetic field can be handled by sandwiching Eq. (14) between two 90-degree coordinate rotations.
Because the matrix is symplectic by construction, it preserves phase-space volume and can be chained with other first-order maps without breaking the symplectic structure.

Reading between the lines

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

If Eq. (14) is right, a Wien filter's optical strength cannot be chosen independently of its spin-rotation angle: specifying $\phi$ fixes $R$, and therefore fixes focusing and dispersion, a design constraint for injectors that need both.
The paper's own validation compares only one dominant focal length (measured 0.78 m versus estimated 0.88 m); a direct measurement of the off-diagonal elements, especially $R_{16}$ and $R_{51}$ whose signs flip with polarity, would be a much sharper test of the full matrix.
The hard-edged idealization ignores fringe fields, but Eq. (14) gives the limiting map that any fringe-field model must approach; comparing it with particle tracking through realistic fields would show how much of the 0.1 m focal-length discrepancy comes from fringing rather than from the matrix itself.
The same parametrization-by-spin-angle structure may extend to other crossed-field devices, such as electrostatic separators or velocity selectors, tying their optics to a physical precession angle rather than to separate field values.

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

Volker Ziemann's note gives the full 6D first-order transfer matrix for an ideal Wien filter, expressed neatly in terms of the spin-rotation angle phi = L/R. That is a legitimate and useful extension of the earlier 4D and 5D work, and the paper is well-written and compact. The result itself looks right, and the matrix is already parametrized by the spin angle, which is a nice observation for CEBAF/MOLLER-type applications. The derivation from the Hamiltonian is transparent and the equations of motion are straightforward.

But the manuscript as submitted has a load-bearing typo. In Eq (3), the square root contains (beta0*delta + 1/beta0 + eEx/(c p0))^2. Expanding that gives a linear term eEx/(beta0 c p0) with a plus sign inside, so the Hamiltonian's linear x-term becomes eB0 x/p0 - eEx/(beta0 c p0). That does not cancel under the stated Wien condition B0 = -E/(beta0 c); it doubles. The next equation, Eq (4), uses the opposite sign, eB0 x/p0 + eEx/(beta0 c p0), which is what you get if the square root had eEx/(c p0) with a minus sign. So the printed Equations (3)-(6) are internally inconsistent, and the central claim does not follow from the equations as written. This is almost certainly a sign typo, but it sits at the root of the derivation and must be corrected.

There are also minor typos in Eq (12)/(13): the x2' line of Eq (13) has cos(phi) x2 where it should be cos(phi) x1'. Eq (12) is a bit redundant but simplifies correctly.

The experimental validation is the weakest part of the paper. They compare a predicted focal length of 0.88 m to a measured f_y = 0.78 m, but the ideal Wien filter focuses in the horizontal plane, not vertical. There is presumably skew and imperfection in the real filter, so the agreement is only circumstantial. A cleaner statement of which plane the measured focal length belongs to, and how much the imperfections matter, would improve the paper. This is not a failure of the theory, just an over-soft validation.

I think the core result is correct and useful, and the sign typo is fixable. This deserves a serious referee, but the referee should require the sign fix and a clarification of the focal-length comparison. I'd give it a conditional accept.

Referee Report

2 major / 7 minor

Summary. The paper derives a six-dimensional first-order transfer matrix for an ideal Wien filter with a vertical magnetic field, starting from the general Hamiltonian of Ref. [4]. After imposing the Wien condition B0 = -E/(beta0 c), the Hamiltonian reduces to a quadratic form describing drift, horizontal focusing, and x-delta dispersion. The transfer matrix is expressed in terms of the spin-rotation angle phi = L/R with R = gamma0 (Brho)/B0, and the authors compare the predicted focusing strength with a measured focal length. The paper also discusses the effect of reversing the field polarity on dispersion and longitudinal coupling.

Significance. If correct, the result provides a compact, parameter-free 6D transfer matrix for ideal Wien filters that can be inserted into beam-optics codes, extending earlier 4D and 5D derivations and connecting the transverse optics to the spin-rotation angle. The derivation is non-circular, the final matrix is symplectic, and the spin-angle parametrization is an elegant and practical addition. The main limitations are the ideal-field assumptions (hard edges, exact Wien condition) and the relatively coarse experimental check, which tests only the focusing block.

major comments (2)

[Section 2, Eq. (3)] The sign of eEx/(c p0) inside the square root is inconsistent with the rest of the derivation. Expanding Eq. (3) as printed gives a linear term eB0x/p0 - eEx/(beta0 c p0) in the Hamiltonian, so the stated Wien condition B0 = -E/(beta0 c) doubles this term instead of canceling it. The expansion in Eq. (4) and the cancellation leading to Eqs. (5)-(6) require the opposite sign, namely (beta0 delta + 1/beta0 - eEx/(c p0))^2 inside the root. Since this cancellation is the linchpin for the transfer matrix Eq. (14), the sign must be corrected for the derivation to be valid.
[Section 4, validation paragraph] The comparison of the predicted focal length f ≈ 0.88 m to the measured f_y = 0.78 m does not test the model as presented. For a vertical magnetic field, Eq. (5) contains x^2/(2R^2) but no y^2 term, so the model predicts horizontal focusing only; the quoted measured focal length is labeled f_y. The text should specify the plane of the measurement and how the reported skew component affects the comparison; without this, the apparent agreement does not provide independent support for the matrix elements.

minor comments (7)

[Section 3, Eq. (12)] In the first equality for z(s), the x'_1 contribution should be R x'_1 / gamma0 (1 - cos(s/R)), not R x'_1 / gamma0 cos(s/R); the second equality is correct.
[Section 3, Eq. (13)] In the equation for x'_2, the coefficient of x1 should be cos(phi) x'_1, not cos(phi) x2; as printed, this line is not the mapping given in Eq. (14).
[Section 3, Eq. (11)] The equation 'dδ'/ds' should read 'dδ/ds' because δ is the dynamical momentum and is not primed.
[Section 2, Eq. (4)] The notation 'δ2/2 (1−β0²)' is confusing; the subsequent simplification to δ²/(2γ0²) is correct, but the intermediate form should be written consistently.
[Section 4, validation] If the measured focal length f_y is actually the horizontal focal length, the subscript y is misleading; please use a notation that makes the plane of focusing explicit.
[Section 5] In the last paragraph, 'can can be found' is a typo for 'can be found'.
[Abstract] The phrase 'from first principles' overstates the starting point, since the derivation begins from the Hamiltonian of Ref. [4]; consider rewording to 'from a standard Hamiltonian'.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the transfer matrix follows from a standard Hamiltonian by parameter-free expansion and integration, with the spin-angle identification and focal-length check providing independent support.

full rationale

The derivation is self-contained after adopting the standard relativistic Hamiltonian in Eq. (1), which is attributed to the author's textbook [4] but is an externally established expression for a charged particle in electromagnetic potentials; using it is not circular. The Wien condition, the expansion to second order, the equations of motion, and the assembled matrix Eq. (14) follow by algebra and integration without fitting any parameter to the target result. The parameter R is defined from the fields, and Eq. (15) independently derives the spin-rotation angle from the Thomas-BMT equation, showing that it equals L/R; this is a connecting identity, not an input used to construct the matrix. The comparison with the measured focal length (0.78 m versus the predicted 0.88 m) is a genuine prediction from the derived matrix, not a fitted constant. The only self-citations, [4] for the formalism and [7] for the measured value, are standard material and an external measurement respectively, so they do not make the argument circular. The paper's own caveats about ideal fields and the absence of fringe fields are accuracy limitations rather than circularity. A separate algebraic sign inconsistency between Eqs. (3) and (4) is a correctness concern outside the scope of circularity analysis.

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

No free parameters are fitted in the derivation; R is defined from physical inputs and checked against the spin-rotation angle. No new physical entities are introduced. The main assumptions are the standard paraxial and ideal-field modeling choices, the prior Hamiltonian formalism, and the Thomas-BMT spin-precession result.

assumptions (4)

domain assumption Paraxial approximation: x' and y' are small, so the square root can be expanded to second order.
Invoked in Section 2 with the phrase 'we tacitly assume the paraxial approximation with small angles x' and y'. This limits the validity to small-angle beams.
domain assumption The general Hamiltonian in Eq (1) from reference [4] is the correct starting point for beam optics in the s-dependent phase space.
Section 2 states 'We base our analysis on the formalism outlined in Section 2 in [4]'. This Hamiltonian is a textbook result, not re-derived in the paper.
domain assumption The Wien filter is ideal: hard-edged fields, no misalignment, no fringe fields, and 1/rho = 0 with the exact Wien condition B0 = -E/(beta0*c).
Section 1 defines 'ideal' and Section 2 assumes the trajectory is straight. The paper explicitly acknowledges at the end that fringe fields are not included.
domain assumption The Thomas-BMT spin precession formula from reference [6] gives the spin-rotation angle as phi = B0*L/(gamma0*(Brho)).
Section 4 uses this to identify phi = L/R. This is a prior published result, not derived in the paper.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Wien-Filter Hamiltonian and Transfer Matrix." pith.science (2026). https://pith.science/paper/AEJSBARM

@misc{pith2026260806052,
  author       = {Pith},
  title        = {Pith review of: Wien-Filter Hamiltonian and Transfer Matrix},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/AEJSBARM}},
  note         = {Machine review of arXiv:2608.06052}
}

read the original abstract

We derive the Hamiltonian and the transfer matrix of a Wien-filter with vertical magnetic field from first principles.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

10 extracted references · 8 canonical work pages

[4]

Ziemann,Hands-On Accelerator Physics Using MATLAB, 2nd edi- tion,CRC Press, Boca Raton, 2025

V. Ziemann,Hands-On Accelerator Physics Using MATLAB, 2nd edi- tion,CRC Press, Boca Raton, 2025

work page 2025
[1]

Salomaa, H

M. Salomaa, H. Enge,Velocity selector for heavy-ion separation,Nuclear Instruments and Methods 145 (1977) 279. 6

work page 1977
[2]

Hurd,Derivation of the first-order transformation matrix for a simple Wien filter and comparison to results of numerical integration,Nuclear Instruments and Methods A258 (1987) 542

J. Hurd,Derivation of the first-order transformation matrix for a simple Wien filter and comparison to results of numerical integration,Nuclear Instruments and Methods A258 (1987) 542

work page 1987
[3]

Brown,A First- and Second-Order Matrix Theory for the Design of Beam Transport Systems and Charged Particle Spectrometers,SLAC-75 (1967)

K. Brown,A First- and Second-Order Matrix Theory for the Design of Beam Transport Systems and Charged Particle Spectrometers,SLAC-75 (1967)

work page 1967
[5]

Bargmann, L

V. Bargmann, L. Michel, and V.L. Telegdi,Precession of the Polarization of Particles Moving in a Homogeneous Electromagnetic Field, Phys. Rev. Lett. 2 (1959) 435

work page 1959
[6]

Tioukine, K

V. Tioukine, K. Aulenbacher,Operation of the MAMI accelerator with a Wien filter based spin rotation system,Nuclear Instruments and Methods A568 (2006) 537

2006
[7]

Ziemann, M

V. Ziemann, M. Bruker,Response-matrix analysis of the injector and measuring the Wien-filter transfer matrices,JLAB-TN-25-054, August 2025

work page 2025
[8]

Moller collaboration,The MOLLER Experiment: An Ultra-Precise Mea- surement of the Weak Mixing Angle Using Møller Scattering, arXiv 1411.4088 (2014)

arXiv 2014

Show all 10 references

[9]

pulsar.nl/gpt/

General Particle Tracer, Pulsar Physics web page athttps://www. pulsar.nl/gpt/
[10]

Ioanoviciu, K

D. Ioanoviciu, K. Tsuno,Third order transfer matrices for real Wien filters with homogeneous main fields, Nuclear Instruments and Methods A614 (2010) 47. 7

2010

Pith tools

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

[1] [4]

Ziemann,Hands-On Accelerator Physics Using MATLAB, 2nd edi- tion,CRC Press, Boca Raton, 2025

V. Ziemann,Hands-On Accelerator Physics Using MATLAB, 2nd edi- tion,CRC Press, Boca Raton, 2025

work page 2025

[2] [1]

Salomaa, H

M. Salomaa, H. Enge,Velocity selector for heavy-ion separation,Nuclear Instruments and Methods 145 (1977) 279. 6

work page 1977

[3] [2]

Hurd,Derivation of the first-order transformation matrix for a simple Wien filter and comparison to results of numerical integration,Nuclear Instruments and Methods A258 (1987) 542

J. Hurd,Derivation of the first-order transformation matrix for a simple Wien filter and comparison to results of numerical integration,Nuclear Instruments and Methods A258 (1987) 542

work page 1987

[4] [3]

Brown,A First- and Second-Order Matrix Theory for the Design of Beam Transport Systems and Charged Particle Spectrometers,SLAC-75 (1967)

K. Brown,A First- and Second-Order Matrix Theory for the Design of Beam Transport Systems and Charged Particle Spectrometers,SLAC-75 (1967)

work page 1967

[5] [5]

Bargmann, L

V. Bargmann, L. Michel, and V.L. Telegdi,Precession of the Polarization of Particles Moving in a Homogeneous Electromagnetic Field, Phys. Rev. Lett. 2 (1959) 435

work page 1959

[6] [6]

Tioukine, K

V. Tioukine, K. Aulenbacher,Operation of the MAMI accelerator with a Wien filter based spin rotation system,Nuclear Instruments and Methods A568 (2006) 537

2006

[7] [7]

Ziemann, M

V. Ziemann, M. Bruker,Response-matrix analysis of the injector and measuring the Wien-filter transfer matrices,JLAB-TN-25-054, August 2025

work page 2025

[8] [8]

Moller collaboration,The MOLLER Experiment: An Ultra-Precise Mea- surement of the Weak Mixing Angle Using Møller Scattering, arXiv 1411.4088 (2014)

arXiv 2014