Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

The van Cittert-Zernike theorem for electromagnetic fields

Open Access Open Access

Abstract

The van Cittert-Zernike theorem, well known for the scalar optical fields, is generalized for the case of vector electromagnetic fields. The deduced theorem shows that the degree of coherence of the electromagnetic field produced by the completely incoherent vector source increases on propagation whereas the degree of polarization remains unchanged. The possible application of the deduced theorem is illustrated by an example of optical simulation of partially coherent and partially polarized secondary source with the controlled statistical properties.

©2009 Optical Society of America

1. Introduction

One of the central results of optical coherence theory is the van Cittert-Zernike theorem which shows that even completely incoherent source can generate a partially coherent optical field in the process of its propagation [15]. This theorem is formulated for a scalar optical field and hence does not take into account the vector nature of the light. At the same time, in the past few years an ever-growing attention is attracted to the theory of coherence of vector electromagnetic fields or, for brevity, vector coherence theory [6]. That is why it seems to be quite natural to try to generalize the van Cittert-Zernike theorem for the vectorial case. Such a generalization is just an object of the present paper.

The paper is organized as follows. In Section 2 we shortly review the basic concepts of the vector coherence theory. When doing so, we give preference to the definition of vector coherence proposed in Refs. [7,8] as the most consistent one (see the corresponding justification in Ref. [6]). In Section 3 we deduce the vector version of the van Cittert-Zernike theorem. An example of possible application of the deduced theorem is given in Section 4. Finally, the main results are summarized in Section 5.

2. Background

We consider a stochastic stationary electromagnetic field propagating within a narrow solid angle around the z axis. According to the coherence theory in the space-frequency domain [4,5] such a field at any typical point x = (x, y) in some plane z = const and at any frequency v can be represented by the statistical ensemble

{E(x,v)}={Ex(x,v)Ey(x,v)},

where Ex and Ey are the Cartesian components of the electric field vector E, written here as a row matrix. For the sake of simplicity later on we will omit the explicit dependence of the considered quantities on frequency v. The second order correlation properties of the electromagnetic field (1) can be completely characterized by the so-called cross-spectral density matrix [5]

W(x1,x2)=[Wij(x1,x2)]=Ei*(x1)Ej(x2)(i,j=x,y),

where the asterisk and the angle brackets denote the complex conjugate and the ensemble average, respectively. For x 1 = x 2 = x, the diagonal element Wii becomes the spectral density, or power spectrum, of field component

Si(x)=Wii(x,x)=Ei(x)2,

so that the spectral density of the hole field may be expressed as

S(x)=E(x)E(x)=Si(x)+Sj(x)=TrW(x,x),

where the dagger denotes the Hermitian conjugate and Tr stands for matrix trace. It may be shown [5] that

Wij(x1,x2)Si(x1)Sj(x2).

Hence, one can define the quantitative measure of the correlation between the vector components Ei and Ej as the normalized absolute value of the cross-correlation function Wij (x 1,x 2 ) , i.e.,

μij(x1,x2)=Wij(x1,x2)Si(x1)Sj(x2),

which will be referred to as the coefficient of correlation. It is obvious that

0μij(x1,x2)1.

The vector components Ei and Ej are fully correlated when ∣μij (x 1, x 2)∣ = 1 and they are fully uncorrelated when ∣μij(x 1, x 2)∣ = 0. In classical coherence theory the coefficient of correlation given by Eq. (6), for j = i, is known as the spectral degree of coherence of a scalar wave field. This fact can serve as a key point for defining the degree of coherence of the vector electromagnetic field.

In Ref. [8] the degree of coherence for the vectorial case is defined as the normalized Frobenius (or Euclidean) norm of the cross-spectral density matrix W(x 1,x 2), i.e.,

μ˜(x1,x2)=(i,jWij(x1,x2)2i,jSi(x1)Sj(x2))1/2=(Tr[W(x1,x2)W(x1,x2)]TrW(x1,x1)TrW(x2,x2))1/2.

On making use of definition (6), one can rewrite Eq. (8) as follows:

μ˜(x1,x2)=(i,jμij(x1,x2)2Si(x1)Sj(x2)i,jSi(x1)Sj(x2))1/2.

