SUPPLEMENTARY INFORMATION

Size: px
Start display at page:

Download "SUPPLEMENTARY INFORMATION"

Transcription

1 Supplementary Information Compact spectrometer based on a disordered photonic chip Brandon Redding, Seng Fatt Liew, Raktim Sarma, Hui Cao* Department of Applied Physics, Yale University, New Haven, CT 06520, USA * hui.cao@yale.edu 1. Spectrum reconstruction algorithm The intensity distribution on the detector array is I = T S, where T is the transmission matrix and S is the input spectrum. After calibrating T, the input spectrum can be reconstructed simply by multiplying the measured I by the inverse of the transmission matrix: S = T -1 I. However, the spectra reconstructed in this manner are susceptible to experimental noise, as shown in Fig. S1(a). To improve the accuracy of the reconstructed spectrum, we used a nonlinear optimization algorithm based on simulated annealing to find the spectrum S which minimized I TS 2. We used the spectrum obtained from the inversion procedure as the initial guess for S and used the simulated annealing algorithm to minimize I TS 2. The final reconstructed spectrum is shown in Figure S1(b). Figure S1 Spectrum reconstruction. The spectrum for a probe signal at λ=1512 nm reconstructed using the matrix inversion procedure (a) and the simulated annealing algorithm (b). In the presence of experimental noise, the simulated annealing algorithm improves the quality of the reconstructed spectrum. 2. Measuring arbitrary spectra with the random spectrometer In addition to accurately identifying the wavelength of narrowband input, the random spectrometer can also measure spectra consisting of multiple wavelengths or even a continuous broad band. Using the random spectrometer shown in Fig. 1, we synthesized various probe spectra as the weighted sum of speckle patterns recorded at individual wavelengths and applied the same reconstruction algorithm as described in Section 1. In Fig. 2(e), we present a reconstructed spectrum of 3 narrow lines with different intensity. The random spectrometer was able to accurately identify the wavelength and relative amplitude of the three lines. Figure 2(f) showed the reconstructed spectrum of a broad continuous spectrum. In this case, since the NATURE PHOTONICS 1

2 spectral feature was broad, we sampled the transmission matrix at the wavelength interval of 1 nm, yielding a transmission matrix with 25 spectral channels and 25 spatial channels (detectors). This result demonstrated that the random spectrometer was able to accurately reconstruct the broadband spectra. 3. Spectral range of the random spectrometer Multiple scattering in a random structure occurs over an extremely broad range of frequencies. This enables the same random spectrometer to operate at different spectral regions simply by using the transmission matrix calibrated for the desired spectral region. The spectral resolution will change for different spectral regions since the transport mean free path of the scattering structure is wavelength dependent. However, the first Mie resonance of the scatterer used in the random spectrometer (75 nm radius air hole in a silicon membrane) is at the wavelength of 440 nm, much shorter than the probe wavelength (~1500 nm). Figure S2 plots the calculated scattering cross section of the air cylinder in the wavelength range of 1250 nm to 1750 nm. The scattering cross section changes by ~50% over the wavelength range of 500 nm. Since the mean free path scales with the scattering cross section, it varies gradually in the wavelength range of 1250 nm to 1750 nm. Therefore, the random spectrometer can operate over a very large spectral range with only a modest change in resolution. Figure S2 Broadband scattering. Calculated scattering cross section of a 75 nm radius air cylinder embedded in a silicon membrane. The scattering cross section decreases gradually from λ=1250 nm to 1750 nm, allowing the random spectrometer to provide similar performance and resolution over a large spectral range. 4. Thermal stability of the random spectrometer One advantage of fabricating the random structure on-chip is that the speckle pattern is very stable, since the scatterers (air holes etched into the silicon membrane) do not move. However, temperature variation will introduce a change in the refractive index of the Si layer and alter the speckle pattern. To estimate the change in temperature that can be tolerated for spectrum reconstruction, we performed numerical simulation of the 20 μm radius random spectrometer using the FDFD method. Specifically, we simulated the speckle pattern generated for a fixed input wavelength (λ = nm) at different temperatures by changing the refractive index of the Si according to the thermo-optic coefficient, dn/dt K -1 [S1]. Figure S3(a) plots the intensity distribution over the detectors as a function of temperature. The speckle pattern gradually decorrelates as the temperature changes. We then reconstructed the input spectra from the speckle patterns using the transmission matrix obtained at one specific temperature. As shown in Fig. S3(b), up to ±4 o K change of temperature, the input wavelength can be recovered 2 NATURE PHOTONICS

3 SUPPLEMENTARY INFORMATION to the accuracy of the spectral resolution (~ 0.5 nm). To achieve better thermal stability, on-chip temperature stabilization techniques, which have been developed for microcavity-based devices, may be used. Alternatively the spectrometer can be calibrated at multiple temperatures, allowing the appropriate transmission matrix to be selected in software as the temperature changes. Figure S3 Modeling thermal stability of the 20 m radius random spectrometer. (a) Calculated intensity of light reaching 20 output waveguides as a function of temperature change for an input signal at λ= nm. As the temperature varies, the speckle pattern gradually decorrelates. (b) Reconstructed spectra for an input light of wavelength λ = nm at different temperatures. When the change of temperature exceeds 4 o K, the recovered wavelength deviates from the input wavelength by 0.5 nm, which is approximately the spectral resolution. 5. Numerical simulation of optical transmission In the random spectrometer, we surrounded the disordered photonic structure with the photonic crystal boundary to reduce the insertion loss. The triangular lattice of air holes that form the boundary were designed to have a 2D photonic bandgap (PBG) for TE polarized light in the wavelength range of 1478 nm 1560 nm. The photonic band structure was calculated with the plane wave expansion method [S2]. To compare the optical transmission with and without the photonic crystal boundary, we performed 2D numerical simulations using the finite-difference frequency-domain method. As shown in Fig. S4, the random spectrometer has the semicircular shape. The input waveguide injects light from the center of the semicircle, and the output waveguides are distributed along the circumference. The density and size of air holes are identical to those of the fabricated spectrometer. The effective index of refraction of the silicon layer is 2.86, it is obtained by matching the PBG of a triangle lattice of air holes (with the same density and size of air holes as in the random array) in the 220 nm silicon layer on top of silica to that of the approximate 2D structure. In Fig. S4(a), the boundary in between the waveguides are open and light can leak out. The boundary at the flat base of the semicircular area containing the scattering medium is also open. In Fig. S4(b), the boundary reflects light back into the random medium. For simplicity, the photonic crystal lattice is replaced by a perfect magnetic conductor at the boundary, which also reflects the TE polarized light from all directions. We summed the flux of light reaching the detectors indicated in green in Fig. S4(a,b) and divided by the input flux to obtain the transmission. For the open random structure of radius 25 μm, the transmission is merely 21% in the wavelength range of 1500 nm to 1525 nm. With the reflecting boundary, the transmission NATURE PHOTONICS 3

