Beer's law with a polychromatic light source
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:


W.A. Shurcliff, Polarized Light, Harvard U. Press, Cambridge, MA ©1962 &amp; 1966.
D.S. Kliger, J.W. Lewis, C.E. Randall, Polarized Light in Optics and Spectroscopy, Academic Press, Boston ©1990.
Mueller calculus article in wikipedia and references therein.
Jones calculus article in wikipedia and references therein.