As can be seen from Eq. (9), μ˜ (x 1,x 2) represents the power-weighted root-mean-square value of the correlation coefficients for all pairs (i, j). Equation (9), together with Eq. (7), states that μ˜ (x 1,x 2) =1 if, and only if, ∣μij (x 1, x 2)∣2 =1 for any pair (i, j), and μ˜ (x 1,x 2) = 0 if, and only if, ∣μij (x 1, x 2)∣2 = 0 for any pair (i, j). In other words, the vector electromagnetic field is completely coherent when its components are fully self- and cross-correlated and it is completely incoherent when its components are fully self- and cross-uncorrelated.

As well known, the polarization is another manifestation of the correlation properties of a stochastic electromagnetic field that can be characterized by the cross-spectral density matrix taken for x 1 = x 2 = x. When

W(x,x)=Wxx(x,x)[1001],

the field is said to be completely unpolarized, and when

DetW(x,x)=0,

where Det stands for determinant, the field is said to be completely polarized. The polarization can be evaluated quantitatively by the degree of polarization defined as [3–5]

P(x)=(14DetW(x,x)[TrW(x,x)]2)1/2.

Of course, for the completely polarized field P(x) =1 and for completely unpolarized field P(x) = 0. It may be shown that the polarization properties of a stochastic electromagnetic field can be characterized by the degree of coherence μ˜ (x 1,x 2) as well [7]. On making use of the relation

2DetW(x,x)=[TrW(x,x)]2Tr[W(x,x)W(x,x)],

one finds

P(x)=2μ˜2(x,x)1.

Equation (14) reveals an important fact: the completely coherent electromagnetic field is necessarily completely polarized. Indeed, for completely coherent field the equality μ˜ (x 1,x 2) = 1 holds for each pair (x 1,x 2), including the case x 1 = x 2 = x, and hence P(x) = 1. Equation (14) shows also that the minimum possible value of μ˜ (x, x) is 1/√2.

3. Vector generalization of the van Cittert-Zernike theorem

Now we will consider the propagation of the electromagnetic field given by the statistical ensemble (1) in free space between two parallel planes normal to the z axis. Let Ei (z)(x′) and Ei (0)(x′) be the fluctuating orthogonal field components (i = x, y) in the source plane z = 0 and any transverse plane z = constant > 0, respectively. Applying the paraxial approximation, it may be shown (see [4]) that these components are related by the following formula:

Ei(z)(x)=1iλzexp(i2πλz)(z=0)Ei(0)(x)exp[iπλz(xx)]dx,

where λ is the wave length. On substituting for Ei (z)(x′) from Eq. (15) into Eq. (2), one finds the following expression that describes the propagation of the cross-spectral density matrix W in free space:

W(z)(x1,x2)=1(λz)2(z=0)W(0)(x1,x2)exp{iπλz[(x2x2)2(x1x1)2]}dx1dx2.

On making use of Eq. (16) and definitions (8) and (12), one can determine the degree of coherence and the degree of polarization of the field in the plane z = constant. It is obvious that these quantities will, in general, change on propagation. Such changes are known as correlation-induced spectral changes [9]. Below we will examine the correlation-induced spectral changes of the electromagnetic field generated by the completely incoherent source.

Let us consider that the field in the source plane z = 0 is completely incoherent. The orthogonal components of such a field are statistically independent and self-uncorrelated. Hence, the completely incoherent source can be characterized by the cross-spectral density matrix W (0)(x 1′,x 2′) with the diagonal elements

Wii(0)(x1,x2)=Ei(0)*(x1)Ei(0)(x2)=ηSi(0)(x1)Si(0)(x2)δ(x1x2)

and off-diagonal elements

Wij(0)(x1,x2)=Ei(0)*(x1)Ej(0)(x2)=Ei(0)*(x1)Ej(0)(x2)=0,

where η is a positive constant, δ(.) is the two-dimensional Dirac delta function, and the components Ei (0) are assumed, without loss of generality, to be the fluctuations with zero mean. On substituting for W(x 1′,x 2′) from Eqs. (17) and (18) into Eq. (16) and making use of the filtering property of delta function, we find

Wii(z)(x1,x2)=η(λz)2exp[iπλz(x12x22)](z=0)Si(0)(x)exp(i2πλzx·Δx)dx

