Hi everyone,

I'm trying to reproduce a plot of water defect modes in a one-dimensional photonic biosensor using MATLAB, but I'm not getting the sharp transmission peak that I expected based on a research paper. Instead, I'm seeing a peak of only around 0.1748, whereas the research paper's plot shows a sharp peak with a transmittance close to 1.

Background: The design I'm working on is for a photonic biosensor intended to detect several types of waterborne bacteria. The defect mode should show a sharp resonance peak when pure water is used as the defect layer in the photonic crystal (PhC) structure. I'm using the transfer matrix method to model the transmittance.

However, I seem to be missing something or making an error in my code that's preventing me from getting the correct results. I’ve tried adjusting parameters like the defect layer thickness and the design wavelength, but it still doesn’t match the expected result. I’m new to this, so any help would be appreciated!

I've attached both my MATLAB code and the original plot from the research paper. Could anyone please help me figure out what I might be doing wrong?