Convergence of experimental multiple Rayleigh matches to. peak L- and M- photopigment sensitivity estimates

Size: px
Start display at page:

Download "Convergence of experimental multiple Rayleigh matches to. peak L- and M- photopigment sensitivity estimates"

Transcription

1 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 1 VIÉNOT F, SERREAULT L., ARDO FERNANDEZ., 2006, Convergence of experimental multiple Rayleigh matches to peak L- an M- photopigment sensitivity estimates, Visual Neuroscience, 23: Convergence of experimental multiple Rayleigh matches to peak L- an M- photopigment sensitivity estimates FRANÇOISE VIÉNOT 1, LUCILE SERREAULT 2 an EDRO ARDO FERNANDEZ 3 1 Muséum National Histoire Naturelle, CRCDG, aris, France 2 Université aris-su, MST Optique physiologique, Optique e contact et Optométrie, Orsay, France 3 Universia e Extremaura, Departamento e Física, Baajoz, Spain Corresponing author Françoise Viénot Muséum National Histoire Naturelle Centre e Recherches et e Restauration es Documents Graphiques 36 rue Geoffroy Saint-Hilaire F aris France vienot@mnhn.fr Tel: Fax: Abstract Shift in the wavelength of peak sensitivity of the cone photopigments is a major cause of inter-iniviual variations in the Rayleigh match. Normal colour observers performe

2 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 2 multiple Rayleigh matches, i.e. a series of seven Rayleigh-type colour matches using various sets of test an primary lights, in orer to erive estimates of the M- an L- photopigment wavelengths of maximum absorbance. We preicte the ratio of incient raiances r / g for a range of plausible values in terms of wavelength of the peak sensitivity of the M-cone an L- cone. An algorithm allowe us to ajust the values of the peak wavelength of the cone photopigments in orer to minimize the sum of square ifferences between preicte an real results of matches. To create caniate M-cone an L-cone photopigments that best preict the set of equation values, we use the low ensity spectral absorbance curves of the M-cone an L-cone photopigments respectively, as tabulate by Stockman an Sharpe 2000, an shifte each template along a logarithmic wavelength scale. For all observers, the system of 7 Rayleigh-like equations converges to a unique solution. M- an L- cone photopigments are estimate separately. The istribution of the wavelength of maximum sensitivity of the photopigments inclues a cluster. Only a few women s results lie outsie the cluster. The choice of the template has a consierable influence on the convergence of the algorithm. Keywors Multiple Rayleigh colour matches, Cone photopigments, Spectral absorbance template, Colorimetry, Human colour vision Introuction Normal colour vision is trichromatic. Colour matches are explaine by the absorption of photons in three families of photopigments: the S-cone photopigment, the M-cone photopigment, the L-cone photopigment. During the past 20 years, a polymorphism of the M- cone an L-cone photopigments has been shown in humans, by molecular genetics an psychophysical stuies see a review in Sharpe et al., 1999, table 1.3. Its consequence on

3 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 3 metameric colour-matches Neitz & Jacobs, 1990 an the spectral luminosity function Carroll et al., 2000 has been moelle an verifie. Shift in the wavelength of peak sensitivity of the cone photopigments is recognise as a major cause of inter-iniviual variations in the Rayleigh matches. Neitz & Jacobs 1990, comparing two ifferent colour-matches from colour normal male observers an consiering together finings from molecular biology, suggeste that there are iscrete variations in both mile an long wavelength cone pigments. He an Shevell 1994, measuring the ratio of two particular Rayleigh-type matches the Dual-Rayleigh, investigate iniviual ifferences in cone photopigments of normal trichromats. Specifically, they showe that the range of the results for 17 observers was too large to be explaine by only iniviual ifferences in photopigment optical ensity an pre-receptoral filtering. Their results were accounte for quantitatively by a small ifference 3-5 nm in the maximum of the L-cone photopigment. He an Shevell left one question unsolve: Which photopigments varies in lambamax? Rayleigh matches are constraine by the spectral positioning both of the M-cone an L- cone photopigments. In this stuy, we aress the inverse problem. We explore a screening metho that specifically estimates the M-cone an L-cone photopigment wavelengths of maximum absorbance in the absence of known genotype. We propose to ask observers to make multiple colour matches, i.e. a series of Rayleigh-type colour matches using various sets of test an primaries, in orer to reveal variations in M-cone an L-cone photopigments. We investigate the convergence an iscuss the accuracy of the metho. Methos Rayleigh- type colour matches The Rayleigh match is a match between a near-monochromatic yellow light an an aitive mixture of near-monochromatic green an re lights. The classical proceure makes use of

4 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb nm, 670 nm an 545 nm raiations. The test an the mixture are presente to the observer on a bipartite fiel. A match is achieve when the same number of photons is absorbe on the left an right half-fiels, in each family of cones. For Rayleigh-like matches where the S- cones are insensitive to the primary an test lights, this necessary conition can be expresse using two equations. We use the term cone funamental to mean the spectral sensitivity of the cones measure at the cornea. Any inert pigment inclue in the eye acts as a filter that transmits a fraction τ of the light at wavelength. It reuces the sensitivity of the eye compare to the sensitivity of the cones. Thus the number of photons absorbe by the L-cones, for example, can be calculate either as the prouct of the raiance incient at the cornea by the L-cone funamental l or as the prouct of the raiance by the spectral transmission factor τ of the inert pigments an by the cone spectral sensitivity. For a Rayleigh-like match, given the raiances y, g, r of the lights an the sensitivity l y, l g, l r an m y, m g, m r of the L- an M-cones to the lights, for L-cones for M-cones l = rr + l gg l yy eq. 1 m = rr + mgg myy eq. 2 With matrices, the equations can be written in terms of ratios of raiances as m l r r lg r m g g y y = ly my eq.3 When l y, l g, l r an m y, m g, m r are known, the system has a unique solution for the ratios of raiances, provie the matrix is non-singular. In most real experiments, test an primary lights are not monochromatic, an the spectral raiance of the light shoul be multiplie by the spectral response l an m of the cone funamentals so as to accurately express the total energy absorbe by the cones. Then, equations 1 an 2 are written as

