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

Chaotic synchronization of a distant star-type laser network with multiple optical injections

Open Access Open Access

Abstract

A novel multi-injection module (MIM) is introduced into a typical distant star-type laser network, which is composed of a hub semiconductor laser node (H-SLN), star semiconductor laser nodes (S-SLNs) and tens of kilometers of fiber links. The chaotic synchronization of this distant network is investigated both experimentally and theoretically. As a result of using the MIM, a significantly low correlation (about 0.2) is successfully achieved between the H-SLN and S-SLNs in different clusters. This correlation is much lower than in previously reported results. Even when the fiber length is extended to 80 kilometers a low correlation (about 0.18) between the H-SLN and S-SLNs in different clusters is also obtained. Moreover, the dependence of chaotic synchronization on the operating conditions, such as the injection power, frequency detuning, and frequency mismatch between arbitrary nodes are examined. Lastly, using a theoretical model, we discuss the broad conditions for achieving chaotic synchronization among S-SLNs in the same cluster, and analyze the effect of the MIM branch number on chaotic synchronization.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Optical chaos has become a hot topic due to such excellent features as wide bandwidth, high complexity, and relatively easy implementation [14]. Novel applications of optical chaos have been found in secure communication [5], radar [6], physical random number generation (PRN) [7], time domain reflectometry [8], information processing [9] and even neural network [1012]. Especially, optical chaos in semiconductor lasers (SLs) have been investigated extensively. Synchronization between chaotic SLs is a necessary condition to achieve secure communication and secure key distribution [1,2,5], and various coupling synchronization regimes have been proposed [5,1323]. Moreover, cluster synchronization properties [2427] in networks with coupled SLs have caught more and more attention. According to the experimental phenomenon of cluster synchronization in coupled laser networks in [25], synchronization is achieved between lasers in the same clusters while there is no synchronization between lasers from different clusters. Moreover, the influence of network topology structures on synchronization is still needed to be further investigated. Xu et al. theoretically investigated the cluster synchronization of mutually coupled vertical-cavity surface-emitting lasers in networks with symmetric topology [26]. The performance of cluster synchronization in semiconductor lasers networks of 12 different topologies is studied both analytically and numerically [27].

Taking realistic conditions into account, the star-type laser network topology structure [2836] is extremely typical for more investigation. In a star-type laser network, side nodes (the star SLs nodes, S-SLN) within the same cluster are linked to a common central node (the hub SL node, H-SLN) in the other cluster. The chaotic synchronization between the nodes of a star-type laser network opens a route for a multi-directional communication scheme and advanced secure key distribution network [3136]. For instance, Bourmpos et al. investigated the parameter sensitivity of the synchronization in a star network for S-SLNs asymmetrically coupled to an H-SLN [33], and found the critical intrinsic and operational variables for zero-lag synchronization between side nodes [34]. Xiang et al. studied theoretically the synchronization regime of a star laser network with heterogeneous time delays [35]. Argyris et al. reported an experimental study of the ultra-high quality chaotic synchronization between the S-SLNs in a star-type SL network [36]. These systems could be classed as single-injection star-type SL (SIST-SL) networks, where the light of the H-SLN is injected into the corresponding S-SLNs through a single-coupling link, and the S-SLNs are classed into a cluster while the H-SLN formed a unique cluster. All the above results are very valuable for guiding the implementation of an SL network in the real world. However, when considering practical circumstances, both a high synchronization among S-SLNs and a low correlation between the H-SLN and S-SLNs are desirable. The former offers the basis for secure-oriented applications among S-SLNs, while the latter could ensure the system’s security [21-22, 30 and 32]. For a distant SL network, the long-distance fiber link between the H-SLN and S-SLNs could offer a viable attack avenue for eavesdroppers when the output of the Hub node is highly similar to that of the S-SLNs [22,30]. Unfortunately, in most SIST-SL networks the correlation between the H-SLN and S-SLNs is significantly high, typically about 0.7 [20,21]. Based on these considerations, it is necessary to explore new ways of reducing the correlation between the H-SLN and S-SLNs.

In this work, we propose a distant multi-injection star-type SL (MIST-SL) network, with a novel multi-injection module (MIM) to obtain multiple optical injections from the H-SLN to the S-SLNs. We build both the experimental setup and a theoretical model to analyze the performance of the MIST-SL network in terms of the temporal waveforms, the spectral profiles, and the cross-correlation properties between nodes in the same cluster and nodes among different clusters. This MIST-SL network owns some advantages compared with previous schemes in [21,22,30], such as all-fiber construction, the long-distance malleability, and easy implementability. In this work, we first describe the configuration of the MIST-SL network and its required synchronization properties. Secondly, we show the experimental results for tests of cluster synchronization among nodes. Thirdly, we explore and discuss the theoretical model of the MIST-SL network. Finally, we draw a number of conclusions.

2. System configuration

Figure 1 is the experimental setup of the MIST-SL network with a Hub SL node in one cluster (Cluster H, CH) and n Star SL nodes within the other cluster (Cluster S, CS). Numbers of 1550 nm InGaAsP/InP DFB-SLs and related optical components constitute the S-SLNs, and are linked to the H-SLN with tens of kilometers of single-mode fiber spools. In the H-SLN section, a chaotic signal is preliminarily generated in the SL via an external optical feedback. Next, the chaotic signal is transmitted through a multi-injection module (MIM) where the signal is divided into m parts. Each part of the signal passes through a non-isometric fiber branch, and the parts are combined again as a newly restructured signal at the MIM exit port. After that, the restructured signal is amplified by the EDFA and injected simultaneously into n S-SLNs through a 1: n optical coupler. Thus, the multiple optical injection from the H-SLN to S-SLNs is formed. It should be noted that the multi-injection here is different from that of previous star-type SL networks, where the optical signal from the H-SLN is directly injected into the S-SLNs in a single fiber link. In a MIST-SL network, the MIM is the key for the achieving a particularly low correlation among different clusters, and the high synchronization between the S-SLNs within the same cluster.

 figure: Fig. 1.

