Beer's law with a polychromatic light source

As is well known, Beer’s law (aka the Beer-Lambert law, among others) is predicated on a number of conditions, one of which is that the light source be monochromatic or, as a practical approximation, that the spectral bandwidth of the light source be much smaller than the spectral bandwidths of whatever absorbers are being investigated. For present purposes, a polychromatic light source may be considered to have spectral bandwidth, $$\ce S_{source}$$, that is not negligible relative to that of the relevant absorbers. The question posed below will deal with what happens when the light source is decisively polychromatic.

Assume a cell path length, b, of 1 cm, and assume a single absorber, having triangular molar absorptivity spectral profile (Fig. 1): The peak molar absorptivity, $$\ce \epsilon_{450 nm}$$, is 100 L/mole cm at 450 nm. Assume the absorber concentration, c, ranges from 0 to 0.01 M. In the figure, the absorber’s molar absorptivity spectral profile is modeled as an isosceles triangle. In reality, absorber profiles are often modeled [1, p. 164] as Gaussian, Lorentzian, Voigt, or whatever the underlying physics requires. Here, however, an isosceles triangle is used for mathematical simplicity.

The figure also shows a narrower red isosceles triangle, centered (and having unity transmittance) at 425 nm and having full width at half maximum height (FWHM) equal to 20 nm. This red triangle is the transmittance spectral profile [2, equation 3-60, p. 70] of the light that exits a monochromator, assuming white light of uniform (unity) intensity is incident upon the monochromator and also assuming both 10 nm/mm reciprocal linear dispersion (RLD) and 2 mm slits ($$\ce W_{slits}$$). With white light incident on the entrance slits of the monochromator, the exiting light has $$\ce S_{source}$$ = RLD x $$\ce W_{slits}$$ = 20 nm. Note that $$\ce S_{source}$$ is defined here as the FWHM of the spectral profile: this is not the only possible such definition, but it is satisfactory for present purposes. Before proceeding, note that a 0.1 M monochromator typically has something like RLD = 4 nm/mm and 2 mm slits are relatively huge, as slit widths go, so the light source’s profile is polychromatic (and not good at all). Indeed, the FWHM is 20 nm, while that of the absorber is 50 nm.

As always, $$\ce T(\lambda) = I_{transmitted}(\lambda)/I_{incident}(\lambda)$$ and $$\ce A(\lambda) = -log[T(\lambda)]$$. If $$\ce c = c_{max} = 0.01 M$$ and b = 1 cm, then the next figure shows the absorbance (black) and transmittance (red) spectra from 380 to 520 nm (Fig. 2): If the light source was monochromatic, at 425 nm, a Beer’s law plot would be a straight line with slope = $$\ce \epsilon_{425 nm}$$ x b = 50 L/mole cm x 1 cm = $$50 M^{-1}$$. For absorber concentrations from 0 to 0.01 M, $$\ce A_{max} (425 nm)$$ = 0.5 at $$\ce c_{max} = 0.01 M$$.

However, for the light source profile in Fig. 1, the Beer’s law plot will curve downward, as undergrad textbooks like to show in cartoon form, and the absorbance at $$\ce c_{max}$$ will be less than 0.5.

Question: For the above scenario, with the specified polychromatic source profile in Fig. 1, what is the absorbance at $$\ce c_{max} = 0.01 M$$?

Hints: 1) Forget about undergrad textbooks for this. 2) There is no added noise or error. 3) The light source cannot be modeled as a pair of monochromatic lines.

Assume the photodetector responsivity is flat from 200 to 700 nm and dark response is zero. There is no detector or preamplifier noise, no digitization error, and no source shot noise. This is idealized in every reasonable way except for the polychromatic source, which, of course, no one would deliberately use.

References:

1. T.C. O’Haver, “Effect of the Source/Absorber Width Ratio on the Signal-to-Noise Ratio of Dispersive Absorption Spectrometry”, Analytical Chemistry, 63 (1991) 164-169.

2. J. D. Ingle, Jr. and S. R. Crouch, Spectrochemical Analysis, Prentice Hall, Englewood Cliffs, NJ, ©1988, pp. 70-74.

