Drude-Lorentz and Debye-Lorentz models for the dielectric constant of metals and water

Computes the complex permittivity of many metals and water with input wavelength

Description

This code computes the complex dielectric constant (i.e. relative permittivity) and the refractive index of various metals using either the Lorentz-Drude (LD) or the Drude model (D) as a function of input light wavelength.

Additionally, it may compute the dielectric constant of pure water using a Debye-Lorentz model.

Choice of materials: silver, aluminum, gold, copper, chromium, nickel, tungsten, titanium, beryllium, palladium, platinum, triply distilled water.

Reference: Bora Ung and Yunlong Sheng, Interference of surface waves in a metallic nanoslit, Optics Express 15, 1182-1190 (2007).

• Required Products: MATLAB
• MATLAB release: MATLAB 7.4 (R2007a)

[Bora Ung]

The efficiency can be improved with the following changes:

replace line 183 with: epsilon_r_L = zeros(size(lambda));
replace lines 185-187 with:
epsilon_r_L = epsilon_r_L + (f(k)*omegap^2)*...
    (((omega(k)^2)*ones(size(lambda)) - omegalight.^2) -...
                i*Gamma(k)*omegalight).^(-1);
eliminate line 190

Hope this file is useful to people as it was to me.

[Duane Hanselman]

Rather than tell use what to do to improve the file, why not UPDATE your submission with your latest work? (Do not submit the revised file as a new submission.)

[Bora Ung]

Thanks Duane for the info. I have already done that a few days ago when I realized that it is possible to make changes. Now there IS NO NEED to make any changes.

[Jin Xu]

thanks a lot.

[Hyunwoo Kihm]

Thanks for valuable information with proper references. Weighting factor of various oscillator made the model fitting almost perfect.

[j kim]

I think, it's good program!!!

[JCBod Bodoque]

Wonderful file! Do you have any work of your own that I can cite for these fits?

[Bora Ung]

As a note, this matlab code was created for and first used in [B. Ung and Y. Sheng, Optics Express, vol.15, pp. 1182-1190 (2007)]

[Noha Anous]

Thanks for uploading such a nice function.
But I have a question, from where are the model equations taken?
are there any approximations?
Because I tried some values and didn't get them the same as I have them in tables found in some papers.
Thanks a lot

[Dominic Bosomtwi]

Can someone explain to me how to use the code?
Thanks a lot.

[chao ding]

Thanks! It's very useful for me.

[Rene]

Thank you, it is useful for me too.

[Johannes Srajer]