Fig. 1. Experimental setup of the MIST-SL network. H-SLN: Hub SL node. S-SLN: Star SL node. OC: optical circulator. PC: polarization controller. VA: variable attenuator. OI: optical isolator. MIM: multi-injection module. EDFA: erbium-doped fiber amplifier. PD: photodetector. OSA: optical spectrum analyzer. ESA: electronic spectrum analyzer. OPM: optical power meter.

Download Full Size | PDF

In the S-SLNs sections, all the SL chips have similar internal parameters and matched characteristics because they are selected from a common wafer. Series of variable attenuators (VAs) and polarization controllers (PCs) are used for adjusting the optical injection power (κ) and optical polarization of the H-SLN. Additionally, optical isolators (OIs) are used to prevent unwanted reflection disturbances. Moreover, a series of 10:90 optical couplers (OCs) are also inserted into the setup to sample light for monitoring and testing. Finally, test equipment is employed to measure the chaotic signals both optically and electronically. This equipment consists of an optical power meter (OPM), optical spectrum analyzer (Anritsu MS9740A, 0.02nm resolution), electronic spectrum analyzer (Keysight N9010A, 26.5GHz) and wideband digital oscilloscope (Keysight DSO91204A, 12GHz bandwidth). All SLs are driven by low noise SL controllers. In the experiment, m is set as 3 and n is 2 due to the limitations of laboratory equipment. The length of three fiber branches are set as about 1000 mm, 1409 mm and 1613 mm, respectively. The chip bias currents are fixed at 22 mA for the H-SLN, 11.38 mA for S-SLN1, and 11.32 mA for S-SLN2. The injection power κ from the H-SLN to each S-SLN is initially set at 170 µw. The chip temperatures are stabilized at 21.02°C for the H-SLN, 24.33°C for S-SLN1, and 24.68°C for S-SLN2. Under these circumstances, the laser free-running wavelengths are measured as 1553.483 nm for the H-SLN and 1553.514 nm for the S-SLNs. As a result, we obtain a positive frequency detuning ΔfH-S=3.85 GHz between the H-SLN and S-SLNs. These parameters are optimally chosen for chaotic synchronization among the SL nodes.

3. Experimental results

Figure 2 shows the experimentally recorded chaotic temporal waveforms and frequency spectra of lasers in the MIST-SL network with 10 kilometer fiber spools. In Figs. 2(a1)–2(a3), we see that the waveforms of the S-SLN1 and S-SLN2 in the CS are highly similar to each other, whereas the waveform of the H-SLN in the CH is quite different. Next, considering the spectra (both optical and electrical) in Figs. 2(b1)–2(b3), 2(c1)–2(c3), the spectra of the S-SLN1 and S-SLN2 are obviously similar, but the spectra of the H-SLN is different from that of the S-SLNs. For the low frequency part of RF spectra, we see that the signal in Fig. 2(c3) is significantly higher than in Figs. 2(c1) and 2(c2). This low frequency enhancement corresponds to the large temporal fluctuations in Fig. 2(a3), and results in the H-SLN temporal waveform being quite different from the S-SLN temporal waveforms. This indicates that there is significant synchronization between lasers within the same cluster and a low correlation among different clusters. Moreover, the bandwidths of S-SLN chaos signals can be further enhanced with careful optimization of laser operations parameters [20].

 figure: Fig. 2.

Fig. 2. The temporal waveforms (the left column), optical spectra (the middle column), and RF spectra (the right column) of the S-SLN1 (the first row), the S-SLN2 (the second row) and the H-SLN (the third row). The time delay between temporal waveforms has been compensated, and the total length of temporal waveforms is 2000 ns. The red curves in the RF spectra are the background noise floor.

Download Full Size | PDF

In order to assess the correlation degree between two temporal waveforms, the well-known cross-correlation (CC) coefficient is introduced by the following equation [2,3]:

$$C{C_{ij}}(\Delta \textrm{t}) = \frac{{\langle [{P_i}(t + \Delta t) - \langle {P_i}(t + \Delta t)\rangle ][{P_j}(t) - \langle {P_j}(t)\rangle ]\rangle }}{{{{\{ \langle {{[{P_i}(t + \Delta t) - \langle {P_i}(t + \Delta t)\rangle ]}^2}\rangle \langle {{[{P_j}(t) - \langle {P_j}(t)\rangle ]}^2}\rangle \} }^{1/2}}}}$$
where the subscripts i and j represent H-SLN, S-SLNs, respectively. Angle brackets represent time average, Δt is the time shift, P is the instance amplitude of temporal waveforms.

Figure 3 shows the correlation plots and calculated CC coefficient curves between arbitrary nodes of the MIST-SL network. From Figs. 3(a1)–3(a2), we can see that the max value of the CC coefficient (CCmax) reaches at 0.92, indicating good synchronization is achieved in the CS which is composed of S-SLN1 and S-SLN2. In contrast, in Figs. 3(b1), 3(b2), 3(c1), and 3(c2), the correlations between the H-SLN and S-SLNs are much lower, and the CCmax is calculated as 0.17 for Fig. 3(b2) and as 0.16 for Fig. 3(c2). That is to say, the correlations between the CH and the CS stay at a low level. These CCmax values are much lower than that of the single-injection scheme, where the typical CCmax values between the H-SLN and S-SLNs are about 0.7 [20,21]. The main general attack scenario for this type of network is that an eavesdropper might extract a clue by comparatively analyzing the signals of the H-SLN and S-SLNs through intercepting fiber links [22,30]. However, the correlation between different clusters in the MIST-SL network is so low that no useful clue can be extracted. Thus, this attack scenario can be eliminated.

 figure: Fig. 3.

