Coherent beam combination of fiber lasers with a strongly confined tapered self-imaging waveguide: theoretical modeling and simulation

Rumao Tao; Xiaolin Wang; Hu Xiao; Pu Zhou; Lei Si

doi:10.1364/PRJ.1.000186

1. INTRODUCTION

Due to high conversion efficiency, excellent beam quality, convenient heat management, and compact configuration of fiber lasers [1], fiber lasers have wide applications, i.e., in industrial processing and optical communication. With the development of high-power laser diode pump technology and double-clad fiber production crafts, the output power of fiber lasers has been increasing rapidly in recent years. However, due to nonlinear effect, facet fracture, and thermal lens [2], the ultimate output power of a single fiber laser cannot increase unrestricted. A promising approach to overcome this difficulty is coherent beam combination (CBC) of multiple fiber lasers, which can achieve a high-power laser beam with good beam quality while maintaining excellent heat-managing capability [3–10]. To date, most of the prevalent architectures for CBC of fiber lasers have involved free-space phased arrays that incorporate multiple tiled emitters. A unity fill factor is required for free-space phased array architecture to ensure the far-field beam quality [4–10]. However, this is too difficult for practical engineering and much of the power in the central lobe is diverted to side lobes. Alternative solutions based on filled-aperture designs, such as those based on diffractive optical elements (DOEs) [11], coherent polarization beam combination (CPBC) [12], and reimaging assisted phased arrays (REAPAR) [13,14], offer much better combining efficiencies and thus reduce side lobes in the far field compared with the aforementioned architectures. The REAPAR architecture, based on the self-imaging properties of waveguides, was first developed by Christensen and Koski for CBC in planar waveguides [13]. Recently it was reported that, based on fiber laser array, the REAPAR technology has been extended to two-dimensional (2D) waveguides and successfully implemented at high power levels [14]. One potential complication during the implementation of the aforementioned configuration is the very close packing of the laser array, which will be on the order of hundreds of micrometers. Tapered waveguide is a special kind of waveguide which, expanding the launch aperture while maintaining the output aperture, can be used in CBC and provide convenient operation with large input beam spacing [13]. However, to the best of our knowledge, the theoretical model to analyze the tapered-waveguide-based CBC of fiber lasers has not been reported.

In the present paper, we discuss the self-imaging properties of fiber lasers in a strongly confined tapered waveguide (SCTW) and its potential application in CBC of fiber lasers. Analytical formulas are derived for the positions, amplitudes, and phases of the $N$ self-images at the output of a SCTW. A finite difference beam propagation method (FDBPM) is used to verify numerically the analytical expressions. Base on the analytical results and the FDBPM, the tapered-waveguide-based CBC of the fiber laser array is investigated, and the CBC results based on a tapered waveguide are compared with that based on a nontapered one. The combining performance of the system is also investigated, followed by some instructive discussions.

2. THEORETICAL MODEL

A. System Introduction

The CBC system based on a tapered waveguide is depicted in Fig. 1. The beams emerging from amplified seed source pass through phase modulators, and then through fiber amplifiers. The amplified beams are coupled into the transport fibers and form a fiber array, which passes through a microlens array. The microlens array is used to focus the beams into the combining tapered waveguide. A small portion of the combined output laser is picked off with a high-reflector and detected with a photodetector. The photodetector output is sent to the controller and is used as a closed-loop feedback signal to control the phase modulators. Then the amplified beams are coherently combined into one single beam.

Fig. 1. Schematic diagram of the waveguide-based CBC system.

Download Full Size | PDF

B. Imaging of a Fiber Laser through a SCTW

For the sake of simplicity, a tapered waveguide of homogeneous refractive index with symmetric cladding (as shown in Fig. 2) is considered, which provides one-dimensional self-imaging. $W_{1}$ is the input width, $W_{2}$ is the output width, and $θ$ is the taper angle between the $z$ axis and the tapered side of the waveguide (if $W_{1} > W_{2}$ , $θ$ is positive; on the contrary, $θ$ is negative). Figure 3 shows the amplitude distributions in SCTW in the case of a Gaussian beam input (which is a good approximation for fiber lasers). In the numerical simulation, the origin of the $x$ axis is taken to be at the center, which is taken to be on the one side as shown in Fig. 2 in derivation. We denote the characteristic distance of propagation, where the input amplitude profile is reconstructed (with some aberrations) in a loss-free, straight waveguide, as $L_{0}$ . It can be seen from Fig. 3 that, if the laser is incident at the offset position, the self-image can only happen at $z = L_{0}$ ; however, self-image can be achieved at $z \leq L_{0}$ when the laser is incident at the axially aligned position. It is advantageous when the fiber laser is incident at the axially aligned position since it can shorten the imaging length of the waveguide, which can save space and make the application more compact, so we mainly focus our attention on the latter situation.