5 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 5 = + y y l g g l r r l eq.4 = + y y m g g m r r m eq.5 In orer to simplify further calculations, it is possible to normalize the raiance to unity as r0, g0 y0 an multiply by an ajusting factor α, β, or γ. = + γ β α y l g l r l eq.6 = + γ β α y m g m r m eq.7 Written with matrices in terms of ratios of multiplying factor, it is r m r l 0 0 g m g l 0 0 γ β α γ = y m y l 0 0 eq.8 Again, when l an m are known, the system has a unique solution for the ratios of multiplying factors α/γ an β/γ, as long as the matrix is non-singular. The inverse problem consists of recovering the sensitivities of the cones from the raiances obtaine for one Rayleigh match. In this case, the system is uner-etermine. We woul expect that with several Rayleigh-like matches, the inverse problem, i.e. recovering the sensitivities of the cones from multiple Rayleigh match settings woul become over-etermine. Moelling of small variations in L-cone an M-cone photopigments. We nee to reconstruct the cone funamentals, i.e. the spectral sensitivity of the cones at the entrance to the eye. Cone funamentals epen on - the low ensity spectral absorbance of the photopigment, - the effective optical ensity of the photopigment in the outer segment of the cones, - the macular pigment optical ensity, - the lens optical ensity.

6 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 6 It has been shown empirically that the low ensity spectral absorbance of the photopigment can be represente by a unique template or nomogram Dartnall, 1953, when plotte on a logarithmic scale versus some function of wavelength: either the logarithmic value of the wavelength Stockman & Sharpe, 2000, or the normalize wavelength Govarovskii, 2000, or the fourth root of wavelength Barlow, 1982; Dartnall et al., 1983, or the logarithmic value of the wave number Baylor et al. 1987, or the normalize frequency Lamb, 1995; Carroll et al., The nomogram rule greatly facilitates the moelling since the low ensity spectral absorbance of an unknown photopigment can be simply calculate from the wavelength of peak sensitivity uniquely. Here, to create caniate M-cone an L-cone photopigments, we use the low ensity spectral absorbance curves of the M-cone an L-cone photopigments respectively, in log quantal units an in 5 nm steps, as propose by Stockman an Sharpe Then we shifte each template along a logarithmic wavelength scale, interpolate in 1 nm steps using the cubic spline metho, transforme to a linear scale, an obtaine the low ensity spectral absorbance curves A 0 of each caniate photopigment. Increasing the ensity of the photopigment in the cone outer segments broaens the banwith of the relative spectral absorptance. We calculate the cone absorptance spectrum, i.e. the fraction F q of the incient photons absorbe at the retinal level in cones of effective optical ensity D as F q DA0 = 1 10 eq.9 The lens an the macular pigment absorb light as two filters in front of the photoreceptors. So at the corneal level the cone funamental sensitivity is obtaine by multiplying the cone absorptance spectrum F q by the transmittance of the macular pigment τ mac an by the transmittance of the lens τ lens l q = F q,l-cones τ mac τ lens eq.10

7 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 7 m q = F q,m-cones τ mac τ lens eq.11 which can be converte then to energy units as l an m in orer to be inclue in a prouct with the spectral raiance of the incient light measure in energy. roceure As explaine previously, we expect that with several Rayleigh-like matches, it will be possible to recover the sensitivities of the cones from the raiances of the matches. The avantage of moelling the spectral sensitivities of the cone funamentals using a template for the photopigment sensitivity is that all l values, as all m values, are linke to 2 inepenent variables the wavelength of the peak sensitivity of the M-cone photopigment, the wavelength of the peak sensitivity of the L-cone photopigment, an 3 parameters the effective optical ensity of the photopigments, the lens optical ensity, the macular pigment optical ensity. Thus, we propose the following proceure: 1. First, the observer performe several real matches an obtaine several values for r / y an g / y, from which we use only the ratio of raiances r / g 2. Secon, we preicte the ratio of raiances r / g for a range of plausible values in terms of wavelength of the peak sensitivity of the M-cone an L-cone photopigment from 515 to 535 nm an from 540 to 570 nm respectively, in 1 nm steps. We coul simplify the proceure by fixing at typical values the effective optical ensity of the photopigment in the outer segment of the cones, the lens optical an the macular pigment optical ensity. 3. Finally, we programme in Matlab an algorithm to minimise the sum of square ifferences between observe an preicte values of the ratio of raiances r / g.

8 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 8 Experimental setup Observers were aske to perform colour-matches between a yellow test light an an aitive mixture of re an green lights. We chose 5 wavelengths in the mile- an long-wave range of the visible spectrum. Among the 10 possible layouts that coul be achieve with 5 lights, only 7 that yiele photopic luminance levels were assesse, with a yellow light as the test light 1. Y1586nm + G551nm = Y2580nm, at about c.m R1617nm + G551nm = Y2580nm, at about c.m R3639nm + G551nm = Y2580nm, at about c.m R1617nm + G551nm = Y1586nm, at about c.m R3639nm + G551nm = Y1586nm, at about c.m R1617nm + Y2580nm = Y1586nm, at about c.m R3639nm + Y2580nm = Y1586nm, at about c.m -2