Fig. 3. The correlation plots (the upper row) and CC coefficient curves (the lower row) between arbitrary nodes of MIST-SL network. a1-a2: the correlation between the S-SLN1 and S-SLN2 in the CS; b1-b2: the correlation between the S-SLN1 in the CS and H-SLN in the CH; c1-c2: the correlation between the S-SLN2 in the CS and H-SLN in the CH.

Download Full Size | PDF

Figure 4 presents the CCmax evolution of the MIST-SL system and the SIST-SL system under different injection powers κ and frequency detuning ΔfH-S between H-SLN and S-SLNs, to show the specific operation conditions for an SL network exhibiting the desired synchronous behavior. Firstly, for the MIST-SL situation, the CCmax values rise with the increase in injection power κ in Fig. 4(a1). The CCmax values between Star nodes (labeled by black squares) increase dramatically in the ‘I’ region (20µw<κ<150µw). The CCmax values between the Hub node and the Star nodes, on the other hand, grow slowly. Next, in the ‘II’ region (κ>150µw), CCmax between the Star nodes reaches levels over 0.9. Conversely, the CCmax values between the Hub node and the Star nodes retain a remarkably low level (less than 0.2). CCmax for different ΔfH-S values is given in Fig. 4(b1). In the whole investigated ΔfH-S range (from -20GHz to 20GHz), CCmax between the Star nodes reaches high levels in the “II” region (-12GHz<ΔfH-S<5GHz), while CCmax between the Hub node and the Star nodes always stay at a low level (less than 0.2). In comparison, Figs. 4(a2) and 4(b2) gives the evolution of CCmax between arbitrary nodes of the SIST-SL network, with similar operating conditions to those of the MIST-SL network. In Figs. 4(a2) and 4(b2), we see that CCmax between the Hub node and the Star nodes is over 0.6 in the “II” regions, indicating a strong correlation between the Hub node and Star nodes of the SIST-SL network. Therefore, the MIST-SL network can achieve good chaotic synchronization among the Star nodes within the same cluster and significantly low correlation between the Hub node in the CH and Star nodes in the CS for a relatively wide range of κ and ΔfH-S.

 figure: Fig. 4.

Fig. 4. The CCmax evolution under different injection powers κ (left column) and frequency detuning ΔfH-S (left column) between H-SLN and S-SLNs of the SL networks. The upper row is the MIST-SL case and the lower row is the SIST-SL case. The ΔfH-S is 3.85 GHz (a1, a2), the injection power is 170µw (b1, b2). The black squares denote the CCmax values between S-SLN1 and S-SLN2 in the CS; the red diamonds denote the CCmax between the H-SLN in the CH and S-SLN1 in the CS; the blue triangles denote the CCmax between the H-SLN in the CH and S-SLN2 in the CS.

Download Full Size | PDF

Next, it is necessary to consider the influence of parameter mismatch since it is inevitable in reality. Figure 5 shows the effects of ΔfS-S mismatch between S-SLN1 and S-SLN2 on chaotic synchronization. Here, ΔfS-S is introduced by fine tuning the wavelength of S-SLN1 (through changing the temperature of the laser chip) and fixing the wavelength of S-SLN2. Then, the synchronization evolution could be investigated between S-SLN1, S-SLN2 and H-SLN. According to Figs. 5(a1), 5(b1), and 5(c1), CCmax between the Star nodes in the CS is decreased sharply as the ΔfS-S mismatch increases. It demonstrates the high synchronization sensitivity of parameter mismatch in the clusters of the MIST-SL network. Figures 5(a1), 5(b1) and 5(c1) also show evidence that CCmax between nodes in different clusters preserves a low level (around 0.2) over the entire test ΔfS-S range. The blue triangles always keeps constant and a low level since the wavelengths of S-SLN2 and H-SLN are not change during the entire test ΔfS-S range. This feature promotes the security of chaos communication among the Star nodes within the same cluster [22,30]. The CCmax evolution in the SIST-SL network is very different from that of the MIST-SL network. As shown in Figs. 5(a2), 5(b2) and 5(c2), CCmax between the Hub node and Star nodes fluctuates significantly, and stays at a high level under specific ΔfS-S mismatch conditions, particularly for the positive ΔfS-S mismatch cases in Figs. 5(b2) and 5(c2).

 figure: Fig. 5.

Fig. 5. CCmax between arbitrary nodes as a function of ΔfS-S between S-SLN1 and S-SLN2 for different injection power levels. The injection power values are successively 50µw (a1), 100µw (b1), 140µw (c1), 100µw (a2), 150µw (b2), and 170µw (c2), respectively. The upper row is the MIST-SL case and the lower row is the SIST-SL case. The black squares denote the CCmax values between S-SLN1 and S-SLN2 in the CS; the red diamonds denote CCmax between the H-SLN in the CH and S-SLN1 in the CS; the blue triangles denote CCmax between the H-SLN in the CH and S-SLN2 in the CS.

Download Full Size | PDF

In order to illustrate the malleability of the transmission distance, Fig. 6 shows the chaotic synchronization among nodes in the cluster CH and CS of MIST-SL network with 80 kilometer fiber spools. The temporal waveform of the H-SLN (Fig. 6(a3)) is distinct from that of the S-SLNs (Fig. 6(a1) and Fig. 6(a2)). Accordingly, the remarkably low CCmax (0.16 and 0.18) between H-SLN of the CH and S-SLNs of the CS is achieved. In contrast, the chaotic synchronization among the Star nodes within the same cluster is demonstrated by the high similarity between the temporal waveforms (Fig. 6(a1) and Fig. 6(a2)), the spindly correlation plot (Fig. 6(b1)) and the calculated CCmax≈0.81. The synchronization among the Star nodes in the CS is distorted by the accumulated fiber dispersion and other nonlinearities [22,37]. Thanks to various advanced compensation devices, such as the dispersion compensating fiber (DCF), fiber Bragg grating (FBG), etc., the implementation of a long-distance MIST-SL network of over 100km is also very possible.

 figure: Fig. 6.

