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).

• 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)]