and

Wij(z)(x1,x2)=0,

where Δx = x 2 - x 1. As can be seen from Eq. (19), each field component reaching the plane z = 0 is not already self-uncorrelated. If we recall the definition (6) and use Eq. (19), we obtain the following formula for the spectral degree of coherence of the field component:

μii(z)(x1,x2)=(z=0)Si(0)(x)exp(i2πλzx·Δx)dx(z=0)Si(0)(x)dx.

Taking into account that the integral in the numerator of Eq. (21) represents the Fourier transform operator with the spatial frequency Δx/λz , we come to the following conclusion: the spectral degree of coherence of the field component in the far zone of the completely incoherent source is equal to the normalized Fourier transform of power spectrum of corresponding field component across the source. In the scalar coherence theory this result is known as the van Cittert-Zernike theorem.

To generalize the van Cittert-Zernike theorem for the vectorial case, we substitute first from Eqs. (19) and (20) into definition of the degree of coherence of the electromagnetic field Eq. (8), i.e.,

μ˜(z)(x1,x2)=(i(z=0)Si(0)(x)exp(i2πλzx·Δx)dx2i,j(z=0)Si(0)(x)dx(z=0)Sj(0)(x)dx)1/2.

It may be readily shown that

i,j(z=0)Si(0)(x)dx(z=0)Si(0)(x)dx=[(z=0)S(0)(x)dx]2,

where S (0)(x′) is the power spectrum of the source defined by Eq. (4). On substituting from Eq. (23) into Eq. (22), we obtain finally

μ˜(z)(x1,x2)=(i(z=0)Si(0)(x)exp(i2πλzx·Δx)dx2)1/2(z=0)S(0)(x)dx.

Equation (24) represents the vector generalization of the van Cittert-Zernike theorem, which can be formulated as follows: the spectral degree of coherence of the electromagnetic field in the far zone of the completely incoherent source is equal to the normalized root-mean-square absolute value of the Fourier transforms of power spectra of field components across the source.

Now we will examine the polarization properties of the electromagnetic field in the far zone of the incoherent source. On making use of definition (12), or alternatively Eq. (14), we obtain the following formula for the degree of polarization of the field in the far zone of the incoherent source:

P(z)(x)=4α24α+1,

where

α=(z=0)Sx(0)(x)dx(z=0)S(0)(x)dx.

Equations (25) and (26) show that the degree of polarization of the completely incoherent electromagnetic field does not change on propagation and depends only on the energy distribution between the field components in the source plane. In particular, for completely unpolarized source, when Sx (0)(x′)=Sy (0)(x′), the equality P (z)(x) = 0 takes place, i.e., the generated field remains completely unpolarized on propagation. On the contrary, whenSx (0)(x′) = 0 or Sy (0)(x′) = 0, it takes place the equality P (z)(x) = 1, i.e., the generated field remains completely polarized on propagation. In general case the vector version of the van Cittert-Zernike theorem given by Eq. (24) must necessarily be supplemented with Eq. (25).

4. Example of theorem application

As an example of possible applications of the deduced theorem we consider the problem of optical simulation of partially coherent and partially polarized secondary source with the controlled statistical properties. Such a simulation may be realized by means of the modified Mach-Zehnder interferometer sketched schematically in Fig. 1.

 figure: Fig. 1

Fig. 1 Modified Mach-Zehnder interferometer: BS - beam splitter; M - mirror; P - polarizer; OW - optical wedge.

Download Full Size | PDF

Let us consider that the completely incoherent circular source of radius R and of uniform power spectrum S (0) 0is placed in the input plane (primary source plane) of the interferometer. Let polarizers P1 and P2 be chosen to transmit only x component and only y component of the incident field, respectively, and let optical wedges OW1 and OW2 be chosen to attenuate the power spectra of the optical fields at the opposite arms of the interferometer in the ratio α: (1 - α) with α being the desired coefficient of attenuation. Then, assuming that the distance z = z 1 + z 2 + z 3 + z 4 between the primary source plane and the secondary source plane is large enough, the degree of coherence of the generated secondary source can be calculated using Eq. (24), i.e.,

μ˜(z)(x1,x2)=α2+(1α)2πR2(xR)exp(i2πλzx·Δx)dx.

