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

Phase control scheme of the coherent beam combining system for generating perfect vectorial vortex beams assisted by a Dammann vortex grating

Open Access Open Access

Abstract

Based on Dammann vortex grating and adaptive gain stochastic parallel gradient descent algorithm, we theoretically proposed a phase control technology scheme of the coherent beam combining system for generating perfect vectorial vortex beams (VVBs). The simulated results demonstrate that the discrete phase locking for different types of VVBs (including vortex beams, vector beams, and generalized VVBs) can be successfully realized. The intensity distributions, polarization orientation, Pancharatnam phases, and beam widths of different |Hm,n〉 states with the obtained discrete phase distribution further prove that the generated beams are perfect VVBs. Subsequently, the phase aberration residual for different VVBs is evaluated using the normalized phase cosine distance function, and their values range from 0.01 to 0.08, which indicates the obtained discrete phase distribution is close to the ideal phase distribution. In addition, benefitting from the high bandwidth of involved devices in the proposed scheme, the influence of dynamic phase noise can be negligible. The proposed method could be beneficial to realize and switch flexible perfect VVBs in further applications.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Vectorial vortex beams (VVBs) are a special type of structured beams with both helical phase and spatially non-uniform polarization distributions [13], which possess simultaneously the spin angular momentum and orbit angular momentum (OAM), and support more encode channels. Perfect VVBs independent of the polarization orders and topological Pancharatnam charges are perfect candidates for spin-orbit interaction [4], free space communication [5], quantum cryptography [6], and optical storage [7]. Therefore, various methods for generating perfect VVBs by employing the spatial light modulators (SLMs) or metasurfaces have been proposed successively [814]. Unfortunately, lower damage thresholds and modulation speed of these devices have difficulty in achieving high power and switchable perfect VVBs.

As a widely studied technology, coherent beam combining (CBC) provides a promising solution for greatly scaling up the output power and maintaining good beam quality [1517]. More importantly, taking the advantage of flexible and independent regulation of amplitude, phase, and polarization of each beamlet in the laser array, the theoretical or experimental studies for vortex beams [1822], vector beams [23,24], and Bessel-Gaussian beam [25], have been carried out. In our previous theoretical work, we proposed an approach to generate perfect VVBs based on CBC with a specially designed radial phase-locked Gaussian laser array [26].

To achieve the desired structured beams based on the CBC, the prerequisite is that the discrete phases of laser array should be arranged and locked in a specific distribution. In addressing this challenge, several approaches have been proposed [2729]. In 2009, Hou et al. proposed that phase information can be extracted for the intensity distributions at the non-focal plane, avoiding the problem of non-uniqueness of phase information at the focal plane. With an optimization algorithm, the desired vortex beams can be generated based on the CBC [27]. Based on the heterodyne method, the phase of each beamlet can be measured and modulated, and the generation of vortex beams with topological charges ranging from 0 to 5 has been realized experimentally [28]. Recently, utilizing the unique advantages of deep learning (DL) in big data and image processing, the phase control method based on DL has been extensively investigated. In 2021, Tünnermann proposed a deep reinforcement learning model to achieve the discrete vortex phase locking in a CBC system by choosing an appropriate reward function [29]. Hou et al. developed a DL-assisted, two-stage phase control method for generating the mode-programmable vortex beams. First, a DL network is trained by a large number of intensity images on the non-focal plane or diffracted by a designed Dammann vortex grating (DVG) [30,31], the phase errors are preliminarily compensated by the trained DL network. Subsequently, the phase errors are further eliminated by an optimization algorithm, thereby achieving precise phase control of the vortex beams. However, these reported methods mostly are focused on the generation of vortex beams, the phase control methods for other types of structured beams however have received less attention.

In this work, based on the DVG and adaptive gain stochastic parallel gradient descent (SPGD) algorithm, a phase control technology scheme for generating perfect VVBs in the CBC system is proposed. For different types of VVBs (including vortex beams, vector beams, and generalized VVBs), the corresponding discrete phase distributions are optimized. With the obtained discrete phase distributions, we verify whether the generated beams are perfect VVBs or not by analyzing different parameters, such as intensity distributions, polarization orientation, Pancharatnam phases, and beam widths. In addition, the phase aberration residuals for different perfect VVBs are evaluated. Finally, based on the proposed method and devices involved, the influence of dynamic phase noise on the optimization results is discussed preliminarily. The developed phase control scheme provides a potential solution for generating high power perfect VVBs in the CBC system, making the application of high power perfect VVBs possible.

2. Principle and method

DVG is a special type of Dammann grating that combines Dammann grating with a helical phase structure, which can be used to generate vortex arrays or detect vortex beams. Its mathematical expression can be written as [32,33]:

$$G = \sum\limits_{n = 1}^{n = m} {\textrm{exp} (\textrm{i}\frac{{2\mathrm{\pi }({a_n}x + {b_n}y)}}{T} + \textrm{i}{l_n}\varphi )}$$
where ln represents the topological charge of the nth vortex beam, T is the period of the grating, and an and bn are the diffraction orders of the nth vortex beam in the x and y directions. In previous research, the number of beamlets is set as M = 26. According to Ref. [31], the number of coherently synthesized vortex beam m = 2[M/3] + 1 ≈ 17, and the topological charge range is from −8 to 8. Therefore, a DVG with a 5 × 5 array is generated, as shown in Fig. 1(a). Figure 1(b) depicts the phase distribution at the central region of DVG. The designed DVG’s corresponding far field diffraction distribution is shown in Fig. 1(c). It can be seen that different locations correspond to vortex beams with different topological charges, where the center position is the 0-order vortex beam, which is exactly a solid point. Figure 1(d) marks the corresponding positions of vortex beams with specific topological charges.

 figure: Fig. 1.

Fig. 1. (a) The phase distribution, (b) phase distribution at the central region, (c) far field diffraction distribution, and (d) corresponding positions of vortex beams with different topological charges of designed DVG of a 5 × 5 array

Download Full Size | PDF

Figure 2 depicts the diffraction distributions of single vortex beams or composite vortex beams with different topological charges passing through a DVG, and the insets show the shape of beam spots at the corresponding diffraction order. Obviously, it can be seen that for a single vortex beam, different topological charges correspond to different positions. While for a composite vortex beam, it is considered that two vortex beams with different topological charges after interacting with a DVG, respectively, producing two Gaussian shaped beam spots at their corresponding positions [34].

 figure: Fig. 2.

Fig. 2. The diffraction distributions of single vortex or composite vortex beams with different topological charges passing through a DVG. (a) l = −2, (b) l = 2, (c) l1= −2, l2= 2, (d) l1= −1, l2= 2. The insets depict the shape of beam spots at the corresponding diffraction order.

Download Full Size | PDF

In previous research [26], by employing a specially designed radial phase-locked Gaussian laser array composed of two discrete vortex arrays with right-handed (RH) and left-handed (LH) circularly polarized states in turn, perfect VVBs based on CBC can be obtained. The principle underlying this research is that VVBs can be considered the superposition of RH and LH circularly polarized vortices, written as a Jones matrix in the form of [35]

$$|{{\psi_{m,n}}} \rangle = \frac{{\psi _R^m\textrm{exp} ({\textrm{i}m\varphi } )}}{{\sqrt 2 }}\left[ {\begin{array}{{c}} 1\\ \textrm{i} \end{array}} \right] + \frac{{\psi _L^n\textrm{exp} ({\textrm{i}n\varphi } )}}{{\sqrt 2 }}\left[ {\begin{array}{{c}} 1\\ { - \textrm{i}} \end{array}} \right]$$
where $\psi _R^m$ and $\psi _L^n$ are RH and LH circularly polarized vortices’ amplitude coefficients, φ is the azimuthal angle of the polar coordinate system. When VVBs pass through a quarter wave plate with a fast axis angle of - 45 degrees, the beam distributions can be represented as
$${E_{m,n}} = \left[ {\begin{array}{cc} \textrm{1}&\textrm{i}\\ \textrm{i}&\textrm{1} \end{array}} \right]\left| {{\psi _{m,n}}} \right\rangle = \sqrt 2 \psi _R^m\textrm{exp} \left( {\textrm{i}m\varphi + {{\textrm{i}\pi } / 2}} \right)\left[ {\begin{array}{c} 0\\ 1 \end{array}} \right] + \sqrt 2 \psi _L^n\textrm{exp} \left( {\textrm{i}n\varphi } \right)\left[ {\begin{array}{c} 1\\ 0 \end{array}} \right]$$

The RH and LH circularly polarized vortex beams are converted into linearly polarized vortex beams with polarization directions of x and y, respectively.

Based on the above discussion and in combination with the characteristics of DVG, we proposed a phase control method of the CBC system for generating perfect VVBs, as schematically shown in Fig. 3. The laser is split into multiple channels by a splitter. In each channel, the power is amplified by a main amplifier, and the output phase is modulated by an electro-optic phase modulator. The output laser is emitted to free space from a special collimation array, which is composed of two Gaussian laser arrays with RH and LH circularly polarized states in turn. After passing through the spectroscope, the most of emitted laser array is reflected into a dynamometer, and a small portion of the transmitted beam, below the damage threshold of the SLM or DVG mask, will be used for phase control. Then, a quarter wave plate is utilized to convert RH and LH circularly polarized beams into linearly polarized beams in the x and y polarization directions, subsequently which are separated after passing through a polarization beam splitter. By using a half wave plate, the polarization direction of the reflected beam is adjusted to coincide with the polarization direction of the transmitted beam. Then two beams can coaxially transmit after passing through a beam splitter. A 4f imaging system with an adjustable diameter circular aperture pinhole is used to filter unwanted sidelobes. The two beams are then incident on an SLM or a DVG mask loaded with the phase of the designed DVG, which is on the front focal plane of the lens. Finally, a photodetector (PD) array in the back focal plane of the lens is used to receive the diffraction beam intensity. The size of the PD array is designed to correspond to the diffraction order of the DVG. In addition, a pinhole is placed before each PD to receive only the beam intensity in the central region. The detailed description of the simulation is given in Supplement 1, S1.

 figure: Fig. 3.

Fig. 3. Schematic illustration of phase control of CBC system to tailor perfect VVBs with a designed DVG. QWP: quarter wave plate, M: mirror, HWP: half wave plate, PBS: polarization beam splitter, BS: beam splitter, SLM: spatial light modulator.