9 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 9 Fig. 1. Layout of the three-channel optical system use for multiple Rayleigh matches. The moifie Maxwellian view layout with iffusers inserte before the Lummer cube allows enlarging the exit beam at the entrance of the pupil so that the observer coul see the 10 eg matching fiel in nearly natural viewing. One neutral wege allows ajusting the amount of "Re" light in the Re an Green mixture. Another neutral wege allows ajusting the amount of "Yellow" light. A three-channel optical system Fig. 1 was use an calibrate. A moifie Maxwellian view layout iffusers were inserte before the Lummer cube allowe enlarging the exit beam at the entrance of the pupil so that the observer coul see the 10 eg matching fiel in nearly natural viewing. Compare to the stanar 2 match, the 10 eg match offers some comfort to the observer an improves the repeatability of his matches. A im white surroun allowe for stabilising the state of aaptation of the observer. The quasi-monochromatic test an primary lights were obtaine using interference filters. One neutral wege allowe ajusting the amount of "Re" light an controlling the hue of the Re an Green mixture. Another neutral wege allowe ajusting the amount of "Yellow" light an controlling the brightness equilibrium of the Rayleigh match. Spectroraiometric calibration was achieve in situ. Subjects Fourteen colour normal observers 5 men, 9 women were aske to perform 7 Rayleigh-like matches 6 repetitions for 12 observers, 15 repetitions for observers FV an F. Their colour vision was screene using Ishihara plates, D15 an Nagel anomaloscope Table 1. Anomaloscope settings were normal for all but two observers CD an who set the regreen mixture scale of the anomaloscope at the borer between normal an protanomalous ranges. A euteranomalous observer serve as a control subject.

10 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 10 Table 1. Description of the group of observers. Normal setting on the anomaloscope is An anomaloscope setting at 46 yiels an anomalous quotient equal to 0.74 an efines the borer between normal an protanomalous observers. Known inheritance is reporte in the text. Wavelengths of M-cone an L-cone peak sensitivity estimates are obtaine using tabulate photopigment spectral absorbance. Observer Sex Age Ishihara anel Nagel Known M-cone L-cone D15 anomaloscope inheritance max max CD F 28 Normal Normal 46 YES DO M Normal Normal Normal FR F 58 Normal Normal Normal YES FV F 58 Normal Normal LS F 23 Normal Normal MH F 23 Normal Normal NB H 27 Deutan Deutan NG F 21 Normal + Normal YES NK M 28 Normal Normal Normal M 21 Normal Normal F M 30 Normal Normal Normal RD M 59 Normal Normal Normal SJV F 27 Normal + Normal & 45 SV F 23 Normal Normal Normal VB F 44 Normal Normal Normal YES

11 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 11 Results The algorithm allowe us to ajust the moel parameters of the putative factors of variation in orer to minimize the sum of square ifferences between preicte an real results of matches. rovisionally, we fixe at typical values Stockman & Sharpe, 2000 the effective optical ensity of the photopigment in the outer segment of the cones 0.38 for a 10 eg fiel, the lens optical ensity at 400 nm an the macular pigment optical ensity at 460 nm for a 10 eg fiel. Fig. 2. 3D mesh plot of the sum of square ifferences between observe an preicte values of the ratio of raiances r / g for observer F.

12 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 12 Fig. 2 shows a 3D mesh plot of the sum of square ifferences between observe an preicte values of the ratio of raiances r / g for one observer. We notice the elongate shape of the hollow region. Fig. 3 isplays the contours of the sum of square ifferences between observe an preicte values of the ratio of raiances r / g for this observer. Here, a minimum is visible at Optimum max,m-cones = nm Optimum max,l-cones = nm Fig. 3. Contours of the sum of square ifferences between observe an preicte values of the ratio of raiances r / g for observer F.

13 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 13 For all observers we foun qualitatively similar contours as in figure 3, with the elongate shape of the hollow space. Nevertheless, the minimisation proceure yiels a unique solution for the wavelength of the peak sensitivities of the M- an L-cone photopigments. Fig. 4. Ajuste wavelengths of the L-cone an M-cone peak sensitivity of all observers. Square symbols for male observers. Diamon symbols for female observers. Only results of observers whose initials are mentione in the text are labelle. Thirteen observers lie within the cluster efine by a circle.

14 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 14 In fig. 4 we present the results of all observers on one graph. Remarkably, observer CD shows a result positione about the upper left corner of the graph, which means that the optimum peak wavelength for the M- an L-photopigments are close to each other. All other normal observers results lie in a cluster about 8 nm wie. As an aitional check of the multiple Rayleigh match proceure, a euteranomalous observer NB on the Nagel anomaloscope was teste. Since the caniate abnormal photopigment can be erive from the M- or the L-template, the solution epens on whether the anomalous photopigment is base on the M- or the L- template. Assuming a small ifference in wavelength of peak sensitivity between the anomalous M-cone photopigment an the normal L-cone photopigment DeMarco et al., 1992, we erive the anomalous photopigment spectral absorbance from the L-cone pigment template. Then, the euteranomalous observer s multiple Rayleigh match settings coul be preicte with a pair of L-cone max = nm an intermeiate cone photopigments max = nm. Discussion Experimental results showe that multiple Rayleigh matches allow ientifying the pair of photopigments that satisfy the iniviual settings. We also examine whether the istribution of the optimum wavelength of peak sensitivity reveals iniviual ifferences in M- an L- cone photopigments. After iscussing the sources of uncertainty, we will concentrate on iniviual ifferences. Sources of uncertainty The raiometric calibration of the equipment, the assume input values of the fixe parameters an the choice of photopigment spectral absorbance templates, are sources of uncertainty which may be of importance consiering the small number of Rayleigh-like equations an the limite repeatability of the observer s settings.

