Abstract

Abstract

**Abstract**:

**Abstract**
The size of the pupil has a large effect on visual function, and pupil size depends mainly on the adapting luminance, modulated by other factors. Over the last century, a number of formulas have been proposed to describe this dependence. Here we review seven published formulas and develop a new unified formula that incorporates the effects of luminance, size of the adapting field, age of the observer, and whether one or both eyes are adapted. We provide interactive demonstrations and software implementations of the unified formula.

Introduction

In the human eye, pupil diameter ranges between approximately 2 and 8 mm. Pupil diameter has a large effect on the optical transfer function of the eye, as illustrated in Figure 1. Pupil size also has a direct effect on depth of field, as well as on retinal illuminance, which in turn influences contrast sensitivity. For these and other reasons, it is often desirable to estimate the pupil size for a given set of conditions when actual measurements are unavailable. Based on available data, various formulas have been proposed to predict pupil diameter. However, none of the existing formulas incorporate the combined effects of the observer's age, the size of the adapting field, and monocular versus binocular stimulation, all of which are known to have a significant impact on pupil size. The purpose of this note is to review these formulas, convert them all to a common format and common units, and propose a unified formula that includes the effects of luminance, age, monocular adaptation, and field size.

Figure 1

Figure 1

We review the proposed formulas in chronological order, but reserve the studies of age and monocular viewing for the end. We note that many investigators have reported large individual differences (Spring & Stiles, 1948), which may in part account for the discrepancies among formulas. Another basis for disagreements is that in some earlier studies, the roles of age, field size, and monocular adaptation appear not to have been understood and may not have been reported.

Studies generally report the size of the eye's entrance pupil, i.e., the virtual image of the physical pupil as seen through the cornea. We follow that practice here. Retinal illuminance is proportional to the area of the entrance pupil (Atchison & Smith, 2000). We also follow the standard practice of specifying stimulus intensities by photopic luminance, even at light levels in the scotopic range. The studies we cite used a variety of luminance units; we use cd m

^{−2}. Table A1 in the Appendix 4 shows conversion factors for other units.Formulas

Holladay (1926 )

Holladay (1926) collected data from three observers of unknown ages binocularly viewing a large adapting area (the interior of an adaptation chamber). He summarized the results with a formula “for a probable average healthy young eye,” (Holladay, 1926, pp. 309–310): where

*M*is luminance in millilamberts. Converting to cd m^{−2}, we have This formula is shown in Figure 2. Clearly it fails for high luminances, but Holladay's data did not go beyond about 600 cd m^{−2}.Figure 2

Figure 2

Crawford (1936 )

Crawford (1936) collected data from 10 subjects binocularly viewing a 55° diameter adapting field. No information is given regarding the age of the observers. He summarized his data with the equation where

*B*is luminance in cd/ft^{2}and log indicates the logarithm of base 10. Converting to luminance in cd m^{−2}, we have This formula is pictured in Figure 3. Crawford notes that his formula is as much as 2 mm below the results of Holladay (1926), Reeves (1918), and others and ascribes this to individual differences.Figure 3

Figure 3

Moon and Spencer (1944 )

Moon and Spencer (1944) considered the data of various authors, but settled on Blanchard (1918), Reeves (1918), and Crawford (1936) as the most reliable. They acknowledged the prior formula of Crawford (but misquoted it), and offered a modification which they asserted provided “an approximate average of all the data,” (Crawford, 1936, p. 321): where

*M*is in millilamberts. Converting to cd m^{−2}, we get which we approximate by This formula is shown in Figure 4. It was also adopted without attribution by Le Grand (1968) (who replaces 4.9 with 5), and appears to be given incorrectly in Wyszecki and Stiles (1982) [Equation 1(2.4.5), p. 106].Figure 4

Figure 4

De Groot and Gebhard (1952 )

De Groot and Gebhard (1952) computed a mean of data sets from eight authors, weighting each set by the number of observers. They noted that previous formulas tended to assume a lower asymptote of 2 mm, not necessarily demanded by the data or physiological constraints. They provided a formula which does not asymptote at 2 mm and was a better fit to their mean than the formula of Moon and Spencer (1944), where M is in millilamberts. Converting to cd m

^{−2}, and expressing as a formula for*D*, we have as shown in Figure 5.Figure 5

