OSA's Digital Library

Applied Optics

Applied Optics


  • Editor: Joseph N. Mait
  • Vol. 48, Iss. 32 — Nov. 10, 2009
  • pp: 6132–6142

Digital simulation of scalar optical diffraction: revisiting chirp function sampling criteria and consequences

David G. Voelz and Michael C. Roggemann  »View Author Affiliations

Applied Optics, Vol. 48, Issue 32, pp. 6132-6142 (2009)

View Full Text Article

Enhanced HTML    Acrobat PDF (719 KB)

Browse Journals / Lookup Meetings

Browse by Journal and Year


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools



Accurate simulation of scalar optical diffraction requires consideration of the sampling requirement for the phase chirp function that appears in the Fresnel diffraction expression. We describe three sampling regimes for FFT-based propagation approaches: ideally sampled, oversampled, and undersampled. Ideal sampling, where the chirp and its FFT both have values that match analytic chirp expressions, usually provides the most accurate results but can be difficult to realize in practical simulations. Under- or oversampling leads to a reduction in the available source plane support size, the available source bandwidth, or the available observation support size, depending on the approach and simulation scenario. We discuss three Fresnel propagation approaches: the impulse response/transfer function (angular spectrum) method, the single FFT (direct) method, and the two-step method. With illustrations and simulation examples we show the form of the sampled chirp functions and their discrete transforms, common relationships between the three methods under ideal sampling conditions, and define conditions and consequences to be considered when using nonideal sampling. The analysis is extended to describe the sampling limitations for the more exact Rayleigh–Sommerfeld diffraction solution.

© 2009 Optical Society of America

OCIS Codes
(050.1940) Diffraction and gratings : Diffraction
(070.2580) Fourier optics and signal processing : Paraxial wave optics
(050.1755) Diffraction and gratings : Computational electromagnetic methods
(070.7345) Fourier optics and signal processing : Wave propagation

ToC Category:
Fourier Optics and Signal Processing

Original Manuscript: June 19, 2009
Manuscript Accepted: August 17, 2009
Published: November 2, 2009

David G. Voelz and Michael C. Roggemann, "Digital simulation of scalar optical diffraction: revisiting chirp function sampling criteria and consequences," Appl. Opt. 48, 6132-6142 (2009)

Sort:  Author  |  Year  |  Journal  |  Reset  


  1. J. W. Goodman, Introduction to Fourier Optics, 3rd ed.(Roberts and Company, 2005).
  2. A. Papoulis, “Pulse compression, fiber communications and diffraction, a unified approach,” J. Opt. Soc. Am. A 11, 3-13(1994). [CrossRef]
  3. F. Gori, “Fresnel transform and sampling theorem,” Opt. Commun. 39, 293-297 (1981). [CrossRef]
  4. D. Mas, J. Garcia, C. Ferreira, L. M. Bernardo, and F. Marinho, “Fast algorithms for free-space diffraction patterns calculations,” Opt. Commun. 164, 233-245 (1999). [CrossRef]
  5. L. Onural, “Sampling of the diffraction field,” Appl. Opt. 39, 5929-5935 (2000). [CrossRef]
  6. R. Rao, “Statistics of the fractal structure and phase singularity of a plane light wave propagation in atmospheric turbulence,” Appl. Opt. 47, 269-276 (2008). [CrossRef] [PubMed]
  7. D. P. Kelly and B. M. Hennelly, “Analytical and numerical analysis of linear optical systems,” Opt. Eng. 45, 088201 (2006). [CrossRef]
  8. A. Stern and B. Javidi, “Analysis of practical sampling and reconstruction from Fresnel fields,” Opt. Eng. 43, 239-250(2004). [CrossRef]
  9. A. Brodzik, “On the Fourier transform of finite chirps,” IEEE Signal Process. Lett. 13, 541-544 (2006). [CrossRef]
  10. S. M. Flatte and J. S. Gerber, “Irradiance-variance behavior by numerical simulation for plane-wave and spherical-wave optical propagation through strong turbulence,” J. Opt. Soc. Am. A 17, 1092-1097 (2000). [CrossRef]
  11. A. Belmonte, “Feasibilty study for the simulation of beam propagation: consideration of coherent lidar performance,” Appl. Opt. 39, 5426-5445 (2000). [CrossRef]
  12. B. M. Hennelly and J. T. Sheridan, “Generalizing, optimizing, and inventing numerical algorithms for the fractional Fourier, Fresnel, and linear canonical transforms,” J. Opt. Soc. Am. A 22, 917-927 (2005). [CrossRef]
  13. Computation time for the FFT algorithm is minimized when N is a power of 2, although depending on the value, other lengths can be computed nearly as fast. The value of N=250 was chosen here to simplify our illustration.
  14. The size is chosen so an odd number of samples spans the aperture and one sample can be placed on the optical axis. In this case 51 samples×0.2 cm=10.2cm width.
  15. The analytic Fourier transform of the chirp q(x,y)=exp(j(k/2z)(x2+y2)) is Q(fX,fY)=jλz exp(−jπλ(fX2+fY2)). In starting with a sampled version of Q, the inverse FFT of Q typically needs to be multiplied by N2ΔfX2 to reproduce the correct magnitude for the sampled version of q (assumes 1/N2 scaling on inverse 2D FFT).
  16. C. Rydberg and J. Bengtsson, “Efficient numerical representation of the optical field for the propagation of partially coherent radiaction with a specified spatial and temporal coherence function,” J. Opt. Soc. Am. A 23, 1616-1625 (2006). [CrossRef]
  17. X. Deng, B. Bihari, J. Gan, F. Zhao, and R. T. Chen, “Fast algorithm for chirp transforms with zooming-in ability and its applications,” J. Opt. Soc. Am. A 17, 762-771 (2000). [CrossRef]
  18. L. Onural, “Exact analysis of the effects of sampling of the scalar diffraction field,” J. Opt. Soc. Am. A 24, 359-367 (2007). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited