## Abstract

A dispersive model of a gradient refractive index (GRIN) lens is introduced based on the idea of iso-dispersive contours. These contours have constant Abbe number and their shape is related to iso-indicial contours of the monochromatic geometry-invariant GRIN lens (GIGL) model. The chromatic GIGL model predicts the dispersion throughout the GRIN structure by using the dispersion curves of the surface and the center of the lens. The analytical approach for paraxial ray tracing and the monochromatic aberration calculations used in the GIGL model is employed here to derive closed-form expressions for the axial and lateral color coefficients of the lens. Expressions for equivalent refractive index and the equivalent Abbe number of the homogeneous equivalent lens are also presented and new aspects of the chromatic aberration change due to aging are discussed. The key derivations and explanations of the GRIN lens optical properties are accompanied with numerical examples for the human and animal eye GRIN lenses.

© 2012 Optical Society of America

## 1. Introduction

The experimental studies have shown a spatial change in chromatic dispersion of the gradient index (GRIN) eye lens. In spite of a noticeable variety in human GRIN lens models, modeling the dispersion of the GRIN lens has been considered only in two studies, [1] and [2]. However, the existing chromatic models of the GRIN lens do not have enough flexibility to be consistent with the experimental data that shows different amount of dispersion at the surface and the center of the GRIN lens. To understand the origin of chromatic effects in the GRIN lens and the corresponding aberrations arising from the spatial change of dispersion within the lens, one should consider the chromatic contribution of different layers in the GRIN structure. For doing this, we need at least to develop a paraxial ray tracing method (preferably analytical) for predicting the ray path within GRIN structure at different wavelengths. The paraxial ray-tracing will be the basis for calculating chromatic aberration of the GRIN structure, provided that the dispersion model for the surface and the center of the GRIN lens is given by experimental measurements.

In this paper we employ the geometry-invariant GRIN lens (GIGL) monochromatic model [3] and introduce wavelength dependence of the refractive index. This allows us to obtain a chromatic model matching experimental data on dispersion of the GRIN lens as well as to retain all properties of the GIGL mode including the analytical description for paraxial ray-tracing. Figure 1 depicts the GIGL geometry and its interior geometry-invariant iso-indicial contours, where *R* is the external radius of curvature and *k* is the external conic constant of the lens surfaces and subscripts ‘*a*’ and ‘*p*’ designate the anterior and the posterior surfaces, respectively, so that *T _{a}* and

*T*are the axial thicknesses of the anterior and the posterior hemispheres, respectively. The sign convention used here is such that the radii and thicknesses are considered to be always positive as in our previous work [3]; however to be consistent with the optical design sign convention, we introduce ‘–’ in front of

_{p}*T*and

_{a}*R*.

_{p}## 2. Dispersion model

#### 2.1. Characteristics of dispersive GRIN lens

In the geometry-invariant GRIN lens model, the refractive index distribution is based on the power law profile, which was originally proposed by Pierscionek [4] and later supported by several studies [5–9]. In the GIGL model [3], this GRIN profile is presented as:

where the parameter*p*in the exponent is employed to account for age-related dependence of the GRIN lens,

*ζ*is the normalized distance from the center of the lens, and

*n*and

_{c}*n*are the refractive indices at the center and at the surface of the GRIN lens, respectively. Here,

_{s}*ζ*varies between −1 to +1 to cover both anterior and posterior hemispheres of the lens. Note that the exponent

*p*is not limited to integer numbers, so to avoid complex numbers we use the form (

*ζ*

^{2})

*. In monochromatic aberration studies,*

^{p}*n*and

_{c}*n*are measured at a certain wavelength (typically at

_{s}*λ*= 555 nm) for which the eye shows its highest sensitivity. It is well known that the GRIN lens is a dispersive medium with different dispersion characteristics at the center and at the surface [10,11]. In other words, the GRIN refractive index profile is different at different wavelengths along the

_{m}*Z*axis. We assume that these axial GRIN distribution profiles follow the power law with their own wavelength specific

*n*(

_{c}*λ*),

*n*(

_{s}*λ*) and