Fig. 2. Schematic diagram of a tapered waveguide.

Download Full Size | PDF

Fig. 3. Amplitude distribution in the SCTW for Gaussian beam propagation. (a) Offset input Gaussian beam and (b) axially aligned input Gaussian beam.

Download Full Size | PDF

C. Imaging of a Tapered Waveguide

The theory is developed for arbitrary symmetric input light distribution and for a 2D tapered waveguide. However, with the effective index method, which is reasonably accurate for strongly guided modes [15], the practical 3D problem can be reduced to two dimensions. The present results are relevant for square and rectangular self-imaging produced by 3D waveguides of square or rectangular cross sections [16,17]. The application of waveguide to coherent combining of lasers is the reverse use of the self-imaging effects [13,14]. We consider the imaging of a symmetric beam through a 2D tapered waveguide to simplify the discussion. For a strongly confined waveguide, guided modes are almost completely confined within the waveguide and their lateral mode profiles contain an integer number of half-periods within the waveguide [17,18]. Defining the incident plane as the original plane and using a spatial Fourier decomposition, the incident light distribution $f (x, 0)$ can be rewritten as a superposition of the infinite number of strongly even guided eigenmodes $E_{i} (x, 0)$ with the coefficients $a_{i} (0)$ :

f (x, 0) = \sum_{i = 0}^{\infty} a_{2 i} (0) E_{2 i} (x, 0), a_{i} (0) = \frac{2}{W_{1}} \int_{0}^{W_{1}} f (x) E_{i}^{*} (x, 0) d x .

The symmetry properties of the input light distribution and the axially aligned incident properties are used to obtain Eq. (1). The strongly guided eigenmodes at plane

z

have the form of

E_{i} (x, z) = \sin [π (i + 1) \frac{x}{W (z)}], i = 0, 1, 2, \dots,

where

i

is the mode number and

W (z)

is the active width of the waveguide.

W (z)

is equal to the physical thickness

W_{1} - 2 z \tan θ

, which is slightly corrected at both sides by the Goos–Hahnchen penetration depth in the cladding of index

n_{c} < n

[15,16,18]:

W (z) = (W_{1} - 2 z \tan θ) + 2 (\frac{λ}{2 π}) {(\frac{n_{c}}{n})}^{2 σ} {(n^{2} - n_{c}^{2})}^{- 1 / 2},

where

σ = 0

for TE polarization,

σ = 1

for TM.

n

is the index of the waveguide.

The transverse propagation constants are

k_{x i} (z) = (i + 1) π / W (z),

and the corresponding longitudinal propagation constants are

β_{i}^{2} (z) = n^{2} k_{0}^{2} - k_{x i}^{2} (z) .

With paraxial approximation, one can obtain

β_{i} (z) ≅ n k_{0} - k_{x i}^{2} (z) / 2 n k_{0},

or

β_{0} (z) ≅ n k_{0} - Δ β_{02} (z) / 8, Δ β_{02} (z) = β_{0} (z) - β_{2} (z) ≅ \frac{8 π^{2}}{2 n k_{0} W^{2} (z)},

β_{2 i} (z) ≅ β_{0} (z) - i (i + 1) \frac{Δ β_{02} (z)}{2},

where

k_{0}

is the propagation constant in vacuum,

k_{0} = (2 π / λ)

.

According to Eq. (6c), we can conclude that if the two lowest-order modes (zero-order and second-order) are coupled, all the other higher modes will fulfill the couple condition. If we set $L = L_{0}$ as the shortest coupling length (self-imaging length), then the following relation can be established as

\frac{8 π^{2}}{2 n k_{0}} \int_{0}^{L_{0}} \frac{d z}{W^{2} (z)} = 2 π .

Substituting Eq. (3) into Eq. (7), self-imaging length can be expressed as

L_{0} = \frac{W_{0}}{\frac{λ}{n W_{0}} + 2 \tan θ},

where

W_{0} = W (0)

.

The output field $f_{out} (x, L_{N})$ at plane $z = L_{N}$ (the subscript $N$ denotes the number of the self-images, which means the number of the fiber lasers that can be combined by the waveguide when the waveguide is used for CBC) is

f_{out} (x, L_{N}) = \sum_{i = 0}^{\infty} a_{2 i} (L_{N}) E_{2 i} (x, L_{N}) \exp (- j \int_{0}^{L_{N}} β_{2 i} (z) d z) .

Substituting Eqs. (5) and (2) into Eq. (9), one can achieve

f_{out} (x, L_{N}) = \sum_{i = 0}^{\infty} a_{2 i} (L_{N}) A_{2 i} E_{2 i} (x, L_{N}),

