Normal vector method for the RCWA with automated vector field generation
|
|
- Coleen Roberts
- 5 years ago
- Views:
Transcription
1 P. Götz, T. Schuster, K. Frenner, S. Rafler, W. Osten Normal vector method for the RCWA with automated vector field generation Stuttgart, May 2011 Institute of Applied Optics, University of Stuttgart, University of Stuttgart, Pfaffenwaldring 9, Stuttgart/ Germany Abstract Early formulations of the RCWA yield, implicated by the erroneous application of factorization rules to discrete Fourier transformations, poor convergence in certain cases. An explanation for this finding and an approach to overcome the problem for crossed gratings was first given by Li [J. Opt. Soc. Am. A 13, 1870 (1996) and 14, 2758 (1997)]. A further improvement was achieved by Schuster et al. [J. Opt. Soc. Am. A 24, 2880 (2007)], using a structure dependent normal vector (NV) field. While it is trivial to create those NV fields for simple geometrical shapes, to our knowledge an appropriate algorithm for arbitrary shapes does not exist, yet. In this work we present such an algorithm. Keywords Gratings Diffraction theory This paper was published in Optics Express and is made available as an electronic reprint with the permission of OSA. The paper can be found at the following URL on the OSA website: Systematic or multiple reproduction or distribution to multiple locations via electronic or other means is prohibited and is subject to penalties under law. Preprint Series Issue No Stuttgart Research Centre for Simulation Technology (SRC SimTech) SimTech Cluster of Excellence Pfaffenwaldring 7a Stuttgart publications@simtech.uni-stuttgart.de
2 Normal vector method for the RCWA with automated vector field generation Peter Götz, Thomas Schuster, Karsten Frenner, Stephan Rafler, Wolfgang Osten Institut für Technische Optik, Universität Stuttgart, Germany Corresponding author: Abstract: Early formulations of the RCWA yield, implicated by the erroneous application of factorization rules to discrete Fourier transformations, poor convergence in certain cases. An explanation for this finding and an approach to overcome the problem for crossed gratings was first given by Li [J. Opt. Soc. Am. A 13, 1870 (1996) and 14, 2758 (1997)]. A further improvement was achieved by Schuster et al. [J. Opt. Soc. Am. A 24, 2880 (2007)], using a structure dependent normal vector (NV) field. While it is trivial to create those NV fields for simple geometrical shapes, to our knowledge an appropriate algorithm for arbitrary shapes does not exist, yet. In this work we present such an algorithm Optical Society of America OCIS codes: ( ) Gratings; ( ) Diffraction theory References and links 1. M. G. Moharam and T. K. Gaylord, Rigorous coupled-wave analysis of planar-grating diffraction, J. Opt. Soc. Am. 71, (1981). 2. L. Li and C. W. Haggans, Convergence of the coupled-wave method for metallic lamellar diffraction gratings, J. Opt. Soc. Am. A 10, (1993). 3. P. Lalanne and M. Morris, Highly improved convergence of the coupled-wave method for TM polarization, J. Opt. Soc. Am. A 13, (1996). 4. G. Granet and B. Guizal, Efficient implemenation of the coupled-wave method for metallic lamellar grating in TM polarization, J. Opt. Soc. Am. A 13, (1996). 5. L. Li, Use of Fourier series in the analysis of discontinuous periodic structures, J. Opt. Soc. Am. A 13, (1996). 6. P. Lalanne, Improved formulation of the coupled-wave method for two-dimensional gratings, J. Opt. Soc. Am. A 14, (1997). 7. L. Li, New formulation of the Fourier modal method for crossed surface-relief gratings, J. Opt. Soc. Am. A 14, (1997). 8. T. Schuster, J. Ruoff, N. Kerwien, S. Rafler, and W. Osten, Normal vector method for convergence improvement using the RCWA for crossed gratings, J. Opt. Soc. Am. A 24, (2007). 9. M. Nevière and E. Popov, Light Propagation in Periodic Media Differential Theory and Design (Marcel Dekker, Inc., New York, 2003). 10. E. Popov and M. Nevière, Maxwell equations in Fourier space: fast-converging formulation for diffraction by arbitrary shaped, periodic, anisotropic media, J. Opt. Soc. Am. A 18, (2001). 11. D. Shepard, A two-dimensional interpolation function for irregularly-spaced data, pp (ACM National Conference, 1968). 12. R. Franke, Scattered Data Interpolation: Tests of Some Methods, Math. Comput. 38, (1982). 13. M. G. Moharam, T. K. Gaylord, E. B. Grann and D. A. Pommet, Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings, J. Opt. Soc. Am. A 12, (1995). 14. R. Bräuer and O. Bryngdahl, Electromagnetic diffraction analysis of two-dimensional gratings Opt. Commun. 100, 1 5 (1993). 15. P. Götz, Simulation von Lichtbeugung mit der Normalenvektormethode, Diploma thesis, Universität Stuttgart, Germany (2008). (C) 2008 OSA 27 October 2008 / Vol. 16, No. 22 / OPTICS EXPRESS 17295
3 1. Introduction For the numerical simulation of grating diffraction rigorous coupled-wave analysis (RCWA) is a commonly used method. After its publication in 1981 by Moharam and Gaylord [1] it turned out that its results for metallic gratings in TM polarization suffered from slow convergence [2] due to numerical problems. While empirical solutions came from Lalanne [3] as well as from Guizal and Granet [4], Li finally gave a theoretical explanation, introducing his factorization rules [5]. Since then different approaches were pursued to apply these rules to crossed gratings in order to improve convergence [6, 7]. One of the most recent improvements for crossed gratings was presented by Schuster et al. [8]. They adopted the idea of the Fast Fourier Factorization by Nevière and Popov [9, 10] and developed a normal vector (NV) method which is an optimized solution for the correct application of Li s factorization rules to crossed gratings. This method, however, requires a structure dependent NV field. The requirements for this vector field will be explained in the next Section (for a detailed explanation see [8]). Afterward, Section 3 and 4 will explain how the algorithm works. Results of the algorithm will be shown in Section Background The RCWA is based on a time-harmonic version of Maxwell s Equations to solve the diffraction problem: E = iωµ 0 H (1) H = iωε 0 εe (2) To solve these equations the RCWA uses pseudo Fourier expansions as an ansatz for E, H and ε to transform this system of PDEs to a system of ODEs, which, in turn, can be reduced to an eigenvalue problem. Important for the following considerations is the calculation of the product ε 0 εe. This product, the dielectric displacement D = ε 0 εe, (3) transforms into a convolution in Fourier space. At material boundaries the normal component of Eq. (3) in contrast to the tangential component is a product of complementary discontinuous functions ε and E, as is generally known from electrodynamics. In truncated Fourier space, as it is used in numerical calculations, this product cannot be calculated by the usual convolution [D ]=ε 0 ε [E ], (4) as this leads to unnecessary oscillations in the dielectric displacement D as shown in Fig. 4 of [5]. In Eq. (4) single brackets refer to vertical vectors containing all Fourier coefficients of the specified function and ε denotes the Toeplitz matrix of the Fourier coefficients of ε with the entries ε jk = ε j k. While Lalanne gave strong empirical evidence [3] that another rule must be used for the normal component of the electrical field, Li later gave a theoretical explanation [5] why the so-called inverse rule [D ]=ε 0 1 ε 1 [E ], (5) must be used for complementary discontinuous functions instead. Consequently the normal and tangential components must be treated differently. While Li used a method of zigzag lines [7] and got an improved convergence, an even further improvement in convergence could be achieved by the normal vector method of Schuster et al. [8]. (C) 2008 OSA 27 October 2008 / Vol. 16, No. 22 / OPTICS EXPRESS 17296
4 This method basically projects the electrical field onto normal vectors to get the tangential and normal components of the electrical field. Until now, however, there was no way to calculate those normal vector fields for completely arbitrary structures. A solution for such an algorithm will be given in the following. 3. Automatic generation of normal vector fields The algorithm is based on the fact that normal vectors are only defined on the boundaries. Except for their orientation (inwards or outwards) vectors are unambiguously defined here. Since the products Nx 2, N2 y and N xn y are transformed into Fourier space, as can be seen in Eq. (8) of [8], the criterion for the vectors in the remaining space is smooth transitions between NVs. An interpolation fulfills this requirement. Generating the vector field is consequently done in two steps: 1. Calculating the normal vectors on the boundaries, called boundary vectors 2. Interpolating the vectors in the remaining space These two steps will be explained further in the following two subsections Calculating the boundary vectors A bitmap serves as a description for the crossed grating structure, where different materials are indicated by different color indices. Determining the boundary vectors is done by calculating the gradient of this bitmap. The pixelation yields a staircase-like modulation of the boundary vectors (Fig. 1a). A blurring of the bitmap by a Gaussian filter smooths the pixelation and (a) (b) (c) (d) Fig. 1. (a) Staircase effect of the gradient (b) Gradient after blurring the bitmap (c) Result of combining (a) and (b) (d) Optional change of the preferred direction hence avoids the modulation (Fig. 1b). However, the filter also smooths the distinct jumps in the bitmap and hence results in a wide band of normal vectors, which is undesired. The gradient of the original image Fig. 1a is only 0 in a narrow band around the boundary and thus can be used as a mask for Fig. 1b. The result of this procedure is a narrow band of proper boundary vectors as can be seen in Fig. 1c. After the boundary vectors are calculated the orientation of the vectors can be given a preferred orientation as it is shown in Fig. 1d. In this way a degree of freedom is introduced which causes differences in both the appearance of the NV fields and the end results of the RCWA. The influence of a preferred orientation is investigated in Section Calculating the remaining vectors Interpolation is done using inverse distance weighting [11, 12]. This is an interpolation method suited for irregularly spaced nodes which can be implemented to work very fast as will be (C) 2008 OSA 27 October 2008 / Vol. 16, No. 22 / OPTICS EXPRESS 17297
5 shown in the next Section. Using this interpolation method, a vector N at position x can be calculated as and be normalized as N (x)= n b i=1 N i N(x)= N (x) N (x). x x i 2 (6) N i denotes the boundary vectors, x i their positions and n b their count. Eq. (6) must be evaluated for every x except for those where boundary vectors already exist. 4. Accelerating the algorithm Using this implementation of the algorithm one realizes that the time complexity is O(R 3 ), where R R is the resolution. With a resolution of e.g this leads to very long computation durations. To gain an advantage for the NV method against other methods the complexity can be reduced to O(R 2 ) using a progressive refinement algorithm. This algorithm makes sure that the number of considered boundary vectors in Eq. (6) becomes smaller with every iteration. The idea is, that the influence of distant vectors is negligible. Hence, they can be omitted completely. Fig. 2 shows how the progressive refinement algorithm works. Only vectors from previous Fig. 2. Progressive refinement algorithm. P is the example vector to be calculated. steps within the gray square region around a point P are included in the computation. Vectors to be computed in the current iteration are drawn in red and the exemplarily considered point P is marked by a red ring. The pre-computed boundary vectors are drawn in black. Vectors computed in earlier refinement steps are called raster vectors and are drawn in green. The set of raster vectors gradually becomes finer until the desired resolution is reached. To get periodic boundary conditions, in the first iteration (Fig. 2a) the region of considered NVs covers the structure itself and its four neighboring structures. Only the boundary vectors (C) 2008 OSA 27 October 2008 / Vol. 16, No. 22 / OPTICS EXPRESS 17298
6 are used as an input for the interpolation in this step, since no raster vectors exist, yet. In the second iteration (Fig. 2b) the region size is identical to the size of the unit cell. Now the red vectors are computed from the boundary vectors in the surrounding square and the raster vector (only a single one at this point). In the third iteration (Fig. 2c) the region shrinks to a quarter of the structure size, leaving only five boundary vectors to consider plus four raster vectors from the iteration before in the corners. Another quartering is done in the last iteration leaving only one boundary vector to consider plus the four raster vectors in the corners. (Fig. 2d). Of course, in practice the sampling of the initial boundary vectors and the final raster size are chosen to be finer, which means a much higher number of iteration steps. As mentioned this algorithm features a complexity of only O(R 2 ) which gives the needed performance boost to calculate NV fields in a fraction of the time the RCWA takes. 5. Results In the following we present the simulation results of two example structures. The NV-method presented here is compared to the zigzag method of Li [7] as well as to an implementation which ignores Li s rules which we refer to as the original formulation. Our implementation is based on [13] which is a more modern formulation of the authors of [1] and at the same time one of the most appreciated publications on RCWA. Although 2D periodic gratings are not explicitly mentioned there, its consideration in RCWA is straightforward as pointed out first by Bräuer and Bryngdahl [14]. For both examples the simulation parameters from Table 1 are taken. The first example is a Table 1. Simulation parameters Parameter Value Wavelength 500 nm Side length of structure 1000nm Thickness of grating 50 nm Refractive index of substrate, superstrate 1.5, 1 Refractive index of white, gray, light gray region 1, i, 1.5 Direction of incidence perpendicular Polarization in y direction C-shaped structure. Fig. 3 shows the two variants with and without preferred direction of the NVs. Plotting the corresponding diffraction efficiencies as a function of the retained Fourier modes, commonly referred to as truncation order (Fig. 4), one can see, that the NV method has a considerably faster convergence than both the original formulation and the zigzag method. It also shows that both variants, although having a completely different look, yield almost the same convergence rate. The second example is a structure which is built out of two intersecting circles. In contrast to the C-shaped structure which is hand drafted it features a certain regularity. However, in order to make sure that the method succeeds not only due to symmetry, odd values for the circle parameters were chosen. For convenience these data are reproduced in Table 2 with the origin being in the upper left corner of the unit cell. NV fields for this structure are shown in Fig. 5. Again, the NV method shows better convergence than the other methods (Fig. 6). Furthermore it is obvious that none of the two NV variants are superior to the other. A general rule of thumb cannot be made for which of both variants must be used. It is also impossible to make rules by arguing with the number of singularities in the NV field or with the distance of singularities from boundaries. Instead, there is very strong evidence that the proper (C) 2008 OSA 27 October 2008 / Vol. 16, No. 22 / OPTICS EXPRESS 17299
7 (a) (b) Fig. 3. NV field for the C-shaped structure (a) without preferred direction, (b) with preferred direction to the right. (Note that both vector fields have in fact a higher resolution than shown here) Diffraction Efficiency (Reflexion) Original formulation 0.1 Li s zigzag method NV method without preferred direction NV method with preferred direction Truncation Order Fig. 4. Convergence plot for the C-shaped structure (a) (b) Fig. 5. NV field for intersecting circles (a) without preferred direction, (b) with preferred direction (C) 2008 OSA 27 October 2008 / Vol. 16, No. 22 / OPTICS EXPRESS 17300
8 Diffraction Efficiency (Reflexion) Original formulation Li s zigzag method NV method without preferred direction NV method with preferred direction Truncation Order Fig. 6. Convergence plot for the structure with intersecting circles Table 2. Data of the intersecting circle structure Parameter Value center of left circle [x, y]=[402.4nm, 503.4nm] radius of left circle 260.7nm center of right circle [x, y]=[640.9nm, 514.2nm] radius of right circle 266.9nm representation of N 2 x, N2 y and N xn y as truncated Fourier series plays an important role. To be precise, it is important that the deviations of the truncated Fourier series of N 2 x, N 2 y and N x N y from the exact values are as small as possible at the material boundary. A detailed investigation of this matter is presented in [15] but omitted here for the sake of brevity. 6. Conclusion With the presented algorithm it is now possible to calculate NV fields for arbitrary structures, e.g. concave shapes or structures with several arbitrarily connected materials. While in this work we presented two variants of the algorithm, both yielding similar results, in practice one of these variants can be chosen. In conjunction with this algorithm the NV method with its excellent convergence is a very attractive alternative to Li s zigzag method for crossed gratings. (C) 2008 OSA 27 October 2008 / Vol. 16, No. 22 / OPTICS EXPRESS 17301
annual report 2005 / 2006 INSTITUT FÜR TECHNISCHE OPTIK UNIVERSITÄT STUTTGART
annual report 2005 / 2006 INSTITUT FÜR TECHNISCHE OPTIK UNIVERSITÄT STUTTGART 33 High Resolution Metrology and Simulation Simulation of diffraction at large structures using the fieldstitching method Scatterometry
More informationRicWaA v RicWaA User s Guide. Introduction. Features. Study with an example.
RicWaA User s Guide RicWaA v1.0 2010 Introduction. RicWaA is a MATLAB based Rigorous coupled-wave Analysis (RCWA) code. Built with the object-oriented programming of MATLAB, RicWaA provides a very simple
More informationCoupling 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 informationDiffraction Gratings as Anti Reflective Coatings Noah Gilbert. University of Arizona ngilbert .arizona.edu Phone: (520)
Diffraction Gratings as Anti Reflective Coatings Noah Gilbert University of Arizona Email: ngilbertemail.arizona.edu Phone: (520)304 4864 Abstract: Diffraction gratings with sub wavelength spatial frequencies
More informationAnalysis of stacked rotated gratings
Brigham Young University BYU ScholarsArchive All Faculty Publications 2007-03-10 Analysis of stacked rotated gratings Gregory P. Nordin nordin@byu.edu P. C. Deguzman See next page for additional authors
More informationA 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 informationTilted Wave Interferometer Improved Measurement Uncertainty. Germany. ABSTRACT
URN (Paper): urn:nbn:de:gbv:ilm1-2014iwk-118:5 58 th ILMENAU SCIENTIFIC COLLOQUIUM Technische Universität Ilmenau, 08 12 September 2014 URN: urn:nbn:gbv:ilm1-2014iwk:3 Tilted Wave Interferometer Improved
More informationInfluence of the Aspect Ratio of Tabular Grains on the Light Scattering
000 International Symposium on Silver 000 Halide International TechnologySymposium on Silver Halide Technology Copyright 000, IS&T Influence of the Aspect Ratio of Tabular Grains on the Light Scattering
More informationDigital Image Processing. Prof. P. K. Biswas. Department of Electronic & Electrical Communication Engineering
Digital Image Processing Prof. P. K. Biswas Department of Electronic & Electrical Communication Engineering Indian Institute of Technology, Kharagpur Lecture - 21 Image Enhancement Frequency Domain Processing
More informationA new Eulerian computational method for the propagation of short acoustic and electromagnetic pulses
A new Eulerian computational method for the propagation of short acoustic and electromagnetic pulses J. Steinhoff, M. Fan & L. Wang. Abstract A new method is described to compute short acoustic or electromagnetic
More informationVector-based synthesis of finite aperiodic. to as subwavelength diffractive optical elements
Prather et al. Vol. 15, No. 6/June 1998/J. Opt. Soc. Am. A 1599 Vector-based synthesis of finite aperiodic subwavelength diffractive optical elements Dennis W. Prather and Joseph N. Mait U.S. Army Research
More informationStratified Volume Diffractive Optical Elements as High Efficiency Gratings
Brigham Young University BYU ScholarsArchive All Faculty Publications 1999-05-01 Stratified Volume Diffractive Optical Elements as High Efficiency Gratings Gregory P. Nordin nordin@byu.edu D. M. Chambers
More informationCS 565 Computer Vision. Nazar Khan PUCIT Lectures 15 and 16: Optic Flow
CS 565 Computer Vision Nazar Khan PUCIT Lectures 15 and 16: Optic Flow Introduction Basic Problem given: image sequence f(x, y, z), where (x, y) specifies the location and z denotes time wanted: displacement
More informationspecular 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 informationFilters. Advanced and Special Topics: Filters. Filters
Filters Advanced and Special Topics: Filters Dr. Edmund Lam Department of Electrical and Electronic Engineering The University of Hong Kong ELEC4245: Digital Image Processing (Second Semester, 2016 17)
More informationCrossed Grating. Diffraction gratings design and analysis by Modal and Chandezon methods. User Manual
Crossed Grating Diffraction gratings design and analysis by Modal and Chandezon methods 2018 User Manual Introduction... - 2 - C-method restrictions... - 3 - Fourier modal method restrictions... - 3 -
More informationELECTROMAGNETIC diffraction by perfectly conducting
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 47, NO. 11, NOVEMBER 1999 1697 Oblique Scattering by a Pair of Conducting Half Planes: TM Case Jong-Won Yu and Noh-Hoon Myung Abstract An exact series
More information2.5D Forward Solver To Model Scattering Of Long Dielectric Cylinders In An Active Millimeter Wave Imaging System
2.5D Forward Solver To Model Scattering Of Long Dielectric Cylinders In An Active Millimeter Wave Imaging System S. Van den Bulcke and A. Franchois Department of Information Technology Ghent University,
More informationCalculation of the diffraction efficiency on concave gratings based on Fresnel Kirchhoff s diffraction formula
Calculation of the diffraction efficiency on concave gratings based on Fresnel Kirchhoff s diffraction formula Yuanshen Huang, 1, * Ting Li, 1 Banglian Xu, 1 Ruijin Hong, 1 Chunxian Tao, 1 Jinzhong Ling,
More informationExperimental reconstruction of a highly reflecting fiber Bragg grating by using spectral regularization and inverse scattering
3284 J. Opt. Soc. Am. A/ Vol. 24, No. 10/ October 2007 Rosenthal et al. Experimental reconstruction of a highly reflecting fiber Bragg grating by using spectral regularization and inverse scattering Amir
More informationScattering reduction effect with overcoated rough surfaces: theory and experiment
Scattering reduction effect with overcoated rough surfaces: theory and experiment H. Giovannini and C. Amra We show that a scattering reduction effect is obtained by coating a rough surface with an antireflection
More informationEnhanced optical absorptance of metals using interferometric femtosecond ablation
Enhanced optical absorptance of metals using interferometric femtosecond ablation K. Paivasaari, J. J. J. Kaakkunen, M. Kuittinen and T. Jaaskelainen Department of Physics and Mathematics, University of
More informationPhysical & Electromagnetic Optics: Diffraction Gratings
31/05/2018 Physical & Electromagnetic Optics: Diffraction Gratings Optical Engineering Prof. Elias N. Glytsis School of Electrical & Computer Engineering National Technical University of Athens Multiple
More informationSupplementary 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 informationSelf-formation, Development and Reproduction of the Artificial System
Solid State Phenomena Vols. 97-98 (4) pp 77-84 (4) Trans Tech Publications, Switzerland Journal doi:.48/www.scientific.net/ssp.97-98.77 Citation (to be inserted by the publisher) Copyright by Trans Tech
More informationSystematic design process for slanted graing couplers,
Brigham Young University BYU ScholarsArchive All Faculty Publications 2006-08-20 Systematic design process for slanted graing couplers, Gregory P. Nordin nordin@byu.edu J. Jiang See next page for additional
More informationContinued Fraction Absorbing Boundary Conditions for Transient Elastic Wave Propagation Modeling
Continued Fraction Absorbing Boundary Conditions for Transient Elastic Wave Propagation Modeling Md Anwar Zahid and Murthy N. Guddati 1 Summary This paper presents a novel absorbing boundary condition
More informationBlood vessel tracking in retinal images
Y. Jiang, A. Bainbridge-Smith, A. B. Morris, Blood Vessel Tracking in Retinal Images, Proceedings of Image and Vision Computing New Zealand 2007, pp. 126 131, Hamilton, New Zealand, December 2007. Blood
More informationA 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 informationModifications of detour phase computer-generated holograms
Modifications of detour phase computer-generated holograms Uriel Levy, Emanuel Marom, and David Mendlovic The detour phase method for the design of computer-generated holograms can be modified to achieve
More informationAn Intuitive Explanation of Fourier Theory
An Intuitive Explanation of Fourier Theory Steven Lehar slehar@cns.bu.edu Fourier theory is pretty complicated mathematically. But there are some beautifully simple holistic concepts behind Fourier theory
More informationModule 1 Lecture Notes 2. Optimization Problem and Model Formulation
Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization
More informationAdvanced 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 informationAngular sensitivities of volume gratings for. of a substrate-mode optical interconnect
Angular sensitivities of volume gratings for substrate-mode optical interconnects Shun-Der Wu, homas K. Gaylord, Elias N. Glytsis, and Yu-Ming Wu he angular sensitivities of slanted volume gratings (VGs)
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 WRI C225 Lecture 04 130131 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Histogram Equalization Image Filtering Linear
More informationLecture 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 informationApplication of Tatian s Method to Slanted-Edge MTF Measurement
Application of s Method to Slanted-Edge MTF Measurement Peter D. Burns Eastman Kodak Company, Rochester, NY USA 465-95 ABSTRACT The 33 method for the measurement of the spatial frequency response () of
More informationDOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS
DOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS RAYLEIGH-SOMMERFELD DIFFRACTION RECTANGULAR APERTURES Ian Cooper School of Physics, University of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECTORY FOR MATLAB
More informationOptimised corrections for finite-difference modelling in two dimensions
Optimized corrections for 2D FD modelling Optimised corrections for finite-difference modelling in two dimensions Peter M. Manning and Gary F. Margrave ABSTRACT Finite-difference two-dimensional correction
More informationBiometrics Technology: Image Processing & Pattern Recognition (by Dr. Dickson Tong)
Biometrics Technology: Image Processing & Pattern Recognition (by Dr. Dickson Tong) References: [1] http://homepages.inf.ed.ac.uk/rbf/hipr2/index.htm [2] http://www.cs.wisc.edu/~dyer/cs540/notes/vision.html
More informationCHAPTER 6 PERCEPTUAL ORGANIZATION BASED ON TEMPORAL DYNAMICS
CHAPTER 6 PERCEPTUAL ORGANIZATION BASED ON TEMPORAL DYNAMICS This chapter presents a computational model for perceptual organization. A figure-ground segregation network is proposed based on a novel boundary
More informationDiffraction 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 informationAdvanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude
Advanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude A. Migukin *, V. atkovnik and J. Astola Department of Signal Processing, Tampere University of Technology,
More informationOptimization of anisotropically etched silicon surface-relief gratings for substrate-mode optical interconnects
Optimization of anisotropically etched silicon surface-relief gratings for substrate-mode optical interconnects Shun-Der Wu, Thomas K. Gaylord, Jonathan S. Maikisch, and Elias N. Glytsis The optimum profiles
More informationAnalysis of the Gaussian Beam on a Corrugated Dielectric Interface
(An ISO 3297: 27 Certified Organization) Vol. 3, Issue 9, September 214 Analysis of the Gaussian Beam on a Corrugated Dielectric Interface Mohamed B. El_Mashade 1, Adel Shaaban 2 Department of Electrical
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationEfficient computation of source magnetic scalar potential
Adv. Radio Sci., 4, 59 63, 2006 Author(s) 2006. This work is licensed under a Creative Commons License. Advances in Radio Science Efficient computation of source magnetic scalar potential W. Hafla, A.
More informationPolarizing properties of embedded symmetric trilayer stacks under conditions of frustrated total internal reflection
University of New Orleans ScholarWorks@UNO Electrical Engineering Faculty Publications Department of Electrical Engineering 3-1-2006 Polarizing properties of embedded symmetric trilayer stacks under conditions
More informationAnalysis of Slanted Gratings for Lightguide Coupling
Analysis of Slanted Gratings for Lightguide Coupling Abstract Slanted gratings are commonly used for coupling light into optical lightguides due to their high efficiency in a certain diffraction order.
More informationOptical Design with Zemax
Optical Design with Zemax Lecture 7: Optimization I 2012-12-11 Herbert Gross Winter term 2012 www.iap.uni-jena.de Time schedule 2 1 16.10. Introduction Introduction, Zemax interface, menues, file handling,
More informationMotion Estimation. There are three main types (or applications) of motion estimation:
Members: D91922016 朱威達 R93922010 林聖凱 R93922044 謝俊瑋 Motion Estimation There are three main types (or applications) of motion estimation: Parametric motion (image alignment) The main idea of parametric motion
More informationExperiments with Edge Detection using One-dimensional Surface Fitting
Experiments with Edge Detection using One-dimensional Surface Fitting Gabor Terei, Jorge Luis Nunes e Silva Brito The Ohio State University, Department of Geodetic Science and Surveying 1958 Neil Avenue,
More informationSpatial Enhancement Definition
Spatial Enhancement Nickolas Faust The Electro- Optics, Environment, and Materials Laboratory Georgia Tech Research Institute Georgia Institute of Technology Definition Spectral enhancement relies on changing
More informationModeling Diffractive Effects due to Micro-lens Arrays on Liquid Crystal Panels in Projectors
Modeling Diffractive Effects due to Micro-lens Arrays on Liquid Crystal Panels in Projectors Synopsys, Inc., 400 Executive Boulevard, Suite 100, Ossining, NY, USA 10562 ABSTRACT The components in optical
More informationA Diagonal Split-cell Model for the High-order Symplectic FDTD Scheme
PIERS ONLINE, VOL. 2, NO. 6, 2006 715 A Diagonal Split-cell Model for the High-order Symplectic FDTD Scheme Wei Sha, Xianliang Wu, and Mingsheng Chen Key Laboratory of Intelligent Computing & Signal Processing
More informationComputer-originated planar holographic optical elements
Computer-originated planar holographic optical elements Silviu Reinhorn, Yaakov Amitai, and Albert A. Friesem We present novel, to our knowledge, methods for the analytical design and recording of planar
More informationSELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND
Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana
More informationDesign of paraxial diffractive elements with the CAD system DIGIOPT
Design of paraxial diffractive elements with the CAD system DIGIOPT Harald Aagedal, Thomas Beth, Heiko Schwarzer, Stephan Teiwes Institute of Algorithms and Cognitive Systems, University of Karlsruhe Am
More informationEffective 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 informationWhite-light interference microscopy: minimization of spurious diffraction effects by geometric phase-shifting
White-light interference microscopy: minimization of spurious diffraction effects by geometric phase-shifting Maitreyee Roy 1, *, Joanna Schmit 2 and Parameswaran Hariharan 1 1 School of Physics, University
More informationUltrasonic Multi-Skip Tomography for Pipe Inspection
18 th World Conference on Non destructive Testing, 16-2 April 212, Durban, South Africa Ultrasonic Multi-Skip Tomography for Pipe Inspection Arno VOLKER 1, Rik VOS 1 Alan HUNTER 1 1 TNO, Stieltjesweg 1,
More informationResearch Collection. Localisation of Acoustic Emission in Reinforced Concrete using Heterogeneous Velocity Models. Conference Paper.
Research Collection Conference Paper Localisation of Acoustic Emission in Reinforced Concrete using Heterogeneous Velocity Models Author(s): Gollob, Stephan; Vogel, Thomas Publication Date: 2014 Permanent
More informationCHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT
CHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT 9.1 Introduction In the previous chapters the inpainting was considered as an iterative algorithm. PDE based method uses iterations to converge
More informationHOUGH TRANSFORM CS 6350 C V
HOUGH TRANSFORM CS 6350 C V HOUGH TRANSFORM The problem: Given a set of points in 2-D, find if a sub-set of these points, fall on a LINE. Hough Transform One powerful global method for detecting edges
More informationDiffraction and Interference of Plane Light Waves
1 Diffraction and Interference of Plane Light Waves Introduction In this experiment you will become familiar with diffraction patterns created when a beam of light scatters from objects placed in its path.
More informationTexture. Frequency Descriptors. Frequency Descriptors. Frequency Descriptors. Frequency Descriptors. Frequency Descriptors
Texture The most fundamental question is: How can we measure texture, i.e., how can we quantitatively distinguish between different textures? Of course it is not enough to look at the intensity of individual
More informationNew Edge-Enhanced Error Diffusion Algorithm Based on the Error Sum Criterion
New Edge-Enhanced Error Diffusion Algorithm Based on the Error Sum Criterion Jae Ho Kim* Tae Il Chung Hyung Soon Kim* Kyung Sik Son* Pusan National University Image and Communication Laboratory San 3,
More informationLet s review the four equations we now call Maxwell s equations. (Gauss s law for magnetism) (Faraday s law)
Electromagnetic Waves Let s review the four equations we now call Maxwell s equations. E da= B d A= Q encl ε E B d l = ( ic + ε ) encl (Gauss s law) (Gauss s law for magnetism) dφ µ (Ampere s law) dt dφ
More informationSIMULATION OF ARTIFICIAL SYSTEMS BEHAVIOR IN PARAMETRIC EIGHT-DIMENSIONAL SPACE
78 Proceedings of the 4 th International Conference on Informatics and Information Technology SIMULATION OF ARTIFICIAL SYSTEMS BEHAVIOR IN PARAMETRIC EIGHT-DIMENSIONAL SPACE D. Ulbikiene, J. Ulbikas, K.
More informationNEAR-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 informationModeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics
Modeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics Anthony L. Bonomo *1 and Marcia J. Isakson 1 1 Applied
More informationPROCEEDINGS OF SPIE. Simulation of lateral color for a hybrid refractive-diffractive eyepiece by field tracing methods
PROCEEDINGS OF SPIE SPIEDigitalLibrary.org/conference-proceedings-of-spie Simulation of lateral color for a hybrid refractive-diffractive eyepiece by field tracing methods D. Batte, M. Kuhn, F. Wyrowski
More informationSupplementary 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 informationModeling Skills Thermal Analysis J.E. Akin, Rice University
Introduction Modeling Skills Thermal Analysis J.E. Akin, Rice University Most finite element analysis tasks involve utilizing commercial software, for which you do not have the source code. Thus, you need
More informationRobust Ring Detection In Phase Correlation Surfaces
Griffith Research Online https://research-repository.griffith.edu.au Robust Ring Detection In Phase Correlation Surfaces Author Gonzalez, Ruben Published 2013 Conference Title 2013 International Conference
More informationHigh Frequency Wave Scattering
High Frequency Wave Scattering University of Reading March 21st, 2006 - Scattering theory What is the effect of obstacles or inhomogeneities on an incident wave? - Scattering theory What is the effect
More informationMesh Based Interpolative Coding (MBIC)
Mesh Based Interpolative Coding (MBIC) Eckhart Baum, Joachim Speidel Institut für Nachrichtenübertragung, University of Stuttgart An alternative method to H.6 encoding of moving images at bit rates below
More informationPhilip E. Plantz. Application Note. SL-AN-08 Revision C. Provided By: Microtrac, Inc. Particle Size Measuring Instrumentation
A Conceptual, Non-Mathematical Explanation on the Use of Refractive Index in Laser Particle Size Measurement (Understanding the concept of refractive index and Mie Scattering in Microtrac Instruments and
More informationKwabena Arthur. at the. June2017. Certified by: George Barbastathis Professor Thesis Supervisor
Simulation of X-Ray Phase Imaging on Integrated Circuits by Kwabena Arthur Submitted to the Department of Mechanical Engineering in Partial Fulfillment of the Requirements for the Degree of Bachelor of
More informationRecent advances in Metamodel of Optimal Prognosis. Lectures. Thomas Most & Johannes Will
Lectures Recent advances in Metamodel of Optimal Prognosis Thomas Most & Johannes Will presented at the Weimar Optimization and Stochastic Days 2010 Source: www.dynardo.de/en/library Recent advances in
More informationNull test for a highly paraboloidal mirror
Null test for a highly paraboloidal mirror Taehee Kim, James H. Burge, Yunwoo Lee, and Sungsik Kim A circular null computer-generated hologram CGH was used to test a highly paraboloidal mirror diameter,
More informationComputation of Three-Dimensional Electromagnetic Fields for an Augmented Reality Environment
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Computation of Three-Dimensional Electromagnetic Fields for an Augmented Reality Environment André Buchau 1 * and Wolfgang M. Rucker
More informationaxis, and wavelength tuning is achieved by translating the grating along a scan direction parallel to the x
Exponential-Grating Monochromator Kenneth C. Johnson, October 0, 08 Abstract A monochromator optical design is described, which comprises a grazing-incidence reflection and two grazing-incidence mirrors,
More informationModel-Based Segmentation of Impression Marks
Model-Based Segmentation of Impression Marks Christoph Brein Institut für Mess- und Regelungstechnik, Universität Karlsruhe (TH), D-76128 Karlsruhe, Germany ABSTRACT Impression marks are commonly found
More informationAPPLYING EXTRAPOLATION AND INTERPOLATION METHODS TO MEASURED AND SIMULATED HRTF DATA USING SPHERICAL HARMONIC DECOMPOSITION.
APPLYING EXTRAPOLATION AND INTERPOLATION METHODS TO MEASURED AND SIMULATED HRTF DATA USING SPHERICAL HARMONIC DECOMPOSITION Martin Pollow Institute of Technical Acoustics RWTH Aachen University Neustraße
More informationRANS Based Analysis of Roll Damping Moments at Bilge Keels
RANS Based Analysis of Roll Damping Moments at Bilge Keels Florian Kluwe (kluwe@tu-harburg.de), Daniel Schmode, Gerhard Jensen Introduction The simulation of ship motions in seaways gets increasing relevance
More information1. 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 informationIntroduction to Digital Image Processing
Fall 2005 Image Enhancement in the Spatial Domain: Histograms, Arithmetic/Logic Operators, Basics of Spatial Filtering, Smoothing Spatial Filters Tuesday, February 7 2006, Overview (1): Before We Begin
More informationFree-Form Shape Optimization using CAD Models
Free-Form Shape Optimization using CAD Models D. Baumgärtner 1, M. Breitenberger 1, K.-U. Bletzinger 1 1 Lehrstuhl für Statik, Technische Universität München (TUM), Arcisstraße 21, D-80333 München 1 Motivation
More informationCh. 22 Properties of Light HW# 1, 5, 7, 9, 11, 15, 19, 22, 29, 37, 38
Ch. 22 Properties of Light HW# 1, 5, 7, 9, 11, 15, 19, 22, 29, 37, 38 Brief History of the Nature of Light Up until 19 th century, light was modeled as a stream of particles. Newton was a proponent of
More informationFast Algorithm for Matrix-Vector Multiply of Asymmetric Multilevel Block-Toeplitz Matrices
" Fast Algorithm for Matrix-Vector Multiply of Asymmetric Multilevel Block-Toeplitz Matrices B. E. Barrowes, F. L. Teixeira, and J. A. Kong Research Laboratory of Electronics, MIT, Cambridge, MA 02139-4307
More informationHomework 4: Clustering, Recommenders, Dim. Reduction, ML and Graph Mining (due November 19 th, 2014, 2:30pm, in class hard-copy please)
Virginia Tech. Computer Science CS 5614 (Big) Data Management Systems Fall 2014, Prakash Homework 4: Clustering, Recommenders, Dim. Reduction, ML and Graph Mining (due November 19 th, 2014, 2:30pm, in
More informationHolographic Elements in Solar Concentrator and Collection Systems
Holographic Elements in Solar Concentrator and Collection Systems Raymond K. Kostuk,2, Jose Castro, Brian Myer 2, Deming Zhang and Glenn Rosenberg 3 Electrical and Computer Engineering, Department University
More informationLocally Weighted Least Squares Regression for Image Denoising, Reconstruction and Up-sampling
Locally Weighted Least Squares Regression for Image Denoising, Reconstruction and Up-sampling Moritz Baecher May 15, 29 1 Introduction Edge-preserving smoothing and super-resolution are classic and important
More informationFast marching methods
1 Fast marching methods Lecture 3 Alexander & Michael Bronstein tosca.cs.technion.ac.il/book Numerical geometry of non-rigid shapes Stanford University, Winter 2009 Metric discretization 2 Approach I:
More informationBragg properties of efficient surface relief gratings in the resonance domain
Optics Communications 235 (24) 261 267 www.elsevier.com/locate/optcom Bragg properties of efficient surface relief gratings in the resonance domain M.A. Golub *, A.A. Friesem, L. Eisen Department of Physics
More informationHigh-resolution 3D profilometry with binary phase-shifting methods
High-resolution 3D profilometry with binary phase-shifting methods Song Zhang Department of Mechanical Engineering, Iowa State University, Ames, Iowa 511, USA (song@iastate.edu) Received 11 November 21;
More informationf. (5.3.1) So, the higher frequency means the lower wavelength. Visible part of light spectrum covers the range of wavelengths from
Lecture 5-3 Interference and Diffraction of EM Waves During our previous lectures we have been talking about electromagnetic (EM) waves. As we know, harmonic waves of any type represent periodic process
More informationx' = c 1 x + c 2 y + c 3 xy + c 4 y' = c 5 x + c 6 y + c 7 xy + c 8
1. Explain about gray level interpolation. The distortion correction equations yield non integer values for x' and y'. Because the distorted image g is digital, its pixel values are defined only at integer
More informationVivekananda. Collegee of Engineering & Technology. Question and Answers on 10CS762 /10IS762 UNIT- 5 : IMAGE ENHANCEMENT.
Vivekananda Collegee of Engineering & Technology Question and Answers on 10CS762 /10IS762 UNIT- 5 : IMAGE ENHANCEMENT Dept. Prepared by Harivinod N Assistant Professor, of Computer Science and Engineering,
More information