Three-channel robust optical encryption via engineering coherence Stokes vector of partially coherent light

Optical encryption strategies utilizing fully coherent light have been widely explored but often face challenges such as speckle noise and beam instabilities. In this work, we introduce a novel protocol for multi-channel optical information encoding and encryption using vectorial spatial coherence engineering of a partially coherent light beam. By characterizing the beam’s spatial coherence structure with a 2 × 2 coherence matrix, we demonstrate independent control over the three components of the coherence Stokes vector. This allows for three-channel optical information encoding and encryption, with applications in color image representation. Unlike existing methods based on fully coherent light modulations, our approach utilizes a two-point dependent coherence Stokes vector, proving resilient to random noise in experimental scenarios. Our findings provide a robust foundation for higher-dimensional optical encoding and encryption, addressing limitations associated with partially coherent light in complex environments

Conversely, partially coherent light has been shown to effectively suppress speckle noise and reduce beam wanders and scintillations during light-matter interactions due to its low coherence feature [23][24][25][26].The rapid development of theory and experiment in optical coherence structure engineering has extended the application of partially coherent light to optical encryption [27,28].Unlike conventional methods that control the deterministic qualities of fully coherent light, optical coherence structure engineering involves the modulation of second-order statistical properties of random light [29].By tailoring the spatial coherence structure, it has been demonstrated that a partially coherent beam can exhibit extraordinary propagation features, including self-shaping, self-focusing, self-splitting, and self-reconstruction in free space [30].The optical coherence structure is now considered an efficient and unique degree of freedom for partially coherent light, finding applications in sub-Rayleigh optical imaging [31,32], robust information transmission [33,34], particle trapping [35], beam shaping [36], and optical measurement and sensing [37][38][39].
Recently, an optical encryption protocol leveraging coherence structure engineering of a scalar partially coherent light beam was proposed [28].In the protocol, the second-order spatial coherence structure of the beam serves as the information carrier for optical encryption, in contrast to deterministic characteristics like amplitude, phase, and polarization of fully coherent optical fields.It was demonstrated that coherence-based optical encryption holds an advantage over encryption methods based on fully coherent light modulation in terms of security.Additionally, the encoded information in the spatial coherence structure is more robust against complex environmental disturbances, such as optical turbulence.However, this method relies on engineering a scalar spatial coherence structure, allowing only single-channel optical encryption.Thus, it faces challenges in achieving multi-channel optical information encoding and encryption.Furthermore, the encryption of color information, crucial in applications such as vision and target recognition, remains unaddressed by coherence-based technology.
In this study, we present a protocol for multi-channel optical information encoding and encryption using vectorial spatial coherence engineering of a partially coherent vector beam.The vectorial spatial coherence structure for a beam is characterized by a 2 × 2 coherence matrix.We demonstrate that, although the elements in the coher- ence matrix are correlated, all three components in the three-dimensional coherence Stokes vector [40], derived from the coherence matrix, are uncorrelated and can be independently controlled.Thus, the coherence Stokes vector for a partially coherent light beam can be regarded as the information carrier, allowing for three-channel optical information encoding and encryption.Associating the three independent components of the coherence Stokes vector with three primary color channels-red, green, and blue-we illustrate that the coherence Stokes vector can be flexibly utilized for color image encoding and encryption.It is noteworthy that, unlike the polarization Stokes vector obtained from the polarization matrix, the coherence Stokes vector is a two-point dependent quantity, describing the second-order statistical properties of the partially coherent vector beam.Through our experiments, we demonstrate that this two-point dependent coherence Stokes vector, and consequently the encoded three-channel information, remains robust even when the partially coherent beam is subjected to random noise of varying strength.Our results pave the way for higherdimensional optical encoding and encryption in complex environments with partially coherent light.