4 increases to 60%. In this 2D simulation, the out-of-plane scattering was not taken into account. Nonetheless, the result illustrates that the introduction of reflecting borders can significantly improve the transmission through a scattering medium. Fig. S4 Numerical simulation of optical transmission through a random spectrometer. (a) Schematic of a semicircular random medium with open boundary. The semicircle has a radius of 25 μm. The random medium consists of a random array of air holes with a radius of 75 nm. PML: Perfectly matched layer, PMC: Perfect magnetic conductor. (b) Schematic of the same random medium with reflecting boundary. The photonic crystal structure that formed the boundary in the fabricated spectrometer was replaced with PMC to reduce the computational load. In both (a) and (b), the optical signal enters the random medium through an input waveguide at the center of the semicircle marked by the black arrow. The output waveguides, distributed along the circumference, delivers the transmitted light to the detectors marked by the green line. (c) Calculated amplitude of the H z field for the TE polarized light at λ = 1500 nm diffusing through the scattering medium shown in (a). Without the reflecting boundary, most of the input light leaves the scattering medium through the base of the semicircle before reaching the output waveguides at the circumference. (d) Calculated amplitude of the H z field for the TE polarized light at λ = 1500 nm diffusing through the random medium shown in (b). The reflecting boundary efficiently channels the input light to the output waveguides. The scale bars indicate 10 μm. 6. Estimation of out-of-plane scattering loss Since 3D FDFD simulations of the on-chip spectrometers are computationally heavy, we estimated the relative out-of-plane leakage of the completely random structure, the photonic amorphous structure, and the golden-angle spiral lattice from 2D simulations using the method in Ref. [S3]. For the TE polarized light, we first calculated the in-plane electric field components (E x and E y ) inside the three scattering structures. Then we conducted a 2D spatial Fourier transform of E x and E y and compared the fraction of wavevectors present within the light cone (these wavevectors can leak out of the plane). Since the field distribution obtained with the 2D simulation which neglects out-of-plane scattering differs from the actual field distribution in the presence of out-of-plane scattering, we do not expect to obtain quantitatively accurate 4 NATURE PHOTONICS

5 SUPPLEMENTARY INFORMATION estimations of the out-of-plane scattering loss, but we can make qualitative comparisons of the three scattering media. The spectrometers in the simulation had the geometry shown in Fig. S4(b) with a radius of 20 μm. Although the simulation included the output waveguides, we extracted the fields (E x and E y ) across the scattering structures within the semicircle to estimate the out-of-plane loss. We then computed ( F(E x ) 2 + F(E y ) 2 ), where F represents the 2D Fourier transform. There are two light cones, one for air above the silicon layer ( k > ω/c, ω is the angular frequency and c is the speed of light in air), the other for silica underneath ( k > n s ω/c, n s = 1.5 is the refractive index of silica). Since the light cone for the silica is larger than that of the air, we computed the fraction of the Fourier components within the silica light cone. We found that for the random, amorphous, and spiral structures, 0.82%, 0.59%, and 0.49% of the electric field intensities were within the light cone, respectively. Hence, the random structure is more likely to scatter light into the light cone, leading to larger out-of-plane leakage than the amorphous and spiral structures. 7. Characterization of amorphous and spiral spectrometers In addition to reducing the out-of-plane scattering loss, the spectrometers based on the photonic amorphous structure and the golden-angle spiral lattice had similar performance to the random spectrometer. We performed the same calibration procedure for transmission matrices of the amorphous and spiral structures, and characterized their operation as a spectrometer. Specifically, the transmission matrices for the 20 m radius spectrometers were recorded for TE polarized light in the wavelength range of 1500 nm to 1525 nm in steps of 0.25 nm. We tested their ability to reconstruct a series of narrow lines. As shown in Fig. S5(a,c), both amorphous and spiral spectrometers accurately recovered the spectral lines over 25 nm bandwidth. We also tested their spectral resolution by measuring two closely spaced lines. Figure S5(b,d) exhibited a wavelength resolution of 0.75 nm around = 1512 nm for both spectrometers, which is comparable to the random spectrometer. Figure S5 Characterization of amorphous and spiral based spectrometers. (a,c) Reconstructed spectra for a series of narrow spectral probes across the 25 nm bandwidth using the spectrometers based on the photonic amorphous structure (a) and the golden-angle spiral lattice (c). The black dotted lines mark the center wavelength of each probe line. (b,c) Reconstructed spectrum (blue line) of two narrow spectral lines separated by 0.75 nm using the amorphous (b) and spiral (d) spectrometers. The red dotted lines mark the center wavelengths of the probe lines. NATURE PHOTONICS 5

6 References [S1] Komma, J., Schwarz, C., Hofmann, G., Heinert, D., Nawrodt, R., Thermo-optic coefficient of silicon at 1550 nm and cryogenic temepratures, Applied Physics Letters, 101, (2012). [S2] Johnson, S. G. & Joannopoulos, J. D. Block-iterative frequency-domain methods for Maxwell s equations in a planewave basis. Optics Express 8, (2001). [S3] Boriskina, S., Gopinath, A., Dal Negro, L. Optical gap formation and localization properties of optical modes in deterministic aperiodic photonic structures. Optics Express 16, (2008). 6 NATURE PHOTONICS

Optimization of metallic biperiodic photonic crystals. Application to compact directive antennas

Optimization of metallic biperiodic photonic crystals. Application to compact directive antennas Optimization of metallic biperiodic photonic crystals Application to compact directive antennas Nicolas Guérin Computational Optic Groups (COG) IFH, ETHZ, http://alphard.ethz.ch Keywords: 3D modeling,