Fig. 6. The experimental results for the MIST-SL network with 80 kilometer fiber spools. a1-a3: the recorded temporal waveforms of S-SLN1, S-SLN2 and H-SLN, respectively. b1-a3: the correlation plots between S-SLN1 and S-SLN2 in the CS, between S-SLN1 in the CS and H-SLN in the CH, and between S-SLN2 in the CS and H-SLN in the CH, respectively.

Download Full Size | PDF

4. Theoretical model and related discussion

The SL dynamics can be described by the famous Lang-Kobayashi (L-K) equations [38]. The model of the MIST-SL network is then established by modifying the L-K equations with multi-optical injection terms. For the H-SLN of the CH:

$$\begin{array}{ll} \frac{{d{E_{HL}}(t )}}{{dt}} &= \frac{1}{2}({1 + i{\alpha_{HL}}} )\left( {{G_{HL}} - \frac{1}{{{\tau_P}}}} \right){E_{HL}}(t )+ {K_{HL}}{E_{HL}}({t - {\tau_{HL}}} )exp ({ - i2\pi {f_{HL}}{\tau_{HL}}} )\\ &\mathop {}\limits_{}^{} + \sqrt {2\beta {N_{HL}}} {\chi _{HL}} \end{array} . $$
$$\frac{{dNHL(t )}}{{dt}} = \frac{I}{e} - \frac{{NHL(t )}}{{{\tau _e}}} - GHL{|{EHL(t )} |^2}.$$
$${E_{H - SLN}} = 1/\sqrt m \times \left\{ \begin{array}{l} {}^{}{}^{}{E_{HL}}(t - {\tau_1}\textrm{ - }{\tau_0})\exp ( - i2\pi {f_{HL}}({\tau_1} + {\tau_0}))\\ + {E_{HL}}(t - {\tau_2}\textrm{ - }{\tau_0})\exp ( - i2\pi {f_{HL}}({\tau_2} + {\tau_0}))\\ + {E_{HL}}(t - {\tau_3}\textrm{ - }{\tau_0})\exp ( - i2\pi {f_{HL}}({\tau_3} + {\tau_0}))\\ \cdots \\ + {E_{HL}}(t - {\tau_m}\textrm{ - }{\tau_0})\exp ( - i2\pi {f_{HL}}({\tau_m} + {\tau_0})) \end{array} \right\}.$$
For the S-SLNs of the CS:
$$\begin{array}{ll} \frac{{d{E_{SLj}}(t)}}{{dt}} &= \frac{1}{2}(1 + i{\alpha _{SLj}})({G_{SLj}} - \frac{1}{{{\tau _p}}}){E_{SLj}}(t) + {K_{H - }}_S{E_{H - SLN}}exp ({i2\pi \Delta {f_{H - S}}t} )\\ & \mathop {}\limits_{}^{} + \sqrt {2\beta {N_{SLj}}} {\chi _{SLj}} \end{array} . $$
$$\frac{{d{N_{^{SLj}}}(t)}}{{dt}} = \frac{I}{q} - \frac{{{N_{SLj}}(t)}}{{{\tau _n}}} - {G_{SLj}}{|{{E_{SLj}}(t)} |^2}.$$
where E is the slowly varying optical field, and N represents the carrier number. The subscripts HL and SLj stand for the Hub laser and for the jth Star laser (j=1, 2, …, n) in the MIST-SL network. n is the number of Star nodes in CS. The multi-optical injection term EH-SLN is expressed as Eq. (4) to characterize the MIM in the experimental setup. The subscript m is the number of injection branches. The τ1, τ2, τ3 and τm are the time delays of each branch. τ0 is the fiber time delay from MIM to the Star lasers (S-SLNs). Next, in Eq. (5), the EH-SLN inject simultaneously into the S-SLNs. KH-S and ΔfH-S represent the injection strength and frequency detuning from the H-SLN to S-SLNs. For the convenience of calculation, the internal parameters of the SL lasers are assumed to be identical. τ0 sets as 0ns since the output of MIM is simultaneously injected into to S-SLNs. Then, we could focus on the effect of MIM on the chaotic synchronization. Next, the parameter values are: the linewidth enhancement factor α=3, the photon lifetime τp=2 ps, the carrier lifetime τn=2 ns, the electronic charge q=1.602×10−19 C, the normalized Gaussian white noise χ, the spontaneous emission rate β=1.1×103 s−1, the gain coefficient G = g(N(t)-N0)/(1+s|E(t)|2), where the differential gain coefficient s=1.5×10−8 ps−1, the transparent carrier number N0=1.5×108, and the gain saturation factor ρ=1×10−7. The time delay of each branch is τ1=4.9 ns, τ2=6.9 ns, τ3=7.9 ns. It is imaginable that the correlation between the H-SLN and S-SLNs should not be offset but be strengthened, if the time delays of fiber branches have integer multiple relationship with the external cavity feedback delay τHL. Therefore, in principle, the non-convention relationship between times delays in our setup is necessary. The feedback time delay is τHL=5.5 ns, the feedback strength is KHL=2.5 ns−1, and the pump current is I=29 mA. For simplicity, the free running wavelength of the lasers is set as λ=1550 nm (with a corresponding central frequency fHL = 1.94×1014Hz), and the value of KH-S is variable and is set to the optimal value. The software Matlab and fourth-order Runge-Kutta algorithms are employed to numerically simulate the dynamics of the SL nodes and the performance of the MIST-SL network.