Figure 5

Stanley and Davies (1995 )

Crawford (1936) and Bouma (1965) demonstrated that pupil size was dependent not on luminance alone, but approximately on the product of luminance and adapting field size. This observation appears as early as Aubert (1876) [quoted in Schweitzer (1956)], but it did not appear in formulas for pupil size until the work of Stanley and Davies (1995), who measured pupil diameter for six observers as a function of luminance and adapting field sizes ranging in diameter from 0.4 to 25.4°. They observed a strong effect of field size, which they suggested might account for the large discrepancies in earlier published results. Evidently unaware of the earlier observations of Crawford and Bouma, they noted that data for all sizes superimposed if the results were plotted as a function of the product of luminance and adapting area in degrees squared (

*corneal flux density*, cd m^{−2}deg^{2}). They fit the combined data with the following function where*a*is the area in deg^{2}. Two examples of the function, for circular fields with diameters of 0.4 and 25.4°, are shown in Figure 6. Those two sizes are the limits of the range tested by Stanley and Davies, who caution against extrapolation outside those bounds. However Bouma (1965) appears to show integration out to much larger diameters, so we have not restricted the diameter in the formula. Atchison et al. (2011) provide a more recent confirmation of the dependence of pupil diameter on corneal flux density.Figure 6

Figure 6

Barten (1999 )

Barten (1999) adopted the formula of Le Grand (1968), itself borrowed from the formula of Moon and Spencer (1944), but inserted a term to modulate luminance by the field area, as in the formula of Stanley and Davies (1995), but here based on the work of Bouma (1965). Barten's formula is given by Barten's formula is plotted in Figure 7, along with that of Stanley and Davies for the same area, to show the similarity.

Figure 7

Figure 7

Blackie and Howland (1999 )

As part of a larger project of modeling the emmetropization of the eye, Blackie and Howland (1999) fit the pupil size data of Flamant (1948) with a formula defined by where

*R*is the log of luminance in cd cm^{−2}. Converting to*L*in cd m^{−2}, and simplifying, we have the formula given below and plotted in Figure 8.Figure 8

Figure 8

Winn, Whitaker, Elliott, and Phillips (1994 )

Winn et al. (1994) measured pupil sizes for 91 subjects ranging in age from 17 to 83 years with normal healthy eyes. Adaptation was monocular, and the adaptation field was a 10° diameter circular disk. Accommodation was controlled and minimized. Five luminance levels between 9 and 4400 cd m

^{−2}were used. For each luminance, pupil diameter was plotted against age, and a linear regression was fit to each set, as shown in Figure 9. The estimated slopes and intercepts are shown in Table 1, along with*r*values.^{2}Figure 9

Figure 9

Table 1

Table 1

cd m^{−2} | Slope | Intercept | r ^{2} |

9 | −0.043 | 8.046 | 0.557 |

44 | −0.04 | 7.413 | 0.486 |

220 | −0.032 | 6.275 | 0.377 |

1100 | −0.02 | 4.854 | 0.226 |

4400 | −0.015 | 4.07 | 0.214 |

The slope values are plotted in Figure 10 as a function of log luminance. In order to create a formula based on the Winn results, we have first fit the slope values as a function of log luminance with a cubic polynomial. This is then used to construct a slope function

*W*(_{S}*L*). It specifies the proportional variation of pupil diameter with age at a given luminance. It is given by In this expression, we restrict the luminance to the range explored by Winn et al., and as a result the function is flat at high and low luminances. This function is plotted as the red curve in Figure 10.Figure 10

Figure 10

We also fit the intercepts as a function of log luminance with a cubic polynomial. We use this to construct a function

*W*that returns an intercept for a given luminance, This function is plotted along with the data in Figure 11._{I}Figure 11

Figure 11

Figure 12

Figure 12

Developing the unified formula

Monocular effect

In pupil experiments it is common to stimulate only one eye and measure pupil size in the other. The effect will generally be different than if both eyes view the same adapting light; a useful pupil formula should predict the difference. Blanchard (1918) and Reeves (1918) provide data comparing monocular and binocular adaptation in a single observer (both present identical figures). ten Doesschate and Alpern (1967) provide similar data for seven observers. Adaptation of one eye produces a larger pupil than adaptation of both eyes. Data from one observer are shown in Figure 13, which illustrates that the effect is as much as 1 mm, but diminishes at the extremes of high and low luminance.