*p*(

*λ*). In view of the well-known concept of iso-indicial contours in the GRIN lens, we propose the idea of iso-dispersive contours. We consider the lens structure consisting of very thin shells whit a constant chromatic dispersion (constant Abbe number). An iso-dispersive contour is the interface between the two adjacent shells.

There are several papers on suggesting theoretically sound equations for chromatic dispersion. Atchison *et al.* [12] have studied and reviewed experimental and theoretical data on ocular media dispersion and found the Cauchy’s equation as the best fit:

*λ*is the main wavelength and

_{m}*n*is the refractive index at

_{m}*λ*.

_{m}*n*,

_{m}*n*

_{λ}_{2}and

*n*

_{λ}_{4}could be found by fitting Eq. (3) to a given dispersion data. Considering the experimental error in one of the dispersion measurements (Ref. [10]) used in the present paper, having only three terms in Eq. (3) provides an acceptable fit. In our experience, employing higher terms of Cauchy’s equation gives the fitted curve a freedom to follow the noise in the data and even to develop a minimum, which does not correspond to a theoretically valid dispersion function.

Equation (1) can be rewritten to describe a dispersive GRIN medium as

*n*and

_{c}*n*are respectively the refractive indices at the center and at the surface of the GRIN lens at the main wavelength

_{s}*λ*, since Eq. (4) would be reduced to Eq. (1) where

_{m}*λ*=

*λ*. Therefore any characteristic of the lens (

_{m}*e.g.*the optical power equations and the third-order aberration representations) defined for the monochromatic GIGL model [3, 13] will remain unchanged when using this representation.

To our knowledge this is the first attempt to represent the dispersive nature of the GRIN lens structure in terms of wavelength-dependent *n _{center}*(

*λ*),

*n*(

_{surface}*λ*), and

*p*(

*λ*). In the following section we show two numerical examples that emphasize the advantage of using this model when experimental dispersion data for the center and the surface of the lens is available.

#### 2.2. Numerical examples

Equation (4) describes the gradual change in dispersion from the surface to the center of the lens. There are not so many published experimental data on the spatial change in dispersion of the GRIN lens in the human. These available data are limited to the dispersion curves at the center and the surface of the lens. This does not provide enough information to determine the rate of change in GRIN profile at any specific wavelength, which corresponds to the exponent *p*(*λ*). Due to the lack of experimental data on wavelength dependence of the exponent *p*(*λ*), for now, we limit our examples to the case of constant *p*. However the presented model is capable of taking into account the wavelength dependence of *p*(*λ*) when more complete data become available. In the following, two sets of data are used to determine all the coefficients in the dispersive GRIN structure given by Eq. (4).

Palmer and Sivak [10] have done a series of measurements on a 70 year old eye for the wavelength range from 410 nm to 680 nm. It is worth mentioning that Palmer and Sivak have stated that the material of the lens has not been altogether homogeneous and there have been inclusions of lower refractive index. Due to this, their data may not be as reliable as other data from human lenses. Figure 2 depicts their data and the least square fit using Eqs. (5) and (6). The fitting equations are

In another work, Sivak and Mandelman [11] measured the GRIN lens dispersion of several subjects from 16-year old to 78-year old. They have provided averaged data for the dispersion curves of the center and the surface of the lens. Figure 3 shows the dispersion curves from the least square fit using Eqs. (5) and (6) to Sivak and Mandelman’s data averaged over 6 to 9 separate eyes given by

In addition to this, Fig. 4 demonstrates the dispersion curves across the lens provided by Eq. (4) for *p* equal to 2.0 and 5.0. This range is chosen in relation to the fitting results by Navarro *et al.* [8], where the case *p* = 2.0 could be considered as an extreme minimum for the refractive index profile of a very young healthy eye, whereas *p* = 5.0 corresponds to an aged eye. The intermediate curves between the lens center and surface show internal dispersion of the lens material predicted by Eq. (4) for *ζ* equal to 0.1, 0.2, . . . , 0.9. The spacing between these curves indicates the power profile of the GRIN structure (characterized by exponent *p*).

Another way of distinguishing two cases (*p* = 2.0 and *p* = 5.0) is examining the refractive index profiles along the lens optical axis; these profiles are shown in Fig. 5. Figure 5 illustrates an age-related alteration originally suggested by Pierscionek [4] in the *lens paradox* explanation.