To illustrate the characteristic of larger MIST-SL network, we theoretically simulate the synchronization when n is set as 2 (the left part) and 5 (the right part) in Fig. 7. One can see that the temporal waveforms of S-SLN1 and S-SLN2 (Figs. 7(a1) and 7(a2)) are nearly identical while the temporal waveforms of H-SLN (Fig. 7(a3)) are very different from that of S-SLNs. As a result, the CCmax between S-SLN1 and S-SLN2 (Fig. 7(b1)) is almost equal to 1, and the CCmax between H-SLN and S-SLNs are about 0.13 (Figs. 7(b2) and 7(b3)). The simulation results are confirmed with the experiment observations in Fig. 2 and Fig. 6. Moreover, for the n = 5 case (Figs. 7(c)–7(d)), the CCmax between laser nodes show similar property with that of n = 2 case. Based on the simulation and experiment, it could be reasonable to suppose that the MIST-SL network has similar performance for larger scale networks.

 figure: Fig. 7.

Fig. 7. The theoretical simulation of the MIST-SL network when the number of S-SLNs is set as 2 (the left part) and 5 (the right part). a1-a3: the temporal waveforms of S-SLN1, S-SLN2 and H-SLN, respectively. b1-b3: the correlation plots between S-SLN1 and S-SLN2 in the CS, between SLN1 in the CS and H-SLN in the CH, and between S-SLN2 in the CS and H-SLN in the CH, respectively. c1-c6: the temporal waveforms of S-SLN1, S-SLN2, S-SLN3, S-SLN4, S-SLN5 and H-SLN, respectively. d1-d6: the correlation plots between S-SLN1 and S-SLN2, between S-SLN2 and S-SLN3, between S-SLN3 and S-SLN4, between S-SLN4 and S-SLN5, between S-SLN5 and S-SLN1, between SLN5 and H-SLN, respectively.

Download Full Size | PDF

To analyze the general correlation feature among nodes of different clusters in a star-type SL network, Fig. 8 shows the simulated CCmax evolution maps in parameter spaces of KH-S and ΔfH-S. For the SIST-SL network case, as shown in Fig. 8(a1), the high synchronization between Star nodes in the CS (CCmax >0.99) region covers the right half of the parameter space. Additionally, in Figs. 8(b1) and 8(c1), the correlation between the Hub node and Star nodes in different clusters increases rapidly with increasing KH-S, and leads to a high correlation (CCmax >0.5) region marked by the white dash lines. We find that the high correlation regions in Figs. 8(b1) and 8(c1) almost coincide with the high synchronization region in Fig. 8(a1). Regrettably, it reveals that the appropriate KH-S and ΔfH-S range that satisfies the high synchronization between the Star nodes in the CS and the low correlation between the Hub node of CH and the Star nodes of CS in the SIST-SL network is very narrow. In comparison, for the MIST-SL network case, the high synchronization region almost covers the right half of the parameter space in Fig. 8(a2). Interestingly, in Figs. 8(b2) and 8(c2), the high correlation (CCmax >0.5) region disappears from the whole investigated parameter space. It means that the whole right half parameter space has the appropriate parameter range for the requirements of low correlation between the Hub node of CH and the Star nodes of CS in the MIST-SL network, along with the high synchronization quality between Star nodes of CS.

 figure: Fig. 8.

Fig. 8. The simulated evolution maps of CCmax in the parameter space of KH-S and ΔfH-S. The upper row is the SIST-SL case and the lower row is the MIST-SL case. The left column: the CCmax maps between S-SLN1 and S-SLN2 in the CS; the middle column: the CCmax maps between the H-SLN in the CH and S-SLN1 in the CS; the last column: the CCmax maps between the H-SLN in the CH and S-SLN2 in the CS. The white dash line marks the boundary with condition CCmax=0.99, while the black dash line marks the boundary with condition CCmax=0.5.

Download Full Size | PDF

Moreover, the effect of MIM branch number on the correlation among clusters is investigated in Fig. 9, where the number of the MIM branch varies from 1 to 6. As shown in Fig. 9(a), CCmax between Star nodes in CS suffer rare affection with the increasing number of the MIM branch. For a different MIM branch case, CCmax gets very high (CCmax>0.99) once KH-S>45ns−1. Next, in Figs. 9(b) and 9(c), CCmax monotonically decreases with the increase of MIM branch number under the same KH-S condition. The more branches are used in the MIM, the lower the CCmax received. Therefore, multi-injection is necessary for the simultaneous realization of low correlation between different clusters and high synchronization within the same cluster. There are several ways to further improve the system synchronization performance, such as manufacturing lasers with better intrinsic parameter consistency, better long-distance fiber pairs, and choosing consistent active devices in system etc.

 figure: Fig. 9.

Fig. 9. Simulated CCmax evolution between arbitrary lasers of the MIST-SL network with different MIM branch numbers as m=1, 2, 3, 4, 5, 6 respectively, a: the CCmax between S-SLN1 and S-SLN2 in the CS; b: CCmax between H-SLN in the CH and S-SLN1 in the CS; c: CCmax between H-SLN in the CH and S-SLN2 in the CS.

Download Full Size | PDF

5. Conclusions

In this paper, we experimentally and theoretically investigate the chaotic synchronization performance of a MIST-SL network, where a novel multi-injection module is introduced and lasers are linked with long distance fiber. Significant low correlation (corresponding CCmax<0.2) is successfully achieved between the H-SLN and S-SLNs in different clusters. We extensively investigate the effect of injection power and frequency detuning, as well as the frequency mismatch between laser nodes, and determine the specific operating conditions for desired synchronous outputs. Moreover, in the 80 kilometer distant MIST-SL network, high synchronization between S-SLNs in the same cluster and dramatically low correlation (CCmax≈0.18) between the H-SLN and S-SLNs in different clusters are also achieved, showing the distance malleability of MIST-SL network. Additionally, the theoretical simulations reveal the broad operation parameter range for the correlation suppression the Hub node and Star nodes among different clusters, and the effectiveness of the multi-injection module. From a broader perspective, the MIM could be transplanted into other topological types of laser networks, such as mesh, ring, or other hybrid types, and offers promise for various high applications, such as giant neurons networks, long-distance chaos communication networks or secure key distribution networks.

