File: bragg.ctl

package info (click to toggle)
mpb 1.12.0-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 4,040 kB
  • sloc: ansic: 13,389; javascript: 9,901; makefile: 214; lisp: 44; sh: 4
file content (25 lines) | stat: -rw-r--r-- 770 bytes parent folder | download | duplicates (5)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
; Compute the bands at the X point for a quarter-wave stack Bragg
; mirror (this is the point that defines the band gap edges).

; the high and low indices:
(define-param n-lo 1.0)
(define-param n-hi 3.0)

(define-param w-hi (/ n-lo (+ n-hi n-lo))) ; a quarter-wave stack

(set! geometry-lattice (make lattice (size 1 no-size no-size))) ; 1d cell
(set! default-material (make dielectric (index n-lo)))
(set! geometry 
      (list
       (make cylinder 
	 (material (make dielectric (index n-hi)))
	 (center 0 0 0) (axis 1 0 0) 
	 (radius infinity) (height w-hi))))

(define-param kx 0.5)
(set! k-points (list (vector3 kx 0 0)))

(set-param! resolution 32)
(set-param! num-bands 8)

(run-tm output-hfield-y) ; note that TM and TE bands are degenerate, so we only need TM