OSA's Digital Library

Journal of the Optical Society of America A

Journal of the Optical Society of America A


  • Vol. 15, Iss. 5 — May. 1, 1998
  • pp: 1320–1328

Reducing canonical diffraction problems to singularity-free one-dimensional integrals

G. W. Forbes and A. A. Asatryan  »View Author Affiliations

JOSA A, Vol. 15, Issue 5, pp. 1320-1328 (1998)

View Full Text Article

Enhanced HTML    Acrobat PDF (496 KB)

Browse Journals / Lookup Meetings

Browse by Journal and Year


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools



The oscillatory integrands of the Kirchhoff and the Rayleigh–Sommerfeld diffraction solutions mean that these two-dimensional integrals typically lead to challenging computations. By adoption of the Kirchhoff boundary conditions, the domain of the integrals is reduced to cover only the aperture. For perfect spherical (both diverging and focused) and plane incident fields, closed forms are derived for vector potentials that allow each of these solutions to be further simplified to just a one-dimensional, singularity-free integral around the aperture rim. The results offer easy numerical access to exact—although, given the approximate boundary conditions, not rigorous—solutions to important diffraction problems. They are derived by generalization of a standard theorem to extend previous results to the case of focused fields and the Rayleigh–Sommerfeld solutions.

© 1998 Optical Society of America

OCIS Codes
(000.3860) General : Mathematical methods in physics
(050.1940) Diffraction and gratings : Diffraction
(050.1960) Diffraction and gratings : Diffraction theory
(260.1960) Physical optics : Diffraction theory
(350.7420) Other areas of optics : Waves

Original Manuscript: July 25, 1997
Revised Manuscript: December 16, 1997
Manuscript Accepted: December 19, 1997
Published: May 1, 1998

G. W. Forbes and A. A. Asatryan, "Reducing canonical diffraction problems to singularity-free one-dimensional integrals," J. Opt. Soc. Am. A 15, 1320-1328 (1998)

Sort:  Author  |  Year  |  Journal  |  Reset  


  1. See the review given in A. Rubinowicz, “The Miyamoto–Wolf diffraction wave,” Prog. Opt. 4, 201–240 (1965), and the more historical account given in A. Rubinowicz, “Thomas Young and the theory of diffraction,” Nature 180, 160–162 (1957). [CrossRef]
  2. T. Gravelsaeter, J. J. Stamnes, “Diffraction by circular apertures. 1. Method of linear phase and amplitude approximation,” Appl. Opt. 21, 3644–3651 (1982). [CrossRef] [PubMed]
  3. W. B. Gordon, “Vector potentials and physical optics,” J. Math. Phys. 16, 448–454 (1975). [CrossRef]
  4. J. S. Asvestas, “Line integrals and physical optics. Part I. The transformation of the solid-angle surface integral to a line integral,” J. Opt. Soc. Am. A 2, 891–895 (1985), and Part II. The conversion of the Kirchhoff surface integral to a line integral,” 2, 896–902 (1985). [CrossRef]
  5. See W. Kaplan, Advanced Calculus (Addison-Wesley, Reading, Mass., 1973), p. 349. A variant of Green’s theorem leads to the more standard, but far less convenient, volume integral for this potential; see, e.g., W. K. H. Panofsky, M. Phillips, Classical Electricity and Magnetism (Addison-Wesley, London, 1964), p. 2, and G. Arfken, Mathematical Methods for Physicists (Academic, Orlando, Fla., 1985), Sec. 1.15.
  6. E. W. Marchand, E. Wolf, “Boundary diffraction wave in the domain of the Rayleigh–Kirchhoff diffraction theory,” J. Opt. Soc. Am. 52, 761–767 (1962). [CrossRef]
  7. For example, the result deserves to be in books such as in J. D. Jackson, Classical Electrodynamics (Wiley, New York, 1975), Chap. 5; G. A. Korn, T. M. Korn, Mathematical Handbook for Scientists and Engineers (McGraw-Hill, New York, 1968), p. 165; and the second and third references given in Ref. 5.
  8. See, for example, M. Born, E. Wolf, Principles of Optics (Pergamon, Elmsford, New York, 1980), Sec. 8.3.
  9. This result is given by A. Rubinowicz, “Die Beugungswelle in der Kirchhoffschen Theorie der Beugungserscheinungen,” Ann. Phys. 53, 257–278 (1917). See the integrals presented after Fig. 1 in this reference. [CrossRef]
  10. By use of Eqs. (3.2) and (2.7) as described, it is found that D(r)=-(1-σ) δO×δS16π2×∫01 ik exp[ik(a2-bτ)/S(τ)]S3(τ)τdτ, where S(τ)=(a2-2bτ+c2τ2)1/2,a=|rO-rS|,b=(rO-rS)·(r-rS), and c=|r-rS|. Again, this integral can be done in closed form: The primitive is just (b2-a2c2)-1 exp[ik(a2-bτ)/S(τ)]. Equation (3.10) now follows.
  11. A good first step is to note that the apparent r dependence can be reduced by use of (r-rO)×(r-rS)=(r-rO)×(rO-rS).
  12. In fact, if this line of singularities cuts a surface of integration, the associated flux contribution is readily shown to mean that, instead of vanishing, the net flux is 2π/|rS-rO|. It follows that the last of the three terms inside the braces of each of Eqs. (3.8), (3.11), and (3.12) corresponds precisely to the discontinuous geometric field component discussed in Section 1. This observation clarifies the link to the results given in K. Miyamoto, E. Wolf, “Generalization of the Maggi–Rubinowicz theory of the boundary diffraction wave. Part I,” J. Opt. Soc. Am. 52, 615–625 (1962), and “Generalization of the Maggi–Rubinowicz theory of the boundary diffraction wave. Part II,” 52, 626–637 (1962). (As indicated previously in the text, the Miyamoto–Wolf potential becomes singular when the observation point falls on the shadow boundary.) [CrossRef]
  13. See, for example, L. Mandel, E. Wolf, Optical Coherence and Quantum Optics (Cambridge U. Press, Cambridge, UK, 1995), Section 3.2.5.
  14. C. R. Schultheisz, “Numerical solution of the Huygens–Fresnel–Kirchhoff diffraction of spherical waves by a circu-lar aperture,” J. Opt. Soc. Am. A 11, 774–778 (1994). [CrossRef]
  15. In the case of the plane-wave results, recall that the incident field has already been scaled before taking the source point to infinity, so it is necessary to divide the potential given in Eq. (3.20), for example, by only exp(iku·rO), and this means that the phase factor now becomes exp(ik2L∞). The same is true of Eq. (4.4), and these results are in keeping with Eqs. (5.1) and (5.3). For the focused field, exp(-ik|rO-rS|)/(4π|rO-rS|) is to be divided out for points well before focus, but its complex conjugate is to be used far beyond the focal plane; the phase is small near the focus, so there is nothing to be gained there in this way.

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.


Fig. 1 Fig. 2

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited