ModelBased Tabu Search Algorithm for FreeSpace Optical Communication with a Novel Parallel Wavefront Correction System
 Author: Li Zhaokun, Zhao Xiaohui, Cao Jingtai, Liu Wei
 Publish: Journal of the Optical Society of Korea Volume 19, Issue1, p45~54, 25 Feb 2015

ABSTRACT
In this study, a novel parallel wavefront correction system architecture is proposed, and a modelbased tabu search (MBTS) algorithm is introduced for this new system to compensate wavefront aberration caused by atmospheric turbulence in a freespace optical (FSO) communication system. The algorithm flowchart is presented, and a simple hypothetical design for the parallel correction system with multiple adaptive optical (AO) subsystems is given. The simulated performance of MBTS for an AOFSO system is analyzed. The results indicate that the proposed algorithm offers better performance in wavefront aberration compensation, coupling efficiency, and convergence speed than a stochastic parallel gradient descent (SPGD) algorithm.

KEYWORD
Parallel correcting system architecture , Model based tabu search algorithm , Freespace optical communications , Atmospheric turbulence

I. INTRODUCTION
Recently freespace optical (FSO) communication has attracted more and more attention as it becomes widely used among the telecommunication community for both ground and spacebased wireless links and “lastmile” applications [1], due to its unregulated spectrum, high potential bandwidth, relatively low power requirement, low bit error rate, and ease of redeployment. However, phase disturbances from atmospheric turbulence along propagation paths, manifesting as intensity fluctuation (scintillation), beam wandering, and beam broadening at the receiver all lead to significant decrease of coupling efficiency [2], which seriously influences the stability and reliability of FSO communication systems [3].
An adaptive optical (AO) system is an effective method to improve laserbeam quality by correcting any wavefront aberration, and has already allowed great achievements [49]. Generally, in a conventional AO system a ShackHartmann wavefront sensor (SH sensor) [10] measures the optical phase deviations of an incoming wavefront, and a deformable mirror (DM) is used to compensate for the phase distortion. Then the DM generates a wavefront phase to compensate for the phase aberration based on phaseconjugation theory [11, 12].
Under strong scintillation, sensorless AO systems are proposed to compensate for the wavefront aberration. Some very effective blind control optimized algorithms, such as stochastic parallel gradient descent (SPGD) and simulated annealing (SA) have been proposed to improve the performance of the FSO system [13]. In this study, a new method called the modelbased tabu search (MBTS) algorithm is proposed to offer better performance in atmospheric turbulence compensation. Tabu search (TS) is a metaheuristic randomsearch algorithm [14] that starts from an initial solution and selects some specified directions for probing. To avoid being trapped in a local minimum, TS uses a flexible method to record the process of optimization and selection to guide the next search step, while a Zernike modelbased (MB) method in an AO system can decrease the dimension of the search space, relieve calculation complexity and accelerate convergence speed. Modelbased tabu search (MBTS), the combination of MB and TS approaches, will benefit from these two methods and perform better in an FSO system, since it can effectively improve coupling efficiency at the receiver and reduce energy loss in the laser propagation path. Most important is that MBTS uses very few iterations, far less than does SPGD (which usually needs hundreds of iterations). Moreover, to further improve the performance of the MBTS algorithm, a novel parallel correcting system is proposed. The newly invented system can increase the speed of correction of the wavefront aberration by simply combining several uniform AO subsystems together, each subsystem correcting certain orders of aberration by MBTS. The working principle of this new system is explained further.
This paper is organized as follows: Section 2 provides models of an FSO communication system, sensorless AO system, and DM. Section 3 gives an analysis of MBTS algorithm and its working principles in the FSO communication system, as well as a simple and novel hypothetical design for a parallel correcting system with multiple AO subsystems. In Section 4, some simulations are carried out to show the performance of the proposed MBTS in a sensorless AOFSO system, in comparison with SPGD algorithm. Finally, some conclusions from this study are given in Section 5.
II. METHODS
2.1. FSO Communication System with Model  Based Tabu Search Algorithm
A laser communication system in atmospheric turbulence works poorly, because an aberrated wavefront leads to energy loss in receiver. In an ideal case when the system does not introduce any aberration to the laser signal, all fluorescence is emitted from a point, so that the laser is well coupled to a single mode fiber (see Fig. 1 (a)). When atmospheric turbulence is present, the intensity in the focal plane depends on the wavefront of incident light, which spreads widely over the focal plane of lens due to the wavefront aberration; thus the coupling efficiency at the receiver will degrade, and the bit error rate will increase (Fig. 1 (b)) [15].
The simple schematic diagram of an FSO experimental setup is illustrated in Fig. 2. A laser from a distant source passes through an atmospheric disturbances box, in which the wavefront is disturbed by the atmospheric turbulence device. In a conventional AO system, when the laser enters it is reflected by a deformable mirror (DM) and directed to a wavefront sensor (WFS). The measured signal is sent to a controller to reconstruct the wavefront from the WFS measurement, and to compare the reconstructed wavefront to a reference plane wavefront for computation of the control signal for the DM. Consequently, the DM deforms its surface to counteract the wavefront aberration. In this way the laser wavefront is compensated, and thus the coupling efficiency in the communication receiver increases. However, in a recently proposed sensorless AO system, a CCD camera replaces the WFS. This camera measures the intensity at the focus plane to obtain a performance metric and generate a good control signal vector. In addition, a good algorithm should also provide a better solution in fewer iterations for a realtime system in practice.
Our goal is to find a control signal so that the DM can generate the best surface deformation to compensate for the wavefront aberration, so that in turn more laser energy is coupled into the single mode fiber.
2.2. Wavefront Compensation Device
In general, a DM with more actuators could generate a more precise correcting phase to compensate for lowerorder but highervalue aberration. However, a DM with too many actuators and too large a receiving aperture is not useful in laser transmission, because it requires reshaping of the laser beam to a large diameter to match the receiving aperture. In our analysis, a DM with 32 actuators is a suitable selection that generates an acceptable result and has no excessive requirement for beam reshaping. The normalized layout of the 32element DM actuators is shown in Fig. 3.
We approximate the DM influence function by a Gaussian mdel as follows:
where
ω is the coupling coefficient determined by the sizes of the electrode actuators and the DM, (x_{j} ,y_{j} ) are the coordinates of the center of thej th actuator,d is the interval between the adjacent actuators, andα is the Gaussian index. The phase compensationϕ (x ,y ) generated by the deformable mirror is given bywhere
v_{j} is thej th voltage of the actuators, and (x ,y ) is the location of each point on the DM. The numerical relationship between the aberration generated by the DM and the voltagesv_{j} applied to the actuators is linear.2.3. Modelbased Compensation in Adaptive Optics
From the conventional point of view, we often pay lot of attention to the control voltages applied to the actuators, but problems occur due to a solution space of high dimension, which can result in slow convergence, large computational burden, and unsatisfied convergence value of the control algorithm, since increasing the number of actuators brings a geometric increase in the number of calculations of the control solution. This is important for realtime performance of the system. A new, alternative approach is to use a Zernike model to fit induced wavefront aberration and to indirectly get appropriate voltages, instead of directly using a DM model. No matter how many actuators exist, we only need to obtain the linear relationship between the control voltages and Zernike coefficients, so that the number of calculations will decrease, which is better for an FSO system under highspeed workload.
We assume the wavefront aberration
ϕ and correcting aberrationφ can be represented by a series ofM andN Zernike functions, respectively. In most cases,M >N . The definitions ofϕ andφ areIn other words, the aberration of the incident wavefront
φ can be denoted by the vectorW , whose elements are the coefficientsw_{i} . Similarly, the correcting aberrationφ and the residual aberrationℜ can be represented by the vectorsU andW with elementsu_{i} andw_{i} respectively andℜ _{i} =w_{i} u_{i} [16]. To get the bestU through the minimization of var(ℜ )is our goal for the sensorless AO system.In practice, the correcting aberration
φ is excited by the actuator voltages from the DM, which is approximated by Eq. (2). The voltagev_{j} applied to thej th actuator determines the correcting aberrationφ . Ideally, the wavefront aberration fitting by the Zernike coefficients (Eq. (4)) and the phase aberration generated by the DM (Eq. (2)) should be identical, even if this cannot be completely realized. From Eq. (2) and (4) we obtainwhere
and
S_{i,k} is thek th effect of the voltage on thei th sampled value, andZ_{i,k} is thek th Zernike coefficient for thei th sampled value. Here we havem sampled values. IfU is obtained, the vectorA iswhere
S ^{†} is the pseudoinverse ofS . Obviously it can be calculated offline and is simplified; what we need to do is simply find the bestU = (u _{1}^{*},u _{2}^{*},…u_{N} ^{*}).2.4. Parallel Wavefront Correcting System
The sensorless AO system described above cannot take the speed of the proposed algorithm to the extreme, so here we put forward a new system architecture that combines many uniform AO subsystems, to increase the processing rate.
We outline a simple conception of this equipment with four AO systems based on the algorithm proposed in this paper, as examplified in Fig. 4. The incident light is divided into some number of parts by beam splitters. Each part of light passes through a sensorless AO subsystem and is corrected by the MBTS algorithm: AO subsystem I gets the best solution (
u _{1}^{*},u _{2}^{*}), AO subsystem II gets the best solution (u _{3}^{*},u _{4}^{*}), and so on. A master control system obtains the best voltage solutions based on the best Zernike mode fitting solutions fromu _{1}^{*} tou _{7}^{*} and Eq. (6) (mentioned in Section 3.1). In our opinion, the light should not be divided into too many parts, to avoid difficulty in measurement due to weak light intensity at the CCD camera. In our analysis, each AO subsystem deals with two orders, and the optimum number of parallel AO systems isceil (N/2), which is the smallest integer bigger than N/2 (where N is the order of aberration to be compensated); it is sufficient and easily selected, based on the specific circumstances.2.5. Tabu Search Algorithm
The tabu search (TS) algorithm is an effective way to find an optimal or nearly optimal global solution, which can guide the search process to avoid a local optimum and find the global optimum. The TS algorithm begins with an initial solution
S , finds a set of neighborhoodsN (S ), selects a new solutionS ^{*} fromN (S ) according to some rules, transforms fromS toS ^{*}, and repeats this process until the termination condition is met. To avoid meaningless cycling and becoming trapped in a local optimal solution, a tabu list of lengthL is introduced in this algorithm. The elements in the list are the solutions that have been selected fromN (S ). In other words, they are the latestL statuses to have been selected; in future processes, we cannot select these solutions for a certain number of iterations.In general, the TS algorithm selects a new solution according to a rule such as this: If the new
S ^{*} is in the tabu list, whether it should be selected or not is based on its corresponding fitness. If the fitness is better than the current best value, then it can be selected directly; otherwise it is disregarded, andS ^{*} is replaced by the secondbest solution withinN (S ). Hence we find whether it is in the tabu list or not, and this cycle repeats.The detailed process is described as follows:
(1) Select an initial solution. A good initial solution is very import to finding a global optimal solution. Here it is the Zernike coefficients that generate the correcting aberration
φ . The proper approach is to know the environmental conditions, record the phase aberration for as long as possible, and use the expectation value as the initial condition. In the simulation, we take zero as a suitable starting point, assuming that the processing rate is fast compared to the DM’s frame rate(2) Find an appropriate
N (S ). To reduce the number of calculations, we decompose the solution into different groups, i.e. we take coefficientsv _{1} andv _{2} as a group (v _{1},v _{2}) and further groups (v _{3},v _{4}), (v _{5},v _{6}),…, respectively, such that the search occurs in a twodimensional space (see Fig. 5). Here we chooseN (S ) as a group of eight points adjacent to^{s} (shown in Fig. 6). Each point inN (S ) is a candidate solution that may be selected. During the search, the fitness is the variance of the residual wavefront phase var((ℜ )).(3) Select a new solution
S ^{*} fromN (S ) by the abovementioned rules, then append it to the tabu list (delete the first record if the list has no place for the new solutionS ^{*}, hence the deleted record could be selected again in a future search), and replaceS withS ^{*}.(4) If the termination condition is satisfied, we have obtained the best solution; otherwise, go to step (2) and loop.
The flow chart is shown as Fig. 7.
2.6. Performance Metric of our FSO System
The goal of an adaptive optical system is to minimize residual phase aberrations after an incoming wave passes the deformable mirror. This corresponds to a maximization of the Strehl ratio (ST). In this paper, we use ST as the system performance metric
J [17],where
A (r ) is the (complex) optical field in the focal plane, andr _{0 }is the desired onaxis location of the center of the fiber end within this plane. In this study, a CCD camera is used to measure the optical power emerging from the hole; the measured quantity is thus proportional to the integral of the intensity within the pinhole with radius^{a} , and we have an experimental performance metric [1820]For direct optimization of the power coupled into a fiber, the system performance metric
J is proportional to the absolute coupling efficiency. In case of a singlemode fiber, the theoretical value is proportional to the overlap integral of the optical fieldA (r ) and the mode profileM _{0}(r ) in the focal plane; thus the performance metricJ is given byAssume that the wavefront phase aberration satisfies a Gaussian distribution, ST can be estimated by variance RMS^{2} as follow
With increasing coupling efficiency, more energy is coupled into the singlemode fiber. When RMS^{2} is close to 0, we can get a simpler formula:
In practice, the pixel size of the CCD camera approximately equals the fiber diameter, so ST is [18]
where
A (i ) is the gray value of thei th pixel, andP is the number of pixels.III. RESULTS
Based on the previous analysis, 12 Zernike modes and the corresponding correcting aberration are given in Fig. 8.
Define the fitting error as fitting capability for the wavefront aberration of DM [18]
where RMS_{C} represents the RMS value of the wavefront phase aberration before correction and RMS_{R} the value of the residual wavefront phase aberration after correction. The smaller E_{fe} is, the more precise a wavefront aberration the DM offers. The fitting capabilities of the DM for each Zernike order are shown in Fig. 9.
From Fig. 9 it is obvious that the DM has good fitting capability for the first 7 orders of wavefront aberrations. Generally, lowerorder aberrations have the most important effect on intensity distribution and beam quality; therefore, we take only the first 7 orders into account. Note that a DM with higher spatial resolution and more actuators can accurately match more types of aberration modes. In our simulations, we introduce an incident wavefront with lower coupling efficiency (shown in Fig. 10).
Our main target is to find the best solution group (
u _{1}^{*},u _{2}^{*}) to fit the Zernike coefficients (v _{1},v _{2}) with the MBTS algorithm. We select the initial solutions as (0, 0), and the simulation results are shown in Fig. 11. From Fig. 11(a) we know that the searching path is an uneven curve that converges to the final point, namely the best solution(u _{1}^{*},u _{2}^{*}). We also find that the final solution group (u _{1}^{*},u _{2}^{*}) is (1.25, 0.75) from Fig. 11(c) and (d), that the RMS of the residual wavefront phase is about 0.2232 from Fig. 11(b), and that the corresponding coupling efficiency is about 80%, which is a significant improvement. In addition, the number of iterations for this search is about 35 to 40. However, we can use another method to compute the best solution (u _{1}^{*},u _{2}^{*}). The tabu list can record the latest local optimuma selected by the MBTS process. We also can calculate the mean values of the records in the current tabu list during each iteration; the result can be considered as the current best solution in the latest iteration. The results are shown in Fig. 10(e). Note that there are no improvements in finding a better final result, except that the curves are smoother. The residual wavefront phase is shown in Fig. 11(f).In Fig. 12, we compare the effects on this search process of starting with the different initial solutions (0, 0), (1, 1), (2, 2), and (4, 4) respectively. From Fig. 12 we find that different initial solutions lead to different search rates: The initial solution (4, 4) offers a faster search process (about 30 iterations), (2, 2) takes second place (35 iterations), and (1, 1) and (0, 0) do less well (about 40 iterations). The choice should be based on some practical experience and several tests.
We also perform a fitting simulation of the remaining Zernike coefficients from
v _{3} tov _{7}. The result is shown in Fig. 13. Here the initial solution is (4, 4), and the final solutions from Fig. 13 areu _{3}^{*} = 0.5,u _{4}^{*} = 0.35,u _{5}^{*} = 0.35,u _{6}^{*} = 0.3, andu _{7}^{*} = 0.3. The final residual wavefront aberration phase is shown in Fig. 14 with the corresponding RMS value of 0.1379. We can see that the coupling efficiency is improved significantly, from 42.81% to about 99%.The variation of coupling efficiency with iteration number is shown in Fig. 15. We select the four different initial solutions mentioned above and find that the final coupling efficiency reaches about 99% in about 30 iterations.
We also perform some experiments to test the SPGD algorithm; the results before and after correction are shown in Figs. 16 and 17 [21]. In these experiments, we know that the MBTS algorithm is much faster than the conventional SPGD algorithm, usually requiring hundreds of iterations (Fig. 16), while the final coupling efficiency for both algorithms is similar. But the iteration number is largely inferior to that for MBTS even in theory, since the highdimensional search space of SPGD reduces the searching rate, while in the MBTS algorithm we project a space with high dimension to a space with lower dimension and decouple them using a new system architecture. Thus the number of iterations will be sharply reduced, but only if the initial point is appropriately selected.
In Fig. 17, from the images received by CCD we see that the initial coupling efficiency before correction is very low, so that light scintillation is serious, but after correction the coupling efficiency is greatly improved. The corresponding results are in Table 1, and are inferior to the simulation results mainly because of the poor stability of the DM. This means that the DM should generate the identical corrected phase as by the same voltage applied to the actuators, in theory, but in practice it is limited by the skillful processing.
We know that the singleshot correction assumption is a concern for the algorithm convergence speed and atmospheric coherenttime. The SPGD is a typical example in [13], with about 500 iterations in experiment, corresponding to 50 ms. The operation rate is limited by the MEMs DM (78 kHz) to satisfy the requirements of the AO system correction. Therefore, the corresponding convergence time of the MBTS is 5 ms to 6.25 ms with the new system model, but with some other algorithms (SPGD being the typical one), the time is tens of milliseconds. In the simulation, the time spent by the SPGD is about 90 ms when the coupling efficiency reaches 0.8; in our experiment, it takes 5060 ms.
It is obvious that the SPGD works far less nicely than the MBTS. However, such a fast search (about 40 iterations) with the MBTS algorithm is at the cost of the construction of complicated equipment, with many sensorless AO subsystems working simultaneously. It places stringent demands on the experiment, and requires synchronization. This is the first such idea and system proposed in this research area, that many wavefront correctors may work in parallel to achieve better performance.
IV. CONCLUSION
In this study we propose a novel algorithm, MBTS, to compensate the wavefront aberration from atmospheric turbulence in an FSO communication system based on a new parallel correcting system architecture. The new idea of the proposed algorithm and its results in an FSO system are described in detail. Simulations indicate that our proposed MBTS can offer a faster search process during optimization (dozens of iterations) and better final coupling efficiency, compared to the SPGD algorithm. The architecture of a novel and interesting parallel correcting system with multiple adaptive optical (AO) subsystems is also given.