More information

NEAR-IR BROADBAND POLARIZER DESIGN BASED ON PHOTONIC CRYSTALS

NEAR-IR BROADBAND POLARIZER DESIGN BASED ON PHOTONIC CRYSTALS U.P.B. Sci. Bull., Series A, Vol. 77, Iss. 3, 2015 ISSN 1223-7027 NEAR-IR BROADBAND POLARIZER DESIGN BASED ON PHOTONIC CRYSTALS Bogdan Stefaniţă CALIN 1, Liliana PREDA 2 We have successfully designed a

More information

Reproducing the hierarchy of disorder for Morpho-inspired, broad-angle color reflection

Reproducing the hierarchy of disorder for Morpho-inspired, broad-angle color reflection Supplementary Information for Reproducing the hierarchy of disorder for Morpho-inspired, broad-angle color reflection Bokwang Song 1, Villads Egede Johansen 2,3, Ole Sigmund 3 and Jung H. Shin 4,1,* 1

More information

SUPPLEMENTARY INFORMATION DOI: /NPHOTON

SUPPLEMENTARY INFORMATION DOI: /NPHOTON DOI:.38/NPHOTON.2.85 This supplement has two parts. In part A, we provide the rigorous details of the wavefront correction algorithm and show numerical simulations and experimental data for the cases of

More information

Modeling of Elliptical Air Hole PCF for Lower Dispersion

Modeling of Elliptical Air Hole PCF for Lower Dispersion Advance in Electronic and Electric Engineering. ISSN 2231-1297, Volume 3, Number 4 (2013), pp. 439-446 Research India Publications http://www.ripublication.com/aeee.htm Modeling of Elliptical Air Hole

More information

Light diffraction from colloidal crystals with low dielectric constant modulation: Simulations using single-scattering theory

Light diffraction from colloidal crystals with low dielectric constant modulation: Simulations using single-scattering theory PHYSICAL REVIEW B 77, 23544 28 Light diffraction from colloidal crystals with low dielectric constant modulation: Simulations using single-scattering theory Alexander Tikhonov, Rob D. Coalson, and Sanford

More information

Chapter 7. Widely Tunable Monolithic Laser Diodes

Chapter 7. Widely Tunable Monolithic Laser Diodes Chapter 7 Widely Tunable Monolithic Laser Diodes We have seen in Chapters 4 and 5 that the continuous tuning range λ is limited by λ/λ n/n g, where n is the index change and n g the group index of the

More information

Diffraction. Single-slit diffraction. Diffraction by a circular aperture. Chapter 38. In the forward direction, the intensity is maximal.

Diffraction. Single-slit diffraction. Diffraction by a circular aperture. Chapter 38. In the forward direction, the intensity is maximal. Diffraction Chapter 38 Huygens construction may be used to find the wave observed on the downstream side of an aperture of any shape. Diffraction The interference pattern encodes the shape as a Fourier

More information

Layered media and photonic crystals. Cord Arnold / Anne L Huillier

Layered media and photonic crystals. Cord Arnold / Anne L Huillier Layered media and photonic crystals Cord Arnold / Anne L Huillier Definition A photonic crystal is a periodic arrangement of a dielectric material that exhibits strong interaction with light Variation

More information

Experimental Observation of Invariance of Spectral Degree of Coherence. with Change in Bandwidth of Light

Experimental Observation of Invariance of Spectral Degree of Coherence. with Change in Bandwidth of Light Experimental Observation of Invariance of Spectral Degree of Coherence with Change in Bandwidth of Light Bhaskar Kanseri* and Hem Chandra Kandpal Optical Radiation Standards, National Physical Laboratory,

More information

A Direct Simulation-Based Study of Radiance in a Dynamic Ocean

A Direct Simulation-Based Study of Radiance in a Dynamic Ocean 1 DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. A Direct Simulation-Based Study of Radiance in a Dynamic Ocean LONG-TERM GOALS Dick K.P. Yue Center for Ocean Engineering

More information

Spectrophotometric Methods of Refractive Indices Measurement

Spectrophotometric Methods of Refractive Indices Measurement Application Note Glass, ceramics, optics Spectrophotometric Methods of Refractive Indices Measurement Measuring the refractive index of single crystal optical materials using two methods Authors N.S. Kozlova

More information

Achieve more with light.

Achieve more with light. Achieve more with light. Comprehensive suite of leading photonic design tools. Component Design Multiphysics Component Design Lumerical s highly integrated suite of component design tools is purposebuilt

More information

specular diffuse reflection.

specular diffuse reflection. Lesson 8 Light and Optics The Nature of Light Properties of Light: Reflection Refraction Interference Diffraction Polarization Dispersion and Prisms Total Internal Reflection Huygens s Principle The Nature

More information

Effective Medium Theory, Rough Surfaces, and Moth s Eyes

Effective Medium Theory, Rough Surfaces, and Moth s Eyes Effective Medium Theory, Rough Surfaces, and Moth s Eyes R. Steven Turley, David Allred, Anthony Willey, Joseph Muhlestein, and Zephne Larsen Brigham Young University, Provo, Utah Abstract Optics in the

More information

Advanced modelling of gratings in VirtualLab software. Site Zhang, development engineer Lignt Trans

Advanced modelling of gratings in VirtualLab software. Site Zhang, development engineer Lignt Trans Advanced modelling of gratings in VirtualLab software Site Zhang, development engineer Lignt Trans 1 2 3 4 Content Grating Order Analyzer Rigorous Simulation of Holographic Generated Volume Grating Coupled

More information

E x Direction of Propagation. y B y

E x Direction of Propagation. y B y x E x Direction of Propagation k z z y B y An electromagnetic wave is a travelling wave which has time varying electric and magnetic fields which are perpendicular to each other and the direction of propagation,

More information

WORCESTER POLYTECHNIC INSTITUTE

WORCESTER POLYTECHNIC INSTITUTE WORCESTER POLYTECHNIC INSTITUTE MECHANICAL ENGINEERING DEPARTMENT Optical Metrology and NDT ME-593L, C 2018 Introduction: Wave Optics January 2018 Wave optics: coherence Temporal coherence Review interference