## 3. Chromatic aberration

#### 3.1. Chromatic coefficients

Using the idea of thin, iso-dispersive shells, we shall describe the chromatic effects occurring within the bulk of the lens. For this we need to revisit the definition of the axial and transverse chromatic aberrations in terms of paraxial optics.

The axial color coefficient for a single refractive surface *C _{L}* [14] is given by

*n*is the refractive index of the medium at the spectral line

_{d}*d*(

*λ*= 587.6 nm) before the surface,

_{d}*y*and

*u*are the height and the angle of the incident ray at the surface, respectively,

*r*is the surface radius of curvature, and

*V*and

_{d}*V*′

*are the Abbe numbers of the media respectively before and after the surface. The Abbe number is defined as where*

_{d}*n*and

_{F}*n*are the refractive indices of the medium at the spectral line

_{c}*F*(

*λ*= 486.1 nm) and the spectral line

_{F}*C*(

*λ*= 656.3 nm), respectively.

_{C}By adding the *C _{L}* coefficients of all surfaces in an optical system the total axial color coefficient is obtained. Using this coefficient, the axial distance between the image position at

*F*and

*C*wavelengths can be calculated as longitudinal axial chromatic aberration

*δ*,

_{AX}*n*is the refractive index at

_{di}*d*line in the last medium and

*u*is the refracted ray angle at the image plane.

_{i}In a similar way, the lateral color coefficient *C _{T}* is given by

*y*and

_{c}*u*are the height and the angle of the incident chief (principal) ray at a refractive surface, respectively. The total coefficient

_{c}*C*can be used in calculating the transverse lateral chromatic aberration

_{T}*δ*, which corresponds to the vertical distance between the images at

_{TLC}*F*and

*C*wavelengths:

The total axial and lateral color coefficients of a GRIN lens are the key to understanding its chromatic behavior in the eye. Using the paraxial ray tracing as in our derivation of the monochromatic aberrations in the GIGL model [3], the total chromatic coefficients can be obtained for the dispersive GRIN structure. Following the analytical ray tracing through the GIGL model, we shall derive the expressions for the axial and lateral color coefficients. First we define the Abbe number of the iso-dispersive shells Using Eq. (4)

*ζ*corresponds to the normalized distance from the lens center to any point in the lens. The distance

*ζ*along the optical axis is simply

*z*/

*T*and

_{a}*z*/

*T*for the anterior and posterior hemispheres, respectively. For a constant Abbe number

_{p}*V*Eq. (16) defines the corresponding constant

_{d}*ζ*of the iso-dispersive contour.

Equation (16) can be described in terms of axial distance z for paraxial ray tracing. Employing Eqs. (11), (4), and (16) we find the contribution of a thin layer with the thickness of *δz* within the GRIN structure in the total axial color coefficient of the lens as

*y*(

*z*) is the height of the ray inside the GRIN lens and

*R*is the radius of curvature of the surface; the analytical expression for

*y*(

*z*) in the GIGL model is derived in Ref. [3]. By summing up the contributions from all infinitely thin layers in the GRIN structure (where

*R*equals to

*R*and

_{a}*R*in the anterior and posterior hemispheres, respectively) and including the contributions from the external surfaces of the lens we find the total axial color coefficient of the lens,

_{p}*y*

_{0}and

*u*are respectively the height and angle of the marginal ray at the anterior surface,

_{a}*y*(

*T*) and

_{p}*u*(

*T*) are respectively the height and the angle of the marginal ray just before the posterior surface,

_{p}*n*and

_{aqu}*V*are respectively the refractive index and the Abbe number of the medium before the lens at the

_{aqu}*d*line, and

*n*and

_{vit}*V*correspond to the medium after the lens. It is worth noticing that the ray tracing for these calculations should be done at the

_{vit}*d*line.

Similarly the coefficients for the lateral color of the lens is given by

*y*(

_{c}*z*) is the chief ray height inside the lens. Assuming the aperture stop is located at the anterior surface of the lens we obtain the total lateral color coefficient of the lens

*u*is the angle of the chief ray at the anterior surface, and

_{ca}*y*(

_{c}*T*) and

_{p}*u*(

_{c}*T*) are respectively the height and the angle of the chief ray just before the posterior surface of the lens.