Funding

National Natural Science Foundation of China (11474233, 61875168); China Postdoctoral Science Foundation (2017M612885); Chongqing Postdoctoral Science Foundation Special Funded project (Xm2017008); Fundamental Research Funds for the Central Universities (XDJK2019B059).

Disclosures

The authors declare no conflicts of interest.

References

1. M. C. Soriano, J. García-Ojalvo, C. R. Mirasso, and I. Fischer, “Complex photonics: Dynamics and applications of delay-coupled semiconductors lasers,” Rev. Mod. Phys. 85(1), 421–470 (2013). [CrossRef]  

2. M. Sciamanna and K. A. Shore, “Physics and applications of laser diode chaos,” Nat. Photonics 9(3), 151–162 (2015). [CrossRef]  

3. J. G. Wu, S. W. Huang, Y. J. Huang, H. Zhou, J. H. Yang, J. M. Liu, M. B. Yu, G. Q. Lo, D. L. Kwong, S. K. Duan, and C. W. Wong, “Mesoscopic chaos mediated by Drude electron-hole plasma in silicon optomechanical oscillators,” Nat. Commun. 8(1), 15570 (2017). [CrossRef]  

4. N. Jiang, Y. J. Wang, A. K. Zhao, S. Q. Liu, Y. Q. Zhang, L. Chen, B. C. Li, and K. Qiu, “Simultaneous bandwidth-enhanced and time delay signature-suppressed chaos generation in semiconductor laser subject to feedback from parallel coupling ring resonators,” Opt. Express 28(2), 1999–2009 (2020). [CrossRef]  

5. A. Argyris, D. Syvridis, L. Larger, V. Annovazzi-Lodi, P. Colet, I. Fischer, J. García-Ojalvo, C. R. Mirasso, L. Pesquera, and K. A. Shore, “Chaos-based communications at high bit rates using commercial fibre-optic links,” Nature 438(7066), 343–346 (2005). [CrossRef]  

6. F. Y. Lin and J. M. Liu, “Chaotic radar using nonlinear laser dynamics,” IEEE J. Quantum Electron. 40(6), 815–820 (2004). [CrossRef]  

7. A. Uchida, K. Amano, M. Inoue, K. Hirano, S. Naito, H. Someya, I. Oowada, T. Kurashige, M. Shiki, S. Yoshimori, K. Yoshimura, and P. Davis, “Fast physical random bit generation with chaotic semiconductor lasers,” Nat. Photonics 2(12), 728–732 (2008). [CrossRef]  

8. Y. Wang, B. Wang, and A. Wang, “Chaotic correlation optical time domain reflectometer utilizing laser diode,” IEEE Photonics Technol. Lett. 20(19), 1636–1638 (2008). [CrossRef]  

9. L. Appeltant, M. C. Soriano, G. Van der Sande, J. Danckaert, J. Dambre, B. Schrauwen, C. R. Mirasso, and I. Fischer, “Information processing using a single dynamical node as complex system,” Nat. Commun. 2(1), 468 (2011). [CrossRef]  

10. S. Y. Xiang, Y. H. Zhang, J. K. Gong, X. X. Guo, L. Lin, and Y. Hao, “STDP-based unsupervised spike pattern learning in a photonic spiking neural network with VCSELs and VCSOAs,” IEEE J. Sel. Top. Quantum Electron. 25(6), 1–9 (2019). [CrossRef]  

11. S. Y. Xiang, Z. X. Ren, Y. H. Zhang, Z. W. Song, and Y. Hao, “All-optics neuromorphic XOR operation with inhibitory dynamics of a single photonic spiking neuron based on a VCSEL-SA,” Opt. Lett. 45(5), 1104–1107 (2020). [CrossRef]  

12. S. Y. Xiang, Z. X. Ren, Z. W. Song, Y. H. Zhang, X. X. Guo, G. Q. Han, and Y. Hao, “Computing primitive of fully VCSEL-based all-optical spiking neural network for supervised learning and pattern classification,” IEEE Trans. Neural Netw. Learn. Syst. 25(6), 1700109 (2020). [CrossRef]  

13. T. Heil, I. Fischer, W. Elsasser, J. Muler, and C. R. Mirasso, “Chaos synchronization and spontaneous symmetry-breaking in symmetrically delay-coupled semiconductor lasers,” Phys. Rev. Lett. 86(5), 795–798 (2001). [CrossRef]  

14. I. Fischer, R. Vicente, J. M. Buldu, M. Peil, C. R. Mirasso, M. C. Torrent, and J. Garcia-Ojalvo, “Zero-lag long-range synchronization via dynamical relaying,” Phys. Rev. Lett. 97(12), 123902 (2006). [CrossRef]  

15. W. L. Zhang, W. Pan, B. Luo, X. H. Zou, M. Y. Wang, and Z. Zhou, “Chaos synchronization communication using extremely unsymmetrical bidirectional injections,” Opt. Lett. 33(3), 237–239 (2008). [CrossRef]  

16. Y. Hong, M. W. Lee, J. Paul, P. S. Spencer, and K. A. Shore, “GHz bandwidth message transmission using chaotic vertical-cavity surface emitting lasers,” J. Lightwave Technol. 27(22), 5099–5105 (2009). [CrossRef]  

17. N. Jiang, W. Pan, L. S. Yan, B. Luo, W. L. Zhang, S. Y. Xiang, L. Yang, and D. Zheng, “Chaos synchronization and communication in mutually coupled semiconductor lasers driven by a third laser,” J. Lightwave Technol. 28(13), 1978–1986 (2010). [CrossRef]  

18. R. Lavrov, M. Jacquot, and L. Larger, “Nonlocal nonlinear electro-optic phase dynamics demonstrating 10 Gb/s chaos communications,” IEEE J. Quantum Electron. 46(10), 1430–1435 (2010). [CrossRef]  