• Generally this sort of extended mathematical analysis is useless. As an analytical chemist I'd create a calibration curve. If it is non-linear then you need more standards to get a good fit. – MaxW Jun 21 at 21:39
• Sure! Using a polychromatic source obviously hurts sensitivity, peak resolution and linearity, but how much? – Ed V Jun 21 at 21:40
• Dear Ed, this is a very interesting question...is this a exam question you used to write for graduate students? Write the answer yourself. The quality of discussions is decreasing esp. in chemistry. This happens in every forum with time. – M. Farooq Jun 21 at 23:20
• @Tyberius I already know the answer (I think!) and thought it would be fun to see what would happen if I asked it here. In real research, it often happens that what is very well understood is right next door to something that was never even suspected. So it is just a little puzzle for anyone who might want to play with it and maybe someone will gain a small pleasure from a bit of unorthodox thinking about our old friend Beer's law. – Ed V Jun 21 at 23:27
• @M. Farooq It is not an exam question for the students I used to teach in the grad Spectroanalytical Chemistry course. I will post my answer in due course, but I am genuinely curious to see if an answer from someone else gets posted. And my reputation is now so large (!) that maybe I should put a bounty on the question! Anyway, your comment is on target and there does seem to be quite a few questions with either no answer or no accepted answer. I have no idea why that is the case. – Ed V Jun 21 at 23:33

There are several ways to answer the question, so here are my three. The first answer involves using a simple integer wavelengths approximation, i.e., noticing that Figs. 1 and 2 have 1 nm point spacing, and then using a spreadsheet to perform a summation. The answer explained in this part also serves to define the terms used in the second answer, which involves analytic calculation of the relevant overlap integrals. Then the third answer, not discussed in detail, uses a Mueller calculus-based method. Before proceeding to the explanations, the numerical answer to the question is $$A_{maximum} = 0.46951302113867$$ at $$C_{max} = 0.01 M$$ and Fig. 4 shows the absorbance as a function of concentration: The absorbance is 93.9% of the ideal value (i.e., 0.5) at $$C_{max}$$, despite the source being significantly polychromatic, and the Beer’s law plot is convex, but perhaps less than might have been expected intuitively.

All three of the answers and explanations that follow are based on the model shown in Fig. 4 below: As per the question’s assumed conditions, the light intensity incident at the monochromator is unity, i.e., 1, with fictive units of Watts. Note that units may be ignored since they will divide out when transmittance is computed later. The 50:50 beam splitter sends half of the light it receives to the sample (upper) beam path and the other half of the light it receives goes to the reference (lower) beam path. Let $$I_{incident}(\lambda )$$ denote the light intensity incident upon the absorber and $$I_{sample}(\lambda )$$ denote the attenuated light intensity that gets through the absorber. Then $$I_{sample}(\lambda )$$ is the incident light intensity on the upper photodetector and transimpedance amplifier (PD-TIA).

From Fig. 4, it is clear that $$I_{reference}(\lambda )$$ = $$I_{incident}(\lambda )$$, assuming the mirror is ideal. This is simply the unattenuated incident light intensity on the lower PD-TIA. The photodiodes were stated to have flat responsivities, i.e., spectrally constant responsivities, so their photocurrent responses are independent of wavelength. Hence, the responsivities, in this idealized scenario, cancel out when transmittance is computed later. Likewise for the constant transimpedances, which are assumed equal. Hence, it suffices to simply ‘add up’ the light intensities, for each color of light, incident at each PD-TIA. Now to the answers and explanations.

Answer 1: As per Figs. 1 and 2, assume a wavelength spacing of 1 nm. There is no light outside the region 405 to 445 nm and wavelengths are effectively “discretized”, i.e., are treated as discrete integer values. Therefore, the number of discrete wavelengths, where light intensity is non-zero, is just 39, i.e., from 406 to 444 nm, in 1 nm increments. At 405 nm, $$I_{incident}(\lambda )$$ is 0 and it increases linearly to 0.5 (at 425 nm). It then decreases linearly to 0 (at 445 nm). The absorbance, $$A_{absorber}(\lambda)$$, increases linearly from 0.1 (at 405 nm) to 0.9 (at 445 nm). Hence, the transmittance through the absorber, $$T_{absorber}(\lambda)$$, decreases non-linearly from about 0.7943 (at 405 nm) to about 0.1259 (at 445 nm). As always, $$T = 10^{-A}$$, so there is no difficulty computing $$T_{absorber}(\lambda )$$ at any of the relevant wavelengths.

Accordingly, $$I_{sample}(\lambda )$$ = $$I_{incident}(\lambda )$$ x $$T_{absorber}(\lambda )$$ and $$I_{reference}(\lambda )$$ = $$I_{incident}(\lambda )$$ x $$T_{reference}(\lambda )$$ = $$I_{incident}(\lambda )$$, since $$T_{reference}(\lambda )$$ = 1. All of the necessary calculations are shown in the spreadsheet screenshot in Fig. 5 below: In Fig. 5, column A contains the wavelengths and column B is a unitless “x” variable that is relative to 405 nm, i.e., at 405 nm, x = 0 and so forth. Column C contains the incident light intensities, column D has the associated absorbances and column E has the corresponding transmittances through the absorber. Column F has the unity transmittances for the reference beam path, column G has the $$I_{sample}(\lambda )$$ values and column H has the $$I_{reference}(\lambda )$$ values. Adding up the 41 values in columns G and H yields total intensities: 3.39164 and 10.000, respectively. Then T is simply 0.339164 and $$A_{maximum} = 0.469590$$. This is close to the true value, as computed in the second answer, despite the rather coarse “discretization” of the wavelengths.