Calculating the integral in Eq. (27) (see, e.g. [3]), we find

μ˜(z)(x1,x2)=2α22α+12J1(2πRλzΔx)(2πRλzΔx),

where J 1(.) is the Bessel function of the first kind and the first order.

The function 2J 1(ξ)/ξ decreases steadily from the value unity when ξ = 0 to the value zero when ξ = 0.61λz/R . Correspondingly, the radius of the area where the degree of coherence (28) exceeds 0.88 of its maximum value is approximately given by

rcoh=0.16λz2R.

As well known, this radius can serves as the numerical measure of coherence. Hence, changing the distance z (really z 4), one can control the coherence of the produced vector field over a rather wide range. On the other hand, varying the attenuation coefficient α, one can control the degree of polarization given by Eq. (25) within the range from zero to unity.

5. Conclusion

We have generalized the well known van Cittert-Zernike theorem for the case of vector electromagnetic fields. The deduced theorem shows that the degree of coherence of the electromagnetic field produced by the completely incoherent vector source increases on propagation whereas the degree of polarization remains unchanged. The possible application of the deduced theorem has been illustrated by an example of optical simulation of partially coherent and partially polarized secondary source with the controlled statistical properties.

Acknowledgments

The authors gratefully acknowledge the financial support from the National Council for Science and Technology (CONACYT) of Mexico and from the Autonomous University of Puebla (BUAP) under projects 55108 and OSA-EXT09-G, respectively. A.S. Ostrovsky thanks also the National Institute of Astrophysics, Optics and Electronics, where he realizes presently his sabbatical leave.

References and links

1. P.H. van Cittert, “Die Wahrscheinliche Schwingungsverteilung in Einer von Einer Lichtquelle Direkt Oder Mittels Einer Linse Beleuchteten Ebene,” Physica 1, 201–210 (1934). [CrossRef]  

2. F. Zernike, “The concept of degree of coherence and its application to optical problems,” Physica 5, 785–795 (1938). [CrossRef]  

3. M. Born and E. Wolf, Principles of Optics (Cambridge University Press, Cambridge, UK, 1997).

4. L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, Cambridge, UK, 1995).

5. E. Wolf, Introduction to the Theory of Coherence and Polarization of Light (Cambridge University Press, Cambridge, UK,2007).

6. A. S. Ostrovsky, P. Martínez-Vara, M. A. Olvera-Santamaría, and G. Martínez-Niconoff, “Vector coherence theory: An overview of Basic concepts and definitions,” in Recent Research Developments in Optics, S.G. Pandalai, ed. (Research Signpost, Trivandrum, Kerala, India, to be published).

7. J. Tervo, T. Setälä, and A. T. Friberg, “Degree of coherence for electromagnetic fields,” Opt. Express 11, 1137–1143 (2003). [CrossRef]   [PubMed]  

8. J. Tervo, T. Setälä, and A. T. Friberg, “Theory of partially coherent electromagnetic fields in the space-frequency domain,” J. Opt. Soc. Am. A 21, 2205–2215 (2004). [CrossRef]  

9. E. Wolf, “Correlation-induced changes in the degree of polarization, the degree of coherence, and the spectrum of random electromagnetic beams on propagation,” Opt. Lett. 28, 1078–1080 (2003). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (1)

Fig. 1
Fig. 1 Modified Mach-Zehnder interferometer: BS - beam splitter; M - mirror; P - polarizer; OW - optical wedge.

Equations (29)

Equations on this page are rendered with MathJax. Learn more.