Coherence Stokes vector and three-channel optical encryption
The second-order statistical properties of a partially coherent vector beam can be characterized by a 2 × 2 coherence matrix [41].In the space-frequency domain, the coher- ence matrix for a partially coherent vector beam, propagating along z axis, is specified by a cross-spectral density matrix, i.e., in which the elements are expressed as with (α, β) ∈ (x, y) .Above, r 1 and r 2 are two arbitrary position vectors in the transverse plane, E x (r) and E y (r) are the x and y components for the random field realizations, the asterisk denotes the complex conjugate, and the angle brackets denote the ensemble average over the field realizations.
According to the generalized van Cittert-Zernike theorem, a partially coherent beam source can be generated by propagating a spatially incoherent source.Therefore, the coherence matrix for a partially coherent beam can be expanded as [42,43] where is the polarization matrix of the spatially incoherent source with E x (v) and E y (v) being the x and y components of its field realization ( v being the position vector in the inco- herent beam source plane), and H x (r, v) and H y (r, v) are the response functions, for the x and y field components, of the optical system between the incoherent source and the partially coherent beam source.In principle, the response functions can be any wellbehaving functions, while p αβ (v) must satisfy the nonnegative condition, i.e., p xx (v) ≥ 0 , p yy (v) ≥ 0 , and From Eq. (3), it is found that by the generalized van Cittert-Zernike theorem, the polarization matrix of an incoherent source can be encoded into the coherence matrix, which is also named as vectorial spatial coherence structure of a partially coherent beam source.This finding has been used in synthesizing the partially coherent vector beams with nonconventional spatial coherence structures [44][45][46][47].In this work, we use this relation to realize the multi-channel optical information encoding and encryption.The diagonal elements p xx (v) and p yy (v) in the polarization matrix , and thus W xx (r 1 , r 2 ) and W yy (r 1 , r 2 ) in the coherence matrix, can be controlled independently.However, the anti- diagonal elements p xy (v) and p yx (v) [ W xy (r 1 , r 2 ) and W yx (r 1 , r 2 ) ] are connected tightly to the diagonal elements p xx (v) and p yy (v) [ W xx (r 1 , r 2 ) and W yy (r 1 , r 2 ) ]. Therefore, the anti- diagonal elements in the polarization and coherence matrices cannot be used to encode additional independent information as the diagonal elements. (1) To separate all the independent degrees of freedom in this encoding, we adopt the famous polarization Poincare sphere, in which the three orthogonal coordinates, S 1 (v) , S 2 (v) , and S 3 (v) can be controlled independently [48].These three orthogonal coordi- nates form the polarization Stokes vector with its elements being determined totally by the polarization matrix p(v) , i.e., Above are three Pauli spin matrices and tr denotes the matrix trace.As the polarization matrix is Hermitian, the elements in the polarization Stokes vector are real quantities.Consequently, the polarization Stokes vector can be employed to encode spatial information with real values across three distinct channels.
Akin to the polarization Stokes vector, the coherence Stokes vector for the partially coherent vector beam can be defined as [40] where the three elements are obtained from the coherence matrix, i.e., The coherence Stokes vector is a two-position dependent vector that describes the second-order statistical properties of the partially coherent vector beam [49,50].The three elements in the coherence Stokes vector have clear physical interpretations [51] such that S 1 (r 1 , r 2 ) , S 2 (r 1 , r 2 ) , and S 3 (r 1 , r 2 ) describe the differences of the electric-field x-and y-component, +π/4 -and −π/4-linearly polarized component, and right-and left-circularly polarized component correlations at points r 1 and r 2 , respectively.
Taking Eq. ( 6) into Eq.( 3) and letting the response functions H x (r, v) = H y (r, v) = H(r, v) , we obtain a generalized van Cittert-Zernike theorem for the Stokes vectors, i.e., [52] Based on the above relation, we find that the polarization Stokes vector, carrying three-channel optical information such as three independent optical images, can be encoded into the coherence Stokes vector of a partially coherent beam.Furthermore, the response function H (r, v) can contribute additional encryption keys during the encod- ing process.Consequently, the initially described three-channel optical encoding system transforms into a more comprehensive three-channel optical encryption protocol.
Figure 1 shows the schematic of our three-channel optical information encoding and encryption protocol based on the vectorial spatial coherence structure engineering of ( 5) a partially coherent light beam.The proposed protocol consists of three stages.At the first stage, three-channel optical information (plaintext) S 1 (v) , S 2 (v) , S 3 (v) , modulated in a fully coherent structured light beam, is encoded into the three two-point coherence Stokes vector elements S 1 (r 1 , r 2 ) , S 2 (r 1 , r 2 ) , and S 3 (r 1 , r 2 ) of a partially coherent ran- dom light beam via a generalized van Cittert-Zernike theorem [c.f., Eq. ( 10)].In Fig. 1a, the dynamic diffuser is used to transform a fully coherent light beam into a spatially incoherent beam, while the encoding system H is used to transform the spatially incoherent beam into the partially coherent random beam.At the second stage, the coherence Stokes vector (ciphertext) of the partially coherent beam is measured through a generalized Hanbury Brown and Twiss experiment [53,54].The elements in the coherence Stokes vector are complex in general, since the coherence matrix only follows the quasi-Hermitian property, i.e., W † (r 1 , r 2 ) = W(r 2 , r 1 ) , where the dagger denotes the Hermitian adjoint.Finally, the three-channel optical information is decoded by using the measured complex coherence Stokes vector and the inverse response function H −1 embedded with the correct decryption keys.As a specific example, we demonstrate in Fig. 1b the theoretical simulation for encoding and decoding three different images with engineering and measuring the polarization and coherence Stokes vectors of optical beams.In the simulation, the amplitude images (letters ' A' , 'B' , and 'C') is firstly loaded by the three elements of the polarization Stokes vector of a fully coherent structured beam by controlling the amplitude and phase of its two orthogonal field components E x (v) and E y (v) .Through adding random phase into the electric field, E x (v) and E y (v) become the components of the random field reali- zation for an spatially incoherent source.We then let the incoherent beam propagates through an optical system with its response function having Fourier transform form, i.e., , where A 0 is an constant, f = 250 mm is the focal dis- tance of the Fourier transform thin lens, and = 532 nm is the wavelength.In the out- put plane, the spatially incoherent beam becomes a partially coherent beam.We collect the field realizations E x (r) and E y (r) of the partially coherent beam in the output plane.The coherence Stokes vector [S 1 (r 1 , r 2 ), S 2 (r 1 , r 2 ), S 3 (r 1 , r 2 )] can be calculated by Eqs. ( 2) and ( 9).The simulated real and imaginary parts of the coherence Stokes vector elements are shown in Fig. 1b.Finally, we find the encoded three optical images are well recovered by inversely Fourier transform of the measured coherence Stokes vector elements.
One remarkable application of our three-channel optical encoding protocol is to realize the color image encoding and encryption.In optics, all colorful image can be decomposed into three primary-color-channels, i.e., red (R), green (G), and blue (B) channels.Therefore, as shown in Fig. 1c, in our protocol, we let S 1 (v) , S 2 (v) , and S 3 (v) carry infor- mation of our lab logo in R, G, B channels, respectively.By the generalized van Cittert-Zernike theorem, the colorful image can be encoded into the coherence Stokes vector of a partially coherent light beam.As shown in Fig. 1c, through measuring the coherence Stokes vector and using the inverse Fourier transform, the lab logo can be well recovered with correct colors.
It is notable that different from the traditional encryption protocols, the three-channel plaintext information in our optical-coherence-based one is encoded into the second-order field correlations rather than into the deterministic characteristics of a fully coherent light beam.Since it requires an entire statistical ensemble of a large number of field realizations to reconstruct the second-order field correlations, the multi-channel information concealed in the ciphertext of the optical-coherence-based protocol is more difficult to compromise, compared to the protocols based on fully coherent light modulation.In addition, it has been proved that the second-order field correlations are much more robust than the deterministic characteristics when the beam suffers disturbances or noisy environments, leading to a robust multi-channel optical encryption.

