Drude model in Transfer Matrix Method?

时间：03-30 整理：3721RD 点击：

Base on described theory about TMM for photonic crystal J.B. Pendry, J. Mod. Opt.41, 209 (1994), Prof. Pendry also uploaded a program using Fortran language [attachment] to calculate transmission/reflection coefficients of complex structure (2D dimensions).
The problem is, i want to perform this code to calculate T&R spectra adag_v1_0.gzfrom cut-wire pair structure.
This is my subroutine to define CWP structure, but it provide wrong results when i compare with CST program:

Code Fortran - [expand]
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
subroutine DEFSTR(EPSPOT,IXMAX,IYMAX,IZMAX,EPS,EPSREF,EPSSUB)
c
      integer IXMAX,IYMAX,IZMAX,IX,IY,IZ
      complex*16 EPSPOT(IXMAX,IYMAX,IZMAX),EPS,EPSREF,EPSSUB
 
      do 5 IX = 1,IXMAX
         do 5 IY = 1,IYMAX
            do 5 IZ = 1,IZMAX
               if ((IZ .EQ. 2) .OR. (IZ .EQ. 4)) then
                   if ((IY .GE. 3).AND.(IY .LE. 16) .AND.
     &                (IX .GE. 4) .AND. (IX .LE. 6)) then
                        EPSPOT(IX,IY,IZ) = EPS
                  else
                        EPSPOT(IX,IY,IZ) = EPSREF
                  end if
               elseif ((IZ .EQ. 3)) then
                   EPSPOT(IX,IY,IZ) = EPSSUB
               else
                        EPSPOT(IX,IY,IZ) = EPSREF
                end if
 5    continue
      return
      end

Code Fortran - [expand]
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23	subroutine DEFSTR(EPSPOT,IXMAX,IYMAX,IZMAX,EPS,EPSREF,EPSSUB) c integer IXMAX,IYMAX,IZMAX,IX,IY,IZ complex*16 EPSPOT(IXMAX,IYMAX,IZMAX),EPS,EPSREF,EPSSUB do 5 IX = 1,IXMAX do 5 IY = 1,IYMAX do 5 IZ = 1,IZMAX if ((IZ .EQ. 2) .OR. (IZ .EQ. 4)) then if ((IY .GE. 3).AND.(IY .LE. 16) .AND. & (IX .GE. 4) .AND. (IX .LE. 6)) then EPSPOT(IX,IY,IZ) = EPS else EPSPOT(IX,IY,IZ) = EPSREF end if elseif ((IZ .EQ. 3)) then EPSPOT(IX,IY,IZ) = EPSSUB else EPSPOT(IX,IY,IZ) = EPSREF end if 5 continue return end

and the second question is, when i try to use Drude's model for epsilon of metal part, then T and R become infinitive.
Anyone know this method, help me plz.

上一篇：Length and Width Parameter in LineCalc
下一篇：最后一页

model Drude Method 相关文章：

栏目分类