Dear Meep Users
1) I am trying to reproduce the results from the paper
Mateus et al IEEE PHOTONICS TECHNOLOGY LETTERS, V. 16, No. 2, pp. 518-520
Feb 2004. But, could not be able to do so. Please suggest modifications in
my approach. Please see attached .ctl file
2) Is there any way we can include 'gain' material and use it for lasing
phenomenan? Also, can it support saturation?
Thanks for time
Regards
Rakesh
(reset-meep)
(define-param lam 1.35) ; free space wavelength (microns)
; Some parameters to describe the geometry:
(define-param len .525) ; length of square
(define-param thick 0.46) ; thickness of square
(define-param gap .175) ; gap between the two squares
(define-param dpml 1) ; PML thickness
(define-param pad 1) ; the distance from pml
; The lattice dimensions
(define sx .7) ; size of cell in x direction
(define sy 6) ; size of cell in y direction
(set! geometry-lattice
(make lattice
(size sx sy no-size)
)
)
(set-param! resolution 50)
(define-param no-geometry? false) ; if true, no geometry
; Set Geometry
(set! geometry
(if no-geometry?
(list
(make block
(center 0 0)
(size len thick)
(material air)
)
)
(list
(make block
(center 0 0)
(size len thick)
(material (make dielectric (epsilon 12.816)))
)
(make block
(center 0 .68)
(size infinity .9)
(material (make dielectric (epsilon 2.1)))
)
(make block
(center 0 1.68)
(size infinity 1.22)
(material (make dielectric (epsilon 12.25)))
)
)
)
)
(set-param! k-point (vector3 0 0)) ; set periodic boundary condition
(set! pml-layers
(list
(make pml
(thickness dpml)
(direction Y)
)
)
)
; Set Source
(define fcen (/ 1 lam)) ; source center frequency
(define df .412) ; free space pulse width (in frequency)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen)
(fwidth df)))
(component Ex)
(center 0 (+ (/ sy -2) dpml pad))
(size sx 0)
)
)
)
(define-param nfreq 80) ; number of frequencies at which to compute flux
(define trans ; transmitted flux
(add-flux fcen df nfreq
(make flux-region
(center 0 1.68)
(size sx 0)
)
)
)
(define refl ; reflected flux
(add-flux fcen df nfreq
(make flux-region
(center 0 (+ (/ Sy -2) pad))
(size sx 0)
)
)
)
(run-sources+
(stop-when-fields-decayed
50 Ex
(vector3 0 1.68)
1e-3
)
(at-beginning output-epsilon)
)
(display-fluxes trans refl)
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