More information

Compact Multilayer Film Structure for Angle Insensitive. Color Filtering

Compact Multilayer Film Structure for Angle Insensitive. Color Filtering 1 Compact Multilayer Film Structure for Angle Insensitive Color Filtering Chenying Yang, Weidong Shen*, Yueguang Zhang, Kan Li, Xu Fang, Xing Zhang, and Xu Liu * E-mail: adongszju@hotmail.com

More information

Chapter 82 Example and Supplementary Problems

Chapter 82 Example and Supplementary Problems Chapter 82 Example and Supplementary Problems Nature of Polarized Light: 1) A partially polarized beam is composed of 2.5W/m 2 of polarized and 4.0W/m 2 of unpolarized light. Determine the degree of polarization

More information

Plane wave in free space Exercise no. 1

Plane wave in free space Exercise no. 1 Plane wave in free space Exercise no. 1 The exercise is focused on numerical modeling of plane wave propagation in ANSYS HFSS. Following aims should be met: 1. A numerical model of a plane wave propagating

More information

Diffuse Optical Tomography, Inverse Problems, and Optimization. Mary Katherine Huffman. Undergraduate Research Fall 2011 Spring 2012

Diffuse Optical Tomography, Inverse Problems, and Optimization. Mary Katherine Huffman. Undergraduate Research Fall 2011 Spring 2012 Diffuse Optical Tomography, Inverse Problems, and Optimization Mary Katherine Huffman Undergraduate Research Fall 11 Spring 12 1. Introduction. This paper discusses research conducted in order to investigate

More information

6-1 LECTURE #6: OPTICAL PROPERTIES OF SOLIDS. Basic question: How do solids interact with light? The answers are linked to:

6-1 LECTURE #6: OPTICAL PROPERTIES OF SOLIDS. Basic question: How do solids interact with light? The answers are linked to: LECTURE #6: OPTICAL PROPERTIES OF SOLIDS Basic question: How do solids interact with light? The answers are linked to: Properties of light inside a solid Mechanisms behind light reflection, absorption

More information

1. Polarization effects in optical spectra of photonic crystals

1. Polarization effects in optical spectra of photonic crystals Speech for JASS 05. April 2005. Samusev A. 1. Polarization effects in optical spectra of photonic crystals Good afternoon. I would like to introduce myself. My name is Anton Samusev. I m a student of Saint

More information

Simulation of Photo-Sensitive Devices with FDTD Method. Copyright 2008 Crosslight Software Inc.

Simulation of Photo-Sensitive Devices with FDTD Method. Copyright 2008 Crosslight Software Inc. Simulation of Photo-Sensitive Devices with FDTD Method Copyright 2008 Crosslight Software Inc. www.crosslight.com What is FDTD method? FDTD=Finite Difference Time Domain FDTD method solves Maxwell s equations

More information

Silicon Avalanche Photodiodes in Dynamic Light Scattering

Silicon Avalanche Photodiodes in Dynamic Light Scattering Silicon Avalanche Photodiodes in Dynamic Light Scattering August 2016 Introduction This application note describes the use of the ID100 single photon counting detector for the measurement of light scattered

More information

Leakage loss and bandgap analysis in air-core photonic bandgap fiber for nonsilica glasses

Leakage loss and bandgap analysis in air-core photonic bandgap fiber for nonsilica glasses Leakage loss and bandgap analysis in air-core photonic bandgap fiber for nonsilica glasses Jonathan Hu and Curtis R. Menyuk University of Maryland Baltimore County, 5200 Westland Blvd., Baltimore, MD 21227

More information

A Graphical User Interface (GUI) for Two-Dimensional Electromagnetic Scattering Problems

A Graphical User Interface (GUI) for Two-Dimensional Electromagnetic Scattering Problems A Graphical User Interface (GUI) for Two-Dimensional Electromagnetic Scattering Problems Veysel Demir vdemir@olemiss.edu Mohamed Al Sharkawy malshark@olemiss.edu Atef Z. Elsherbeni atef@olemiss.edu Abstract

More information

Supporting Information: Highly tunable elastic dielectric metasurface lenses

Supporting Information: Highly tunable elastic dielectric metasurface lenses Supporting Information: Highly tunable elastic dielectric metasurface lenses Seyedeh Mahsa Kamali, Ehsan Arbabi, Amir Arbabi, u Horie, and Andrei Faraon SUPPLEMENTAR NOTE : SAMPLING FREQUENC OF THE PHASE

More information

ratio of the volume under the 2D MTF of a lens to the volume under the 2D MTF of a diffraction limited

ratio of the volume under the 2D MTF of a lens to the volume under the 2D MTF of a diffraction limited SUPPLEMENTARY FIGURES.9 Strehl ratio (a.u.).5 Singlet Doublet 2 Incident angle (degree) 3 Supplementary Figure. Strehl ratio of the singlet and doublet metasurface lenses. Strehl ratio is the ratio of

More information

Accuracy requirements in the mechanical assembly of photonic crystals. Martin Deterre Corey Fucetola Sebastien Uzel

Accuracy requirements in the mechanical assembly of photonic crystals. Martin Deterre Corey Fucetola Sebastien Uzel Accuracy requirements in the mechanical assembly of photonic crystals Martin Deterre Corey Fucetola Sebastien Uzel Agenda Introduction to photonic crystals: theory, background, applications Photonic crystal

More information

High spatial resolution measurement of volume holographic gratings

High spatial resolution measurement of volume holographic gratings High spatial resolution measurement of volume holographic gratings Gregory J. Steckman, Frank Havermeyer Ondax, Inc., 8 E. Duarte Rd., Monrovia, CA, USA 9116 ABSTRACT The conventional approach for measuring

More information

19 Total internal reflection (TIR) and evanescent

19 Total internal reflection (TIR) and evanescent 19 Total internal reflection (TIR) and evanescent waves Consider a TE- or TM-polarized wave (or a superposition) incident on an interface at x =0surface as depicted in the margin at an incidence angle.

More information

WAVELENGTH MANAGEMENT