[FIG. 1.] Effect of aberrant wavefront on intensity: (a) laser with plane wavefront, (b) laser with aberrant wavefront.

[FIG. 2.] FSO system with sensorless AO.

[FIG. 3.] Layout of actuators (filled circles) for a 32element deformable mirror. The layout of the actuators is a square arrangement.

[]

[]

[]

[]

[]

[]

[]

[FIG. 4.] Parallel processing system with multiple sensorless AO subsystems.

[FIG. 5.] Twodimensional search space.

[FIG. 6.] The eight adjacent groups for N(S).

[FIG. 7.] Flow chart of the TS algorithm.

[]

[]

[]

[]

[]

[]

[FIG. 8.] Set of aberration modes: (a) 12 Zernike modes (314), (b) correcting wavefront (314).

[]

[FIG. 9.] Fitting capabilities of the DM for each Zernike order.

[FIG. 10.] Introduced wavefront aberration: RMS 0.8484, coupling efficiency 42.81%.

[FIG. 11.] Compensation by MBTS algorithm: (a) all points explored during searching, (b) RMS values of residual wavefront for different iterations, (c) solutions of u1* for different iterations, (d) solutions of u2* for different iterations, (e) mean values of current tabu list, (f) residual wavefront after compensation.

[FIG. 12.] Comparison of different initial solutions: initial solution (a) (0, 0), (b) (1,1), (c) (2, 2), and (d) (4,4).

[FIG. 13.] Solutions for u3*  u7* .

[FIG. 14.] Final residual wavefront aberration phase.

[FIG. 15.] Coupling efficiency based on MBTS.

[FIG. 16.] Coupling efficiency based on SPGD.

[FIG. 17.] Images captured by CCD (a) before and (b) after correction.

[TABLE 1.] Experimental results using SPGD