15 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 15 We checke that propagation of error of the spectral position of the maximum transmittance of the filters as little as 1 nm results in a shift of 1 or 2 nm of the L-cone photopigment sensitivity. Although the matches were easy to perform, an although the observations were reproucible, we face several ifficulties to erive iniviual factors of variation. At first, we ignore four factors of variation: the macular pigment that has little effect on large fiels beyon 551 nm, the peak optical ensity of the visual pigments, the peak optical ensity of the macular pigment, an the optical ensity of the lens. In a previous stuy Viénot, 2001, we showe that, with careful control of the colour-matching functions, the preiction of all factors of iniviual variation can be precise. Here, since our observers performe 7 equations, we moifie the software in orer to ajust, in aition to the wavelength of peak sensitivity of the photopigments, the optical ensity of the lens, or the effective optical ensity of the M-cone photopigment, or the effective optical ensity of the L-cone photopigment. Inee, ifferent estimations of the effective optical ensity of the M-cone an L-cone photopigments have alreay been accounte for in the literature Smith et al., Three observers were oler than 58 years FR, FV, RD. Observer FV lies at the borer of the cluster, the other two insie. We verifie that multiplying the lens optical ensity by 1.56 as for a 60 years ol observer okorny, Smith & Lutze, 1976 change only moerately the spectral position of the M-cone photopigment less than 1.8 nm an the L-cone photopigment less than 0.7 nm. In aition, it is known that a change of effective optical ensity of the photopigment affects notably the ajuste ratio of raiances in Rayleigh-type matches. In the case of the ual match which resembles equation #5 in our experimental stuy, He an Shevell 1994 calculate that the effect of a 5 nm shift in the L-cone funamental on the ratio of two particular Rayleigh-type matches is about as large as the effect of the normal extreme

16 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 16 variation in optical ensity an pre-receptoral filtering. Thomas an Mollon 2004 insiste on the impossibility to ientify a unique pair of photopigments that satisfy a given Rayleigh match because the match is etermine by the max of the photopigment an also its optical ensity. Here, we evaluate the observers who mae the largest number of repetitions. As anticipate, optimizing the effective optical ensity of a single photopigment results in a noticeable spectral shift up to 9 nm for the M-cone photopigment an up to 2.5 nm for the L- cone photopigment. Although ajusting aitional parameters allows reuction of the resiuals of the fit, such results must be consiere with caution as the system of 7 equations with 3 or more unknowns becomes only slightly over-etermine with respect to the variability of the observer s experimental settings. The choice of the template has a consierable influence on the final convergence of the algorithm. Alternatively, we moele the optical ensity spectrum t x of the ilute photopigment using analytical expressions. From Stockman an Sharpe 2000, we use log10t x = a + bx 2 + cx 4 + x 6 + ex 8 + fx 10 + gx 12 + hx 14 eq.12 where is the wavelength, max is the wavelength at peak sensitivity a, b, c,, e, f, g, h are the constant terms propose by Stockman an Sharpe 2000 x = log10 *558 / max. From Lamb 1995, we use S x = expa A x + expb B x + expc C x + D 1 where a, b, c, A, B, C, D are the constant terms propos by Lamb 1995 x = max /.

17 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 17 Then the photopigment sensitivity only epens on the wavelength at peak sensitivity max. Surprisingly, the spectral positions of the M-cone an the L-cone photopigment maximal absorbance collapse to the same peak wavelengths for 3 female observers CD, MH, FV when using the Lamb s template an for two observers CD, MH when using Stockman & Sharpe s template. For the remaining 11 observers, the algorithm yiels shifte wavelengths of the M-cone an L-cone photopigment peak sensitivity, epening upon the choice of template Fig. 5. Respective ranges are [519.5 : 524.3] an [555.4 : 561.9] with the tabulate templates, [529.3 : 538.8] an [556.7 : 567.8] with Stockman an Sharpe s polynomial template, [523.7 : 533.7] an [547.5 : 559] with Lamb s analytical template. Such isagreements reflect variations in the esign of the template, mainly in the positioning of the maximum of the template with respect to the slope at long wavelengths. The template that serves to moel photopigment absorbance is as critical to change the slope at long wavelengths as the cone effective optical ensity.

18 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 18 Fig. 5. Ajuste wavelengths of the L-cone an M-cone peak sensitivity are higher using Stockman an Sharpe s polynomial template than using their tabulate template. Figure 5 also shows that the range of M-cone an L-cone photopigment sensitivity is more restricte with the tabulate template than with Stockman an Sharpe s polynomial template. Thus, an accurate template is necessary for the multiple Rayleigh match proceure to converge. In the ignorance of a template base on soun theoretical principles, our ecision to relate the M-cone pigment an the L-cone pigment to inepenent patterns is probably wise.

19 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 19 Iniviual ifferences For all observers, the system of 7 equations converges to a unique solution. The istribution of the wavelength of maximum sensitivity of the photopigments inclues a cluster with 13 observers as shown in figure 4. Sources of uncertainty make it ifficult to rely on the absolute location of the wavelengths of peak sensitivity of the photopigments. Nevertheless, compare to the ata reporte in the literature, from for the normal M-cone photopigment to nm for the hybri M-cone photopigment an from for the hybri L-cone photopigment to nm for the normal L-cone photopigment after Sharpe et al, 1999, the multiple Rayleigh matches yiel somewhat shorter peak wavelengths for the M-cone photopigment: to nm an a quite plausible range for the L-cone photopigment: to nm. The case of observer CD nees to be explaine. rotanomaly coul be suspecte from the spectral position of her photopigments. However it is not obvious from her anomaloscope settings. We note that CD an have anomaloscope settings at the borer of the normal range but only CD s point is outsie of the cluster. Four observers are aware of eficiencies in their family. Observer CD has a colour blin uncle, from the mother s sie, an observer NG is aware of some colour vision efect for her father, with no more information, so some extra anomalous pigment might well be operative. Nevertheless, only CD is out of the cluster. Observer VB who has a euteranopic brother, an observer FR whose father an two sons are euteranopes lie insie the cluster. As expecte they are not carriers of an anomalous photopigment. Most instructive is the istribution between men an women within the cluster. All men s multiple Rayleigh matches yiel well separate L-cone an M-cone photopigments. For a few women FV an NG, an possibly MH, the minimum ifference between observe an preicte results occurs when both L an M optimum peak wavelengths move towar each

