The quickest and simplest way to do this I think is either with finite-difference time-domain (FDTD) or finite-difference frequency-domain (FDFD). You can find notes and code snippets for both of these all over the internet. For FDTD, try here:

http://emlab.utep.edu/ee5390fdtd.htm

For FDFD, see Lectures 6 to 14 here:

http://emlab.utep.edu/ee5390cem.htm

There is also an older set of notes and codes for FDFD here under "Short Courses":

http://emlab.utep.edu/academics.htm

If I had to start from scratch and get the quickest possible answer, I would use FDFD. The materials above explain a 2D FDFD. You can use that as a 1D code just by setting the size of the grid in one of the directions to 1.

Hope this helps!

Hi Bhuvneshwer

A possible alternative to also consider as well as Raymond's suggestions is the transfer matrix method (TMM) and the related formalism of the scattering matrix method (SMM).

I've successfully implemented a TMM/SMM in MatLab which is able to calculate fields for finite 1D photonic crystals and other arbitrary slab stacks. It is possible to use the TMM/SMM with infinite PhCs which is often referred to as "rigourous coupled wave analysis" (RCWA) methods, or you could investigate using a large finite structure to approximate an infinite one. The advantages for me of this method (I have no experience of implementing any other method) were that I found it straightforward to follow and MatLab obviously is built to handle matrices well!

I'm making the assumption that you're refer to stacked slabs when you say 1D PhCs, rather than another configuration (e.g. line gratings or similar), which requires the RCWA side to the TMM/SMM approach.

There's lots of information available on the TMM/SMM which is easy to come by. To start, I found "Optical waves in layered media" by Pochi Yeh a very good text on the topic. There's a freely available piece of powerful software called CAMFR which implements this model in the context you're interested in and runs using a Python interface, available at http://camfr.sourceforge.net/ . Its development is described well in this thesis: http://photonics.intec.ugent.be/publications/PhD.asp?ID=104 , which I found extremely useful.

Best of luck!

Chris Lewis can you send me the matlab program using TMM at [email protected]

TMM is pretty straightforward to implement, but a bit more complicated if you wish to visualize the fields. The same goes for RCWA, method of lines, and other similar methods. I believe the CAMFR program is based on the method of lines. If you can't find codes for them, I cover the methods in detail here:

http://emlab.utep.edu/ee5390cem.htm

The block diagrams are as close to code as it gets.

Hi Bhuvneshwer, 

The algorithm to plot the fields in layered structures is described here

http://journals.aps.org/prb/abstract/10.1103/PhysRevB.90.045414

Good luck in field plotting

Hello Sergey A. Dyakov,

Thank you for your reply. But I don't have excess PRL's paper. So, kindly request to please send me this paper through email: [email protected] along with program or algorithms.

Can anyone send me the code of effective refractive calculation through Matlab ?

Thanks