Download Full Size | PDF

According to Eq. (2), vortex beams and vector beams are two special types of VVBs. The vortex beam is composed of RH and LH circularly polarized vortex beams with the same topological charges, the vector beam is composed of RH and LH circularly polarized vortex beams with the opposite topological charges, and the generalized VVBs are decomposed into RH and LH circularly polarized vortex beams with the different topological charges. Therefore, using our proposed approach and based on Fig. 2, after diffraction of the designed DVG, the vortex beams will form a Gaussian shaped spot at the corresponding position, while the vector beams will generate double positional symmetric Gaussian shaped spots, and the VVBs will produce two Gauss shaped spots at the asymmetric position. Therefore, on the one hand, classification of different types of structured light fields can be achieved based on the proposed method. On the other hand, based on the intensity signals detected by the PD array, the phase controller continuously and independently adjusts the phase modulator in each laser module using an optimization algorithm, thereby regulating the output phase of each laser module. Through continuous feedback iterations, phase locking of the desired VVBs can be achieved.

3. Results and discussions

We choose the adaptive gain SPGD algorithm as the optimization algorithm [36], and the cost function is defined as the total intensity in the target area

$$f = \int_0^{2\mathrm{\pi }} {\int_0^{{R_f}} {{I_{\textrm{mea }}}} } r\; \textrm{d}r\textrm{d}\varphi$$
where Rf is the radius of the pinhole before each photodetector, and Imea represents the intensity distribution. Through continuous iteration by employing the adaptive gain SPGD algorithm, the cost function f gradually evolves to its maximum.

In the numerical simulation, the wavelength, waist width, and aperture diameter of each laser array element are set to be λ = 1064 nm, ω0 = ρ0= 0.35 mm, the radius of collimation array, the number of beamlets, and the focal length of lens are set as R1 = 3.5 mm, M = 48, and f = 500 mm. To obtain sufficient statistics for evaluating optimization effects, 20 realizations for each type of VVBs have been independently simulated. Meanwhile, all initial discrete phase and the phase perturbation amplitude are set as 0 and 0.2, respectively. We choose the maximum detected intensity by PD and use it to normalize the cost function.

For the convenience of expression, the unified representation of the VVBs is |Hm,n〉 state (where m and n are non-zero integers). Clearly, the case of m = n represents the vortex beams with the topological charge of m, m = -n accounts for the vector beams with polarization order of |m|, and |m| ≠ | n | stands for the generalized VVBs.

To quantitatively characterize the changing trend of the purity of vortex mode during the iteration process, we adopt the concept of OAM spectrum, which is formed by the normalized energy weight of OAM mode Pl. Normalized energy weight Pl can be represented as [37]

$$\begin{aligned} {P_l} &= \frac{{{p_l}}}{{\sum\limits_{l ={-} \infty }^\infty {{p_l}} }}\\ {p_l} &= \int_0^{{R_f}} {\left\langle {{{|{{a_l}(\rho ,z)} |}^2}} \right\rangle \rho \textrm{d}\rho } \end{aligned}$$

Firstly, taking |H2,2〉 (namely, the vortex beam with the topological charge of 2) as the optimized target, the distribution of discrete phases using the adaptive gain SPGD algorithm is optimized, and the results are shown in Fig. 4. Due to the randomness of phase perturbations, each iteration is different, but its overall change trend is consistent. That is, at the beginning of each iteration, the cost function changes slowly; when the number of iterations reaches a certain number, the cost function increases rapidly; when the number of iterations exceeds 500, the cost function hardly changes significantly. Figure 4(b)-(e) and (f)-(i) show the diffraction distributions and normalized OAM spectrum for iterations of 1, 150, 300, and 600, respectively. By comparison, it is found that at the beginning of the iteration, due to the phase distribution being almost zero, the Gaussian spot is located at the 0-order diffraction position, which results in the OAM spectrum being almost concentrated at the 0-order diffraction; When the number of iterations increases to 150, part of the energy has been transferred to the 2-order diffraction, and the corresponding OAM spectral intensity is 0.4611; When the iterations number increases to 300, Gaussian spot also appears at the 2-order diffraction position, and the corresponding OAM spectral intensity increases to 0.8529; When the iterations number is 600, the diffraction distribution is basically consistent with that in Fig. 2(a), and the OAM spectral intensity increases to 0.9898, which means that the discrete phase locking for vortex beam with the topological charge of 2 has been achieved through algorithm iteration.

 figure: Fig. 4.

Fig. 4. Optimized results for desired |H2,2〉 state: (a) The convergence curves of the cost function with the number of iterations, (b)–(e) and (f)–(i) are the diffraction distributions and normalized OAM spectrum with iterations of 1, 150, 300 and 600, respectively. The insets show the shape of beam spots at the corresponding diffraction order.

Download Full Size | PDF

Similarly, we choose the |H−2,−2〉 (namely, the vortex beam with the topological charge of −2) state as the optimization target, the optimized results are shown in Fig. 5. The convergence curves change trend of the cost function are consistent with that of |H2,2〉 state. With the increase of iterations number, the intensity of the -2-order diffraction position gradually increases from zero to a Gaussian spot, and the OAM spectral intensity also increases from 0 to 0.9914, which has also been proven that the discrete phase locking for vortex beam with the topological charge of −2 has been successfully achieved.

 figure: Fig. 5.

