Home Forums FDTD Nonlinear Photonic Band Gap Simulation

Nonlinear Photonic Band Gap Simulation

Profile Photo
(joined November 2018)
Participant

Hello,

I am creating a nonlinear 3rd order photonic crystal out of silicon, with n2 value of 3*10^-18 through which I calculated the susceptibility using the equation n(2)=3Chi(3)/(4n^2). The dielectric constant I used was 15.813. The crystal has air holes with n=1. However, when I ran the PWE band solver simulation, there were no band gaps found, even though the contrast between the two refractive indices is great. Can you assist me with what might possibly yield band gaps?

I attached the file to this post.

Responses (1):

  • #50258
    Profile Photo
    Malik Javaid
    Participant

    Hello,

    I am creating a nonlinear 3rd order photonic crystal out of silicon, with n2 value of 3*10^-18 through which I calculated the susceptibility using the equation n(2)=3Chi(3)/(4n^2). The dielectric constant I used was 15.813. The crystal has air holes with n=1. However, when I ran the PWE band solver simulation, there were no band gaps found, even though the contrast between the two refractive indices is great. Can you assist me with what might possibly yield band gaps?

    I attached the file to this post.

  • #50269
    Profile Photo
    Scott Newman
    Participant

    Hello Malik,

    Run your PWEM calculation using the TM polarization and you will find that there is 2 bandgaps. The one I believe you are looking for is the one from [w/2pic] 0.517826 to 0.547916.

    I realize that the reference you are likely using probably says you should be looking at the TE polarization. It is important to note that there are different conventions for TE and TM. For a 2D simulation with infinite Y, OptiFDTD uses TE (Ey, Hx, Hz) and TM (Hy, Ex, Ez). While this may differ from some references it is inline with many waveguide based reference.

    On another note, the PWEM calculation in OptiFDTD does not take into consideration the power of the source so the bandgap you are calculating is based on the linear permitivity of 15.813 that you entered.

    Scott

Login You must be logged in to reply to this topic.