_{p}At this stage we have all equations ready for analysis of chromatic effects in the GRIN lens.

#### 3.2. Numerical example

Using Eqs. (18) and (20) we can calculate the chromatic effects arising from a typical GRIN lens. To use these two equations, one needs the dispersion data for the GRIN lens as well as the dispersion data for the medium surrounding the lens. For the media before and after the lens, the aqueous and vitreous, respectively, we have used combined data provided by Atchison and Smith in terms of coefficients for Cauchy equation presented in Table 5 in Ref. [12].

### 3.2.1. Sivak and Mandelman’s experimental data

As an example, Table 1 presents a typical GRIN lens geometry and its surrounding media with dispersive characteristics. Since the eye GRIN lens receives a converging beam we shall take into account its convergence by assuming a typical corneal shape with the anterior and posterior radii of curvature of 7.8 mm and 6.7 mm, respectively, the corneal thickness of 0.5 mm, and the 3.5 mm axial distance to the GRIN lens. This corresponds to the marginal ray angle *u _{a}* = −0.036352 rad at the pupil such that

*y*

_{0}= 1 mm. The aperture stop (iris) is located just before the anterior surface of the lens, so the height of the chief ray at the anterior surface of the lens is zero. For the object located at infinity and the full field of view of 2 deg the chief ray angle at the anterior surface of the lens is

*u*= 0.015000 rad.

_{ca}Table 2 presents the chromatic coefficients of this typical GRIN lens for each part of the lens using the dispersion data from Sivak and Mandelman [11] in Eqs. (9) and (10). Here to describe the dispersive GRIN medium in accordance with Eq. (11) we assumed *p* = 4.0, which impacts the magnitude of the coefficients *C _{L}* and

*C*for the GRIN structure.

_{T}Axial color coefficients in Table 2 show an interesting compensation effect in axial chromatic aberration of the lens. The negative *C _{L}* arising from the peripheral surfaces is to some extent corrected by the positive

*C*originating from the GRIN structure of the lens. To understand the origin of the sign in coefficient

_{L}*C*we revisit Eq. (11). Assuming

_{L}*u*≪

*y*/

*r*and noting that 0 <

*y*, it is clear that the radius of curvature

*r*and the variation of the quantity (

*n*− 1)/(

_{d}*n*) play the main role in establishing the sign of the axial color coefficient.

_{d}V_{d}Figure 6 illustrates the variation of the quantity (*n _{d}* − 1)/(

*n*) inside the lens is positive for the anterior hemisphere and negative for the posterior hemisphere. On the other hand, the radius of curvature

_{d}V_{d}*r*for the iso-indicial contours in the GRIN structure is also positive in the anterior and negative in the posterior hemispheres, thus the overall sign of the axial color coefficient of the GRIN structure is positive in both hemispheres.

In contrast, the quantity (*n _{d}* − 1)/(

*n*) in the medium before and after the lens (0.004979 and 0.004888 respectively) is less than that of the surface of the lens (0.010015). This leads to the negative axial color coefficients at both anterior and posterior lens surfaces, see Table 2. A similar approach could be used to understand the change of sign in lateral color coefficients of the lens, although for the posterior surface we have an additional contribution to the coefficient from the non-zero height of the chief ray.

_{d}V_{d}To get a feeling about the magnitude of calculated axial color coefficient, we also calculate the axial chromatic aberration *δ _{AX}* in the whole eye using this typical lens. We assume that the lens in the eye is not tilted and then the total axial color coefficient of the eye is the sum of the axial color coefficients of the cornea and the GRIN lens. We assume that the optical power of an average eye is 60 D and the corneal refractive index and Abbe number are

*n*= 1.3677 and

_{d}*V*= 55.48 [12], which corresponds to the axial color coefficient of −0.000843 for the cornea. The sum of the axial coefficients of the cornea and the lens is −0.001718, and from Eq. (13) we find that the axial chromatic aberration of this eye is 2.09 D. This amount of chromatic aberration is within the expected range of 1.87 ± 0.26 D found in the study by Gilmartin and Hogan [15], which used a similar spectral range of 488 nm and 633 nm.

