The deformation, flicker, and drift of a light field owing to complex media such as a turbulent atmosphere have limited its practical applications. Thus, research on invariants in randomly fluctuated light fields has garnered considerable attention in recent years. Coherence is a statistical property of light, while its full and quantitative characterization is challenging. Herein, we successfully realize the orthogonal modal decomposition of partially coherent beams and introduce the application of coherence entropy as a global coherence characteristic of such randomly fluctuated light fields. It is demonstrated that coherence entropy remains consistent during propagation in a unitary system by unraveling complex channels. As representative examples, we study the robustness of coherence entropy for partially coherent beams as they propagate through deformed optical systems and turbulent media. Coherence entropy is anticipated to serve as a key metric for evaluating the propagation of partially coherent beams in complex channels. This study paves the way for a broader application scope of a customized low-coherence light field through nonideal optical systems and complex media.

- Advanced Photonics
- Vol. 6, Issue 4, 046002 (2024)
Abstract
Keywords
Video Introduction to the Article
1 Introduction
In contrast to amplitude, phase, and polarization, coherence is a statistical property of photons that characterizes their phase space occupation.1
However, the full characterization of incoherent light fields is challenging. A four-dimensional cross-spectral density or mutual coherence function is required to describe the coherence properties in the quasi-monochromatic case.1
In this work, we introduce the utilization of coherence entropy, a global coherence characteristic of randomly fluctuated light fields, in the analysis of partially coherent beams generated from orthogonal modes.26
Sign up for Advanced Photonics TOC. Get the latest issue of Advanced Photonics delivered right to you!Sign up now
2 Theory
A partially coherent light field, which is a randomly fluctuated electric field, is the statistical average of instantaneous electric fields. The cross-spectral density function aids in analyzing second-order statistics,2 such as the average intensity and degree of coherence when it propagates through complex media. In practical applications, a single slice of four-dimensional cross-spectral density can be measured.23
We now consider a one-sided unitary optical system with a transmission matrix (which acts only on spatial distributions), after which the cross-spectral density of the output partially coherent light is expressed as
To quantitatively evaluate the changes of mode-weights in complex media, the coherence entropy was introduced, and defined as26
The coherence entropy ranges from a minimum value of 0, indicating a fully ordered state where only one mode possesses a weight of 1, to a maximum value of 1, representing the most disordered state characterized by multiple modes with equal weights and no bias. Thus, the value of coherence entropy is contingent upon the relative distribution of mode-weights. It can be utilized to evaluate the robustness of partially coherent light transmission in optical systems and to detect variations in the light field induced by the system. Entropy is a fundamental concept that quantifies the level of disorder. Recent research has explored the connection between information entropy and the spatial coherence width of a partially coherent beam.16 The presence of disordered phases in a dielectric nanofin array leads to varying sizes of instantaneous intensity speckles. Consequently, a greater amount of information entropy is associated with reduced spatial coherence. Similarly, in the context of modal decomposition, the coherence entropy tends to increase as the coherence width decreases.
3 Results
To experimentally investigate the coherence entropy during propagation, multimode ptychography was employed to reconstruct the full information of the cross-spectral density function, that is .31 represents the measured mixed modes that are not necessarily orthogonal, and denotes the total number, which is not necessarily equal to . To further calculate the mode-weight for the chosen orthogonal basis, e.g., HG basis, each measured mode was decomposed as . Thus, the cross-spectral density in Eq. (2) was reorganized as
Therefore, the mode-weight in Eq. (2) was calculated using
Figure 1 shows a schematic of orthogonal modal decomposition for a GSM beam (see Appendix), whose coherence depends on the beam spot size on the rotating ground glass disk (RGGD). As shown in Part-3, the mixed modes were reconstructed according to the multimode ptychography iterative engine, whose update function follows Eq. (10). Examples of the reconstructed mixed modes for a low-coherence GSM beam and the diffraction patterns are shown in Fig. 1. Considering a GSM beam and two mode-superposed partially coherent light sources as examples, Fig. S1 in the Supplementary Material shows that the calculated mode-weight of an orthogonal basis (shown in bars) agrees with that actually used for superposition. If only the traditional two-dimensional focal cross-spectral density is measured, the mode-weights cannot be obtained, except for beams with known models and analytical solutions. However, the four-dimensional cross-spectral density measurement results of our modal decomposition method yield complete information and exhibit a high flexibility.
Figure 1.Schematic diagram of incoherent modal decomposition. Part-1: generation of a partially coherent beam with adjustable coherence width. The coherence width can be modified by shifting the spherical lens before the RGGD back and forth. Part-2: intensity monitor. The GSM beam propagates through a spherical lens system (
To demonstrate the recovery of mode-weights in deformed systems and the significance of the measurement scheme, we examined a GSM beam passing through a cylindrical lens system. The schematic of this system is shown in Fig. 1(a). In contrast to the GSM beam passing through a spherical lens, the distribution of focal intensity has been elongated after it passes through the cylindrical lens [Figs. 2(a) and 2(b)], and the deformation also occurs on the cross-spectral density distribution. Here, “” represents the transmission matrix of the cylindrical lens system. If original symmetric HG modes [Fig. 2(c)] defined according to the beam width and coherence width along the short axis are used, the mode-weights of the focused elongated beam [Fig. 2(e)] are not consistent with those of the source GSM beam [Fig. 2(d)]. This is in good agreement with the theoretical analysis (see Supplementary Material). However, when we applied the transmission matrix on the source orthogonal basis to obtain a new spatial basis [Fig. 2(f)], the mode-weight distribution can be recovered.37,38 We can capitalize on such flexibility of the measurement scheme for broader applications. Figures 2(g)–2(j) show the mode-weight adjustment of a GSM beam with various degrees of coherence. Here, represents coherence length and represents the beam width. The “before” shown with orange bars means that the mode-weights were calculated according to a symmetrically distributed HG basis, and the “after” (green bars) means that the mode-weights were adjusted by applying the transmission matrix on the source orthogonal basis. After basis adjustment, the mode-weight can be corrected to be consistent with that of the source plane for both high- and low-coherence cases.
Figure 2.Basis adjustment for modal decomposition. (a) A GSM beam in the source plane and (b) the elongated beam spot after going through a cylindrical lens system. Using the symmetric HG basis, mode-weight results in panels (d) and (e) are calculated for (a) and (b), respectively. (f) Basis after applying the transmission matrix on panel (c), which helps recover the mode-weights. (g)–(j) Mode-weight distributions of a GSM beam propagating through a cylindrical lens for various coherent widths. The orange bars correspond to the mode-weight distribution before adjusting the orthogonal basis, and the green bars show the mode-weights after adjusting the orthogonal basis. Here, only the main part (
To further illustrate that coherence entropy can be recovered in deformed systems, theoretical mode-weights and their corresponding coherence entropies were simulated for a GSM beam passing through a cylindrical lens system, as shown in Table 1. The row labeled “GSM source” shows the dependence of coherence entropy on coherence width. The lower the coherence width, the higher the coherence entropy. After passing through the cylindrical lens system, if the spatial basis is not suitable, the computed coherence entropy will deviate from the coherence entropy of the source. The last row shows the results obtained after basis adjusting by applying the transmission matrix on the source basis. It shows that the coherence entropy can be recovered.
GSM source | 0.1094 | 0.2657 | 0.5142 | 0.7777 |
With inappropriate basis | 0.3615 | 0.4424 | 0.5816 | 0.7645 |
After basis adjustment | 0.1092 | 0.2640 | 0.5139 | 0.7777 |
Table 1. Theoretical simulation of GSM coherence entropy before and after basis adjustment.
In comparison to static systems, the complexity of the transmission matrix in turbulent media is heightened. When a GSM beam passes through the turbulent atmosphere, the intensity captured by the camera drifts and deforms constantly. Figure 3 shows the mode-weight results of a GSM beam with three coherence widths passing through a turbulent atmosphere generated using a temperature-adjustable hot plate. The low-turbulence atmosphere in Figs. 3(a)–3(c) corresponds to room temperature, and the high-turbulence atmosphere in Figs. 3(d)–3(f) corresponds to 200°C. The coherence length was determined to be for high coherence, for medium coherence, and for low coherence (see Appendix). Figures 3(g)–3(i) show the difference between the mode-weight matrices of these two turbulence cases. As the coherence width decreases, the mode-weight profile expands, owing to the broadening of the beam and degradation in coherence induced by turbulence. However, a comparison of the first and second rows shows that the fluctuation of mode-weights decreases with the coherence.
Figure 3.Coherence entropy of a GSM beam passing through the turbulent atmosphere. The mode-weights for a GSM beam propagating through (a)–(c) low turbulence (
Figures 3(j)–3(l) provide a quantitative analysis of the propagation robustness by assessing the coherence entropy of a GSM beam as it propagates through turbulent atmosphere. Three coherence and three turbulence cases were compared. For each parameter setting, 10 sets of data were analyzed. The bars in blue, yellow, and red correspond to the mean values of coherence entropy, with the hot plate set at room temperature (low turbulence), 100°C (medium turbulence), and 200°C (high turbulence). The corresponding beam drift in Fig. S5 in the Supplementary Material shows the effect of the turbulent atmosphere generated using the hot plate. For the high-coherence case, the beam may drift out of the region defined by the beam width at a turbulence of 200°C. As shown in Figs. 3(j)–3(l), for the low-turbulence case (blue bars), as the coherence decreases, the height of the blue bar increases, which indicates an increase in coherence entropy. For medium-turbulence (yellow bars) and high-turbulence (red bars) cases, the coherence entropy is higher than that for the low-turbulence case. This is attributed to the broadening of mode-weights. Although the coherence entropy increases as the turbulence becomes stronger, for a low-coherence beam, the coherence entropy is nearly constant for various turbulence cases. Furthermore, the standard deviation of coherent entropy is represented by the black lines on the bars. The stronger the turbulence, the larger the standard deviation. Compared with the high- and medium-coherence cases, the coherence entropy for a low-coherence beam is more robust to the turbulent atmosphere.
As illustrated in previous research,39 for a Schell-model light source passing through a random medium (linear, isotropic, and homogeneous), the spectral density is the convolution of the three terms, written as . Here, means convolution, is the spectral density of the coherent portion at a distance , and and are the Fourier transforms of the source degree of coherence and second-order complex phase correlation function of the random medium, respectively. The size of kernels of and can be estimated from Figs. 3(c) and 3(d), respectively. The kernel of can be neglected when the size is considerably smaller than the kernel . It corresponds to a light source with lower coherence, where the impact of turbulence is mitigated to some extent due to the discrepancy between the broadening caused by turbulence and the profile of the mode-weights. Similar results can be seen in Figs. 2(g)–2(j), and the mode-weights of low-coherence light are more robust to complex systems than that of high-coherence light before basis adjustment. However, a trade-off must be considered in this context: as coherence decreases, the number of modes in the light fields increases, leading to increased complexity in modal decomposition experiments.
When the turbulent transmission matrix is determined or measured, the examination of mode-weights and coherence entropy in turbulent media becomes more flexible through the use of basis adjustment, facilitating applications in information coding and robustness evaluation. Additionally, to showcase the versatility of the proposed modal decomposition approach across different bases, experiments involving orbital angular momentum (OAM) multiplexing,18,40,41 and decomposition were conducted. With traditional demultiplexing methods, the broadening of the OAM spectrum induced by propagation through complex media has not been sufficiently addressed. To solve this problem, we demonstrate the improvement in OAM-based optical encryption and decryption through complex media using basis adjustment (Fig. 4). The multifocal array consisting of various integer vortices42 is widely used for OAM demultiplexing. A comparison of the regular OAM demultiplexing in free space based on incoherent decomposition with Laguerre–Gaussian basis (see Appendix) and a multifocal array is given in Figs. 4(a) and 4(b). The phase distortion induced by turbulent media [see Fig. 4(c)] was simulated using Zernike polynomials.43 Consequently, the OAM demultiplexing through turbulent media became chaotic for both cases, as shown in Figs. 4(d)–4(f). The application of the unitary transmission matrix incorporating the perturbation phase on the source orthogonal bases results in the preservation of orthogonality in the output bases, as demonstrated in Fig. S6 in the Supplementary Material. Subsequently, OAM demultiplexing after propagation through turbulent media can be achieved through incoherent decomposition. For example, an input binary number series “001101011-011100111” was decrypted based on a threshold of 0.25 when reading the normalized mode-weight distribution in Fig. 4(f). Similarly, a string of clock-style numbers “215006” was encrypted using OAM multiplexing; its coding correspondence is shown in Fig. 4(g). Traditional multifocal arrays cannot decrypt such information [see Fig. 4(h) and corresponding simulation results in Fig. S7 in the Supplementary Material]. However, the decryption can be successfully performed after applying basis adjustment in incoherent decomposition.
Figure 4.Performance improvement in OAM-based optical encryption and decryption through turbulent media using basis adjustment. OAM demultiplexing in free-space using (a) incoherent decomposition and (b) traditional multifocal array consisting of different integer vortices. (c) Wavefront distortion induced by the turbulent media. Chaotic OAM demultiplexing through turbulent media using (d) orthogonal modal decomposition and (e) traditional multifocal array. (f) Corrected OAM demultiplexing through turbulent media using mode adjustment for incoherent decomposition. (g) Encryption of a clock-style number string “215006.” (h) Experimental demonstration showing that the correct decryption may be achieved after applying mode adjustment in incoherent decomposition, while the traditional method yields highly inaccurate results.
Furthermore, the coherence entropy proposed in this study defined using an OAM mode-weight matrix may prove instrumental in OAM-based optical communications through complex media.44,45 In addition, the evaluation of random fluctuations induced by the light source or optical system is no longer limited to the beam shape or other two-dimensional attributes. To demonstrate this, we performed a channel robustness experiment based on the conservation of coherence entropy (Fig. 5). Mode-weights depicted in Fig. 5(a) and Figs. 5(b) and 5(c) were calculated for a focal system in the absence and presence of turbulent media, respectively. In Figs. 5(a) and 5(b), the decomposition bases utilized on the output plane were determined by applying only the transmission matrix of focused propagation on the Laguerre–Gaussian bases. Conversely, in Fig. 5(c), the turbulence phase was taken into account in the transmission matrix and basis calculation. The inset intensity patterns were captured at a focal plane that showed the distortion effects owing to the turbulent media. After applying basis adjustment, the mode-weights were corrected from Figs. 5(b) and 5(c); the reconstructed intensity is shown in Fig. 5(c) (see corresponding HG decomposition in Fig. S8 in the Supplementary Material). In Fig. 5(d), the measured coherence entropy without turbulence is highlighted using a black line. A blue dashed line and orange circles correspond to the coherence entropy through turbulent media and corrected coherence entropy after basis adjustment, respectively. A simple demonstration of encoding, which leverages the conservation coherence entropy, is shown in Figs. 5(e)–5(g). The correspondence between colors and the values of coherence entropy is shown in Fig. 5(d). Given that the intensity is highly chaotic after a light beam passes through turbulent media, the input values of coherence entropy drastically increased. Consequently, the color information of architectural images was scrambled [Fig. 5(f)]. However, owing to the conservation of coherence entropy, the correct image could be reconstructed, as shown in Fig. 5(g). The retrieval of information benefits from the precise feedback of the channel transmission matrix. Consequently, when the system encounters unforeseen perturbations, the coherence entropy will fluctuate, indicating the instability of the system.
Figure 5.Channel robustness evaluation using the coherence entropy. Mode-weight distribution (a) in free space and (b) through turbulent media. Inset intensity patterns show the distortion effect caused by turbulent media. The mode-weights and corresponding reconstructed intensity are shown in panel (c).
4 Discussions and Conclusions
This study examines the coherence entropy as defined by the mode-weights of orthogonal bases within partially coherent light fields. The use of orthogonal modal decomposition allows for an exploration of the propagation robustness of coherence entropy within a unitary system, encompassing both static and turbulent media. Our findings demonstrate the recovery of coherence entropy through the unraveling of complex channels and the selection of appropriate bases. The recovery of both mode-weight and coherence entropy is closely linked to the principle of energy conservation. Therefore, in the event that absorption is distributed evenly among all modes, the relative distribution of mode-weights and corresponding coherence entropy will remain unchanged. Conversely, if absorption is not uniform, the coherent entropy will be altered. Given the susceptibility of coherent entropy to changes in mode weights and basis selection, this metric can be utilized to evaluate the robustness of the optical system. The existing approach for detecting coherence entropy is reliant on the chosen basis, and the development of a basis-independent detection method could enhance the applicability of coherence entropy in various fields.
However, it is noteworthy that the coherence entropy, as defined by the mode-weights of a nonorthogonal basis, may not necessarily be retrievable. An interesting example is the van Cittert–Zernike theorem (no energy loss), wherein the modes undergo a Fourier transform, i.e., a typical unitary channel. We discuss an incoherent source with equal mode-weights, i.e., . Following the Fourier theorem, on the output plane, a partially coherent light beam can be obtained (nonequal mode-weights). This result may indicate a decreased entropy in a unitary system. However, each mode on the output plane corresponds to a plane wave with a certain tilted phase, which is not an orthogonal mode. Thus, the entropy defined according to these mode-weights is not included in our analysis. The analysis of coherence entropy can be extended to a vector beam with diagonal beam coherence-polarization matrix based on the assumption that the two initial polarization components share the same set of orthogonal spatial basis (see Supplementary Material). The modal decomposition of diagonal coherence-polarization matrix can be simplified as that of a scalar cross-spectral density defined with .33 Thus, the overall mode-weights and corresponding coherence entropy for such a vector partially coherent beam in a unitary system can be reconstructed by selecting appropriate decomposition bases.
Owing to the random fluctuations in turbulence, the coherence entropy carried by high-coherence light increases as it propagates through these media. The variation in coherence entropy increases with the turbulence strength. Without basis adjustment, light fields with lower coherence are more robust to turbulence in terms of coherence entropy. Turbulence can also be described as a set of phase screens and propagation (slicing up a thick medium), both of which are unitary. The product of unitary matrices is also unitary. Then, the application of the transmission matrix to the selection of orthogonal bases offers a more flexible and accurate approach to studying the robustness of coherence entropy in turbulent media. The coherence entropy, as a global characteristic, is expected to become a key metric for evaluating the propagation of partially coherent beams through complex media with broader practical application prospects.
5 Appendix: Methods
5.1 Theoretical Calculation of Mode-Weights
According to the literature,1,2 the cross-spectral density function may be decomposed as . Here, represents the mode-weights of the mode, represents eigenfunctions orthogonal to each other, and is the number of orthogonal eigenfunctions. Considering a GSM beam as an example, by solving the homogenous Fredholm integral equation , the cross-spectral density may be decomposed into the incoherent superposition of HG modes36 and
5.2 Experimental Generation of a Partially Coherent Light Source
As shown in Fig. 1, Part-1 shows the generation of a GSM beam according to the adjustable coherent width. The coherent light emitted from a semiconductor laser (VentusHR, wavelength ) is expanded by a beam expander (BE) and then focused by a spherical lens (focal length ) onto the RGGD. The scattered light is collimated by the second spherical lens () and shaped by a Gaussian amplitude filter (GAF, 1-mm beam waist). After the GAF, a GSM beam with an adjustable coherence is obtained. The coherence width may be modified by shifting the first spherical lens before rotating the RGGD back and forth. The larger the beam spot, the lower the coherence. For convenience, the GAF can also be replaced with a spatial light modulator or digital micromirror device. Part-2 shows the intensity monitor of GSM propagating through spherical lens system (). The first camera (CCD1, ECO445 SVS-VISTEK, , ) is located on the focal plane of the spherical lens. The coherence widths in Fig. 3 correspond to (high coherence), (medium coherence), and (low coherence).
5.3 Experimental Measurement of Mode-Weights
In Fig. 1, Part-3 is used to measure the mixed modes of the GSM beam. The GSM beam that is to be measured illuminates a high-scattering object (1951 USAF resolution test chart), and the second camera (CCD2, Electron-Multiplying CCD, iXon Life Oxford, , ) is placed at a distance of 146 mm after the object to take the diffraction patterns. To obtain a higher volume of information for reconstruction, a two-dimensional translation stage (CONEX-MFACC Newport) is used to scan the USAF for ptychographic data acquisition. The focal spot is approximately wide, and the step size is . Finally, the mixed modes and sample information are reconstructed using a multimode ptychography engine.31 The cross-spectral density and mode-weights may then be calculated further.
For the partially coherent beam with a known model, e.g., GSM, the mode-weights can be calculated using the analytical solution, e.g., Eq. (7). Here, we introduced a model-independent mode-weight measurement method that is suitable for more general partially coherent beams. The characteristics of the partially coherent beam are evident once the propagation or scattering occurs. Therefore, coherent diffraction imaging is used to reconstruct the complete four-dimensional cross-spectral density function of a partially coherent beam. To improve the redundancy of information in the diffraction pattern, which is the key to the convergence of iterative algorithm, a high-scattering object (USAF) is placed between the beam and camera (see Fig. 1). In the conventional coherent diffraction imaging (CDI) device, the beam acts as a probe, and the object is the target to be reconstructed. Here, although the target of the measurement is the light source, we can still use the procedure of CDI. Given that the oversampling ratio should be larger than 2,47 one diffraction pattern is not sufficient to reconstruct the mixed modes of the GSM. Therefore, a ptychographic scanning is performed.
The purpose of the CDI iterative algorithm is to minimize the difference of collected intensity and calculated average-intensity or the output field right behind the object plane. For instance, in the coherence case (i.e., single-mode case), one must solve . Here, is the probe, and is the object. The partially coherent beam is the superposition of multimodes, and thus, the problem should be transformed into minimizing collected intensity and summation of calculated intensity of each mode. Subsequently, each single diffraction field can be renewed with as48
Regarding the error of modal decomposition, the upper limit of the sum should tend to infinity. For the partially coherent beam with a known model and analytical solution of mode-weights or the mode-superposed beam, the orthogonal basis may be easily determined. However, for a totally random partially coherent beam, the basis is difficult to choose, and the error of finite decomposition modes increases. Nevertheless, we still emphasize that the scheme of the proposed modal decomposition is flexible, especially when dealing with beams superimposed by orthogonal modes.
Xingyuan Lu is a postdoctoral fellow at Soochow University, collaborating with Prof. Chengliang Zhao. She received her BSc and PhD degrees from Soochow University in 2014 and 2021, respectively. From 2019 to 2021, she worked as a visiting PhD with Prof. John (Jianwei) Miao at the University of California, Los Angeles. Her main research interests are light manipulation, optical field with topology/vortex, and coherent diffraction imaging.
Zhuoyi Wang received his MSc degree in the School of Physical Science and Technology at Soochow University, China, in 2022. He is now a PhD candidate in the School of Physical Science and Technology at Soochow University, under the supervision of Prof. Chengliang Zhao. His research is mainly focused on the development of statistical optics and orbital angular momentum with their applications.
Qiwen Zhan received his PhD in electrical and computer engineering from the University of Minnesota. He is currently the principal investigator of Nano-photonics Research Group at the University of Shanghai for Science and Technology.
Yangjian Cai is a professor at the School of Physics and Electronics, Shandong Normal University, Jinan, China. He received his BSc in physics at Zhejiang University, his PhD in physics at Zhejiang University, and his PhD in electromagnetic theory at Royal Institute of Technology. His research interests include light manipulation and applications, optical coherence and polarization, atmospheric optics, laser physics, and optical imaging.
Chengliang Zhao is a professor at the School of Physical Science and Technology, Soochow University, China. He received his PhD in physics from Zhejiang University. His research interests include coherent optics, diffractive imaging, phase retrieval, and optical tweezers.
References
[1] E. Wolf. Introduction to the Theory of Coherence and Polarization of Light(2007).
[4] J. W. Goodman. Statistical Optics(2015).
[8] A. M. Paniagua-Diaz et al. Blind ghost imaging. Optica, 6, 460-464(2019).
[27] K. Kim, D. Y. Park, J. G. Kim. Entropic measure of global coherence of a fluctuating field. J. Korean Phys. Soc., 35, 186-189(1999).
[30] C. Okoro et al. Demonstration of an optical-coherence converter. Optica, 4, 1052-1058(2017).
[32] F. Gori et al. Beam coherence-polarization matrix. Pure Appl. Opt., 7, 941(1998).
[35] B. L. Moisewitsch. Integral Equations(2011).
[38] M. Mazilu et al. Random super-prism wavelength meter. Opt. Lett., 39, 96-99(2014).

Set citation alerts for the article
Please enter your email address