Equation (10b) can be calculated recursively:

A_{0} = \exp [- j [n k_{0} - \frac{π^{2}}{2 n k_{0} W_{0} (W_{0} - 2 \tan θ L_{N})}] L_{N}],

Then the coefficient $B_{2 i}$ is introduced, which is defined as

B_{2 i} = \frac{1}{C} \sum_{q = 0}^{N - 1} \exp (j Φ_{2 i, q}),

with

Φ_{2 i, q} = - π (2 i + 1) \frac{x_{q}}{W (L_{N})} - φ_{q},

According to Eqs. (12) and (13d),

A_{0} = B_{0}

.

Using Eqs. (13a)–(13c), one can obtain

Φ_{2 i, q - 1} = Φ_{2 (i - 1), q} + π \frac{1}{N} (2 i) .

Recalling the $N$ periodicity of the summands, the coefficients $B_{2 i}$ can be written as

B_{2 i} = B_{2 (i - 1)} \exp (j π \frac{1}{N} 2 i) .

Comparing Eq. (11) with Eq. (15), only if one of the following relations is fulfilled:

\frac{L_{N}}{(\frac{n W_{0}}{λ}) (W_{0} - 2 \tan θ L_{N})} = \frac{1}{N},

or

L_{N} = \frac{n W_{0}^{2}}{N λ + 2 n \tan θ W_{0}}, W (L_{N}) = W_{0} - \frac{2 n \tan θ W_{0}^{2}}{N λ + 2 n \tan θ W_{0}},

can the following relationship be established:

A_{2 i} = B_{2 i} = \frac{1}{C} \sum_{q = 0}^{N - 1} \exp (- j π (2 i + 1) \frac{x_{q}}{W (L_{N})} - j φ_{q}) .

Replacing the summation index $q$ of the sum with $N - q - 1$ , the above expression can be written as

A_{2 i} = B_{2 i} = \frac{1}{C} \sum_{q = 0}^{N - 1} \exp (j π (2 i + 1) \frac{x_{q}}{W (L_{N})} - j φ_{q}) .

Inserting Eq. (15) into Eq. (9), one can obtain

A_{2 i} E_{2 i} (x, L_{N}) = \frac{1}{2 j} {A_{2 i} \exp [j π (2 i + 1) \frac{x}{W (L_{N})}] - A_{2 i} \exp [- j π (2 i + 1) \frac{x}{W (L_{N})}]} .

Then

A_{2 i} E_{2 i} (x, L_{N}) = \frac{1}{C} \sum_{q = 0}^{N - 1} \frac{1}{2 j} {\exp [j π (2 i + 1) \frac{x - x_{q}}{W (L_{N})}] - \exp [- j π (2 i + 1) \frac{x - x_{q}}{W (L_{N})}]} \exp (- j φ_{q}) . = \frac{1}{C} \sum_{q = 0}^{N - 1} E_{2 i} (x - x_{q}, L_{N}) \exp (- j φ_{q}) .

The output field distribution can be rewritten as

f_{out} (x, L_{N}) = \frac{1}{C} \sum_{i = 0}^{\infty} a_{2 i} (L_{N}) \sum_{q = 0}^{N - 1} E_{2 i} (x - x_{q}, L_{N}) \exp (- j φ_{q}), = \frac{1}{C} \sum_{q = 0}^{N - 1} \exp (- j φ_{q}) \sum_{i = 0}^{\infty} a_{2 i} (L_{N}) E_{2 i} (x - x_{q}, L_{N}) .

According to Eq. (1), Eq. (20) can be rewritten as

f_{out} (x, L_{N}) = \frac{1}{C} \sum_{q = 0}^{N - 1} f (x - x_{q}, L_{N}) \exp (- j φ_{q}) .

According to Eq. (21), a sum of $N$ images with equal amplitudes numbered by $q = 0, 1, \dots, N - 1$ are obtained. Due to the fact that the waveguide is tapered, the images at the output of the waveguide ( $f (x, L_{N})$ ) are different from the initial input field ( $f (x, 0)$ ). Setting $θ = 0$ , the relations for nontapered waveguides can be obtained. Equations (3), (8), (13), and (16) are the main results obtained in this manuscript, which are the essential condition to realize the SCTW-based CBC. Combined with a FDBPM, the waveguide-based CBC of fiber lasers can be studied numerically.

3. NUMERICAL RESULTS