_{d}In addition to the lens dispersion data, Sivak and Mandelman have also provided the dispersion of the capsule of the lens. To study the effect of the lens capsule dispersion on the lens axial color coefficient we have performed a numerical fit to Sivak and Mandelmans measurements using Eq. (3) for the capsule

### 3.2.2. Palmer and Sivak’s experimental data

Employing the same lens geometry and the same dispersive characteristics of the surrounding media, we shall calculate the axial color coefficients of the lens using the dispersion data from the study by Palmer and Sivak of a 70 year old eye [10]. In this study the age of the eye is well defined, so we consider the exponent *p* = 5.5 based on the equation suggested by Navarro *et al.* [8] for connecting the exponent *p* to the age of the eye. Table 4 presents the color coefficients of the lens and Fig. 7 depicts the quantity (*n _{d}* − 1)/(

*n*) as a function of the normalized distance

_{d}V_{d}*ζ*from the lens center.

In this example we can also see the compensation of chromatic aberrations in the lens, however in comparison to Table 2 the compensation happens in a different way. Table 4 indicates that the GRIN structure of the lens produces a negative axial color coefficient, while the axial color coefficients of the external surfaces are positive. This is due to the different rate of growth of the quantity (*n _{d}* − 1)/(

*n*) presented in Fig. 7, and also due to the different amount of the quantity (

_{d}V_{d}*n*− 1)/(

_{d}*n*) at the lens surface, 0.003328, which is less than the surrounding media. As a result, the total axial color coefficient of the lens is still negative (−0.000537) in comparison with Table 2, but its absolute value is reduced by 20%. Using the same characteristics for the cornea, as in Section 3.2.1, the chromatic aberration of the whole eye becomes 1.68 D. This amount of chromatic aberration is within the expected range of 1.87 ± 0.26 D found in the study by Gilmartin and Hogan [15] using a similar spectral range.

_{d}V_{d}### 3.2.3. Theoretical equations

In addition to the experimental data used in our two examples, it is worth revisiting the lens dispersion equations suggested by Atchison and Smith (Table 5 in Ref. [12]), which is based on the combined theoretical data from Le Grand [17] and Navarro *et al.* [18]. One can notice that the two equations describing the center and the surface dispersion can be reduced to the simple expression *n _{center}*(

*λ*) ≃

*n*(

_{surface}*λ*) × 1.014430. For such a simple connection between the center and the surface dispersion curves, the quantity (

*n*− 1)/(

_{d}*n*) becomes independent of

_{d}V_{d}*ζ*terms and remains constant (0.006127) inside the GRIN lens. This becomes clear when we rewrite the quantity (

*n*− 1)/(

_{d}*n*) as (

_{d}V_{d}*n*−

_{F}*n*)/

_{c}*n*. The significance of the constancy of the quantity (

_{d}*n*− 1)/(

_{d}*n*) is that the GRIN structure is free from both chromatic effects (axial and lateral color) and only the external surfaces of the lens contribute to the chromatic aberrations. Table 5 shows this fact numerically under the same conditions considered in Table 1.

_{d}V_{d}We would like to emphasize that for a non-dispersive GRIN structure as well as for a homogeneous refractive index lens (featuring in Le Grand [17] and Navarro *et al.* [18] models) one has only to worry about the dispersion description of the external surfaces of the lens. However the experimental data [10, 11] suggest that the chromatic contribution of the GRIN structure is comparable to that of the external surfaces, and furthermore plays an important role in the chromatic aberration compensation of the whole eye.

It is worth mentioning that for achieving an achromatic lens (*i.e.* ∑*C _{L}* = 0), the axial color of the external surfaces is to be compensated by that from the GRIN structure, so the latter should not be zero. As an example of such aberration balancing, the lens in Section 3.2 can be made achromatic simply by adjusting one of the center dispersion coefficients,

*n*

_{c}_{λ}_{2}, in Eq. (9) from 7256.06 to 5098.80. This corresponds to the change in the Abbe number of the lens center from 33.71 to 51.51, using Eq. (4). The example illustrates that the proposed dispersion model can also help in developing bio-inspired lens designs, where the achromatic correction in individual components is needed.

## 4. Equivalent Abbe number approximation

