Journal of Optometry Journal of Optometry
J Optom 2018;11:10-20 - Vol. 11 Num.1 DOI: 10.1016/j.optom.2017.04.003
Original Article
Extended depth of focus contact lenses vs. two commercial multifocals: Part 1. Optical performance evaluation via computed through-focus retinal image quality metrics
Lentes de contacto de profundidad de foco extendido vs. dos lentes multifocales comerciales: Parte 1. Evaluación del rendimiento óptico del enfoque computarizado mediante métricas de calidad de imagen retinianas
Ravi C. Bakarajua,b,, , Klaus Ehrmanna,b, Arthur Hoa,b
a Brien Holden Vision Institute, Sydney, Australia
b School of Optometry and Vision Science, University of New South Wales, Sydney, Australia
Received 01 November 2016, Accepted 21 April 2017

To compare the computed optical performance of prototype lenses designed using deliberate manipulation of higher-order spherical aberrations to extend depth-of-focus (EDOF) with two commercial multifocals.


Emmetropic, presbyopic, schematic eyes were coupled with prototype EDOF and commercial multifocal lenses (Acuvue Oasys for presbyopia, AOP, Johnson & Johnson & Air Optix Aqua multifocal, AOMF, Alcon). For each test configuration, the through-focus retinal image quality (TFRIQ) values were computed over 21 vergences, ranging from −0.50 to 2.00D, in 0.125D steps. Analysis was performed considering eyes with three different inherent aberration profiles: five different pupils and five different lens decentration levels.


Except the LOW design, the AOP lenses offered ‘bifocal’ like TFRIQ performance. Lens performance was relatively independent to pupil and aberrations but not centration. Contrastingly, AOMF demonstrated distance centric performance, most dominant in LOW followed by MED and HIGH designs. AOMF lenses were the most sensitive to pupil, aberrations and centration. The prototypes demonstrated a ‘lift-off’ in the TFRIQ performance, particularly at intermediate and near, without trading performance at distance. When compared with AOP and AOMF, EDOF lenses demonstrated reduced sensitivity to pupil, aberrations and centration.


With the through focus retinal image quality as the gauge of optical performance, we demonstrated that the prototype EDOF designs were less susceptible to variations in pupil, inherent ocular aberrations and decentration, compared to the commercial designs. To ascertain whether these incremental improvements translate to a clinically palpable outcome requires investigation through human trials.


Comparar el rendimiento óptico computado de los prototipos de lentes de contacto diseñadas mediante manipulación deliberada de aberraciones esféricas de alto orden para conseguir la profundidad de foco extendida (EDOF) con dos lentes multifocales comerciales.


Se adaptaron lentes de contacto con prototipo EDOF, y lentes multifocales comerciales a ojos esquemáticos présbitas emétropes (Acuvue Oasys para presbicia, AOP, Johnson & Johnson & Air Optix Aqua multifocal, AOMF, Alcon). Para cada configuración de la prueba se obtuvieron los valores de la calidad de imagen de la retina a través del foco (TFRIQ) en 21 vergencias, que oscilaron entre -0,5 y 2D, en pasos de 0,125D. Se llevó a cabo el análisis considerando los ojos con tres perfiles de aberración diferentes: cinco pupilas diferentes y cinco niveles de descentramiento de lentes distintos.


Exceptuando el diseño LOW, las lentes AOP aportaron un rendimiento TFRIQ de tipo ‘bifocal’. El desempeño de las lentes fue relativamente independiente de la pupila y aberraciones, pero no del descentramiento. Por contra, AOMF demostró un rendimiento céntrico en visión lejana, más dominante en el diseño LOW, seguido de los diseños MED y HIGH. Las lentes AOMF fueron más sensibles a la pupila, aberraciones y centrado. Los prototipos demostraron un efecto de mejora en el rendimiento TFRIQ, particularmente en visión intermedia y próxima, sin pérdida de rendimiento en visión lejana. Al compararlas con las lentes AOP y AOMF, las lentes EDOF demostraron una reducción de la sensibilidad al factor pupila, las aberraciones y el centrado.


Con la calidad de imagen de la retina a través del foco, como medida del rendimiento óptico, demostramos que los diseños del prototipo EDOF fueron menos susceptibles a las variaciones de la pupila, las aberraciones oculares inherentes, y al descentramiento, en comparación a los diseños comerciales. Se precisa más investigación futura, mediante ensayos en humanos, para verificar si estas mejoras incrementales se trasladan a resultados clínicamente evaluables.

Presbyopia, Higher order aberrations, Extended depth of focus contact lenses, Ray tracing simulation, Fourier Optics
Palabras clave
Presbicia, Aberraciones de alto orden, Lentes de contacto de profundidad de campo ampliada, Simulación de trazado de rayos, Óptica de Fourier

Presbyopia currently affects more than 20% of the world's population.1 Whether it is due to the cosmetic advantages that contact lenses (CLs) offer over spectacles,2 a more demanding lifestyle in this digital era,3 increased practitioner confidence to fit newer multifocal CLs (MFCLs) lenses,4 or a combination of the above, uptake of presbyopic MFCLs appears to have gained momentum in recent years.5–7

Current simultaneous image MFCLs can be broadly classified into three types: center-distance, center-near and concentric ring bifocals. The power distribution as a function of the optic zone diameter in the former two designs is monotonic (i.e. there is a smooth transition of power from distance to near zone or vice versa); while the latter design uses a non-monotonic power distribution, which is often classified as a zonal refractive design.8,9 The merits and demerits of visual performance with such conventional optical designs have been thoroughly explored.4,10–23 Recently, our group have pilot tested extended depth of focus (EDOF) designs, achieved through manipulation of higher order spherical aberration terms using non-linear, optimization algorithms that aimed to minimize the merit function defined by through focus retinal image quality (TFRIQ).24 Preliminary investigations on the visual performance in presbyopes fitted with EDOF CLs after a short term (1h) exposure have been reported.25,26 However, a detailed theoretical investigation on how the performance of the proposed EDOF CLs differ from the commercial designs like center-near designs and concentric bifocals has not been reported.