To verify the formulas derived in Section 2, numerical simulation examples are presented. We take the values determined by numerical techniques as numerically exact [19] and use them in the comparison to justify our formulas by numerical examples. A FDBPM [20,21] is used to simulate the light propagation in the tapered waveguide section. The computational core of our simulation program is described in [20,21], which has been improved based on the method described in [22–25] and referenced therein. The simulation results obtained with the FDBPM can provide reliable accuracy [26,27] for comparison. In the FDBPM simulation, we choose the discretization sizes $Δ x = 0.02 μm$ . A relative error $ζ$ is defined to characterize the length error of our analytical model, which is expressed as

ζ = \frac{| L_{N_simulation} - L_{N_analytical} |}{L_{N_simulation}} \times 100 %,

where

L_{N_analytical}

is the waveguide length predicted by our analytical model and

L_{N_simulation}

is the waveguide length obtained through FDBPM. The

L_{N_simulation}

is obtained by maximizing the power in the bucket (PIB) criterion, which indicates how much energy is transferred to the main lob and is defined as

P I B = \frac{\int_{- h}^{h} \int_{- h}^{h} 〈 I (r⃗, z) 〉 d x d y}{\int_{- \infty}^{\infty} \int_{- \infty}^{\infty} 〈 I (r⃗, z) 〉 d x d y},

where

h

is the radius of the bucket. With the

L_{N_analytical}

predicted by Eq. (16), we just need to analyze the simulation results near

L_{N_analytical}

to find more accurate results,

L_{N_simulation}

, which are time-saving and useful to reduce the cost in practical application.

We consider a SCTW with $W_{1} = 50 μm$ , $\tan θ = 0.0103$ , and $n = 1.45$ . The waveguide is situated in air, i.e., $n_{c} = 1.00$ ; $λ = 1.00 μm$ is taken here. Gaussian beams with waist width of 5 μm are launched at $z = 0$ . The relative small size of the waveguide and beam size compared to that used in the real experiments are due to the limitation of computation resources. Figure 4 shows the field distribution of the waveguide by the FDBPM simulation with “transparent” boundary conditions [28,29]. Equation (7) predicts that the self-imaging length is $L_{0} = 1465.5 μm$ and the self-imaging length given by the FDBPM simulation results is $L_{0} = 1466 μm$ [see Fig. 4(c)], with $ζ = 0.03 %$ . Choosing $N = 2$ with incident position and phase relation indicated by Eq. (13), we obtained $L_{2} = 1047.1 μm$ , and the self-imaging length given by the FDBPM simulation results is $L_{2} = 1048 μm$ [see Fig. 4(d)], with $ζ = 0.08 %$ . Figures 4(e)–4(h) give the field distribution at the simulated length. From Figs. 4(e) and 4(g), we can obtain that the self-imaging of the input laser beam is obtained with some degree of compression; from Figs. 4(f) and 4(h), the combining of the two Gaussian beams is achieved, which proves the accuracy of the incident position and phase relation indicated by Eq. (13). One can obtain from Figs. 4(g)–4(h) that, within the main lobe area, the phase difference between our analytical model (plane) and the simulation is negligible (less than $λ / 10$ ). It can be also found from Figs. 4(a) and 4(b) that the self-imaging properties of the tapered waveguide are distorted as indicated in Eq. (21), and the combined performance is degraded.

Fig. 4. Field distribution at the output of the SCTW. (a), (c), (e), and (g) are for the self-imaging of the laser beam; (b), (d), (f), and (h) are for the combination application.

Download Full Size | PDF

For perfect self-imaging, the field of the self-imaging beam is the same as that of the input beam from Eqs. (8) and (21). However, it reveals in Figs. 4(e)–4(h) that the fields of the laser beam are distorted. The phase distributions of the self-imaged beam fluctuate within a narrow range while the phase distributions of the input beam are plane. This can be understood easily from the analytical derivation. For complete mode decomposition in Eq. (8), we need an infinite number of guided modes totally confined within the waveguide. The resulting approximation errors have been investigated and are negligible for nontapered, strongly confined structures [16]. However, due to the taper of the waveguide, the confining capacity of the SCWT becomes weakened and more high-order modes leak out of the waveguide, which result in larger approximation errors and imperfect combining (imaging) performance.

To further study the effects quantitatively, we simulated the combining of two fiber lasers in SCWT with different taper angles. The parameters are listed in Table 1 and the results are shown in Fig. 5. It is shown that the distortion is negligible when $θ = 0$ and perfect combining performance can be achieved. The distortion of the field at the output of the SCWT increases with the taper angle, and the power in the main lobe is reduced to 76% when $θ = 2.8 °$ ( $\tan θ = 0.05$ ). This is due to the fact that more high-order modes leak out with larger taper angle, which results in deterioration of the combining (imaging) performance. To ensure more than 90% of the energy contained in the main lobe, the taper angle must be set to less than 1.4° ( $\tan θ < 0.025$ ). Within the main lobe area, the phase distribution for the tapered waveguide is nearly plane and the phase fluctuation is less than $λ / 10$ when $θ < 1.4 °$ . From Figs. 4 and 5, one sees that our analytical formulas agree well with the FDBPM simulation results when $θ < 1.4 °$ and the deviation becomes larger as the taper angle increases. Our analytical model is accurate for small taper angles, i.e., $θ < 1.4 °$ , and can be used to design the SCWT and the fiber laser array in waveguide-based CBC systems, which is useful in practical application.

