materials.scm

; Materials Library

; default unit length is 1 um
(define um-scale 1.0)

; conversion factor for eV to 1/um [=1/hc]
(define eV-um-scale (/ um-scale 1.23984193))

;------------------------------------------------------------------
; crystaline silicon (c-Si) from A. Deinega et al., J. Optical Society of America A, Vol. 28, No. 5, pp. 770-77, 2011
; based on experimental data for intrinsic silicon at T=300K from M.A. Green and M. Keevers, Progress in Photovoltaics, Vol. 3, pp. 189-92, 1995
; wavelength range: 0.4 - 1.0 um

(define cSi-frq1 (/ 3.64 um-scale))
(define cSi-gam1 0)
(define cSi-sig1 8)

(define cSi-frq2 (/ 2.76 um-scale))
(define cSi-gam2 (/ (* 2 0.063) um-scale))
(define cSi-sig2 2.85)

(define cSi-frq3 (/ 1.73 um-scale))
(define cSi-gam3 (/ (* 2 2.5) um-scale))
(define cSi-sig3 -0.107)

(define cSi (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency cSi-frq1) (gamma cSi-gam1) (sigma cSi-sig1))
  (make lorentzian-susceptibility
    (frequency cSi-frq2) (gamma cSi-gam2) (sigma cSi-sig2))
  (make lorentzian-susceptibility
    (frequency cSi-frq3) (gamma cSi-gam3) (sigma cSi-sig3)))))

;------------------------------------------------------------------
; amorphous silicon (a-Si) from Horiba Technical Note 08: Lorentz Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Lorentz_Dispersion_Model.pdf
; wavelength range: 0.21 - 0.83 um

(define aSi-frq1 (/ (* 0.315481407124682 um-scale)))
(define aSi-gam1 (/ (* 0.645751005208333 um-scale)))
(define aSi-sig1 14.571)

(define aSi (make medium (epsilon 3.109)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency aSi-frq1) (gamma aSi-gam1) (sigma aSi-sig1)))))

;------------------------------------------------------------------
; hydrogenated amorphous silicon (a-Si:H) from Horiba Technical Note 08: Lorentz Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Lorentz_Dispersion_Model.pdf
; wavelength range: 0.21 - 0.83 um

(define aSi-H-frq1 (/ (* 0.334189199460916 um-scale)))
(define aSi-H-gam1 (/ (* 0.579365387850467 um-scale)))
(define aSi-H-sig1 12.31)

(define aSi-H (make medium (epsilon 3.22)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency aSi-H-frq1) (gamma aSi-H-gam1) (sigma aSi-H-sig1)))))

;------------------------------------------------------------------
; indium tin oxide (ITO) from Horiba Technical Note 08: Lorentz Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Lorentz_Dispersion_Model.pdf
; wavelength range: 0.21 - 0.83 um

(define ITO-frq1 (/ (* 0.182329695588235 um-scale)))
(define ITO-gam1 (/ (* 1.94637665620094 um-scale)))
(define ITO-sig1 2.5)

(define ITO (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency ITO-frq1) (gamma ITO-gam1) (sigma ITO-sig1)))))

;------------------------------------------------------------------
; alumina (Al2O3) from Horiba Technical Note 08: Lorentz Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Lorentz_Dispersion_Model.pdf
; wavelength range: 0.21 - 2.07 um

(define Al2O3-frq1 (/ (* 0.101476668030774 um-scale)))
(define Al2O3-gam1 0)
(define Al2O3-sig1 1.52)

(define Al2O3 (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency Al2O3-frq1) (gamma Al2O3-gam1) (sigma Al2O3-sig1)))))

;------------------------------------------------------------------
; aluminum nitride (AlN) from Horiba Technical Note 08: Lorentz Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Lorentz_Dispersion_Model.pdf
; wavelength range: 0.26 - 1.65 um

(define AlN-frq1 (/ (* 0.139058089950651 um-scale)))
(define AlN-gam1 0)
(define AlN-sig1 3.306)

(define AlN (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency AlN-frq1) (gamma AlN-gam1) (sigma AlN-sig1)))))

;------------------------------------------------------------------
; aluminum arsenide (AlAs) from R.E. Fern and A. Onton, J. Applied Physics, Vol. 42, pp. 3499-500, 1971
; ref: https://refractiveindex.info/?shelf=main&book=AlAs&page=Fern
; wavelength range: 0.56 - 2.2 um

(define AlAs-frq1 (/ (* 0.2822 um-scale)))
(define AlAs-gam1 0)
(define AlAs-sig1 6.0840)

(define AlAs-frq2 (/ (* 27.62 um-scale)))
(define AlAs-gam2 0)
(define AlAs-sig2 1.900)

(define AlAs (make medium (epsilon 2.0792)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency AlAs-frq1) (gamma AlAs-gam1) (sigma AlAs-sig1))
  (make lorentzian-susceptibility
    (frequency AlAs-frq2) (gamma AlAs-gam2) (sigma AlAs-sig2)))))

;------------------------------------------------------------------
; borosilicate glass (BK7) from SCHOTT Zemax catalog 2017-01-20b
; ref: https://refractiveindex.info/?shelf=glass&book=BK7&page=SCHOTT
; wavelength range: 0.3 - 2.5 um

(define BK7-frq1 (/ (* 0.07746417668832478 um-scale)))
(define BK7-gam1 0)
(define BK7-sig1 1.03961212)
(define BK7-frq2 (/ (* 0.14148467902921502 um-scale)))
(define BK7-gam2 0)
(define BK7-sig2 0.231792344)
(define BK7-frq3 (/ (* 10.176475470417055 um-scale)))
(define BK7-gam3 0)
(define BK7-sig3 1.01046945)

(define BK7 (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency BK7-frq1) (gamma BK7-gam1) (sigma BK7-sig1))
  (make lorentzian-susceptibility
    (frequency BK7-frq2) (gamma BK7-gam2) (sigma BK7-sig2))
  (make lorentzian-susceptibility
    (frequency BK7-frq3) (gamma BK7-gam3) (sigma BK7-sig3)))))

;------------------------------------------------------------------
; fused quartz (silica) from I.H. Malitson, J. Optical Society of America, Vol. 55, pp. 1205-9, 1965
; ref: https://refractiveindex.info/?shelf=glass&book=fused_silica&page=Malitson
; wavelength range: 0.21 - 6.7 um

(define fused-quartz-frq1 (/ (* 0.0684043 um-scale)))
(define fused-quartz-gam1 0)
(define fused-quartz-sig1 0.696166300)
(define fused-quartz-frq2 (/ (* 0.1162414 um-scale)))
(define fused-quartz-gam2 0)
(define fused-quartz-sig2 0.407942600)
(define fused-quartz-frq3 (/ (* 9.896161 um-scale)))
(define fused-quartz-gam3 0)
(define fused-quartz-sig3 0.897479400)

(define fused-quartz (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency fused-quartz-frq1) (gamma fused-quartz-gam1) (sigma fused-quartz-sig1))
  (make lorentzian-susceptibility
    (frequency fused-quartz-frq2) (gamma fused-quartz-gam2) (sigma fused-quartz-sig2))
  (make lorentzian-susceptibility
    (frequency fused-quartz-frq3) (gamma fused-quartz-gam3) (sigma fused-quartz-sig3)))))

;------------------------------------------------------------------
; gallium arsenide (GaAs) from T. Skauli et al., J. Applied Physics, Vol. 94, pp. 6447-55, 2003
; ref: https://refractiveindex.info/?shelf=main&book=GaAs&page=Skauli
; wavelength range: 0.97 - 17 um

(define GaAs-frq1 (/ (* 0.4431307 um-scale)))
(define GaAs-gam1 0)
(define GaAs-sig1 5.466742)
(define GaAs-frq2 (/ (* 0.8746453 um-scale)))
(define GaAs-gam2 0)
(define GaAs-sig2 0.02429960)
(define GaAs-frq3 (/ (* 36.9166 um-scale)))
(define GaAs-gam3 0)
(define GaAs-sig3 1.957522)

(define GaAs (make medium (epsilon 5.372514)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency GaAs-frq1) (gamma GaAs-gam1) (sigma GaAs-sig1))
  (make lorentzian-susceptibility
    (frequency GaAs-frq2) (gamma GaAs-gam2) (sigma GaAs-sig2))
  (make lorentzian-susceptibility
    (frequency GaAs-frq3) (gamma GaAs-gam3) (sigma GaAs-sig3)))))

;------------------------------------------------------------------
; silicon nitride (Si3N4) from H. R. Philipp, J. Electrochemical Society, Vol. 120, pp. 295-300, 1973
; ref: https://refractiveindex.info/?shelf=main&book=Si3N4&page=Philipp
; wavelength range: 0.207 - 1.24 um

(define Si3N4-VISNIR-frq1 (/ (* 0.13967 um-scale)))
(define Si3N4-VISNIR-gam1 0)
(define Si3N4-VISNIR-sig1 2.8939)

(define Si3N4-VISNIR (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency Si3N4-VISNIR-frq1) (gamma Si3N4-VISNIR-gam1) (sigma Si3N4-VISNIR-sig1)))))

