Brendel–Bormann oscillator model explained

The Brendel–Bormann oscillator model is a mathematical formula for the frequency dependence of the complex-valued relative permittivity, sometimes referred to as the dielectric function. The model has been used to fit to the complex refractive index of materials with absorption lineshapes exhibiting non-Lorentzian broadening, such as metals and amorphous insulators, across broad spectral ranges, typically near-ultraviolet, visible, and infrared frequencies. The dispersion relation bears the names of R. Brendel and D. Bormann, who derived the model in 1992, despite first being applied to optical constants in the literature by Andrei M. Efimov and E. G. Makarova in 1983. Around that time, several other researchers also independently discovered the model. The Brendel-Bormann oscillator model is aphysical because it does not satisfy the Kramers–Kronig relations. The model is non-causal, due to a singularity at zero frequency, and non-Hermitian. These drawbacks inspired J. Orosco and C. F. M. Coimbra to develop a similar, causal oscillator model.

Mathematical formulation

The general form of an oscillator model is given by

\varepsilon(\omega) = \varepsilon _{infty} + \sum _{j} \chi _{j} where\varepsilon is the relative permittivity,\varepsilon _{infty} is the value of the relative permittivity at infinite frequency,\omega is the angular frequency,\chi _{j} is the contribution from the j th absorption mechanism oscillator.\left(\chi ^{L} \right) and Gaussian oscillator \left(\chi ^{G} \right) , given by\chi ^{L} _{j} (\omega; \omega _{0,j} ) = (s _{j})/(\omega ^{2} _{0 ,j} - \omega ^{2} - i \Gamma _{j} \omega) \chi ^{G} _{j} (\omega) = (1)/(\sqrt{ 2 \pi) \sigma _{j} } \expwheres _{j} is the Lorentzian strength of the j th oscillator,\omega _{0,j} is the Lorentzian resonant frequency of the j th oscillator,\Gamma _{j} is the Lorentzian broadening of the j th oscillator,\sigma _{j} is the Gaussian broadening of the j th oscillator.

The Brendel-Bormann oscillator

\left(\chi ^{BB} \right) is obtained from the convolution of the two aforementioned oscillators in the manner of\chi ^{BB} _{j} (\omega) = \int ^{inf ty} _{-inf ty} \chi ^{G} _{j} (x - \omega _{0,j} ) \chi ^{L} _{j} (\omega; x) dx ,

which yields

\chi ^{BB} _{j} (\omega) = (i \sqrt{\pi)/(s) _{j} } \left[w\left(\frac{a_{j}(\omega) - \omega_{0,j}}{\sqrt{2}\sigma_{j}} \right) + w\left(\frac{a_{j}(\omega) + \omega_{0,j}}{\sqrt{2}\sigma_{j}} \right) \right]wherew(z) is the Faddeeva function,a _{j} = \sqrt{\omega ^{2} + i \Gamma _{j} \omega} .

The square root in the definition of

a _{j} must be taken such that its imaginary component is positive. This is achieved by:\Re\left( a _{j} \right) = \omega \sqrt{ (\sqrt{ 1+ \left( \Gamma _{j} /\omega \right) ^{2})/(+1 }{ 2 }}) \Im\left( a _{j} \right) = \omega \sqrt{ (\sqrt{ 1+ \left( \Gamma _{j} /\omega \right) ^{2})/(-1 }{ 2 }})

See also

• Cauchy equation
• Sellmeier equation
• Forouhi–Bloomer model
• Tauc–Lorentz model
• Lorentz oscillator model