There are several proposed methods that enable vision scientists and lens designers to gauge the optical performance of simultaneous image MFCLs.27 These methods include schematic optical modeling using ray tracing and in vitro testing of lenses on a bench-top, life-sized, physical model eye.27,28 In this paper, the first of two parts, we report a theoretical evaluation of the optical performance of a proposed EDOF CLs designed to correct presbyopia. The optical performance of the proposed designs is then compared to the modeled performance of two commercial multifocal alternatives: ACUVUE® OASYS® for Presbyopia (AOP, Johnson & Johnson, USA) and AIR OPTIX® Aqua Multifocal (AOMF, Alcon, USA). To emulate real world performance, we aim to evaluate the TFRIQ calculations on a schematic model eye at various pupils, inherent aberrations and decentration levels. In part 2 of the paper, we report data on presbyopic visual performance with proposed EDOF CLs compared to commercial alternatives, obtained from dispensing clinical studies.

MethodsLens designs (power profiles)

Fig. 1 shows the power profiles of all the study lenses, described in detail elsewhere.9,26,29 The AOMF lenses have a monotonic power distribution across the optic zone radius. On the other hand, the AOP lens is a concentric-ring zonal bifocal. In contrast to commercials, the EDOF CLs have series of smooth, non-monotonic, aperiodic, power variation across the optic zone that was designed using deliberate manipulation of multiple HOSA terms. Here, the objective was to achieve an extension of depth of focus that would facilitate a balanced visual performance across a range of predefined dioptric distances. The profiles in each column of Fig. 1 are designed to satisfy presbyopes with LOW, MED and HIGH-add requirements.

Figure 1.

Power profiles of the study lenses. Units of the color scale is in Diopters (D). L: Low; M: Medium; H: High add power.

Transformation of lens designs (power profiles) into Zemax compatible surfaces

The optical designs of the study lenses (described in Fig. 1) were obtained from the physical lens measurements performed using NIMO TR1504 (Lambda-X, Nivelles, Belgium).9 To model the commercial and prototype designs, the optical designs were transformed to appropriate surface shapes in Zemax ray tracing software (Optic Studio 16 SP1, professional version) before proceeding with the image quality computations. This transformation was accomplished using a custom-written application that routed Zemax through Matlab (MathWorks® 2015b, USA) via dynamic data exchange interface. This program offers the user a graphical user interface to elect a base curve (defined as a standard conic section), central thickness, refractive index, and a target power profile. To achieve the target power profile, the in-built damped least-squared optimization routine of Zemax was employed to facilitate the optimization of the front surface of the correcting CL, which was described as a ‘super conic’ surface. During the optimization routine, the sagittal power across the optic zone diameter was obtained as an interim output. Whenever a local minimum was achieved, the front surface was tested for its accuracy by comparing the target and output power profiles using the summated square error metric. The optimization was deemed complete when the lenses achieved a summated square error ≤0.05D, calculated over the entire optic zone diameter.

Modeling presbyopic schematic eyes

The Navarro–Escudero30 eye was modified to mimic presbyopic eyes with no residual accommodation. These modifications were based on unaccommodated emmetropic eyes, refractive error dependent models, corneal and lens shapes and concomitant changes with the age, and in vitro lenticular refractive index (RI) measurements.27,28,31–43 The effect of residual accommodation was not taken into account. Further details of the model are reported elsewhere.27,33,44 All the ray-tracing routines were performed in Zemax (Optic Studio 16 SP1, professional version 2016) executed through Matlab (MathWorks 2015b) via dynamic data exchange interface.

Ray tracing through combination of contact lens and model eyes

Soft CLs wrap the eye in its flexed state,45,46 and consequently, the back surface of a correcting lens would then be aligned with cornea, provided the fit is acceptable. Hence, we chose to design the lenses in their flexed state, with a base curve radius of 7.75mm, asphericity equal to −0.25 and refractive index of 1.42, over an 8mm optic zone. Only the optical zone of the CLs was modeled, the peripheral carrier zone of the CL was ignored. The CLs chosen for this modeling exercise had emmetropic distance prescription whose central thickness was fixed at 0.100mm.

At first, we modified the baseline Navarro–Escudero schematic model eye to produce an aberration free model. In the interest of brevity, only an emmetropic, presbyopic model eye was considered in this experiment. However, these same methods could also be used to test presbyopic eyes with myopia or hyperopia. To avoid unnecessary complexity, tear film dynamics expected in a real eye, were ignored. The two commercial MFCLs and the prototype EDOF CLs (each of their respective LOW, MED and HIGH designs), described in earlier sections, were used to correct the emmetropic presbyopic schematic model eye. All the test lenses were well centered on the corneal surface, and ‘into-the-eye’ ray tracing method was adapted to analyze the optical performance of the test lens+model eye combinations.44 A monochromatic reference wavelength of 0.589μm was chosen for all calculations. The pupil surface was made the aperture stop, and a 512×512-ray grid on the entrance pupil was used in all calculations. The ray tracing was performed in a real, robust ray aiming, cache mode as described elsewhere.33,44 The entrance pupil was modeled to emulate the effects of the human Stiles–Crawford effect using a simple Gaussian function with a standard deviation of −0.115, as described elsewhere.47,48

Modeling effects of different pupils, inherent aberrations and decentration levels

To mimic the diversity observed in normal eyes and generalize our results to the population, we have modeled three levels of inherent aberration profiles. The distribution of inherent spherical aberration in the relaxed eyes is well characterized.49–61 The inherent spherical aberrations were modeled based on 5th, 50th and 95th percentiles of the normative distribution, which translate to approximately −0.10μm, 0.10μm and +0.30μm at 6mm pupil.44 In addition to spherical aberration, normal eyes have considerable amounts of on-axis, asymmetric, higher order aberrations (HOA, e.g. horizontal/vertical coma and trefoil).49 Therefore, to further align our models to the real world observations, we supplemented the aberration models with mean values of coma [horizontal: +0.143μm and VERTICAL: +0.09μm] and trefoil [horizontal: +0.106μm and vertical: +0.081μm] at 6mm pupil from published data.49 These three levels of inherent aberration models are referred here and elsewhere in the manuscript as Models 1, 2 and 3. The inherent ocular aberrations once modeled, were dialed into the CL+eye system via customized Matlab routines. For each test configuration, the TFRIQ metrics were computed over target vergences, ranging from −0.50D to 2.00D in 0.125D steps. To make the modeling more realistic, the entire analysis was performed for five different pupils, ranging from 3 to 5mm, in 0.5mm steps. To further mimic reality, in addition to the well-centered CL on eye scenario, the TFRIQ values were also computed for four preset CL vertical decentrations: 0.125mm, 0.250mm, 0.375mm and 0.500mm at 3.5mm pupil. When studying the effects of lens decentration, the asymmetric HOA were ignored from the three models.