Fig. 5. Optimized results for desired |H−2,−2〉 state: (a) The convergence curves of the cost function with the number of iterations, (b)–(e) and (f)–(i) are the diffraction distributions and normalized OAM spectrum with iterations of 1, 150, 300 and 600, respectively. The insets show the shape of beam spots at the corresponding diffraction order.

Download Full Size | PDF

According to the principle of beam splitting and combing depicted in Fig. 3, if the incident beam is a vector beam, two Gaussian spots will be generated at the symmetrical positions about 0-order diffraction (as shown in Fig. 2(c)). Here, taking |H2,−2〉 (corresponding to the vector beam with polarization order of 2) state as the optimization target, the cost function is defined as the summation of the intensity at the -2-order and 2-order diffraction. Figure 6 gives the optimized results for |H2,−2〉 state. Clearly, the convergence speed is similar to that of a vortex beam. At the beginning of the iteration, there is almost no intensity in the central region of the - 2 and 2-order diffraction, and the corresponding OAM spectral intensity is also about 0; When the iterations number increases to 150, a small amount of energy appears in the central region of the - 2 and 2-order diffraction, but the distribution is relatively dispersed, with corresponding OAM spectral intensities of 0.0331 and 0.0351, respectively; When the iterations number is 300, the energy at the −2 and 2-order diffraction gradually concentrates towards the central region, and the OAM spectral intensity also increases to 0.3520 and 0.3894, respectively; When the iterations number reaches 600, the energy at the −2 and 2-order diffraction points is almost concentrated in the central region, and the OAM spectral intensities also reache 0.4996 and 0.4868, which is very close to the ideal case of 0.5 and 0.5. This indicates we have also successfully realized phase locking for the desired vector beams with the proposed method.

 figure: Fig. 6.

Fig. 6. Optimized results for desired |H2,−2〉 state: (a) The convergence curves of the cost function with the number of iterations, (b)–(e) and (f)–(i) are the diffraction distributions and normalized OAM spectrum with iterations of 1, 150, 300 and 600, respectively. The insets show the shape of beam spots at the corresponding diffraction order.

Download Full Size | PDF

Similarly, for generalized VVBs, after DVG diffraction, there are also two Gaussian spots, which only appear at the two asymmetric positions with respect to 0-order diffraction, as shown in Fig. 2(d). Figure 7 displays the optimized results for the |H−1,2〉 state, and the cost function is defined as the summation of the intensity at the -1-order and 2-order diffraction. The convergence curves of the cost function, diffraction distribution, and normalized OAM spectrum are consistent with those of vector beams. Finally, when the iterations number is 600, the OAM spectral intensities of the −1 and 2-order are 0.5214 and 0.4714, respectively, which are close to ideal cases of 0.5 and 0.5. Therefore, using the proposed method, the discrete phase locking for desired generalized VVBs can also be achieved.

 figure: Fig. 7.

Fig. 7. Optimized results for desired |H−1,2〉 state: (a) The convergence curves of the cost function with the number of iterations, (b)–(e) and (f)–(i) are the diffraction distributions and normalized OAM spectrum with iterations of 1, 150, 300 and 600, respectively. The insets show the shape of beam spots at the corresponding diffraction order.

Download Full Size | PDF

In addition, the discrete phase locking for desired higher order generalized VVBs (such as |H−2,4〉 state) is also numerically investigated, as shown in Fig. 8. For higher order vortices, to obtain more phase information, the circular aperture pinhole of the 4f imaging system is appropriately expanded (ω0 is set as 0.7 mm), so that it can be seen that the intensity distribution area on each diffraction order is also appropriately expanded, but this has no impact on the convergence curve. After 600 iterations, the OAM spectral intensities of the −2 and 4-order are 0.4345 and 0.4379, respectively, which are slightly lower than the |H−1,2〉 state, but the convergence effect is still acceptable. Furthermore, in order to verify the universality of the proposed method, more examples of higher order VVBs are given in Supplement 1, S2.

 figure: Fig. 8.

Fig. 8. Optimized results for desired |H−2,4〉 state: (a) The convergence curves of the cost function with the number of iterations, (b)–(e) and (f)–(i) are the diffraction distributions and normalized OAM spectrum with iterations of 1, 150, 300 and 600, respectively. The insets show the shape of beam spots at the corresponding diffraction order.

Download Full Size | PDF

Based on the discrete phase information obtained from the optimization algorithm, the intensity distribution, polarization orientation, and Pancharatnam phases of |H2,2〉, |H−2,−2〉, |H−2,2〉, |H−1,2〉 and |H−2,4〉 states are numerically simulated [35,38], as shown in Fig. 9. According to the definition of VVBs, the polarization order and topological Pancharatnam charge are given as P = (m - n)/2 and l0 = (m + n)/2. For different |Hm,n〉 states, the polarization orientations separately rotate 0, 0, −4π, −3π, and −6π in a full circle, which correspond to the polarization order P of 0, 0, 1, −2, −1.5, and −3. While the Pancharatnam phases spirally vary 4π, −4π, 0, π, and 2π, which means that the topological Pancharatnam charges l0 are exactly equal to 2, −2, 0, 0.5, and 1. In addition, according to Eq. (12) of Ref. [39], the beam width ω of different |Hm,n〉 states are calculated and labeled in Fig. 9. Obviously, the beam width for different |Hm,n〉 states is basically the same, with the maximum variation not exceeding 1.1%, which indicates the beam width is independent of polarization order and topological Pancharatnam charge of VVBs, namely, the generated VVBs is perfect.

 figure: Fig. 9.