WAVELENGTH MANAGEMENT BEAM DIAGNOS TICS SPECIAL PRODUCTS OEM DETECTORS THZ DETECTORS PHOTO DETECTORS HIGH POWER SOLUTIONS POWER DETECTORS ENERGY DETECTORS MONITORS Camera Accessories WAVELENGTH MANAGEMENT UV CONVERTERS UV Converters

More information

Supporting Information

Supporting Information Supporting Information Min et al. 10.1073/pnas.1701092114 UV-Cross-Linking Silk Fibroin Using Stilbene Stilbene chromophore and its derivatives have been used as photoreactive building blocks or dopants

More information

ENHANCEMENT OF DIFFUSERS BRDF ACCURACY

ENHANCEMENT OF DIFFUSERS BRDF ACCURACY ENHANCEMENT OF DIFFUSERS BRDF ACCURACY Grégory Bazalgette Courrèges-Lacoste (1), Hedser van Brug (1) and Gerard Otter (1) (1) TNO Science and Industry, Opto-Mechanical Instrumentation Space, P.O.Box 155,

More information

Thin film solar cell simulations with FDTD

Thin film solar cell simulations with FDTD Thin film solar cell simulations with FDTD Matthew Mishrikey, Prof. Ch. Hafner (IFH) Dr. P. Losio (Oerlikon Solar) 5 th Workshop on Numerical Methods for Optical Nano Structures July 7 th, 2009 Problem

More information

Engineered Diffusers Intensity vs Irradiance

Engineered Diffusers Intensity vs Irradiance Engineered Diffusers Intensity vs Irradiance Engineered Diffusers are specified by their divergence angle and intensity profile. The divergence angle usually is given as the width of the intensity distribution

More information

Dynamical Theory of X-Ray Diffraction

Dynamical Theory of X-Ray Diffraction Dynamical Theory of X-Ray Diffraction ANDRE AUTHIER Universite P. et M. Curie, Paris OXFORD UNIVERSITY PRESS Contents I Background and basic results 1 1 Historical developments 3 1.1 Prologue 3 1.2 The

More information

Mu lt i s p e c t r a l

Mu lt i s p e c t r a l Viewing Angle Analyser Revolutionary system for full spectral and polarization measurement in the entire viewing angle EZContrastMS80 & EZContrastMS88 ADVANCED LIGHT ANALYSIS by Field iris Fourier plane

More information

Lecture 7 Notes: 07 / 11. Reflection and refraction

Lecture 7 Notes: 07 / 11. Reflection and refraction Lecture 7 Notes: 07 / 11 Reflection and refraction When an electromagnetic wave, such as light, encounters the surface of a medium, some of it is reflected off the surface, while some crosses the boundary

More information

RSoft Product Applications

RSoft Product Applications RSoft Product Applications Complete design solutions for photonic components, circuits and systems synopsys.com/optical-solutions/rsoft RSoft Photonic Design Software Applications Synopsys RSoft products

More information

Supplementary Figure 1 Optimum transmissive mask design for shaping an incident light to a desired

Supplementary Figure 1 Optimum transmissive mask design for shaping an incident light to a desired Supplementary Figure 1 Optimum transmissive mask design for shaping an incident light to a desired tangential form. (a) The light from the sources and scatterers in the half space (1) passes through the

More information

Diffraction Efficiency

Diffraction Efficiency Diffraction Efficiency Turan Erdogan Gratings are based on diffraction and interference: Diffraction gratings can be understood using the optical principles of diffraction and interference. When light

More information

Analysis of waveguide-splitter-junction in highindex Silicon-On-Insulator waveguides

Analysis of waveguide-splitter-junction in highindex Silicon-On-Insulator waveguides Analysis of waveguide-splitter-junction in highindex Silicon-On-Insulator waveguides Damian Goldring *, Evgeny Alperovich *, Uriel Levy ** and David Mendlovic * * Faculty of Engineering, Tel Aviv University,

More information

Structural color printing based on plasmonic. metasurfaces of perfect light absorption

Structural color printing based on plasmonic. metasurfaces of perfect light absorption Supplementary Information Structural color printing based on plasmonic metasurfaces of perfect light absorption Fei Cheng 1, Jie Gao 1,*, Ting S. Luk 2, and Xiaodong Yang 1,* 1 Department of Mechanical

More information

Supplementary materials of Multispectral imaging using a single bucket detector

Supplementary materials of Multispectral imaging using a single bucket detector Supplementary materials of Multispectral imaging using a single bucket detector Liheng Bian 1, Jinli Suo 1,, Guohai Situ 2, Ziwei Li 1, Jingtao Fan 1, Feng Chen 1 and Qionghai Dai 1 1 Department of Automation,

More information

Strong angular dispersion using higher bands of planar silicon photonic crystals

Strong angular dispersion using higher bands of planar silicon photonic crystals Strong angular dispersion using higher bands of planar silicon photonic crystals Babak Momeni, * Maysamreza Chamanzar, Ehsan Shah Hosseini, Murtaza Askari, Mohammad Soltani, and Ali Adibi School of Electrical

More information

Coupling of surface roughness to the performance of computer-generated holograms

Coupling of surface roughness to the performance of computer-generated holograms Coupling of surface roughness to the performance of computer-generated holograms Ping Zhou* and Jim Burge College of Optical Sciences, University of Arizona, Tucson, Arizona 85721, USA *Corresponding author:

More information

Class 11 Introduction to Surface BRDF and Atmospheric Scattering. Class 12/13 - Measurements of Surface BRDF and Atmospheric Scattering

Class 11 Introduction to Surface BRDF and Atmospheric Scattering. Class 12/13 - Measurements of Surface BRDF and Atmospheric Scattering University of Maryland Baltimore County - UMBC Phys650 - Special Topics in Experimental Atmospheric Physics (Spring 2009) J. V. Martins and M. H. Tabacniks http://userpages.umbc.edu/~martins/phys650/ Class

More information

Multi-GPU Scaling of Direct Sparse Linear System Solver for Finite-Difference Frequency-Domain Photonic Simulation

Multi-GPU Scaling of Direct Sparse Linear System Solver for Finite-Difference Frequency-Domain Photonic Simulation Multi-GPU Scaling of Direct Sparse Linear System Solver for Finite-Difference Frequency-Domain Photonic Simulation 1 Cheng-Han Du* I-Hsin Chung** Weichung Wang* * I n s t i t u t e o f A p p l i e d M