Measures of the optical performance

For each configuration, the exit pupil raw wavefront (512×512 array) was numerically propagated to the Fourier plane to yield a point spread function. The modulation transfer function was computed as the real part of optical transfer function, obtained after Fourier transformation of the point spread function. The real part of the optical transfer function weighted with the neural contrast sensitivity function was calculated, within spatial frequency bounds of 0–30cycles/°. The ratio of the volumes under the product of the optical transfer function and the contrast sensitivity function calculated with a given correction to the diffraction limited case was defined as the retinal image quality (RIQ), defined as visual Strehl ratio (VSOTF), described in detail elsewhere.27,44,62–64 A natural log transformation was applied to the RIQ computations. The lower the log transformed RIQ, the poorer the optical performance of the system.

ResultsAnalysis of overall optical performance

Figs. 2–4 are presented as a montage of 9 graphs assembled into a 3-by-3 matrix. These graphs illustrate the computed TFRIQ values as a function of the following test variables: pupil (Fig. 2), inherent aberrations (Fig. 3) and lens decentration (Fig. 4), respectively.

Figure 2.

Log transformed through password-focus retinal image quality values computed as a function of 5 pupil diameters: from 3 to 5mm in 0.5mm steps. This data represents computations performed on the schematic eye embedded with Model 2 inherent aberrations. The correcting lens was well-centered over the cornea. A red dashed line across subgraphs indicates the ideal threshold RIQ. The bounded area under the perpendicular lines drawn from X-axis (absolute value) is inscribed at the corner of each sub-graphs indicate the variance of TFRIQ between the five pupil diameters.

Figure 3.

Log transformed through-focus retinal image quality values computed as a function of 3 inherent aberration models: Models 1–3. This data represents computations performed on the schematic eye with 4mm pupil. The correcting lens was well-centered over the cornea. A red dashed line across sub-graphs indicates the ideal threshold RIQ. The bounded area under the perpendicular lines drawn from X-axis (absolute value) is inscribed at the corner of each sub-graphs indicate the variance of TFRIQ between the three aberration models.

Figure 4.

Log transformed through-focus retinal image quality values computed as a function of 5 preset lens centration positions: 0.0–0.50mm in 0.125mm steps. This data represents computations performed on the schematic eye embedded with Model 2 inherent aberrations (excluding asymmetric HOA) at 3.5mm pupil diameter. A red dashed line across sub-graphs indicates the ideal threshold RIQ. The bounded area under the perpendicular lines drawn from X-axis (absolute value) is inscribed at the corner of each sub-graphs indicate the variance of TFRIQ between the five decentration levels.

The rows of the 3-by-3 matrix represent computed performance for AOP, AOMF and EDOF lenses; while the columns represent data for lens types designed to cater for participants with LOW, MED and HIGH add requirements, respectively. The section in ‘white’ ribbons on each of the sub-graphs of Figs. 2–4, highlight the intended range of vision correction for those specific study lenses. The lengths of ‘white’ ribbons for LOW, MED and HIGH correspond to vergence ranges: 0–0.75D, 0–1.50D and 0–2.00D, respectively.

Intended bifocality, multifocality or extension of depth of focus can only be achieved at an expense of loss of RIQ at an alternate distance. Successful design of an ‘ideal’ simultaneous image multifocals for a presbyope attempts to minimize the compromise at any given distance. An eye (free of HOA), when imposed with approximately ∼0.375D of defocus at 4mm pupil would yield a computed log transformed RIQ of −1.50. Given that most people can clinically tolerate up to approximately 0.25D blur without significantly affecting their visual performance, a log transformed RIQ value of −1.50 log units was selected to be the ideal threshold for this simulation. A red dotted line plotted along the horizontal dimension (anchored at the proposed ideal threshold) is offered on every graph to enable the reader to differentiate the performance between the individual lens types. To further facilitate overall appreciation of the TFRIQ, in each of the sub-graphs of Figs. 2–4, the standard deviation of the TFRIQ measure, as function of the defined variable is shown with a ‘stem’ plot (perpendicular lines drawn from X-axis). In addition, for each test case, the area (absolute value) bounding the ‘stem’ plot calculated over the defined white ribbon length (AUR), was inscribed onto the top-right section of the each sub-graph. For any given variable (e.g. pupil size, inherent aberrations and lens decentration), the lower the AUR, the greater the insensitivity of TFRIQ to the changes in that particular variable (i.e. more consistent the optical performance of the lens).

To demonstrate how the optical performance varies as a function of three independent variables, one could populate an exhaustive number of permutations of possible performance plots involving the proposed variables: pupil, inherent aberrations and lens centration. In the interest of brevity, we choose to only represent three distinct examples that closely represent population data. Firstly, to demonstrate effects of change in pupil on the through focus performance, we locked the remaining two variables. Hence, data in Fig. 2 represents the test case, when the correcting CL is well-centered on the schematic eye embedded with population mean aberration levels defined by Model 2 (which includes both symmetric and asymmetric HOA). Secondly, to demonstrate the effects of change in inherent aberration profiles of the corrected eyes on the through focus performance, we froze the pupil to 4mm and again considered a well centered correcting CL scenario (Fig. 3). Lastly, to demonstrate the effects of lens decentration on performance (Fig. 4), we choose to present data for 3.5mm pupil on an eye with aberrations defined by Model 2, however the asymmetric HOA were ignored for this step.