It is a common practice, in the reduced eye models, to replace the GRIN lens by a simple lens providing spatially constant refractive index and Abbe number. Typically the characteristics of the eye lens in the reduced eye models are defined based on the optical power and axial chromatic aberration of the whole eye. In this section a different approach is proposed to calculate the equivalent Abbe number of the lens just using the experimental dispersion curves of the surface layer (cortex) and the center (core) of the GRIN lens. This approach may help in improving the reconstruction methods of a subject-specific GRIN lens.

The GIGL model provides a convenient equation for the optical power of the lens derived from a thin lens approximation as

*n*provides the equivalent optical power for the lens. Equating the right hand sides of Eqs. (22) and (23) and solving for

_{eqv}*n*result

_{eqv}The definition of *n _{eqv}* by Eq. (24) is not limited to the main wavelength and can be expanded to the equivalent powers of other wavelengths as

*p*with wavelength

*λ*is negligible the equivalent Abbe number can be rewritten as

*V*and

_{dc}*V*are the Abbe numbers of the center and the surface of the GRIN lens, respectively. As an example, Table 1 presents the calculated lens equivalent refractive index

_{ds}*n*and the equivalent Abbe number

_{eqv}*V*using the dispersion data from Sivak and Mandelman [11], and Palmer and Sivak [10] for

_{eqv}*p*= 3.0 (representing a 34-year old eye according to the study by Navarro

*et al.*[8]).

Employing the equivalent refractive indices and Abbe numbers listed in Table 6, one can calculate the axial and lateral color coefficients of the equivalent lens. These coefficients are approximate, while the exact color coefficients derived in Eqs. (18) and (20). The comparison between the axial and lateral color of a typical GRIN lens and its equivalent refractive index lens is presented in Table 7. Here, the angles and the height of the incident rays, and the geometry of the lens and its surrounding media are the same as in the previous numerical examples (See Section 3.2.1), except for the exponent *p*, here *p* = 3.0.

Table 7 demonstrate a good agreement between the exact and equivalent color coefficients, which indicates practical advantages of this approach in finding the exact equivalent optical characteristics. However, replacing a complicated GRIN structure with a homogeneous material will not transfer all optical characteristics of the lens simultaneously. Here, the equivalent refractive index calculation aims at preserving the optical power of the lens, yet other optical characteristics of the lens need their own equivalent refractive index. As an example, we consider the optical path length (OPL) in the GRIN lens. The axial OPL in GIGL [3] is derived as

which defines the quantity (2*n*+

_{c}p*n*)/(2

_{s}*p*+ 1) as the equivalent refractive index in the OPL calculation of the lens. This OPL equivalent refractive index known as the average refractive index [19] is different from the one calculated in Eq. (24).

## 5. Axial Chromatic aberration and aging

There are two independent studies, which found no significant changes in the magnitude of chromatic aberration with age [20, 21], but some earlier studies claim that the magnitude of chromatic aberration decreases with age [22, 23]. Considering the growth of the GRIN lens, it is well known that for an unaccommodated lens, the external surfaces become more curved with aging [24]. The lens considered in Section 3.2.2 shows positive axial color coefficients at its surfaces (Table 4). Aging lens increases the optical power of its surfaces, thus this lens surfaces will show larger positive axial color coefficients, see Eq. (11). In addition to this, the chromatic effect arising from the GRIN structure should be taken into account. To examine the role of the GRIN structure in the age-related chromatic effect individually, we keep the geometry of the lens in Section 3.2.2 unchanged, only adjust the age-related exponent *p* to three different age groups. Following the study by Navarro *et al.* [8] we selected three age groups (20, 40, and 60 year old) with corresponding value for *p* and calculated the resultant axial color coefficients for each group in Table 8.

Table 8 indicates that the magnitude of the axial color coefficient originated from the GRIN structure of the lens decreases due to an increase in the exponent *p*. Looking at the sign of the coefficients of the lens surfaces, it is evident that the total axial color coefficient of the lens goes toward zero or even positive amounts. As mentioned in Section 3.2.1, a typical cornea shows a negative axial color coefficient. Due to aging, the surfaces of the cornea get slightly more curved [25], which increases the magnitude of this negative corneal axial color coefficient. Considering these two effects for the GRIN lens and the cornea, one could argue that the eye, as a whole, might show a constant or a better chromatic performance with aging.

