1. Introduction
In recent years, the development of coherent light sources at a wavelength of about 2.3–2.5 μm has been attracting increasing interest [
1,
2,
3,
4,
5,
6,
7]. One of the main reasons is that the earth’s atmosphere is highly transparent in this spectral region. Hazardous gases such as HF, CH
4, and CO have pronounced absorption lines near 2.3 µm, which opens up an opportunity for using 2.3 μm lasers for remote sensing. In addition, laser sources at this wavelength are in demand for LiDARs and other applications. Many real applications require high power (tens of watts and even higher), which is a challenge for modern laser systems, but works in this direction are actively underway based on solid-state lasers [
1,
2,
3] and fiber lasers and amplifiers [
4,
5,
6,
7]. Bulk lasers can provide higher output powers on average, but fiber lasers are generally characterized by high beam quality, alignment-free operation, and good heat dissipation. In addition, some fiber laser systems are capable of producing high average power too (~100 W and higher).
A very promising way to develop laser sources operating beyond 2.1 μm is based on using tellurite (tellurium dioxide-based) glass fibers [
8]. Silica glasses commonly used for near-IR fiber devices have large losses beyond 2.1 µm, so their use is very limited at longer wavelengths (although not excluded [
9]). Tellurite glasses are transparent in the range of ~0.5–5 µm, have excellent optical and physicochemical properties [
8,
10,
11,
12,
13] and fairly high damage thresholds (i.e., they can resist high peak intensities [
14]), can be doped with high concentrations of active ions [
15], and are produced with low background losses [
8]. Lasers at 2.3 µm based on tellurite fibers doped with trivalent thulium ions have been successfully demonstrated [
5,
6,
7], but their experimental powers still do not exceed several hundred mW [
7]. Another promising approach to generate light at 2.3 µm is to use stimulated Raman scattering (SRS) in undoped tellurite fibers pumped at 2 µm. This approach was theoretically explored in [
16], where a continuous wave (CW) Raman laser with a power of ~160 W at a wavelength of 2.35 µm based on a single-mode fiber with a pump power of 300 W at 2 μm was demonstrated numerically. However, single-mode fiber lasers at 2 μm with a power of 300 W are non-standard and rather complex devices, manufactured only in individual scientific laboratories. In addition, it remains an open question whether a single-mode tellurite fiber can resist such CW powers.
In this work, for obtaining high-power radiation at 2.3 μm in a Tm-doped tellurite multicore fiber amplifier, we propose to simultaneously use a unique combination of cascade lasing at two radiative transitions together with SRS processes. This is inherent specifically for zinc–tellurite (TZ, TeO
2–ZnO-based) glasses doped with trivalent thulium ions, which, in principle, makes it possible to obtain the efficiency of conversion of the pump at 793 nm to a signal at 2.3 µm above the Stokes limit (0.793 µm/2.3 µm = 34.5%). A simplified diagram of laser levels of trivalent thulium in a TZ glass matrix is shown in
Figure 1a [
15]. Previously, it was shown that the cascade scheme of dual-band operation at the
3H
4 →
3H
5 and
3F
4 →
3H
6 transitions with a pump wave at 793 nm can significantly increase the efficiency for a wave at 2.3 μm compared to operation only at the
3H
4 →
3H
5 transition, since it helps to overcome the bottleneck problem related to a long lifetime of
3F
4 [
5]. It is also known that Raman gain is very large for TZ glasses, with the maximum at about 750 cm
–1 (22.5 THz), which is ~60 times higher than the maximum for silica glass [
15]. Therefore, additional amplification of a signal at 2.3 μm due to the SRS of the wave at ~2 μm is also possible. The scheme of the SRS process is shown in
Figure 1b. In this case, for the Raman amplification at 2.3 μm, the wavelength of 1.96 μm displaced by 22.5 THz is most suitable. The Raman gain of the TZ glass in comparison with the Raman gain of the silica glass is plotted in
Figure 1c based on the data presented in [
17]. Note that in [
5], the parameters of the system were such that the SRS processes for them were negligible.
The concept of a master-oscillator power amplifier (MOPA) considered in this paper is often used for high-power laser systems [
18,
19]. At the same time, special fiber designs, e.g., multimode large mode area (LMA) fibers operating in quasi-single mode, are frequently developed for high-power amplifiers [
18,
20,
21].
It is also possible to use taper fibers with a thin single-mode core at the input and a thick core at the output [
22,
23]. With adiabatic signal amplification, the output beam quality can be very high, while the scattering into higher modes is rather small. One can also use multicore fibers with
N cores when each core is a single-mode one [
24,
25,
26]. In the case of evanescent-field coupling, there are
N supermodes in the system [
25,
27]. It was demonstrated that the use of an out-of-phase supermode of a multicore fiber with a symmetrical arrangement of axially symmetric cores, for which the spatial phases in neighboring cores differ by π, has a number of advantages over a “common” in-phase supermode, for which the phases are the same in all cores [
28,
29,
30,
31]. This issue was discussed in detail in [
28,
29,
30,
31], here we only note that the in-phase mode is unstable at high powers, while the out-of-phase mode is stable. Moreover, it has been predicted theoretically that the self-focusing limit for the beam power (for continuous medium) can be overcome using the out-of-phase mode, while this is impossible for the in-phase mode due to its instability [
28,
29]. As shown experimentally and numerically, with selective excitation of the out-of-phase supermode during amplification, the mode composition at the output does not deteriorate [
30]. Moreover, the out-of-phase mode overlaps better with the doped core, so it has a higher gain (this statement will be demonstrated below). Multicore fibers can be fabricated from various glasses, including tellurite ones. A rare-earth ion-doped multicore fiber was experimentally reported in [
32], and a passive seven-core tellurite fiber was demonstrated, for example, in [
33]. Although the manufacture of high-quality multicore tellurite fibers is still a challenge today, in principle, technologies can be improved to the required level.
Here, we propose and numerically investigate high-power amplifiers at a wavelength of 2.3 μm based on a special double-clad TZ fiber with
N = 10 step-index cores doped with trivalent thulium ions symmetrically arranged along the ring. A multicore fiber geometry leads to a decrease in peak intensities of laser beams relative to a single-mode fiber geometry and helps to prevent potential fiber damage for high-power laser amplifiers. It will be demonstrated that the use of the out-of-phase supermode is more efficient than the in-phase supermode. We propose to use a low-power seed signal at 2.3 μm and an additional low-power seed signal at ~2 μm with a standard 300 W diode pump at 793 nm. Two seeding low-power waves at 2.3 μm and 1.96 μm (with a power of about 1 mW) can be generated simultaneously in laser elements based on various matrices doped with trivalent thulium ions [
5,
6,
7]. Note that at present, relatively cheap commercially available laser diodes provide powers up to several hundred watts at a wavelength of 793 nm. In this case, the efficiency of converting the pump power to a wave at 2.3 μm, as will be shown below, can reach 50% due to the combination of laser and SRS processes. The scheme of a cascade laser amplifier, with one wave used for additional amplification of another due to SRS, has not been studied before to the best of our knowledge.
2. Materials and Methods
We considered a special design of a double-clad TZ glass multicore fiber shown in
Figure 2a. Ten step-index cores doped with trivalent thulium ions are arranged in a ring. The diameter of each core is
d = 7 µm, and the distance between the centers of neighboring cores is
L = 1.2 ×
d = 8.4 µm. The refractive index of the first undoped cladding is chosen to be 0.1 less than the refractive index of the core glass (
Figure 2b), which can be achieved by a slight change in the content of zinc oxide in TZ glass compositions [
5]. The parameters of each core are chosen so that it is single mode at signal wavelengths. A diode pump at 793 nm is launched into the first undoped cladding, which has a diameter of
D = 200 μm. This diameter
D coincides with the typical diameter of the fiber core, which outputs a diode pump with a power of hundreds of watts. In this case, the numerical aperture between undoped cladding and cladding 2 of the TZ multicore fiber should correspond to the numerical aperture of the fiber that outputs the pump (typically NA = 0.22). The multicore fiber also has the second cladding (its diameter is not important in the calculations).
We found supermode structures for the special multicore TZ fiber using finite element modeling, as in our previous work [
34].
Figure 3 shows examples of all ten supermodes calculated at a wavelength of 1.96 µm. Mode 1 (in-phase supermode) has the highest effective refractive index, and mode 10 (out-of-phase supermode) has the lowest effective refractive index. The pair of modes 2 and 3 is degenerate in the effective refractive index. The pairs of modes 4, 5 as well as 6, 7, and 8, 9 are also degenerate. The in-phase and out-of-phase supermodes have the most symmetrical structures. For them, |
Ex| distributions are the same for each core. For other modes, |
Ex| distributions are not the same in different cores; moreover, in some cores, the fields are close to zero or small. Therefore, for modes 2–9, the gain coefficients provided by trivalent thulium ions are obviously lower due to the smaller overlap with the doped region than for modes 1 and 10. Moreover, the preparation of seed signals for amplifiers in modes 2–9 is more difficult than in modes 1 and 10. In-phase and out-of-phase supermodes can be prepared using, for example, a special beam shaper (systems based on spatial light modulators are commonly used for such purposes) [
30,
35]. Further, we consider only in-phase and out-of-phase modes.
We calculated effective mode areas for in-phase and out-of-phase supermodes as functions of wavelength. The results are plotted in
Figure 4. The mode spatial structures at 1.96 µm and 2.3 µm which are of primary interest are shown in the insets in
Figure 4.
Effective supermode areas are calculated as
where
Sz is the
z-component of the Poynting vector.
A schematic diagram of the considered multicore fiber laser amplifier seeded by a dual-wavelength source at 1.96 µm and 2.3 µm with a high-power diode pump at 793 nm is shown in
Figure 5. To excite a required supermode, a beam shaper, for example, based on a spatial light modulator can be used [
30].
To simulate a dual-wavelength laser amplifier with allowance for the SRS process between waves at 1.96 and 2.3 µm, we generalized the system of rate equations and the power evolution equations presented in [
5,
15]. The SRS was considered in a standard way [
27]. We analyzed CW waves, but this model can be valid for describing pulsed amplifiers with signal durations noticeable longer than the total lifetime of level
3F
4. In this case, all powers refer to peak values. This approach may be justified when average powers should be reduced. The rate equations for the population densities
n1,
n2,
n3, and
n4 (normalized to the concentration of trivalent thulium ions
NTm) in the steady state are [
5,
15]:
where
t is time,
τj is the total lifetime of level
j,
τ4R and
τ4NR are the radiative and non-radiative lifetimes of level 4,
Wjk are the stimulated rates from level
j to level
k,
β4j is the branching ratio from level 4 to level
j, and
KCR is the cross-relaxation rate (see the scheme in
Figure 1a) [
5,
15]. Stimulated absorption and emission rates for a pump wave are:
where
c is the speed of light,
h is the Planck’s constant,
σ14 and
σ41 are the absorption and emission cross sections at a pump wavelength
λp,
Pp is the wave power,
Adoped is the area of doped cores, and
Γp is the overlap integral of the pump intensity distribution with the doped regions estimated as
Γp =
Adoped/(π
D2/4). Stimulated emission and absorption rates for signal waves at
λ1 = 1.96 µm and
λ2 = 2.3 µm with powers
P1 and
P2, respectively, are:
where
Γ1,2 are the overlap integrals of the signal waves with the doped regions.
The SRS-modified equations for the power evolution along the
z-axis are:
where
α1,
α2, and
αp are the linear losses at wavelengths
λ1,
λ2, and
λp, respectively;
gR12 is the Raman gain coefficient for the TZ fiber calculated as
Here,
gR(TZ@1µm) is the maximum Raman gain for bulk TZ glass for the pump wave at 1 µm [
17,
27]. For most of the simulations, we set optical losses of 0.3 dB/km for laser waves, which is easily achievable for modern single-core tellurite fibers and, in principle, can be achieved for multicore tellurite fibers. At the very end of the section Results, we will also analyze the impact of higher losses (up to 5 dB/m) that can be caused by light scattering at the core-cladding interfaces, bending losses, and so on.
The system of Equations (9)–(11), with allowance for the rate Equations (2)–(5) in the steady state, was solved numerically by the classical Runge–Kutta method (with adaptive step size). The used parameters are listed in
Table 1. We neglected the spontaneous Raman scattering to the anti-Stokes wave and the second-order Raman scattering to the Stokes wave. At room temperature, the process of scattering to the anti-Stokes component is weak. Moreover, its expected wavelength of 1.7 µm is strongly absorbed by Tm ions. The 2nd order Stokes Raman wave at the expected wavelength of about 2.8 µm is strongly absorbed by hydroxyl groups in zinc–tellurite glass.
3. Results
Let us first revisit the main idea of the paper: in high-power dual-wavelength fiber amplifiers based on a special TZ fiber doped with trivalent thulium ions, the interaction between waves at 1.96 μm and 2.3 μm due to SRS leads to a significant increase in the efficiency of power conversion from the diode pump at 793 nm to the wave at 2.3 μm. Moreover, the conversion efficiency is higher for the out-of-phase supermode than for the in-phase one because the effective area of the out-of-phase supermode is smaller than that of the in-phase one (for example, the effective mode area is 430 µm
2 for out-of-phase and 620 µm
2 for in-phase supermodes at 2.3 µm); the out-of-phase supermode is better localized near the cores. This is due to the fact that the topological structures of the field distributions of in-phase and out-of-phase modes are different. There are lines of zero field amplitude between the cores for the out-of-phase mode (due to equal amplitudes and phase difference by π at these lines), so most of its power is contained in the cores. On the contrary, there is significant field amplitude in the regions between the cores for the in-phase mode, as clearly seen in
Figure 4. The overlap integrals with doped cores for the out-of-phase supermode are higher (see
Table 1); hence, the laser gain coefficients are higher too. In addition, the Raman gain is also higher for the out-of-phase supermode (see
Table 1).
Figure 6 (left column) shows the results of calculations performed for the out-of-phase supermode, and
Figure 6 (right column) for the in-phase one. Hereinafter, the seed signal powers at 1.96 µm and at 2.3 mW are 10 mW and 1 mW, respectively. We checked that the exact values of low-power seed signals do not affect the output powers of 100-W class signals. The pump power was 200 W, and the concentration of active ions was
NTm = 4 × 10
20 cm
−3. We purposely switched on and off the SRS in the calculations to reveal its contribution. The results obtained for the complete model are plotted as solid curves and for the model without SRS as dashed curves.
Figure 6a,b show the evolution of the wave power at 2.3 μm, and
Figure 6c,d at 1.96 μm; whereas
Figure 6e,f show the evolution of the pump wave power.
At short distances from the fiber input end, there is cascading laser amplification of both signal waves, and the contribution of the SRS is small. The SRS effect becomes significant when the power in the wave at 1.96 µm is sufficiently large and almost reaches its maximum. This occurs approximately when most of the pump power has been absorbed (vertical line through
Figure 6a,c,e). After that, the wave at 1.96 µm ceases to be amplified. In the complete model with SRS, the power is further efficiently transferred from the wave at 1.96 µm to the wave at 2.3 µm. In the model without SRS, the powers in both signal waves reach a maximum at close points
z and then decay due to losses; there is no further interaction between them. The wave at 1.96 µm decays faster since it is partially absorbed from the ground state (at the
3H
6 →
3F
4 transition). The wave at 2.3 µm is not absorbed at the
3H
6 →
3F
4 transition. At practically zero pump powers, the levels above
3F
4 are not populated; therefore, the wave at 2.3 µm cannot be absorbed from excited states. Note that the power distribution in the pump wave is practically the same both in the model with and without SRS. In the case with allowance for SRS, the maximum power at 2.3 μm is 86 W in the out-of-phase supermode and 66 W in the in-phase one. For the out-of-phase supermode, the diode pump power conversion efficiency is 43%, exceeding the Stokes limit of 34.5%. Further study and optimization allowed us to find parameters of the system for which it is possible to achieve higher values of powers and conversion efficiencies.
Next, we studied the dependence of the power in two signal waves on two variables: the diode pump power and the fiber length for out-of-phase (
Figure 7, left-hand panels) and in-phase (
Figure 7, right-hand panels) supermodes with and without SRS (at
NTm = 4 × 10
20 cm
−3). The power curves at 2.3 µm are plotted in the top row (
Figure 7a–d) and at 1.96 µm in the bottom row (
Figure 7e–h).
It can be seen that at pump power < 100 W, the wave at 2.3 μm is weakly amplified, while the wave at 1.96 μm is significantly enhanced. This is explained by the cross-relaxation effect, which is rather strong at the chosen concentration of active ions. Its action leads to a decrease in the population at the 3H4 level and a decrease in the laser gain at the 3H4 → 3H5 transition, as well as an increase in the population at the 3F4 level and a significant increase in the laser gain at the 3H4 → 3H6 transition. Thus, at a pump power of 100 W, the maximum wave power at 1.96 μm is 60 W (the conversion efficiency reaches 60%) in the out-of-phase mode, which significantly exceeds its Stokes limit (0.793/1.96 = 40%). In this case, the red dashed-dotted curves in the lower panels demonstrate the optimal lengths of the wave at 1.96 μm depending on the pump power (i.e., the z value at which the maximum P1 is reached at a certain Pp). These dependences are rather weak and are determined by the pump absorption length, which changes weakly for the considered parameters. The blue-dotted curves in all plots show the optimal fiber lengths for maximizing wave power at 2.3 µm as a function of pump power (i.e., the z value at which P2 is maximized at a certain Pp).
In the case of switching off the SRS process, the optimal lengths are quite close for both waves since amplification occurs only due to cascade laser processes, which cease when the pump is absorbed. In the case of an allowance for SRS, the optimal lengths maximizing wave power at 2.3 μm strongly depend on the pump power. The higher the pump power, the shorter the optimal length. This has a very simple explanation: the higher the pump power, the higher the power in the wave at 1.96 µm, which means that the faster it amplifies the wave at 2.3 µm due to SRS. The optimal length for a wave at 2.3 µm is reached when the gain becomes equal to the optical loss. With the chosen TZ fiber parameters, this occurs for
P1 = 9.5 W for the out-of-phase mode and
P1 = 13.5 W for the in-phase mode (i.e., the blue-dotted curves in
Figure 7e,g show the level lines corresponding to these values). In addition, as
Pp increases, the population of the
3H
6 level increases, and the laser effects are also enhanced for the wave at 2.3 μm. At any pump power, the maximum power achieved in the wave at 2.3 μm is higher for the out-of-phase supermode than for the in-phase one.
Next, we made optimizations for the wave at 2.3 µm for various concentrations of active ions. The ratio of the power in this wave to the pump power (
P2/
Pp) as a function of two variables, pump power and fiber length, is presented in
Figure 8.
The following features are visible here: At low concentrations of trivalent thulium ions, amplification of the wave at 2.3 μm with a noticeable efficiency (>10%) begins at a pump power of about 10 W (
Figure 8a), while with increasing concentration, a higher pump power is required (
Figure 8b–d). At the highest considered value of
NTm, the cross-relaxation is so strong that the laser amplification of the wave at 2.3 μm is small, and notable values of its power are achieved mainly due to the SRS process at high powers and long fiber lengths (
Figure 8d). The optimum for
P2 is achieved at intermediate concentrations of
NTm =24 × 10
20 cm
−3 at high pump power (
Pp > 150 W) (
Figure 8b,c). In this case, both effects are important: laser amplification and SRS. The maximum conversion efficiency reaches 50%, which exceeds the Stokes limit by 15%.
Figure 8.
(a) Ratio of the power in the wave at 2.3 µm to the pump power (P2/Pp) as a function of two variables: pump power and fiber length calculated for NTm = 1 × 1020 cm−3 (a); NTm = 2 × 1020 cm−3 (b); NTm = 4 × 1020 cm−3 (c); and NTm = 8 × 1020 cm−3 (d). Dotted blue curves in all panels demonstrate the optimal fiber length maximizing P2 for certain pump power.
Figure 8.
(a) Ratio of the power in the wave at 2.3 µm to the pump power (P2/Pp) as a function of two variables: pump power and fiber length calculated for NTm = 1 × 1020 cm−3 (a); NTm = 2 × 1020 cm−3 (b); NTm = 4 × 1020 cm−3 (c); and NTm = 8 × 1020 cm−3 (d). Dotted blue curves in all panels demonstrate the optimal fiber length maximizing P2 for certain pump power.
Next, for clarity, we plotted the dependence of the wave power at 2.3 µm on the pump power for optimal concentrations:
NTm = 2 × 10
20 cm
−3 (
Figure 9a) and
NTm = 4 × 10
20 cm
−3 (
Figure 9b). Each curve in
Figure 9a,b corresponds to a certain length L = 1, 2, …, 7 m and, in fact, reproduces the data from
Figure 8b,c. It is clearly seen that at the most optimal lengths of 3–4 m and at high powers (>150 W), these dependences are close to linear functions.
Finally, we studied the effect of increased fiber losses on the quantitative characteristics of laser amplifiers. For all simulations presented above, we set optimistic losses of 0.4 dB/m at 793 nm and 0.3 dB/m at 1.96 and 2.3 µm. Here, we purposely increased losses at laser wavelengths to take into account additional imperfections (including bending losses). The simulated output powers at 2.3 µm vs. the fiber length for optimal concentration
NTm = 4 × 10
20 cm
−3 calculated for out-of-phase supermodes with SRS and without SRS are plotted in
Figure 10a,b, respectively, for different losses in the 0.3–5 dB/m range. It is not surprising that the higher the loss, the lower the output power and the shorter the optimal fiber length. More interesting is that the relative contribution of the Raman gain is much higher for relatively low losses. Comparing
Figure 10a,b, it can be seen that at losses up to 2 dB/m, the output power with SRS is significantly higher than without it, while at losses ≥3 dB/m, the Raman contribution becomes relatively small. Therefore, for the experimental implementation of the proposed amplifier, one should aim to manufacture a high-quality multicore tellurite fiber.
4. Discussion and Conclusions
To conclude, we have shown numerically that a trivalent-thulium-ion-doped TZ multicore fiber with coupled axially symmetrical cores arranged in a symmetrical ring structure is a promising solution for diode-pumped optical amplifiers of 100-W power class at 2.3 μm, which is required for important applications. We studied the use of dual-wavelength seed signals with a power of order 1 mW at 1.96 μm and 2.3 μm and a commercially available diode pump at 793 nm with a power of up to 300 W. We demonstrated the possibility of attaining up to 50% conversion efficiency from the pump power to the wave at 2.3 µm which exceeds the Stokes limit by 15% due to a combination of laser amplification and SRS between signal waves (for losses of 0.3 dB/m). The proposed scheme can provide significantly higher efficiency in comparison to trivalent-thulium-ion-doped fiber lasers and amplifiers that do not benefit from additional Raman energy transfer [
4,
5,
6,
7]. Note that the use of a fiber based on TZ glass is substantial, since its Raman shift of 22.5 THz coincides with the difference in the frequencies for two different radiative transitions of the trivalent thulium ion. For other glasses that are transparent in this range, the Raman shift is smaller (~17 THz for ZBLAN [
36], 10 THz for arsenic [
37], ~7 THz for arsenic triselenide [
37]), which is not suitable for implementing the concept of combining cascade laser amplification and SRS. We emphasize that the use of the out-of-phase supermode is important since its overlap integral with doped cores is higher compared to the commonly used in-phase supermode, whereas the out-of-phase supermode provides a higher total gain. In addition, it was previously shown that the out-of-phase supermode is stable at high powers, while the in-phase one is unstable [
28,
29,
30]. As was shown in [
30], small perturbations of the ideal transverse fiber structure led to a slight deformation of the out-of-phase mode, but do not lead to qualitative changes. Therefore, when deviating from the ideal geometry, only minor quantitative changes in the amplifier parameters are expected.
Note that the proposed design is especially good for high pump powers > 150 W. It provides high total power with an allowable load on each core. In this case, all powers refer to peak values. At lower target peak powers, a similar optimization can be performed using fewer cores, thereby reducing the effective mode area and increasing the Raman gain. One can also use a smaller cladding diameter to increase the overlap integral of the pump with doped cores. We considered CW waves, but the developed approach is valid for pulsed amplifiers with signal durations notably longer than the total lifetime of level 3F4 (3 ms).