Figure 13

Figure 13

We have found that monocular and binocular curves (like those in Figure 13) can be roughly superimposed by a horizontal shift on a log axis. To illustrate, we plot in Figure 14 the result of shifting the monocular data for all eight available data sets. In each case we optimized the shift by minimizing the area (shown in yellow) between the curves in the range of luminances where they overlap. The average optimal shift was −0.994 log units, or a factor of 0.1015.

Figure 14

Figure 14

In light of the success of the shifting operation, we can model the monocular effect by simply assuming an attenuation of luminance by a factor of 0.1 when monocular viewing is used. This result is new and surprising. It also accounts almost perfectly for the results of Bartleson (1968). We formalize this with a function expressing attenuation as a function

*M*(*e*) of number of eyes*e*,Effective corneal flux density

We have seen that pupil diameter is to a first approximation dependent on the product of luminance and adapting field area, or corneal flux density (Atchison et al., 2011; Aubert, 1876; Bouma, 1965; Crawford, 1936; Stanley & Davies, 1995). In the previous section we noted that corneal flux density is effectively attenuated by a factor of 10 if only one eye is adapted. We introduce the concept of

*effective corneal flux density*to describe the quantity that effectively controls pupil diameter, equal to the product of luminance, area, and the monocular effect,*F*=*LaM*(*e*).Age slope

The data of Winn et al. (1994) show that pupil diameter varies with age, and the slope function