;------------------------------------------------------------------
; silicon nitride (Si3N4) from K. Luke, et. al., Optics Letters, Vol. 40, pp. 4823-4826, 2015
; ref: https://refractiveindex.info/?shelf=main&book=Si3N4&page=Luke
; wavelength range: 0.310 - 5.504 um

(define Si3N4-NIR-frq1 (/ (* 0.1353406 um-scale)))
(define Si3N4-NIR-gam1 0)
(define Si3N4-NIR-sig1 3.0249)
(define Si3N4-NIR-frq2 (/ (* 1239.842 um-scale)))
(define Si3N4-NIR-gam2 0)
(define Si3N4-NIR-sig2 40314)

(define Si3N4-NIR (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency Si3N4-NIR-frq1) (gamma Si3N4-NIR-gam1) (sigma Si3N4-NIR-sig1))
  (make lorentzian-susceptibility
    (frequency Si3N4-NIR-frq2) (gamma Si3N4-NIR-gam2) (sigma Si3N4-NIR-sig2)))))

;------------------------------------------------------------------
; elemental metals from A.D. Rakic et al., Applied Optics, Vol. 37, No. 22, pp. 5271-83, 1998
; wavelength range: 0.2 - 12.4 um

; silver (Ag)

(define Ag-plasma-frq (* 9.01 eV-um-scale))

(define Ag-f0 0.845)
(define Ag-frq0 1e-10)
(define Ag-gam0 (* 0.048 eV-um-scale))
(define Ag-sig0 (/ (* Ag-f0 (sqr Ag-plasma-frq)) (sqr Ag-frq0)))

(define Ag-f1 0.065)
(define Ag-frq1 (* 0.816 eV-um-scale)) ; 1.519 um
(define Ag-gam1 (* 3.886 eV-um-scale))
(define Ag-sig1 (/ (* Ag-f1 (sqr Ag-plasma-frq)) (sqr Ag-frq1)))

(define Ag-f2 0.124)
(define Ag-frq2 (* 4.481 eV-um-scale)) ; 0.273 um
(define Ag-gam2 (* 0.452 eV-um-scale))
(define Ag-sig2 (/ (* Ag-f2 (sqr Ag-plasma-frq)) (sqr Ag-frq2)))

(define Ag-f3 0.011)
(define Ag-frq3 (* 8.185 eV-um-scale)) ; 0.152 um
(define Ag-gam3 (* 0.065 eV-um-scale))
(define Ag-sig3 (/ (* Ag-f3 (sqr Ag-plasma-frq)) (sqr Ag-frq3)))

(define Ag-f4 0.840)
(define Ag-frq4 (* 9.083 eV-um-scale)) ; 0.137 um
(define Ag-gam4 (* 0.916 eV-um-scale))
(define Ag-sig4 (/ (* Ag-f4 (sqr Ag-plasma-frq)) (sqr Ag-frq4)))

(define Ag-f5 5.646)
(define Ag-frq5 (* 20.29 eV-um-scale)) ; 0.061 um
(define Ag-gam5 (* 2.419 eV-um-scale))
(define Ag-sig5 (/ (* Ag-f5 (sqr Ag-plasma-frq)) (sqr Ag-frq5)))

(define Ag (make medium (epsilon 1.0)
  (E-susceptibilities
     (make drude-susceptibility
       (frequency Ag-frq0) (gamma Ag-gam0) (sigma Ag-sig0))
     (make lorentzian-susceptibility
       (frequency Ag-frq1) (gamma Ag-gam1) (sigma Ag-sig1))
     (make lorentzian-susceptibility
       (frequency Ag-frq2) (gamma Ag-gam2) (sigma Ag-sig2))
     (make lorentzian-susceptibility
       (frequency Ag-frq3) (gamma Ag-gam3) (sigma Ag-sig3))
     (make lorentzian-susceptibility
       (frequency Ag-frq4) (gamma Ag-gam4) (sigma Ag-sig4))
     (make lorentzian-susceptibility
       (frequency Ag-frq5) (gamma Ag-gam5) (sigma Ag-sig5)))))

;------------------------------------------------------------------
; gold (Au)

(define Au-plasma-frq (* 9.03 eV-um-scale))

(define Au-f0 0.760)
(define Au-frq0 1e-10)
(define Au-gam0 (* 0.053 eV-um-scale))
(define Au-sig0 (/ (* Au-f0 (sqr Au-plasma-frq)) (sqr Au-frq0)))

(define Au-f1 0.024)
(define Au-frq1 (* 0.415 eV-um-scale)) ; 2.988 um
(define Au-gam1 (* 0.241 eV-um-scale))
(define Au-sig1 (/ (* Au-f1 (sqr Au-plasma-frq)) (sqr Au-frq1)))

(define Au-f2 0.010)
(define Au-frq2 (* 0.830 eV-um-scale)) ; 1.494 um
(define Au-gam2 (* 0.345 eV-um-scale))
(define Au-sig2 (/ (* Au-f2 (sqr Au-plasma-frq)) (sqr Au-frq2)))

(define Au-f3 0.071)
(define Au-frq3 (* 2.969 eV-um-scale)) ; 0.418 um
(define Au-gam3 (* 0.870 eV-um-scale))
(define Au-sig3 (/ (* Au-f3 (sqr Au-plasma-frq)) (sqr Au-frq3)))

(define Au-f4 0.601)
(define Au-frq4 (* 4.304 eV-um-scale)) ; 0.288 um
(define Au-gam4 (* 2.494 eV-um-scale))
(define Au-sig4 (/ (* Au-f4 (sqr Au-plasma-frq)) (sqr Au-frq4)))

(define Au-f5 4.384)
(define Au-frq5 (* 13.32 eV-um-scale)) ; 0.093 um
(define Au-gam5 (* 2.214 eV-um-scale))
(define Au-sig5 (/ (* Au-f5 (sqr Au-plasma-frq)) (sqr Au-frq5)))

(define Au (make medium (epsilon 1.0)
  (E-susceptibilities
     (make drude-susceptibility
       (frequency Au-frq0) (gamma Au-gam0) (sigma Au-sig0))
     (make lorentzian-susceptibility
       (frequency Au-frq1) (gamma Au-gam1) (sigma Au-sig1))
     (make lorentzian-susceptibility
       (frequency Au-frq2) (gamma Au-gam2) (sigma Au-sig2))
     (make lorentzian-susceptibility
       (frequency Au-frq3) (gamma Au-gam3) (sigma Au-sig3))
     (make lorentzian-susceptibility
       (frequency Au-frq4) (gamma Au-gam4) (sigma Au-sig4))
     (make lorentzian-susceptibility
       (frequency Au-frq5) (gamma Au-gam5) (sigma Au-sig5)))))

;------------------------------------------------------------------
; copper (Cu)

(define Cu-plasma-frq (* 10.83 eV-um-scale))

(define Cu-f0 0.575)
(define Cu-frq0 1e-10)
(define Cu-gam0 (* 0.030 eV-um-scale))
(define Cu-sig0 (/ (* Cu-f0 (sqr Cu-plasma-frq)) (sqr Cu-frq0)))

(define Cu-f1 0.061)
(define Cu-frq1 (* 0.291 eV-um-scale)) ; 4.261 um
(define Cu-gam1 (* 0.378 eV-um-scale))
(define Cu-sig1 (/ (* Cu-f1 (sqr Cu-plasma-frq)) (sqr Cu-frq1)))

(define Cu-f2 0.104)
(define Cu-frq2 (* 2.957 eV-um-scale)) ; 0.419 um
(define Cu-gam2 (* 1.056 eV-um-scale))
(define Cu-sig2 (/ (* Cu-f2 (sqr Cu-plasma-frq)) (sqr Cu-frq2)))

(define Cu-f3 0.723)
(define Cu-frq3 (* 5.300 eV-um-scale)) ; 0.234 um
(define Cu-gam3 (* 3.213 eV-um-scale))
(define Cu-sig3 (/ (* Cu-f3 (sqr Cu-plasma-frq)) (sqr Cu-frq3)))

(define Cu-f4 0.638)
(define Cu-frq4 (* 11.18 eV-um-scale)) ; 0.111 um
(define Cu-gam4 (* 4.305 eV-um-scale))
(define Cu-sig4 (/ (* Cu-f4 (sqr Cu-plasma-frq)) (sqr Cu-frq4)))

(define Cu (make medium (epsilon 1.0)
  (E-susceptibilities
     (make drude-susceptibility
       (frequency Cu-frq0) (gamma Cu-gam0) (sigma Cu-sig0))
     (make lorentzian-susceptibility
       (frequency Cu-frq1) (gamma Cu-gam1) (sigma Cu-sig1))
     (make lorentzian-susceptibility
       (frequency Cu-frq2) (gamma Cu-gam2) (sigma Cu-sig2))
     (make lorentzian-susceptibility
       (frequency Cu-frq3) (gamma Cu-gam3) (sigma Cu-sig3))
     (make lorentzian-susceptibility
       (frequency Cu-frq4) (gamma Cu-gam4) (sigma Cu-sig4)))))

;------------------------------------------------------------------
; aluminum (Al)

(define Al-plasma-frq (* 14.98 eV-um-scale))

(define Al-f0 0.523)
(define Al-frq0 1e-10)
(define Al-gam0 (* 0.047 eV-um-scale))
(define Al-sig0 (/ (* Al-f0 (sqr Al-plasma-frq)) (sqr Al-frq0)))