More information

Full-Colour, Computational Ghost Video. Miles Padgett Kelvin Chair of Natural Philosophy

Full-Colour, Computational Ghost Video. Miles Padgett Kelvin Chair of Natural Philosophy Full-Colour, Computational Ghost Video Miles Padgett Kelvin Chair of Natural Philosophy A Quantum Ghost Imager! Generate random photon pairs, illuminate both object and camera SPDC CCD Identical copies

More information

Supporting Information

Supporting Information Supporting Information Chameleon-Inspired Mechanochromic Photonic Films Composed of Nonclose-Packed Colloidal Arrays Gun Ho Lee, Tae Min Choi, Bomi Kim, Sang Hoon Han, Jung Min Lee, and Shin-Hyun Kim *,

More information

Chapter 2: Wave Optics

Chapter 2: Wave Optics Chapter : Wave Optics P-1. We can write a plane wave with the z axis taken in the direction of the wave vector k as u(,) r t Acos tkzarg( A) As c /, T 1/ and k / we can rewrite the plane wave as t z u(,)

More information

Supplementary Figure 1: Schematic of the nanorod-scattered wave along the +z. direction.

Supplementary Figure 1: Schematic of the nanorod-scattered wave along the +z. direction. Supplementary Figure 1: Schematic of the nanorod-scattered wave along the +z direction. Supplementary Figure 2: The nanorod functions as a half-wave plate. The fast axis of the waveplate is parallel to

More information

Mie scattering off plasmonic nanoparticle

Mie scattering off plasmonic nanoparticle Mie scattering off plasmonic nanoparticle Model documentation COMSOL 2009 Version: COMSOL 3.5a1 (build 3.5.0.608) Contents I. Model Overview II. Model Navigator III. Options and settings IV. Geometry modeling

More information

Cambridge University Press Fundamentals of Photonic Crystal Guiding Maksim Skorobogatiy and Jianke Yang Excerpt More information

Cambridge University Press Fundamentals of Photonic Crystal Guiding Maksim Skorobogatiy and Jianke Yang Excerpt More information 1 troduction When thinking about traditional optical materials one invokes a notion of homogeneous media, where imperfections or variations in the material properties are minimal on the length scale of

More information

Deterministic microlens diffuser for Lambertian scatter

Deterministic microlens diffuser for Lambertian scatter Deterministic microlens diffuser for Lambertian scatter Tasso R. M. Sales, Donald J. Schertler, and Stephen Chakmakjian RPC Photonics, Inc. 330 Clay Road, Rochester, New York 14623 Phone: 585-272-2840

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION Supplementary Information S1. Galvo-dithered direct laser writing (GD-DLW) method In order to rectify the asymmetry of the fabrication voxel of the direct laser writing (DLW) method, we have developed

More information

Supporting information for: A highly directional room-temperature single. photon device

Supporting information for: A highly directional room-temperature single. photon device Supporting information for: A highly directional room-temperature single photon device Nitzan Livneh,, Moshe G. Harats,, Daniel Istrati, Hagai S. Eisenberg, and Ronen Rapaport,, Applied Physics Department,

More information

Interference with polarized light