Answer 2: From the first answer, it is clear how the relevant overlap integrals should be calculated. First, x is the independent (unitless) wavelength variable and it is continuous from 0 to 40. As x increases from 0 to 20, $$I_{incident}(x)$$ increases from 0 to 0.5, so $$I_{incident}(x)$$ = x/40. As x increases from 20 to 40, $$I_{incident}(x)$$ decreases from 0.5 to 0, so $$I_{incident}(x)$$ = 1- (x/40). For the absorbances, $$A_{absorber}(x)$$ = 0.1 + 0.02x for x = 0 to 40. Hence, $$T_{absorber}(x)$$ = $$10^{0.1 + 0.02x}$$ and $$I_{sample}(x)$$ = $$I_{incident}(x)$$ x $$10^{0.1 + 0.02x}$$.

For the sample beam path, the total intensity, $$I_{sample}^{total}$$, is the following overlap integral: $$I_{sample}^{total} = \int\limits_0^{40} {I_{sample}(x)\ \mathrm dx}$$

But this is simply the sum of two definite integrals: $$\int\limits_0^{20} {(x/40)\ 10^{0.1 + 0.02x}\ \mathrm dx} + \int\limits_{20}^{40} {[1-(x/40)] 10^{0.1 + 0.02x}\ \mathrm dx}$$

Using the fact that $$10^{-x} = e^{-xln(10)}$$, and defining $$a \equiv -0.02 ln(10)$$, then results in $$\frac{e^{5a}}{40} \int\limits_{0}^{20} {xe^{ax} \mathrm dx} + e^{5a} \int\limits_{20}^{40} {e^{ax} \mathrm dx} - \frac{e^{5a}}{40} \int\limits_{20}^{40} {xe^{ax} \mathrm dx}$$

The integrals are easily evaluated, resulting in $$\frac{e^{5a}}{40a^2} {[(20a - 1) e^{20a} + 1]} + \frac{ e^{5a}}{a}{[e^{40a} – e^{20a}]} - \frac{e^{5a}}{40a^2} {[(40a - 1) e^{40a}] - [(20a - 1) e^{20a}]}$$

With $$a \equiv -0.02 ln(10)$$, $$I_{sample}^{total}$$ = 3.39224318670537.

For the reference beam path, the integrals are easier because $$T_{reference}(x)$$ = 1. Hence $$I_{reference}^{total} = \int\limits_0^{20} {(x/40) \ \mathrm dx} + \int\limits_{20}^{40} {[1-(x/40)] \ \mathrm dx} = 5 + 5 = 10$$ Alternatively, this is simply the area of the triangle: 0.5 x 40 x 0.5 = 10. Therefore, $$T = I_{sample}^{total}\ / \ I_{reference}^{total} = 3.39224318670537\ /\ 10 = 0.339224318670537$$. Then $$A_{maximum} = 0.46951302113867$$.

The analytical method worked, as expected, but the system was so heavily idealized that it was essentially a “spherical cow” and they don’t milk. Any realistic complications, e.g., even simply using a tungsten bulb light source, make the analytic method harder and it does not take much to make the equations all but analytically intractable. Then numerical methods may be employed, but the difficulty then mostly lies in getting the right model equations: correctly evaluating equations is not necessarily the same as evaluating correct equations.

What about convolution? Convolution, per se, is unnecessary in this heavily idealized problem, for the reasons discussed below. First, note that the spectral profile of $$I_{incident}(\lambda )$$ is the same as the monochromator’s transmittance profile (shown in Fig. 1), except for having a maximum value of 0.5 rather than unity. The spectral profile is symmetric. Second, the absorber’s transmittance profile is as shown in Fig. 2 and it is also symmetric. Third, the convolution of two real functions involves reversing one of the functions and then computing the overlap integrals as a function of a shift parameter, i.e., a relative displacement of the two functions. Since the two functions to be convolved are both symmetric, nothing needs to be reversed, so the convolution reduces to computing the overlap integrals as a function of a shift parameter, $$\tau$$. Zero displacement, i.e., $$\tau = 0$$, can be conveniently defined as when the two functions are perfectly aligned.

