Beer–Lambert law

The law was discovered by Pierre Bouguer before 1729, while looking at red wine, during a brief vacation in Alentejo, Portugal.[1] It is often attributed to Johann Heinrich Lambert, who cited Bouguer's Essai d'optique sur la gradation de la lumière (Claude Jombert, Paris, 1729) – and even quoted from it – in his Photometria in 1760.[2] Lambert's law stated that the loss of light intensity when it propagates in a medium is directly proportional to intensity and path length. Much later, the German scientist August Beer discovered another attenuation relation in 1852. Beer's law stated that the transmittance of a solution remains constant if the product of concentration and path length stays constant.[3] The modern derivation of the Beer–Lambert law combines the two laws and correlates the absorbance, which is the negative decadic logarithm of the transmittance, to both the concentrations of the attenuating species and the thickness of the material sample.[4] The first modern formulation was given possibly by Robert Luther and Andreas Nikolopulos in 1913.[5]

A common and practical expression of the Beer–Lambert law relates the optical attenuation of a physical material containing a single attenuating species of uniform concentration to the optical path length through the sample and absorptivity of the species. This expression is:

$${\displaystyle A=\varepsilon \ell c}$$

Where

• ${\displaystyle A}$ is the absorbance
• ${\displaystyle \varepsilon }$ is the molar attenuation coefficient or absorptivity of the attenuating species
• ${\displaystyle \ell }$ is the optical path length in cm
• ${\displaystyle c}$ is the concentration of the attenuating species

A more general form of the Beer–Lambert law states that, for ${\displaystyle N}$ attenuating species in the material sample,

$${\displaystyle T=e^{-\sum _{i=1}^{N}\sigma _{i}\int _{0}^{\ell }n_{i}(z)\mathrm {d} z}=10^{-\sum _{i=1}^{N}\varepsilon _{i}\int _{0}^{\ell }c_{i}(z)\mathrm {d} z},}$$

or equivalently that

$${\displaystyle \tau =\sum _{i=1}^{N}\tau _{i}=\sum _{i=1}^{N}\sigma _{i}\int _{0}^{\ell }n_{i}(z)\,\mathrm {d} z,}$$

$${\displaystyle A=\sum _{i=1}^{N}A_{i}=\sum _{i=1}^{N}\varepsilon _{i}\int _{0}^{\ell }c_{i}(z)\,\mathrm {d} z,}$$

where

• ${\displaystyle \sigma _{i}}$ is the attenuation cross section of the attenuating species ${\displaystyle i}$ in the material sample;
• ${\displaystyle n_{i}}$ is the number density of the attenuating species ${\displaystyle i}$ in the material sample;
• ${\displaystyle \varepsilon _{i}}$is the molar attenuation coefficient or absorptivity of the attenuating species ${\displaystyle i}$ in the material sample;
• ${\displaystyle c_{i}}$ is the amount concentration of the attenuating species ${\displaystyle i}$ in the material sample;
• ${\displaystyle \ell }$ is the path length of the beam of light through the material sample.

In the above equations, the transmittance ${\displaystyle T}$ of material sample is related to its optical depth ${\displaystyle {\tau }}$ and to its absorbance A by the following definition

$${\displaystyle T={\frac {\Phi _{\mathrm {e} }^{\mathrm {t} }}{\Phi _{\mathrm {e} }^{\mathrm {i} }}}=e^{-\tau }=10^{-A},}$$

where

• ${\displaystyle \Phi _{\mathrm {e} }^{\mathrm {t} }}$ is the radiant flux transmitted by that material sample;
• ${\displaystyle \Phi _{\mathrm {e} }^{\mathrm {i} }}$is the radiant flux received by that material sample.

Attenuation cross section and molar attenuation coefficient are related by

$${\displaystyle \varepsilon _{i}={\frac {\mathrm {N_{A}} }{\ln {10}}}\,\sigma _{i},}$$

and number density and amount concentration by

$${\displaystyle c_{i}={\frac {n_{i}}{\mathrm {N_{A}} }},}$$

where ${\displaystyle \mathrm {N_{A}} }$ is the Avogadro constant.

In case of uniform attenuation, these relations become[6]

$${\displaystyle T=e^{-\ell \sum _{i=1}^{N}\sigma _{i}n_{i}}=10^{-\ell \sum _{i=1}^{N}\varepsilon _{i}c_{i}},}$$

