-
Notifications
You must be signed in to change notification settings - Fork 591
/
metasurface_lens_farfield.ctl
70 lines (50 loc) · 2.53 KB
/
metasurface_lens_farfield.ctl
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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
;; far fields at focal-length distance of a metalens with binary-phase grating unit cell from the Meep tutorial
(set-param! resolution 50) ; pixels/μm
(define-param dpml 1.0) ; PML thickness
(define-param dsub 2.0) ; substrate thickness
(define-param dpad 2.0) ; padding between grating and PML
(define-param gp 0.3) ; grating period
(define-param gh 1.8) ; grating height
(define-param focal-length 200) ; focal length of metalens
(define-param spot-length 100) ; far field line length
(define-param ff-res 10) ; far field resolution (points/μm)
; list of grating duty cycles
(define-param gdc-list (list '()))
; # of cells
(define num-cells (length gdc-list))
; return gdc of nth cell
(define gdc-cell (lambda (n) (list-ref gdc-list n)))
(define sx (+ dpml dsub gh dpad dpml))
(define sy (* num-cells gp))
(define cell (make lattice (size sx sy no-size)))
(set! geometry-lattice cell)
(define boundary-layers (list (make pml (thickness dpml) (direction X))))
(set! pml-layers boundary-layers)
(define-param lcen 0.5) ; center wavelength
(define fcen (/ lcen)) ; center frequency
(define df (* 0.2 fcen)) ; frequency width
(define pulse-src (list (make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez)
(center (+ (* -0.5 sx) dpml (* 0.5 dsub)) 0 0)
(size 0 sy 0))))
(set! sources pulse-src)
(set! k-point (vector3 0 0 0))
(define glass (make medium (index 1.5)))
(set! geometry (list (make block
(material glass)
(size (+ dpml dsub) infinity infinity)
(center (+ (* -0.5 sx) (* 0.5 (+ dpml dsub))) 0 0))))
(set! geometry (append geometry
(map (lambda (n)
(make block
(material glass)
(size gh (* (gdc-cell n) gp) infinity)
(center (+ (* -0.5 sx) dpml dsub (* 0.5 gh)) (+ (* -0.5 sy) (* (+ n 0.5) gp)) 0)))
(arith-sequence 0 1 num-cells))))
(define symm (list (make mirror-sym (direction Y))))
(set! symmetries symm)
(define mon-pt (vector3 (- (* 0.5 sx) dpml (* 0.5 dpad)) 0 0))
(define n2f-obj (add-near2far fcen 0 1 (make near2far-region (center mon-pt) (size 0 sy 0))))
(run-sources+ 500)
(output-farfields n2f-obj (string-append "numcells-" (number->string num-cells)) (volume (center focal-length 0 0) (size spot-length 0 0)) ff-res)