The information required to facilitate comparison of three study lenses was condensed into Figs. 2–4. For each figure, at first, the results are reported from the top-left, sub-graph. Later, to gauge performance of one lens type against the other, results are compared against their corresponding counterparts based on add design, along the same column. Later, this procedure is repeated across the other two columns of the figure. Two descriptive metrics are used as tools to compare results within in sub-graphs and also across the columns (add type) and rows (lens-types); (a) the TFRIQ performance comparison with the ideal scenario (red-dashed line); and (b) the standard deviation plot of the TRIQ as function of a specific test variable (numerical value, AUR). As this simulation experiment treats emmetropic presbyopia, plano distance prescription was chosen for all lenses. However, it is noticeable from Fig. 1 that the AOMF designs, particularly LOW and MED have excessive negative power that could potentially make its ‘effective’ distance prescription in the myopic direction. Therefore, to make the comparisons between the study lenses reasonable and more importantly to account for a potential over refraction with lenses that employ monotonic power profiles, the TFRIQ performance obtained with AOMF LOW and MED designs were shifted to the right by 0.25D.

As seen in Fig. 2, the TFRIQ calculations obtained with the AOP LOW design satisfies the ideal threshold of −1.5 log units for all pupils. This observation also holds partially true for the EDOF LOW design, at least for pupils ≤4mm. However, the actual RIQ values obtained with the EDOF LOW seem to be higher than those obtained with the AOP LOW. With an appropriate offset, as seen in the highlighted ‘white’ sections of the sub-graphs, AOMF LOW satisfies the threshold requirements. Considering the variability of TFRIQ values with change in pupil, the AOP LOW design seems to the best performer (AUR=1.13), followed by the EDOF LOW (AUR=2.09) and the AOMF LOW (AUR=3.53) designs. Although not critical for a LOW design, it is worthwhile to note that the TFRIQ curves obtained with the AOMF LOW and EDOF LOW seem to have a ‘lift-off’ on the right hand side of the TFRIQ curve. This would theoretically translate to better performance for near vergences. Exactly the opposite is observed with the AOP LOW design, where a dip on the right hand side of TFRIQ curve is observed.

Considering the MED designs (i.e. second column of Fig. 2), we notice that the TFRIQ curves of the AOP MED design becomes bimodal (two distinct peaks) in nature, which is a pattern clearly appreciated for pupils less than 4mm. This would theoretically translate to improved optical performance at vergence of 1.50 D, which comes at an expense of poor performance at vergences of about 0.75D. However, the bimodal pattern of the TFRIQ curves seem to fade away for pupils >4mm. Unlike the AOP MED design, the EDOF MED sits close to the ideal threshold requirement, without inducing bimodality, at least for pupils less than 4mm. For pupils >4mm, the TFRIQ curve produces an infliction point at a vergence range of 1.00D but it is not distinct enough to create a ‘bimodal’ pattern observed with AOP MED. However, it is important to notice that the near peak obtained with the EDOF MED slips marginally lower than its AOP MED counterpart. The AOMF MED design fails to fully satisfy the requirements of an ideal scenario at all pupils. Again, based on the variability measure, the AOP MED design (AUR=3.15) outperforms the EDOF MED (3.96) and AOMF MED (6.26) designs.

Next, we examine the results obtained with the HIGH designs (i.e. third column of Fig. 2). These graphs clearly reveal that the previous trends with MED designs are reproduced here. The bimodal nature of AOP HIGH in this case is evident for all pupils. The EDOF HIGH design remains close to the ideal case scenario, but that statement does not hold true for pupils greater than 4mm. Again, for pupils >4mm, the EDOF HIGH produces an appreciable infliction point, at about 1.00D target vergence, however it does not present characteristics of a bimodal pattern seen with AOP HIGH. The performance of the AOMF HIGH design seems to be largely pupil dependent and its performance approaches ideal threshold only for specific pupil×vergence pairs. For example, a favorable pattern is observed when considering near vergence with 3mm pupil or far vergence with 5mm pupil. However, the converse pattern emerges when performance is computed at far vergence with 5mm pupil or at near vergence with 3mm pupil. Again, the AOP HIGH design demonstrates least variability in the shape of TFRIQ as a function of change in pupil size (AUR=5.96), closely followed by the EDOF HIGH (AUR=6.96) with the AOMF HIGH (AUR=8.55) design a distant third. The performance of AOMF design does not exhibit bimodal nature at any of the pupils tested.

Fig. 3 showcases the TFRIQ performance values computed for three different inherent aberration models at 4mm pupil, with a well-centered correcting CL. Here, the AOP LOW design not only satisfies the ideal criterion but also demonstrates minimal variance in TFRIQ values obtained with the three inherent aberration models (AUR=1.41). The same observation holds true for the EDOF LOW design (AUR=1.37). Amongst the LOW designs, AOMF LOW (AUR=1.47) demonstrates greatest variance in TFRIQ values between the three inherent aberration models but it is not distinctly different from the performance of AOP or EDOF designs.

The AOP MED holds to the pattern of the AOP LOW, producing a minimal variance in the TFRIQ response within the three aberration models (AUR=2.10), and also meets the ideal threshold requirements for two out of three SA models. However, both these test cases with AOP MED produce a valley in the through focus output. Again, the AOMF MED design (AUR=3.85) produces greatest variance in TFRIQ, as a function of the three aberration models. It just meets the threshold requirement for Model 1, provided an additional offset of 0.25D is applied to the TFRIQ curves. However, this observation does not hold true for results obtained with Models 2 and 3. Unlike the AOP and AOMF designs, the EDOF MED design remains closer to the ideal threshold requirement. This holds true for all three models. Further, it also demonstrates smaller variance across the three aberration models (AUR=2.84) than AOMF design.

The patterns produced by AOP HIGH, AOMF HIGH and EDOF HIGH are very similar to those obtained by their respective MED design counterparts. The TFRIQ curve obtained with the AOP HIGH design is again distinctly bimodal, with the maximum peak-to-trough variation of approximately 1.5 log units. The amplitude of the bimodal curve does not seem to be strongly associated with the three inherent aberration models. The EDOF HIGH design (AUR=4.29) produces least variability between the TFRIQ values obtained with three different inherent aberration levels, closely followed by AOP HIGH (AUR=4.52) with AOMF HIGH (AUR=5.48).