or equivalently

$${\displaystyle \tau =\ell \sum _{i=1}^{N}\sigma _{i}n_{i},}$$

$${\displaystyle A=\ell \sum _{i=1}^{N}\varepsilon _{i}c_{i}.}$$

Cases of non-uniform attenuation occur in atmospheric science applications and radiation shielding theory for instance.

The law tends to break down at very high concentrations, especially if the material is highly scattering. Absorbance within range of 0.2 to 0.5 is ideal to maintain linearity in the Beer–Lambart law. If the radiation is especially intense, nonlinear optical processes can also cause variances. The main reason, however, is that the concentration dependence is in general non-linear and Beer's law is valid only under certain conditions as shown by derivation below. For strong oscillators and at high concentrations the deviations are stronger. If the molecules are closer to each other interactions can set in. These interactions can be roughly divided into physical and chemical interactions. Physical interaction do not alter the polarizability of the molecules as long as the interaction is not so strong that light and molecular quantum state intermix (strong coupling), but cause the attenuation cross sections to be non-additive via electromagnetic coupling. Chemical interactions in contrast change the polarizability and thus absorption.

Expression with attenuation coefficient

The Beer–Lambert law can be expressed in terms of attenuation coefficient, but in this case is better called Lambert's law since amount concentration, from Beer's law, is hidden inside the attenuation coefficient. The (Napierian) attenuation coefficient ${\displaystyle \mu }$ and the decadic attenuation coefficient ${\displaystyle \mu _{10}=\mu /\ln 10}$ of a material sample are related to its number densities and amount concentrations as

$${\displaystyle \mu (z)=\sum _{i=1}^{N}\mu _{i}(z)=\sum _{i=1}^{N}\sigma _{i}n_{i}(z),}$$

$${\displaystyle \mu _{10}(z)=\sum _{i=1}^{N}\mu _{10,i}(z)=\sum _{i=1}^{N}\varepsilon _{i}c_{i}(z)}$$

respectively, by definition of attenuation cross section and molar attenuation coefficient. Then the Beer–Lambert law becomes

$${\displaystyle T=e^{-\int _{0}^{\ell }\mu (z)\mathrm {d} z}=10^{-\int _{0}^{\ell }\mu _{10}(z)\mathrm {d} z},}$$

and

$${\displaystyle \tau =\int _{0}^{\ell }\mu (z)\,\mathrm {d} z,}$$

$${\displaystyle A=\int _{0}^{\ell }\mu _{10}(z)\,\mathrm {d} z.}$$

In case of uniform attenuation, these relations become

$${\displaystyle T=e^{-\mu \ell }=10^{-\mu _{10}\ell },}$$

or equivalently

$${\displaystyle \tau =\mu \ell ,}$$

$${\displaystyle A=\mu _{10}\ell .}$$

In many cases, the attenuation coefficient does not vary with ${\displaystyle z}$, in which case one does not have to perform an integral and can express the law as:

$${\displaystyle I(z)=I_{0}e^{-\mu z}}$$

where the attenuation is usually an addition of absorption coefficient ${\displaystyle \alpha }$ (creation of electron-hole pairs) or scattering (for example Rayleigh scattering if the scattering centers are much smaller than the incident wavelength).[7] Also note that for some systems we can put ${\displaystyle 1/\lambda }$ (1 over inelastic mean free path) in place of ${\displaystyle \mu }$.[8]

Assume that a beam of light enters a material sample. Define z as an axis parallel to the direction of the beam. Divide the material sample into thin slices, perpendicular to the beam of light, with thickness dz sufficiently small that one particle in a slice cannot obscure another particle in the same slice when viewed along the z direction. The radiant flux of the light that emerges from a slice is reduced, compared to that of the light that entered, by dΦ_e(z) = −μ(z)Φ_e(z) dz, where μ is the (Napierian) attenuation coefficient, which yields the following first-order linear ODE:

$${\displaystyle {\frac {\mathrm {d} \Phi _{\mathrm {e} }(z)}{\mathrm {d} z}}=-\mu (z)\Phi _{\mathrm {e} }(z).}$$

The attenuation is caused by the photons that did not make it to the other side of the slice because of scattering or absorption. The solution to this differential equation is obtained by multiplying the integrating factor

