Issue 
J. Eur. Opt. SocietyRapid Publ.
Volume 18, Number 2, 2022



Article Number  12  
Number of page(s)  4  
DOI  https://doi.org/10.1051/jeos/2022012  
Published online  06 December 2022 
Research Article
Singular value representation of the coherence Poincaré sphere
^{1}
Institute of Photonics, University of Eastern Finland, P.O. Box 111, 80101 Joensuu, Finland
^{2}
Department of Physics, University of Miami, Coral Gables, FL 33146, USA
^{*} Corresponding author: jyrki.laatikainen@uef.fi
Received:
1
February
2022
Accepted:
7
November
2022
The socalled coherence Poincaré sphere was recently introduced for geometrical visualization of the state of twopoint spatial coherence of a random electromagnetic beam. The formalism and its interpretation strongly utilized a specific decomposition of the Gram matrix of the crossspectral density (CSD) matrix. In this work, we show that the interpretation of the coherence Poincaré sphere is obtained exclusively and straightforwardly via the singular value decomposition of the CSD matrix.
Key words: Electromagnetic beams / Optical coherence / Polarization
© The Author(s), published by EDP Sciences, 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Optical coherence theory deals with the analysis and consequences of randomness in optical fields which gives rise to partial (spatial and temporal) coherence of light [1, 2]. A topic of substantial interest in the last two decades has been the coherence properties of vectorial light [3], and among the very recent results is the geometrical representation called the coherence Poincaré sphere [4, 5]. This formalism displays the spatial coherence of a random electromagnetic beam similarly as the traditional Poincaré sphere [6, 7] depicts the beam’s polarization characteristics, and it is the first graphical representation of electromagnetic twopoint coherence. Illustrative geometrical interpretations of this kind are often found to be extremely useful in physics as evidenced by the conventional Poincaré sphere in polarization optics and the Bloch sphere in quantum mechanics [8]. In optics, the Poincaré sphere and its variants have found important applications, e.g., in the context of full Poincaré beams [9], orbital angular momentum [10], higherorder polarization states [11, 12], vector fields [13], and scalar twobeam interference [14].
In the previous works [4, 5], the derivation and interpretation of the coherence Poincaré sphere was based on decomposing the Gram matrix of the crossspectral density (CSD) into two parts in full analogy to the division of the polarization matrix [1] into parts corresponding to a completely unpolarized and fully polarized beams. The present work complements and extends these earlier studies by employing the singular value decomposition (SVD) of the CSD matrix, which was utilized in [5] and more extensively studied in [15]. More precisely, we derive the formalism of the coherence Poincaré sphere using the SVD exclusively, and show that this approach leads to a physical interpretation for the sphere as a geometric representation of the intertwined coherence and polarization information conveyed by the singular values and vectors of the CSD matrix.
2 Discussion
Consider a random, polychromatic, and statistically stationary electromagnetic beam field. The spatial coherence properties of the field at two positions r _{1} and r _{2} on a transversal plane with respect to the propagation direction are described in the spacefrequency domain by the CSD matrix [1, 2, 6, 16](1)taken as an average over an ensemble of monochromatic (transverse, two component) electric field realizations E(r, ω) at angular frequency ω. The angle brackets, asterisk, and superscript T stand for the ensemble average, complex conjugate, and matrix transpose, respectively. The SVD of the CSD matrix is written as [5–7](2)where the dagger denotes Hermitian conjugation and W _{12} = W . From now on we do not explicitly show the frequency dependence of various quantities. In equation (2), and are unitary matrices, and D = diag[ν _{+},ν _{−}], with ν _{+} and ν _{−} representing the singular values of the CSD. The complex unit vectors and in the columns of U and V are the left and right singular vectors of the CSD matrix obeying(3) (4)
The singular values are real and satisfy . Furthermore, their squares coincide with the eigenvalues of the Gram matrices of the CSD and its Hermitian adjoint given, respectively, by [4, 5](5) (6)
where Ω _{12} = Ω , Ω _{21} = Ω , and we employed the quasiHermiticity property . We note that the Gram matrices contain secondorder coherence information only. The vectors and fulfil the eigenvalue equations(7) (8)
Expressions for can therefore be obtained from the characteristic equation , where σ _{0} is the 2×2 unit matrix and det denotes the determinant. The characteristic equation can be written as(9)where tr stands for the trace. The (squared) singular values are then obtained as(10)where(11)is bounded as 0 ≤ P _{Ω}(r _{1},r _{2}) ≤ 1.
Next, we highlight some central properties of Ω _{12} and Ω _{21}. Firstly, they are Hermitian and nonnegative definite matrices that satisfy the conditions , , , det Ω _{12} = det Ω _{21} ≥ 0, and contain nonnegative diagonal entries. Secondly, the coherence information in matrix Ω _{12} is generally different from that in Ω _{21} due to the quasiHermiticity of the CSD matrix as is evident from equations (5) and (6). The mathematical properties of Ω _{12} and Ω _{21} are similar to those of the polarization matrix and they can formally be decomposed into two parts, one of which is proportional to the identity matrix and the other has zero determinant. This division is analogous to the decomposition of the polarization matrix into parts corresponding to a completely unpolarized beam and a fully polarized beam [1]. The earlier works concerning the coherence Poincaré sphere [4, 5] were extensively based on this division. Here we follow a different procedure and interpret the sphere using the SVD of the CSD only.
We proceed by defining the Stokes parameters of Ω _{12} as(12)where σ _{ j }, with j ∈ (1, 2, 3), are the Pauli spin matrices [1]. We remark that analogous definitions hold naturally for Ω _{21}. These parameters are realvalued and contain information on the twopoint spatial coherence of the beam. The Stokes parameters can be normalized as(13)with S _{0}(r) = trW(r, r, ω) being the spectral density of the beam. The parameters q _{ j }(r _{1},r _{2}), j = 1, 2, 3, obey the quadratic equation(14)where(15)is the electromagnetic degree of coherence [3, 17]. This degree is bounded as 0 ≤ μ(r _{1},r _{2}) ≤ 1, with the lower and upper bounds corresponding to complete incoherence and full coherence of the beam at points r _{1} and r _{2}, respectively.
Next we approach the geometric interpretation of the coherence Poincaré sphere in equation (14) via the singular value decomposition of the CSD. We first define the coherence Poincaré vector(16)that displays the spatial coherence information in Ω _{12} as points on or within a unit sphere in the (q _{1}, q _{2}, q _{3}) space. We remark that since the information content of Ω _{12} is in general different from that of Ω _{21}, two coherence Poincaré vectors q _{12} = q(r _{1},r _{2}) and q _{21} = q(r _{2},r _{1}) are required to display the spatial coherence of the beam. Analytical expressions of these vectors are obtained from the SVD of the CSD matrix as is shown below. For this purpose we recall the unitarity conditions U ^{†} U = V ^{†} V = σ _{0}, which together with the SVD and equations (5) and (6) yield(17) (18)
Furthermore, unitarity of U and V implies that and similarly for . These together with equations (10), (17), and (18) result in(19) (20)where we have written P _{Ω} = P _{Ω}(r _{1},r _{2}). It is important to note that these two expressions coincide with the decompositions in equation (3) of [4] and equation (8) of [5] which constituted the starting point of the mentioned works without a reference to the SVD. Combining the equations above with the definitions in equations (12), (13), (15), and (16), we find that(21) (22)where μ = μ These expressions provide the singularvalue interpretation of the two coherence Poincaré vectors that represent the state of spatial coherence of a partially coherent and partially polarized electromagnetic beam. Both vectors have the same length, q _{12} = q _{21} = P _{Ω} μ ^{2}, and their directions are specified by the vectors and . In addition, we note that the equalities and are obtained from equation (10), and by using them together with equations (11) and (15) we see that . Hence, the length of the coherence Poincaré vectors can be viewed as the intensitynormalized distance between the squared singular values and of the CSD matrix, and their directions are specified by the singular vectors and corresponding to the larger singular value ν _{+}.
Next, we elucidate the physical meaning of the formalism. Firstly, for a completely coherent beam the degree of coherence equals unity, which yields [4] and hence the vectors and are unitlength vectors. Fully coherent beams are thus located on the surface of a unit sphere in the space. Secondly, the origin is preserved for beams with or . The former includes the socalled pure unpolarized beams [18] and beams that can be transformed into such by a suitable unitary operation [5]. The latter naturally means that the beam is spatially fully incoherent.
For a fully polarized but spatially partially coherent beam holds [4], and the lengths of the coherence Poincaré vectors depend only on the degree of coherence as . In addition, we note that the CSD matrix of a beam with an arbitrary state of full polarization can be written as , where is a correlation function over an ensemble of random scalars , and , , are the deterministic Jones vectors that specify the polarization state of the beam at positions . As a consequence, the singular vectors are of the form and . This implies that the coherence Poincaré vectors are expressible as and , where represent the polarization Poincaré vectors at with , , being the (normalized) polarization Stokes parameters, . In other words, for a fully polarized beam the directions of the coherence Poincaré vectors depict the state of polarization as in the context of the polarization Poincaré sphere. The coordinate axes in the space represent the states of , , , righthand, and lefthand circular polarization whereas elsewhere the beam is elliptically polarized. The vector points out the polarization state of the beam at and does so at . Furthermore, if the state of polarization is uniform across the beam, these vectors converge into a single coherence Poincaré vector whose orientation expresses the uniform polarization state and length displays the squared degree of coherence at a pair of points. Finally, we observe that in a single point, , the quantities reduce to the polarization Stokes parameters , , and the vector is the polarization Poincaré vector. Consequently, the traditional polarization Poincaré sphere [6, 7] is encountered when the formalism is applied at a single point. The various reductions described above are illustrated graphically in Figure 1.
Figure 1 Poincaré sphere of electromagnetic spatial coherence. For an arbitrary partially polarized and partially coherent beam two coherence Poincaré vectors and with the same length but different orientations are generally needed (Sphere A). If the beam is fully polarized the lengths of these vectors display the degree of coherence, , and their directions show the polarization state of the beam at points and , respectively (Sphere B). If the beam is uniformly polarized, these vectors coincide and a single coherence Poincaré vector is sufficient to represent the beam, with its length again showing the degree of coherence and the direction specifying the polarization state (Sphere C). At a single point the formalism reduces to the traditional polarization Poincaré sphere where the distance from the origin is given by the degree of polarization [1] (Sphere D). 
3 Conclusion
In summary, as an extension to the previous works [4, 5], we have derived the concept of the Poincaré sphere of electromagnetic twopoint spatial coherence exclusively from the point of view of the singular value decomposition of the crossspectral density matrix. The interpretation of the concept for an arbitrary partially polarized, partially spatially coherent beam follows directly from this approach; the state of coherence of the beam is depicted by two coherence Poincaré vectors whose lengths are defined by the normalized distance between the (squared) singular values of the CSD matrix and orientations are determined by the singular vectors related to the larger singular value. Furthermore, we highlighted the interpretation of this construction for fully polarized beams for which the coherence and polarization characteristics are closely linked, and noted that at a single point the formalism coincides with the traditional polarization Poincaré sphere.
Acknowledgments
This work was supported by the Academy of Finland [projects 308393, 310511, 320166 (PREIN)]. OK thanks the Joensuu University Foundation for financial support.
Competing interests
The authors declare that they have no competing interests.
References
 Mandel L., Wolf E. (1995) Optical coherence and quantum optics, Cambridge University Press, Cambridge. [CrossRef] [Google Scholar]
 Korotkova O. (2014) Random light beams: theory and applications, CRC Press, Boca Raton. [Google Scholar]
 Friberg A.T., Setälä T. (2016) Electromagnetic theory of optical coherence [Invited], J. Opt. Soc. Am. A 33, 2431. [NASA ADS] [CrossRef] [Google Scholar]
 Laatikainen J., Friberg A.T., Korotkova O., Setälä T. (2021) Poincaré sphere of electromagnetic spatial coherence, Opt. Lett. 46, 2143. [NASA ADS] [CrossRef] [Google Scholar]
 Laatikainen J., Friberg A.T., Korotkova O., Setälä T. (2022) Coherence Poincaré sphere of partially polarized optical beams, Phys. Rev. A 105, 033506. [NASA ADS] [CrossRef] [Google Scholar]
 Brosseau C. (1998) Fundamentals of polarized light: a statistical optics approach, Wiley, New York. [Google Scholar]
 Gil J.J., Ossikovski R. (2016) Polarized light and the Mueller matrix approach, CRC Press, Boca Raton. [Google Scholar]
 Nielsen A., Chuang I.L. (2010) Quantum computation and quantum information, Cambridge University Press, Cambridge. [Google Scholar]
 Beckley A.M., Brown T.G., Alonso M.A. (2010) Full Poincaré beams, Opt. Express 18, 10777. [CrossRef] [Google Scholar]
 Padgett M.J., Courtial J. (1999) Poincaré sphere equivalent for light beams containing orbital angular momentum, Opt. Lett. 24, 439. [Google Scholar]
 Milione G., Sztul H.I., Nolan D.A., Alfano R.R. (2011) Higherorder Poincaré sphere, Stokes parameters, and the angular momentum of light, Phys. Rev. Lett. 107, 053601. [NASA ADS] [CrossRef] [Google Scholar]
 Zang X., Bautista G., Turquet L., Setälä T., Kauranen M., Turunen J. (2021) Efficient hybridmode excitation in plasmonic nanoantennas by tightly focused higherorder vector beams, J. Opt. Soc. Am. B. 38, 521. [CrossRef] [Google Scholar]
 Ren Z.C., Kong L.J., Li S.M., Qian S.X., Li Y., Tu C., Wang H.T. (2015) Generalized Poincaré sphere, Opt. Express 23, 26586. [CrossRef] [Google Scholar]
 Halder A., Norrman A., Friberg A.T. (2021) Poincaré sphere representation of scalar twobeam interference under spatial unitary transformations, Opt. Lett. 46, 5619. [NASA ADS] [CrossRef] [Google Scholar]
 Luo M., Laatikainen J., Friberg A.T., Korotkova O., Setälä T. (2022) Singularvalue decomposition and electromagnetic coherence of optical beams, Opt. Lett. 47, 5337. [NASA ADS] [CrossRef] [Google Scholar]
 Tervo J., Setälä T., Friberg A.T. (2004) Theory of partially coherent electromagnetic fields in the space–frequency domain, J. Opt. Soc. Am. A 21, 2205. [NASA ADS] [CrossRef] [Google Scholar]
 Tervo J., Setälä T., Friberg A.T. (2003) Degree of coherence for electromagnetic fields, Opt. Express 11, 1137. [NASA ADS] [CrossRef] [Google Scholar]
 Gori F., Tervo J., Turunen J. (2009) Correlation matrices of completely unpolarized beams, Opt. Lett. 34, 1447. [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Figure 1 Poincaré sphere of electromagnetic spatial coherence. For an arbitrary partially polarized and partially coherent beam two coherence Poincaré vectors and with the same length but different orientations are generally needed (Sphere A). If the beam is fully polarized the lengths of these vectors display the degree of coherence, , and their directions show the polarization state of the beam at points and , respectively (Sphere B). If the beam is uniformly polarized, these vectors coincide and a single coherence Poincaré vector is sufficient to represent the beam, with its length again showing the degree of coherence and the direction specifying the polarization state (Sphere C). At a single point the formalism reduces to the traditional polarization Poincaré sphere where the distance from the origin is given by the degree of polarization [1] (Sphere D). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.