20 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 20 other. Observer NG was recruite because she coul rea the theoretically conceale igits of plate 45 of the Ishihara album without any hesitancy but observer SJV who coul rea conceale igits of plates 5 an 45 has well separate photopigments. We note that for FV an MH, algorithms using a unique analytical template for L-cone an M-cone photopigments converge to the same peak wavelength estimates for the two photopigments. Shifts of L-cone an M-cone optimum peak wavelength seem to be correlate, as if the pigments were attracte by a unique generic photopigment. Several explanations can be given. One possibility is that the variability of the matches is too high to allow the fit to converge. Alternatively, the Rayleigh-type colour matches are highly correlate an the equations cannot be regare as inepenent. Or, the assumption that a unique L-cone photopigment an a unique M-cone photopigment oversimplifies the analysis. Although genetics has shown that multiple pigments are not rare in the colour normal population: 8% for males Sjoberg et al., 1998, the pattern is expecte to be complex with women since about 15% of them inherit an X-chromosome carrying an abnormal opsin gene array from one parent Sharpe et al., Therefore, the female pattern is complex. The fact that in our reuce sample, only females ten to have close L-cone an M-cone photopigments might reflect multiple pigments. Conclusion In conclusion, with simple assumptions, multiple Rayleigh matches from a colour normal observer converge to a pair of L-cone an M-cone photopigments that satisfy the iniviual s settings. The locations of the wavelength of peak sensitivity of the L- an M- cone photopigments are preicte separately which represents an avantage for multiple Rayleigh matches versus ual Rayleigh matches. In our sample, the istribution of optimum preicte wavelengths of peak sensitivity inclues a cluster. The case of one observer is ambiguous. Nevertheless, in the face of the sources of error an uncertainty, it remains ifficult to

21 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 21 accurately preict the L-cone an M-cone photopigment wavelengths of maximum absorbance. Accurate tabulate templates are essential to ensure convergence to precise positioning of iniviual M-cone an L-cone photopigments.

22 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 22 Acknowlegments We thank colleagues an stuents who volunteer as observers. We acknowlege the contribution of Estelle Michon from the MST Optométrie at Université Orsay to a preliminary phase of this stuy, Ken Knoblauch an the reviewers for their valuable comments.

23 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 23 References Barlow, H.B What causes trichromacy? A theoretical analysis using comb-filtere spectra. Vision Research, 22, Baylor, D.A., Nunn, B. J. & Schnapf, J.L Spectral sensitivity of cones of the monkey Macaca fascicularis. Journal of hysiology, 390, Carroll, J., McMahon, C., Neitz, M. & Neitz, J Flicker-photometric electroretinogram estimates of L:M cone photoreceptor ratio in men with photopigment spectra erive from genetics. Journal of the Optical Society of America A, 17, Dartnall, H.J.A The interpretation of spectral sensitivity curves. British Meical Bulletin, 9, Dartnall, H.J.A., Bowmaker, J.K., & Mollon, J.D Human visual pigments: microspectrophotometric results from the eyes of seven persons. roceeings of the Royal Society, Lonon, Series B, 220, DeMarco,., okorny, J., & Smith, V.C Full spectrum cone sensitivity functions for X-chromosome linke anomalous trichromats. Journal of the Optical Society of America A, 9, Govarovskii V.I., Fyhrquist N., Reuter T., Kuzmin D.G. & Donner K In search of the visual pigment template. Vis Neurosci. 17: He, J.C. & Shevell, S.K Iniviual ifferences in cone photopigments of normal trichromats measure by ual Rayleigh-type color matches. Vision Research, 34, Lamb, T. D hotoreceptor spectral sensitivities: common shape in the longwavelength region. Vision Research, 35, Neitz, J. & Jacobs, G.H olymorphism in normal human color vision an its mechanisms. Vision Research, 30,

24 23/04/ :27igments e cônes_vns _vienotserraultaro_siteweb 24 okorny J., Smith V.C. & Lutze M Aging of the human lens. Appl. Optics 26: Sharpe, L. T., Stockman, A., Jägle, H., & Nathans, J Opsin genes, cone photopigments, color vision an colorblinness. In Color vision: From Genes to erception, e. Gegenfurtner K.R. & Sharpe L. T., pp Cambrige, UK: Cambrige University ress. Sjoberg, S.A., Neitz M., Baling S.D. & Neitz J L-cone pigment genes expresse in normal colour vision. Vision Research, 38: Smith V.C., okorny J. & DeMarco Full-spectrum cone sensitivity functions for X-chromosome-linke anomalous trichromats. Journal of the Optical Society of America A, 9, Smith V.C., okorny J. & Starr S.J Variability of color mixture ata I. Interobserver variability in the unit coorinates. Vision Research, 16, Stockman, A. & Sharpe, L. T Spectral sensitivities of the mile- an longwavelength sensitive cones erive from measurements in observers of known genotype. Vision Research, 40, Thomas.B.M. & Mollon J.D Moelling the Rayleigh match. Visual Neuroscience, 21, Viénot F Retinal istribution of the macular pigment an the cone effective optical ensity from colour matches of real observers. Color Research an Application, 26, S264- S268.

Modelling the Rayleigh match

Modelling the Rayleigh match Visual Neuroscience (2004), 21, 477 482. Printed in the USA. Copyright 2004 Cambridge University Press 0952-5238004 $16.00 DOI: 10.10170S095252380421344X Modelling the Rayleigh match P.B.M. THOMAS and

More information

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

OPTO 5320 VISION SCIENCE I

OPTO 5320 VISION SCIENCE I OPTO 5320 VISION SCIENCE I Monocular Sensory Processes of Vision: Color Vision Mechanisms of Color Processing VI. Retinal fundamentals A. Retinal fundamentals and cone photopigments B. Properties of cone

More information

MATH , 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208

