Optical coherence theory typically deals with the average properties of randomly fluctuating fields. However, in some circumstances the averaging process can mask important physical aspects of the field propagation. We derive a new method of simulating partially coherent fields of nearly arbitrary spatial and temporal coherence. These simulations produce the expected coherence properties when averaged over sufficently long time intervals. Examples of numerous fields are given, and an analytic formula for the intensity fluctuations of the field is given. The method is applied to the propagation of partially coherent fields through random phase screens.
©2006 Optical Society of America
Optical coherence theory, which has its origins in studies of the coherence of sunlight , is now a well-established discipline in optical science and numerous books [2, 3] deal with the general theory and applications. Optical coherence theory typically deals with the average properties of randomly fluctuating fields. However, the averaging process can in some significant instances mask important physical aspects of the field behavior. For instance, there has been much evidence that partially coherent fields are less susceptible to turbulence degradation than their fully coherent counterparts (see, for instance, [4, 5, 6, 7, 8], and earlier references therein). This has lead to the possibility of using partially coherent fields as sources in free-space optical communications. Such problems, however, involve (at least) three significant time scales: the coherence time of the field, the rate at which turbulence changes with time (the Greenwood frequency), and the data communication rate. Coherence theory calculations involve a long time average over all three scales, masking any possible issues that might arise over finite intervals of time. An example of one such issue is beam wander in turbulence .
Other important field properties can be ‘hidden’ by the averaging process. For instance, in recent years there has been much interest in beams possessing orbital angular momentum, or optical vortices . Such vortices are characterized by an intensity null at their center and a helical phase front, and they are stable under amplitude and phase perturbations of the field. When the vortex field is partially coherent, or its coherence is reduced on propagation through turbulence, the vortex position fluctuates and no point in space, on average, possesses an intensity null. It can be said that the vortex is hidden; it has been shown that some of its behavior is preserved in the correlation properties of the field [11, 12]. At any instant of time, however, the vortex is present in the field.
In this paper we derive a method of numerically generating realizations of partially coherent fields of nearly arbitrary spatial and temporal coherence. The spectral properties, spatial coherence properties, and intensity profile of the field can be freely and independently chosen. The method is an extension and reimagining of a technique  used to study the invariance properties of random fields in dispersive media; it may also be considered an extension of techniques for studying random electrical noise .
The paper is organized as follows. In section 2 the partially coherent field generator is introduced. In section 3 the intensity fluctuations of the generated field are determined. In section 4 a number of examples of partially coherent fields are described, and their propagation through random phase screens is considered. Section 5 presents concluding remarks.
2. Generating realizations of fields of arbitrary spatial and temporal coherence
The partially coherent field generator is introduced as a generalization of the technique used in Ref.  to study the propagation of partially coherent fields in dispersive media; in that reference, only temporal coherence properties were considered. As the goal of this paper is to generate a realization of a field with specified spatial and temporal coherence properties, we first introduce a model of a partially coherent field which can be implemented numerically; we then determine how the parameters of this model are related to the average properties of the partially coherent field.
We consider a source of partially coherent radiation in the plane z = 0, which emits optical pulses of fixed spatial and temporal shape at random times. We initially restrict ourselves to a time interval t ∈ [-T/2,T/2]. Assuming that the pulses are emitted independently of one another, the probability p(N) that N are emitted in this interval is dictated by Poissonian statistics, i.e.
Let us assume that N pulses are emitted in this interval. The field of these N pulses is then given by
where Λ(r,t) is the field amplitude of a single pulse in the plane z = 0 at transverse position r = (x,y), aside from a linear phase term, tj is the time of emission of the jth pulse and K j is the angle of inclination of the jth pulse. The time of arrival is assumed to be a random variable uniformly distributed throughout the interval, and the angle of inclination is a random variable whose probability distribution P(K) is for now unspecified. This representation of the field is very similar to that used in Ref. ; however, the introduction of the angle of inclination K j and its probability distribution allows us to control the spatial coherence properties of the field as well as the temporal properties.
Once the pulse shape Λ(r,t), the average number of pulses N̅ and the probability distribution P(K) are specified, realizations of the field can be generated. We now consider how these quantities are related to the average properties of the field, such as the mutual coherence function.
The mutual coherence function of a statistically stationary field V(r,t) is defined as
where τ = t 2 -t 1 and the angle brackets denote ensemble averaging. It is to be noted that this ensemble average is equivalent to three independent averages: the average over the arrival times tj of the pulses, the average over the inclination factors K j , and the average over the number of pulses per interval N. The instantaneous form of this function for our field of N pulses is
We first wish to evaluate the functional form of this quantity when we average over the arrival times ti . To do so, we express Λ(r,t) in terms of its temporal Fourier transform, i.e.
Our expression (4) may then be written in the form
The time average of a function of ti and tj over the interval can be written as
We may therefore write
where sinc[x] = sin[x]/x. We may separate our expression (6) into two distinct sums: one for which i = j and one for which i ≠ j. We may then write
At this point, the angle brackets on the exponentials refer only to averaging over the inclination factors K j . We now use our probability density function to evaluate this average. We have
where P̃(r) represents the two-dimensional Fourier transform of the probability density function, defined by
The sums may then be evaluated, and we find that
At this point the integrals are completely independent of the individual realizations of pulses, i.e. the particular values of tj and K j for each pulse, and the only random variable remaining is the number of pulses in the interval. We may average over this quantity as well, to get the mutual coherence function as
where p(N) is the Poisson distribution. We need the well-established results
The mutual coherence function may then be written as
Defining the average rate of pulse emission as η= N̅/T, we may rewrite Eq. (18) in the form
Letting the measurement interval T → ∞, the functions f(ω)T and g(ω)T reduce to
where δ(ω) is the Dirac delta function and δ (e) is the even half-delta function, defined such that
By use of these results, the mutual coherence function becomes
The latter term is the DC-contribution to the field. If we are considering sufficiently narrowband optical signals, it may be safely neglected.
One more simplification will be convenient. We assume that the spatial and temporal parts of the field factorize, i.e. that
We may then write
If we consider this field in the space-frequency domain, and consider instead the cross-spectral density W(r 1,r 2,ω) of the field, we have
Equation (27) is the main result of this paper. It demonstrates that the cross-spectral density of our realization of pulses will have a spectrum | (ω)|2, an average field profile Θ*(r 1) and a spectral degree of coherence P̃(r 2 - r 1). These three functions can be chosen independently of one another, and we may therefore construct a realization of a partially coherent field which has quite general spatial and temporal coherence.
It is to be noted that our results are not completely general, i.e. there exist partially coherent fields which cannot be realized by this random pulse technique. In particular, there exist fields whose spatial and temporal coherence properties are not factorizable, and there also exist fields whose spectral degree of coherence is frequency dependent. Nevertheless, our method provides an excellent tool for creating realizations of quite arbitrary spatial and temporal coherence.
3. Intensity fluctuations of partially coherent field realizations
We have as yet only considered the second-order coherence properties of the optical field. It is worthwhile, however, to investigate the fourth-order coherence properties achievable by our simulation method as well, as intensity fluctuations of an optical signal (in particular, its scintillation index ) play an important role in optical communications.
The scintillation index is defined as [15, Sec. 1.7]
where I(r,t) = 〈|V(r,t)|2〉 is the intensity of the field which is on average independent of time. The quantity 〈I(r, t)〉2 can be derived from the results of the previous section, so we focus on the quantity 〈I(r,t)2〉. Using our collection of pulses over a finite interval T again, we have
where VN (r,t) is defined by Eq. (2). Looking at the temporal Fourier decomposition of this equation, we find that
There are N 4 terms to the summation, but most of them result in either a zero-frequency contribution which will be neglected or a negative-frequency term which is identically zero. The only non-zero terms are those for which
These averages can be calculated as in the previous section. We have
Similarly, for the spatial average, we have
Our expression for the intensity fluctuations takes on the form
We again factorize our field into a temporal and spatial part. Noting that the frequency dependence of the first term is a convolution, and the second part is simply 〈I(r,t)〉2, we may write
On substitution into the scintillation index, we readily find that
In the limit η → ∞, i.e. a large rate of pulse emission, we find that σI ~ 1. This result is consistent with a light field said to be chaotic or Gaussian [16, chapter 3]. Our simulation method therefore cannot produce a pure coherent laser field. However, partially coherent fields derived from laser light, for instance by passing coherent light through a rotating ground glass plate, are well-known to be chaotic  and our simulation method will reliably simulate many partially coherent field configurations.
Although = 1 is the ideal limit for our simulation method, for small values of η (small rate of pulse emission) the value of can be much higher. We can estimate how large η must be to reach the ideal limit by considering Gaussian pulses,
On substitution into our formula (38), we readily find that
We can approach the ideal value of ~ 1 by taking ησt to be sufficiently large.
We now demonstrate the use of the field generator by numerical examples.
4.1. Gaussian Schell-model fields
We consider fields of Gaussian intensity profile and Gaussian spatial correlation, known as Gaussian Schell-model fields [3, Sec. 5.2.2]. The mutual coherence function of such fields may be written as
where I(r) is the average field intensity,
σI being the beam width, μ(r 2 -r 1) is the spatial correlation function (equivalent to the spectral degree of coherence in the frequency domain),
σg being the correlation length, and γ(τ) is the temporal coherence function, to be taken as Lorentzian or Gaussian. On comparison with Eq. (27), it can be seen that our field generator should generate a Gaussian Schell-model field if we take and take | (ω)|2 to be the temporal Fourier transform of γ(τ). Figure 1 illustrates the intensity of the field generated by our simulation method for several realizations, with σI = 2 cm, σg = 1 cm, and Gaussian spectrum of center frequency 1 × 1015 Hz and 1% bandwidth. The average pulse rate is taken to be 5 pulses/cycle. The pictures show the gradual evolution of the field in time; the frames are each separated by 5 periods at the center frequency.
To be a valid technique for generating realizations, the field must possess the proper prescribed average properties. Figure 2 illustrates the average intensity of the field, taken over 50 instantaneous values of the field each separated by 20 periods at the center frequency. Part (a) shows the cross-section of the beam, while (b) shows the cross-section of the beam along the line y = 0. The ideal Gaussian is shown as a dashed line, and it can be seen that there is excellent agreement. Convergence could be further improved by extending the duration of the time average.
The spatial correlation properties can also be numerically calculated to test the technique. Figure 3 shows the average spatial correlation properties, taken over 50 instantaneous values of the field each separated by 20 periods at the center frequency. The spatial correlation properties were calculated at points +x, -x along the line y = 0. The circles represent the results generated from our realization, while the dashed line represents the ideal Gaussian Schell-model case. Again it can be seen that there is excellent agreement.
The average temporal correlation properties of the field can also be numerically calculated. Figure 4 shows the complex degree of coherence γ(τ) calculated at the center of a coherent Gaussian beam, for a Gaussian and Lorentian lineshape. Again there is excellent agreement between the results of the simulation and the expected average behavior.
4.2. Propagation through random phase screens
We now consider the propagation of our partially coherent field realizations through random phase screens. Such phase screens are extensively used to model the effects of beam propagation through weak turbulence . The field must be propagated from the source plane to the plane of the phase screen, and then propagated beyond it to a detector plane. Though such Fresnel-type propagation has of course been done many times, is it uncommon to see it in the time domain, so we briefly summarize the technique. We may use Huygens’ principle to express the propagation of the field V(r′,t) from the plane z = 0 in the form
where V 0 (r′, t) is the field in the source plane z = 0,R= |r - r′| is the source point-field point distance, and the integration is carried out over the source plane. For quasi-monochromatic fields such as we are considering, we may separate the source field into the components
where ω 0 is the central frequency of oscillation of the source and F 0(r,t) is the slowly-varying piece of the source field. We may then write
We are also primarily interested in paraxial fields, for which
The field may then be simplified to the form
Since the function F 0 is assumed to be slowly-varying with respect to time and paraxial, the spatial argument which depends on the tranverse coordinates will be negligible compared to the other terms. We may simplify our field calculation to the form
where t̂ = t - z/c. The field in any plane of constant z is therefore simply the time-shifted Fresnel transform of the field in the source plane, and can be calculated straightforwardly with a fast Fourier transform.
We study the response of the partially coherent field on propagation through a Gaussian random phase screen with power spectrum
where L 0 is the outer scale of the phase screen. Such a screen is a poor model for quantitative studies of atmospheric turbulence (in such cases more sophisticated spectral models exist ), but will suffice for a brief illustration of partially coherent field propagation through random media. The screens are generated by the method described in Ref. . We take as screen parameters = 1, L 0 = 15cm and the screen is located at z = 1.0km. A spatially coherent field with σI = 2cm is propagated through this screen and to a detector at z = 2km; the intensity at the detector is shown in Fig. 5(a). It can be seen that the spot has been aberrated by the phase screen and, perhaps more important, is no longer centered on the center of the detector plane. This is an example of beam wander; it is to be noted that if the beam wanders sufficiently from the axis, it may not illuminate the detector at all, resulting in information loss (a ‘miss’). Temporal fluctuations of the coherent field do not fix the problem; as illustrated in Fig. 5(b), over a 200 cycle time average of the field, the spot shape and position remain essentially unchanged.
We now consider the propagation of a partially coherent field with σg = 1cm through the same phase screen. As illustrated in Fig. 6(a), this field is also significantly distorted by the phase screen. However, some field intensity is still present in the center of the detector plane. When a long time average is taken, as shown in Fig. 6(b), we see that there is still a tendency for the field to ‘wander’ from the center of the screen, but an appreciable amount of field intensity remains near the center of the detector plane.
These results suggest the following picture for the improved propagation characteristics of partially coherent fields in turbulence, as illustrated in Fig. 7. A coherent laser essentially propagates its energy through a single coherent mode, which is subject to distortion on propagation through the inhomogeneous medium. The signal may or may not arrive fully at the detector, and will have suffered degradation due to diffraction. The partially coherent beam sends its energy through multiple (independent) modes, each of which propagates differently in the turbulent medium. Although any individual mode may not arrive cleanly at the detector, it is likely that at least one will make the journey. For our particular method of generating realizations, the modes are the individual pulses, each of which has a different ‘angle of attack’ (value of K j ) into the phase screen.
We have derived a new method of generating realizations of partially coherent fields of nearly arbitrary spatial and temporal coherence. The method was demonstrated for Gaussian Schell-model fields of different degrees of spatial and temporal coherence, and the numerically calculated average properties agree well with the predicted analytical results. The method was also used to gain some insight into the propagation of partially coherent beams through random phase screens, which could be used as a rough model for atmospheric turbulence. It is expected that the ability to study the behavior of fluctuating fields on many time scales (center frequency, bandwidth, rate of turbulence fluctuations, detector response) will provide useful insight into the behaviors of many coherence-related phenomena.
It is to be noted that this method can also be used as an educational tool. Coherence theory is often difficult for the beginner to understand. This method of generating realizations can be used as an additional visual/conceptual aid for the student of optical coherence, as well as a technique for generating unique homework assignments - each student could be given a realization of a field and asked to determine its spatial and temporal properties.
The research was supported by the US Air Force Office of Scientific Research under grant FA 9550-05-1-0288.
References and links
01. È Verdet, Leçons d’Optique Physique, vol. 1, (L’Imprimierie Impériale, Paris, 1869).
02. M. Born and E. Wolf, Principles of Optics (Cambridge University Press, Cambridge, 1999, 7th (expanded) edition).
03. L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, Cambridge, 1995).
04. J. Wu, “Propagation of a Gaussian-Schell beam through turbulent media,” J. Mod. Opt. 37, 671–684 (1990). [CrossRef]
05. G. Gbur and E. Wolf, “Spreading of partially coherent beams in random media,” J. Opt. Soc. Am. A 19, 1592–1598 (2002). [CrossRef]
06. J. Wu and A. D. Boardman, “Coherence length of a Gaussian Schell-model beam and atmospheric turbulence,” J. Mod. Opt. 38, 1355–1363 (1991). [CrossRef]
07. J. C. Ricklin and F. M. Davidson, “Atmospheric turbulence effects on a partially coherent Gaussian Beam: implications for free-space laser communication,” J. Opt. Soc. Am. A 19, 1794–1802 (2002). [CrossRef]
08. O. Korotkova, L. C. Andrews, and R. L. Phillips, “A Model for a Partially Coherent Gaussian Beam in Atmospheric Turbulence with Application in Lasercom,” Opt. Eng. 43, 330–341 (2004). [CrossRef]
09. L.C. Andrews and R.L. Phillips, Laser Beam Propagation through Random Media (SPIE Press, Bellingham, Washington, 1998).
10. M.S. Soskin and M.V. Vasnetsov, Singular Optics, in Progress in Optics, ed. E. Wolf (Elsevier, Amsterdam), vol. 42 (2001), p. 219–276.
11. G. Gbur, T.D. Visser, and E. Wolf, “‘Hidden’ singularities in partially coherent wavefields”, J. Opt. A 6, S239–S242 (2004). [CrossRef]
12. G. Gbur and T.D. Visser, “Phase singularities and coherence vortices in linear optical systems,” Opt. Commun. 259, 428–435 (2006). [CrossRef]
13. W. Wang and E. Wolf, “Invariance properties of random pulses and of other random fields in dispersive media,” Phys. Rev. E 52, 5532–5539 (1995). [CrossRef]
14. S.O. Rice, Mathematical analysis of random noise, in Selected Papers on Noise and Stochastic Processes, ed. N. Wax (Dover, NY, 1954).
15. L.C. Andrew and R.L. Phillips, Laser Beam Scintillation with Applications (SPIE Press, Bellingham, Washington, 2001). [CrossRef]
16. R. Loudon, The Quantum Theory of Light (Oxford University Press, Oxford, 1983, 2nd edition).
17. F.T. Arecchi, E. Gatti, and A. Sona, “Time distribution of photons from coherent and Gaussian sources,” Phys. Lett. 20, 27–29 (1966). [CrossRef]
18. D.L. Knepp, “Multiple phase-screen calculation of the temporal behavior of stochastic waves,” Proc. IEEE 71, 722–737 (1983). [CrossRef]