Fig. 9. Intensity distribution, polarization orientation, and Pancharatnam phase of different |Hm,n〉 states according to discrete phase distribution obtained from the proposed method. (a) m = n = 2; (b) m = n = −2; (c) -m = n = 2; (d) m = −1, n = 2; (e) m = −2, n = 4.

Download Full Size | PDF

However, compared to the Ref. [26], the intensity distribution, polarization distribution, and phase distribution of the generated VVBs are not so perfect, because there is still a certain deviation between the discrete phase distribution obtained from the proposed method and the ideal phase distribution. Here, we adopted the Normalized Phase Cosine Distance (NPCD) proposed by Zuo et al. to evaluate the phase aberration residual, and the NPCD is defined as [40]

$$NPCD = \frac{1}{2} - \frac{1}{{2(M - 1)}}\sum\limits_{m = 1}^{M - 1} {\cos } ({\varphi_{it\textrm{ }}^m - \varphi_{id}^m} )$$
where ${\varphi _{it}}$ and ${\varphi _{id}}$ are the discrete phase distribution obtained from the proposed method and ideal phase distribution, respectively. The NPCD value range is 0 to 1, and the closer the NPCD is to 0, the closer the obtained discrete phase distribution is to the ideal phase distribution. According to Eq. (6), the NPCD for different |Hm,n〉 states are calculated at initial iteration and 1000 iterations, as shown in Table. 1. Since the phase distribution of the initial iteration is set to 0, the NPCD for all states is 0.5. After 1000 iterations, the NPCD of |H2,2〉, |H−2,2〉, and |H−1,2〉 states decrease to about 0.01, indicating that the discrete phase distribution obtained from the proposed method is very close to the ideal phase distribution. while the NPCD of |H−2,4〉 state only drops to 0.0855, slightly higher than that of other states, which proves that there is a slight deviation between the obtained discrete phase distribution and the ideal phase distribution for higher order generalized VVBs. This is also consistent with the results characterized by parameters such as OAM spectral intensity and intensity distribution. The NPCD can be further reduced by selecting appropriate parameters, such as phase perturbation amplitude. However, a smaller disturbance intensity will result in the longer convergence time, so it is necessary to tradeoff of convergence speed and convergence effect. The detailed description of the influence of phase perturbation amplitude is given in Supplement 1, S3.

Tables Icon

Table 1. NPCD for different |Hm,n〉 states at initial iteration and 1000 iterations

In practical application, the changes of refractive index induced by thermal effects in the optical fiber amplifiers, mechanical vibration of the system, and the surrounding environment both lead to the dynamical phase noise, which result in affecting the convergence speed and effectiveness. Benefiting from the high rate of electro-optic phase modulators, phase controller and high bandwidth of PD, the frequency of a single iteration approximates 100 megahertz [31,41]. For the proposed optical system, it takes 500 iterations to achieve a good convergence effect, so the bandwidth of the SPGD algorithm is about 200 kHz. According to the current experimental measurement results of optical fiber amplifier systems, the frequency of dynamical phase noise of fiber amplifiers is approximately between 10 Hz and 10 kHz [4244]. When the bandwidth of the SPGD algorithm is far greater than the frequency of the dynamic phase noise, that is, the optimization of the SPGD algorithm has been completed, and the dynamic phase noise has not yet changed. This means that the dynamic phase noise is only equivalent to an initial phase noise and does not affect the optimization results of the SPGD algorithm.

4. Conclusions

In conclusion, based on the DVG and adaptive gain SPGD algorithm, the phase control technology scheme of the CBC system for generating perfect VVBs is proposed. By analyzing the convergence curves, diffraction distributions, and normalized OAM spectrum, the discrete phase locking of different types of VVBs (including vortex beams, vector beams, and generalized VVBs) is successfully realized. Based on obtained discrete phase distribution, the intensity distribution, polarization orientation, Pancharatnam phases, and beam widths of different |Hm,n〉 states are numerically simulated, and the results further confirm that the generated beams are perfect VVBs. In addition, the difference between the obtained discrete phase distribution and ideal phase distribution for different |Hm,n〉 states is evaluated using the NPCD. The NPCDs will decrease from 0.5 to between 0.01 and 0.08, which means the obtained discrete phase distribution is close to the ideal phase distribution. Finally, the influence of dynamic phase noise on the SPGD algorithm is analyzed. The results show that due to the high bandwidth of the devices used in the proposed phase control scheme, the bandwidth of the SPGD algorithm can approximate 200 kHz, so the impact of dynamic phase noise is negligible to the CBC system. Our work enriches the phase control scheme for generating structured beams in the CBC system and enables further applications of a high power perfect VVBs.

Funding

Youth Innovation Promotion Association XIOPM-CAS (XIOPMQCH2021003); National Natural Science Foundation of China (62005310, 62171443); Key Research and Development Projects of Shaanxi Province (2021GY-298, 2023-ZDLGY-37).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Supplemental document