19. A. Jafari, H. Sedghi, K. Mabhouti, and S. Behnia, “Slave-master dynamics of semiconductor laser with short external cavity,” Opt. Commun. 284(12), 3018–3029 (2011). [CrossRef]  

20. J. G. Wu, Z. M. Wu, G. Q. Xia, T. Deng, X. D. Lin, X. Tang, and G. Y. Feng, “Isochronous synchronization between chaotic semiconductor lasers over 40-km fiber links,” IEEE Photonics Technol. Lett. 23(24), 1854–1856 (2011). [CrossRef]  

21. J. G. Wu, Z. M. Wu, X. Tang, L. Fan, W. Deng, and G. Q. Xia, “Experimental demonstration of LD-based bidirectional fiber-optic chaos communication,” IEEE Photonics Technol. Lett. 25(6), 587–590 (2013). [CrossRef]  

22. J. G. Wu, Z. M. Wu, Y. R. Liu, L. Fan, X. Tang, and G. Q. Xia, “Simulation of bidirectional long-distance chaos communication performance in a novel fiber-optic chaos synchronization system,” J. Lightwave Technol. 31(3), 461–467 (2013). [CrossRef]  

23. S. Y. Xiang, Y. N. Han, H. N. Wang, A. J. Wen, and Y. Hao, “Zero-lag chaos synchronization properties in a hierarchical tree-type network consisting of mutually coupled semiconductor lasers,” Nonlinear Dyn. 99(4), 2893–2906 (2020). [CrossRef]  

24. L. Y. Zhang, W. Pan, L. S. Yan, B. Luo, X. H. Zou, and M. F. Xu, “Cluster synchronization of coupled semiconductor lasers network with complex topology,” IEEE J. Sel. Top. Quantum Electron. 25(6), 1501007 (2019). [CrossRef]  

25. M. Nixon, M. Friedman, E. Ronen, A. A. Friesem, N. Davidson, and I. Kanter, “Synchronized cluster formation in coupled laser networks,” Phys. Rev. Lett. 106(22), 223901 (2011). [CrossRef]  

26. M. F. Xu, W. Pan, S. Y. Xiang, and L. Y. Zhang, “Cluster synchronization in symmetric VCSELs networks with variable-polarization optical feedback,” Opt. Express 26(8), 10754–10761 (2018). [CrossRef]  

27. Y. A. Han, S. Y. Xiang, and L. Y. Zhang, “Cluster synchronization in mutually-coupled semiconductor laser networks with different topologies,” Opt. Commun. 445, 262–267 (2019). [CrossRef]  

28. W. L. Zhang, W. Pan, B. Luo, X. H. Zou, and M. Y. Wang, “One-to-many and many-to-one optical chaos communications using semiconductor lasers,” IEEE Photonics Technol. Lett. 20(9), 712–714 (2008). [CrossRef]  

29. J. Zamora-Munt, C. Masoller, J. Garcia-Ojalvo, and R. Roy, “Crowd synchrony and quorum sensing in delay-coupled lasers,” Phys. Rev. Lett. 105(26), 264101 (2010). [CrossRef]  

30. P. Li, J. G. Wu, Z. M. Wu, X. D. Lin, D. Deng, Y. R. Liu, and G. Q. Xia, “Bidirectional chaos communication between two outer semiconductor lasers coupled mutually with a central semiconductor laser,” Opt. Express 19(24), 23921–23931 (2011). [CrossRef]  

31. Y. Aviad, I. Reidler, M. Zigzag, M. Rosenbluh, and I. Kanter, “Synchronization in small networks of time-delay coupled chaotic diode lasers,” Opt. Express 20(4), 4352–4359 (2012). [CrossRef]  

32. H. Aida, M. Arahata, H. Okumura, H. Koizumi, A. Uchida, K. Yoshimura, J. Muramatsu, and P. Davis, “Experiment on synchronization of semiconductor lasers by common injection of constant-amplitude random-phase light,” Opt. Express 20(11), 11813–11829 (2012). [CrossRef]  

33. M. Bourmpos, A. Argyris, and D. Syvridis, “Sensitivity analysis of a star optical network based on mutually coupled semiconductor lasers,” J. Lightwave Technol. 30(16), 2618–2624 (2012). [CrossRef]  

34. M. Bourmpos, A. Argyris, and D. Syvridis, “Analysis of the bubbling effect in synchronized networks with semiconductor lasers,” IEEE Photonics Technol. Lett. 25(9), 817–820 (2013). [CrossRef]  

35. S. Y. Xiang, A. J. Wen, and W. Pan, “Synchronization regime of star-type laser network with heterogeneous coupling delays,” IEEE Photonics Technol. Lett. 28(18), 1988–1991 (2016). [CrossRef]  

36. A. Argyris, E. Pikasis, and D. Syvridis, “Highly correlated chaotic emission from bidirectionally coupled semiconductor lasers,” IEEE Photonics Technol. Lett. 28(17), 1819–1822 (2016). [CrossRef]  

37. G. P. Agrawal, Nonlinear Fiber Optics, San Diego, (Academic, CA, 1989).