Fig. 4 depicts the TFRIQ values computed at various preset correcting CL centration positions. The performance values were obtained over a 3.5mm pupil diameter, on a schematic eye embedded with a population mean of on-axis spherical aberration (i.e. Model 2 without asymmetric HOA terms). The AOP and EDOF LOW also remains close to the ideal case for most CL decentration levels tested. Conversely, AOMF LOW offers a TFRIQ function that seems to be slightly sensitive to lens decentration compared to AOP and EDOF. The AOP LOW demonstrates performance that is most insensitive to lens centration (AUR=0.35), closely followed by the EDOF LOW design (AUR=0.61). The sensitivity of the performance of AOMF LOW design (AUR=1.51) to lens centration is approximately 3 times that obtained with the AOP and EDOF LOW designs.

Unlike the AOP LOW design, the AOP MED performance is susceptible to lens decentration. A 0.50mm decentration produced a 1.5 log unit change in the RIQ at 1.50D target vergence. On the other hand, the EDOF MED design demonstrated robust performance with lens decentration. The AOMF MED was the most vulnerable of the three designs, a 0.50mm decentration produced a RIQ change of about 2.5 log units at 1.50D target vergence. The variance in TFRIQ as a function of lens decentration levels obtained with AOMF MED (4.79) is approximately 2–3 times more than those obtained with the EDOF MED design AOP MED (2.19) and (AUR=1.35), respectively.

Next, we examine the results of HIGH designs (i.e. column 3 of Fig. 4). Here, the AOMF HIGH lens is the poorest performer with regard to the sensitivity of the TFRIQ curve as a function of lens decentration. A 0.50mm decentration sharply upset the shape of the entire TFRIQ curve, which seems to affect target vergences encompassing 1.00 to 2.00D. Any further decentration would further accentuates the gap between the obtained TFRIQ and the ideal response, pushing the lens to behave like a single vision lens. The AOP HIGH design also appears to be sensitive to lens centration. However, when compared to AOMF HIGH, the effect of decentration on the degradation of the TFRIQ shape with AOP HIGH seems to limit itself to only distance target vergences. The performance of the EDOF HIGH lens seems to be the most robust of the three HIGH designs. The variance in the TFRIQ values as a function of lens centration, obtained with AOMF HIGH (AUR=9.64) are approximately 3 times of that obtained with AOP HIGH (3.69) and EDOF HIGH (3.09).


The impetus for the development of our current theoretical model was to determine the functional capacity of lenses designed with multiple higher order spherical aberrations under varying pupil size, inherent ocular aberrations and lens decentration.27,28,44,55,65 Performance of the EDOF lenses was compared with that of commercial lenses designed with monotonic or concentric power distribution across the optic zone diameter.9,29

To accomplish this task, various scenarios were simulated to reflect reality. Undoubtedly, pupil responses to ambient lighting assert undue influence when evaluating presbyopic correction. Computing the optical performance at various pupils, ranging from 3 to 5mm, should encompass most daily living conditions experienced by presbyopes. Knowing that lens alignment with the pupil center or the corneal geometric center is imperfect, we included a number of preset lens decentered positions. Our previous work indicates that eyes having the same refractive prescriptions but different inherent spherical aberration perform differently when fitted with identical designs.44 We not only took these learnings into account for our current theoretical modeling, but also expanded the concept to include asymmetric aberrations, described here.49

Marsack et al. have proposed numerous (31) retinal image quality metrics to quantify the visual performance of a human eye.63 Their work suggests that visual Strehl ratio in the Fourier domain is the best surrogate of visual performance as it is designed to capture both optical and neural components of the visual processing. For this reason, we have chosen to use this proposed metric for the current modeling exercise. Further, to be able to palpate smaller differences between the lens types, the retinal image quality scale has been log transformed as suggested by Ravikumar et al.66

When an optical system, free of HOA, is imposed with defocus, every location in the entrance pupil has an identical power, which would need only one spherical refraction to correct all pupils. However, introducing spherical aberration forces the best focus (i.e. the refractive correction needed) to become pupil dependent. This observation becomes relevant to those designs that employ monotonic power profiles to ameliorate presbyopic symptoms, particularly the AOMF design in the current experiment. In such instances, an over refraction may become necessary to gauge the ‘true’ through-focus performance. To facilitate a fair comparison of the through-focus performance measures obtained with AOMF designs with the rest of the study lenses, a deliberate offset was made to ‘mimic’ an over refraction of +0.25D for the LOW and MED designs.

Summarizing the results for LOW designs, the performance of the AOP LOW and EDOF LOW stand out for the following reasons: (a) close to the ideal threshold; (b) relatively insensitive to variation in pupil; (c) relatively invariant across different inherent aberrations and (d) performance relatively invariant to lens centration. However, the actual RIQ values obtained by the EDOF LOW designs were marginally higher than AOP LOW. The AOMF LOW design met the threshold requirements, however its performance was twice as more sensitive to pupil than AOP and EDOF counterparts. Clearly, AOMF LOW's performance was distance-centric (as it had highest RIQ values at the distance vergences amongst all the lenses), which made it more akin to single vision than a MFCL. The increased sensitivity of AOMF LOW's TFRIQ to inherent aberration levels and lens centration, further separates AOMF design as a presbyopic lens.

If one were to speculate and extrapolate these optical performance attributes to real world, the AOMF LOW would offer excellent distance performance, like any other single vision, but may not offer much help at near, especially for ‘genuine’ presbyopes, unless they are prescribed with an over-plus refraction. The pupil dependent changes observed suggest inadequate aid at near vergences exacerbated under lower illumination conditions. Further, the small variability in performance with different inherent aberration profiles may be polarizing when applied to a population. In other words, participants with favorable inherent aberration profiles (i.e. low inherent spherical aberration) would benefit from the lenses while those with substantially positive inherent spherical aberration may experience an effect similar to wearing a single vision lens.