(define Al-f1 0.227)
(define Al-frq1 (* 0.162 eV-um-scale)) ; 7.654 um
(define Al-gam1 (* 0.333 eV-um-scale))
(define Al-sig1 (/ (* Al-f1 (sqr Al-plasma-frq)) (sqr Al-frq1)))

(define Al-f2 0.050)
(define Al-frq2 (* 1.544 eV-um-scale)) ; 0.803 um
(define Al-gam2 (* 0.312 eV-um-scale))
(define Al-sig2 (/ (* Al-f2 (sqr Al-plasma-frq)) (sqr Al-frq2)))

(define Al-f3 0.166)
(define Al-frq3 (* 1.808 eV-um-scale)) ; 0.686 um
(define Al-gam3 (* 1.351 eV-um-scale))
(define Al-sig3 (/ (* Al-f3 (sqr Al-plasma-frq)) (sqr Al-frq3)))

(define Al-f4 0.030)
(define Al-frq4 (* 3.473 eV-um-scale)) ; 0.357 um
(define Al-gam4 (* 3.382 eV-um-scale))
(define Al-sig4 (/ (* Al-f4 (sqr Al-plasma-frq)) (sqr Al-frq4)))

(define Al (make medium (epsilon 1.0)
  (E-susceptibilities
     (make drude-susceptibility
       (frequency Al-frq0) (gamma Al-gam0) (sigma Al-sig0))
     (make lorentzian-susceptibility
       (frequency Al-frq1) (gamma Al-gam1) (sigma Al-sig1))
     (make lorentzian-susceptibility
       (frequency Al-frq2) (gamma Al-gam2) (sigma Al-sig2))
     (make lorentzian-susceptibility
       (frequency Al-frq3) (gamma Al-gam3) (sigma Al-sig3))
     (make lorentzian-susceptibility
       (frequency Al-frq4) (gamma Al-gam4) (sigma Al-sig4)))))

;------------------------------------------------------------------
; beryllium (Be)

(define Be-plasma-frq (* 18.51 eV-um-scale))

(define Be-f0 0.084)
(define Be-frq0 1e-10)
(define Be-gam0 (* 0.035 eV-um-scale))
(define Be-sig0 (/ (* Be-f0 (sqr Be-plasma-frq)) (sqr Be-frq0)))

(define Be-f1 0.031)
(define Be-frq1 (* 0.100 eV-um-scale)) ; 12.398 um
(define Be-gam1 (* 1.664 eV-um-scale))
(define Be-sig1 (/ (* Be-f1 (sqr Be-plasma-frq)) (sqr Be-frq1)))

(define Be-f2 0.140)
(define Be-frq2 (* 1.032 eV-um-scale)) ; 1.201 um
(define Be-gam2 (* 3.395 eV-um-scale))
(define Be-sig2 (/ (* Be-f2 (sqr Be-plasma-frq)) (sqr Be-frq2)))

(define Be-f3 0.530)
(define Be-frq3 (* 3.183 eV-um-scale)) ; 0.390 um
(define Be-gam3 (* 4.454 eV-um-scale))
(define Be-sig3 (/ (* Be-f3 (sqr Be-plasma-frq)) (sqr Be-frq3)))

(define Be-f4 0.130)
(define Be-frq4 (* 4.604 eV-um-scale)) ; 0.269 um
(define Be-gam4 (* 1.802 eV-um-scale))
(define Be-sig4 (/ (* Be-f4 (sqr Be-plasma-frq)) (sqr Be-frq4)))

(define Be (make medium (epsilon 1.0)
  (E-susceptibilities
     (make drude-susceptibility
       (frequency Be-frq0) (gamma Be-gam0) (sigma Be-sig0))
     (make lorentzian-susceptibility
       (frequency Be-frq1) (gamma Be-gam1) (sigma Be-sig1))
     (make lorentzian-susceptibility
       (frequency Be-frq2) (gamma Be-gam2) (sigma Be-sig2))
     (make lorentzian-susceptibility
       (frequency Be-frq3) (gamma Be-gam3) (sigma Be-sig3))
     (make lorentzian-susceptibility
       (frequency Be-frq4) (gamma Be-gam4) (sigma Be-sig4)))))

;------------------------------------------------------------------
; chromium (Cr)

(define Cr-plasma-frq (* 10.75 eV-um-scale))

(define Cr-f0 0.168)
(define Cr-frq0 1e-10)
(define Cr-gam0 (* 0.047 eV-um-scale))
(define Cr-sig0 (/ (* Cr-f0 (sqr Cr-plasma-frq)) (sqr Cr-frq0)))

(define Cr-f1 0.151)
(define Cr-frq1 (* 0.121 eV-um-scale)) ; 10.247 um
(define Cr-gam1 (* 3.175 eV-um-scale))
(define Cr-sig1 (/ (* Cr-f1 (sqr Cr-plasma-frq)) (sqr Cr-frq1)))

(define Cr-f2 0.150)
(define Cr-frq2 (* 0.543 eV-um-scale)) ; 2.283 um
(define Cr-gam2 (* 1.305 eV-um-scale))
(define Cr-sig2 (/ (* Cr-f2 (sqr Cr-plasma-frq)) (sqr Cr-frq2)))

(define Cr-f3 1.149)
(define Cr-frq3 (* 1.970 eV-um-scale)) ; 0.629 um
(define Cr-gam3 (* 2.676 eV-um-scale))
(define Cr-sig3 (/ (* Cr-f3 (sqr Cr-plasma-frq)) (sqr Cr-frq3)))

(define Cr-f4 0.825)
(define Cr-frq4 (* 8.775 eV-um-scale)) ; 0.141 um
(define Cr-gam4 (* 1.335 eV-um-scale))
(define Cr-sig4 (/ (* Cr-f4 (sqr Cr-plasma-frq)) (sqr Cr-frq4)))

(define Cr (make medium (epsilon 1.0)
  (E-susceptibilities
     (make drude-susceptibility
       (frequency Cr-frq0) (gamma Cr-gam0) (sigma Cr-sig0))
     (make lorentzian-susceptibility
       (frequency Cr-frq1) (gamma Cr-gam1) (sigma Cr-sig1))
     (make lorentzian-susceptibility
       (frequency Cr-frq2) (gamma Cr-gam2) (sigma Cr-sig2))
     (make lorentzian-susceptibility
       (frequency Cr-frq3) (gamma Cr-gam3) (sigma Cr-sig3))
     (make lorentzian-susceptibility
       (frequency Cr-frq4) (gamma Cr-gam4) (sigma Cr-sig4)))))

;------------------------------------------------------------------
; nickel (Ni)

(define Ni-plasma-frq (* 15.92 eV-um-scale))

(define Ni-f0 0.096)
(define Ni-frq0 1e-10)
(define Ni-gam0 (* 0.048 eV-um-scale))
(define Ni-sig0 (/ (* Ni-f0 (sqr Ni-plasma-frq)) (sqr Ni-frq0)))

(define Ni-f1 0.100)
(define Ni-frq1 (* 0.174 eV-um-scale)) ; 7.126 um
(define Ni-gam1 (* 4.511 eV-um-scale))
(define Ni-sig1 (/ (* Ni-f1 (sqr Ni-plasma-frq)) (sqr Ni-frq1)))

(define Ni-f2 0.135)
(define Ni-frq2 (* 0.582 eV-um-scale)) ; 2.130 um
(define Ni-gam2 (* 1.334 eV-um-scale))
(define Ni-sig2 (/ (* Ni-f2 (sqr Ni-plasma-frq)) (sqr Ni-frq2)))

(define Ni-f3 0.106)
(define Ni-frq3 (* 1.597 eV-um-scale)) ; 0.776 um
(define Ni-gam3 (* 2.178 eV-um-scale))
(define Ni-sig3 (/ (* Ni-f3 (sqr Ni-plasma-frq)) (sqr Ni-frq3)))

(define Ni-f4 0.729)
(define Ni-frq4 (* 6.089 eV-um-scale)) ; 0.204 um
(define Ni-gam4 (* 6.292 eV-um-scale))
(define Ni-sig4 (/ (* Ni-f4 (sqr Ni-plasma-frq)) (sqr Ni-frq4)))

(define Ni (make medium (epsilon 1.0)
  (E-susceptibilities
     (make drude-susceptibility
       (frequency Ni-frq0) (gamma Ni-gam0) (sigma Ni-sig0))
     (make lorentzian-susceptibility
       (frequency Ni-frq1) (gamma Ni-gam1) (sigma Ni-sig1))
     (make lorentzian-susceptibility
       (frequency Ni-frq2) (gamma Ni-gam2) (sigma Ni-sig2))
     (make lorentzian-susceptibility
       (frequency Ni-frq3) (gamma Ni-gam3) (sigma Ni-sig3))
     (make lorentzian-susceptibility
       (frequency Ni-frq4) (gamma Ni-gam4) (sigma Ni-sig4)))))

;------------------------------------------------------------------
; palladium (Pd)

(define Pd-plasma-frq (* 9.72 eV-um-scale))

(define Pd-f0 0.330)
(define Pd-frq0 1e-10)
(define Pd-gam0 (* 0.008 eV-um-scale))
(define Pd-sig0 (/ (* Pd-f0 (sqr Pd-plasma-frq)) (sqr Pd-frq0)))