MATH , 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208 MATH 321-03, 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208 Instructor: Brent Deschamp Email: brent.eschamp@ssmt.eu Office: McLaury 316B Phone:

More information

Impurities in inelastic Maxwell models

Impurities in inelastic Maxwell models Impurities in inelastic Maxwell moels Vicente Garzó Departamento e Física, Universia e Extremaura, E-671-Baajoz, Spain Abstract. Transport properties of impurities immerse in a granular gas unergoing homogenous

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

Final Exam Study Guide and Practice Problems Solutions

Final Exam Study Guide and Practice Problems Solutions Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

Code_Aster. Detection of the singularities and computation of a card of size of elements

Code_Aster. Detection of the singularities and computation of a card of size of elements Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : Josselin DLMAS Clé : R4.0.04 Révision : 9755 Detection of the singularities an computation of a car of size

More information

Code_Aster. Detection of the singularities and calculation of a map of size of elements

Code_Aster. Detection of the singularities and calculation of a map of size of elements Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : DLMAS Josselin Clé : R4.0.04 Révision : Detection of the singularities an calculation of a map of size of

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

The Press-Schechter mass function

The Press-Schechter mass function The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Sensors & Transducers 2015 by IFSA Publishing, S. L.

Sensors & Transducers 2015 by IFSA Publishing, S. L. Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity

More information

23 Implicit differentiation

23 Implicit differentiation 23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For

More information

NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS

NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS ANZIAM J. 46(2004), 1 16 NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS A. C. MCINTOSH 1,R.O.WEBER 2 ang.n.mercer 2 (Receive 14 August, 2002) Abstract This

More information

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer Key Engineering Materials Online: 4-8-5 I: 66-9795, Vols. 7-73, pp 38-33 oi:.48/www.scientific.net/kem.7-73.38 4 rans ech ublications, witzerlan Citation & Copyright (to be inserte by the publisher imulation

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

We G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies

We G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables* 51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

Quantile function expansion using regularly varying functions

Quantile function expansion using regularly varying functions Quantile function expansion using regularly varying functions arxiv:705.09494v [math.st] 9 Aug 07 Thomas Fung a, an Eugene Seneta b a Department of Statistics, Macquarie University, NSW 09, Australia b

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

Three-Dimensional Modeling of Green Sand and Squeeze Molding Simulation Yuuka Ito 1,a and Yasuhiro Maeda 2,b*

Three-Dimensional Modeling of Green Sand and Squeeze Molding Simulation Yuuka Ito 1,a and Yasuhiro Maeda 2,b* Materials Science Forum Submitte: 2017-08-13 ISSN: 1662-9752, Vol. 925, pp 473-480 Revise: 2017-12-09 oi:10.4028/www.scientific.net/msf.925.473 Accepte: 2018-01-12 2018 Trans Tech Publications, Switzerlan

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Extinction, σ/area. Energy (ev) D = 20 nm. t = 1.5 t 0. t = t 0

Extinction, σ/area. Energy (ev) D = 20 nm. t = 1.5 t 0. t = t 0 Extinction, σ/area 1.5 1.0 t = t 0 t = 0.7 t 0 t = t 0 t = 1.3 t 0 t = 1.5 t 0 0.7 0.9 1.1 Energy (ev) = 20 nm t 1.3 Supplementary Figure 1: Plasmon epenence on isk thickness. We show classical calculations

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

Role of parameters in the stochastic dynamics of a stick-slip oscillator

Role of parameters in the stochastic dynamics of a stick-slip oscillator Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

arxiv: v1 [hep-ex] 4 Sep 2018 Simone Ragoni, for the ALICE Collaboration

arxiv: v1 [hep-ex] 4 Sep 2018 Simone Ragoni, for the ALICE Collaboration Prouction of pions, kaons an protons in Xe Xe collisions at s =. ev arxiv:09.0v [hep-ex] Sep 0, for the ALICE Collaboration Università i Bologna an INFN (Bologna) E-mail: simone.ragoni@cern.ch In late

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

3.2 Shot peening - modeling 3 PROCEEDINGS

3.2 Shot peening - modeling 3 PROCEEDINGS 3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:

More information

Some vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10

Some vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10 Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the

More information

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Human color vision and the unsaturated blue color of the daytime sky

Human color vision and the unsaturated blue color of the daytime sky Human color vision and the unsaturated blue color of the daytime sky Glenn S. Smith a) School of Electrical and Computer Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0250 Received

More information

Improved Geoid Model for Syria Using the Local Gravimetric and GPS Measurements 1

Improved Geoid Model for Syria Using the Local Gravimetric and GPS Measurements 1 Improve Geoi Moel for Syria Using the Local Gravimetric an GPS Measurements 1 Ria Al-Masri 2 Abstract "The objective of this paper is to iscuss recent research one by the author to evelop a new improve

More information

Solutions to Math 41 Second Exam November 4, 2010

Solutions to Math 41 Second Exam November 4, 2010 Solutions to Math 41 Secon Exam November 4, 2010 1. (13 points) Differentiate, using the metho of your choice. (a) p(t) = ln(sec t + tan t) + log 2 (2 + t) (4 points) Using the rule for the erivative of

More information

Tutorial on Maximum Likelyhood Estimation: Parametric Density Estimation

Tutorial on Maximum Likelyhood Estimation: Parametric Density Estimation Tutorial on Maximum Likelyhoo Estimation: Parametric Density Estimation Suhir B Kylasa 03/13/2014 1 Motivation Suppose one wishes to etermine just how biase an unfair coin is. Call the probability of tossing

More information

Cascaded redundancy reduction

Cascaded redundancy reduction Network: Comput. Neural Syst. 9 (1998) 73 84. Printe in the UK PII: S0954-898X(98)88342-5 Cascae reunancy reuction Virginia R e Sa an Geoffrey E Hinton Department of Computer Science, University of Toronto,

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

