K. HICKE et al. PHYSICAL REVIEW E 83, 056211 (2011)Compar<strong>in</strong>g the stability <strong>of</strong> the <strong>synchronization</strong> manifoldwithout delay mismatch [Eq. (8)] for L = K <strong>and</strong> with delaymismatch [Eq. (17)], we see that the variational equationsare identical, <strong>and</strong> thus similar ECMs [i.e., solutions S(t) withsimilar frequency <strong>and</strong> carrier density] <strong>of</strong> both systems willhave comparable stability properties. Thus we can suspect thatan asymmetrically placed mirror does not have a large effecton the <strong>synchronization</strong> properties <strong>of</strong> the system, provided thatL = K. When simulat<strong>in</strong>g the system numerically, we <strong>in</strong>deedf<strong>in</strong>d identical <strong>synchronization</strong> <strong>of</strong> the two semiconductor lasers,with a time lag correspond<strong>in</strong>g to the difference <strong>in</strong> propagationtime between the lasers <strong>and</strong> the mirror, as shown <strong>in</strong> Fig. 7(upper panel). Also, <strong>in</strong> the cross-correlation function, the ma<strong>in</strong>correlation peak is shifted by τ (see Fig. 7, lower panel) asexpected.How does the peak height, i.e., the lag correlation, dependon the delay mismatch? In Fig. 8, the upper red curve depictsthe height <strong>of</strong> the correlation peak as a function <strong>of</strong> τ. Formost <strong>of</strong> the allowed range 0 τ τ, the correlation peakis large <strong>and</strong> the <strong>synchronization</strong> quality is high. The position <strong>of</strong>the mirror does not have a large effect on the <strong>synchronization</strong>level.If we consider both a mismatch <strong>in</strong> the transmission <strong>and</strong>reflection <strong>of</strong> the mirror <strong>and</strong> a delay mismatch (Fig. 9), we f<strong>in</strong>dthat the area around the diagonal L = K synchronizes well,but the correlation breaks down even for a small mismatch <strong>in</strong>the coupl<strong>in</strong>g, <strong>in</strong> contrast to the case τ = 0 (Fig. 2). In otherwords, if the mirror is asymmetrically placed, the system ismuch more sensitive to a coupl<strong>in</strong>g mismatch, as was impliedby Eq. (16) <strong>and</strong> discussed above.B. Dynamical regimesAlthough the delay mismatch does not affect the quality <strong>of</strong>the (lag) <strong>synchronization</strong> <strong>of</strong> the coupled lasers, it does havea significant effect on the dynamics <strong>in</strong> the <strong>synchronization</strong>FIG. 8. (Color onl<strong>in</strong>e) Cross-correlation C(t)atshiftt = τ<strong>of</strong> the output <strong>in</strong>tensities <strong>of</strong> both lasers (upper red l<strong>in</strong>e), <strong>and</strong> the peakat t = 2τ <strong>of</strong> the autocorrelation function <strong>of</strong> one laser (lower bluel<strong>in</strong>e), versus the delay mismatch parameter τ. The peaks co<strong>in</strong>cide:higher regularity (peak <strong>in</strong> the autocorrelation) means stronger crosscorrelation(peak <strong>in</strong> the correlation coefficient). Parameters areτ = 1000, α = 4.0, p = 1.0, μ = 0.26, T = 200, L = K = 0.1,β = 10 −5 .FIG. 9. (Color onl<strong>in</strong>e) Cross-correlation C(t)atshiftt = τ<strong>of</strong> the field <strong>in</strong>tensities <strong>of</strong> two coupled lasers <strong>in</strong> a configuration,accord<strong>in</strong>g to Fig. 6, vsL <strong>and</strong> K. In comparison with Fig. 2, thesystem is more sensitive to a mismatch <strong>of</strong> the coupl<strong>in</strong>g L − K. τ =20. Other parameters are τ = 1000, μ = 0.26, α = 4.0, p = 1.0,T = 200, β = 10 −5 .manifold itself. When synchronized (with a time lag), thetwo subsystems behave like one laser, subject to two equallystrong feedbacks with different delays. This <strong>in</strong>troduction <strong>of</strong> asecond delay changes the alignment <strong>of</strong> the ECMs. The ECMansatz E(t) = Ae iωt , n(t) = n leads, after elim<strong>in</strong>at<strong>in</strong>g A, toa transcendental frequency equation <strong>and</strong> an equation for thecarrier density,ω =−2K √ 1 + α 2 s<strong>in</strong>(arctan α − ωτ) cos(ωτ), (18)pμ − 2K[cos(ωτ) cos(ωτ)]n = . (19)1 + μThe modes <strong>of</strong> this system are located <strong>in</strong>side the area <strong>of</strong> theellipse that is the solution space for a laser subject to only onefeedback with a strength <strong>of</strong> 2K. We note that if r τ = τ 1 /τ 2 ∈Q, the modes lie on a closed curve, <strong>and</strong> irrational ratios r τresult <strong>in</strong> a more complicated mode spectrum. Tronciu et al.[12], who <strong>in</strong>vestigated the correspond<strong>in</strong>g case <strong>of</strong> lasers subjectto feedback from an <strong>in</strong>tegrated double cavity, found chaoticbehavior for lower feedback strengths than needed for thecase <strong>of</strong> s<strong>in</strong>gle feedback. The more complex mode alignmentwas postulated as a cause for a more complex dynamics. Thecomplex alignment <strong>of</strong> the modes for τ ≠ 0 only dependson the delay mismatch, not on the delay time τ. S<strong>in</strong>ce thedynamics is organized by the mode spectrum, it is strongly<strong>in</strong>fluenced by the changes <strong>in</strong> the mode structure.For a small mismatch, the double delay setup exhibits afilter<strong>in</strong>g effect similar to a Michelson <strong>in</strong>terferometer [51],s<strong>in</strong>ce the frequencies for which ωτ ≈ (2m + 1) π 2 , m ∈ Zare filtered out. This effect results <strong>in</strong> the formation <strong>of</strong> “modeisl<strong>and</strong>s.” Close to simple rational ratios <strong>of</strong> r τ , we observesimilar phenomena.In Fig. 10, exemplary phase-space portraits are shown fordifferent values <strong>of</strong> the delay mismatch τ for a fixed me<strong>and</strong>elay time τ = 1000. The areas where the transversely stablemodes are located are similar for all delay mismatches τ.Dueto the different mode spectrum, the dynamics is significantlyaltered compared to the case τ = 0.A semiconductor laser subject to two self-feedbacks withdifferent strength <strong>and</strong> delay has been studied by Liu <strong>and</strong>Ohtsubo [52], who showed that one can stabilize the dynamicsto fixed po<strong>in</strong>ts or limit cycles. This effect was found strongestfor unequal feedback strengths. Rogister et al. [53] showed056211-6
MISMATCH AND SYNCHRONIZATION: INFLUENCE OF ... PHYSICAL REVIEW E 83, 056211 (2011)FIG. 11. (Color onl<strong>in</strong>e) Cross-correlation coefficient (upper redl<strong>in</strong>es), <strong>and</strong> auto-correlation at t = 2τ (lower blue l<strong>in</strong>es), vs the delaymismatch parameter τ for rational ratios <strong>of</strong> the delay times τ 1 <strong>and</strong>τ 2 . Magnifications <strong>of</strong> plot <strong>in</strong> Fig. 8. (a)r τ = τ 1τ 2= 1/1, (b) r τ = 2/1,(c) r τ = 3/1, (d) τ ≈ τ. Parameters as <strong>in</strong> Fig. 8.FIG. 10. (Color onl<strong>in</strong>e) Dynamics <strong>of</strong> the symmetrized solution(E S ,n S )<strong>in</strong>the(ω,n S ) phase space for different values <strong>of</strong> the delaymismatch parameter τ after t = 5 × 10 4 . For a def<strong>in</strong>ition <strong>of</strong> ω<strong>and</strong> n S , see Fig. 3. The ECMs for a laser subject to two delayedfeedbacks with different delay times are aligned <strong>in</strong>side the area <strong>of</strong>the mode ellipse for the case where the delays are equal. Green (lightgray) circles are transversely stable modes; red (black) triangles aretransversely unstable modes <strong>and</strong> antimodes. The it<strong>in</strong>eracy is drawn<strong>in</strong> blue (dark gray) <strong>and</strong> marked with an arrow if stabilized. The me<strong>and</strong>elay time τ is fixed at τ = 1000. The other parameters are α = 4.0,p = 1.0 (except lower right), μ = 0.26, K = 0.1. Note that the lowerright plot has a different y scale than the others.that chaotic dynamics <strong>and</strong> LFF can be suppressed for a s<strong>in</strong>glelaser, subject to two different feedbacks, by suppress<strong>in</strong>g theantimodes that are responsible for power dropouts [44]. Thisstabilization <strong>of</strong> the dynamics occurs ma<strong>in</strong>ly for a short secondfeedback. Increas<strong>in</strong>g the second feedback strength from a lowlevel to the magnitude <strong>of</strong> the first one results <strong>in</strong> a bifurcationcascade <strong>in</strong> the laser, lead<strong>in</strong>g to several dynamical regimes<strong>in</strong>clud<strong>in</strong>g stable behavior.In a configuration <strong>of</strong> two lasers coupled via a semitransparentmirror, we observe similar changes <strong>in</strong> the dynamics.We calculated the secondary peak <strong>of</strong> the autocorrelation <strong>of</strong>each laser at t = 2τ to <strong>in</strong>vestigate the regularity <strong>of</strong> thetime series. The height <strong>of</strong> the autocorrelation peak exhibitssignificant extrema (see Fig. 8, lower blue l<strong>in</strong>e) for certa<strong>in</strong>delay mismatches. In particular, we f<strong>in</strong>d dips <strong>and</strong> peaks <strong>in</strong>the vic<strong>in</strong>ity <strong>of</strong> simple rational values <strong>of</strong> r τ <strong>and</strong> also aroundr τ = 1(τ ≈ 0) <strong>and</strong> r τ →∞ (τ ≈ τ). At these po<strong>in</strong>ts,naturally the cross correlation also exhibits dips <strong>and</strong> peaks(see Fig. 11, upper red l<strong>in</strong>e). These are primarily caused bya change <strong>in</strong> the underly<strong>in</strong>g dynamical state <strong>of</strong> the system<strong>and</strong> not by a change <strong>of</strong> the <strong>synchronization</strong> quality. Loweredcross-correlation values have been tested to be caused by asmall signal-to-noise ratio when the dynamics is stabilized to afixed po<strong>in</strong>t. Delay mismatches that lead to a stabilization <strong>of</strong> thedynamics to a limit cycle result <strong>in</strong> a peak <strong>in</strong> the autocorrelation<strong>and</strong> cross correlation due to the coherent behavior.To further <strong>in</strong>vestigate the changes <strong>in</strong> the laser dynamics forvary<strong>in</strong>g delay mismatch, we compute a bifurcation diagram. Inthe diagram, which is shown <strong>in</strong> Fig. 12, several <strong>in</strong>tensity maxima<strong>and</strong> m<strong>in</strong>ima extracted from a time trace <strong>of</strong> each laser areplotted versus the delay mismatch τ. The dynamics changessignificantly for vary<strong>in</strong>g delay mismatch, when the delays areclose to a simple ratio, e.g., r τ = 1/1, 2/1, 3/2, .... EspeciallyFIG. 12. (Color onl<strong>in</strong>e) Bifurcation diagram <strong>of</strong> the <strong>in</strong>tensityextrema <strong>of</strong> the two lasers vs the delay mismatch parameter τ. Uppergreen po<strong>in</strong>ts are maxima; lower red po<strong>in</strong>ts are m<strong>in</strong>ima. Vary<strong>in</strong>g thedelay mismatch results <strong>in</strong> radical changes <strong>in</strong> the dynamics, which isevident, e.g., by the changes <strong>in</strong> the variance <strong>of</strong> the <strong>in</strong>tensity outputs.Magnifications are shown for a small delay mismatch τ < 15 (lowerleft), a delay times ratio <strong>of</strong> r τ ≈ 2/1 (lower middle), <strong>and</strong> one veryshort delay τ 1 ≫ τ 2 (lower right). Other parameters are τ = 1000,α = 4.0, μ = 0.26, p = 1.0, T = 200, β = 10 −5 .056211-7