|
|
1.Introduction.Let me begin by saying that I have known Jim for 33 years. This is what he looked like when I joined the Optical Sciences Center (now College of Optical Sciences) in the Fall of 1988. You can tell that he hasn’t changed a bit. In contrast, see how the rest of us have aged; this is what I looked like in those days! For reasons that will become clear in a few minutes, whenever I think of Jim Wyant, the following story pops into my head. A destitute man happens to meet the Lord Almighty. He asks: “God, what’s a million years to you?” The Lord responds: “Oh, nothing; it’s just like a second to me.” The poor man asks: “God, what’s a million dollars to you?” The Lord says: “Nothing; it’s like a penny to me.” The man finally says: “Dear Lord, may I have a penny?” To which the Lord responds: “Absolutely; please wait a second!” ![]() Now, here’s the background story. There was a time that I thought I had many rich and generous friends. So, I had this idea of asking each of them to donate half a million dollars to our College, so that we could create endowed scholarships in their own names. I asked Jim to accompany me to the office of the University’s vice-president to see if she would agree to match such donations by covering the students’ tuition and other expenses. To our delight, vice-president Leslie Tolbert loved the idea and consented to our request. Moreover, she allowed us to have as many as 30 such endowed scholarships with matching provided from her office’s budget. On our way back from the VP’s office, I turned to Jim and said: “Jim, why don’t you become a donor and give us the first half-million dollars to start this Friends of Tucson Optics (FoTO) scholarship program. To which Jim responded: “You’ll never get a penny from me!” Fortunately for our College and for our students (then and forever after), we didn’t have to wait a million years for Jim to change his mind! In less than two weeks, Jim had decided to endow the first FoTO scholarship—to be named after his beloved deceased wife, Louise. We went on to receive a few more gifts from generous individuals and established a couple more endowed FoTO scholarships. A few years later, Jim announced that he would donate another 10 million dollars to the College, provided that the College would use the money for a four-to-one matching of $100,000 gifts. That is when the money started pouring in and, in a matter of two to three years, we had all thirty FoTO scholarships fully endowed, thanks to Jim’s generosity. As for my rich and generous friends that I had initially intended to rely upon, it turned out that my rich friends weren’t generous, and my generous friends weren’t rich! With regard to my presentation today, I wanted to do something that touched upon the educational aspects of our College. As you all know, Jim is a gifted teacher and, to my eternal consternation (!), our students have always had a much higher opinion of him as an instructor than they’ve had of the rest of us on the faculty. So, I thought I should contribute a paper to this symposium that might have some educational value. That is when I hit upon the unifying theme of Fourier transformation. Our colleague, Prof. Sasian, edited a special edition of Applied Optics last year. This was a commemorative issue in honor of the renaming of our College to James C. Wyant College of Optical Sciences. My plan was to submit a paper entitled “Ubiquity of Fourier Transformation in Optical Sciences” to that special issue of Applied Optics, then present some of its results at this Tribute Symposium to Dean Wyant, which was initially scheduled for last August. Unfortunately, the pandemic intervened and the Tribute was postponed until this week. That Applied Optics paper1 has now been in the public domain for more than a year and, in any event, there isn’t much time available today for me to discuss the technical aspects of the paper in any detail. So, I’ll just mention the kind of topics that are covered in the paper, hoping that you will find them sufficiently interesting to look up the details on your own. The published version of this presentation, however, will contain a few more applications of Fourier transformation that were not covered in the Applied Optics paper. The organization of this paper is as follows. In Sec.2, I describe the central limit theorem of probability and statistics, which is the fundamental mathematical argument as to why many natural phenomena appear to have a Gaussian probability distribution. The comb function, comprising an infinite array of identical, equi-spaced delta-functions, is the subject of Sec.3, where I use two different methods to show that comb(s) is the Fourier transform of comb(x). Section 4 describes how one may use the comb function to derive the Fourier series coefficients of periodic functions. Another application of the comb function in the context of the Shannon-Nyquist sampling theorem is the subject of Sec.5. Finally, in Sec.6, I use a Fourier-transform-based method to solve Maxwell’s equations of classical electrodynamics and derive the radiation field of a magnetic dipole oscillator. 2.The central limit theorem.It is a well-known fact that the Fourier transform of the unit rectangular function, rect(x), is sinc(s) = sin(πs)/(πs).2 Successive convolutions of rect(x) with itself yield tri(x) = rect(x) * rect(x), whose Fourier transform is sinc2(s), and hump(x) = rect(x) * rect(x) * rect(x), whose Fourier transform is sinc3(s).† Clearly, the N-fold convolution of rect(x) with itself yields a function of total width N, unit area, and Fourier transform sincN(s), whose profile along the x-axis approaches that of a Gaussian function as N continues to rise. Considering that sin(πs) = πs – (πs)3/3! + ···, it is seen that The method described in the preceding paragraph is quite general and can be applied to the convolution of a large number of functions that may or may not be identical, say, f1(x) * f2(x) * ··· * fN(x), provided that the individual functions satisfy certain constraints. As an example, let us consider a large number of independent random variables, say, x1, x2, · · ·, xN, having respective probability distributions p1(x), p2(x), ···, pw(x). In general, pn(x) ≤ 0 and Next, we define an average value Denoting the so-called characteristic function2-4 of the probability density It is now easy to verify that, for the N independent random variables under consideration, The inverse Fourier transform of Eq.(8) now yields the probability density function of xsum as This is the main conclusion of the central limit theorem of probability and statistics, explaining why so many natural phenomena, each arising as the superposition of a large number of more or less independent random variables, tend to have a Gaussian probability distribution. 3.The comb function.The function comb(x) is a periodic function of the real variable x, where each period consists of a single delta-function; that is, Now, in the limit when N → ∞, the function sin[π(2N + 1)s] oscillates rapidly with a changing s, and one must carefully examine its behavior in the vicinity of sn = nπ (i.e., around the points where s is an integer multiple of π), because the denominator sin(πs) as well as the numerator are very close to zero around such points. It is not difficult to see that, in the immediate vicinity of each sn, the function approaches (2N + 1)sinc[(2N + 1)(s – sn)], where, as before, the sinc function is defined as sinc(s) = sin(πs)/(πs). Thus, in the limit when N → ∞, the function appearing on the right-hand side of Eq.(1) approaches comb(s). An alternative way of demonstrating the same thing involves raising the function cos(πx) to the power 2N, where N is a positive integer, then allowing N to go to infinity. In the immediate vicinity of xn = n, where n is a positive, zero, or negative integer, the function cos2(πx) is very close to 1, which makes its Nth power also close to 1. However, as soon as x moves away from xn, cos2N(πx) drops below zero and its Nth power rapidly approaches zero as N → ∞. Thus, in cos2N (πx), we have a periodic function consisting if narrow bumps of height 1, centered at xn = n, which approaches zero everywhere else when N → ∞. The area under each bump is readily found using the method of integration by parts, as follows: Considering that Next, we consider the Fourier transform of This is an array of δ-functions centered at s = 0 and extending (symmetrically) over integer values of s from –N to N. The magnitude of the largest δ-function, located at s = 0, is The adjacent δ-functions located at s = ±1 are slightly smaller, having magnitude N/(N + 1). The next adjacent δ-functions located at s = ±2 are a bit smaller, at N(N – 1)/(N + 1) (N + 2), and so on. It should thus be clear that, as N → ∞ and the function (2M)!! cos2N(πx)/(2N – 1)!! approaches comb(x), its Fourier transform, a tapered array of 2N + 1 delta-functions centered at s = 0, approaches comb(s). 4.The Fourier series.Many textbooks on Fourier transform theory begin by introducing the Fourier series, associated with periodic functions, as a first step toward discussing the Fourier integral for an arbitrary function f(x). However, the reverse approach of starting with the Fourier integral Thus, the Fourier transform of a periodic function consists of an array of equi-spaced δ-functions (located at sn = n/p), where the amplitude of the nth delta-function is given by The inverse Fourier transform now yields the Fourier series of the periodic function, as follows: In the special case when fp(x) is real-valued, Eq.(15) yields In the above equation, In this way, the Fourier series coefficients of a periodic function f(x) of period p, are obtained as the normalized sampled values p–1Fp(n/p) of the Fourier transform Fp(s) of a single period fp (x), in accordance with Eq.(16). Alternatively, if f(x) happens to be real-valued, the Fourier coefficients are derived from the sine and cosine transforms of fp(x), in accordance with Eqs.(17). 5.The Shannon-Nyquist Sampling theorem.In modern digital signal processing, electronic communication, and information storage, analog waveforms are sampled at regular intervals, then converted to a digital format.2,4 The Fourier transform F(s) of the waveform (or signal) under consideration, namely, f(t), is said to be bandlimited if F(s) turns out to be negligible outside the frequency interval (–w, w). Under such circumstances, multiplying f(t) with the sampling function p–1comb(t/p), namely, an array of unit δ-functions at regular intervals Δt = p, creates a sampled function p–1comb(t/p)/(t), whose Fourier transform, comb(ps) * F(s), may be written as During the reconstruction of the original signal f(t) from its uniformly-spaced samples f(np), the discretized function p–1comb(t/p)/(t) is sent through a linear, shift-invariant, low-pass filter, whose transfer function H(s) = ρ rect(ps) multiplies the Fourier transform of the input waveform. Assuming the sampling rate satisfies the aforesaid Nyquist criterion, the individual copies of F(s) that are present within the spectrum of the input waveform will be well-separated from each other, thus allowing the filter to extract the original (unshifted) copy of F(s) and perfectly reproduce the signal f(t) at its output. The inverse Fourier transform of H(s) is the so-called impulse-response of the filter, which is presently given by h(t) = sinc(t/p). The output f(t) of the filter, being the convolution between the input waveform and the impulse-response, is thus seen to be The function sinc(t/p), whose properly shifted and scaled copies reconstruct the original signal f(t) in accordance with Eq.(18), is known as the interpolation function.2 6.Radiation by an oscillating magnetic point-dipole.Consider a magnetic point-dipole Fourier transforming the above magnetization profile in four-dimensional spacetime (r, t), we find The electric current-density associated with our point-dipole is As for the scalar potential ψ(r, t) produced by the dipole (again, in the Lorenz gauge), we note that, since magnetic dipoles do not possess an electric charge, that is, ρ(r, t) = 0, their scalar potential Our next task is to determine the spatio-temporal distribution A(r, t) of the vector potential by inverse Fourier transforming This is neither the retarded nor the advanced vector potential at this point. To arrive at the retarded potential, we need to include an appropriate superposition of plane-waves that exist, independently of our oscillating dipole, in free space. Such plane-waves have a temporal frequency ω = ω0 and reside on a spherical surface of radius k = ω0/ω in the k-space.1 It is not difficult to guess that the desired distribution in the (k, ω) space of this vacuum vector potential should be Here, δ(ω ± ω0) ensure that the vacuum potential oscillates at the frequency ω0; similarly, the confinement of the k-vectors to the requisite spherical surface is guaranteed by δ[k – (ω0/c)]. The term Adding the vacuum potential of Eq.(24) to that of the magnetic dipole given by Eq.(22), we finally arrive at the retarded vector potential of the oscillating point-dipole, namely, This is the exact solution of Maxwell’s equations for the vector potential produced by our magnetic point-dipole. Given that the corresponding scalar potential ψ(r, t) equals zero, it is easy to verify the satisfaction of the Lorenz gauge condition ∇ · A + c–2 ∂ψ/∂t = 0.5 The electromagnetic fields may now be derived using The rate of flow of electromagnetic energy is given by the time-averaged Poynting vector, namely, Here, As expected,5 the radiated power P is seen to be proportional to the square of the dipole moment m0 and the fourth power of its oscillation frequency ω0. A personal note of gratitude. I am grateful to Dean James Wyant for, among many other things, providing an environment in which my colleagues and I could thrive, engaging in teaching and research without feeling pressured, and without being burdened by excessive bureaucratic concerns. It has been an honor to be able to call Jim my colleague (and my friend) for the past 33 years. ReferencesM. Mansuripur,
“Ubiquity of Fourier transformation in Optical Sciences,”
Applied Optics, 59 G47
–G63
(2020). https://doi.org/10.1364/AO.390342 Google Scholar
R. N. Bracewell, The Fourier Transform and Its Applications, 3McGraw-Hill, New York
(1999). Google Scholar
A. Papoulis, Probability, Random Variables, and Stochastic Processes, 2McGraw-Hill, New York
(1984). Google Scholar
M. Mansuripur, Mathematical Methods in Science and Engineering: Applications in Optics and Photonics, Cognella Academic Publishing, San Diego
(2020). Google Scholar
J. D. Jackson, Classical Electrodynamics, 3Wiley, New York
(1999). Google Scholar
Notes[1] Here, we are using the standard definition of rect(x) = 1 when |x| ≤ ½, and 0 otherwise. Similarly, tri(x) = 1 – |x| when |x| ≤ 1, and 0 otherwise.2 The name used for the function hump(x) = ¾ – x2 when |x| ≤ ½, and [2] The convention used in this paper with regard to the magnetic fields is B = μ0H + M; consequently, the magnetization M has the units of the B-field. This means that our magnetic dipole |