(define Pd-f1 0.649)
(define Pd-frq1 (* 0.336 eV-um-scale)) ; 3.690 um
(define Pd-gam1 (* 2.950 eV-um-scale))
(define Pd-sig1 (/ (* Pd-f1 (sqr Pd-plasma-frq)) (sqr Pd-frq1)))

(define Pd-f2 0.121)
(define Pd-frq2 (* 0.501 eV-um-scale)) ; 2.475 um
(define Pd-gam2 (* 0.555 eV-um-scale))
(define Pd-sig2 (/ (* Pd-f2 (sqr Pd-plasma-frq)) (sqr Pd-frq2)))

(define Pd-f3 0.638)
(define Pd-frq3 (* 1.659 eV-um-scale)) ; 0.747 um
(define Pd-gam3 (* 4.621 eV-um-scale))
(define Pd-sig3 (/ (* Pd-f3 (sqr Pd-plasma-frq)) (sqr Pd-frq3)))

(define Pd-f4 0.453)
(define Pd-frq4 (* 5.715 eV-um-scale)) ; 0.217 um
(define Pd-gam4 (* 3.236 eV-um-scale))
(define Pd-sig4 (/ (* Pd-f4 (sqr Pd-plasma-frq)) (sqr Pd-frq4)))

(define Pd (make medium (epsilon 1.0)
  (E-susceptibilities
     (make drude-susceptibility
       (frequency Pd-frq0) (gamma Pd-gam0) (sigma Pd-sig0))
     (make lorentzian-susceptibility
       (frequency Pd-frq1) (gamma Pd-gam1) (sigma Pd-sig1))
     (make lorentzian-susceptibility
       (frequency Pd-frq2) (gamma Pd-gam2) (sigma Pd-sig2))
     (make lorentzian-susceptibility
       (frequency Pd-frq3) (gamma Pd-gam3) (sigma Pd-sig3))
     (make lorentzian-susceptibility
       (frequency Pd-frq4) (gamma Pd-gam4) (sigma Pd-sig4)))))

;------------------------------------------------------------------
; platinum (Pt)

(define Pt-plasma-frq (* 9.59 eV-um-scale))

(define Pt-f0 0.333)
(define Pt-frq0 1e-10)
(define Pt-gam0 (* 0.080 eV-um-scale))
(define Pt-sig0 (/ (* Pt-f0 (sqr Pt-plasma-frq)) (sqr Pt-frq0)))

(define Pt-f1 0.191)
(define Pt-frq1 (* 0.780 eV-um-scale)) ; 1.590 um
(define Pt-gam1 (* 0.517 eV-um-scale))
(define Pt-sig1 (/ (* Pt-f1 (sqr Pt-plasma-frq)) (sqr Pt-frq1)))

(define Pt-f2 0.659)
(define Pt-frq2 (* 1.314 eV-um-scale)) ; 0.944 um
(define Pt-gam2 (* 1.838 eV-um-scale))
(define Pt-sig2 (/ (* Pt-f2 (sqr Pt-plasma-frq)) (sqr Pt-frq2)))

(define Pt-f3 0.547)
(define Pt-frq3 (* 3.141 eV-um-scale)) ; 0.395 um
(define Pt-gam3 (* 3.668 eV-um-scale))
(define Pt-sig3 (/ (* Pt-f3 (sqr Pt-plasma-frq)) (sqr Pt-frq3)))

(define Pt-f4 3.576)
(define Pt-frq4 (* 9.249 eV-um-scale)) ; 0.134 um
(define Pt-gam4 (* 8.517 eV-um-scale))
(define Pt-sig4 (/ (* Pt-f4 (sqr Pt-plasma-frq)) (sqr Pt-frq4)))

(define Pt (make medium (epsilon 1.0)
  (E-susceptibilities
     (make drude-susceptibility
       (frequency Pt-frq0) (gamma Pt-gam0) (sigma Pt-sig0))
     (make lorentzian-susceptibility
       (frequency Pt-frq1) (gamma Pt-gam1) (sigma Pt-sig1))
     (make lorentzian-susceptibility
       (frequency Pt-frq2) (gamma Pt-gam2) (sigma Pt-sig2))
     (make lorentzian-susceptibility
       (frequency Pt-frq3) (gamma Pt-gam3) (sigma Pt-sig3))
     (make lorentzian-susceptibility
       (frequency Pt-frq4) (gamma Pt-gam4) (sigma Pt-sig4)))))

;------------------------------------------------------------------
; titanium (Ti)

(define Ti-plasma-frq (* 7.29 eV-um-scale))

(define Ti-f0 0.148)
(define Ti-frq0 1e-10)
(define Ti-gam0 (* 0.082 eV-um-scale))
(define Ti-sig0 (/ (* Ti-f0 (sqr Ti-plasma-frq)) (sqr Ti-frq0)))

(define Ti-f1 0.899)
(define Ti-frq1 (* 0.777 eV-um-scale)) ; 1.596 um
(define Ti-gam1 (* 2.276 eV-um-scale))
(define Ti-sig1 (/ (* Ti-f1 (sqr Ti-plasma-frq)) (sqr Ti-frq1)))

(define Ti-f2 0.393)
(define Ti-frq2 (* 1.545 eV-um-scale)) ; 0.802 um
(define Ti-gam2 (* 2.518 eV-um-scale))
(define Ti-sig2 (/ (* Ti-f2 (sqr Ti-plasma-frq)) (sqr Ti-frq2)))

(define Ti-f3 0.187)
(define Ti-frq3 (* 2.509 eV-um-scale)) ; 0.494 um
(define Ti-gam3 (* 1.663 eV-um-scale))
(define Ti-sig3 (/ (* Ti-f3 (sqr Ti-plasma-frq)) (sqr Ti-frq3)))

(define Ti-f4 0.001)
(define Ti-frq4 (* 19.43 eV-um-scale)) ; 0.064 um
(define Ti-gam4 (* 1.762 eV-um-scale))
(define Ti-sig4 (/ (* Ti-f4 (sqr Ti-plasma-frq)) (sqr Ti-frq4)))

(define Ti (make medium (epsilon 1.0)
  (E-susceptibilities
     (make drude-susceptibility
       (frequency Ti-frq0) (gamma Ti-gam0) (sigma Ti-sig0))
     (make lorentzian-susceptibility
       (frequency Ti-frq1) (gamma Ti-gam1) (sigma Ti-sig1))
     (make lorentzian-susceptibility
       (frequency Ti-frq2) (gamma Ti-gam2) (sigma Ti-sig2))
     (make lorentzian-susceptibility
       (frequency Ti-frq3) (gamma Ti-gam3) (sigma Ti-sig3))
     (make lorentzian-susceptibility
       (frequency Ti-frq4) (gamma Ti-gam4) (sigma Ti-sig4)))))

;------------------------------------------------------------------
; tungsten (W)

(define W-plasma-frq (* 13.22 eV-um-scale))

(define W-f0 0.206)
(define W-frq0 1e-10)
(define W-gam0 (* 0.064 eV-um-scale))
(define W-sig0 (/ (* W-f0 (sqr W-plasma-frq)) (sqr W-frq0)))

(define W-f1 0.054)
(define W-frq1 (* 1.004 eV-um-scale)) ; 1.235 um
(define W-gam1 (* 0.530 eV-um-scale))
(define W-sig1 (/ (* W-f1 (sqr W-plasma-frq)) (sqr W-frq1)))

(define W-f2 0.166)
(define W-frq2 (* 1.917 eV-um-scale)) ; 0.647 um
(define W-gam2 (* 1.281 eV-um-scale))
(define W-sig2 (/ (* W-f2 (sqr W-plasma-frq)) (sqr W-frq2)))

(define W-f3 0.706)
(define W-frq3 (* 3.580 eV-um-scale)) ; 0.346 um
(define W-gam3 (* 3.332 eV-um-scale))
(define W-sig3 (/ (* W-f3 (sqr W-plasma-frq)) (sqr W-frq3)))

(define W-f4 2.590)
(define W-frq4 (* 7.498 eV-um-scale)) ; 0.165 um
(define W-gam4 (* 5.836 eV-um-scale))
(define W-sig4 (/ (* W-f4 (sqr W-plasma-frq)) (sqr W-frq4)))

(define W (make medium (epsilon 1.0)
  (E-susceptibilities
     (make drude-susceptibility
       (frequency W-frq0) (gamma W-gam0) (sigma W-sig0))
     (make lorentzian-susceptibility
       (frequency W-frq1) (gamma W-gam1) (sigma W-sig1))
     (make lorentzian-susceptibility
       (frequency W-frq2) (gamma W-gam2) (sigma W-sig2))
     (make lorentzian-susceptibility
       (frequency W-frq3) (gamma W-gam3) (sigma W-sig3))
     (make lorentzian-susceptibility
       (frequency W-frq4) (gamma W-gam4) (sigma W-sig4)))))

;------------------------------------------------------------------
; metals from D. Barchiesi and T. Grosges, J. Nanophotonics, Vol. 8, 08996, 2015
; wavelength range: 0.4 - 0.8 um

; gold (Au)
; fit to P.B. Johnson and R.W. Christy, Physical Review B, Vol. 6, pp. 4370-9, 1972

(define Au-JC-visible-frq0 (/ (* 0.139779231751333 um-scale)))
(define Au-JC-visible-gam0 (/ (* 26.1269913352870 um-scale)))
(define Au-JC-visible-sig0 1)