Journal of Engineering Science and Technology Review 7 (1) (2014) Research Article

Journal of Engineering Science and Technology Review 7 (1) (2014) Research Article Jestr Journal of Engineering Science an Technology Review 7 () (04) 8 Research Article JOURNA OF Engineering Science an Technology Review www.jestr.org Research on the Measurement Error of MWIR Average

More information

Make graph of g by adding c to the y-values. on the graph of f by c. multiplying the y-values. even-degree polynomial. graph goes up on both sides

Make graph of g by adding c to the y-values. on the graph of f by c. multiplying the y-values. even-degree polynomial. graph goes up on both sides Reference 1: Transformations of Graphs an En Behavior of Polynomial Graphs Transformations of graphs aitive constant constant on the outsie g(x) = + c Make graph of g by aing c to the y-values on the graph

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

Non-Linear Bayesian CBRN Source Term Estimation

Non-Linear Bayesian CBRN Source Term Estimation Non-Linear Bayesian CBRN Source Term Estimation Peter Robins Hazar Assessment, Simulation an Preiction Group Dstl Porton Down, UK. probins@stl.gov.uk Paul Thomas Hazar Assessment, Simulation an Preiction

More information

Summary: Differentiation

Summary: Differentiation Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin

More information

Problem Set 2: Solutions

Problem Set 2: Solutions UNIVERSITY OF ALABAMA Department of Physics an Astronomy PH 102 / LeClair Summer II 2010 Problem Set 2: Solutions 1. The en of a charge rubber ro will attract small pellets of Styrofoam that, having mae

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Model for Dopant and Impurity Segregation During Vapor Phase Growth

Model for Dopant and Impurity Segregation During Vapor Phase Growth Mat. Res. Soc. Symp. Proc. Vol. 648, P3.11.1-7 2001 Materials Research Society Moel for Dopant an Impurity Segregation During Vapor Phase Growth Craig B. Arnol an Michael J. Aziz Division of Engineering

More information

Wiener Deconvolution: Theoretical Basis

Wiener Deconvolution: Theoretical Basis Wiener Deconvolution: Theoretical Basis The Wiener Deconvolution is a technique use to obtain the phase-velocity ispersion curve an the attenuation coefficients, by a two-stations metho, from two pre-processe

More information

Mathematical Review Problems

Mathematical Review Problems Fall 6 Louis Scuiero Mathematical Review Problems I. Polynomial Equations an Graphs (Barrante--Chap. ). First egree equation an graph y f() x mx b where m is the slope of the line an b is the line's intercept

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica

More information

Modeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions. Paul M. Mwasame, Norman J. Wagner, Antony N.

Modeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions. Paul M. Mwasame, Norman J. Wagner, Antony N. Submitte to the Journal of Rheology Moeling the effects of polyispersity on the viscosity of noncolloial har sphere suspensions Paul M. Mwasame, Norman J. Wagner, Antony N. Beris a) epartment of Chemical

More information

Computing Derivatives J. Douglas Child, Ph.D. Rollins College Winter Park, FL

Computing Derivatives J. Douglas Child, Ph.D. Rollins College Winter Park, FL Computing Derivatives by J. Douglas Chil, Ph.D. Rollins College Winter Park, FL ii Computing Inefinite Integrals Important notice regaring book materials Texas Instruments makes no warranty, either express

More information

Lecture Notes: March C.D. Lin Attosecond X-ray pulses issues:

Lecture Notes: March C.D. Lin Attosecond X-ray pulses issues: Lecture Notes: March 2003-- C.D. Lin Attosecon X-ray pulses issues: 1. Generation: Nee short pulses (less than 7 fs) to generate HHG HHG in the frequency omain HHG in the time omain Issues of attosecon

More information

Moist Component Potential Vorticity

Moist Component Potential Vorticity 166 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 60 Moist Component Potential Vorticity R. MCTAGGART-COWAN, J.R.GYAKUM, AND M. K. YAU Department of Atmospheric an Oceanic Sciences, McGill University, Montreal,

More information

Resilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study

Resilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study Resilient Moulus Preiction Moel for Fine-Graine Soils in Ohio: Preliminary Stuy by Teruhisa Masaa: Associate Professor, Civil Engineering Department Ohio University, Athens, OH 4570 Tel: (740) 59-474 Fax:

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

Characterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy

Characterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy JOURNAL OF APPLIED PHYSICS VOLUME 85, NUMBER 1 15 JUNE 1999 Characterization of lea zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy Haifeng Wang, Wenhua Jiang, a) an Wenwu

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

Average value of position for the anharmonic oscillator: Classical versus quantum results

Average value of position for the anharmonic oscillator: Classical versus quantum results verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive

More information

Integration Review. May 11, 2013

Integration Review. May 11, 2013 Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.

More information

The new concepts of measurement error s regularities and effect characteristics

The new concepts of measurement error s regularities and effect characteristics The new concepts of measurement error s regularities an effect characteristics Ye Xiaoming[1,] Liu Haibo [3,,] Ling Mo[3] Xiao Xuebin [5] [1] School of Geoesy an Geomatics, Wuhan University, Wuhan, Hubei,

More information

Differentiability, Computing Derivatives, Trig Review. Goals:

Differentiability, Computing Derivatives, Trig Review. Goals: Secants vs. Derivatives - Unit #3 : Goals: Differentiability, Computing Derivatives, Trig Review Determine when a function is ifferentiable at a point Relate the erivative graph to the the graph of an

More information

Entanglement is not very useful for estimating multiple phases

Entanglement is not very useful for estimating multiple phases PHYSICAL REVIEW A 70, 032310 (2004) Entanglement is not very useful for estimating multiple phases Manuel A. Ballester* Department of Mathematics, University of Utrecht, Box 80010, 3508 TA Utrecht, The

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction

FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