Interference with polarized light Interference with polarized light Summary of the previous lecture (see lecture 3 - slides 12 to 25) With polarized light E 1 et E 2 are complex amplitudes: E 1 + E 2 e iϕ 2 = E 1 2 + E 2 2 + 2 Re(E 1 *

More information

Powerful features (1)

Powerful features (1) HFSS Overview Powerful features (1) Tangential Vector Finite Elements Provides only correct physical solutions with no spurious modes Transfinite Element Method Adaptive Meshing r E = t E γ i i ( x, y,

More information

Imaging and Deconvolution

Imaging and Deconvolution Imaging and Deconvolution Urvashi Rau National Radio Astronomy Observatory, Socorro, NM, USA The van-cittert Zernike theorem Ei E V ij u, v = I l, m e sky j 2 i ul vm dldm 2D Fourier transform : Image

More information

SILICON PHOTONICS WAVEGUIDE AND ITS FIBER INTERCONNECT TECHNOLOGY. Jeong Hwan Song

SILICON PHOTONICS WAVEGUIDE AND ITS FIBER INTERCONNECT TECHNOLOGY. Jeong Hwan Song SILICON PHOTONICS WAVEGUIDE AND ITS FIBER INTERCONNECT TECHNOLOGY Jeong Hwan Song CONTENTS Introduction of light waveguides Principals Types / materials Si photonics Interface design between optical fiber

More information

APSS Apollo Application Note on Photonic Crystal Fibers (PCFs)

APSS Apollo Application Note on Photonic Crystal Fibers (PCFs) APSS Apollo Application Note on Photonic Crystal Fibers (PCFs) Design and simulation APN-APSS-PCF Apollo Inc. 1057 Main Street West Hamilton, Ontario L8S 1B7 Canada Tel: (905)-524-3030 Fax: (905)-524-3050

More information

Spectrographs. C. A. Griffith, Class Notes, PTYS 521, 2016 Not for distribution.

Spectrographs. C. A. Griffith, Class Notes, PTYS 521, 2016 Not for distribution. Spectrographs C A Griffith, Class Notes, PTYS 521, 2016 Not for distribution 1 Spectrographs and their characteristics A spectrograph is an instrument that disperses light into a frequency spectrum, which

More information

Chap. 4. Jones Matrix Method

Chap. 4. Jones Matrix Method Chap. 4. Jones Matrix Method 4.1. Jones Matrix Formulation - For an incident light with a polarization state described by the Jones vector - Decompose the light into a linear combination of the "fast"

More information

Visible-frequency dielectric metasurfaces for multi-wavelength achromatic and highly-dispersive holograms

Visible-frequency dielectric metasurfaces for multi-wavelength achromatic and highly-dispersive holograms Supporting Materials Visible-frequency dielectric metasurfaces for multi-wavelength achromatic and highly-dispersive holograms Bo Wang,, Fengliang Dong,, Qi-Tong Li, Dong Yang, Chengwei Sun, Jianjun Chen,,

More information

dq dt I = Irradiance or Light Intensity is Flux Φ per area A (W/m 2 ) Φ =

dq dt I = Irradiance or Light Intensity is Flux Φ per area A (W/m 2 ) Φ = Radiometry (From Intro to Optics, Pedrotti -4) Radiometry is measurement of Emag radiation (light) Consider a small spherical source Total energy radiating from the body over some time is Q total Radiant

More information

IMPLEMENTATION OF ANALYTICAL (MATLAB) AND NUMERICAL (HFSS) SOLUTIONS ADVANCED ELECTROMAGNETIC THEORY SOHAIB SAADAT AFRIDI HAMMAD BUTT ZUNNURAIN AHMAD

IMPLEMENTATION OF ANALYTICAL (MATLAB) AND NUMERICAL (HFSS) SOLUTIONS ADVANCED ELECTROMAGNETIC THEORY SOHAIB SAADAT AFRIDI HAMMAD BUTT ZUNNURAIN AHMAD STUDY OF SCATTERING & RESULTANT RADIATION PATTERN: INFINITE LINE CURRENT SOURCE POSITIONED HORIZONTALLY OVER A PERFECTLY CONDUCTING INFINITE GROUND PLANE IMPLEMENTATION OF ANALYTICAL (MATLAB) AND NUMERICAL

More information

WAVELENGTH MANAGEMENT

WAVELENGTH MANAGEMENT Camera Accessories WAVELENGTH MANAGEMENT UV CONVERTERS UV Converters take advantage of a phenomenon called fluorescence to extend the performance range of the Beamage beam profiling camera to ultraviolet

More information

Supporting Information for Azimuthal Polarization Filtering for Accurate, Precise, and Robust Single-Molecule Localization Microscopy

Supporting Information for Azimuthal Polarization Filtering for Accurate, Precise, and Robust Single-Molecule Localization Microscopy Nano Letters Supporting Information for Azimuthal Polarization Filtering for Accurate, Precise, and Robust Single-Molecule Localization Microscopy Matthew D. Lew, and W. E. Moerner *, Departments of Chemistry

More information

Optics Vac Work MT 2008

Optics Vac Work MT 2008 Optics Vac Work MT 2008 1. Explain what is meant by the Fraunhofer condition for diffraction. [4] An aperture lies in the plane z = 0 and has amplitude transmission function T(y) independent of x. It is

More information

2011 Optical Science & Engineering PhD Qualifying Examination Optical Sciences Track: Advanced Optics Time allowed: 90 minutes

2011 Optical Science & Engineering PhD Qualifying Examination Optical Sciences Track: Advanced Optics Time allowed: 90 minutes 2011 Optical Science & Engineering PhD Qualifying Examination Optical Sciences Track: Advanced Optics Time allowed: 90 minutes Answer all four questions. All questions count equally. 3(a) A linearly polarized

More information

X-Ray fluorescence and Raman spectroscopy

X-Ray fluorescence and Raman spectroscopy X-Ray fluorescence and Raman spectroscopy Advanced physics laboratory (nd part) 4CFU Catalini Letizia, De Angelis Giulia Vittoria, Piselli Verdiana Abstract In this paper we report about two different

More information

Michelson Interferometer

Michelson Interferometer Michelson Interferometer The Michelson interferometer uses the interference of two reflected waves The third, beamsplitting, mirror is partially reflecting ( half silvered, except it s a thin Aluminum

More information

Optical Photon Processes

Optical Photon Processes Optical Photon Processes GEANT4 is an effective and comprehensive tool capable of realistically modeling the optics of scintillation and Cerenkov detectors and their associated light guides. This is founded

More information

Physics 202, Lecture 23

Physics 202, Lecture 23 Physics 202, Lecture 23 Today s Topics Lights and Laws of Geometric Optics Nature of Light Reflection and Refraction Law of Reflection Law of Refraction Index of Reflection, Snell s Law Total Internal

More information

Instruction manual for T3DS calculator software. Analyzer for terahertz spectra and imaging data. Release 2.4

Instruction manual for T3DS calculator software. Analyzer for terahertz spectra and imaging data. Release 2.4 Instruction manual for T3DS calculator software Release 2.4 T3DS calculator v2.4 16/02/2018 www.batop.de1 Table of contents 0. Preliminary remarks...3 1. Analyzing material properties...4 1.1 Loading data...4

More information

An Introduction to the Finite Difference Time Domain (FDTD) Method & EMPIRE XCcel

An Introduction to the Finite Difference Time Domain (FDTD) Method & EMPIRE XCcel An Introduction to the Finite Difference Time Domain (FDTD) Method & EMPIRE XCcel Simulation Model definition for FDTD DUT Port Simulation Box Graded Mesh six Boundary Conditions 1 FDTD Basics: Field components

More information

9. Polarizers. Index of. Coefficient of Material Wavelength ( ) Brewster angle refraction (n)

9. Polarizers. Index of. Coefficient of Material Wavelength ( ) Brewster angle refraction (n) 9. Polarizers All polarized light is to some degree elliptical in nature. Basic states of polarization like linear and circular are actually special cases of elliptically polarized light which is defined

More information

Modelling of standard and specialty fibre-based systems using finite element methods

Modelling of standard and specialty fibre-based systems using finite element methods Modelling of standard and specialty fibre-based systems using finite element methods Natascia Castagna* a, Jacques Morel a, Luc Testa b, Sven Burger c,d a Federal Institute of Metrology METAS, Lindenweg

More information

Outline. Darren Wang ADS Momentum P2

Outline. Darren Wang ADS Momentum P2 Outline Momentum Basics: Microstrip Meander Line Momentum RF Mode: RFIC Launch Designing with Momentum: Via Fed Patch Antenna Momentum Techniques: 3dB Splitter Look-alike Momentum Optimization: 3 GHz Band

More information

A Single Grating-lens Focusing Two Orthogonally Polarized Beams in Opposite Direction

A Single Grating-lens Focusing Two Orthogonally Polarized Beams in Opposite Direction , pp.41-45 http://dx.doi.org/10.14257/astl.2016.140.08 A Single Grating-lens Focusing Two Orthogonally Polarized Beams in Opposite Direction Seung Dae Lee 1 1* Dept. of Electronic Engineering, Namseoul

More information

Grating Spectrometer GRATING.TEX KB

Grating Spectrometer GRATING.TEX KB UNIVERSITÄTK OSNABRÜCK 1 Grating Spectrometer GRATING.TEX KB 20020119 KLAUS BETZLER 1,FACHBEREICH PHYSIK, UNIVERSITÄT OSNABRÜCK This short lecture note recalls some of the well-known properties of spectrometers.

More information

Optical simulations within and beyond the paraxial limit

Optical simulations within and beyond the paraxial limit Optical simulations within and beyond the paraxial limit Daniel Brown, Charlotte Bond and Andreas Freise University of Birmingham 1 Simulating realistic optics We need to know how to accurately calculate

More information

Fig The light rays that exit the prism enter longitudinally into an astronomical telescope adjusted for infinite distance.

Fig The light rays that exit the prism enter longitudinally into an astronomical telescope adjusted for infinite distance. Romanian Master of Physics 07 Problem I Reflection and refraction of light A. An interesting prism The main section of a glass prism, situated in air n '.00, has the form of a rhomb with. A thin yellow

More information

ECE 595, Section 10 Numerical Simulations Lecture 33: Introduction to Finite- Difference Time-Domain Simulations. Prof. Peter Bermel April 3, 2013

ECE 595, Section 10 Numerical Simulations Lecture 33: Introduction to Finite- Difference Time-Domain Simulations. Prof. Peter Bermel April 3, 2013 ECE 595, Section 10 Numerical Simulations Lecture 33: Introduction to Finite- Difference Time-Domain Simulations Prof. Peter Bermel April 3, 2013 Recap from Monday Numerical ODE solvers Initial value problems

More information

A Direct Simulation-Based Study of Radiance in a Dynamic Ocean

A Direct Simulation-Based Study of Radiance in a Dynamic Ocean A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Dick K.P. Yue Center for Ocean Engineering Massachusetts Institute of Technology Room 5-321, 77 Massachusetts Ave, Cambridge, MA 02139 phone:

More information

CRISM (Compact Reconnaissance Imaging Spectrometer for Mars) on MRO. Calibration Upgrade, version 2 to 3

CRISM (Compact Reconnaissance Imaging Spectrometer for Mars) on MRO. Calibration Upgrade, version 2 to 3 CRISM (Compact Reconnaissance Imaging Spectrometer for Mars) on MRO Calibration Upgrade, version 2 to 3 Dave Humm Applied Physics Laboratory, Laurel, MD 20723 18 March 2012 1 Calibration Overview 2 Simplified

More information

Magnetics. Introduction to Filtering using the Fourier Transform. Chuck Connor, Laura Connor. Potential Fields. Magnetics.

Magnetics. Introduction to Filtering using the Fourier Transform. Chuck Connor, Laura Connor. Potential Fields. Magnetics. Introduction to using the Chuck Connor, Laura Connor Potential Fields for this week Nm Sm Schematic Earth dipolar magnetic field. The field lines placed in the page plane are drawn as thick lines, those

More information

Design of Electromagnetic Test Sites

Design of Electromagnetic Test Sites Sensor and Simulation Notes Note 533 3 August 2008 Design of Electromagnetic Test Sites Carl E. Baum University of New Mexico Department of Electrical and Computer Engineering Albuquerque New Mexico 87131

More information

ECEG105/ECEU646 Optics for Engineers Course Notes Part 5: Polarization

ECEG105/ECEU646 Optics for Engineers Course Notes Part 5: Polarization ECEG105/ECEU646 Optics for Engineers Course Notes Part 5: Polarization Prof. Charles A. DiMarzio Northeastern University Fall 2008 Sept 2008 11270-05-1 Wave Nature of Light Failure of Raytracing Zero-λ

More information

Lecture PowerPoints. Chapter 24 Physics: Principles with Applications, 7 th edition Giancoli

Lecture PowerPoints. Chapter 24 Physics: Principles with Applications, 7 th edition Giancoli Lecture PowerPoints Chapter 24 Physics: Principles with Applications, 7 th edition Giancoli This work is protected by United States copyright laws and is provided solely for the use of instructors in teaching

More information

Progress of the Thomson Scattering Experiment on HSX

Progress of the Thomson Scattering Experiment on HSX Progress of the Thomson Scattering Experiment on HSX K. Zhai, F.S.B. Anderson, D.T. Anderson HSX Plasma Laboratory, UW-Madison Bill Mason PSL, UW-Madison, The Thomson scattering system being constructed

More information

Chapter 35 &36 Physical Optics

Chapter 35 &36 Physical Optics Chapter 35 &36 Physical Optics Physical Optics Phase Difference & Coherence Thin Film Interference 2-Slit Interference Single Slit Interference Diffraction Patterns Diffraction Grating Diffraction & Resolution

More information

LASer Cavity Analysis and Design

LASer Cavity Analysis and Design The unique combination of simulation tools for LASer Cavity Analysis and Design During the last 15 years LASCAD has become industry-leading so ware for LASer Cavity Analysis and Design. The feedback from

More information

Instruction Manual for Modelling of Trap Detectors. Helsinki University of Technology. Metrology Research Institute. Instruction Manual for

Instruction Manual for Modelling of Trap Detectors. Helsinki University of Technology. Metrology Research Institute. Instruction Manual for Page 1 (10) Helsinki University of Technology Metrology Research Institute Instruction Manual for Modelling Version: 2.2 Date of Issue: December 21, 2005 Page 2 (10) Table of Contents 1. INTRODUCTION 3

More information

Iterative procedure for in-situ EUV optical testing with an incoherent source

Iterative procedure for in-situ EUV optical testing with an incoherent source APS/123-QED Iterative procedure for in-situ EUV optical testing with an incoherent source Ryan Miyakawa and Patrick Naulleau Lawrence Berkeley National Laboratory, Berkeley, CA 94720 Avideh Zakhor Dept.

More information

CMSC427 Shading Intro. Credit: slides from Dr. Zwicker

CMSC427 Shading Intro. Credit: slides from Dr. Zwicker CMSC427 Shading Intro Credit: slides from Dr. Zwicker 2 Today Shading Introduction Radiometry & BRDFs Local shading models Light sources Shading strategies Shading Compute interaction of light with surfaces

More information