pith. sign in

arxiv: 2605.01249 · v1 · submitted 2026-05-02 · ⚛️ physics.optics

A New Perspective on Matrix Representation of Paraxial Geometric Optics using Two Kinds of Three-Matrix Decompositions of the 2times 2 Special-Linear-Group Matrices

Satoshi Itoh This is my paper

Pith reviewed 2026-05-09 18:44 UTC · model grok-4.3

classification ⚛️ physics.optics
keywords ABCD matrixparaxial opticsmatrix decompositionoptical designrefraction surfacesspecial linear groupgeometric opticsray transfer
0
0 comments X

The pith

ABCD matrices in paraxial optics admit two distinct three-matrix decompositions connected by a transformation that raises or lowers the number of refraction surfaces while leaving the overall paraxial ray transfer unchanged.

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

The paper establishes that any 2x2 ABCD matrix, which encodes the paraxial behavior of a rotationally symmetric optical system and possesses exactly three real degrees of freedom, can be factored into a product of three elementary matrices in each of two different patterns. One pattern corresponds to a particular ordering of thin-lens powers and free-space propagations; the second pattern uses a different ordering. The authors derive an explicit algebraic map that converts any factorization of the first kind into one of the second kind, and vice versa. Because the map preserves the product exactly, the composite paraxial imaging properties remain identical even though the number of individual refraction surfaces changes. This supplies a designer with a family of physically distinct layouts that all satisfy the same paraxial specification, thereby opening a route to choose the layout that later minimizes aberrations when non-paraxial effects are taken into account.

Core claim

We propose two kinds of three-matrix decomposition of ABCD matrices by focusing on the fact that the ABCD matrices have three real-number degrees of freedom. In addition, we formulate a transformation between the two kinds of decomposition for a single matrix, which can increase or decrease the number of refraction surfaces in the optical configuration while keeping the paraxial specifications fixed.

What carries the argument

A pair of three-matrix factorizations of an SL(2,R) matrix together with the algebraic transformation that interconverts them, each factorization representing a distinct sequence of thin-lens powers and propagation distances whose product recovers the original ABCD matrix.

If this is right

  • Any paraxial specification can be realized by optical trains containing different numbers of refraction surfaces.
  • The transformation supplies an exact algebraic way to move between these trains without disturbing the overall ray-transfer matrix.
  • Designers can therefore enumerate candidate layouts that share identical paraxial focal length, magnification, and principal-plane locations.
  • The same paraxial matrix can be used as the starting point for subsequent non-paraxial optimization of each candidate layout.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The method could be inserted into automated design loops that iterate over surface count as a discrete parameter before continuous optimization of curvatures and spacings begins.
  • One could test the practical utility by starting from the ABCD matrix of a simple doublet, generating the transformed four-surface version, and comparing the Seidel aberration coefficients of both realizations.
  • The same algebraic structure might extend to systems with thick lenses or GRIN media once the appropriate elementary matrices are identified.

Load-bearing premise

That the extra freedom to change surface count through these decompositions will translate into measurable improvements when real aberrations are minimized in actual multi-element designs.

What would settle it

Take any concrete ABCD matrix that arises from a known thin-lens system, apply the claimed transformation to produce a new decomposition with a different surface count, then recompute the composite ABCD matrix from the new sequence and verify whether it equals the original matrix to machine precision.

Figures

Figures reproduced from arXiv: 2605.01249 by Satoshi Itoh.

Figure 1
Figure 1. Figure 1: Schematic of the + decomposition 𝑀 = 𝑇 (𝑎)𝑅(𝛽)𝑇 (𝑐). This diagram conceptually shows an optical configuration. In this diagram, the ray goes from right to left. The horizontal lines express the light transition. The vertical line indicates refraction. The symbols 𝛽, 𝑎, and 𝑐 denote the parameters of the operator matrices. 3.2. H Decomposition When 𝐶 ≠ 0 or when the ABCD matix 𝑀 is a refraction matrix, we c… view at source ↗
Figure 2
Figure 2. Figure 2: Schematic of the H decomposition 𝑀 = 𝑅(𝛼)𝑇 (𝑏)𝑅(𝛾). This diagram conceptually shows an optical configuration. In this diagram, the ray goes from right to left. The horizontal lines express the light transition. The vertical line indicates refraction. The symbols 𝑏, 𝛼, and 𝛾 denote the parameters of the operator matrices. angles. Thus, the matrices of [Eq. (11)] expresses a special type of optical conjugate… view at source ↗
Figure 3
Figure 3. Figure 3: Schematic for how the + ↔ H transformation increases or decreases the number of the refraction surfaces in an optical configuration. These diagrams conceptually show optical configurations. In these diagrams, the ray goes from right to left. The horizontal lines express the light transition. The vertical line indicates refraction. The + ↔ H transformation interchanges the + part within the dashed closed cu… view at source ↗
read the original abstract