It is worth mentioning that employing the Sivak and Mandelman’s data in the same approach provides an increase in the total axial color coefficient, which does not support the argument above. Since the dispersion data from Sivak and Mandelman is averaged over 6 to 9 subjects with different ages, some valuable information of each individual lens might be missing. This highlights the need of more chromatic measurements on the GRIN lens and its surrounding media for a more confident conclusion.

## 6. Discussion and conclusion

In this study the potential of the existing GIGL monochromatic model [3] has been realized and demonstrated by introducing a flexible chromatic GRIN lens model. The model is examined with different data fittings to gain better understanding of the GRIN lens chromatic aberration behavior of the lens. The advantages of the provided equations are not limited to this data fitting only, since the model can also be used to help explaining other chromatic characteristics of the eye for example spherochromatism. The latter can be calculated using the spherical aberration coefficients for different color wavelength.

We would like to emphasize that the provided model should be regarded as a tool, which can be employed in experimental data fitting and further analysis for better understanding of the chromatic nature of the eye. The scope of the chromatic model is not limited to the human eye and can be applied to animal GRIN lenses, where unusual chromatic behavior might happen (*e.g.* [26,27]). Table 9 provides the axial color coefficients of a semi-spherical octopus eye lens calculated based on the measurements from Jagger and Sands [27].

The proposed dispersive GIGL model completes the description of the monochromatic GIGL model. The concept of iso-dispersive contours introduced in this paper is a unique feature of the dispersive GIGL model. The optically-friendly geometry of the GIGL model featuring iso-dispersive contours (constant Abbe number *V _{d}*) supports the calculations for the contribution of the individual layers to the total chromatic effects of the lens. Thus the model gives a new insight in the derivation of the equivalent Abbe number based on experimental dispersion curves for the center and the surface of the lens. The proposed model offers easy analysis of equivalent refractive index and OPL equivalent index. One could use the model to investigate chromatic aberration evolution in the aging human eye, in particular the effect of the power constant p on the color coefficients. The model also predicts the amount of induced defocus (change in optical power with wavelength) when using different wavelengths,

*e.g.*adaptive optic system with wavefront sensor operating at one wavelength and the retinal imaging at another wavelength, although retinal absorption effects still have to be included in this case.

## Acknowledgments

The authors would like to thank Prof. Chris Dainty for his valuable comments. This research was supported by Science Foundation Ireland under grant 07/IN.1/1906.

## References and links

**1. **H. L. Liou and N. A. Brennan, “Anatomically accurate, finite model eye for optical modeling,” J. Opt. Soc. Am. A **14**, 1684–1695 (1997). [CrossRef]

**2. **J. A. Díaz, C. Pizarro, and J. Arasa, “Single dispersive gradient-index profile for the aging human lens,” J. Opt. Soc. Am. A **25**, 250–261 (2008). [CrossRef]

**3. **M. Bahrami and A. V. Goncharov, “Geometry-invariant GRIN lens: analytical ray tracing” J. Biomed. Opt. **17**, 055001 (2012). [CrossRef] [PubMed]

**4. **B. K. Pierscionek, “Presbyopia - effect of refractive index,” Clin. Exp. Optom. **73**, 23–30 (1990). [CrossRef]

**5. **G. Smith, D. A. Atchison, and B. K. Pierscionek, “Modeling the power of the aging human eye,” J. Opt. Soc. Am. A **9**, 2111–2117 (1992). [CrossRef] [PubMed]

**6. **G. Smith and B. K. Pierscionek, “The optical structure of the lens and its contribution to the refractive status of the eye,” Ophthalmic Physiol. Opt. **18**, 21–29 (1998). [CrossRef] [PubMed]

**7. **C. E. Jones, D. A. Atchison, R. Meder, and J. M. Pope, “Refractive index distribution and optical properties of the isolated human lens measured using magnetic resonance imaging (MRI),” Vision Res. **45**, 2352–2366 (2005). [CrossRef] [PubMed]

**8. **R. Navarro, F. Palos, and L. González, “Adaptive model of the gradient index of the human lens. I. formulation and model of aging ex vivo lenses,” J. Opt. Soc. Am. A **24**, 2175–2185 (2007). [CrossRef]

