The Forouhi–Bloomer model is a mathematical formula for the frequency dependence of the complex-valued refractive index. The model can be used to fit the refractive index of amorphous and crystalline semiconductor and dielectric materials at energies near and greater than their optical band gap. The dispersion relation bears the names of Rahim Forouhi and Iris Bloomer, who created the model and interpreted the physical significance of its parameters. The model is aphysical due to its incorrect asymptotic behavior and non-Hermitian character. These shortcomings inspired modified versions of the model as well as development of the Tauc–Lorentz model.
Mathematical formulation The complex refractive index is given by
n ~ ( E ) = n ( E ) + i κ ( E ) {\displaystyle {\tilde {n}}(E)=n(E)+i\kappa (E)}
where
n {\displaystyle n} is the real component of the complex refractive index, commonly called the refractive index,
κ {\displaystyle \kappa } is the imaginary component of the complex refractive index, commonly called the extinction coefficient,
E {\displaystyle E} is the photon energy (related to the angular frequency by E = ℏ ω {\displaystyle E=\hbar \omega } ). The real and imaginary components of the refractive index are related to one another through the Kramers-Kronig relations. Forouhi and Bloomer derived a formula for κ ( E ) {\displaystyle \kappa (E)} for amorphous materials. The formula and complementary Kramers–Kronig integral are given by
κ ( E ) = A ( E − E g ) 2 E 2 − B E + C {\displaystyle \kappa (E)={\frac {A\left(E-E_{g}\right)^{2}}{E^{2}-BE+C}}}
n ( E ) = n ∞ + 1 π P ∫ − ∞ ∞ κ ( ξ ) − κ ∞ ξ − E d ξ {\displaystyle n(E)=n_{\infty }+{\frac {1}{\pi }}{\mathcal {P}}\int _{-\infty }^{\infty }{\frac {\kappa (\xi )-\kappa _{\infty }}{\xi -E}}d\xi }
where
E g {\displaystyle E_{g}} is the bandgap of the material,
A {\displaystyle A} , B {\displaystyle B} , C {\displaystyle C} , and n ∞ {\displaystyle n_{\infty }} are fitting parameters,
P {\displaystyle {\mathcal {P}}} denotes the Cauchy principal value,
κ ∞ = lim E → ∞ κ ( E ) = A {\displaystyle \kappa _{\infty }=\lim _{E\rightarrow \infty }\kappa (E)=A} .
A {\displaystyle A} , B {\displaystyle B} , and C {\displaystyle C} are subject to the constraints A > 0 {\displaystyle A>0} , B > 0 {\displaystyle B>0} , C > 0 {\displaystyle C>0} , and 4 C − B 2 > 0 {\displaystyle 4C-B^{2}>0} . Evaluating the Kramers-Kronig integral,
n ( E ) = n ∞ + B 0 E + C 0 E 2 − B E + C {\displaystyle n(E)=n_{\infty }+{\frac {B_{0}E+C_{0}}{E^{2}-BE+C}}}
where
Q = 1 2 4 C − B 2 {\displaystyle Q={\frac {1}{2}}{\sqrt {4C-B^{2}}}} ,
B 0 = A Q ( − 1 2 B 2 + E g B − E g 2 + C ) {\displaystyle B_{0}={\frac {A}{Q}}\left(-{\frac {1}{2}}B^{2}+E_{g}B-E_{g}^{2}+C\right)} ,
C 0 = A Q ( 1 2 B ( E g 2 + C ) − 2 E g C ) {\displaystyle C_{0}={\frac {A}{Q}}\left({\frac {1}{2}}B\left(E_{g}^{2}+C\right)-2E_{g}C\right)} . The Forouhi–Bloomer model for crystalline materials is similar to that of amorphous materials. The formulas for n ( E ) {\displaystyle n(E)} and κ ( E ) {\displaystyle \kappa (E)} are given by
n ( E ) = n ∞ + ∑ j B 0 , j E + C 0 , j E 2 − B j E + C j {\displaystyle n(E)=n_{\infty }+\sum _{j}{\frac {B_{0,j}E+C_{0,j}}{E^{2}-B_{j}E+C_{j}}}} .
κ ( E ) = ( E − E g ) 2 ∑ j A j E 2 − B j E + C j {\displaystyle \kappa (E)=\left(E-E_{g}\right)^{2}\sum _{j}{\frac {A_{j}}{E^{2}-B_{j}E+C_{j}}}} . where all variables are defined similarly to the amorphous case, but with unique values for each value of the summation index j {\displaystyle j} . Thus, the model for amorphous materials is a special case of the model for crystalline materials when the sum is over a single term only.
References
See also Cauchy equation Sellmeier equation Lorentz oscillator model Tauc–Lorentz model Brendel–Bormann oscillator model