The next step is to compute the convolution function either analytically, if possible or desired, or numerically, if necessary or convenient. Then the convolution function would be evaluated at $$\tau = 25$$ (nm), since this corresponds to the 25 nm shift between the two spectral profiles and $$\tau$$ is integer valued, i.e., data values are available every nm, as shown in Fig. 5. With finer wavelength spacing, the shift would scale up accordingly. If the convolution is computed analytically, then $$\tau$$ is a real number.

For the data in Fig. 5, the numerically computed convolution is shown in Fig. 6: From Fig. 6, the value of the convolution at 425 nm is 3.39164, as in Fig. 5. If the two spectral profiles do not overlap, then the value of the convolution is 10, as expected. But note that all that was needed was the single overlap integral at the known $$\tau = 25$$ (nm), and that value is 3.39164, as shown above. For the reference beam path, the overlap integral was 10, as shown in the trivially simple calculation above.

This problem is so heavily idealized that there is no necessary reason to compute the convolution function. In a more realistic situation, the computation would be harder and even numerical evaluation of the convolution function might be fruitless.

Answer 3: Figure 4 appears to just be a schematic block diagram, but it is actually a screenshot of a Mueller optical calculus-based computer simulation program. Although discussion of the fundamental principles of the several optical calculuses is beyond the scope of this answer, both Mueller calculus and Jones calculus were invented in the early 1940’s and there is a considerable body of information about them in the literature, e.g., [1-4].

In the simulation depicted in Fig. 4, stacks of Stokes vectors (one per discrete color) are passed from one optical component to the next optical component(s), e.g., from the lamp block to the monochromator block. Each component has a resident stack of matrices, with components as appropriate. When an optical component receives an incoming stack of optical vectors, each vector is multiplied by the corresponding matrix and the resulting stack of processed optical vectors is sent out to the next optical component (or components) “downstream”. All the stacks of optical vectors are received, processed and then sent on in a single computer simulation step.

The plot shown in Fig. 3 was obtained by using 1,001 concentrations: 0 to 0.01 M, in increments of 10 µM. The wavelength range was from 405 to 445 nm, in 1 pm increments. Thus, with every one of the 1,001 simulation steps, stacks of 40,001 Stokes vectors (used in Mueller calculus) were passed, starting with the lamp block and ending with the sample and reference PD-TIA blocks. Both of those blocks, for every step, computed 40,001 intensities, multiplied them by responsivities, summed the resulting photocurrents and then multiplied by the transimpedances.

In the present case, the responsivities were simply constant default values (0.1 A/W) and the transimpedances were $$1 k\Omega$$. For every step of the simulation, T was computed as the quotient of the real number outputs of the two PD-TIA blocks and $$A_{maximum} = -log(T) = 0.46951302122$$. Compared with the exact answer, this differs in the ninth decimal place. As well, there was no need to solve any equations or use a spreadsheet: no wheel had to be re-invented.

Another of the advantages of using optical calculus-based spectrometric simulations is that they can easily accommodate non-idealities, such as realistic lamp spectral profiles and photodetector responsivities, imperfect beam splitters and mirrors, polarized light, optical activity, circular dichroism, stray light, light source modulation and so on. In fact, any linear optical effect can be modeled with an optical calculus. In addition, noise effects are easily included, e.g., detector/preamplifier noises and light source noises, including low frequency noises.

Summary: Answer 2 provided an exact solution to the posed question, but it was a bit tedious to obtain and would not be useful if the system was altered in any significant way, e.g., if the photodiode responsivities were those of a real red-sensitive photodiode. Answer 1 was adequately accurate, despite its simplicity, but not particularly amenable to generalizations. Answer 3 gave an accurate answer and made it simple to get the plot shown in Fig. 3. There is no significant difficulty in generalizing to more realistic scenarios.

Lastly, it is clear that using the decisively polychromatic source in Fig. 1 caused a relatively small reduction in $$A_{maximum}$$ and the Beer’s law plot is convex, but not badly so. The deleterious effect on spectral resolution was not investigated here, but could easily be done, perhaps as another question.

References:

1. W.A. Shurcliff, Polarized Light, Harvard U. Press, Cambridge, MA ©1962 & 1966.
2. D.S. Kliger, J.W. Lewis, C.E. Randall, Polarized Light in Optics and Spectroscopy, Academic Press, Boston ©1990.
3. Mueller calculus article in wikipedia and references therein.
4. Jones calculus article in wikipedia and references therein.
• Why use Mueller/Jones matrices when you are not overtly dealing with a polarisation problem ? The response you measure $I_{sample}$ will be the convolution of your instrument slit width (monochromator transmission in fig 1 of the question) with the molecular transmission profile as you scan the former over the latter. – porphyrin Jun 26 at 8:02
• Your answer 2 seems to be a straight integration not convolution. – porphyrin Jun 26 at 12:13
• There is no need it is a well understood process. – porphyrin Jun 26 at 14:12