We require decomposition methods for the ABCD-matrix formulation in rotationally symmetric paraxial geometric optics when designing a multi-component optical system from a given single paraxial specification (represented by an ABCD matrix) to optimize non-paraxial specifications (e.g., optical aberrations). In this study, we propose two kinds of three-matrix decomposition of ABCD matrices by focusing on the fact that the ABCD matrices have three real-number degrees of freedom. In addition, we formulate a transformation between the two kinds of decomposition for a single matrix, which can increase or decrease the number of refraction surfaces in the optical configuration while keeping the paraxial specifications fixed. This nature is useful for the optical design of multi-component systems with optimized non-paraxial characteristics.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript proposes two distinct three-matrix decompositions of 2×2 SL(2,R) ABCD matrices in paraxial geometric optics, each parameterized by the three real degrees of freedom of the group. It further derives an explicit algebraic transformation between the two decompositions for any given matrix, which is claimed to permit increasing or decreasing the number of refraction surfaces in an equivalent optical configuration while exactly preserving the overall paraxial ray-transfer matrix. The authors argue that this property is useful for designing multi-component systems with improved non-paraxial performance such as reduced optical aberrations.

Significance. If the decompositions and the transformation are algebraically correct, the work supplies a systematic, parameter-counting method for inserting or removing thin-lens or surface elements without disturbing the target ABCD matrix. This could, in principle, give designers additional degrees of freedom to minimize Seidel or higher-order aberrations while holding paraxial specifications fixed. The approach is grounded in the standard matrix formalism of paraxial optics and does not introduce new physical assumptions.

major comments (2)
  1. [Abstract] Abstract and concluding section: the central utility claim—that the transformation 'can increase or decrease the number of refraction surfaces … while keeping the paraxial specifications fixed' and thereby optimizes 'non-paraxial characteristics'—is asserted without any concrete numerical example, ray-trace verification, or aberration calculation demonstrating that the added surfaces actually reduce aberrations for a fixed ABCD matrix.
  2. [Decomposition sections] The manuscript supplies the algebraic parametrizations of the two decompositions and the explicit map between them, yet contains no verification (analytic identity or numerical check) that the product of the three matrices in each decomposition recovers the original ABCD matrix for arbitrary parameter values.
minor comments (1)
  1. Notation for the two decomposition families is introduced without a compact summary table comparing the parameter ranges or the explicit matrix elements of each kind.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and constructive feedback on our manuscript. The comments identify key areas where additional verification and demonstration would strengthen the presentation. We respond to each major comment below and outline the revisions we will implement.

read point-by-point responses

  1. Referee: [Abstract] Abstract and concluding section: the central utility claim—that the transformation 'can increase or decrease the number of refraction surfaces … while keeping the paraxial specifications fixed' and thereby optimizes 'non-paraxial characteristics'—is asserted without any concrete numerical example, ray-trace verification, or aberration calculation demonstrating that the added surfaces actually reduce aberrations for a fixed ABCD matrix.

    Authors: We agree that the utility claim benefits from concrete illustration. The manuscript emphasizes the algebraic construction that enables varying the number of surfaces while exactly preserving the ABCD matrix. In the revised version we will add a specific numerical example: we select a representative ABCD matrix, apply the transformation to obtain equivalent decompositions with increased and decreased numbers of refraction surfaces, confirm the paraxial matrix is unchanged, and compute the resulting Seidel aberration coefficients to demonstrate how the added degrees of freedom can be used to reduce aberrations. revision: yes

  2. Referee: [Decomposition sections] The manuscript supplies the algebraic parametrizations of the two decompositions and the explicit map between them, yet contains no verification (analytic identity or numerical check) that the product of the three matrices in each decomposition recovers the original ABCD matrix for arbitrary parameter values.

    Authors: This observation is correct. Although the parametrizations and the transformation map are derived from the three real degrees of freedom of SL(2,R), an explicit verification step is missing. In the revision we will insert both an analytic identity proving that the product of the three matrices equals the original ABCD matrix for general parameter values and a numerical check performed on several arbitrary matrices to confirm the decompositions. revision: yes