*W*provides a means of adjusting computed values as a function of age. However, the function is deficient in two ways. First, it is limited to a modest range of luminance (9 to 4400 cd m_{S}^{−2}), and second, it is defined only for monocular viewing and a 10° diameter adapting field. We generalize the function in the following ways. First, we assume that the slope is a function of effective corneal flux density, rather than luminance alone. Second, we observe that the reason the slope magnitude declines with age is because the maximum pupil diameter, and thus the total range of diameters, declines with age. If this is so, we should find a simple relation between the slope and the computed pupil diameter at the same luminance. In Figure 15a we show the fit of a linear model (*r*^{2}= 0.986) between the slope estimates from Winn et al. and the pupil diameters computed from the formula of Stanley and Davies for the corresponding effective corneal flux density (using Winn's conditions of monocular viewing and 10° diameter adapting field).Figure 15

Figure 15

Using that linear relation, we can create a new age slope function that is a linear transformation of the pupil diameter as computed by the formula of Stanley and Davies The age slope function (for the conditions of Winn et al.) is plotted in Figure 15b, along with the values estimated by Winn. The agreement is evident. But note that the age slope function extends beyond the limits explored by Winn, as we intended.

Age effect

The age slope function describes the rate of change of pupil diameter with age (mm/year), as a function of effective corneal flux density. We construct an age effect function that describes the change in pupil diameter as a function of luminance, area, age, a reference age, and number of eyes, This function will be used to adjust a reference formula for pupil diameter as a function of age and luminance. It applies only to ages between 20 and 83, the range investigated by Winn et al. (1994). In Appendix 1 we discuss an extension of the formula to ages below 20. The reference age is the age for which the reference function is defined. Ideally it would be the mean age of the group of observers on whose data the reference formula was based.

Unified formula

In addition to the well known effect of luminance, we have seen that there are systematic effects of age (Winn et al., 1994), field size, and number of eyes adapted (Blanchard, 1918; Reeves, 1918; ten Doesschate & Alpern, 1967). Here we construct a unified formula that includes all four effects where

*a*is field area in deg^{2},*y*is age in years,*y*is the reference age, and_{0}*e*is the number of eyes (one or two). Since we are using the formula of Stanley and Davies as the reference formula, the reference age should be the mean age of the population of observers used by Stanley and Davies. If we recall the definition of effective corneal flux density*F*=*LaM*(*e*), we can rewrite the unified formula in this expanded but simplified form that makes clear the role of the Stanley Davies formula (Equation 10) and the age effectEstimating the reference age

As noted above, the reference age

*y*should be the mean age of the population of observers used by Stanley and Davies. This mean age is unknown, but we can estimate it by comparing the unified formula with the data of Winn et al. (1994) as shown in Figure 9 above. We have plotted predictions from our unified formula (red) with the reference age parameter_{0}*y*estimated by minimizing the mean squared error between the predictions and the data. The estimated value is_{0}*y*= 28.58. The agreement between our unified formula (red) and the Winn estimates (blue) is striking and provides additional confirmation of the robustness of our approach._{0}The unified formula is shown in Figure 16 along with the other formulas developed above. We omit the formula of LeGrand, as it is nearly identical to Moon & Spencer (1944). In Figure 16a we have used parameters of field area

*a*= 900*π*(diameter = 60°), age*y*= 30, and eyes*e*= 2. This approximates the curve of Moon and Spencer, even though the unified formula is based on that of Stanley and Davies. Likewise in Figure 16b we show that when the field diameter and number of eyes are set to match those of Winn et al. (10° and monocular viewing), the unified formula provides an excellent fit to their curve over the range they measured.Figure 16

Figure 16

Demonstration and calculator

To assist the reader in visualizing the effects of luminance, age, field size, and binocularity on pupil diameter, and to allow calculation of specific values, we present a demonstration below in Figure 17. The demonstration makes use of the free Wolfram CDF player available at http://www.wolfram.com/cdf-player/. The demonstration allows the reader to set the values of the various parameters. A single value calculator is shown in Figure 18.

Figure 17

Figure 17

Figure 18

Figure 18

Retinal illuminance versus luminance

Equipped with a formula for pupil diameter as a function of luminance, we can easily compute the expected retinal illuminance corresponding to a given luminance. Of course, the relationship will depend somewhat upon age and other parameters, as shown in the demonstration in Figure 19. The near linearity of these log-log curves is striking (

*r*^{2}> 0.99 for all ages). If the relationship was in fact linear, it would imply a linear relationship between log luminance and log pupil diameter. This is not quite so, though the error is usually less than 1 mm.Figure 19

Figure 19

Discussion

Unified formula

The development of our unified formula follows a logical argument that we summarize here.

- At a given age, the pupil diameter is some function of the effective corneal flux density, given by
*F*=*LaM*(*e*). Evidence for this are the studies on field size (Atchison et al., 2011; Aubert, 1876; Bouma, 1965; Crawford, 1936; Stanley & Davies, 1995), the data on monocular adaptation (Bartleson, 1968; Blanchard, 1918; Reeves, 1918; ten Doesschate & Alpern, 1967), and our analysis of those data (Figure 13). - At the reference age, we take that function to be the formula of Stanley and Davies (1995). The rationale is that Stanley and Davies used relatively well-defined and modern methods and six observers, and their formula is simple and consistent with other historical formulas.
- At any other age, pupil diameter is adjusted by a linear function of age relative to the reference age. Supporting this are the data of Winn et al. (1994) showing a linear variation with age at a range of luminances.
- The slope of the linear age adjustment at any luminance is itself a linear function of the reference pupil diameter (pupil diameter at the reference age for that luminance). The rationale is that the slope depends on the available range of pupil diameter variation, and that in turn is dependent upon the reference pupil diameter.

Variability

Our formula describes the mean pupil diameter for a given condition, but it should be understood that there may be considerable variability around that mean. We identify four possible sources of variability: between the two eyes, in a single eye over a short interval of time (seconds), in a single eye over a long interval of time (days), and between different observers.

Typically the two pupils have approximately equal diameters, regardless of whether one or both eyes are illuminated (Loewenfeld & Lowenstein, 1999). Large discrepancies (>0.4 mm) are called anisocoria, and may be diagnostic of neural pathology. Ettinger, Wyatt, and London (1991) found mean absolute differences of between 0.12 mm (at 343 cd m

^{−2}) and 0.36 mm (in darkness).Careful observation shows that pupil size under prolonged constant illumination varies continuously over time. The most significant variation is generally known as

*pupillary unrest*, or sometimes*hippus*. It is a spectrally broadband low frequency random fluctuation with a bandpass spectrum ranging from 0.02 to 2.0 Hz and an amplitude around 0.25 mm, somewhat dependent on the mean (Stark, Campbell, & Atwood, 1958). It is synchronous in the two eyes and thus does not contribute to anisocoria. It never occurs in total darkness, and its magnitude is greatest at the midrange of pupil diameters (Loewenfeld & Lowenstein, 1999). Pupillary unrest complicates the measurement of steady state pupil diameter that our formula is designed to compute and may contribute to some of the variability noted in the literature. Methods using still photography will produce variable results depending on the momentary phase of the waveform; much better are modern methods that continuously monitor pupil diameter over an interval of time and then average out the unrest (Winn et al., 1994).Variability between observers is in part due to variations in age, and that component is dealt with by our unified formula. However, considerable variability remains among observers of the same age. We have estimated this variability from the data of Winn et al. (1994), extracted from their Figure 2 as described above (Figure 9). We estimated the standard deviation of the departures of the points from the best-fitting linear function of age. The result, expressed as a function of adapting luminance, is shown in Figure 20. The values range between about 1 and 0.6 mm, and are smallest at the largest luminances.

Figure 20

Figure 20

Other influences on pupil size

Our formula assumes that the observer fixates the center of the adapting field, consistent with the instructions given in the experiments of Winn et al. (1994). Crawford (1936) compared pupil sizes produced by a glare source located at eccentricities of 0 to 56°, finding an exponential decline of as much as 1.73 mm with eccentricity. We have found no studies other than Crawford's that systematically explore the effect of the placement of a uniform adapting field within the visual field, but his results do suggest the possibility of an effect, and suggest that our formula should be used with caution when the adapting field is not uniform or not centered on fixation.

While we have dwelt on the influence of adapting field luminance, it is important to acknowledge that there are many other factors that influence pupil size. Among these are accommodation (Kasthurirangan & Glasser, 2005; Marg & Morgan, 1949), mental activity (Hess & Polt, 1964), emotional arousal (Bradley, Miccoli, Escrig, & Lang, 2008), contrast (Barbur, 2004), detection (Privitera, Renninger, Carney, Klein, & Aguilar, 2010), recognition (Heaver & Hutton, 2011), and attention (Hoeks & Levelt, 1993). However, the cognitive effects are quite small, and all of these responses are largely transient, meaning the pupil diameter returns near to its prevailing value after a number of seconds.

Area summation and ipRGCs

Recently, intrinsically photosensitive retinal ganglion cells (ipRGC) have been discovered in the primate retina (Dacey et al., 2005). These cells express the photopigment melanopsin, with an action spectrum peaking at about 483 nm. These cells have been shown to exert significant control over the steady state pupil response in both macaque and human (Gamlin et al., 2007; McDougal & Gamlin, 2010). Thus while we have considered photopic luminance as the controlling variable, the action spectrum of the actual controlling quantity is likely to be a more complex mixture of rod, cone, and ipRGC sensitivities. The discovery of ipRGCs may also explain one mystery of the pupillary reflex. Our unified formula exhibits perfect summation of luminance over space (see formula of Stanley and Davies above), yet typical retinal ganglion cells have relatively small summation areas. In contrast, ipRGCs have very large receptive fields (Dacey et al., 2005) more consistent with the pupillary reflex.

Conclusions

Based on existing data and formulas, we have derived a unified formula for light adapted pupil diameter that includes the effects of luminance, adapting field size, age of the observer, and whether one eye or two are used. Our formula only describes the steady state in what is otherwise a dynamic reflex subject to many influences. Likewise there are large individual differences. With these caveats, we hope that our formula will be of some use in scientific and practical applications.

Acknowledgments

This work was supported by NASA Space Human Factors Engineering (WBS 466199). We thank Larry Thibos for providing the wavefront data used to create Figure 1.

Commercial relationships: none.

Corresponding author: Andrew B. Watson.

Email: andrew.b.watson@nasa.gov

Address: NASA Ames Research Center, Moffett Field, CA, USA.

References

Atchison D.
Smith G.
(2000).

*Optics of the Human Eye*. Oxford: Butterworth-Heinemann.
Atchison D. A.
Girgenti C. C.
Campbell G. M.
Dodds J. P.
Byrnes T. M.
Zele A. J.
(2011). Influence of field size on pupil diameter under photopic and mesopic light levels.

*Clinical and Experimental Optometry*, 94(6), 545–548, http://www.ncbi.nlm.nih.gov/pubmed/21929524. [CrossRef] [PubMed]
Aubert H.
(1876). Physiologische Optik. In Graefe
A.
Saemisch T.
(Eds.),

*Handbuch Der Gesamten Augenheilkunde*(Vol. 2, p. 453), Leipzig: Wilhelm Engelmann.
Barbur J. L.
(2004). Learning from the pupil-studies of basic mechanisms and clinical applications. In Chalupa
L. M.
Werner J. S.
(Eds.),

*The Visual Neurosciences*. (Vol. 1, pp. 641–656). Cambridge, MA: MIT Press.
Barten P. G. J.
(1999).

*Contrast sensitivity of the human eye and its effects on image quality*. Bellingham, WA: SPIE Optical Engineering Press.
Blackie C. A.
Howland H. C.
(1999). An extension of an accommodation and convergence model of emmetropization to include the effects of illumination intensity.

*Ophthalmic and Physiological Optics**,*19(2), 112–125. [CrossRef] [PubMed]
Blanchard J.
(1918). The brightness sensibility of the retina.

*The Physical Review**,*11(2), 81–99. [CrossRef]
Bouma H.
(1965). Receptive systems: Mediating certain light reactions of the pupil of the human eye.

*Philips Research Reports**,*(Suppl. 5).
Bradley M. M.
Miccoli L.
Escrig M. A.
Lang P. J.
(2008). The pupil as a measure of emotional arousal and autonomic activation.

*Psychophysiology*, 45(4), 602–607. [CrossRef] [PubMed]
Crawford B. H.
(1936). The dependence of pupil size upon external light stimulus under static and variable conditions.

*Proceedings of the Royal Society of London*,*Series B, Biological Sciences*, 121(823), 376–395. [CrossRef]
Dacey D. M.
Liao H. W.
Peterson B. B.
Robinson F. R.
Smith V. C.
Pokorny J.
(2005). Melanopsin-expressing ganglion cells in primate retina signal colour and irradiance and project to the LGN.

*Nature**,*433(7027), 749–754. [PubMed] [CrossRef] [PubMed]
De Groot S. G.
Gebhard J. W.
(1952). Pupil size as determined by adapting luminance.

*Journal of the Optical Society of America A**,*42(7), 492–495. [PubMed] [CrossRef]
Ettinger E. R.
Wyatt H. J.
London R.
(1991). Anisocoria. Variation and clinical observation with different conditions of illumination and accommodation.

*Investigative Ophthalmology & Visual Science**,*32(3):501–509. http://www.iovs.org/content/32/3/501. [PubMed] [Article] [PubMed]
Flamant F.
(1948). Variation du diametre de la pupille de l'oeil en fonction de la brilliance.

*Revue d'Optique**,*27, 751–758.
Gamlin P. D.
McDougal D. H.
Pokorny J.
Smith V. C.
Yau K. W.
Dacey D. M.
(2007). Human and macaque pupil responses driven by melanopsin-containing retinal ganglion cells.

*Vision Research**,*47(7), 946–954. [PubMed] [CrossRef] [PubMed]
Heaver B.
Hutton S. B.
(2011). Keeping an eye on the truth? Pupil size changes associated with recognition memory.

*Memory**,*19(4), 398–405. [CrossRef] [PubMed]
Hess E. H.
Polt J. M.
(1964). Pupil size in relation to mental sctivity during simple problem-solving.

*Science**,*143(3611), 1190–1192. [CrossRef] [PubMed]
Hoeks B.
Levelt W.
(1993). Pupillary dilation as a measure of attention: A quantitative system analysis.

*Behavior Research Methods**,*25(1), 16–26. [CrossRef]
Holladay L.
(1926). The fundamentals of glare and visibility.

*Journal of the Optical Society of America*, 12(4), 271–319. [CrossRef]
Kasthurirangan S.
Glasser A.
(2005). Characteristics of pupil responses during far-to-near and near-to-far accommodation.

*Ophthalmic and Physiological Optics**,*25(4), 328–339. [PubMed] [CrossRef] [PubMed]
Kobashi H.
Kamiya K.
Ishikawa H.
Goseki T.
Shimizu K.
(2012). Daytime variations in pupil size under photopic conditions.

*Optometry & Vision Science**,*89(2), 197. [CrossRef]
Le Grand Y.
(1968).

*Light, colour and vision*(English 2nd ed.). Translated by R. W. G. Hunt, J. W. T. Walsh, & F. R. W. Hunt. London: Chapman & Hall.
Loewenfeld I. E.
Lowenstein O.
(1999).

*The pupil: Anatomy, physiology, and clinical applications*. Boston: Butterworth-Heinemann.
MacLachlan C.
Howland H. C.
(2002). Normal values and standard deviations for pupil diameter and interpupillary distance in subjects aged 1 month to 19 years.

*Ophthalmic and Physiological Optics**,*22(3), 175–182. [CrossRef] [PubMed]
Marg E.
Morgan
M. W.Jr.
(1949). The pupillary near reflex; The relation of pupillary diameter to accommodation and the various components of convergence.

*American Journal of Optometry & Archives of American Academy of Optometry*, 26(5), 183–198. [PubMed] [CrossRef]
McDougal D. H.
Gamlin P. D.
(2010). The influence of intrinsically-photosensitive retinal ganglion cells on the spectral sensitivity and response dynamics of the human pupillary light reflex.

*Vision Research**,*50(1), 72–87. [PubMed] [CrossRef] [PubMed]
Moon P.
Spencer D. E.
(1944). On the Stiles-Crawford effect.

*Journal of the Optical Society of America*, 34(6), 319–329, http://www.opticsinfobase.org/abstract.cfm?URI=josa-34-6-319. [CrossRef]
Privitera C. M.
Renninger L. W.
Carney T.
Klein S.
Aguilar M.
(2010). Pupil dilation during visual target detection.

*Journal of Vision*, 10(10):3, 1–14, http://www.journalofvision.org/content/10/10/3, doi:10.1167/10.10.3. [PubMed] [Article] [CrossRef] [PubMed]
Ravikumar S.
Thibos L. N.
Bradley A.
(2008). Calculation of retinal image quality for polychromatic light.

*Journal of the Optical Society of America, A, Optics, Image Science, and Vision**,*25(10), 2395–2407. [PubMed] [CrossRef]
Reeves P.
(1918). Rate of pupillary dilation and contraction.

*Psychological Review*, 25(4), 330–340. [CrossRef]
Robl C.
Sliesoraityte I.
Hillenkamp J.
Prahs P.
Lohmann C. P.
Helbig H.
(2009). Repeated pupil size measurements in refractive surgery candidates.

*Journal of Cataract and Refractive Surgery**,*35(12), 2099–2102. [PubMed] [CrossRef] [PubMed]
Schweitzer N. M. J.
(1956). Threshold measurements on the light reflex of the pupil in the dark adapted eye.

*Documenta Ophthalmologica**,*10(1), 1–78.
Spring K. H.
Stiles W. S.
(1948). Variation of pupil size with change in the angle at which the light stimulus strikes the retina.

*British Journal of Ophthalmology**,*32(6), 340–346. [PubMed] [CrossRef] [PubMed]
Stanley P. A.
Davies A. K.
(1995). The effect of field of view size on steady-state pupil diameter.

*Ophthalmic & Physiological Optics**,*15(6), 601–603. [PubMed] [CrossRef]
Stark L.
Campbell F. W.
Atwood J.
(1958). Pupil unrest: An example of noise in a biological servomechanism.

*Nature**,*182(4639), 857–858. [CrossRef] [PubMed]
Thibos L. N.
Hong X.
Bradley A.
Cheng X.
(2002). Statistical variation of aberration structure and image quality in a normal population of healthy eyes.

*Journal of the Optical Society of America A**,*19(12), 2329–2348. [CrossRef]
Watson A. B.
Ahumada
A. J.Jr.
(2008). Predicting visual acuity from wavefront aberrations.

*Journal of Vision*, 8(4):17, 1–19, http://www.journalofvision.org/8/4/17, doi:10.1167/8.4.17. [PubMed] [Article] [CrossRef] [PubMed]
Watson A. B.
Ahumada
A. J.Jr.
(2012, in press). Modeling acuity for optotypes varying in complexity.

*Journal of Vision*.
Winn B.
Whitaker D.
Elliott D. B.
Phillips N. J.
(1994). Factors affecting light-adapted pupil size in normal human subjects.

*Investigative Ophthalmology & Visual Science*, 35(3):1132–1137, http://www.iovs.org/content/35/3/1132. [PubMed] [Article] [PubMed]
Wyszecki G.
Stiles W. S.
(1982).

*Color science*(2nd ed.). New York: John Wiley and Sons.Appendix 1: Extension to ages below 20 years

The unified formula developed in the body of this paper computes pupil diameter for ages between 20 and 80 years, the limits of Winn et al.'s (1994) data. Below we provide an extension to ages between 1 and 20. This extension rests on a small amount of data and some assumptions, which is why we relegate it to an appendix. But we include it in part to inspire further research to confirm, refute, or refine the formula.

MacLachlan and Howland (2002 )

MacLachlan and Howland (2002) measured pupil size in younger observers aged 1 month to 19 years. Pupils were illuminated by dim ambient illumination of 15.9 lux. Observers were grouped by gender and by age into bins about one year in width. In their Table 1 they published for each group the mean age, mean pupil diameter, standard deviation, and number of observers. They provide functions describing pupil diameter as function of age separately for male and females. However they did not describe the method by which they arrived at the functions, and we sought a single function, so we have fit our own function to the data. The differences between male and female were small and inconsistent so we ignored gender. The standard deviations varied little between groups and averaged 0.93 mm. The number of subjects per group declined strongly with age, ranging from 378 for females of mean age 0.59 to 5 at mean age 18.82; this trend was similar for males. We fit the ensemble group data (mean pupil diameter versus mean age) using a linear model, with weights given by the number of subjects in the group divided by the group variance.

We considered two linear models, a polynomial of age or of log age, and in each case we considered polynomial of various orders. We evaluated the fit of the model using adjusted

*r*^{2}, AIC, and BIC statistics. By all three criteria the best fit was given by a third order polynomial of log age. The resulting fit is shown in Figure A1.We use this fit to define a function describing pupil diameter as a function of age for the conditions of MacLachlan and Howland (2002)

Equivalent age

The age slope function derived above will allow us to adjust pupil diameter for age, but only for ages greater than 20, the lower limit of Winn's data. We know that below 20, the trend with age reverses and pupils become smaller, as shown in Figure 12 (MacLachlan & Howland, 2002). Because data on younger observers were collected at only one low effective corneal flux density (

*F*), we do not know how young pupils react to adapting luminance. We make the assumption that younger pupils will behave in an equivalent fashion to those of an observer at an age which yields the same pupil size at_{MH}*F*._{MH}First we estimate