Fig. 5. Field distributions at the output of the SCTW for different taper angles. (a), (c) Amplitude distribution and (b), (d) Phase distribution.

Download Full Size | PDF

Table 1. Parameters of the Tapered Waveguide

View Table | View all tables in this article

Table 2. Parameters for Setting of $3 \times 3$ Fiber Laser Array

View Table | View all tables in this article

A. Combination Based on SCTW

Based on our analytical formulas and the FDBPM, the combining of $3 \times 3$ fiber laser array (shown in Fig. 6) in the 3D square-cross tapered waveguide is simulated with different $θ$ . The waist diameter of the laser beam is 5 μm. The relative phase and positions in the input $x o y$ plane are shown in Table 2, and the other parameters are taken to be the same as those in Fig. 4. The parameters of the waveguide are listed in Table 3. The fiber lasers can be effectively combined provided the relative phase and incident position indicated in Eq. (13) are fulfilled. From Fig. 7, we can conclude that, with larger taper angle, the output beam is more compressed and more energy is diverted into the side lobes.

Fig. 6. $3 \times 3$ fiber laser array.

Download Full Size | PDF

Table 3. Parameters of the Square-Cross Tapered Waveguide

View Table | View all tables in this article

Fig. 7. Transverse intensity distribution of $3 \times 3$ fiber lasers.

Download Full Size | PDF

B. Combination Based on Nontapered Waveguides and Tapered Waveguides

The combining of $3 \times 3$ fiber laser arrays based on tapered waveguides is compared with that based on nontapered waveguides in this section. The output size of the waveguides is taken to be the same. Beam waist radius of each fiber laser is 2.5 μm with power $P = 0.1 a.u.$ . If a nontapered waveguide is adopted to combine the array laser beams, we can determine that $L_{3} = 540 μm$ . Figure 8 presents the output intensity distribution of $3 \times 3$ fiber lasers based on different waveguides. It reveals that, with the same size output port, the peak energy of the combined beam is higher for CBC based on a tapered waveguide, which is due to the taper of the waveguide.

Fig. 8. Transverse intensity distribution of $3 \times 3$ fiber lasers.

Download Full Size | PDF

The results of two kinds of CBC systems are listed in Table 4 for comparison. It indicates that, compared with planar waveguides, larger input size can be achieved using tapered waveguides, which is meaningful for coupling the fiber laser array into the waveguide but is achieved at the cost of waveguide length. In the practical application, we should choose waveguide based on a practical situation.

Table 4. Results of Two Kinds of CBC Systems

View Table | View all tables in this article

C. Performance Analysis

For applications such as energy transmission, the beam quality of the combined beam is a key parameter to our concern. The well-known parameter $M^{2}$ is taken as the characteristic parameter to analyze the combination performance of the system which, for the 2D case, is defined as [30]

M^{2} = \sqrt{\frac{w^{2} (z) - w_{0}^{2}}{{(λ / π w_{0})}^{2} (z - z_{0})}}

with

w (z) = 2 σ (z),

Defining a fill factor $t = 2 N w_{0} / W_{1}$ to characterize the filling of the input part of the waveguide, we can get the results in Fig. 9. For the sake of simplicity, we consider the 1D array. The parameters of the waveguide are taken to be the same as in Fig. 4, except for the variable parameters. From Fig. 9(a), we can conclude that there exists an optimal $t$ , which corresponds to a best input beam waist. It proves that near-diffraction limited beam quality ( $M^{2} \leq 1.41$ ) is obtained. We define a tolerance to ensure that the value of the $M^{2}$ is less than 1.43, and the tolerance on fill factor is $0.45 < t < 0.67$ . It is shown, by comparing the results in Figs. 4(f) and 9(b), that the phase distortion of the output field can be mitigated by optimization of the fill factor $t$ . The far-field intensity distribution [Fig. 8(c)] is obtained by propagating the near field in free space using the Fast Fourier Transformation method [31,32]. It reveals that although the near field has side lobe, one single beam without side lobes can be achieved in far field by optimization of the fill factor.

Fig. 9. Optimal designation of the system and the results. (a) $M^{2}$ as a function of $t$ , (b) near-field distribution for the optimum $t$ , and (c) far-field intensity distribution for different $t$ .