{ E ( x , v ) } = { E x ( x , v ) E y ( x , v ) } ,
W ( x 1 , x 2 ) = [ W ij ( x 1 , x 2 ) ] = E i * ( x 1 ) E j ( x 2 ) ( i , j = x , y ) ,
S i ( x ) = W ii ( x , x ) = E i ( x ) 2 ,
S ( x ) = E ( x ) E ( x ) = S i ( x ) + S j ( x ) = Tr W ( x , x ) ,
W ij ( x 1 , x 2 ) S i ( x 1 ) S j ( x 2 ) .
μ ij ( x 1 , x 2 ) = W ij ( x 1 , x 2 ) S i ( x 1 ) S j ( x 2 ) ,
0 μ ij ( x 1 , x 2 ) 1 .
μ ˜ ( x 1 , x 2 ) = ( i , j W ij ( x 1 , x 2 ) 2 i , j S i ( x 1 ) S j ( x 2 ) ) 1 / 2 = ( Tr [ W ( x 1 , x 2 ) W ( x 1 , x 2 ) ] Tr W ( x 1 , x 1 ) Tr W ( x 2 , x 2 ) ) 1 / 2 .
μ ˜ ( x 1 , x 2 ) = ( i , j μ ij ( x 1 , x 2 ) 2 S i ( x 1 ) S j ( x 2 ) i , j S i ( x 1 ) S j ( x 2 ) ) 1 / 2 .
W ( x , x ) = W xx ( x , x ) [ 1 0 0 1 ] ,
Det W ( x , x ) = 0 ,
P ( x ) = ( 1 4 Det W ( x , x ) [ Tr W ( x , x ) ] 2 ) 1 / 2 .
2 Det W ( x , x ) = [ Tr W ( x , x ) ] 2 Tr [ W ( x , x ) W ( x , x ) ] ,
P ( x ) = 2 μ ˜ 2 ( x , x ) 1 .
E i ( z ) ( x ) = 1 i λz exp ( i 2 π λ z ) ( z = 0 ) E i ( 0 ) ( x ) exp [ i π λz ( x x ) ] d x ,
W ( z ) ( x 1 , x 2 ) = 1 ( λz ) 2 ( z = 0 ) W ( 0 ) ( x 1 , x 2 ) exp { i π λz [ ( x 2 x 2 ) 2 ( x 1 x 1 ) 2 ] } d x 1 d x 2 .
W ii ( 0 ) ( x 1 , x 2 ) = E i ( 0 ) * ( x 1 ) E i ( 0 ) ( x 2 ) = η S i ( 0 ) ( x 1 ) S i ( 0 ) ( x 2 ) δ ( x 1 x 2 )
W ij ( 0 ) ( x 1 , x 2 ) = E i ( 0 ) * ( x 1 ) E j ( 0 ) ( x 2 ) = E i ( 0 ) * ( x 1 ) E j ( 0 ) ( x 2 ) = 0 ,
W ii ( z ) ( x 1 , x 2 ) = η ( λ z ) 2 exp [ i π λ z ( x 1 2 x 2 2 ) ] ( z = 0 ) S i ( 0 ) ( x ) exp ( i 2 π λz x · Δ x ) d x
W ij ( z ) ( x 1 , x 2 ) = 0 ,
μ ii ( z ) ( x 1 , x 2 ) = ( z = 0 ) S i ( 0 ) ( x ) exp ( i 2 π λz x · Δ x ) d x ( z = 0 ) S i ( 0 ) ( x ) d x .
μ ˜ ( z ) ( x 1 , x 2 ) = ( i ( z = 0 ) S i ( 0 ) ( x ) exp ( i 2 π λz x · Δ x ) d x 2 i , j ( z = 0 ) S i ( 0 ) ( x ) d x ( z = 0 ) S j ( 0 ) ( x ) d x ) 1 / 2 .
i , j ( z = 0 ) S i ( 0 ) ( x ) d x ( z = 0 ) S i ( 0 ) ( x ) d x = [ ( z = 0 ) S ( 0 ) ( x ) d x ] 2 ,
μ ˜ ( z ) ( x 1 , x 2 ) = ( i ( z = 0 ) S i ( 0 ) ( x ) exp ( i 2 π λz x · Δ x ) d x 2 ) 1 / 2 ( z = 0 ) S ( 0 ) ( x ) d x .
P ( z ) ( x ) = 4 α 2 4 α + 1 ,
α = ( z = 0 ) S x ( 0 ) ( x ) d x ( z = 0 ) S ( 0 ) ( x ) d x .
μ ˜ ( z ) ( x 1 , x 2 ) = α 2 + ( 1 α ) 2 π R 2 ( x R ) exp ( i 2 π λz x · Δ x ) d x .
μ ˜ ( z ) ( x 1 , x 2 ) = 2 α 2 2 α + 1 2 J 1 ( 2 π R λ z Δ x ) ( 2 π R λ z Δ x ) ,
r coh = 0.16 λz 2 R .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.