Lecture 04: HFK Propagation Physical Optics II (Optical Sciences 330) (Updated: Friday, April 29, 2005, 8:05 PM) W.J. Dallas

Similar documents
DOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS

Fresnel Diffraction. monchromatic light source

Chapter 3 Optical Systems with Annular Pupils

The condition for maximum intensity by the transmitted light in a plane parallel air film is. For an air film, μ = 1. (2-1)

ANTENNAS. Vector and Scalar Potentials. Maxwell's Equations. D = εe. For a linear, homogeneous, isotropic medium µ and ε are contant.

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

Qualifying Examination Electricity and Magnetism Solutions January 12, 2006

11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.

Multipole Radiation. February 29, The electromagnetic field of an isolated, oscillating source

3. Electromagnetic Waves II

you of a spring. The potential energy for a spring is given by the parabola U( x)

TheWaveandHelmholtzEquations

Question Bank. Section A. is skew-hermitian matrix. is diagonalizable. (, ) , Evaluate (, ) 12 about = 1 and = Find, if

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s

Do not turn over until you are told to do so by the Invigilator.

2 Governing Equations

PES 3950/PHYS 6950: Homework Assignment 6

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

Chapter 16. Fraunhofer Diffraction

$ i. !((( dv vol. Physics 8.02 Quiz One Equations Fall q 1 q 2 r 2 C = 2 C! V 2 = Q 2 2C F = 4!" or. r ˆ = points from source q to observer

THE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2

1 Spherical multipole moments

Stress, Cauchy s equation and the Navier-Stokes equations

Lecture 8 - Gauss s Law

Math 209 Assignment 9 Solutions

Geometry of the homogeneous and isotropic spaces

1 Fundamental Solutions to the Wave Equation

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

Scattering in Three Dimensions

Lecture 2 Date:

Review Notes on Maxwell's Equations

Right-handed screw dislocation in an isotropic solid

Flux. Area Vector. Flux of Electric Field. Gauss s Law

The geometric construction of Ewald sphere and Bragg condition:

Hopefully Helpful Hints for Gauss s Law

Antennas & Propagation

I. CONSTRUCTION OF THE GREEN S FUNCTION

B. Spherical Wave Propagation

Light Time Delay and Apparent Position

Black Body Radiation and Radiometric Parameters:

Math 259 Winter Handout 6: In-class Review for the Cumulative Final Exam

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

sinγ(h y > ) exp(iωt iqx)dωdq


University of Illinois at Chicago Department of Physics. Electricity & Magnetism Qualifying Examination

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II September 15, 2012 Prof. Alan Guth PROBLEM SET 2

Motion in One Dimension

Electromagnetism Physics 15b

Module 9: Electromagnetic Waves-I Lecture 9: Electromagnetic Waves-I

Inverse Square Law and Polarization

Boundary Layers and Singular Perturbation Lectures 16 and 17 Boundary Layers and Singular Perturbation. x% 0 Ω0æ + Kx% 1 Ω0æ + ` : 0. (9.

COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS

Stanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012

A New Approach to General Relativity

4. Electrodynamic fields

Fields and Waves I Spring 2005 Homework 8. Due: 3 May 2005

EM Boundary Value Problems

As is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.

Phys101 Lectures 30, 31. Wave Motion

Double-angle & power-reduction identities. Elementary Functions. Double-angle & power-reduction identities. Double-angle & power-reduction identities

Green s Identities and Green s Functions

Basic Interference and. Classes of of Interferometers

Graphs of Sine and Cosine Functions

Magnetic Field. Conference 6. Physics 102 General Physics II

f(k) e p 2 (k) e iax 2 (k a) r 2 e a x a a 2 + k 2 e a2 x 1 2 H(x) ik p (k) 4 r 3 cos Y 2 = 4

Chapter 22 The Electric Field II: Continuous Charge Distributions

Analytical calculation of the power dissipated in the LHC liner. Stefano De Santis - LBNL and Andrea Mostacci - CERN

, the tangent line is an approximation of the curve (and easier to deal with than the curve).

Homework # 3 Solution Key

Supporting Information Wedge Dyakonov Waves and Dyakonov Plasmons in Topological Insulator Bi 2 Se 3 Probed by Electron Beams

17.1 Electric Potential Energy. Equipotential Lines. PE = energy associated with an arrangement of objects that exert forces on each other

Electrostatic Potential

Objectives: After finishing this unit you should be able to:

Related Rates - the Basics

Waves and Polarization in General

Question 1: The dipole

MONTE CARLO SIMULATION OF FLUID FLOW

Euclidean Figures and Solids without Incircles or Inspheres

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

Vector d is a linear vector function of vector d when the following relationships hold:

PY208 Matter & Interactions Final Exam S2005

Gauss Law. Physics 231 Lecture 2-1

Physics 2B Chapter 22 Notes - Magnetic Field Spring 2018

Chem 453/544 Fall /08/03. Exam #1 Solutions

ME 210 Applied Mathematics for Mechanical Engineers

Chapter 2: Basic Physics and Math Supplements

Physics 235 Chapter 5. Chapter 5 Gravitation

MATH 417 Homework 3 Instructor: D. Cabrera Due June 30. sin θ v x = v r cos θ v θ r. (b) Then use the Cauchy-Riemann equations in polar coordinates

TELE4652 Mobile and Satellite Communications

Tutorial on Strehl ratio, wavefront power series expansion, Zernike polynomials expansion in small aberrated optical systems By Sheng Yuan

Physics 181. Assignment 4

Review: Electrostatics and Magnetostatics

Physics 2212 GH Quiz #2 Solutions Spring 2016

Voltage ( = Electric Potential )

4. Kruskal Coordinates and Penrose Diagrams.

MULTILAYER PERCEPTRONS

Math 2263 Solutions for Spring 2003 Final Exam

CHAPTER 25 ELECTRIC POTENTIAL

221B Lecture Notes Scattering Theory I

5. Plane Electromagnetic Waves

Transcription:

C:\Dallas\0_Couses\0_OpSci_330\0 Lectue Notes\04 HfkPopagation.doc: Page of 9 Lectue 04: HFK Popagation Physical Optics II (Optical Sciences 330) (Updated: Fiday, Apil 9, 005, 8:05 PM) W.J. Dallas The Pinhole Wave The popagation of optical waves in fee space is govened by the Helmholt equation which is also known simply as the wave equation. Fo now, we will concentate on a paticulaly useful kind of popagation, that of waves taveling between paallel planes. Fo this discussion, we will use a heuistic, but highly accuate agument. Plane-Paallel Popagation Figue illustates the popagation, in fee space, between two planes of the wave coming fom a pinhole in a flat opaque sceen. We ll call it a pinhole wave. It is paticulaly impotant because the patten in the second plane epesents the PSF of the point in the fist plane. Moving the point lateally just moves the PSF, i.e., wave popagation in fee space is an LSI system. We can completely descibe the system by knowing only the PSF. The pictue in the second plane of Figue is a stylied epesentation of the eal pat of the pinhole wave. It is lateally tuncated and thesholded at a value of eo in that plane. A moe accuate illustation will come afte we have looked into the mathematical epesentation of the pinhole wave. y X =0 = 0 Figue : A pinhole wave popagating fom the pinhole to a second plane. Why the pinhole wave is not the same as the spheical wave The wave fom a point standing fee in space is a spheical wave. A spheical wave has the complex eik amplitude us ( ) =. Fo example, it dies-off lateally in the x-diection as in the plane of the souce: x it cannot be a pinhole wave; that wave is eo except at the oigin.

C:\Dallas\0_Couses\0_OpSci_330\0 Lectue Notes\04 HfkPopagation.doc: Page of 9 How we modify a spheical wave to suit ou needs We can imagine a pinhole-wave candidate by looking at the spheical wave shown in Figue. In that figue we have also de-emphasied the left-popagating wave. Notice that the wave-fonts intesect the x-y plane at ight angles except at the oigin. This behavio suggests that we use u s. x Taking the deivative, with espect to, gives Figue : Spheical wave intesecting the x-y plane ik ik x + y + ik e e e uph( x, y, ) = = = ik π π x + y + π Ou pinhole wave is pecisely the Huygens wavelet. It is closely elated to the Rayleigh-Sommefeld diffaction equation. Nomaliing the pinhole wave We want to be sue that, at the plane = 0 the pinhole wave is a -D delta function. We need to have that ph (,,0) = δ δ u x y x y

C:\Dallas\0_Couses\0_OpSci_330\0 Lectue Notes\04 HfkPopagation.doc: Page 3 of 9 Checking the value at the oigin diectly is not possible because of the delta function s infinite value at that point. Instead, we use the fact that δ ( x) δ ( y) dxdy = of. π Fee-space popagation as an LSI system If we shift the pinhole in the fist plane lateally we get ph (,,0) = δ ( ) δ ( ) u x x y y x x y y 0 0 0 0. It is this integal that gives us the leading facto In the second plane, this shift has the effect of shifting the pinhole wave in that plane, giving ph (,, ) u x x y y 0 0 0 We have a linea shift-invaiant system with a point spead function (,, ) uph x y 0. This fact gives us yet anothe name fo the pinhole wave: the point spead function of fee space. Fom ou knowledge of LSI systems, we know that the popagation is descibed by the convolution u x, y, = u x, y,0 u x, y,. Witten out, this convolution is 0 ph 0 u x y u x y u x x y y dx dy (,, ) = (,,0) (,, ) 0 0 0 ph 0 0 0 0 0 The wave field eveywhee to the ight of the plane = 0 is detemined by the value of the wave in the plane = 0. The exact Huygens-Fesnel-Kichhoff (HFK) fomula Let s look at two abitay planes, not necessaily having one plane at the oigin. In the fom of a convolution, the HFK fomula fo popagation fom the plane = 0 to the plane = is u x, y, = u x, y, u x, y,. Witten out using the explicit wave fomula, this equation is 0 ph 0 Paaxial popagation ( 0) u( x, y, ) = u( x0, y0, 0) ik π i e ( x x ) + ( y y ) + ( ) 0 0 0 ( ) + ( ) + ( ) ik x x0 y y0 0 ( x x ) + ( y y ) + ( ) 0 0 0 dx dy This fomula is highly accuate, but a little cumbesome. It is a common and useful pactice to examine popagation when the obsevation egion is bounded to small distances fom the -axis. The appoximations that go into that calculation lead to paaxial optics. The appoximations ae done in fou steps stating with the pinhole wave fomula. 0 0

C:\Dallas\0_Couses\0_OpSci_330\0 Lectue Notes\04 HfkPopagation.doc: Page 4 of 9 ik e uph( x, y, ) = ik π

Appoximations C:\Dallas\0_Couses\0_OpSci_330\0 Lectue Notes\04 HfkPopagation.doc: Page 5 of 9. π i ik = We note that so that λ is negligible fo λ π the x-y plane of seveal wavelengths. This tem becomesik ik. >>, i.e., a distance away fom. The tem is the angle that the position vecto to a point on the pinhole wave makes with the -axis. y x Figue 3: The angle to the optical axis Fo small angles, i.e. nea the -axis, = cosθ ik ik e e 3. Nea the -axis, so that. ik 4. The appoximation fo the phase, k, must be done caefully because e is apidly vaying. The paabolic appoximation of is done using the Taylo expansion. We begin with the expession x + = x + y + = + y then use the fist-ode Taylo expansion

/ + ε + ε, ε << to get So that: x+ y πi ik ik ik λ C:\Dallas\0_Couses\0_OpSci_330\0 Lectue Notes\04 HfkPopagation.doc: Page 6 of 9 x + y x + y x + y + fo = +. x y iπ + λ e e e = e e. Taking all of the steps togethe we have The paaxial pinhole-wave fomula Tidying up the equation gives πi λ uph ( x, y, ) = e e iλ ( + ) iπ x y λ u ph π i e ( x, y, ) π λ πi λ x + y iπ λ e. The paaxial HFK fomula fo popagation Inseting this expession into the popagation elation we get the paaxial Huygens-Fesnel-Kichhoff (HFK) popagation equation, π ( 0) λ i iπ ( x x0) + ( y y0) e λ( 0) u x y u x y e dx dy,, =,, 0 0 0 0 0 iλ 0 Thin optical elements We will be woking with complex-amplitude tansmitting thin elements. These elements will be chaacteied by a complex-amplitude tansmittance. This tansmittance descibes the element's action on a geneal wave. If we suppose that the element is located in the plane " = 0", then an incident wave slightly to the left of the element will be conveted to an exiting wave slightly to the ight of the element in the manne: ( 0 + ) ( 0 ) u x, y, = u x, y, t x, y

t(x,y) C:\Dallas\0_Couses\0_OpSci_330\0 Lectue Notes\04 HfkPopagation.doc: Page 7 of 9 y x u(x, y, = 0 - ) u(x, y, = 0 + ) Figue 4: A thin optical element The Neutal-Density Filte The thin neutal-density filte is chaacteied by a constant-amplitude attenuation 0 0 t x,y = A whee A Apetues These ae clea holes in opaque sceens. Thee ae seveal common types. Rectangula apetue x x t ( x, y) = ect ect x0 y0 Finite-width infinite-length slit x t ( x, y) = ect x0 Infinitesimal-width infinite-length slit t x,y = δ x Infinitesimal-width finite-length slit y t ( x, y) = δ ( x) ect y0 Cicula apetue

t x,y = cyl d0 C:\Dallas\0_Couses\0_OpSci_330\0 Lectue Notes\04 HfkPopagation.doc: Page 8 of 9 Pinhole t x,y = δ x δ y Double pinhole x0 x0 t x,y = δ x + δ x + δ y Infinite double slit x0 x0 t ( x,y) = δ x + δ x + Finite double slit x0 x0 y t ( x, y) = δ x + δ x + ect y0

The Plane-Paallel Plate C:\Dallas\0_Couses\0_OpSci_330\0 Lectue Notes\04 HfkPopagation.doc: Page 9 of 9 The thin plane-paallel plate is chaacteied by a constant phase shift ove its entie suface. The expession fo the plate's complex amplitude tansmittance is i 0 t( x,y ) = e φ Whee φ0 is a constant. The Pism A pism is a thin element which has a phase shift that vaies linealy with position. The expession fo the pism's complex amplitude tansmittance is i 0 x + 0 t x,y = e π ξ η y whee ξ 0 and η 0 ae constants. The Lens The lens is a focusing element. It essentially convets one pinhole wave into a second pinhole wave. It is consideed to have, in the paabolic appoximation, a quadatic phase facto. The complex-amplitude tansmittance of a lens located on-axis is iπ ( x +y ) λ f t x, y = e ± powe To be examined in subsequent lectues Diffaction gatings Fesnel one plates Chomatic filtes Polaies