Download Full Size | PDF

Figure 10 gives optimum $t$ for different numbers of lasers. One can see that, with input beam number $N$ increasing, the optimal $t$ increases and the beam quality of the combined beam decreases: the optimum $M^{2}$ is 1.84 for $N = 3$ and 2.26 for $N = 4$ . This is due to the fact that more high-order modes are generated in the waveguide (tapered or nontapered) with larger $N$ , which leak out and result in deterioration of the beam quality.

Fig. 10. $M^{2}$ as a function of $t$ for different beam array.

Download Full Size | PDF

Figure 11 reveals the influence of the tapered angle. The beam waist width is 7.5 μm. It shows that, with larger taper angle, the beam quality of the combined beam degrades and the length of the SCWT becomes short. Although the best beam quality ( $M^{2} = 1.02$ ) can be obtained when nontapered waveguides are used, the length of the waveguide ( $L_{2} = 1823 μm$ ) is one half longer than that of the tapered waveguide ( $M^{2} = 1.43$ , $L_{2} = 1216 μm$ ), which may prevent a compact system configuration. One can conclude that the use of tapered waveguide can result in a compact configuration with moderate degrading of the beam quality. It also shows that the tolerance on the taper angle is $θ < 0.8 °$ .

Fig. 11. Dependence on taper angle. $N = 2$ , $W_{1} = 50 μm$ .

Download Full Size | PDF

D. Discussion

Although the manufacturing of 3D tapered waveguides is not easy, the designation of 2D tapered waveguides is realizable. Figure 12 shows the combination based on a 2D tapered waveguide. The 2D waveguide is approximated by a rectangular cross with the short side (in the $y$ direction) being $W_{1} = 5 μm$ and the long side (in the $x$ direction) being 50 μm. The waveguide is shown in Fig. 12(a). The long side is far longer than the short side and the longer side of the waveguide has no taper while the short side tapers gradually. The parameters are taken as $\tan θ = 0.0103$ . The waist width of two Gaussian beams is 0.5 μm. Figure 12(c) shows the two fiber laser beams combined into a single laser beam. In practice, the aforementioned approach can be realized by using two metal plates with dielectric coating or glass plates with grazing incidence, high-reflection coating.

Fig. 12. Combining based on 2D SCWT. (a) Schematic diagram of a 2D tapered waveguide, (b) intensity distribution of input laser beams, and (c) intensity distribution of output laser beams.

Download Full Size | PDF

The actual order of the waveguide is hundreds of micrometers in the experiments. Hence, modeling the real waveguide requires at least one order of magnitude more computational time and memory than the modeled waveguide, which is beyond the computer resources in our lab for the time being. We are working to improve and optimize the computer programs to study the waveguide with the size in the experiment in the future. The conclusions obtained in the simulations are instructive for real experiments.

4. CONCLUSIONS

In this paper, we have presented a rigorous analysis for the self-imaging properties of fiber lasers in a SCWT and derived some analytical formulas for the positions, amplitudes, and phases at the output of the waveguide. Numerical simulation has been carried out to verify our analytical formulas. The results agree well with the predication of our analytical formulas when the taper angle is less than 1.4° and the phase distortion resulting from the taper of the waveguide is less than $λ / 10$ . The model is useful in designing and performance analyzing in the following ways:

(1) It can predict the combing length (self-imaging length) of the tapered waveguide with error less than 1%, which is meaningful for saving time during designing the waveguide.
(2) It can provide the phase relationships between the input laser beams and incident positions of the laser beams on the input part of the waveguide, which is useful in simulation of waveguide-based CBC based on beam propagation method.
(3) It provides a direct insight of the CBC of laser beams based on tapered waveguide and proves this novel method theoretically.

Based on our analytical formulas and the FDBPM, we simulated the CBC of fiber lasers based on waveguide, and the effects of input beam number, aperture fill factor, and taper angle on the combination performance have been investigated. The engineering realization of such a configuration is discussed. It shows that this novel approach can achieve a high beam quality output beam ( $M^{2} \leq 1.41$ ). If the output size is fixed, with larger taper angle, larger input size can be obtained; if the input size is fixed, shorter length of the waveguide and more compact system configuration can be obtained with larger taper angle. There exists an optimum fill factor $t$ . With input beam number increase, the optimum fill factor will increase and the beam quality will degrade. Although the near field distribution has a multilobe configuration, one single beam can be achieved in the far field. The tolerance of the system on the fill factor and the taper angle is studied, which is $0.45 < t < 0.67$ and $θ < 0.8 °$ . The CBC based on waveguide can be realized in 2D waveguide effectively. In practical application, we should choose and design waveguides based on the practical situation and take beam quality, system compaction, and cost into consideration.