*F*_{MH}. We note that the pupil diameter for a 20 year old from MacLachlan and Howland is about equal to the unified result for a 20 year old with a 20° diameter binocular adapting field of 0.174 cd m^{−2}and so we estimate*F*as 0.174 cd m_{MH}^{−2}× 100*π*deg^{2}= 54.66. Next we compute, for ages*y*from 1 to 19, the age required by*D*at_{U}*F*to yield the same pupil diameter as_{MH}*D*(_{MH}*y*). We discovered that those values were fit exceptionally well by an exponential (constrained to have a value of 20 at 20). The result of the fit is shown in Figure A2.Unified formula for ages below 20

The equivalent age transformation allows us to extend our unified formula to ages below 20 by mapping the younger age to an older age, and then using the standard formula, An illustration of the effect of age from 1 to 80 years is shown in Figure A3. The figure shows the change in pupil diameter from a peak at 20 years for two different adapting luminances. Pupil size declines on either side of the peak at age 20, but does so much more modestly for bright adapting fields.

We include this extension to younger ages in our demonstration and calculator (above) and our simplified formulas (below) but it should be understood that the extension is based on much less data, and many more assumptions than the formula for ages 20 to 80 years.

Figure A1

Figure A1

Figure A2

Figure A2

Figure A3

Figure A3

Appendix 2: Simplified formulas

Our unified formula is defined above through a series of nested expressions, notably Equations 10 and 17–24. However for simplicity of calculation it is possible to produce simplified final expressions for the two cases of age less than and greater than or equal to 20. First we define a term that is the effective corneal flux density, raised to the power 0.41, Then

Appendix 3: Mathematica notebook

An implementation of the unified formula is provided in a Mathematica Notebook called PupilDiameter.nb. It contains functions implementing all of the individual formulas described above, as well as the unified formula. A few plotting functions are also provided to illustrate the formulas. Mathematica is computer system for programming and mathematics (http://www.wolfram.com).

Appendix 4: Units