Synthesis of partially coherent vector beam for three-channel optical encoding
We now turn to discuss how to synthesize the partially coherent vector beam in experiment to realize the three-channel information encoding.In our experiment, the partially coherent vector beam is synthesized with the help of the complex-random-mode representation (RMR), with which the coherence matrix is expanded as the summation of N vector modes [55,56] where E n (r) and n are the vectorial complex random mode and the corresponding modal weight.According to the expression in Eq. ( 11), the partially coherent vector beam can be synthesized by superposing a set of vectorial complex random modes.In the RMR method, the modal weights for all random modes are equal to 1/N and the electric field of the random mode can be written as where τ x (r) and τ y (r) are the common amplitudes for all random modes, T x,n (r) and T y,n (r) are the complex random functions imposed on the deterministic common ampli- tudes, and êx and êy are two Cartesian unit vectors along the x and y directions.Through Eqs. ( 2), (11), and ( 12), it is found that the auto-correlations of the complex random functions T x,n (r) and T y,n (r) give rise to the diagonal elements W xx (r 1 , r 2 ) and W yy (r 1 , r 2 ) of the coherence matrix, while the cross-correlations between T x,n (r) and T y,n (r) intro- duce the anti-diagonal elements W xy (r 1 , r 2 ) and W yx (r 1 , r 2 ).
We now focus on how to encode the auto-and cross-correlations, i.e., all the information of the coherence matrix, into two complex random functions.We first obtain the polarization matrix p(v) for the fully coherent beam by using the three-channel plain- text information S 1 (v) , S 2 (v) , and S 3 (v) .Two complex random functions for the partially coherent beam are then obtained by introducing random functions into the electric field of the fully coherent beam and letting the resulting spatially incoherent beam propagates through the optical system with response function H(r, v) , i.e., where Above, r 1,n (v) , r 2,n (v) , r 3,n (v) , and r 4,n (v) are four independent random functions obey- ing Gaussian statistics with zero-mean and unit-variance, C(v) = p xy (v)/ p xx (v)p yy (v) is a complex function that controls the anti-diagonal elements in the coherence matrix of the synthesized partially coherent beam [57,58].
By adding all the complex random modes with Eq. ( 11), the partially coherent beam, with three-channel plaintext information encoded in its vectorial coherence structure, (11) can be synthesized.The RMR method used in our work is quite similar to the protocol that using the rotating ground-glass disk.Nevertheless, our method is more compact and flexible since the complex random modes E n (r) can be generated and modulated flexibly by a single spatial light modulator (SLM).