ACKNOWLEDGMENTS

This project was supported by Innovation Foundation for Excellent Graduates in National University of Defense Technology under grant B120704.

REFERENCES

1. D. J. Richardson, J. Nilsson, and W. A. Clarkson, “High power fiber laser: current status and future perspectives,” J. Opt. Soc. Am. B 27, B63–B92 (2010). [CrossRef]

2. J. R. Leger, J. Nilsson, J. P. Huignard, A. P. Napartovich, T. M. Shay, and A. Shirakawa, “Laser beam combining and fiber laser systems,” IEEE J. Sel. Top. Quantum Electron. 15, 237–239 (2009). [CrossRef]

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

4. T. M. Shay, J. T. Baker, A. D. Sanchez, C. A. Robin, C. L. Vergien, C. Zerinque, D. Gallant, C. A. Lu, B. Pulford, T. J. Bronder, and A. Lucero, “High-power phase locking of a fiber amplifier array,” Proc. SPIE 7195, 71951M (2009). [CrossRef]

5. B. He, Q. H. Lou, J. Zhou, J. Dong, Y. Wei, D. Xue, Y. Qi, Z. Su, L. Li, and F. Zhang, “High power coherent beam combination from two fiber lasers,” Opt. Express 14, 2721–2726 (2006). [CrossRef]

6. J. Wang, K. Duan, and Y. Wang, “Experimental study of coherent beam combining of two fiber lasers,” Acta Phys. Sin. 57, 5627–5631 (2008).

7. P. Zhou, Z. J. Liu, X. L. Wang, Y. X. Ma, H. T. Ma, X. J. Xu, and S. F. Guo, “Coherent beam combining of fiber amplifiers using stochastic parallel gradient descent algorithm and its application,” IEEE J. Sel. Top. Quantum Electron. 15, 248–256 (2009). [CrossRef]

8. L. Liu, M. A. Vorontsov, E. Polnau, T. Weyrauch, and L. A. Beresnev, “Adaptive phase-locked fiber array with wavefront phase tip-tilt compensation using piezoelectric fiber positioners,” Proc. SPIE 6708, 6708K (2007).

9. J. E. Kansky, C. X. Yu, D. V. Murphy, S. E. J. Shaw, R. C. Lawrence, and C. Higgs, “Beam control of a 2D polarization maintaining fiber optic phased array with high-fiber count,” Proc. SPIE 6306, 63060G (2006). [CrossRef]

10. P. Zhou, Y. Ma, X. Wang, H. Ma, J. Wang, X. Xu, and Z. Liu, “Coherent beam combination of a hexagonal distributed high power fiber amplifier array,” Appl. Opt. 48, 6537–6540 (2009). [CrossRef]

11. E. C. Cheung, J. G. Ho, G. D. Goodno, R. R. Rice, J. Rothenberg, P. Thielen, M. Weber, and M. Wickham, “Diffractive-optics-based beam combination of a phase-locked fiber laser array,” Opt. Lett. 33, 354–356 (2008). [CrossRef]

12. R. Uberna, A. Bratcher, and B. G. Tiemann, “Coherent polarization beam combination,” IEEE J. Quantum Electron. 46, 1191–1196 (2010). [CrossRef]

13. S. E. Christensen and O. Koski, “2-Dimensional waveguide coherent beam combiner,” in Advanced Solid-State Photonics, OSA Technical Digest (CD) (Optical Society of America, 2007), paper WC1.

14. R. Uberna, A. Bratcher, T. G. Alley, A. D. Sanchez, A. S. Flores, and B. Pulford, “Coherent combination of high power fiber amplifiers in a two-dimensional re-imaging waveguide,” Opt. Express 18, 13547–13553 (2010). [CrossRef]

15. W. S. C. Chang, Fundamentals of Guided-Wave Optoelectronic Devices (Cambridge University, 2010).

16. R. Ulrich and T. Kamiya, “Resolution of self-images in planar optical waveguides,” J. Opt. Soc. Am. 68, 583–592 (1978). [CrossRef]

17. M. Bachmann, P. A. Besse, and H. Melchior, “General self-imaging properties in N×N multimode interference couplers including phase relations,” Appl. Opt. 33, 3905–3911 (1994). [CrossRef]

18. S. He, X. Ao, and V. Romanov, “General properties of N×M self-images in a strongly confined rectangular waveguide,” Appl. Opt. 42, 4855–4859 (2003). [CrossRef]

19. H. J. Baker, J. R. Lee, and D. R. Hall, “Self-imaging and high-beam-quality operation in multi-mode planar waveguide optical amplifiers,” Opt. Express 10, 297–302 (2002). [CrossRef]