$${\displaystyle e^{\int _{0}^{z}\mu (z')\mathrm {d} z'}}$$

throughout to obtain

$${\displaystyle {\frac {\mathrm {d} \Phi _{\mathrm {e} }(z)}{\mathrm {d} z}}\,e^{\int _{0}^{z}\mu (z')\mathrm {d} z'}+\mu (z)\Phi _{\mathrm {e} }(z)\,e^{\int _{0}^{z}\mu (z')\mathrm {d} z'}=0,}$$

which simplifies due to the product rule (applied backwards) to

$${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} z}}{\bigl (}\Phi _{\mathrm {e} }(z)\,e^{\int _{0}^{z}\mu (z')\mathrm {d} z'}{\bigr )}=0.}$$

Integrating both sides and solving for Φ_e for a material of real thickness ℓ, with the incident radiant flux upon the slice Φ_eⁱ = Φ_e(0) and the transmitted radiant flux Φ_e^t = Φ_e(ℓ ) gives

$${\displaystyle \Phi _{\mathrm {e} }^{\mathrm {t} }=\Phi _{\mathrm {e} }^{\mathrm {i} }\,e^{-\int _{0}^{\ell }\mu (z)\mathrm {d} z},}$$

and finally

$${\displaystyle T={\frac {\Phi _{\mathrm {e} }^{\mathrm {t} }}{\Phi _{\mathrm {e} }^{\mathrm {i} }}}=e^{-\int _{0}^{\ell }\mu (z)\mathrm {d} z}.}$$

Since the decadic attenuation coefficient μ₁₀ is related to the (Napierian) attenuation coefficient by μ₁₀ = μ/ln 10, one also have

$${\displaystyle T=e^{-\int _{0}^{\ell }\ln {10}\,\mu _{10}(z)\mathrm {d} z}={\bigl (}e^{-\int _{0}^{\ell }\mu _{10}(z)\mathrm {d} z}{\bigr )}^{\ln {10}}=10^{-\int _{0}^{\ell }\mu _{10}(z)\mathrm {d} z}.}$$

To describe the attenuation coefficient in a way independent of the number densities n_i of the N attenuating species of the material sample, one introduces the attenuation cross section σ_i = μ_i(z)/n_i(z). σ_i has the dimension of an area; it expresses the likelihood of interaction between the particles of the beam and the particles of the specie i in the material sample:

$${\displaystyle T=e^{-\sum _{i=1}^{N}\sigma _{i}\int _{0}^{\ell }n_{i}(z)\mathrm {d} z}.}$$

One can also use the molar attenuation coefficients ε_i = (N_A/ln 10)σ_i, where N_A is the Avogadro constant, to describe the attenuation coefficient in a way independent of the amount concentrations c_i(z) = n_i(z)/N_A of the attenuating species of the material sample:

$${\displaystyle {\begin{aligned}T=e^{-\sum _{i=1}^{N}{\frac {\ln {10}}{\mathrm {N_{A}} }}\varepsilon _{i}\int _{0}^{\ell }n_{i}(z)\mathrm {d} z}=\\\left(e^{-\sum _{i=1}^{N}\varepsilon _{i}\int _{0}^{\ell }{\frac {n_{i}(z)}{\mathrm {N_{A}} }}\mathrm {d} z}\right)^{\ln {10}}=10^{-\sum _{i=1}^{N}\varepsilon _{i}\int _{0}^{\ell }c_{i}(z)\mathrm {d} z}.\end{aligned}}}$$

Under certain conditions the Beer–Lambert law fails to maintain a linear relationship between attenuation and concentration of analyte. These deviations are classified into three categories:

1. Real—fundamental deviations due to the limitations of the law itself.
2. Chemical—deviations observed due to specific chemical species of the sample which is being analyzed.
3. Instrument—deviations which occur due to how the attenuation measurements are made.

There are at least six conditions that need to be fulfilled in order for the Beer–Lambert law to be valid. These are:

1. The attenuators must act independently of each other.
2. The attenuating medium must be homogeneous in the interaction volume.
3. The attenuating medium must not scatter the radiation—no turbidity—unless this is accounted for as in DOAS.
4. The incident radiation must consist of parallel rays, each traversing the same length in the absorbing medium.
5. The incident radiation should preferably be monochromatic, or have at least a width that is narrower than that of the attenuating transition. Otherwise a spectrometer as detector for the power is needed instead of a photodiode which cannot discriminate between wavelengths.
6. The incident flux must not influence the atoms or molecules; it should only act as a non-invasive probe of the species under study. In particular, this implies that the light should not cause optical saturation or optical pumping, since such effects will deplete the lower level and possibly give rise to stimulated emission.

If any of these conditions are not fulfilled, there will be deviations from the Beer–Lambert law.

The Beer–Lambert law can be applied to the analysis of a mixture by spectrophotometry, without the need for extensive pre-processing of the sample. An example is the determination of bilirubin in blood plasma samples. The spectrum of pure bilirubin is known, so the molar attenuation coefficient ε is known. Measurements of decadic attenuation coefficient μ₁₀ are made at one wavelength λ that is nearly unique for bilirubin and at a second wavelength in order to correct for possible interferences. The amount concentration c is then given by

$${\displaystyle c={\frac {\mu _{10}(\lambda )}{\varepsilon (\lambda )}}.}$$

For a more complicated example, consider a mixture in solution containing two species at amount concentrations c₁ and c₂. The decadic attenuation coefficient at any wavelength λ is, given by

$${\displaystyle \mu _{10}(\lambda )=\varepsilon _{1}(\lambda )c_{1}+\varepsilon _{2}(\lambda )c_{2}.}$$

Therefore, measurements at two wavelengths yields two equations in two unknowns and will suffice to determine the amount concentrations c₁ and c₂ as long as the molar attenuation coefficient of the two components, ε₁ and ε₂ are known at both wavelengths. This two system equation can be solved using Cramer's rule. In practice it is better to use linear least squares to determine the two amount concentrations from measurements made at more than two wavelengths. Mixtures containing more than two components can be analyzed in the same way, using a minimum of N wavelengths for a mixture containing N components.

The law is used widely in infra-red spectroscopy and near-infrared spectroscopy for analysis of polymer degradation and oxidation (also in biological tissue) as well as to measure the concentration of various compounds in different food samples. The carbonyl group attenuation at about 6 micrometres can be detected quite easily, and degree of oxidation of the polymer calculated.

This law is also applied to describe the attenuation of solar or stellar radiation as it travels through the atmosphere. In this case, there is scattering of radiation as well as absorption. The optical depth for a slant path is τ′ = mτ, where τ refers to a vertical path, m is called the relative airmass, and for a plane-parallel atmosphere it is determined as m = sec θ where θ is the zenith angle corresponding to the given path. The Beer–Lambert law for the atmosphere is usually written

$${\displaystyle T=e^{-m(\tau _{\mathrm {a} }+\tau _{\mathrm {g} }+\tau _{\mathrm {RS} }+\tau _{\mathrm {NO_{2}} }+\tau _{\mathrm {w} }+\tau _{\mathrm {O_{3}} }+\tau _{\mathrm {r} }+\cdots )},}$$

where each τ_x is the optical depth whose subscript identifies the source of the absorption or scattering it describes:

• a refers to aerosols (that absorb and scatter);
• g are uniformly mixed gases (mainly carbon dioxide (CO₂) and molecular oxygen (O₂) which only absorb);
• NO₂ is nitrogen dioxide, mainly due to urban pollution (absorption only);
• RS are effects due to Raman scattering in the atmosphere;
• w is water vapour absorption;
• O₃ is ozone (absorption only);
• r is Rayleigh scattering from molecular oxygen (O₂) and nitrogen (N₂) (responsible for the blue color of the sky);
• the selection of the attenuators which have to be considered depends on the wavelength range and can include various other compounds. This can include tetraoxygen, HONO, formaldehyde, glyoxal, a series of halogen radicals and others.

m is the optical mass or airmass factor, a term approximately equal (for small and moderate values of θ) to 1/cos θ, where θ is the observed object's zenith angle (the angle measured from the direction perpendicular to the Earth's surface at the observation site). This equation can be used to retrieve τ_a, the aerosol optical thickness, which is necessary for the correction of satellite images and also important in accounting for the role of aerosols in climate.

• Beer–Lambert Law Calculator
• Beer–Lambert Law Simpler Explanation