On the other hand, presbyopes wearing the AOP LOW or EDOF LOW designs can expect to achieve improvements at near, even in the low-illumination conditions. Additionally, the minimal variance observed in the performance with different inherent aberrations profiles indicates a more predictable outcome across an unselected population. Whether the improvements at near expected with the AOP LOW or EDOF LOW design come at the expense of significant loss (subjective perception) at any other distance needs to be explored in a clinical trial.

Next, we summarize the results for MED and HIGH designs. Unlike the stand out performance observed with the AOP LOW design, the modeled through focus performance responses with the AOP MED and HIGH designs failed to satisfy the ideal threshold response at most pupil sizes. One may argue that they have partially satisfied the ideal curve requirement, at least when gauged in terms of the range of vergences spanned by the TFRIQ curve (where RIQ>−1.50 log units). However, the biggest caveat in their performance is adherence to a ‘bimodal’ trait, more obvious in HIGH than MED. When transformed to a real world setting, this would result in improved performance at near at an expense of performance at intermediate. Further, the performance of the AOP MED and HIGH designs were found to be sensitive to lens decentration. Despite this potential shortcoming, the performance stands relatively invariant to pupil changes and inherent aberrations; whether these positive attributes put the ‘bimodal’ caveat under a shadow would need further clinical investigation.

The AOMF MED and HIGH designs portray attributes that are similar to AOMF LOW design. When compared with corresponding counterparts of AOP and EDOF, performance of AOMF was not only the most sensitive to pupil, but also to inherent aberrations and lens centration. In contrast, the EDOF MED and HIGH designs not only demonstrated reduced sensitivity to pupil and inherent aberrations, but also remained relatively less sensitive to lens centration, when compared to AOP and AOMF. The EDOF MED and HIGH designs distinguish themselves from their AOP counterparts, by not exhibiting bimodality, thus driving them closer to the ideal response. Again, whether the theoretical improvements observed at intermediate and near vergences translate to participant subjective preference of the EDOF MED and HIGH over AOP and AOMF counterparts needs clinical investigation.

The concept of all refractive, simultaneous-image, presbyopic lenses is underpinned by deliberate introduction of multiple powers over a preferred pupil diameter. This is done in an attempt to offer depth-of-focus to the corrected eye. However, it is of paramount importance to note that the magnitudes and signs of the induced aberrations play a critical role in executing the fine balancing act of maintaining retinal image quality over a critical threshold, across the preferred range of target vergences.44 For a deeper understanding on this subject, let us deliberate two simple illustrations. For a specific pupil diameter, if the higher order aberrations of the correcting lens were similar in magnitude to the inherent ocular aberrations (but opposite in sign) then these aberrations would compensate each other, resulting in a low aberrated optical system. This scenario would unsurprisingly precipitate a much better performance for one specific distance but at the expense of poor performance at other distances.44,67 This can be appreciated when AOMF designs were coupled with schematic eyes embedded with aberrations defined by Model 3. Both AOP and EDOF, which do not rely on a monotonic power distribution, result in a performance that is relatively insensitive to inherent aberrations compared to AOMF.

On the other hand, when the accumulated total aberrations (CL+eye system) go beyond a critical threshold, the retinal image quality at one or more distances over the preferred vergences will be diminished to such an extent that its relative measure with respect to any other distances may become unappreciable. This effect can be appreciated by examining the performance obtained with the AOP HIGH design. As the design was predominantly loaded with power compensation for two vergences (distance and near), the accrual of their total aberrations creates an appreciable massive ‘trough’ in the through focus performance, corresponding to an intermediate vergence.27 This would clinically translate to improved performance at near, however with a significant compromise at intermediate. Also, such high levels of aberrations, concentrated in annular ‘chunks’ of the optic zone, would also result in performance that can become exceedingly sensitive to lens centration. Conversely, concentric designs offer performance that is relatively insensitive to pupil compared with monotonic power profile distributions.

Both of the outcomes (performance akin to single vision and performance which is bimodal), would not alleviate the ‘genuine’ presbyopic symptoms, where visual needs span a continuum of distances from far to near. The current theoretical modeling reveals that approaches which either: utilize a monotonic power profile variation (i.e. predominantly involving primary spherical aberration); or employing a concentric-ring type multifocal (where the distribution of power is allocated in distinct zones) suffer undesirable compromises in the through-focus performance resulting from differences in pupil size, inherent aberrations and lens centration. This modeling demonstrates that utilizing the approach that involves deliberately induced multiple, higher-order spherical aberration terms, resulting in a non-monotonic, aperiodic, refractive power profile design may render the performance of such designs less susceptible to variations in pupil, aberrations and centration. However, these optical performance predictions need to be validated by clinical outcomes, which we intend to do in part 2 of this manuscript.


With the computed through focus retinal image quality as the gauge of optical performance, we have demonstrated that lenses designed with deliberately induced multiple, higher-order, spherical aberration terms theoretically render visual performance less susceptible to variations in pupil, inherent ocular aberrations and lens decentration. Especially, when compared to commercial alternatives that are designed based on monotonic power changes or concentric power distribution across the optic zone radius. Whether these incremental improvements to the through focus retinal image quality and its relative insensitivity toward change in pupil, inherent aberrations and decentration actually translate to a clinically palpable outcome needs thorough investigation via human clinical trials.

Conflicts of interest

This project is entirely funded by the Brien Holden Vision Institute and the Institute has proprietary interests in the intellectual property governing extended depth-of-focus contact lenses (Patent application number WO2014059465 A1) and receives royalties from the sale of center-near contact lenses described in this manuscript.


The authors would like to acknowledge Ms. Perera and Dr. Kim for help with actual lens measurements and Ms. J. Sha for help with post-processing the data. Author extends thanks to Mr. D. Tilia, Dr. P. Erickson and Dr. J. Flanagan for reviewing the manuscript.