Experimental demonstration of three-channel optical encoding and encryption
We next discuss how to experimental realize the three-channel optical encoding and encryption.Figure 2a shows our experimental setup.A fully coherent laser beam with wavelength = 532 nm is emitted from a single-longitudinal-mode laser and is converted into a y-linearly polarized beam by a half-wave plate (HWP).After expanded by a beam splitter (BE 1 ), the beam splits into two parts by a beam-splitter (BS) and then they go into two different arms.The purpose of the Arm 1 in Fig. 2a is to encode the three-channel information into the coherence Stokes vector of a partially coherent light beam, while the purpose of Arm 2 is to create the reference beams for measuring the coherence properties of the synthesized partially coherent beam.
In Arm 1, the incident beam is first modulated by a phase-only SLM.The screen of the SLM is split into two halves that are used to realize the independent modulation of the x and y components for the complex random modes E n (r) .The computer-generated holo- gram (CGH) on the SLM is displayed dynamically to ensure the synthesis of partially coherent light beam.Figure 2b and c show an example for encoding the three-channel target images "1", "2", and "3" (for short, "123") into the partially coherent beam.Figure 2c displays a realization of the corresponding dynamic CGHs loaded in the SLM.The details for creating the CGHs can be found in Supplementary material 1.After modulated by the SLM, the reflected light beam goes into a 4f common-path interferometric arrangement composed by thin lenses L 1 and L 2 with the same focal distance f = 250 mm.We note the distances between the SLM and L 1 , between L 2 and the Ronchi grating (RG) are both f, and the distance between the L 1 and L 2 is 2f.The SLM is located in the input plane (i.e., front focal plane of L 1 ) of the 4f system, while a two-pinhole filter is placed in the spatial-frequency plane (rear focal plane of L 1 ) of the system to select out two +1 diffraction orders from the SLM.The beams from two diffraction orders are converted into x and y polarization states, respectively, by two different HWPs.In the output plane (rear focal plane of L 2 ), two beam modes E x,n (r)ê x and E y,n (r)ê y are synthesized by a RG.By dynamically playing the CGHs on the SLM, the partially coherent beam, with three-channel target images encoded in its coherence Stokes vector, is synthesized.
We employ the Arm 2 in the Fig. 2a to generate coherent reference light beams.The y-polarization light beam incidents onto a neutral density filter (NDF), linear polarizer (LP), and a quarter-wave plate (QWP).The axis of the LP is set to be π/4 with respect to the x axis, while the fast axis of the QWP is set to be parallel or perpendicular to the polarization direction of the beam transmitted from LP.Thus, two π/4 linearly polarized beams with a π/2 phase difference can be generated, respectively.The reference beams are then expanded by a BE 2 to produce the collimated beams of virtually uniform intensity distribution.Two reference beams are interfered with the synthesized partially coherent beam, respectively, by a polarization beam splitter (PBS).The thin lens L 3 , positioned after the RG, is used to create a 4f imaging system for imaging the synthesized partially coherent source into the CCD 1 and CCD 2 .The distances between the RG and L 3 , as well as between L 3 and CCD 1(2) , are both 2f, where f = 250 mm is the focal distance of L 3 .We then record the instantaneous intensities for the x and y field components.The CCD 1 and CCD 2 are controlled by a direct-digital-synthesis (DDS) signal generator, used to simultaneously capture the intensities of the x and y components of the composite random field.With the principle of the generalized Hanbury Brown and Twiss experiment, the real and imaginary parts of the coherence matrix can be obtained by [54] (17) , where the prime and double-prime denote the real and imaginary parts, respectively, I C1 α (r) and I C2 α (r) denote the α component intensities for the fields composed by the par- tially coherent beam and two reference beams, respectively, I α (r) is the α component intensity for the partially coherent field, while I R1 α (r) and I R2 α (r) are the α component intensities for the first and the second reference beams.The angle brackets in Eqs. ( 17) and ( 18) denote the ensemble average.We note that in the experiment, the ensemble average can be replaced with the spatial average over a single field realization if the partially coherent light beam is synthesized by a Fourier or a fractional Fourier transform system [28,59,60], i.e., the response function H(r, v) obeys Fourier or fractional Fourier transform form.The details for the generalized Hanbury Brown and Twiss experiment are shown in Supplementary material 1.After the complex coherence matrix is measured, the coherence Stokes vector can be obtained by using Eq. ( 9).
Figure 2d and e show the measured instantaneous intensities in the case when the target images "123" are encoded in the partially coherent beam, while Fig. 2f shows the three elements of the coherence Stokes vector by the spatial average over these instantaneous intensities.The corresponding simulation results for the coherence Stokes vector are shown in Fig. 2g.It is found that the measured coherence Stokes vector is consistent well with the simulation.We will show in the next section that by the inverse transform of the measured coherence Stokes vector, the encoded three-channel images can be well reconstructed.