(define Au-JC-visible-frq1 (/ (* 0.404064525036786 um-scale)))
(define Au-JC-visible-gam1 (/ (* 1.12834046202759 um-scale)))
(define Au-JC-visible-sig1 2.07118534879440)

(define Au-JC-visible (make medium (epsilon 6.1599)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency Au-JC-visible-frq0) (gamma Au-JC-visible-gam0) (sigma Au-JC-visible-sig0))
   (make lorentzian-susceptibility
     (frequency Au-JC-visible-frq1) (gamma Au-JC-visible-gam1) (sigma Au-JC-visible-sig1)))))

;------------------------------------------------------------------
; gold (Au)
; fit to E.D. Palik, Handbook of Optical Constants, Academic Press, 1985

(define Au-visible-frq0 (/ (* 0.0473629248511456 um-scale)))
(define Au-visible-gam0 (/ (* 0.255476199605166 um-scale)))
(define Au-visible-sig0 1)

(define Au-visible-frq1 (/ (* 0.800619321082804 um-scale)))
(define Au-visible-gam1 (/ (* 0.381870287531951 um-scale)))
(define Au-visible-sig1 -169.060953137985)

(define Au-visible (make medium (epsilon 0.6888)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency Au-visible-frq0) (gamma Au-visible-gam0) (sigma Au-visible-sig0))
   (make lorentzian-susceptibility
     (frequency Au-visible-frq1) (gamma Au-visible-gam1) (sigma Au-visible-sig1)))))

;------------------------------------------------------------------                       
;; UNSTABLE: field divergence may occur

; silver (Au)
; fit to E.D. Palik, Handbook of Optical Constants, Academic Press, 1985

(define Ag-visible-frq0 (/ (* 0.142050162130618 um-scale)))
(define Ag-visible-gam0 (/ (* 18.0357292925015 um-scale)))
(define Ag-visible-sig0 1)

(define Ag-visible-frq1 (/ (* 0.115692151792108 um-scale)))
(define Ag-visible-gam1 (/ (* 0.257794324096575 um-scale)))
(define Ag-visible-sig1 3.74465275944019)

(define Ag-visible (make medium (epsilon 0.0067526)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency Ag-visible-frq0) (gamma Ag-visible-gam0) (sigma Ag-visible-sig0))
   (make lorentzian-susceptibility
     (frequency Ag-visible-frq1) (gamma Ag-visible-gam1) (sigma Ag-visible-sig1)))))

;------------------------------------------------------------------
;; UNSTABLE: field divergence may occur

; aluminum (Al)
; fit to E.D. Palik, Handbook of Optical Constants, Academic Press, 1985

(define Al-visible-frq0 (/ (* 0.0625841659042985 um-scale)))
(define Al-visible-gam0 (/ (* 0.606007002962666 um-scale)))
(define Al-visible-sig0 1)

(define Al-visible-frq1 (/ (* 0.528191199577075 um-scale)))
(define Al-visible-gam1 (/ (* 0.291862527666814 um-scale)))
(define Al-visible-sig1 -44.4456675577921)

(define Al-visible (make medium (epsilon 0.13313)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency Al-visible-frq0) (gamma Al-visible-gam0) (sigma Al-visible-sig0))
   (make lorentzian-susceptibility
     (frequency Al-visible-frq1) (gamma Al-visible-gam1) (sigma Al-visible-sig1)))))

;------------------------------------------------------------------
; chromium (Cr)
; fit to E.D. Palik, Handbook of Optical Constants, Academic Press, 1985

(define Cr-visible-frq0 (/ (* 0.118410119507342 um-scale)))
(define Cr-visible-gam0 (/ (* 0.628596264869804 um-scale)))
(define Cr-visible-sig0 1)

(define Cr-visible-frq1 (/ (* 0.565709598452496 um-scale)))
(define Cr-visible-gam1 (/ (* 0.731117670900812 um-scale)))
(define Cr-visible-sig1 13.2912419951294)

(define Cr-visible (make medium (epsilon 2.7767)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency Cr-visible-frq0) (gamma Cr-visible-gam0) (sigma Cr-visible-sig0))
   (make lorentzian-susceptibility
     (frequency Cr-visible-frq1) (gamma Cr-visible-gam1) (sigma Cr-visible-sig1)))))

;------------------------------------------------------------------
;; UNSTABLE: field divergence may occur

; titanium (Ti)
; fit to E.D. Palik, Handbook of Optical Constants, Academic Press, 1985

(define Ti-visible-frq0 (/ (* 0.101331651921602 um-scale)))
(define Ti-visible-gam0 (/ (* 0.365743382258719 um-scale)))
(define Ti-visible-sig0 1)

(define Ti-visible-frq1 (/ (* 4.56839173979216e-09 um-scale)))
(define Ti-visible-gam1 (/ (* 5.86441957443603e-10  um-scale)))
(define Ti-visible-sig1 54742662.1963414)

(define Ti-visible (make medium (epsilon -5.4742e7)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency Ti-visible-frq0) (gamma Ti-visible-gam0) (sigma Ti-visible-sig0))
   (make lorentzian-susceptibility
     (frequency Ti-visible-frq1) (gamma Ti-visible-gam1) (sigma Ti-visible-sig1)))))

;------------------------------------------------------------------
; aluminum (Al) from Horiba Technical Note 09: Drude Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Drude_Dispersion_Model.pdf
; wavelength range: 0.19 - 0.83 um

(define Al-drude-frq (/ (* 0.0789607648707171 um-scale)))
(define Al-drude-gam (/ (* 1.78138208333333 um-scale)))
(define Al-drude-sig 1)

(define Al-drude (make medium (epsilon 1.0)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency Al-drude-frq) (gamma Al-drude-gam) (sigma Al-drude-sig)))))

;------------------------------------------------------------------
; cobalt (Co) from Horiba Technical Note 09: Drude Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Drude_Dispersion_Model.pdf
; wavelength range: 0.26 - 1.65 um

(define Co-frq (/ (* 0.0789607648707171 um-scale)))
(define Co-gam (/ (* 0.213802712536644 um-scale)))
(define Co-sig 1)

(define Co (make medium (epsilon 3.694)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency Co-frq) (gamma Co-gam) (sigma Co-sig)))))

;------------------------------------------------------------------
;; WARNING: unstable; field divergence may occur

; molybdenum (Mo) from Horiba Technical Note 09: Drude Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Drude_Dispersion_Model.pdf
; wavelength range: 0.25 - 0.83 um

(define Mo-frq (/ (* 0.0620790071099539 um-scale)))
(define Mo-gam (/ (* 0.148359690080172 um-scale)))
(define Mo-sig 1)

(define Mo (make medium (epsilon -1.366)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency Mo-frq) (gamma Mo-gam) (sigma Mo-sig)))))

;------------------------------------------------------------------
; nickel chrome (NiCr) from Horiba Technical Note 09: Drude Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Drude_Dispersion_Model.pdf
; wavelength range: 0.25 - 0.83 um

(define NiCr-frq (/ (* 0.0868845080588648 um-scale)))
(define NiCr-gam (/ (* 0.308418390547264 um-scale)))
(define NiCr-sig 1)

(define NiCr (make medium (epsilon 1.0)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency NiCr-frq) (gamma NiCr-gam) (sigma NiCr-sig)))))

;------------------------------------------------------------------
; nickel iron (NiFe) from Horiba Technical Note 09: Drude Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Drude_Dispersion_Model.pdf
; wavelength range: 0.25 - 0.83 um

(define NiFe-frq (/ (* 0.0838297450980392 um-scale)))
(define NiFe-gam (/ (* 0.259381156903766 um-scale)))
(define NiFe-sig 1)

(define NiFe (make medium (epsilon 1.0)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency NiFe-frq) (gamma NiFe-gam) (sigma NiFe-sig)))))

;------------------------------------------------------------------
; titanium (Ti) from Horiba Technical Note 09: Drude Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Drude_Dispersion_Model.pdf
; wavelength range: 0.21 - 1.24 um

(define Ti-drude-frq (/ (* 0.113746966055046 um-scale)))
(define Ti-drude-gam (/ (* 0.490056098814229 um-scale)))
(define Ti-drude-sig 1)

(define Ti-drude (make medium (epsilon 1.0)
  (E-susceptibilities
   (make drude-susceptibility
     (frequency Ti-drude-frq) (gamma Ti-drude-gam) (sigma Ti-drude-sig)))))

;------------------------------------------------------------------
; silicon nitride (SiN), non-stoichiometric, from Horiba Technical Note 08: Lorentz Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Lorentz_Dispersion_Model.pdf
; wavelength range: 0.21 - 2.07 um

(define SiN-frq1 (/ (* 0.190891752117013 um-scale)))
(define SiN-gam1 (/ (* 3.11518072864322 um-scale)))
(define SiN-sig1 1.2650)

(define SiN (make medium (epsilon 2.320)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency SiN-frq1) (gamma SiN-gam1) (sigma SiN-sig1)))))

;------------------------------------------------------------------
; silicon nitride (Si3N4), stoichiometric, from Horiba Technical Note 08: Lorentz Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Lorentz_Dispersion_Model.pdf
; wavelength range: 0.23 - 0.83 um