See Supplement 1 for supporting content.

References

1. N. Radwell, R. D. Hawley, J. B. Götte, and S. Franke-Arnold, “Achromatic vector vortex beams from a glass cone,” Nat. Commun. 7(1), 10564 (2016). [CrossRef]  

2. Z. Liu, Y. Liu, Y. Ke, Y. Liu, W. Shu, H. Luo, and S. Wen, “Generation of arbitrary vector vortex beams on hybrid-order poincaré sphere,” Photonics Res. 5(1), 15–21 (2017). [CrossRef]  

3. Z. Bomzon, V. Kleiner, and E. Hasman, “Pancharatnam–berry phase in space-variant polarization-state manipulations with subwavelength gratings,” Opt. Lett. 26(18), 1424–1426 (2001). [CrossRef]  

4. M. Li, Y. Cai, S. Yan, Y. Liang, P. Zhang, and B. Yao, “Orbit-induced localized spin angular momentum in strong focusing of optical vectorial vortex beams,” Phys. Rev. A 97(5), 053842 (2018). [CrossRef]  

5. G. Milione, T. A. Nguyen, J. Leach, D. A. Nolan, and R. R. Alfano, “Using the non-separability of vector beams to encode information for optical communication,” Opt. Lett. 40(21), 4887–4890 (2015). [CrossRef]  

6. X. Li, T. H. Lan, C. H. Tien, and M. Gu, “Three-dimensional orientation-unlimited polarization encryption by a single optically configured vectorial beam,” Nat. Commun. 3(1), 998 (2012). [CrossRef]  

7. L. Zeng, Y.H. Ye, M. Dong, W.H. Zhang, E. Z. Li, D.C. Li, D.S. Ding, and B.S. Shi, “Optical memory for arbitrary perfect poincaré states in an atomic ensemble,” Opt. Lett. 48(2), 477–480 (2023). [CrossRef]  

8. P. Li, Y. Zhang, S. Liu, C. Ma, L. Han, H. Cheng, and J. Zhao, “Generation of perfect vectorial vortex beams,” Opt. Lett. 41(10), 2205–2208 (2016). [CrossRef]  

9. D. Li, S. Feng, S. Nie, C. Chang, J. Ma, and C. Yuan, “Generation of arbitrary perfect poincaré beams,” J. Appl. Phys. 125(7), 073105 (2019). [CrossRef]  

10. Y. Bao, J. Ni, and C.-W. Qiu, “A minimalist single-layer metasurface for arbitrary and full control of vector vortex beams,” Adv. Mater. 32(6), 1905659 (2020). [CrossRef]  

11. M. Vergara and C. Iemmi, “Multiple quasi-perfect vector vortex beams with arbitrary 3D position on focus,” Appl. Opt. 61(20), 5926–5933 (2022). [CrossRef]  

12. T. Wang, S. Fu, F. He, and C. Gao, “Generation of perfect polarization vortices using combined gratings in a single spatial light modulator,” Appl. Opt. 56(27), 7567–7571 (2017). [CrossRef]  

13. S. Fu, T. Wang, and C. Gao, “Generating perfect polarization vortices through encoding liquid-crystal display devices,” Appl. Opt. 55(23), 6501–6505 (2016). [CrossRef]  

14. M. Liu, P. Huo, W. Zhu, C. Zhang, S. Zhang, M. Song, S. Zhang, Q. Zhou, L. Chen, H. J. Lezec, A. Agrawal, Y. Lu, and T. Xu, “Broadband generation of perfect poincare beams via dielectric spin-multiplexed metasurface,” Nat. Commun. 12(1), 2230 (2021). [CrossRef]  

15. G. D. Goodno, H. Komine, S. J. McNaught, S. B. Weiss, S. Redmond, W. Long, R. Simpson, E. C. Cheung, D. Howland, P. Epp, M. Weber, M. McClellan, J. Sollee, and H. Injeyan, “Coherent combination of high-power, zigzag slab lasers,” Opt. Lett. 31(9), 1247–1249 (2006). [CrossRef]  

16. T. Fan, “Laser beam combining for high-power, high-radiance sources,” IEEE J. Sel. Top. Quantum Electron. 11(3), 567–577 (2005). [CrossRef]  

17. Z. Liu, X. Jin, R. Su, P. Ma, and P. Zhou, “Development status of high power fiber lasers and their coherent beam combination,” Sci. China Inf. Sci. 62(4), 41301 (2019). [CrossRef]  

18. D. Zhi, T. Hou, P. Ma, Y. Ma, P. Zhou, R. Tao, X. Wang, and L. Si, “Comprehensive investigation on producing high-power orbital angular momentum beams by coherent combining technology,” High Power Laser Sci. Eng. 7, e33 (2019). [CrossRef]  

19. J. Long, T. Hou, Q. Chang, T. Yu, R. Su, P. Ma, Y. Ma, P. Zhou, and L. Si, “Generation of optical vortex lattices by a coherent beam combining system,” Opt. Lett. 46(15), 3665–3668 (2021). [CrossRef]  