Circularity Check

0 steps flagged

Direct algebraic decompositions with no self-referential reductions

full rationale

The paper constructs two explicit three-matrix factorizations of SL(2,R) ABCD matrices by parameterizing the known three real degrees of freedom, then derives an algebraic map between the factorizations that preserves the composite matrix. These steps are purely constructive and rely only on standard matrix algebra and the determinant-1 constraint; no parameters are fitted to data, no predictions are generated from subsets of results, and no self-citations are invoked to justify the decompositions or the transformation. The claimed utility for non-paraxial design is stated as a motivation rather than a derived result, leaving the core derivations self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard fact that ABCD matrices form the group SL(2,R) with determinant 1 and therefore possess three real degrees of freedom; no additional free parameters or invented entities are introduced in the abstract.

axioms (1)
  • standard math ABCD matrices belong to SL(2,R) and therefore have exactly three real degrees of freedom.
    Invoked in the abstract to justify the three-matrix decomposition.

pith-pipeline@v0.9.0 · 5434 in / 1175 out tokens · 28655 ms · 2026-05-09T18:44:43.338942+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

8 extracted references · 8 canonical work pages

  1. [1]

    2017 , publisher=

    Optics , author=. 2017 , publisher=

    work page 2017

  2. [2]

    American Journal of Physics , volume=

    Radiance , author=. American Journal of Physics , volume=. 1963 , publisher=

    work page 1963

  3. [3]

    2003 , publisher=

    Lie Groups, Lie Algebras, and Representations: An Elementary Introduction , author=. 2003 , publisher=

    work page 2003

  4. [4]

    Photoinduced active terahertz metamaterials with nanostructured vanadium dioxide film deposited by sol-gel method , volume =

    Yaxin Zhang and Shen Qiao and Linlin Sun and Qi Wu Shi and Wanxia Huang and Ling Li and Ziqiang Yang , journal =. Photoinduced active terahertz metamaterials with nanostructured vanadium dioxide film deposited by sol-gel method , volume =. 2014 , url =

    work page 2014

  5. [5]

    Climate Change 2007: The Physical Science Basis

    P Forster and V Ramaswamy and P Artaxo and T Bernsten and R Betts and D Fahey and J Haywood and J Lean and D Lowe and G Myhre and J Nganga and R Prinn and G Raga and M Schulz and R V Dorland , title=. Climate Change 2007: The Physical Science Basis. Contribution of Working Group 1 to the Fourth Assesment Report of Intergovernmental Panel on Climate Change...

    work page 2007

  6. [6]

    Multiphoton plasmon-resonance microscopy , volume =

    Dvir Yelin and Dan Oron and Stephan Thiberge and Elisha Moses and Yaron Silberberg , journal =. Multiphoton plasmon-resonance microscopy , volume =. 2003 , url =. doi:10.1364/OE.11.001385 , abstract =

    work page doi:10.1364/oe.11.001385 2003

  7. [7]

    Cluster formation in ferrofluids induced by holographic optical tweezers , volume =

    Jan Masajada and Marcin Bacia and S. Cluster formation in ferrofluids induced by holographic optical tweezers , volume =. Opt. Lett. , keywords =. 2013 , url =. doi:10.1364/OL.38.003910 , abstract =

    work page doi:10.1364/ol.38.003910 2013

  8. [8]

    and Shiri, Ron and Acton, Scott D

    Dean, Bruce H and Aronstein, David L and Smith, Scott J. and Shiri, Ron and Acton, Scott D. , title =. Space Telescopes and Instrumentation I: Optical, Infrared, and Millimeter , year =. doi:10.1117/12.673569 , url =

    work page doi:10.1117/12.673569