Three-channel optical images encoding
According to the relation in Eq. ( 10), once the elements of the coherence Stokes vector S j (r 1 , r 2 ) are determined, all the information (including the shapes and positions) on the encoded three-channel target images can be recovered by using the inverse transform system whose response function is H −1 (r, v) .We first study three-channel optical images encoding with a Fourier transform system.Therefore, by evaluating an inverse Fourier transform of the measured coherence Stokes vector, the three-channel optical images S 1 (v) , S 2 (v) , and S 3 (v) can be reconstructed.In Fig. 3b-d, we show the recovered three-channel images "123" by the inverse Fourier transform of the measured coherence Stokes vector shown in Fig. 2f.The original target images are shown in Fig. 3a.The enlarged images are shown in Fig. 3e-g.We find that, indeed, both the shape and the position of the target images are well recovered from the measured coherence Stokes vector.To verify the flexibility of our method, we change the three-channel target images to the capital letters "ABC" [see in Fig. 3h] and to the "circle, square, and triangle" shapes [see in Fig. 3l], respectively.The experimentally reconstructed images are shown in Fig. 3i-k and m-o.The corresponding simulation results are shown in each inset of the figures.The experimental and simulation results verify that the three-dimensional coherence Stokes vector, indeed, can be viewed as an effective information carrier for realizing three-channel optical encoding. (18) ,

