Abstract
We present the formulation of a paraxial ray transfer or ABCD matrix for onion-type GRIN lenses. In GRIN lenses, each iso-indicial surface (IIS) can be considered a refracting optical surface. If each IIS is a shell or layer, the ABCD matrix of a GRIN lens is computed by multiplying a typically high number of translation and refraction matrices corresponding to the K layers inside the lens. Using a differential approximation for the layer thickness, this matrix product becomes a sum. The elements A, B, C, and D of the approximated GRIN ray transfer matrix can be calculated by integrating the elements of a single-layer matrix. This ABCD matrix differs from a homogeneous lens matrix in only one integration term in element C, corresponding to the GRIN contribution to the lens power. Thus the total GRIN lens power is the sum of the homogeneous lens power and the GRIN contribution, which offers a compact and simple expression for the ABDC matrix. We then apply this formulation to the crystalline lens and implement both numerical and analytical integration procedures to obtain the GRIN lens power. The analytical approximation provides an accurate solution in terms of Gaussian hypergeometric functions. Last, we compare our numerical and analytical procedures with published ABCD matrix methods in the literature, and analyze the effect of the iso-indicial surface's conic constant (Q) and inner curvature gradient (G) on the lens power for different lens models.