2. Model
Our model begins with the scalar paraxial optical wave equation valid for weakly-guiding large mode area fibers given by
Eq. (1):
where
is the scalar field amplitude,
is a suitably chosen wave vector magnitude,
is the transverse Laplacian operator,
is the refractive index distribution, and
is the vacuum optical wave vector magnitude. If we now set
the solution of the eigenproblem given by
Eq. (2)yields the stationary mode field profile
and propagation constant
Additionally we impose the normalization condition on
given by
Eq. (3)where
represents the fiber cross section. We now look for a higher-order mode
satisfying
Eq. (4) such that
Substituting
and invoking
Eq. (2) modified to describe
and
reveals that
which constitutes a somewhat non-standard yet deliberate definition of
The field evolution along the fiber is then given by
Eq. (5)where we have introduced the complex field amplitude
Allowing slow variation of the field amplitudes with
we arrive at a general expression for the scalar field representing the propagation of two transverse modes given by
Eq. (6):
The coefficients
are normalized such that the signal optical intensity field is given by
in which case the power in the modes is given by
To study the effect of thermally induced index variations on the evolution of the coefficients we introduce them perturbatively such that
Furthermore, laser gain is introduced locally so as to increase the optical intensity according to the usual two-level rate equations. We neglect resonantly-induced refractive index changes [
9A. A. Fotiadi, O. L. Antipov, and P. M’egret, “Resonantly induced refractive index changes in Yb-doped fibers: the origin, properties and application for all-fiber coherent beam combining,” in Frontiers in Guided Wave Optics and Optoelectronics, B. Pal, ed. (Intec, 2010), 209–234.
] since we are focusing on operating regimes for which these changes are much smaller than the thermally-induced ones. Invoking
Eq. (1),
Eq. (2),
Eq. (4) and
Eq. (6) leads to the equation for the evolution of the field
where the higher order term in the index perturbation has been neglected. Substituting
Eq. (6) into
Eq. (7) yields
where we have further approximated
This approximation is valid for weakly guiding fibers with little index variation across the core. Multiplying
Eq. (8) by each mode profile
, integrating both sides of the equation over the fiber cross-section and invoking the orthogonality of the two modes yields a coupled set of differential equations for the amplitudes
where the expression for
is given by
Eq. (10)This model thus accounts for both intra and intermodal non-linear phase modulation and gain. The first term on the right hand side of
Eq. (9a)-
(9b) represents the change in optical length of the fiber as its temperature changes as well as the optical gain of the modes. The second term represents the resonant energy transfer between the two as well as inter-modal laser gain. The choice of
as the reference wavevector for the slowly varying envelope approximation may lead to a breakdown in the accuracy of the solutions of
Eq. (9a)-
(9b) in the presence of a large thermal gradient. To mitigate this possibility we apply a variable transformation given by
Eq. (11)to
Eq. (9a)-
(9b) leading to the equations for the transformed amplitudes
where
is chosen to minimize the local self-phase modulation along the fiber taking into account the thermal gradient. Alternatively we could have decreased the longitudinal grid spacing to accurately approximate the increasing strength of the phase fluctuations. This transformation effectively locally optimizes the reference wavevector throughout the fiber length in the presence of thermal gradients.
The deliberate choice to include only two modes in the description of the optical field limits the ability of the model to treat cases where the index profile changes appreciably over the length of the fiber. An arbitrary number of modes may be added to our model in the future to treat such cases. Including the LP
02 mode in particular would enable the model to treat fibers with a thermally-induced shrinking mode field diameter [
10F. Jansen, F. Stutzki, H.-J. Otto, T. Eidam, A. Liem, C. Jauregui, J. Limpert, and A. Tünnermann, “Thermally induced waveguide changes in active fibers,” Opt. Express 20(4), 3997–4008 (2012), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-20-4-3997. [CrossRef] [PubMed]
]. No power loss results from this limitation however because the resonant coupling terms in
Eq. (12a)-
(12b) conserve total power.
To solve E
Eq. (12a)-
(12b) both the temperature and gain distributions throughout the fiber must be known. The thermal distribution evolves according to the time dependent heat equation given by
Eq. (13)with appropriate boundary conditions. Quantum defect heating is assumed to be the sole contributor to the heat load
The gain within the doped region of the fiber is given by
where the population densities
in
Eq. (14) are determined by the equilibrium solution to the two level rate equations and the quantum defect heating is given by
Eq. (15)It is important to note that the population densities and therefore the gain, heating, and temperature vary throughout the fiber as well as in time. We have neglected laser dynamics in our treatment on the grounds that we are primarily considering optical pump and signal intensities far above their respective saturation intensities such that the populations are assumed to reach equilibrium on a time scale smaller than that describing the thermal variations within the fiber.
We now introduce the key approximation of our approach. We treat the spatially-varying temperature distribution as static during the time required for the signal to completely traverse the length of the amplifier. Accordingly the rate equations are solved for the equilibrium population distributions and
Eq. (12a)-
(12b) are integrated for constant thermo-optic and gain profiles. The resulting modal amplitudes along the fiber are used to calculate the heat load which is then used to update the temperature distribution. More precisely we use a Runge-Kutta method to solve
Eq. (12a)-
(12b) at a given time step to obtain the optical intensity distribution throughout the fiber assuming uniform pump intensity within the cladding cross section at each point along the length of the fiber. This is then used to calculate the heating distribution. We then use a Crank Nicolson scheme to solve
Eq. (13) for the temperature distribution at the next time step. All thermal and optical degrees of freedom are represented on a Cartesian grid. In the counter-pumped case we furthermore approximate the total pump absorption to be independent of the relative amplitudes between the two modes. This alleviates the requirement to iterate the optical solution multiple times at each time step to converge on the appropriate residual pump level to ensure the launched pump level is uniform.
3. Modeling results
We model a large mode area photonic crystal fiber (PCF) with a pump cladding diameter of 170 µm, an outer cladding diameter of 400 µm, and doping concentration of 3.5 × 10
25 m
−1 which results in a device length of approximately 1.6 meters. This type of fiber is comparable in mode area and device length to some studied previously [
4T. Eidam, C. Wirth, C. Jauregui, F. Stutzki, F. Jansen, H.-J. Otto, O. Schmidt, T. Schreiber, J. Limpert, and A. Tünnermann, “Experimental observations of the threshold-like onset of mode instabilities in high power fiber amplifiers,” Opt. Express 19(14), 13218–13224 (2011), http://www.opticsinfobase.org/abstract.cfm?URI=oe-19-14-13218. [CrossRef] [PubMed]
,
5,
7C. Jauregui, T. Eidam, H.-J. Otto, F. Stutzki, F. Jansen, J. Limpert, and A. Tünnermann, “Temperature-induced index gratings and their impact on mode instabilities in high-power fiber laser systems,” Opt. Express 20(1), 440–451 (2012), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-20-1-440. [CrossRef] [PubMed]
,
8]. The large mode area results in a relatively long beat length between the two transverse modes. The beat length is the primary length scale that drives the longitudinal grid resolution requirement. The modeled fiber has a 7-cell core design resulting in a doped area that is somewhat smaller than the overall core area. The fiber has a pitch of 25 µm and an air hole diameter of 5.2 µm resulting in a mode field area of 2773 µm
2 for the fundamental LP
01 mode and 2939 µm
2 for the LP
11 mode as shown in
Fig. 1
. The beat length between these two modes in the absence of thermally induced index change is 22.6 mm.
Fig. 1 Photonic crystal fiber design (a), LP01 mode field profile (b) and LP11 field profile (c).
Our model describes a wide range of amplifier configurations. In this initial work we consider only counter-pumped configurations with 47.5 Watts of input power in the LP
01 mode and 2.5 Watts of input power in the LP
11 mode both inputs being monochromatic. To save the need to time-step throughout the entire warm up period, which can be several seconds, we first solve the time-independent heat equation for the thermal distribution assuming all of the power is in the LP
01 mode before adding in the LP
11 power and allowing time to evolve. We have established that this method of initiating the simulation gives the most rapid approach to conditions conducive to equilibrium.
Figure 2
shows results for three different pump powers with convective cooling at the boundary at a uniform 200 W/m
2K. The time step used was 40 µs and the spatial grid spacing was 111 µm in the longitudinal direction and 5.5 µm in the transverse direction. These were varied to ensure the numerical convergence of representative simulations. These values resulted in a spatial numerical grid with 7.9 × 10
7 points.
Fig. 2 Propagation distance dependence of the powers of the LP
01 and LP
11 modes after a time of 60 ms (a, c, e) and the output power in each mode as a function of time (b, d, f) for pump powers of 380, 636, and 952 Watts for a convectively-cooled amplifier. The media show the evolution over time from 0 to 60 ms (a [
Media 1], c [
Media 2], e [
Media3]).
Figure 2(a)-
2(b) describe the behavior of the amplifier when pumped with 380 Watts. As shown in the embedded movie (
Media 1), during the initial transient behavior power transfers from the LP
01 to the LP
11 mode. This coupling requires the phase of the optical intensity to lag the phase of the temperature variations in a particular region of the fiber [
6]. This is made possible by the fact that the optical phase can change much faster than the temperature phase. The optical phase shift at a given point in the fiber relative to the phase at the fiber entrance is accumulated throughout the fiber length. Therefore the relative phase change, and hence the coupling behavior, is strongest at the output end of the amplifier. The direction of the power transfer depends on the size of the phase shift. For shifts less than π radians, power couples from the LP
11 to the LP
01 mode while for shifts greater than π radians, power couples in the reverse direction. Thus for a monotonically increasing optical phase shift, power transfers back and forth between the modes as the local phase passes multiples of π as the signals propagate down the amplifier. According to this picture, once the temperature distribution reaches equilibrium with the optical phase, no coupling will be possible and the LP
01 mode will dominate the output as shown in
Fig. 2(a)-
2(b).
An additional feature is rapid temporary coupling of power between the modes at a spatial frequency corresponding to the inter-modal beat length. This rapid temporary coupling occurs because the effect of the local optical intensity on the temperature is mediated through the heat equation and the laser rate equations rather directly through the intensity. This leads to a coupling constant that reverses sign within one period of the induced grating. Previously published work has shown that no such temporary coupling exists where the index changed is caused purely by the non-linear refractive index [
5]. Another important feature of this rapid coupling is that it apparently does not depend on a particular thermo-optic phase shift because it remains present even for amplifiers below the instability threshold that reach equilibrium as shown in
Fig. 2(a) and
Fig. 4(a)
.
Fig. 4 Propagation distance dependence of the powers of the LP
01 and LP
11 modes after a time of 60 ms (a, c, e) and the output power in each mode as a function of time (b, d, f) for pump powers of 500, 1060, and 1750 Watts for a conductively-cooled amplifier. The media show the evolution over time from 0 to 60 ms (a [
Media 4], c [
Media 5], e [
Media6]).
Figure 2(c)-
2(d) describe the amplifier with 636 Watts of pump. We observe that the beat-length oscillations increase in amplitude and that the LP
11 mode acquires the most power for a period of about 10 ms after the initial transients. After this the modal output powers are approximately equal out to the end of the run. We do not claim that this represents an equilibrium state, but it is clear that the dynamic behavior persists for a longer time than for the case with lower pumping power. This power level represents the approximate instability onset threshold for this fiber under these conditions. We do not attempt to identify precise threshold values due to the multiple very long-time simulations required to do so.
Figure 2(e)-
2(f) describe the amplifier with 952 Watts of pump. Very interestingly, the amplifier appears to reach equilibrium with most of the power in the LP
11 mode. We chose these power levels to demonstrate three different behavior regimes. To accurately determine the rapidity of the onset of the instability we would need to run the simulation for longer times. For all configurations the total amplifier output is steady in time in agreement with experimental observations [
4T. Eidam, C. Wirth, C. Jauregui, F. Stutzki, F. Jansen, H.-J. Otto, O. Schmidt, T. Schreiber, J. Limpert, and A. Tünnermann, “Experimental observations of the threshold-like onset of mode instabilities in high power fiber amplifiers,” Opt. Express 19(14), 13218–13224 (2011), http://www.opticsinfobase.org/abstract.cfm?URI=oe-19-14-13218. [CrossRef] [PubMed]
]. We did not study the convectively-cooled amplifier at pump powers higher than 952 Watts because the maximum temperature of the fiber in this last case approached 500 °C which is probably as high as and perhaps higher than such amplifiers should sustain. We have indeed observed instability in simulations of convectively cooled amplifiers at higher powers that exhibited intolerable fiber core temperatures.
Figure 3
shows the temperature profile at the end of the 60 ms simulation for the 952 Watt pump case. We observe that the variation in the profile due to the inter-modal beating is small compared to the overall temperature profile.
Fig. 3 Temperature profiles at the end of the 60 ms simulation for the convectively-cooled amplifier with 952 Watts pump within planes passing through the center of the fiber core. Temperature is in °C with respect to the temperature of the cooling medium.
More effective cooling enables higher power levels. We next discuss the behavior of amplifiers with conductive cooling that exhibit constant temperature on one edge of the fiber as occurs with fibers wrapped on a cooled mandrel for example. We treat the fiber edge lying at the location lying on the negative
x axis 200 µm from axis in
Fig. 1. This arrangement is consistent with wrapping the fiber on a cooled mandrel. The effect of different cooling symmetries on the amplifier beyond the scope of this paper.
Figure 4(a)-
4(b) describe the behavior of such an amplifier pumped at 500 Watts. One interesting observation is that the modal contents appear to approach constant finite values as shown in
Fig. 4(b). An examination of
Fig. 4(a) as well as
Fig. 2(a) reveals that this steady value depends on the phase of the beat-length coupling at equilibrium. This suggests that fine-tuning of the pump power may enable optimization of the beam quality by ensuring that a trough in this variation occurs for the LP
11 modal power. In this case it is important to note that the relative phase of the two modes at the amplifier input is constrained by the fact that the LP
11 content is likely caused by launching the seed slightly off center. This particular effect may not be observable, however, for fibers exhibiting effective loss discrimination between the modes.
Figure 4(c)-
4(d) describe the behavior of a conductively-cooled amplifier pumped at 1060 Watts. One of the interesting features of this simulation is the presence of an apparent temporary equilibrium with most of the power in the LP
11 mode at the amplifier output followed by a rapid transition to an oscillatory behavior with most of the power in the LP
01 mode. As the pump power is increased, these oscillations increase in amplitude until they cause a transition to dynamic full-depth oscillation between the modes at the output shown in
Fig. 4(e)-
4(f) which describe the amplifier when pumped with 1750 Watts.
We propose that longitudinal heat flow combined with the different heat loads caused by the different transverse modes explain the onset of these oscillations. Without longitudinal thermal diffusion, time dependent instability theoretically should not exist in the absence of external perturbation. This may be seen by realizing that in the absence of such heat flow, an equilibrium solution to
Eq. (12) may be obtained by calculating the equilibrium
x-y temperature distribution at each point along the propagation direction using only the optical fields, and hence the heat load, at that longitudinal position. Furthermore, it has been shown [
7C. Jauregui, T. Eidam, H.-J. Otto, F. Stutzki, F. Jansen, J. Limpert, and A. Tünnermann, “Temperature-induced index gratings and their impact on mode instabilities in high-power fiber laser systems,” Opt. Express 20(1), 440–451 (2012), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-20-1-440. [CrossRef] [PubMed]
] that there is significant difference between the temperature profiles of such amplifiers in the presence and absence of longitudinal heat flow. These two facts support the view that longitudinal heat flow is necessary for dynamic modal instability. Nevertheless, we observed temporal oscillations in the modal content of the amplifier with a broader frequency spread at the highest pumping power with longitudinal heat flow artificially turned off in the simulation. We believe this point requires further investigation.
The next link in our proposed instability feedback loop is the effect of modal content on the quantum defect heating load. We have observed that the LP
11 mode creates less heat per unit length than the LP
01 mode. This can be seen in the movies linked to
Fig. 2 and
Fig. 4 where the temperature is shown to be lowered slightly in the regions along the fiber where the modal content is dominated by LP
11 and elevated where the modal content is dominated by LP
01. The phase of the optical interference pattern at any point along the fiber is determined by the relative phases of the two modes which are determined by the initial phases and the accumulated thermal phase shift from the input end of the fiber up to that point. These observations suggest that the ability of changes in modal content to affect the optical phase shift and the ability of the thermo-optic phase shift to affect the modal content through resonant coupling together give rise to the observed instability. We further note that this instability mechanism does not require the higher-order mode to possess a different optical frequency than the fundamental as is the case for a previously proposed mechanism [
6].
We also observe that the amplifier exhibits two distinct regions along the length of the fiber. Within the first part of the fiber, the amplitudes are stable in time. Beyond a critical length at which the rapid temporary coupling becomes sufficiently strong, the modal contents oscillate in time. Our model predicts that the length over which oscillating behavior occurs within the fiber should increase with increasing pump power.
Figure 5(a)
-
5(c) show the spectra of the modal content of the amplifier output at three different power levels exhibiting three different types of behavior. At the lower power level the output is relatively stable for long times as shown in
Fig. 4(a) with a corresponding peaked spectrum shown in
Fig. 5(a). At an intermediate power level the modal content exhibits relatively regular oscillatory behavior for the longest simulated times. The corresponding spectral peaks shown in
Fig. 5(b) are relatively small due to the limitations in the simulation run time. The time series shown in
Fig. 4(f) representing the highest output power has no dominant spectral components but is rather described by a noisy peak with a tail extending to several kHz as shown in
Fig. 5(c) indicating chaotic oscillatory behavior. This behavior remained unchanged when the cooling point on the fiber was shifted to a point on the
y axis in
Fig. 1. The temperature profile in the conductively-cooled case shows all of the expected features including the overall background profile with local variations imposed by the modal interference pattern as shown in
Fig. 6
. The maximum temperature of 100 °C relative to the heat sink should be well tolerable.
Fig. 5 Logarithmic spectral intensity of the time variations in the simulated output power of the LP
01 mode for pump powers of 500 Watts (a), 1060 Watts (b), and 1750 Watts (c) obtained via the Fourier transform of the time series shown in
Fig. 4(b),
4(d), and
4(f).
Fig. 6 Temperature profiles at the end of the 60 ms simulation for the conductively-cooled amplifier with 1750 Watts pump within planes passing through the center of the fiber core. Temperature is in °C with respect to the temperature of the heat sink.
There are several other dependencies that may be checked with our model that are beyond the scope of this paper. These include fiber dimensions such as length, core diameter, pump cladding diameter, and outer cladding diameter. The dependence on the seeding configuration, both total power level and modal content may also be investigated. The model may also be used to study the differences between co-pumped and counter-pumped configurations as well as the effect of differential passive loss between the modes and tailored dopant profiles within the core. Rather than trying to study the dependencies of the instability on all of these factors at the same time it makes more sense to conduct this type of analysis in the context of overall evaluation of proposed amplifier configurations.
4. Experimental results
We experimentally investigated modal instabilities in an amplifier featuring an SBS-suppressing photonic crystal fiber featuring a segmented acoustically tailored (SAT) core. Details of the core design are provided in [
11]. The SAT fiber allowed us to investigate the modal instabilities using a single-frequency seed, as well as phase-modulated and broadband seed sources. Initially, the amplifier was seeded with the output of an intermediate amplifier chain seeded with a single-longitudinal mode non-planar ring oscillator (NPRO).
The PCF exhibited a nominal core diameter and mode field diameter of 39.5 and 30 µm respectively. Borosilicate glass stress rods within the cladding enabled single-polarization operation. The inner pump cladding and outer cladding had diameters of 329 µm and 580 µm respectively with a pumping numerical aperture of 0.6 enabled by an air-clad structure. This fiber is substantially different to the one analyzed in the previous section. The experimental fiber has significantly shorter beat length, longer device length and larger outer cladding diameter than the modeled fiber. Hence modeling the experimental fiber requires at least 1 order of magnitude more computational time and memory than the modeled fiber. We are working to improve and optimize the computer programs to study such fibers in the future. Thus in the present work we seek to qualitatively compare the amplifier we were able to model to the available experimental amplifier.
Figure 7
depicts the experimental setup. The NPRO provided 10 mW of continuous wave seed power with a nominal linewidth of a few kHz that was subsequently amplified using a four stage amplifier chain the final stage of which consisted of the counter-pumped PCF amplifier. For the single-frequency studies, the phase modulator shown in the
Fig. 7 was not used. Our measurements indicated that ~30 W of the 50 W of third-stage output power was coupled into the PCF core. We used fiber-coupled diode stacks to launch up to about 1 kW of ~976 nm pump into the 10 meter length of PCF. Investigation into the thermal dependence of the phenomena was prompted by a comparison this fiber design to commercially available fibers used in high power fiber amplifiers. The 976 nm pump absorption of this fiber is at least twice that of 25/400 step index or 40/400 photonic crystal fibers. As such, the heat deposition per unit length is much higher in the fiber discussed here relative to other fiber designs.
Fig. 7 Experimental setup of counter-pumped PCF amplifier. A three stage amplifier system was utilized to provide approximately 30 Watts of seed power. An image of the beam was captured using a commercial M2 beam analyzer. The phase modulator was not utilized in the single-frequency experiments.
To investigate the influence of active cooling on modal instability we employed two different cooling arrangements. At first, the fiber was suspended in the air with a coil diameter of 53 cm and allowed to cool through convection. Each end of the fiber was held in a water cooled chuck, which was mounted on a 5-axis translation stage. The beam was analyzed using a commercial M2 measurement device. At power level below 300 Watts, the M2 value was less than 1.3 and the near field image showed a Gaussian-like profile consistent with the fundamental fiber mode (LP01). As the signal power approached 380 W, there was a sudden jump in the M2 value which was measured at 2.2. We have made multiple attempts to improve mode-matching of the seed laser into the photonic crystal fiber core and modify the fiber coiling geometry to increase higher order mode filtering. These exercises had negligible impact on the modal instability threshold.
A comparison of M
2 as a function of signal output power was performed for a fiber amplifier under a second cooling condition. The fiber was conductively cooled with an aluminum ring, as shown in
Fig. 8
. Water-cooled plates were affixed to the aluminum ring to improve cooling of the fiber. Approximately 2.5 meters of fiber, near the seed input end, was taken off the ring and suspended in air.
Fig. 8 Conductively cooled fiber amplifier setup. The seed laser is input from the left, where the fiber is held in a water-cooled chuck. The output end of the fiber departs the aluminum ring on the right and is held in a water-cooled chuck.
Each end of the fiber was held in a water cooled chuck, which was mounted on a 5-axis translation stage. The location of each stage was fixed on the optical table. The section of fiber left in the air was required to allow the fiber to conform to the input and output chuck positions. Improving thermal management of the fiber amplifier by using the cold ring was shown to increase the modal instability threshold beyond 500 Watts. For conductive cooling, the M
2 value was ~1.4 at power levels approaching 500 W.
Figure 9
compares the two cases by showing the measured M
2 value versus the signal power. As shown, improving thermal management of the fiber amplifier can increase the modal instability threshold. We note that without the benefit of the SAT fiber, the investigation of the modal instability in the conductively-cooled single-frequency amplifier would not have been possible due to the onset of SBS. In fact, in an experiment that utilized a reference fiber with similar dimensions to the SAT fiber but with an acoustically homogeneous core, the SBS threshold was encountered at ~150 Watts. Further power scaling of the SAT fiber was hampered by rapid fluctuations in the near field output of the amplifier. Near field images of the optical intensity profile below and above these modal instabilities are shown in
Fig. 10
. For the SAT fiber amplifier described above, these instabilities occurred at ~530 Watts. To rule out linewidth and acoustic tailoring effects on the modal instabilities we assembled an ASE source with an approximate spectral width of 10 nm and used it to seed both the SAT fiber and the reference fiber (i.e. fiber without acoustic tailoring). Since the reference fiber amplifier was SBS limited at ~150 Watts, the ASE source allowed us to reach power levels >500 Watts with this type of fiber. There was little difference in the threshold of the modal instabilities between the SAT and reference fibers.
Fig. 9 Beam quality as a function of output power for the PCF amplifier seeded with a single-frequency source which is (blue) conductively cooled and (red) convectively cooled (suspended in air).
Fig. 10 Near field image of SAT PCF amplifier operating at 300 Watts output power (a). Near field image of SAT PCF amplifier operating at ~530 Watts showing evidence of multi-mode behavior (b).
Additional measurements confirmed that the linewidth of the seed source had little effect on the threshold of the modal instabilities in the SAT fiber. In one study, the linewidth of the ASE source was reduced to 0.5 nm using a notch filter. In another study, the output linewidth from the single-frequency NPRO was broadened by modulating at 1 GHz using the LiNbO3 modulator. The modal instability threshold for the conductively cooled SAT amplifier remained constant at 500 Watts within the few Watt uncertainty of the optical power meter regardless of the linewidth.
Finally, we conducted a study to capture the time-dependent behavior of the output signal. It is known that low-speed cameras such as the one used for the M
2 studies presented above are not fast enough to capture the time dynamics of the energy transfer between the LP
01 and LP
11 modes [
12F. Stutzki, H.-J. Otto, F. Jansen, C. Gaida, C. Jauregui, J. Limpert, and A. Tünnermann, “High-speed modal decomposition of mode instabilities in high-power fiber lasers,” Opt. Lett. 36(23), 4572–4574 (2011), http://www.opticsinfobase.org/ol/abstract.cfm?URI=ol-36-23-4572. [CrossRef] [PubMed]
]. To increase the time resolution we focused a sample of the output beam into a single-mode fiber and monitored the output using a fast photodiode detector. The output from the fiber is proportional to the power in the fundamental mode. The trace of the output power in time is shown in
Fig. 11
exhibiting oscillations with frequencies on the order of several kHz. This trace was obtained near the threshold of 500 Watts and indicates that the amplifier was operating in the intermediate instability regime.
Fig. 11 Time series output (a) and frequency spectrum (b) of the fiber-coupled fast photodetector measuring the output of the PCF amplifier showing periodic oscillatory behavior. The output power was 500 Watts.
We observed threshold-like transverse modal instability in both simulation and experimental results. In each case, active conductive cooling of the fiber increased the threshold. For the conductively-cooled amplifier cases the onset of the instability was manifest by oscillatory behavior in both the experimental and simulated cases. In the experimental case, the power output was not pushed very far past the onset of the oscillatory instability. In the simulated case, however, pushing the power higher resulted in chaotic oscillations characterized by a decaying frequency spectrum. Thus in the conductively cooled case we find three basic stability regimes characterized by stable operation, periodic oscillations, and chaotic oscillations in modal content. In the convectively cooled case it appears possible that the amplifier could reach equilibrium with most of the power exiting the fiber in the LP
11 mode. The different behavior in the convectively and conductively cooled cases is not surprising in light of the drastically different temperature profiles observed in the two cases as shown in
Fig. 3 and
Fig. 6. One important aspect of our simulations is that they assume no external perturbations to the fiber such as those that would be caused by mechanical vibrations and fluctuations in pump and seed powers. Such perturbations may prevent experimental observation of the simulated LP
11 equilibrium state.
5. Mitigation strategies
Our model suggests a few promising approaches for mitigating the onset of this instability. The first approach is to minimize the maximum longitudinal thermal gradient within the core thus reducing longitudinal heat flow. One method of doing this is to modify the fiber design to promote constant gain per unit length throughout the fiber. Gain is a function of doping concentration as well as pump and signal intensities therefore the doping concentration or pump cladding diameter could be modified along the fiber length to achieve this goal. The latter of these two approaches, while perhaps easier to implement, would have brightness conservation implications for the pump. Another method would be to create an opposite gradient on the heat sink. Another approach would be to design the fiber core such that the thermo-optic overlap of the two modes is the same. This would cause the local heating load to be independent of the modal content. These two general approaches are compatible with each other and may be employed together to achieve the best results.
On the other hand, our simulation results suggest that filtering strength requirements for mode filtering fibers [
3] are more stringent as a result of the observed rapid intermodal coupling. Even differential mode losses of several dB/m may not be able to completely suppress higher order modes over lengths of millimeters. For example in the results depicted in
Fig. 4(c)-
4(d), the power in the LP
11 mode rises from nearly zero to about 200 Watts within a few millimeters. Another known approach for suppressing higher order modes is to maximize the overlap of the gain and the fundamental mode while minimizing the overlap with higher-order modes [
2T. Eidam, S. Hädrich, F. Jansen, F. Stutzki, J. Rothhardt, H. Carstens, C. Jauregui, J. Limpert, and A. Tünnermann, “Preferential gain photonic-crystal fiber for mode stabilization at high average powers,” Opt. Express 19(9), 8656–8661 (2011), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-19-9-8656. [CrossRef] [PubMed]
]. This is most directly accomplished by concentrating the doped area in the center of the core. According to our model, this would increase the disparity between the modal thermal self-phase modulation coefficients leading to a more rapid optical phase shift for a given temperature change thus possibly enhancing the efficiency of the proposed instability mechanism.
Modal instability now seems to be the limiting factor for further power scaling of narrow-linewidth large mode area fibers. Just as there have been multiple methods developed for increasing the SBS threshold in these fibers, we anticipate methods will emerge for increasing the modal instability threshold with specially designed transverse and longitudinal thermal profiles assuming a role similar to that of acoustic engineering in SBS-suppressing fibers.