20. R. Scarmozzino, A. Gopinath, R. Pregla, and S. Helfert, “Numerical technique for modeling guided-wave photonic devices,” IEEE J. Sel. Top. Quantum Electron. 6, 150–162 (2000). [CrossRef]

21. R. Scarmozzino and R. M. Osgood Jr., “Comparison of finite-difference and Fourier-transform solutions of the parabolic wave equation with emphasis on integrated-optics applications,” J. Opt. Soc. Am. A 8, 724–731 (1991). [CrossRef]

22. Y. Chung and N. Dagli, “An assessment of finite difference beam propagation method,” IEEE J. Quantum Electron. 26, 1335–1339 (1990). [CrossRef]

23. Y. Chung and N. Dagli, “Modeling of guided-wave optical components with efficient finite-difference beam propagation methods,” in Antennas and Propagation Society International Symposium (IEEE, 1992), pp. 248–251.

24. C. Vassallo, Optical Waveguide Concepts (Elsevier, 1991).

25. C. Vassallo, “Interest of improved three-point formulas for finite-difference modeling of optical devices,” J. Opt. Soc. Am. A 14, 3273–3284 (1997). [CrossRef]

26. K. Kawano and T. Kitoh, Introduction to Optical Waveguide Analysis (Wiley, 2001).

27. C. Chen, Foundations for Guided-Wave Optics (Wiley, 2007).

28. G. R. Hadley, “Transparent boundary condition for the beam propagation method,” Opt. Lett. 16, 624–626 (1991). [CrossRef]

29. G. R. Hadley, “Transparent boundary condition for the beam propagation method,” IEEE J. Quantum Electron. 28, 363–370 (1992). [CrossRef]

30. A. E. Siegman, “How to (maybe) measure laser beam quality,” in Diode Pumped Solid State Lasers: Applications and Issues, M. Dowley, ed., Vol. 17 of OSA Trends in Optics and Photonics (Optical Society of America, 1998), paper MQ1.

31. T. Poon and T. Kim, Engineering Optics with MATLAB (World Scientific, 2006).

32. J. D. Schmidt, Numerical Simulation of Optical Wave Propagation with Examples in MATLAB (SPIE, 2010).

$W_{1} (μm)$	50	50	50	50	50	50
$\tan θ$	0.05	0.025	0.0069	0.0034	0.0017	0
Analytical results $L_{2} (μm)$	394.6	649.5	1219.8	1469.2	1631.2	1833.5
Simulation results $L_{2} (μm)$	395	649.8	1220.2	1469.6	1631.6	1833.8

	1	2	3	4	5	6	7	8	9
Relative phase	0	$- π / 3$	0	$- π / 3$	$- 2 π / 3$	$- π / 3$	0	$- π / 3$	0
Position	$- W_{1} / 3$ , $W_{1} / 3$	0, $W_{1} / 3$	$W_{1} / 3$ , $W_{1} / 3$	$- W_{1} / 3$ , 0	0, 0	$W_{1} / 3$ ,0	$- W_{1} / 3$ , $- W_{1} / 3$	0, $- W_{1} / 3$	$W_{1} / 3$ , $- W_{1} / 3$

	1	2	3	4
$W_{1} (μm)$	50	50	50	50
$\tan θ$	0.0103	0.0121	0.0138	0.0155
Analytical results $L_{3} (μm)$	814.5	769.6	731.6	697.1
Simulation results $L_{3} (μm)$	815	770.2	732.5	698.5

	$W_{1} (μm)$	$W_{2} (μm)$	$L_{3} (μm)$	$I_{peak}$ (a.u.)	$η$ (%) ( $h = 2 μm$ )
Tapered waveguides	50	33	814	1.34	85.19
Nontapered waveguides	33	33	540	0.89	86.07

$W_{1} (μm)$	50	50	50	50	50	50
$\tan θ$	0.05	0.025	0.0069	0.0034	0.0017	0
Analytical results $L_{2} (μm)$	394.6	649.5	1219.8	1469.2	1631.2	1833.5
Simulation results $L_{2} (μm)$	395	649.8	1220.2	1469.6	1631.6	1833.8

Coherent beam combination of fiber lasers with a strongly confined tapered self-imaging waveguide: theoretical modeling and simulation

Abstract

1. INTRODUCTION

2. THEORETICAL MODEL

A. System Introduction

B. Imaging of a Fiber Laser through a SCTW

C. Imaging of a Tapered Waveguide

3. NUMERICAL RESULTS

A. Combination Based on SCTW

B. Combination Based on Nontapered Waveguides and Tapered Waveguides

C. Performance Analysis

D. Discussion

4. CONCLUSIONS

ACKNOWLEDGMENTS

REFERENCES

Cited By

Figures (12)

Tables (4)

Equations (37)

Photonics Research