(define Si3N4-frq1 (/ (* 0.389153148148148 um-scale)))
(define Si3N4-gam1 (/ (* 0.693811936205932 um-scale)))
(define Si3N4-sig1 4.377)

(define Si3N4 (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency Si3N4-frq1) (gamma Si3N4-gam1) (sigma Si3N4-sig1)))))

;------------------------------------------------------------------
; silicon dioxide (SiO2) from Horiba Technical Note 08: Lorentz Dispersion Model
; ref: http://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/Lorentz_Dispersion_Model.pdf
; wavelength range: 0.25 - 1.77 um

(define SiO2-frq1 (/ (* 0.103320160833333 um-scale)))
(define SiO2-gam1 (/ (* 12.3984193000000 um-scale)))
(define SiO2-sig1 1.12)

(define SiO2 (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency SiO2-frq1) (gamma SiO2-gam1) (sigma SiO2-sig1)))))

;------------------------------------------------------------------
; indium phosphide (InP) from Handbook of Optics, 2nd edition, Vol. 2, McGraw-Hill, 1994
; ref: https://refractiveindex.info/?shelf=main&book=InP&page=Pettit
; wavelength range: 0.95 - 10 um

(define InP-frq1 (/ (* 0.6263 um-scale)))
(define InP-gam1 0)
(define InP-sig1 2.316)

(define InP-frq2 (/ (* 32.935 um-scale)))
(define InP-gam2 0)
(define InP-sig2 2.765)

(define InP (make medium (epsilon 7.255)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency InP-frq1) (gamma InP-gam1) (sigma InP-sig1))
  (make lorentzian-susceptibility
    (frequency InP-frq2) (gamma InP-gam2) (sigma InP-sig2)))))

;------------------------------------------------------------------
; germanium (Ge) from N. P. Barnes and M. S. Piltch, J. Optical Society America, Vol. 69, pp. 178-180, 1979
; ref: https://refractiveindex.info/?shelf=main&book=Ge&page=Icenogle
; wavelength range: 2.5 - 12 um

(define Ge-frq1 (/ (* 0.6641159 um-scale)))
(define Ge-gam1 0)
(define Ge-sig1 6.7288)

(define Ge-frq2 (/ (* 62.210127 um-scale)))
(define Ge-gam2 0)
(define Ge-sig2 0.21307)

(define Ge (make medium (epsilon 9.28156)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency Ge-frq1) (gamma Ge-gam1) (sigma Ge-sig1))
  (make lorentzian-susceptibility
    (frequency Ge-frq2) (gamma Ge-gam2) (sigma Ge-sig2)))))

;------------------------------------------------------------------
; silicon (Si) from C. D. Salzberg and J. J. Villa, , J. Optical Society America, Vol. 47, pp. 244-246, 1957
; ref: https://refractiveindex.info/?shelf=main&book=Si&page=Salzberg
; wavelength range: 1.36 - 11 um

(define Si-frq1 (/ (* 0.301516485 um-scale)))
(define Si-gam1 0)
(define Si-sig1 10.6684293)

(define Si-frq2 (/ (* 1.13475115 um-scale)))
(define Si-gam2 0)
(define Si-sig2 0.0030434748)

(define Si-frq3 (/ (* 1104 um-scale)))
(define Si-gam3 0)
(define Si-sig3 1.54133408)

(define Si (make medium (epsilon 9.28156)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency Si-frq1) (gamma Si-gam1) (sigma Si-sig1))
  (make lorentzian-susceptibility
    (frequency Si-frq2) (gamma Si-gam2) (sigma Si-sig2))
  (make lorentzian-susceptibility
    (frequency Si-frq3) (gamma Si-gam3) (sigma Si-sig3)))))

;------------------------------------------------------------------
; poly(methyl methacrylate) (PMMA) from N. Sultanova et al., Acta Physica Polonica A, Vol. 116, pp. 585-7, 2009
; ref: https://refractiveindex.info/?shelf=organic&book=poly%28methyl_methacrylate%29&page=Sultanova
; wavelength range: 0.437 - 1.052 um

(define PMMA-frq1 (/ (* 0.106362587407415 um-scale)))
(define PMMA-gam1 0)
(define PMMA-sig1 1.1819)

(define PMMA (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency PMMA-frq1) (gamma PMMA-gam1) (sigma PMMA-sig1)))))

;------------------------------------------------------------------
; polycarbonate (PC) from N. Sultanova et al., Acta Physica Polonica A, Vol. 116, pp. 585-7, 2009
; ref: https://refractiveindex.info/?shelf=organic&book=polycarbonate&page=Sultanova
; wavelength range: 0.437 - 1.052 um

(define PC-frq1 (/ (* 0.145958898324152 um-scale)))
(define PC-gam1 0)
(define PC-sig1 1.4182)

(define PC (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency PC-frq1) (gamma PC-gam1) (sigma PC-sig1)))))

;------------------------------------------------------------------
; polystyrene (PS) from N. Sultanova et al., Acta Physica Polonica A, Vol. 116, pp. 585-7, 2009
; ref: https://refractiveindex.info/?shelf=organic&book=polystyren&page=Sultanova
; wavelength range: 0.437 - 1.052 um

(define PS-frq1 (/ (* 0.142182980697410 um-scale)))
(define PS-gam1 0)
(define PS-sig1 1.4435)

(define PS (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency PS-frq1) (gamma PS-gam1) (sigma PS-sig1)))))

;------------------------------------------------------------------
; cellulose (CLS) from N. Sultanova et al., Acta Physica Polonica A, Vol. 116, pp. 585-7, 2009
; ref: https://refractiveindex.info/?shelf=organic&book=cellulose&page=Sultanova
; wavelength range: 0.437 - 1.052 um

(define CLS-frq1 (/ (* 0.105294824184287 um-scale)))
(define CLS-gam1 0)
(define CLS-sig1 1.124)

(define CLS (make medium (epsilon 1.0)
 (E-susceptibilities 
  (make lorentzian-susceptibility
    (frequency CLS-frq1) (gamma CLS-gam1) (sigma CLS-sig1)))))

;------------------------------------------------------------------
; barium borate (BaB2O4), beta phase, from G. Tamosauskas et al., Optical Materials Express, Vol. 8, pp. 1410-18, 2018
; ref: https://refractiveindex.info/?shelf=main&book=BaB2O4&page=Tamosauskas-o
; ref: https://refractiveindex.info/?shelf=main&book=BaB2O4&page=Tamosauskas-e
; wavelength range: 0.188 - 5.2 um

;; NOTE: ordinary (o) axes in X and Y, extraordinary (e) axis in Z

(define BaB2O4-frq1 (/ (* 0.06265780079128216 um-scale)))
(define BaB2O4-gam1 0)
(define BaB2O4-sig1 0.90291)

(define BaB2O4-frq2 (/ (* 0.13706202975295528 um-scale)))
(define BaB2O4-gam2 0)
(define BaB2O4-sig2 0.83155)

(define BaB2O4-frq3 (/ (* 7.746612162745725 um-scale)))
(define BaB2O4-gam3 0)
(define BaB2O4-sig3 0.76536)

(define BaB2O4-susc-o
  (list
   (make lorentzian-susceptibility
     (frequency BaB2O4-frq1) (gamma BaB2O4-gam1) (sigma-diag (vector3* BaB2O4-sig1 (vector3 1 1 0))))
   (make lorentzian-susceptibility
     (frequency BaB2O4-frq2) (gamma BaB2O4-gam2) (sigma-diag (vector3* BaB2O4-sig2 (vector3 1 1 0))))
   (make lorentzian-susceptibility
     (frequency BaB2O4-frq3) (gamma BaB2O4-gam3) (sigma-diag (vector3* BaB2O4-sig3 (vector3 1 1 0))))))

(define BaB2O4-frq1 (/ (* 0.0845103543951864 um-scale)))
(define BaB2O4-gam1 0)
(define BaB2O4-sig1 1.151075)

(define BaB2O4-frq2 (/ (* 0.15029970059850417 um-scale)))
(define BaB2O4-gam2 0)
(define BaB2O4-sig2 0.21803)

(define BaB2O4-frq3 (/ (* 16.217274740226856 um-scale)))
(define BaB2O4-gam3 0)
(define BaB2O4-sig3 0.656)

(define BaB2O4-susc-e
  (list
   (make lorentzian-susceptibility
     (frequency BaB2O4-frq1) (gamma BaB2O4-gam1) (sigma-diag (vector3* BaB2O4-sig1 (vector3 0 0 1))))
   (make lorentzian-susceptibility
     (frequency BaB2O4-frq2) (gamma BaB2O4-gam2) (sigma-diag (vector3* BaB2O4-sig2 (vector3 0 0 1))))
   (make lorentzian-susceptibility
     (frequency BaB2O4-frq3) (gamma BaB2O4-gam3) (sigma-diag (vector3* BaB2O4-sig3 (vector3 0 0 1))))))

(define BaB2O4 (make medium
		 (epsilon 1.0)
		 (E-susceptibilities (append BaB2O4-susc-o BaB2O4-susc-e))))