**9. **S. Kasthurirangan, E. L. Markwell, D. A. Atchison, and J. M. Pope, “In vivo study of changes in refractive index distribution in the human crystalline lens with age and accommodation,” Invest. Ophthalmol. Vis. Sci. **49**, 2531–2540 (2008). [CrossRef] [PubMed]

**10. **D. A. Palmer and J. Sivak, “Crystalline lens dispersion,” J. Opt. Soc. Am. A **71**, 780–782 (1981). [CrossRef]

**11. **J. G. Sivak and T. Mandelman, “Chromatic dispersion of the ocular media,” Vision Res. **22**, 997–1003 (1982). [CrossRef] [PubMed]

**12. **D. A. Atchison and G. Smith, “Chromatic dispersions of the ocular media of human eyes,” J. Opt. Soc. Am. A **22**, 29–37 (2005). [CrossRef]

**13. **The geometry-Invariant lens computational code. This is a computable document format (CDF) for the equations presented in Ref. [3]. Our source CDF code can be accessed via Mathematica, the computational software developed by Wolfram Research (Oct. 2011), *http://optics.nuigalway.ie/people/mehdiB/CDF.html*.

**14. **M. J. Kidger, *Fundamental Optical Design* (SPIE Press, 2002).

**15. **B. Gilmartin and R. E. Hogan, “The magnitude of longitudinal chromatic aberration of the human eye between 458 and 633 nm,” Vision Res. **25**, 1747–1755 (1985). [CrossRef] [PubMed]

**16. **R. I. Barraquer, R. Michael, R. Abreu, J. Lamarca, and F. Tresserra, “Human lens capsule thickness as a function of age and location along the sagittal lens perimeter,” Invest. Ophthalmol. Vis. Sci. **47**, 2053–2060 (2006). [CrossRef] [PubMed]

**17. **Y. Le Grand, *Form and Space Vision*, rev. ed., translated by M. Millodot and G. Heath (Indiana University Press, 1967).

**18. **R. Navarro, J. Santamaría, and J. Bescós, “Accommodation dependent model of the human eye with aspherics,” J. Opt. Soc. Am. A **2**, 1273–1281 (1985). [CrossRef] [PubMed]

**19. **S. R. Uhlhorn, D. Borja, F. Manns, and J. M. Parel, “Refractive index measurement of the isolated crystalline lens using optical coherence tomography,” Vision Res. **48**, 2732–2738 (2008). [CrossRef] [PubMed]

**20. **C. Ware, “Human axial chromatic aberration found not to decline with age,” A. Graefes Arch. Klin. Exp. Ophthalmol. **218**, 39–41 (1982). [CrossRef]

**21. **P. A. Howarth, X. X. Zhang, D. L. Still, and L. N. Thibos, “Does the chromatic aberration of the eye vary with age?,” J. Opt. Soc. Am. A **5**, 2087–2096 (1988). [CrossRef] [PubMed]

**22. **M. Millodot, “The influence of age on the chromatic aberration of the eye,” A. Graefes Arch. Klin. Exp. Ophthalmol. **198**, 235–243 (1976). [CrossRef]

**23. **J. A. Mordi and W. K. Adrian, “Influence of age on chromatic aberration of the human eye,” A. Graefes Arch. Klin. Exp. Ophthalmol. **198**, 235–243 (1976). [CrossRef]

**24. **N. Brown, “The change in lens curvature with age,” Exp. Eye Res. **19**, 175–183 (1974). [CrossRef] [PubMed]

**25. **A. V. Goncharov and C. Dainty, “Wide-field schematic eye models with gradient-index lens,” J. Opt. Soc. Am. A **24**, 2157–2174 (2007). [CrossRef]

**26. **R. H. H. Kröger and M. C. W. Campbell, “Dispersion and longitudinal chromatic aberration of the crystalline lens of the African cichlid fish Haplochromis burtoni,” J. Opt. Soc. Am. A **13**, 2341–2347 (1996). [CrossRef]

**27. **W. S. Jagger and P. J. Standsl, “A wide-angle gradient index optical model of the crystalline lens and eye of the octopus,” Vision Res. **39**, 2841–2852 (1999). [CrossRef] [PubMed]