38. R. Lang and K. Kobayashi, “External optical feedback effects on semiconductor injection laser properties,” IEEE J. Quantum Electron. 16(3), 347–355 (1980). [CrossRef]  

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. Experimental setup of the MIST-SL network. H-SLN: Hub SL node. S-SLN: Star SL node. OC: optical circulator. PC: polarization controller. VA: variable attenuator. OI: optical isolator. MIM: multi-injection module. EDFA: erbium-doped fiber amplifier. PD: photodetector. OSA: optical spectrum analyzer. ESA: electronic spectrum analyzer. OPM: optical power meter.
Fig. 2.
Fig. 2. The temporal waveforms (the left column), optical spectra (the middle column), and RF spectra (the right column) of the S-SLN1 (the first row), the S-SLN2 (the second row) and the H-SLN (the third row). The time delay between temporal waveforms has been compensated, and the total length of temporal waveforms is 2000 ns. The red curves in the RF spectra are the background noise floor.
Fig. 3.
Fig. 3. The correlation plots (the upper row) and CC coefficient curves (the lower row) between arbitrary nodes of MIST-SL network. a1-a2: the correlation between the S-SLN1 and S-SLN2 in the CS; b1-b2: the correlation between the S-SLN1 in the CS and H-SLN in the CH; c1-c2: the correlation between the S-SLN2 in the CS and H-SLN in the CH.
Fig. 4.
Fig. 4. The CCmax evolution under different injection powers κ (left column) and frequency detuning ΔfH-S (left column) between H-SLN and S-SLNs of the SL networks. The upper row is the MIST-SL case and the lower row is the SIST-SL case. The ΔfH-S is 3.85 GHz (a1, a2), the injection power is 170µw (b1, b2). The black squares denote the CCmax values between S-SLN1 and S-SLN2 in the CS; the red diamonds denote the CCmax between the H-SLN in the CH and S-SLN1 in the CS; the blue triangles denote the CCmax between the H-SLN in the CH and S-SLN2 in the CS.
Fig. 5.
Fig. 5. CCmax between arbitrary nodes as a function of ΔfS-S between S-SLN1 and S-SLN2 for different injection power levels. The injection power values are successively 50µw (a1), 100µw (b1), 140µw (c1), 100µw (a2), 150µw (b2), and 170µw (c2), respectively. The upper row is the MIST-SL case and the lower row is the SIST-SL case. The black squares denote the CCmax values between S-SLN1 and S-SLN2 in the CS; the red diamonds denote CCmax between the H-SLN in the CH and S-SLN1 in the CS; the blue triangles denote CCmax between the H-SLN in the CH and S-SLN2 in the CS.
Fig. 6.
Fig. 6. The experimental results for the MIST-SL network with 80 kilometer fiber spools. a1-a3: the recorded temporal waveforms of S-SLN1, S-SLN2 and H-SLN, respectively. b1-a3: the correlation plots between S-SLN1 and S-SLN2 in the CS, between S-SLN1 in the CS and H-SLN in the CH, and between S-SLN2 in the CS and H-SLN in the CH, respectively.
Fig. 7.
Fig. 7. The theoretical simulation of the MIST-SL network when the number of S-SLNs is set as 2 (the left part) and 5 (the right part). a1-a3: the temporal waveforms of S-SLN1, S-SLN2 and H-SLN, respectively. b1-b3: the correlation plots between S-SLN1 and S-SLN2 in the CS, between SLN1 in the CS and H-SLN in the CH, and between S-SLN2 in the CS and H-SLN in the CH, respectively. c1-c6: the temporal waveforms of S-SLN1, S-SLN2, S-SLN3, S-SLN4, S-SLN5 and H-SLN, respectively. d1-d6: the correlation plots between S-SLN1 and S-SLN2, between S-SLN2 and S-SLN3, between S-SLN3 and S-SLN4, between S-SLN4 and S-SLN5, between S-SLN5 and S-SLN1, between SLN5 and H-SLN, respectively.
Fig. 8.
Fig. 8. The simulated evolution maps of CCmax in the parameter space of KH-S and ΔfH-S. The upper row is the SIST-SL case and the lower row is the MIST-SL case. The left column: the CCmax maps between S-SLN1 and S-SLN2 in the CS; the middle column: the CCmax maps between the H-SLN in the CH and S-SLN1 in the CS; the last column: the CCmax maps between the H-SLN in the CH and S-SLN2 in the CS. The white dash line marks the boundary with condition CCmax=0.99, while the black dash line marks the boundary with condition CCmax=0.5.
Fig. 9.
Fig. 9. Simulated CCmax evolution between arbitrary lasers of the MIST-SL network with different MIM branch numbers as m=1, 2, 3, 4, 5, 6 respectively, a: the CCmax between S-SLN1 and S-SLN2 in the CS; b: CCmax between H-SLN in the CH and S-SLN1 in the CS; c: CCmax between H-SLN in the CH and S-SLN2 in the CS.

Equations (6)

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

C C i j ( Δ t ) = [ P i ( t + Δ t ) P i ( t + Δ t ) ] [ P j ( t ) P j ( t ) ] { [ P i ( t + Δ t ) P i ( t + Δ t ) ] 2 [ P j ( t ) P j ( t ) ] 2 } 1 / 2
d E H L ( t ) d t = 1 2 ( 1 + i α H L ) ( G H L 1 τ P ) E H L ( t ) + K H L E H L ( t τ H L ) e x p ( i 2 π f H L τ H L ) + 2 β N H L χ H L .
d N H L ( t ) d t = I e N H L ( t ) τ e G H L | E H L ( t ) | 2 .
E H S L N = 1 / m × { E H L ( t τ 1  -  τ 0 ) exp ( i 2 π f H L ( τ 1 + τ 0 ) ) + E H L ( t τ 2  -  τ 0 ) exp ( i 2 π f H L ( τ 2 + τ 0 ) ) + E H L ( t τ 3  -  τ 0 ) exp ( i 2 π f H L ( τ 3 + τ 0 ) ) + E H L ( t τ m  -  τ 0 ) exp ( i 2 π f H L ( τ m + τ 0 ) ) } .
d E S L j ( t ) d t = 1 2 ( 1 + i α S L j ) ( G S L j 1 τ p ) E S L j ( t ) + K H S E H S L N e x p ( i 2 π Δ f H S t ) + 2 β N S L j χ S L j .
d N S L j ( t ) d t = I q N S L j ( t ) τ n G S L j | E S L j ( t ) | 2 .
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.