Beer–Lambert law

A demonstration of the Beer–Lambert law: green laser light in a solution of Rhodamine 6B. The beam radiant power becomes weaker as it passes through solution

The Beer–Lambert law, also known as Beer's law, the Lambert–Beer law, or the Beer–Lambert–Bouguer law relates the attenuation of light to the properties of the material through which the light is travelling. The law is commonly applied to chemical analysis measurements and used in understanding attenuation in physical optics, for photons, neutrons, or rarefied gases. In mathematical physics, this law arises as a solution of the BGK equation.

History

The law was discovered by Pierre Bouguer before 1729, while looking at red wine, during a brief vacation in Alentejo, Portugal. 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. 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, 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. 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.

Mathematical formulation

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:

$A=\varepsilon \ell c$

Where

• $\varepsilon$ is the molar attenuation coefficient or absorptivity of the attenuating species
• $\ell$ is the optical path length
• $c$ is the concentration of the attenuating species

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

$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

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

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

where

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

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

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

where

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

Attenuation cross section and molar attenuation coefficient are related by

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

and number density and amount concentration by

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

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

In case of uniform attenuation, these relations become

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

or equivalently

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

$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 the linearity in 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 $\mu$ and the decadic attenuation coefficient $\mu _{10}=\mu /\ln 10$ of a material sample are related to its number densities and amount concentrations as

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

$\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

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

and

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

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

In case of uniform attenuation, these relations become

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

or equivalently

$\tau =\mu \ell ,$

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

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

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

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

Derivation

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:

${\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

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

throughout to obtain

${\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

${\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

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

and finally

$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

$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:

$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:

${\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}=\\{\Bigl (}e^{-\sum _{i=1}^{N}\varepsilon _{i}\int _{0}^{\ell }{\frac {n_{i}(z)}{\mathrm {N_{A}} }}\mathrm {d} z}{\Bigr )}^{\ln {10}}=10^{-\sum _{i=1}^{N}\varepsilon _{i}\int _{0}^{\ell }c_{i}(z)\mathrm {d} z}.\end{aligned}}$

The above assumption that the attenuation cross sections are additive is generally incorrect since electromagnetic coupling occurs if the distances between the absorbing entities is small.

The derivation of the concentration dependence of the absorbance is based on electromagnetic theory. Accordingly, the macroscopic polarization of a medium $P$ derives from the microscopic dipole moments $p$ in the absence of interaction according to

$P=N\ p\$

where $p$ is the dipole moment and $N$ the number of absorbing entities per unit volume. On the other hand, macroscopic polarization is given by:

$P=(\varepsilon _{r}-1)\cdot \varepsilon _{0}\cdot E$

Here $\varepsilon _{r}$represents the relative dielectric function, $\varepsilon _{0}$ the vacuum permittivity and $E$ the electric field. After equating and solving for the relative dielectric function the result is:

$\varepsilon _{r}=1+{\frac {P}{\varepsilon _{0}\cdot E}}$

If we take into account that the polarizability $\alpha$ is defined by $p=\alpha \cdot E$ and that for the number of absorbers per unit volume $N=N_{A}\cdot c\$holds, it follows that:

$\varepsilon _{r}=1+c{\frac {N_{A}\cdot \alpha }{\varepsilon _{0}}}$

According to Maxwell's wave equation the following relation between the complex dielectric function and the complex index of refraction function holds $\varepsilon _{r}={\hat {n}}^{2}$for isotropic and homogeneous media. Therefore:

${\hat {n}}={\sqrt {1+c{\frac {N_{A}\cdot \alpha }{\varepsilon _{0}}}}}$

The imaginary part of the complex index of refraction is the index of absorption $k$. Employing the imaginary part of the polarizability $\alpha ''$and the approximation $\surd (1+x)\approx 1+x/2$ it follows that:

$k=c{\frac {N_{A}\cdot \alpha ''}{2\varepsilon _{0}}}$

Taking into account the relation between $k$ and $A$, $A=4\pi (\log _{10}e)k\cdot c\cdot d/\lambda$ it eventually follows that

$A={\frac {2\pi (\log _{10}e)N_{A}\alpha ''}{\lambda \cdot \varepsilon _{0}}}\cdot c\cdot d$

As a consequence, the linear relation between concentration and absorbance is generally an approximation, and holds in particular only for small polarisabilities and weak absorptions, i.e. oscillator strengths. If we do not introduce the approximation $\surd (1+x)\approx 1+x/2$, and employ instead the following relation between the imaginary part of the relative dielectric function and index of refraction and absorption $\varepsilon _{r}''=2nk$ it can be seen that the molar attenuation coefficient depends on the index of refraction (which is itself concentration dependent):

$A={\frac {2\pi (\log _{10}e)N_{A}\alpha ''}{n\cdot \lambda \cdot \varepsilon _{0}}}\cdot c\cdot d$

Validity

Under certain conditions 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 Beer–Lambert law to be valid. These are:

1. The attenuators must act independently of each other. Electromagnetic coupling must be excluded.
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 has not a selective wavelength dependence.
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.
7. The wave properties of light must be negligible. In particular interference enhancement or decrease must not occur.

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

The Beer–Lambert law is not compatible with Maxwell's equations. Being strict, the law does not describe the transmittance through a medium, but the propagation within that medium. It can be made compatible with Maxwell's equations if the transmittance of a sample with solute is ratioed against the transmittance of the pure solvent which explains why it works so well in spectrophotometry. As this is not possible for pure media, the uncritical employment of the Beer–Lambert law can easily generate errors of the order of 100% or more. In such cases it is necessary to apply the Transfer-matrix method. A detailed discussion of the incompatibility between the Beer–Lambert law and Maxwell's equations can be found in the review The Bouguer‐Beer‐Lambert law: Shining light on the obscure .

Recently it has also been demonstrated that Beer's law is a limiting law, since the absorbance is only approximately linearly depending on concentration. The reason is that the attenuation coefficient also depends on concentration and density, even in the absence of any interactions. These changes are, however, usually negligible except for high concentrations and large oscillator strength. For high concentrations and/or oscillator strengths, it is the integrated absorbance which is linearly depending on concentration, at least as long as there are no local field effects. If there are local field effects, they can be approximately taken into account by applying the Lorentz-Lorenz relation. In fact, Beer's law, i.e. the concentration dependence of absorbance, can be derived directly from the Lorentz-Lorenz relation (or, equivalently, the Clausius-Mossotti relation). Correspondingly, it can be demonstrated that there is a twin law according to which the change of the refractive index is approximately linear to the molar concentration for diluted solutions. This twin law can also be derived from the Lorentz-Lorenz relation.

Chemical analysis by spectrophotometry

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

$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

$\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.

Application for the atmosphere

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

$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} }+\ldots )},$

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.

See also

• Applied spectroscopy
• Atomic absorption spectroscopy
• Absorption spectroscopy
• Cavity ring-down spectroscopy
• Clausius-Mossotti relation
• Infra-red spectroscopy
• Job plot
• Laser absorption spectrometry
• Lorentz-Lorenz relation
• Logarithm
• Polymer degradation
• Scientific laws named after people
• Quantification of nucleic acids
• Tunable diode laser absorption spectroscopy

External links

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