;------------------------------------------------------------------
; lithium niobate (LiNbO3) from D.E. Zelmon et al., J. Optical Society of America B, Vol. 14, pp. 3319-22, 1997
; ref: https://refractiveindex.info/?shelf=main&book=LiNbO3&page=Zelmon-o
; ref: https://refractiveindex.info/?shelf=main&book=LiNbO3&page=Zelmon-e
; wavelength range: 0.4 - 5.0 um

;; NOTE: ordinary (o) axes in X and Y, extraordinary (e) axis in Z

(define LiNbO3-frq1 (/ (* 0.13281566172707193 um-scale)))
(define LiNbO3-gam1 0)
(define LiNbO3-sig1 2.6734)

(define LiNbO3-frq2 (/ (* 0.24318717071424636 um-scale)))
(define LiNbO3-gam2 0)
(define LiNbO3-sig2 1.2290)

(define LiNbO3-frq3 (/ (* 21.78531615561271 um-scale)))
(define LiNbO3-gam3 0)
(define LiNbO3-sig3 12.614)

(define LiNbO3-susc-o
  (list
   (make lorentzian-susceptibility
     (frequency LiNbO3-frq1) (gamma LiNbO3-gam1) (sigma-diag (vector3* LiNbO3-sig1 (vector3 1 1 0))))
   (make lorentzian-susceptibility
     (frequency LiNbO3-frq2) (gamma LiNbO3-gam2) (sigma-diag (vector3* LiNbO3-sig2 (vector3 1 1 0))))
   (make lorentzian-susceptibility
     (frequency LiNbO3-frq3) (gamma LiNbO3-gam3) (sigma-diag (vector3* LiNbO3-sig3 (vector3 1 1 0))))))

(define LiNbO3-frq1 (/ (* 0.14307340773183533 um-scale)))
(define LiNbO3-gam1 0)
(define LiNbO3-sig1 2.9804)

(define LiNbO3-frq2 (/ (* 0.2580697580112788 um-scale)))
(define LiNbO3-gam2 0)
(define LiNbO3-sig2 0.5981)

(define LiNbO3-frq3 (/ (* 20.39803912144498 um-scale)))
(define LiNbO3-gam3 0)
(define LiNbO3-sig3 8.9543)

(define LiNbO3-susc-e
  (list
   (make lorentzian-susceptibility
     (frequency LiNbO3-frq1) (gamma LiNbO3-gam1) (sigma-diag (vector3* LiNbO3-sig1 (vector3 0 0 1))))
   (make lorentzian-susceptibility
     (frequency LiNbO3-frq2) (gamma LiNbO3-gam2) (sigma-diag (vector3* LiNbO3-sig2 (vector3 0 0 1))))
   (make lorentzian-susceptibility
     (frequency LiNbO3-frq3) (gamma LiNbO3-gam3) (sigma-diag (vector3* LiNbO3-sig3 (vector3 0 0 1))))))

(define LiNbO3 (make medium
		 (epsilon 1.0)
		 (E-susceptibilities (append LiNbO3-susc-o LiNbO3-susc-e))))

;------------------------------------------------------------------
; calcium tungstate (CaWO4) from W.L. Bond, J. Applied Physics, Vol. 36, pp. 1674-77, 1965
; ref: https://refractiveindex.info/?shelf=main&book=CaWO4&page=Bond-o
; ref: https://refractiveindex.info/?shelf=main&book=CaWO4&page=Bond-e
; wavelength range: 0.45 - 4.0 um

;; NOTE: ordinary (o) axes in X and Y, extraordinary (e) axis in Z

(define CaWO4-frq1 (/ (* 0.1347 um-scale)))
(define CaWO4-gam1 0)
(define CaWO4-sig1 2.5493)

(define CaWO4-frq2 (/ (* 10.815 um-scale)))
(define CaWO4-gam2 0)
(define CaWO4-sig2 0.9200)

(define CaWO4-susc-o
  (list
   (make lorentzian-susceptibility
     (frequency CaWO4-frq1) (gamma CaWO4-gam1) (sigma-diag (vector3* CaWO4-sig1 (vector3 1 1 0))))
   (make lorentzian-susceptibility
     (frequency CaWO4-frq2) (gamma CaWO4-gam2) (sigma-diag (vector3* CaWO4-sig2 (vector3 1 1 0))))))

(define CaWO4-frq1 (/ (* 0.1379 um-scale)))
(define CaWO4-gam1 0)
(define CaWO4-sig1 2.6041)

(define CaWO4-frq2 (/ (* 21.371 um-scale)))
(define CaWO4-gam2 0)
(define CaWO4-sig2 4.1237)

(define CaWO4-susc-e
  (list
   (make lorentzian-susceptibility
     (frequency CaWO4-frq1) (gamma CaWO4-gam1) (sigma-diag (vector3* CaWO4-sig1 (vector3 0 0 1))))
   (make lorentzian-susceptibility
     (frequency CaWO4-frq2) (gamma CaWO4-gam2) (sigma-diag (vector3* CaWO4-sig2 (vector3 0 0 1))))))

(define CaWO4 (make medium
		(epsilon 1.0)
		(E-susceptibilities (append CaWO4-susc-o CaWO4-susc-e))))

;------------------------------------------------------------------
; calcium carbonate (CaCO3) from G. Ghosh, Optics Communication, Vol. 163, pp. 95-102, 1999
; ref: https://refractiveindex.info/?shelf=main&book=CaCO3&page=Ghosh-o
; ref: https://refractiveindex.info/?shelf=main&book=CaCO3&page=Ghosh-e
; wavelength range: 0.204 - 2.172 um

;; NOTE: ordinary (o) axes in X and Y, extraordinary (e) axis in Z

(define CaCO3-frq1 (/ (* 0.13940057496294625 um-scale)))
(define CaCO3-gam1 0)
(define CaCO3-sig1 0.96464345)

(define CaCO3-frq2 (/ (* 10.954451150103322 um-scale)))
(define CaCO3-gam2 0)
(define CaCO3-sig2 1.82831454)

(define CaCO3-susc-o
  (list
   (make lorentzian-susceptibility
     (frequency CaCO3-frq1) (gamma CaCO3-gam1) (sigma-diag (vector3* CaCO3-sig1 (vector3 1 1 0))))
   (make lorentzian-susceptibility
     (frequency CaCO3-frq2) (gamma CaCO3-gam2) (sigma-diag (vector3* CaCO3-sig2 (vector3 1 1 0))))))

(define CaCO3-frq1 (/ (* 0.1032906302623815 um-scale)))
(define CaCO3-gam1 0)
(define CaCO3-sig1 0.82427830)

(define CaCO3-frq2 (/ (* 10.954451150103322 um-scale)))
(define CaCO3-gam2 0)
(define CaCO3-sig2 0.14429128)

(define CaCO3-susc-e
  (list
   (make lorentzian-susceptibility
     (frequency CaCO3-frq1) (gamma CaCO3-gam1) (sigma-diag (vector3* CaCO3-sig1 (vector3 0 0 1))))
   (make lorentzian-susceptibility
     (frequency CaCO3-frq2) (gamma CaCO3-gam2) (sigma-diag (vector3* CaCO3-sig2 (vector3 0 0 1))))))

(define CaCO3 (make medium
		(epsilon-diag 1.73358749 1.73358749 1.35859695)
		(E-susceptibilities (append CaCO3-susc-o CaCO3-susc-e))))

;------------------------------------------------------------------
; silicon dioxide (SiO2) from G. Ghosh, Optics Communication, Vol. 163, pp. 95-102, 1999
; ref: https://refractiveindex.info/?shelf=main&book=SiO2&page=Ghosh-o
; ref: https://refractiveindex.info/?shelf=main&book=SiO2&page=Ghosh-e
; wavelength range: 0.198 - 2.0531 um

;; NOTE: ordinary (o) axes in X and Y, extraordinary (e) axis in Z

(define SiO2-frq1 (/ (* 0.10029257051247614 um-scale)))
(define SiO2-gam1 0)
(define SiO2-sig1 1.07044083)

(define SiO2-frq2 (/ (* 10 um-scale)))
(define SiO2-gam2 0)
(define SiO2-sig2 1.10202242)

(define SiO2-susc-o
  (list
   (make lorentzian-susceptibility
     (frequency SiO2-frq1) (gamma SiO2-gam1) (sigma-diag (vector3* SiO2-sig1 (vector3 1 1 0))))
   (make lorentzian-susceptibility
     (frequency SiO2-frq2) (gamma SiO2-gam2) (sigma-diag (vector3* SiO2-sig2 (vector3 1 1 0))))))

(define SiO2-frq1 (/ (* 0.10104546699382412 um-scale)))
(define SiO2-gam1 0)
(define SiO2-sig1 1.09509924)

(define SiO2-frq2 (/ (* 10 um-scale)))
(define SiO2-gam2 0)
(define SiO2-sig2 1.15662475)

(define SiO2-susc-e
  (list
   (make lorentzian-susceptibility
     (frequency SiO2-frq1) (gamma SiO2-gam1) (sigma-diag (vector3* SiO2-sig1 (vector3 0 0 1))))
   (make lorentzian-susceptibility
     (frequency SiO2-frq2) (gamma SiO2-gam2) (sigma-diag (vector3* SiO2-sig2 (vector3 0 0 1))))))

(define SiO2-aniso (make medium
		     (epsilon-diag 1.28604141 1.28604141 1.28851804)
		     (E-susceptibilities (append SiO2-susc-o SiO2-susc-e))))