Mathcad Lecture #5 In-class Worksheet Plotting and Calculus

Mathcad Lecture #5 In-class Worksheet Plotting and Calculus Mathca Lecture #5 In-class Worksheet Plotting an Calculus At the en of this lecture, you shoul be able to: graph expressions, functions, an matrices of ata format graphs with titles, legens, log scales,

More information

SYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d

SYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d SYNTHESIS ON THE ASSESSMENT OF HIPS ONTRATION OEFFIIENT Marius MILEA Technical University Gheorghe Asachi of Iași, Machine Manufacturing an Inustrial Management Faculty, B-ul D Mangeron, nr 59, A, 700050,

More information

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical

More information

Shape Effect on Blind Frequency for Depth Inversion in Pulsed Thermography

Shape Effect on Blind Frequency for Depth Inversion in Pulsed Thermography Shape Effect on Blin Frequency for Depth Inversion in Pulse Thermography M. Genest 1, E. Grinzato 2, P. Bison 2, S. Marinetti 2 C. Ibarra-Castaneo 1, X. Malague 1 1 Electrical an Computing Eng. Dept.,

More information

The maximum sustainable yield of Allee dynamic system

The maximum sustainable yield of Allee dynamic system Ecological Moelling 154 (2002) 1 7 www.elsevier.com/locate/ecolmoel The maximum sustainable yiel of Allee ynamic system Zhen-Shan Lin a, *, Bai-Lian Li b a Department of Geography, Nanjing Normal Uni ersity,

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

CONTROL CHARTS FOR VARIABLES

CONTROL CHARTS FOR VARIABLES UNIT CONTOL CHATS FO VAIABLES Structure.1 Introuction Objectives. Control Chart Technique.3 Control Charts for Variables.4 Control Chart for Mean(-Chart).5 ange Chart (-Chart).6 Stanar Deviation Chart

More information

COUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS

COUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

A new identification method of the supply hole discharge coefficient of gas bearings

A new identification method of the supply hole discharge coefficient of gas bearings Tribology an Design 95 A new ientification metho of the supply hole ischarge coefficient of gas bearings G. Belforte, F. Colombo, T. Raparelli, A. Trivella & V. Viktorov Department of Mechanics, Politecnico

More information

Sources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences

Sources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences Sources an Sinks of Available Potential Energy in a Moist Atmosphere Olivier Pauluis 1 Courant Institute of Mathematical Sciences New York University Submitte to the Journal of the Atmospheric Sciences

More information

A simple tranformation of copulas

A simple tranformation of copulas A simple tranformation of copulas V. Durrleman, A. Nikeghbali & T. Roncalli Groupe e Recherche Opérationnelle Créit Lyonnais France July 31, 2000 Abstract We stuy how copulas properties are moifie after

More information

A USER-FRIENDLY TOOL FOR SIMULATING THE TIME-DEPENDENT FIELD DISTRIBUTION IN PE INSULATION ON THE BASIS OF A PHYSICAL APPROACH

A USER-FRIENDLY TOOL FOR SIMULATING THE TIME-DEPENDENT FIELD DISTRIBUTION IN PE INSULATION ON THE BASIS OF A PHYSICAL APPROACH A USER-FRIENLY TOOL FOR SIMULATING THE TIME-EPENENT FIEL ISTRIBUTION IN PE INSULATION ON THE BASIS OF A PHYSICAL APPROACH Séverine LE ROY, Laboratoire Plasma et Conversion Energie, (France), severine.leroy@laplace.univ-tlse.fr

More information

COLOR SCIENCE. Concepts and Methods, Quantitative Data and Formulae, 2nd Edition. John Wiley & Sons New York Chichester Brisbane Toronto Singapore

COLOR SCIENCE. Concepts and Methods, Quantitative Data and Formulae, 2nd Edition. John Wiley & Sons New York Chichester Brisbane Toronto Singapore COLOR SCIENCE Concepts and Methods, Quantitative Data and Formulae, 2nd Edition GÜNTER WYSZECKI National Research Council, Ottawa, Ontario, Canada W. S. STILES Richmond, Surrey, England t^- n M 1982 A

More information

MATHEMATICAL REPRESENTATION OF REAL SYSTEMS: TWO MODELLING ENVIRONMENTS INVOLVING DIFFERENT LEARNING STRATEGIES C. Fazio, R. M. Sperandeo-Mineo, G.

MATHEMATICAL REPRESENTATION OF REAL SYSTEMS: TWO MODELLING ENVIRONMENTS INVOLVING DIFFERENT LEARNING STRATEGIES C. Fazio, R. M. Sperandeo-Mineo, G. MATHEMATICAL REPRESENTATION OF REAL SYSTEMS: TWO MODELLING ENIRONMENTS INOLING DIFFERENT LEARNING STRATEGIES C. Fazio, R. M. Speraneo-Mineo, G. Tarantino GRIAF (Research Group on Teaching/Learning Physics)

More information

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything

More information

Lecture 16: The chain rule

Lecture 16: The chain rule Lecture 6: The chain rule Nathan Pflueger 6 October 03 Introuction Toay we will a one more rule to our toolbo. This rule concerns functions that are epresse as compositions of functions. The iea of a composition

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

PARALLEL-PLATE CAPACITATOR

PARALLEL-PLATE CAPACITATOR Physics Department Electric an Magnetism Laboratory PARALLEL-PLATE CAPACITATOR 1. Goal. The goal of this practice is the stuy of the electric fiel an electric potential insie a parallelplate capacitor.

More information

Space-time Linear Dispersion Using Coordinate Interleaving

Space-time Linear Dispersion Using Coordinate Interleaving Space-time Linear Dispersion Using Coorinate Interleaving Jinsong Wu an Steven D Blostein Department of Electrical an Computer Engineering Queen s University, Kingston, Ontario, Canaa, K7L3N6 Email: wujs@ieeeorg

More information