Color image encoding
Next, we demonstrate the color image encoding with partially coherent vector beam.The information of three primary-color-channels for a color image is first encoded into three elements of the coherence Stokes vector, respectively, by engineering the dynamic CGHs in the SLM.We then measure the complex coherence Stokes vector by the generalized Hanbury Brown and Twiss experiment.Through the inverse Fourier transform, the three primary-color-channels and therefore the color image can be recovered.In Fig. 4 we display the recovered results of four different target color images.The left panels of Fig. 4 show the target images, i.e., the colorful letters "RGB", a color windmill, a color peking opera mask, and our lab logo.The corresponding recovered images including the separate primary-color-channels are shown in four right panels of Fig. 4. The numerical simulation results are displayed in the inset of each figures.We find the color images is recovered, however with disturbance in the recovered color images.We note the disturbance is due to that in our experiment only a single realization of the measured intensities is used to calculate the complex coherence Stokes vector, i.e., in Eqs. ( 17) and ( 18) the ensemble average is replaced with the spatial average over a single realization.By increasing the number of the realizations in calculating the coherence Stokes vector, the recovered image quality can be further improved.In Fig. 5a-j and k-t, we show the recovered windmill and peking opera mask with different number N of intensities realizations.We can find that the recovered image quality can be improved gradually with the increase of N. The quality of the recovered color images is nearly perfect comparing to the target color images when N = 10 .We remark here that by further increase of the number of realizations used in calculating the coherence Stokes vector, the quality of the recovered image can be further improved.The quality of the recovered image depends also on the complexity of the image itself, i.e., more complex structure of the image, more number of the realizations is required to recover well the encoded image (for example, see in Supplementary Fig. S1 for the recovery of our lab logo).