;------------------------------------------------------------------
; gallium nitride (GaN), alpha phase (wurtzite), from A.S. Barker Jr. and M. Ilegems, Physical Review B, Vol. 7, pp. 743-50, 1973
; ref: https://refractiveindex.info/?shelf=main&book=GaN&page=Barker-o
; ref: https://refractiveindex.info/?shelf=main&book=GaN&page=Barker-e
; wavelength range: 0.35 - 10 um

;; NOTE: ordinary (o) axes in X and Y, extraordinary (e) axis in Z

(define GaN-frq1 (/ (* 0.256 um-scale)))
(define GaN-gam1 0)
(define GaN-sig1 1.75)

(define GaN-frq2 (/ (* 17.86 um-scale)))
(define GaN-gam2 0)
(define GaN-sig2 4.1)

(define GaN-susc-o
  (list
   (make lorentzian-susceptibility
     (frequency GaN-frq1) (gamma GaN-gam1) (sigma-diag (vector3* GaN-sig1 (vector3 1 1 0))))
   (make lorentzian-susceptibility
     (frequency GaN-frq2) (gamma GaN-gam2) (sigma-diag (vector3* GaN-sig2 (vector3 1 1 0))))))

(define GaN-frq1 (/ (* 18.76 um-scale)))
(define GaN-gam1 0)
(define GaN-sig1 5.08)

(define GaN-susc-e
  (list
   (make lorentzian-susceptibility
     (frequency GaN-frq1) (gamma GaN-gam1) (sigma-diag (vector3* GaN-sig1 (vector3 0 0 1))))))

(define GaN (make medium
	      (epsilon-diag 3.6 3.6 5.35)
	      (E-susceptibilities (append GaN-susc-o GaN-susc-e))))

;------------------------------------------------------------------
; aluminum nitride (AlN) from J. Pastrnak and L. Roskovcova, Physica Status Solidi, Vol. 14, K5-8, 1966
; ref: https://refractiveindex.info/?shelf=main&book=AlN&page=Pastrnak-o
; ref: https://refractiveindex.info/?shelf=main&book=AlN&page=Pastrnak-e
; wavelength range: 0.22 - 5 um

;; NOTE: ordinary (o) axes in X and Y, extraordinary (e) axis in Z

(define AlN-frq1 (/ (* 0.1715 um-scale)))
(define AlN-gam1 0)
(define AlN-sig1 1.3786)

(define AlN-frq2 (/ (* 15.03 um-scale)))
(define AlN-gam2 0)
(define AlN-sig2 3.861)

(define AlN-susc-o
  (list
   (make lorentzian-susceptibility
     (frequency AlN-frq1) (gamma AlN-gam1) (sigma-diag (vector3* AlN-sig1 (vector3 1 1 0))))
   (make lorentzian-susceptibility
     (frequency AlN-frq2) (gamma AlN-gam2) (sigma-diag (vector3* AlN-sig2 (vector3 1 1 0))))))

(define AlN-frq1 (/ (* 0.1746 um-scale)))
(define AlN-gam1 0)
(define AlN-sig1 1.6173)

(define AlN-frq2 (/ (* 15.03 um-scale)))
(define AlN-gam2 0)
(define AlN-sig2 4.139)

(define AlN-susc-e
  (list
   (make lorentzian-susceptibility
     (frequency AlN-frq1) (gamma AlN-gam1) (sigma-diag (vector3* AlN-sig1 (vector3 0 0 1))))
   (make lorentzian-susceptibility
     (frequency AlN-frq2) (gamma AlN-gam2) (sigma-diag (vector3* AlN-sig2 (vector3 0 0 1))))))

(define AlN-aniso (make medium
		     (epsilon-diag 3.1399 3.1399 3.0729)
		     (E-susceptibilities (append AlN-susc-o AlN-susc-e))))

;------------------------------------------------------------------
; alumina/sapphire (Al2O3) from I.H. Malitson and M.J. Dodge, J. Optical Society of America, Vol. 62, pp. 1405, 1972
; ref: https://refractiveindex.info/?shelf=main&book=Al2O3&page=Malitson-o
; ref: https://refractiveindex.info/?shelf=main&book=Al2O3&page=Malitson-e
; wavelength range: 0.2 - 5 um

;; NOTE: ordinary (o) axes in X and Y, extraordinary (e) axis in Z

(define Al2O3-frq1 (/ (* 0.0726631 um-scale)))
(define Al2O3-gam1 0)
(define Al2O3-sig1 1.4313493)

(define Al2O3-frq2 (/ (* 0.1193242 um-scale)))
(define Al2O3-gam2 0)
(define Al2O3-sig2 0.65054713)

(define Al2O3-frq3 (/ (* 18.02825 um-scale)))
(define Al2O3-gam3 0)
(define Al2O3-sig3 5.3414021)

(define Al2O3-susc-o
  (list
   (make lorentzian-susceptibility
     (frequency Al2O3-frq1) (gamma Al2O3-gam1) (sigma-diag (vector3* Al2O3-sig1 (vector3 1 1 0))))
   (make lorentzian-susceptibility
     (frequency Al2O3-frq2) (gamma Al2O3-gam2) (sigma-diag (vector3* Al2O3-sig2 (vector3 1 1 0))))
   (make lorentzian-susceptibility
     (frequency Al2O3-frq3) (gamma Al2O3-gam3) (sigma-diag (vector3* Al2O3-sig3 (vector3 1 1 0))))))

(define Al2O3-frq1 (/ (* 0.0740288 um-scale)))
(define Al2O3-gam1 0)
(define Al2O3-sig1 1.5039759)

(define Al2O3-frq2 (/ (* 0.1216529 um-scale)))
(define Al2O3-gam2 0)
(define Al2O3-sig2 0.55069141)

(define Al2O3-frq3 (/ (* 20.072248 um-scale)))
(define Al2O3-gam3 0)
(define Al2O3-sig3 6.5927379)

(define Al2O3-susc-e
  (list
   (make lorentzian-susceptibility
     (frequency Al2O3-frq1) (gamma Al2O3-gam1) (sigma-diag (vector3* Al2O3-sig1 (vector3 0 0 1))))
   (make lorentzian-susceptibility
     (frequency Al2O3-frq2) (gamma Al2O3-gam2) (sigma-diag (vector3* Al2O3-sig2 (vector3 0 0 1))))
   (make lorentzian-susceptibility
     (frequency Al2O3-frq3) (gamma Al2O3-gam3) (sigma-diag (vector3* Al2O3-sig3 (vector3 0 0 1))))))

(define Al2O3-aniso (make medium
		      (epsilon 1.0)
		      (E-susceptibilities (append Al2O3-susc-o Al2O3-susc-e))))

;------------------------------------------------------------------
; yttrium oxide (Y2O3) from Y. Nigara, Japanese J. of Applied Physics, Vol. 7, pp. 404-8, 1968
; ref: https://refractiveindex.info/?shelf=main&book=Y2O3&page=Nigara
; wavelength range: 0.25 - 9.6 um

(define Y2O3-frq1 (/ (* 0.1387 um-scale)))
(define Y2O3-gam1 0)
(define Y2O3-sig1 2.578)

(define Y2O3-frq2 (/ (* 22.936 um-scale)))
(define Y2O3-gam2 0)
(define Y2O3-sig2 3.935)

(define Y2O3 (make medium (epsilon 1.0)
 (E-susceptibilities
  (make lorentzian-susceptibility
    (frequency Y2O3-frq1) (gamma Y2O3-gam1) (sigma Y2O3-sig1))
  (make lorentzian-susceptibility
    (frequency Y2O3-frq2) (gamma Y2O3-gam2) (sigma Y2O3-sig2)))))

;------------------------------------------------------------------
; undoped yttrium aluminum garnet (YAG) from D.E. Zelmon et al., Applied Optics, Vol. 37, 4933-5, 1998
; ref: https://refractiveindex.info/?shelf=main&book=Y3Al5O12&page=Zelmon
; wavelength range: 0.4 - 5.0 um

(define YAG-frq1 (/ (* 0.1088577052853862 um-scale)))
(define YAG-gam1 0)
(define YAG-sig1 2.28200)

(define YAG-frq2 (/ (* 16.814695953242804 um-scale)))
(define YAG-gam2 0)
(define YAG-sig2 3.27644)

(define YAG (make medium (epsilon 1.0)
 (E-susceptibilities
  (make lorentzian-susceptibility
    (frequency YAG-frq1) (gamma YAG-gam1) (sigma YAG-sig1))
  (make lorentzian-susceptibility
    (frequency YAG-frq2) (gamma YAG-gam2) (sigma YAG-sig2)))))

;------------------------------------------------------------------
; cadmium telluride (CdTe) from D.T.F. Marple, J. Applied Physics, Vol. 35, pp. 539-42, 1964
; ref: https://refractiveindex.info/?shelf=main&book=CdTe&page=Marple
; wavelength range: 0.86 - 2.5 um

(define CdTe-frq1 (/ (* 0.6049793384901669 um-scale)))
(define CdTe-gam1 0)
(define CdTe-sig1 1.53)

(define CdTe (make medium (epsilon 5.68)
 (E-susceptibilities
  (make lorentzian-susceptibility
    (frequency CdTe-frq1) (gamma CdTe-gam1) (sigma CdTe-sig1)))))