×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Jobs

Analysis of Optical Properties

Analysis of Optical Properties

Analysis of Optical Properties

(OP)
Hey all, I'm trying to create a method of analysis of the optical properties of layered materials, but I've encountered some difficulty.  We are using the Fresnel equation for Reflectance to calculate the n value, k value, and thickness (d) of a certain material on a substrate.  The formula is thus:

a = n0 - (n1 - ik1)
b = n0 + (n1 - ik1)
c = (n1 - ik1) - (n2 - ik2)
d = (n1 - ik1) + (n2 - ik2)
f = (2(pi)/y)*(n1 - ik1)*(d1)

R = [(a/b) + (c/d)*e^-2i(f)]/[1 + (a/b)*(c/d)*e^-2i(f)]

In this formula, the subscript (0) refers to air, (1) refers to the material being analyzed, and (2) is the substrate, silicon.  We can disregard any polarization in the equation, because we are using light that is normal to the surface.  (y) is wavelength in nanometers, and (i) is imaginary.  (e) is 2.181...

The goal of this analysis is to use measured values of reflectance to calculate n1, k1, and d1 simultaneously.  To generate a system of equations for solution, we assume the values of n1, k1, and d1 are very similar between different wavelengths, and let all other values (which are known) also vary with wavelength.  After tabulating the data for 850 wavelengths, I set up a Matlab program thus:

--Import the data.
--Assert that data fits specific parameters
Define relevant terms.
--Create a loop that analyzes the data 3 rows at a time (to create    the system of equations).
--Call a function file and apply loop to input argument.
--fsolve system of equations.
--End loop.
--Output results for n, k, and d in a new file.

The function file includes the above equation for reflectance, and creates equations with relevant data for the loop.

The problem I am having with this analysis is the resulting data.  None of the values of n, k, or d are appropriate.  As a test, we used materials in which all valued were known, including n1, k1, and d1.  My big question right now is this:  Is there a better way to feed Matlab this data?  Did we make a mistake somewhere?  Would you create a different method?

I appreciate all of your help in this matter, I've been working on the problem for a few weeks now with no improved results.  I am attaching the relevant file, which references both the program and the function file.

RE: Analysis of Optical Properties

Seems like you're trying to reinvent something that was successfully done more than 40 yrs ago, at least, for two layers of thin films, using polarized light, i.e., ellipsometry.

Nonetheless, I don't see how you could possibly get a plausible solution with made up data.  You asserted "we assumed the values of n1, k1, and d1 are very similar between different wavelengths."  So what possible differences could you be measuring?  The only plausible approach is that the different layers are NOT similar, and you are looking at the aggregate, spectrally variegated transmission.  Are you trying to measure some sort of Fabry-Perot effect?

But, it seems to me that there are a number of issues with the overall concept:
>  The measurement would be dominated by the first surface
>  The interaction of adjacent layers with different indices is a non trivial matter
>  The equation ignores the possibility of absorption by the materials
>  The equation seems to ignore the multiple reflections that would generate a Fabry-Perot behavior.

TTFN
FAQ731-376: Eng-Tips.com Forum Policies

RE: Analysis of Optical Properties

(OP)
The idea behind the project is to create a method of analysis for multi-layered thin film structures.  This is a test of data that is already known (not fabricated).  We are analyzing the light across 850 different wavelengths of light (not just a few), so across the data there are significant differences, but between 3 consecutive wavelengths the differences are very small.  We want to re-create a method of approximation, but we're running into the above difficulties.  The materials we are using for the analysis also fit the model above, so we can disregard absorption and Fabry-Perot effect.  I'm not looking for an analysis of my methods, I would like help with the actual Matlab code, which, regardless of how appropriate it's use is, should still work.

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members!


Resources