Robust property
The advantage of our three-channel optical information encoding based on the coherence Stoke vector engineering is that the system is effectively resistant to the disturbances induced by the complex environments, such that when the partially coherent vector beam, with three-channel information encoded in its coherence Stokes vector, suffers random noise, the encoded information can still be well Here, we demonstrate such robustness by introducing white noise in the partially coherent beam.In our experiment, the white noise with different strengths are added into the E x,n (r) and E y,n (r) components of the complex random modes E n (r) .We note the complex random modes E n (r) are used to synthesize the partially coherent beam by the RMR method.Figure 6 shows the experimental results for the reconstructed three-channel optical images "123" and color "windmill" from the partially coherent vector beams with different signal-to-noise ratios (SNRs).The SNR is defined as [61] where I signal and I total denote the intensity of the complex random mode E n (r) with and without introducing the white noise, respectively.The measured instantaneous intensities of E n (r) with SNR = −0.77dB, −6.53 dB, and −13.47 dB for the partially coher- ent vector beams with three-channel images "123" encoded inside are displayed in Fig. 6a, e, and i, respectively.The respective ratios of signal power to total power are 45.60 %, 18.20 %, and 4.30 %, respectively.Figure 6b-d, f-h, and g-l show the experimental results of the recovered three-channel images under the three different SNRs.It can be found that the three-channel images "123" can be reconstructed well in all three cases.The color "windmill" and its three primary-color-channels are recovered in the noisy environments with SNR = −0.17dB, −6.92 dB, and −13.68 dB, respectively.The corresponding signal power ratios to the total power are 49.03 %, 16.89 %, and 4.11 %, respectively.experimental recovered results are shown in Fig. 6m-x.The experimental results are consistent well with our predictions (see in simulation results in Supplementary Fig. S2) that the color images can also be recovered with high quality in the noisy environments.Our experimental and simulation results verified that our protocol is quite insensitive to the environmental noise.Thus, the two-point coherence Stokes vector of a partially coherent light beam can be regarded as a robust information carrier for three-dimensional optical encoding.

Three-channel optical encryption
Finally, our protocol can significantly enhance information security for optical information encoding, specifically enabling three-channel optical encryption by customizing the encoding rules of the response function H(r, v) .We would like to emphasize here that our protocol is adaptable to any encoding system, similar to those utilized in previous optical encryption protocols with fully coherent light.This flexibility arises from the fact that the response function H (r, v) can be designed arbitrarily.To illustrate the feasibil- ity of our three-channel optical encryption, we present an example in which the optical encryption system performs a fractional Fourier transform with a controllable fractional order.The corresponding transfer function is expressed as [13,28]: where A 0 = −icscϕ E / f and ϕ E = p E π/2 with p E being the fractional order of the frac- tional Fourier transform.The value of p E serves as an encryption key for the protocol, with p E ∈ (0, 1] ; setting p E = 1 reduces it to the ordinary Fourier transform.Following from Eq. ( 10), the ciphertext can be expressed as: (20) where S j (v) with j ∈ (1, 2, 3) denote the three-channel plaintext information and r = r 1 − r 2 .The three-channel plaintext information is now encoded into the coher- ence Stokes vector of the partially coherent beam the fractional Fourier transform system.In the experiment, such the partially coherent beam can still be synthesized with the RMR method.We only require to modulate the dynamic CGHs in the SLM by following Eqs.( 13) and ( 14).After the coherence Stokes vector being measured, the encoded three-channel information can be recovered by the inverse fractional Fourier transform, i.e., Above, S R j (v) denotes the recovered three-channel information, S R j (r 1 , r 2 ) denotes the measured coherence Stokes vector, and ϕ R = p R π/2 with p R serving as a decoding key.It is noted that only when the measured S R j (r 1 , r 2 ) matches the ciphertext S j (r 1 , r 2 ) , and the decoding key p R = p E , the encoded information can be correctly recovered.
As a proof-of-principle experiment, the encryption key, i.e., the fractional order of the transform system is set to be p E = 0.80 .Figure 7 depicts our experimental decryp- tion results for the three-channel images "123" and the color "windmill".Figure 7a-c and j-m show the recovered results with correct decryption key, i.e., p R = 0.80 , while Fig. 7d-f and n-q show the recovered results with incorrect decryption key, i.e., p R = 1.00 .It can be inferred from the recovered results that the target images can be well recovered only when the encryption and decoding keys match well.Otherwise, the recovered images cannot be recognized.It is emphasized that our decryption results are quite sensitive to the value of p R (see in Supplementary Fig. S3).We also carry out the experiment by introducing the white noise in the partially coherent (21) 2 ) d 2 r 1 d 2 r 2 .