20. M. Veinhard, S. Bellanger, L. Daniault, I. Fsaifes, J. Bourderionnet, C. Larat, E. Lallier, A. Brignon, and J.C. Chanteloup, “Orbital angular momentum beams generation from 61 channels coherent beam combining femtosecond digital laser,” Opt. Lett. 46(1), 25–28 (2021). [CrossRef]  

21. T. Hou, Q. Chang, J. Long, P. Ma, and P. Zhou, “Design considerations and performance analysis of a fiber laser array system for structuring orbital angular momentum beams: a simulation study,” Opt. Express 30(9), 15279–15309 (2022). [CrossRef]  

22. T. Yu, H. Xia, W. Xie, and Y. Peng, “Orbital angular momentum mode detection of the combined vortex beam generated by coherent combining technology,” Opt. Express 28(24), 35795–35806 (2020). [CrossRef]  

23. R. S. Kurti, K. Halterman, R. K. Shori, and M. J. Wardlaw, “Discrete cylindrical vector beam generation from an array of optical fibers,” Opt. Express 17(16), 13982–13988 (2009). [CrossRef]  

24. P. Ma, P. Zhou, Y. Ma, X. Wang, R. Su, and Z. Liu, “Generation of azimuthally and radially polarized beams by coherent polarization beam combination,” Opt. Lett. 37(13), 2658–2660 (2012). [CrossRef]  

25. T. Yu, H. Xia, W. Xie, G. Xiao, and H. Li, “The generation and verification of bessel-gaussian beam based on coherent beam combining,” Results Phys. 16, 102872 (2020). [CrossRef]  

26. P. Ju, W. Fan, W. Gao, Z. Li, Q. Gao, G. Li, X. Jiang, and T. Zhang, “Generation of perfect vectorial vortex beams by employing coherent beam combining,” Opt. Express 31(7), 11885–11898 (2023). [CrossRef]  

27. T. Hou, Y. Zhang, Q. Chang, P. Ma, R. Su, J. Wu, Y. Ma, and P. Zhou, “High-power vortex beam generation enabled by a phased beam array fed at the nonfocal-plane,” Opt. Express 27(4), 4046–4059 (2019). [CrossRef]  

28. D. R. Gozzard, J. T. Spollard, P. G. Sibley, L. E. Roberts, and D. A. Shaddock, “Optical vortex beams with controllable orbital angular momentum using an optical phased array,” OSA Continuum 3(12), 3399–3406 (2020). [CrossRef]  

29. H. Tünnermann and A. Shirakawa, “Deep reinforcement learning for tiled aperture beam combining in a simulated environment,” J. Physics: Photonics 3(1), 015004 (2021). [CrossRef]  

30. T. Hou, Y. An, Q. Chang, P. Ma, J. Li, L. Huang, D. Zhi, J. Wu, R. Su, Y. Ma, and P. Zhou, “Deep-learning-assisted, two-stage phase control method for high-power mode-programmable orbital angular momentum beam generation,” Photonics Res. 8(5), 715–722 (2020). [CrossRef]  

31. T. Hou, Y. An, Q. Chang, J. Long, P. Ma, J. Li, and P. Zhou, “Deep learning of coherent laser arrays in angular domain for orbital angular momentum beams customization,” IEEE J. Sel. Top. Quantum Electron. 28(4), 1–10 (2022). [CrossRef]  

32. I. Moreno, J. A. Davis, D. M. Cottrell, N. Zhang, and X.-C. Yuan, “Encoding generalized phase functions on dammann gratings,” Opt. Lett. 35(10), 1536–1538 (2010). [CrossRef]  

33. N. Zhang, X. C. Yuan, and R. E. Burge, “Extending the detection range of optical vortices by dammann vortex gratings,” Opt. Lett. 35(20), 3495–3497 (2010). [CrossRef]  

34. S. Fu, S. Zhang, T. Wang, and C. Gao, “Measurement of orbital angular momentum spectra of multiplexing optical vortices,” Opt. Express 24(6), 6240–6248 (2016). [CrossRef]  

35. X. Yi, Y. Liu, X. Ling, X. Zhou, Y. Ke, H. Luo, S. Wen, and D. Fan, “Hybrid-order poincaré sphere,” Phys. Rev. A 91(2), 023801 (2015). [CrossRef]  

36. Q. Hu, L. Zhen, Y. Mao, S. Zhu, X. Zhou, and G. Zhou, “Adaptive stochastic parallel gradient descent approach for efficient fiber coupling,” Opt. Express 28(9), 13141–13154 (2020). [CrossRef]  

37. S. Fu and C. Gao, “Influences of atmospheric turbulence effects on the orbital angular momentum spectra of vortex beams,” Photonics Res. 4(5), B1–B4 (2016). [CrossRef]  

38. G. Milione, H. I. Sztul, D. A. Nolan, and R. R. Alfano, “Higher-order poincaré sphere, stokes parameters, and the angular momentum of light,” Phys. Rev. Lett. 107(5), 053601 (2011). [CrossRef]  

39. J. Pinnell, V. Rodríguez-Fajardo, and A. Forbes, “How perfect are perfect vortex beams?” Opt. Lett. 44(22), 5614–5617 (2019). [CrossRef]  