B.A. Holden,T.R. Fricke,S.M. Ho
Global vision impairment due to uncorrected presbyopia
Arch Ophthalmol, 126 (2008), pp. 1731-1739
D.H. Akerman
40 is the New 20/20—Presbyopia Equals Opportunity
Contact Lens Spectrum, (2010)
N. Agahi,M.G. Parker
Are today's older people more active than their predecessors? Participation in leisure-time activities in Sweden in 1992 and 2002
Ageing Soc, 25 (2005), pp. 925-941
J. Woods,C. Woods,D. Fonn
Visual performance of a multifocal contact lens versus monovision in established presbyopes
Optom Vis Sci, 92 (2014), pp. 175-182
N. Efron,J.J. Nichols,C.A. Woods,P.B. Morgan
Trends in US contact lens prescribing 2002 to 2014
Optom Vis Sci, 92 (2015), pp. 758-767
P.B. Morgan,C.A. Woods,I.G. Tranoudis
International contact lens prescribing in 2014
Contact Lens Spectr, 30 (2015), pp. 28-33
P.B. Morgan,C.A. Woods,I.G. Tranoudis
International contact lens prescribing in 2015
Contact Lens Spectr, 31 (2016), pp. 24-29
R. Montes-Mico,D. Madrid-Costa,A. Dominguez-Vicent,L. Belda-Salmeron,T. Ferrer-Blasco
In vitro power profiles of multifocal simultaneous vision contact lenses
Contact Lens Anterior Eye, 37 (2014), pp. 162-167
E. Kim,R.C. Bakaraju,K. Ehrmann
Power profiles of commercial multifocal soft contact lenses
Optom Vis Sci, (2016),
K. Fisher
Presbyopic visual performance with modified monovision using multifocal soft contact lenses
Int Contact Lens Clin, 24 (1997), pp. 91-100
C.M. Hutnik,D. O’Hagan
Multifocal contact lenses—look again!
Can J Ophthalmol, 32 (1997), pp. 201-205
M. Guillon,C. Maissa,P. Cooper,K. Girard-Claudon,T.R. Poling
Visual performance of a multi-zone bifocal and a progressive multifocal contact lens
CLAO J, 28 (2002), pp. 88-93
D.F. Ezekiel
Soft Multifocal Contact Lens: Google Patents
A.S. Rajagopalan,E.S. Bennett,V. Lakshimanarayanan
Visual performance of subjects wearing presbyopic contact lenses
K. Richdale,G.L. Mitchell,K. Zadnik
Comparison of multifocal and monovision soft contact lens corrections in patients with low-astigmatic presbyopia
B.J. Evans
Monovision: a review
Ophthalmic Physiol Opt, 27 (2007), pp. 417-439
E.B. Papas,T. Decenzo-Verbeten,D. Fonn
Utility of short-term evaluation of presbyopic contact lens performance
Eye Contact Lens, 35 (2009), pp. 144-148
A. Llorente-Guillemot,S. Garcia-Lazaro,T. Ferrer-Blasco,R.J. Perez-Cambrodi,A. Cervino
Visual performance with simultaneous vision multifocal contact lenses
P.R.B. Fernandes,H.I.F. Neves,D.P. Lopes-Ferreira,J.M.M. Jorge,J.M. Gonzalez-Meijome
Adaptation to multifocal and monovision contact lens correction
Optom Vis Sci, 90 (2013), pp. 228-235
D. Madrid-Costa,S. Garcia-Lazaro,C. Albarran-Diego,T. Ferrer-Blasco,R. Montes-Mico
Visual performance of two simultaneous vision multifocal contact lenses
Ophthalmic Physiol Opt, 33 (2013), pp. 51-56
S. Garcia-Lazaro,T. Ferrer-Blasco,D. Madrid-Costa,C. Albarran-Diego,R. Montes-Mico
Visual performance of four simultaneous-image multifocal contact lenses under dim and glare conditions
Eye Contact Lens, 41 (2015), pp. 19-24
D. Madrid-Costa,J. Ruiz-Alcocer,S. Garcia-Lazaro,T. Ferrer-Blasco,R. Montes-Mico
Optical power distribution of refractive and aspheric multifocal contact lenses: effect of pupil size
Contact Lens Anterior Eye, 38 (2015), pp. 317-321
J. Sha,R.C. Bakaraju,D. Tilia
Short-term visual performance of soft multifocal contact lenses for presbyopia
Arq Bras Oftalmol, 79 (2016), pp. 73-77
Bakaraju RC, Ehrmann K, Ho A, Holden BA. Lenses, Devices, Methods and Systems for Refractive Error. USA: Brien Holden Vision Institute.
D. Tilia,R.C. Bakaraju,J. Chung
Short-term visual performance of novel extended depth-of-focus contact lenses
Optom Vis Sci, 93 (2016), pp. 435-444
D. Tilia,A. Munro,J. Chung
Short-term comparison between extended depth-of-focus prototype contact lenses and a commercially-available center-near multifocal
R.C. Bakaraju
Optical Performance of Simultaneous Vision Multifocal Contact Lenses Using Schematic and Physical Eye Models
University of New South Wales, (2010)
R.C. Bakaraju,K. Ehrmann,D. Falk,A. Ho,E. Papas
Optical performance of multifocal soft contact lenses via a single-pass method
Optom Vis Sci, 89 (2012), pp. 1107-1118
E. Kim,R.C. Bakaraju,K. Ehrmann
Reliability of power profiles measured on NIMO TR1504 (Lambda-X) and effects of lens decentration for single vision, bifocal and multifocal contact lenses
J Optom, (2015),
I. Escudero-Sanz,R. Navarro
Off-axis aberrations of a wide-angle schematic eye model
J Opt Soc Am, 16 (1999), pp. 1881-1891
D.A. Atchison
Optical models for human myopic eyes
R.C. Bakaraju,K. Ehrmann,D. Falk,A. Ho,E. Papas
Physical human model eye and methods of its use to analyse optical performance of soft contact lenses
Opt Express, 18 (2010), pp. 16868-16882
R.C. Bakaraju,K. Ehrmann,E. Papas,A. Ho
Finite schematic eye models and their accuracy to in-vivo data
M. Dubbelman,V. Sicam,G. van der Heijde
The shape of the anterior and posterior surface of the aging human cornea
M. Dubbelman,G.L. van der Heijde
The shape of the aging human lens: curvature, equivalent refractive index and the lens paradox
Vis Res, 41 (2001), pp. 1867-1877
M. Dubbelman,G.L. van der Heijde,H.A. Weeber
Change in shape of the aging human crystalline lens with accommodation
R. Navarro
The optical design of the human eye: a critical review
J Optom, 2 (2009), pp. 3-18
R. Navarro,J. Santamaria,J. Bescos
Accommodation-dependent model of the human eye with aspherics
J Opt Soc Am, 2 (1985), pp. 1273-1281
R.A. Schachar
Effect of change in central lens thickness and lens shape on age-related decline in accommodation
J Cataract Refract Surg, 32 (2006), pp. 1897-1898
M. Sheridan,W.A. Douthwaite
Corneal asphericity and refractive error
Ophthalmic Physiol Opt, 9 (1989), pp. 235-238
D. Siedlecki,H. Kasprzak,B. Pierscionek
Schematic eye with a gradient-index lens and aspheric surfaces
Opt Lett, 29 (2004), pp. 1197-1199
G. Smith
The optical properties of the crystalline lens and their significance
Clin Exp Ophthalmol, 86 (2003), pp. 3-18
G. Smith,P. Bedggood,R. Ashman,M. Daaboul,A. Metha
Exploring ocular aberrations with a schematic human eye model
Optom Vis Sci, 85 (2008), pp. 330-340
R.C. Bakaraju,K. Ehrmann,A. Ho,E. Papas
Inherent ocular spherical aberration and multifocal contact lens optical performance
Optom Vis Sci, 87 (2010), pp. 1009-1022
L. Janoff,O.H. Dabezies Jr.
Power change induced by soft contact lens flexure
CLAO J, 9 (1983), pp. 32-38
B.A. Weissman,K.M. Gardner
Power and radius changes induced in soft contact lens systems by flexure
Am J Optom Physiol Opt, 61 (1984), pp. 239-245
R.A. Applegate,V. Lakshminarayanan
Parametric representation of Stiles–Crawford functions: normal variation of peak location and directionality
J Opt Soc Am A, 10 (1993), pp. 1611-1623
D.A. Atchison,A. Joblin,G. Smith
Influence of Stiles–Crawford effect apodization on spatial visual performance
J Opt Soc Am A Opt Image Sci Vis, 15 (1998), pp. 2545-2551
T.O. Salmon,C. van de Pol
Normal-eye Zernike coefficients and root-mean-square wavefront errors
J Cataract Refract Surg, 32 (2006), pp. 2064-2074
S. Amano,Y. Amano,S. Yamagami
Age-related changes in corneal and ocular higher-order wavefront aberrations
Am J Ophthalmol, 137 (2004), pp. 988-992
A. Mathur,D.A. Atchison,W.N. Charman
Effect of accommodation on peripheral ocular aberrations
A. Mathur,D.A. Atchison,W.N. Charman
Myopia and peripheral ocular aberrations
A. Mathur,D.A. Atchison,W.N. Charman
Effects of age on peripheral ocular aberrations
Opt Express, 18 (2010), pp. 5840-5853
A. Mathur,D.A. Atchison,D.H. Scott
Ocular aberrations in the peripheral visual field
Opt Lett, 33 (2008), pp. 863-865
S. Plainis,G. Ntzilepis,D.A. Atchison,W.N. Charman
Through-focus performance with multifocal contact lenses: effect of binocularity, pupil diameter and inherent ocular aberrations
Ophthalmic Physiol Opt, 33 (2013), pp. 42-50
H. Radhakrishnan,W.N. Charman
Age-related changes in ocular aberrations with accommodation
L.N. Thibos
Retinal image quality for virtual eyes generated by a statistical model of ocular wavefront aberrations
Ophthalmic Physiol Opt, 29 (2009), pp. 288-291
X.H. Wan,S.M. Li,Y. Xiong
Ocular monochromatic aberrations in a rural Chinese adult population
Optom Vis Sci, 91 (2014), pp. 68-75
L. Wang,D.D. Koch
Ocular higher-order aberrations in individuals screened for refractive surgery
J Cataract Refract Surg, 29 (2003), pp. 1896-1903
L. Wang,D.D. Koch
Age-related changes in corneal and ocular higher-order aberrations
Am J Ophthalmol, 138 (2004), pp. 897
author reply 897
R.H. Wei,L. Lim,W.K. Chan,D.T. Tan
Higher order ocular aberrations in eyes with myopia in a Chinese population
J Refract Surg, 22 (2006), pp. 695-702
X. Cheng,A. Bradley,L.N. Thibos
Predicting subjective judgment of best focus with objective image quality metrics
J Vis, 4 (2004), pp. 310-321
J.D. Marsack,L.N. Thibos,R.A. Applegate
Metrics of optical quality derived from wave aberrations predict visual performance
J Vis, 4 (2004), pp. 322-328
L.N. Thibos,X. Hong,A. Bradley,R.A. Applegate
Accuracy and precision of objective refraction from wavefront aberrations
J Vis, 4 (2004), pp. 329-351
A. Plakitsi,W.N. Charman
Comparison of the depths of focus with the naked eye and with three types of presbyopic contact lens correction
J Br Contact Lens Assoc, 18 (1995), pp. 119-125
A. Ravikumar,R.A. Applegate,Y. Shi,H.E. Bedell
Six just-noticeable differences in retinal image quality in 1 line of visual acuity: toward quantification of happy versus unhappy patients with 20/20 acuity
J Cataract Refract Surg, 37 (2011), pp. 1523-1529
Y. Benard,N. Lopez-Gil,R. Legras
Optimizing the subjective depth-of-focus with combinations of fourth- and sixth-order spherical aberration
Corresponding author at: Level 5, Rupert Myers Building, North Wing Gate 14, Barker Street, University of New South Wales, Sydney, NSW 2052, Australia. (Ravi C. Bakaraju
Copyright © 2017. Spanish General Council of Optometry
J Optom 2018;11:10-20 - Vol. 11 Num.1 DOI: 10.1016/j.optom.2017.04.003