1 Introduction
In our previous paper[
1] we studied extraction characteristics of one-dimensional Talbot arrays. In this model each emitter is described by a complex field differential equation in the Rigrod approximation. All emitters are coupled through the Talbot cavity which constitutes one-half of the two-point resonator boundary conditions. We showed that the far-field performance of linear arrays is a strong function of the Talbot cavity length displaying multiple maxima and minima as a function of Talbot cavity length. This is true for either internal or external cavity phasing. This complicated behavior is due primarily to the Talbot propagation phase embodied in the interference between emitters. Also, we showed that the best far-field performance was in the half-Talbot plane for external phasing; the quarter-Talbot plane gives almost the same performance[
1]. When a linear stochastic propagation phase for each emitter is included we showed analytically that the standard threshold equation is replaced by a stochastic threshold equation which on the average predicts an increase in threshold gain. Further more, in the linear array development we introduced two measures of performance: a locking probability, and a coherence function. The former shows the number of times the integrator converges for a given rms phase
σ
2. This represents a measure of how often a lasing mode will be established for a certain rms optical path difference between the emitters. The coherence function is a normalized on-axis intensity.
In the following two-dimensional hexagonal case we will see entirely different behavior. There are four major differences. First, for a deterministic hexagonal Talbot cavity the far-field intensity as a function of the Talbot cavity length does not show any complicated structure, such a multiple maxima and minima as in the linear array case. Specifically, the far-field on-axis intensity shows a monotonic decrease as the Talbot cavity length increases. Next, when the stochastic propagation phase is included the coherence function shows an exp(-
σ
2) behavior independent of the Talbot plane. In the linear array the coherence function changes dramatically depending on the Talbot cavity length. Third, again as the Talbot cavity length changes the coherence function, for small
σ
2 remains near 19, the number of emitters in a double-hexagonal array, where as the coherence function varied widely in the linear case[
1]. Fourth, in our previous paper[
1] we showed that the one-dimensional linear array starts in an threshold eigen-state, for example the out-of-phase supermode, and remains in this mode as the gain is increased. However, the hexagonal array threshold mode is not necessarily maintained above threshold. Instead, we find the lasing mode becomes uniform characterized by nearly equal amplitudes and phases for all emitters. These features indicate the dominance of an in-phase lasing mode which is due to rotational symmetry. Another difference between these two cases is in the losses. The losses for a hexagonal array are larger than that a linear array with an identical number of emitters. For example, in a linear arrays there are always just two end emitters, but in a single hexagonal array of say 7 emitters, six are peripheral and contribute to the losses. This leads to a larger threshold and smaller extraction. In contrast, many identical nearest neighbors in the hexagonal array makes it more stable. Specifically, the probability of locking does not decrease as rapidly as a function of the rms phase as that of the linear array. We clearly see this behavior in our simulations. Briefly then, the above-threshold hexagonal array lases in a mode which is almost in-phase and can, therefore, be considered to be naturally phased.
2 Theory
In our previous paper[
1] we presented the theory of coupling between
N emitters in a Talbot cavity when each emitter had a stochastic linear propagation phase due to different fiber lengths and dispersion. Rather than repeat than development we will just present a synopsis.
We consider N gain elements coupled through a Talbot cavity terminated by a monolithic mirror of field reflectivity rt
. The jth gain element has outcoupling field reflectivity rj
and gain gj
. With these assumptions, and in the Rigrod approximation, the jth gain element supports a forward
(z) field and a reverse field (z) both described by
Since the gain is real the phases satisfy
and the constant Cj
is given by
These four equations constitute our basic working equations.
Fig. 1. Double hexagonal array and coordinate nomenclature
At the outcoupling end, z=0, we have
(0)=rj
(0). At the Talbot cavity end of the laser array, z=L, continuity of the electric fields requires that the reverse field is composed as (Li
)=∑Ri,j
(Lj
). Ri,j
is the complex reflection matrix developed by evaluating the overlap integral between the Fresnel propagated electric field of the jth
emitter integrated over the aperture of the of the ith
laser, see the Appendix. These two restrictions constitute the two-point boundary conditions.
When the fields contain a random linear propagation phase the Talbot cavity boundary condition leads to the stochastic threshold equation[
1]
which reduces to the deterministic form[
2
David Mehuys, William Streifer, Robert G. Waarts, and Davie F. Welsh, “Modal analysis of linear Talbot-cavity semiconductor lasers,” Opt. Lett.
16, 823–825(1991). [CrossRef] [PubMed]
] for
σ
2=0 and for
σ
2≠0 this leads to an increased threshold gain.
As a measure of locking we define the coherence function C as
which is related to the far-field on-axis intensity
I(0) by
I(0)=
C∑|
Am
|
2. For completely coherent emitters this function has a value
N equal to the number of emitters. If, however, the phases are random then the average of
eq. (6) takes a simple form when the amplitude spectrum is also flat. With these restrictions using and the previous averaging procedure
eq. (6) averages to
which approaches N exp(-σ
2) for large N.
3 Examples
Our example is a double hexagonal array with 19 elements as shown in
fig. (1) along with our coordinate nomenclature. The period of a typical V-groove fiber holder is
d=150
µ; the emitting aperture
a is 10
µm. The monolithic Talbot mirror is 100% reflecting while outcoupling field reflectivities are
rj
=
r=.8, also we pick a value of
gj
=
g=4 for the gain.
We begin with a short discussion of previous results obtained for the hexagonal array Talbot planes[
3
V. P. Kandidov, A. V. Kondrat’ev, and M. B. Surovitskii, “Collective modes of two-dimensional laser arrays in a Talbot cavity,” Quant. Elect.
28, 692–696(1998). [CrossRef]
]. These authors[
3
V. P. Kandidov, A. V. Kondrat’ev, and M. B. Surovitskii, “Collective modes of two-dimensional laser arrays in a Talbot cavity,” Quant. Elect.
28, 692–696(1998). [CrossRef]
] consider Talbot cavity properties of oblique-angled arrays, in general. In particular, for a hexagonal array they show that the paraxial quadratic propagation phase factor is
. The x-spatial frequency is
νp
=2
πp/d, the y-spatial frequency is
µp,q
=4
π(
q +
p/2)/√3
d,
p, q are integers,
d is the x-period, and
k=2
π/λ. (Note that the linear array results are obtained for
µp,q
=0.) Thus, for the hexagonal array the self-replication Talbot plane is when
. This has the solution
where the quantity in brackets is the Talbot distance for a linear array. The quadratic phase also shows that there is no precise out-of-phase plane similar to the linear case where the phase alternates between (0,
π) for adjacent emitters. Since this phase factor analysis does not yield outstanding fractional Talbot cavity lengths we turn to our threshold results to indicate important intermediate planes, see
eq. (8).
An analysis of the threshold equation,
eq. (8), shows that there is a mode which is analogous to the linear array out-of-phase solution for
z=2/3
zt
. Here, however, the phases do not alternate between (0,
π) but rather have the structure: (0,
π, 0); (
π,
π,
π,
π); (0,
π, 0,
π, 0) corresponding to the rows shown in
fig. (1). We refer to this plane as the out-of-phase Talbot plane. Similar to the linear case we look for an in-phase solution at 1/2 of the out-of-phase Talbot distance at 1/3
zt
. Here, we find that the phases are small and have the structure :(.076, .078, .076); (.079, .086, .086, .079): (.076, .086, .098, .086, .076). Thus, this solution is almost in phase with a very small phase difference over the array. This solution acts in the far-field as an in-phase eigenvalue solution. We will consider these two planes in more detail in a moment.
Threshold is determined by solving
obtained from
eq. (5) by setting
σ
2=0. The threshold gain
ḡi
is related to the eigenvalue
i
by
ḡi
=ln(1/
r)+ln(1/|
i
|) where
r is the field outcoupling reflectivity. In
fig. (2) we plot the three smallest threshold eigenvalue gains, ln(1/|
i
|), along with the largest threshold eigenvalue gain all as functions of the normalized Talbot cavity length
z/zt
for a field reflectivity of
r=.8. This figure does not indicate the threshold supermode phase structure.
Fig. (3) contains this information through the threshold eigenfunction coherence function
C, see
eq. (6), along with a repeat of the threshold eigenvalue gain for all 19 modes.
Fig. (3a) is for in the in-phase plane (1/3
zt
) and
fig, (3b) is for our out-of-phase plane (2/3
zt
). Now we return to
fig. (2). This figure shows that the three lowest threshold eigenvalue gains are within 10% of one another; the only exception is the 1/3 plane which is the in-phase solution. With equal importance,
figs. (3a,b) show that one of the lowest gains is an in-phase solution. In
fig. (3a) the lowest threshold gain is the in-phase mode. However, in
fig. (3b) the lowest threshold gain is the out-of phase mode with a gain of 2.02 while the next highest gain is the in-phase solution with a gain of 2.25 with
C=.33 and 18.6, respectively.. Further more, when this coherence and threshold analysis is repeated for intermediate planes it again shows that one of the lowest two gains is an in-phase solution. Thus, because the threshold gains are nearly equal the lasing mode is not necessarily the threshold mode. In fact, the following simulations show that this is the case. This is not true for the linear array where the threshold eigenvalues are well separated and the lasing mode is the same as the threshold mode. Specifically, in the out-of phase plane the lowest threshold gain is, of course, the out-of phase mode, however, the next higher gain in this plane is about a factor of 4 greater[
1].
Fig. 2. The three lowest threshold eigenvalue gains (colored) and the bounding largest threshold eigenvalue gain (black) as functions of the normalized Talbot length.
Fig. 3. (a) the coherence function (green) and the threshold eigenvalue gain (black) for the 1/3 Talbot plane;(b), the coherence function and the threshold eigenvalue gain for the 2/3 Talbot plane. Both are functions of mode number.
We now turn to the above threshold extraction by solving
eqs. (1) for the 19 coupled amplitude differential equations.
Figs. (4a,b) shows the on-axis intensity, and the Talbot cavity effective reflectivity as functions of the Talbot cavity length normalized to
zt
. The far-field on-axis intensity displays a smooth monotonic decrease as the cavity length increases and the coherence function remains at 19, but is not shown. Thus, the lasing modes are always in-phase affirming our earlier argument that the lasing mode is different from the threshold mode. This behavior is completely different from the on-axis intensity for a linear array which shows a very ragged dependence on the cavity length[
1]. In
fig. (4b) we also show the effective Talbot cavity reflectivity for an individual emitter defined by
Ri
=
Ā
i
(
L)/
(
L). One can easily see that the decreasing behavior is given by a (
z/
zt
)
-1 dependence due to the decreasing energy of the propagated field
Ep
in the Talbot cavity, see the Appendix. This
z/zt
dependence is the origin of the similar decrease seen in
fig. (4a). In this example
z
0≈.01cm,
zt
≈2cm for
d=150
µm,
ω
0=
a/√2=10/√2
µm and λ=1.55
µm. Note that 5%<
Ri
<56% for .2<
z/zt
<.8, and that
Ri
is independent of
i, as it is a ratio.
Fig. 4. (a), the far-field on-axis intensity; (b), the Talbot cavity effective reflectivity as functions of the normalized Talbot cavity length.
Optimum operating conditions are always of interest. For
z=1/3
zt
and
z=2/3
zt
,
fig. (5a) shows the outcoupled far-field on-axis intensity (1-
r
2)
I(0) as a function of the intensity reflectivity
r
2. Here,
I(0)=|∑
Ai
exp(
ϕi
)|
2 where
Ai
and
ϕi
are solutions of the coupled differential equations,
eqs. (1). For a gain of 4 and a fill factor of .047
fig. (5a) clearly shows an optimum outcoupling intensity reflectivity of
=18% and 20%, for
z/zt
=1/3 and
z/zt
=2/3 respectively.
Fig. (5b) shows the forward recirculating power for the central most emitter as a function of the outcoupling reflectivity
r
2. Here, the recirculating powers in the central emitter are similar except the 2/3 plane is reduced. Also, our simulations show that the effective reflectivity
Ri
, as defined above, is locked at 28% for all emitters over the entire range of
r
2 for the 1/3 plane and at
Ri
=7.2% for the out-of-phase plane where
z/zt
=2/3. For all cases the coherence function is a constant 19 indicating that the lasing mode is an in-phase mode. Comparing the graphs in
fig. (5a) we see that the optimum reflectivity is a mild function of Talbot cavity length. However, it is a stronger function of gain. For example, for
g=3,
=23%; and for
g=8, then
=10%.
Fig. 5. (a), the optimum far-field outcoupled on-axis intensity for the the 1/3 plane (red), and the 2/3 plane (black); (b) the forward recirculating power for the central emitter in the 1/3 plane (red), and in the 2/3 plane (black). All graphs are functions of the intensity reflectivity r
2.
Our last topic is the decreased performance when the emitters have a linear stochastic phase[
1]. Briefly, this is modeled by adding the same propagation phase to each emitter before and after each propagation. However, between all emitters the phase ensemble is random specified by some rms phase.
Eqs. (1) are then integrated until convergence or non-convergence for many such ensembles. Thus, we introduce the locking probability as the ratio of the number of times the integrator converges to the total number of tries for a specific rms phase. This then allows an assessment of the lasing performance in a stochastic environment.
Our simulations have shown that the hexagonal array lases in an in-phase mode. Thus, consideration of just one Talbot plane is sufficient and for this we pick our threshold out-of-phase plane located at
z/zt
=2/3. As a function of rms phase,
Fig. (6a) shows the coherence function;
fig. (6b) displays the outcoupled near-field intensity of the central emitter;
fig. (6c) shows the far-field on-axis intensity of the array; and
fig. (6d) shows the probability of locking. The average of the coherence function in
fig. (6a) satisfies
eq. (6) to within 10%. The decrease in the near-field emitter intensity,
fig. (6b) is due to a decrease in extracted energy as the rms phase increases. Thus the on-axis intensity,
fig. (6c), decreases much more rapidly than the coherence function due to the added decrease in the near-field extraction. Finally, the probability of locking remains near unity which means that the Talbot phase dominates. These figures are consistent because for a symmetric phased mode the central emitter near-field intensity of 1.3, shown in
fig. (6b), times 19
2 gives the far-field on-axis intensity of approximately 450, shown in
fig. (6c), all for a zero rms phase. This, again, conforms to an in-phase solution. These results are typical for any Talbot cavity length.
We conclude our analysis with a consideration of the most favorable fill factor by just plotting the amplitude of
Ri,j
for nearest, second-nearest, and third-nearest neighbors. The strong coupling region occurs when these three amplitudes are nearly equal. For the hexagonal array coupling to many neighbors beyond just the nearest neighbors requires the fill factor
f be less than .13 for the out-of phase (2/3
zt
) plane and less than .08 for the in-phase (1/3
zt
) plane. We do not show these graphs since they are similar to those in our previous work[
1].
We have only considered a double hexagonal array of 19 elements and did not present any results for a single array of just 7 elements. The reason is that the 7 element array solutions do not show the symmetry evident in the 19 element array and are more like the linear array results. In other words, the emitters of the 7 element array do not have as many identical neighbors as the 19 element array and thus the solutions are not as appealing.