40. J. Zuo, H. Jia, C. Geng, Q. Bao, F. Zou, Z. Li, J. Jiang, F. Li, B. Li, and X. Li, “Deep learning piston aberration control of fiber laser phased array by spiral phase modulation,” J. Lightwave Technol. 40(12), 3980–3991 (2022). [CrossRef]  

41. T. Hou, Q. Chang, T. Yu, J. Long, H. Chang, P. Ma, R. Su, Y. Ma, and P. Zhou, “Switching the orbital angular momentum state of light with mode sorting assisted coherent laser array system,” Opt. Express 29(9), 13428–13440 (2021). [CrossRef]  

42. S. J. Augst, T. Y. Fan, and A. Sanchez, “Coherent beam combining and phase noise measurements of ytterbium fiber amplifiers,” Opt. Lett. 29(5), 474–476 (2004). [CrossRef]  

43. Q. Chang, T. Hou, J. Long, Y. Deng, H. Chang, P. Ma, R. Su, Y. Ma, and P. Zhou, “Experimental phase stabilization of a 397-channel laser beam array via image processing in dynamic noise environment,” J. Lightwave Technol. 40, 1–7 (2022). [CrossRef]  

44. H. Xiao, X. Wang, Y. Ma, B. He, P. Zhou, J. Zhou, and X. Xu, “Phase noise measurement of high-power fiber amplifiers,” Chin. Opt. Lett. 9(4), 041404 (2011). [CrossRef]  

Supplementary Material (1)

NameDescription
Supplement 1       Detailed description of the simulation, more examples of higer order VVBs, and the influence of phase perturbation amplitude

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

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 (9)

Fig. 1.
Fig. 1. (a) The phase distribution, (b) phase distribution at the central region, (c) far field diffraction distribution, and (d) corresponding positions of vortex beams with different topological charges of designed DVG of a 5 × 5 array
Fig. 2.
Fig. 2. The diffraction distributions of single vortex or composite vortex beams with different topological charges passing through a DVG. (a) l = −2, (b) l = 2, (c) l1= −2, l2= 2, (d) l1= −1, l2= 2. The insets depict the shape of beam spots at the corresponding diffraction order.
Fig. 3.
Fig. 3. Schematic illustration of phase control of CBC system to tailor perfect VVBs with a designed DVG. QWP: quarter wave plate, M: mirror, HWP: half wave plate, PBS: polarization beam splitter, BS: beam splitter, SLM: spatial light modulator.
Fig. 4.
Fig. 4. Optimized results for desired |H2,2〉 state: (a) The convergence curves of the cost function with the number of iterations, (b)–(e) and (f)–(i) are the diffraction distributions and normalized OAM spectrum with iterations of 1, 150, 300 and 600, respectively. The insets show the shape of beam spots at the corresponding diffraction order.
Fig. 5.
Fig. 5. Optimized results for desired |H−2,−2〉 state: (a) The convergence curves of the cost function with the number of iterations, (b)–(e) and (f)–(i) are the diffraction distributions and normalized OAM spectrum with iterations of 1, 150, 300 and 600, respectively. The insets show the shape of beam spots at the corresponding diffraction order.
Fig. 6.
Fig. 6. Optimized results for desired |H2,−2〉 state: (a) The convergence curves of the cost function with the number of iterations, (b)–(e) and (f)–(i) are the diffraction distributions and normalized OAM spectrum with iterations of 1, 150, 300 and 600, respectively. The insets show the shape of beam spots at the corresponding diffraction order.
Fig. 7.
Fig. 7. Optimized results for desired |H−1,2〉 state: (a) The convergence curves of the cost function with the number of iterations, (b)–(e) and (f)–(i) are the diffraction distributions and normalized OAM spectrum with iterations of 1, 150, 300 and 600, respectively. The insets show the shape of beam spots at the corresponding diffraction order.
Fig. 8.
Fig. 8. Optimized results for desired |H−2,4〉 state: (a) The convergence curves of the cost function with the number of iterations, (b)–(e) and (f)–(i) are the diffraction distributions and normalized OAM spectrum with iterations of 1, 150, 300 and 600, respectively. The insets show the shape of beam spots at the corresponding diffraction order.
Fig. 9.
Fig. 9. Intensity distribution, polarization orientation, and Pancharatnam phase of different |Hm,n〉 states according to discrete phase distribution obtained from the proposed method. (a) m = n = 2; (b) m = n = −2; (c) -m = n = 2; (d) m = −1, n = 2; (e) m = −2, n = 4.

Tables (1)

Tables Icon

Table 1. NPCD for different |Hm,n〉 states at initial iteration and 1000 iterations

Equations (6)

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

G = n = 1 n = m exp ( i 2 π ( a n x + b n y ) T + i l n φ )
| ψ m , n = ψ R m exp ( i m φ ) 2 [ 1 i ] + ψ L n exp ( i n φ ) 2 [ 1 i ]
E m , n = [ 1 i i 1 ] | ψ m , n = 2 ψ R m exp ( i m φ + i π / 2 ) [ 0 1 ] + 2 ψ L n exp ( i n φ ) [ 1 0 ]
f = 0 2 π 0 R f I mea  r d r d φ
P l = p l l = p l p l = 0 R f | a l ( ρ , z ) | 2 ρ d ρ
N P C D = 1 2 1 2 ( M 1 ) m = 1 M 1 cos ( φ i t   m φ i d m )
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.