Conclusions
In summary, we introduced a novel approach for three-channel optical encoding and encryption by engineering a single partially coherent vector light beam.The three-channel information carried by the polarization Stokes vector of a fully coherent light beam is initially encoded into the three components of the two-point coherence Stokes vector of a partially coherent beam via the generalized van Cittert-Zernike theorem.Subsequently, the encoded three-channel information is recovered by measuring the coherence Stokes vector through the recently introduced generalized Hanbury Brown-Twiss experiment.The advanced protocol exhibits several key advantages, including increased resilience to disturbances induced by complex environments and enhanced security for the encoding information through the introduction of additional quantities as encryption keys in the encoding system.In addition, we have shown that by associating the three independent components of the coherence Stokes vector with three primary color channels, the coherence Stokes vector can be flexibly utilized for color image encoding and encryption.The feasibility of our protocol has been demonstrated through a proof-of-principle experiment and numerical simulation.Both results affirmed that the encoded three-channel optical images and color images can be successfully decoded or decrypted even in the presence of high levels of random noise during the recovery process.Our findings suggest that partially coherent light engineering may find applications in multi-target image encryption and high-security, multi-channel optical communication in complex environments.

Fig. 1
Fig. 1 Principle of three-channel optical information encoding and encryption based on engineering the coherence Stokes vector of a partially coherent beam.a Schematic of encoding three-channel information S 1 (v) , S 2 (v) , and S 3 (v) into the coherence Stokes vector S(r 1 , r 2 ) = [S 1 (r 1 , r 2 ), S 2 (r 1 , r 2 ), S 3 (r 1 , r 2 )] of a partially coherent beam via the encoding system H.The three-channel information is decoded by the decoding system H −1 .b Three-channel images encoding and decoding.c Color image encoding and decoding

Fig. 2 a
Fig. 2 a Experimental setup for realizing the three-channel optical information encoding and decoding.b The target three-channel optical images.c A realization of the corresponding computer-generated holograms (CGH) loaded in the spatial light modulator (SLM) for generating the components of the field realization.d and e The experimentally measured instantaneous intensities for the x and y field components.f The experimental results of the coherence Stokes vector by spatial average over the instantaneous intensities.g The corresponding simulation results

Fig. 3
Fig. 3 Experimental results of the recovered three-channel images from the measured coherence Stokes vectors.a, h, and l are the target images.b-g The recovered images "123" and their enlarged images.i-k The recovered images "ABC".m-o The recovered images of "circle, square, and triangle" shapes.The corresponding simulation results for the recovered three-channel images are shown in the inset (within the yellow dashed circle)

Fig. 4
Fig. 4 Experimental results of the recovered color images from the measured coherence Stokes f, k, and p are the target images.b-e The recovered colorful letters "RGB" and its the separate primary-color-channels.g-j The recovered color windmill.l-o The recovered color peking opera mask.q-t The recovered lab logo.The corresponding simulation results for the recovered color images are shown in the inset (within the yellow dashed circle)

Fig. 5
Fig. 5 Experimental results of the recovered color images with variable number N of intensities realizations.a-j The recovered color windmill.k-t The recovered color peking opera mask

Fig. 6
Fig. 6 Robustness of the recovered images in noisy environments.a, e, and i The measured instantaneous intensities of E n (r) with SNR = −0.77dB, −6.53 dB, and −13.47 dB for the partially coherent vector beams with three-channel images "123" encoded inside.b-d,f-h, and g-lThe corresponding recovered three-channel images in the noisy environments.m-p, q-t, and u-x The recovered color images in the noise environments with SNR = −0.17dB, −6.92 dB, and −13.68 dB, respective

Fig. 7
Fig. 7 Three-channel optical encryption with coherence Stokes vector.a-c and j-m show the recovered three-channel images and color image with correct decryption key, i.e., p R = 0.80 .d-f and n-q show the recovered three-channel images and color image with incorrect decryption key, p R = 1.00 .g-i and r-u show the recovered images in the noise environments with SNR = −13.40dB and SNR = −13.61dB, respectively