1. Introduction
When light interacts with a medium, acoustic waves are generated on which the light is reflected. This reflected optical wave, in the following referred to as the Brillouin wave, will have a frequency in the order of 10 GHz lower than the incident light for silica based fibers [
1
G. P. Agrawal, Nonlinear Fiber Optics , Third Edition, (Academic Press, San Diego, 2001).
]. This phenomenon, known as stimulated Brillouin scattering (SBS), is the nonlinear optical phenomenon with the lowest threshold for narrow-band light. This low threshold has been discussed in relation to limitations of continuous wave (CW) high power fiber lasers and has thus been an active research subject since the achievement of low loss optical fibers [
2
E. P. Ippen, “Low-power quasi-cw Raman oscillator,” Appl. Phys. Lett.
16, 303–305 (1970).
[CrossRef]
–
11
V. I. Kovalev and R. G. Harrison, “Continuous wave stimulated Brillouin scattering in optical fibers: New results and applications for high power lasers,” in Topical Problems of Nonlinear Wave Physics
, Proc. SPIE
5975, 59750L-1-13 (2006).
]. The current investigation is motivated by the recent use of Raman amplification of a Brillouin fiber laser [
12
J. W. Lou, F. K. Fatemi, and M. Currie, “Brillouin fibre laser enhanced by Raman amplification,” Elect. Lett.
40, 1044–1046 (2004).
[CrossRef]
] and experimental observation and theoretical modeling of passive Q-switching in Raman fiber lasers due to Brillouin self-pulsation [
13
G. Ravet, A.A. Fotiadi, M. Blondel, and P. Mégret, “Passive Q-switching in all-fibre Raman laser with distributed Rayleigh feedback,” Elect. Lett.
40, 528–529 (2004).
[CrossRef]
,
14
A. A. Fotiadi, P. Mégret, and M. Blondel, “Dynamics of a self-Q-switched fiber laser with a Rayleigh-stimulated Brillouin scattering ring mirror,” Opt. Lett.
29, 1078–1080 (2004).
[CrossRef]
[PubMed]
].
When propagating CW light in optical fibers the backward propagating Brillouin wave is amplified due to SBS, where the input CW signal act as a pump. As the power of the input CW signal is increased, the Brillouin gain is increased and consequently the backward propagating power of the Brillouin wave increases to a level where it starts to deplete the input signal. If the input signal power is increased even further, the depletion causes the signal power to drop so that the system eventually enters a mode of self-pulsation. If feedback is introduced to the system ie. through reflections at the fiber ends, the input power threshold of this oscillating behavior is reduced. The threshold condition for SBS and stimulated Raman scattering individually in the absence of feedback was first given by R. G. Smith [
5
R. G. Smith, “Optical Power Handling Capacity of Low Loss Optical Fibers as Deterined by Stimulated Raman and Brillouin Scattering,” App. Opt.
11, 2489–2494 (1972).
[CrossRef]
]. In the case of low reflectivity the SBS threshold was given by A. L. Gaeta and R. W. Boyd [
8
A. L. Gaeta and R. W. Boyd, “Stimulated Brillouin scattering in the presence of external feedback,” J. Nonlinear Opt. Phys. Mater.
1, 581–594 (1992).
[CrossRef]
] including the material response time and for arbitrary external feedback by M. Dämmig
et al. [
9
M. Dämmig, G. Zinner, F. Mitschke, and H. Welling, “Stimulated Brillouin scattering in fibers with and without external feedback,” Phys. Rev. A
48, 3301–3309 (1993).
[CrossRef]
[PubMed]
] not including the material response time.
In the following a threshold condition for oscillations due to SBS of a Raman amplified signal is derived for arbitrary external feedback under assumptions on the depletion due to the SBS. The derivation follow the method used by R. G. Smith [
5
R. G. Smith, “Optical Power Handling Capacity of Low Loss Optical Fibers as Deterined by Stimulated Raman and Brillouin Scattering,” App. Opt.
11, 2489–2494 (1972).
[CrossRef]
] and M. Dämmig
et al. [
9
M. Dämmig, G. Zinner, F. Mitschke, and H. Welling, “Stimulated Brillouin scattering in fibers with and without external feedback,” Phys. Rev. A
48, 3301–3309 (1993).
[CrossRef]
[PubMed]
]. The analytic results are compared to numerical investigations of a discrete and a distributed Raman amplification schemes. From these specific examples the result of depletion on the threshold condition is shown, a point which to the knowledge of the authors has not been addressed previously.
2. Theory
The present work focus on the effects of loss and Raman amplification and in the following an analytical equation for the power threshold, at which self-pulsation occur, is derived. The setup considered here, shown in
Fig. 1, consist of a single mode, inversion symmetric fiber of length
L placed between two reflectors having power reflectivities,
R
+ in the signal input end and
R
- in the signal output end; both reflects the light back into the fiber and thus creates a Fabry-Pérot cavity. The cavity can thus arrive either on purpose or unintentionally i.e. due to fiber splices. It is noted that the length between the two reflectors considered here is in the order of hundreds of meters to kilometers. Consequently the longitudinal modespacing is small (~500 kHz for a 200 m cavity) compared to the linewidth of the Brillouin gain spectrum (in the order of 20–100 MHz for silica based fibers [
1
G. P. Agrawal, Nonlinear Fiber Optics , Third Edition, (Academic Press, San Diego, 2001).
]). Thus the investigation does not involve mode locking due to cavity resonances, but Brillouin self-pulsation due to wave dynamics. For a polarization maintaining single mode fiber, neglecting four wave mixing e.g. due to phase mismatch, neglecting higher order Brillouin scattering, and assuming instantaneous material response the governing equations can be expressed in power as [
1
G. P. Agrawal, Nonlinear Fiber Optics , Third Edition, (Academic Press, San Diego, 2001).
]
where P
±
S=P
±S (z,t), P
±
B=P
±
B (z,t) and P
±
R=P
±
R (z,t) are the powers of the signal, Brillouin and Raman waves, the +(-) superscripts indicates forward(backward) propagating waves, gi
is the Brillouin (i=B) or Raman (i=R) gain coefficient, z is the length, and t is the time. ωS
, ω
B and ωR
are the angular frequencies of the signal, Brillouin and Raman waves respectively. Their ratio takes into account the difference in energy between the Raman wave, the signal, and Brillouin waves. Since the change in frequency of the acoustic wave is usually very small compared to the optical frequency, the ratio between ω
S and ω
B is approximately 1. Furthermore α
S, α
B, and α
R are the attenuation coefficients and since ω
B≈ω
S then α
S≈α
B=α is used in the following. v
g,S,
v
g,B, and v
g,R are the group velocities, which can assume arbitrary values in the analytical approach. The pump for the Raman amplifier is assumed to have a line width broad enough to disregard Brillouin contributions. The assumption of an instantaneous material response is valid for long temporal pulses compared to the response of the material as in the present case.
The boundary conditions (BC)s for the Fabry-Pérot cavity are
where
P
in,±
i
is the signal (
i=S) or Raman (
i=R) input power,
ε
± is the Brillouin seed due to vacuum or thermal lattice fluctuations and
R
±
i
are the cavity end reflectivities for which the superscript +(-) refer to the reflection of the waves in the -(+) into the +(-) direction. The setup is illustrated in
Fig. 1.
Fig. 1. A fiber cavity of length L with power reflectivity R
+
i
in the left end and R
-
i
in the right end creates a cavity. AC Winput signal in this produce a Brillouin wave in the opposite direction. Both waves are amplified by a Raman amplifier. The signal output power, given by P
out
S (t)=(1-R
-
S)P
+
S (L,t), oscillates in case of high signal input power levels.
A way of defining the threshold condition, as it is done in Refs. [
5
R. G. Smith, “Optical Power Handling Capacity of Low Loss Optical Fibers as Deterined by Stimulated Raman and Brillouin Scattering,” App. Opt.
11, 2489–2494 (1972).
[CrossRef]
] and [
9
M. Dämmig, G. Zinner, F. Mitschke, and H. Welling, “Stimulated Brillouin scattering in fibers with and without external feedback,” Phys. Rev. A
48, 3301–3309 (1993).
[CrossRef]
[PubMed]
], is when the input signal power at which the intracavity SBS power in the signal input end in steady state reaches a fraction,
r, of the intracavity signal power at the input end, i.e.
r signifies the limit at which SBS depletion of the signal is nearly negligible. It will shortly be related to an experimentally obtainable parameter involving the low reflectivity limit of the threshold. The self-pulsations occur when the Brillouin wave depletes the signal. The power threshold, Pthr, is then given by the BC’s from insertion of
P
+
S (0)=
R
+
P
-
S (0)+
P
thr in
Eq. 3. It is noted that the threshold power defined here is the experimentally obtainable input signal power thus different from the steady-state intra cavity signal power at the signal input end,
P
+
S (0), used in Ref. [
9
M. Dämmig, G. Zinner, F. Mitschke, and H. Welling, “Stimulated Brillouin scattering in fibers with and without external feedback,” Phys. Rev. A
48, 3301–3309 (1993).
[CrossRef]
[PubMed]
]. The value
r=1 was used by R. G. Smith in Ref. [
5
R. G. Smith, “Optical Power Handling Capacity of Low Loss Optical Fibers as Deterined by Stimulated Raman and Brillouin Scattering,” App. Opt.
11, 2489–2494 (1972).
[CrossRef]
] yielding the estimate for no Raman amplification nor external feedback of
I
thr≈21
α=(
g
B
A
eff), with
I
thr being the threshold intensity and
A
eff the effective area of the signal.
The threshold is defined as the point immediately before the Brillouin wave significantly influence the signal. This implies that steady state is assumed for all the waves and the signal is assumed not to be effected by the Brillouin scattering. Furthermore the Raman pump is assumed undepleted. This gives, see Appendix A,
where G
±, P̃±
i
and J
± have been introduced for convenience and will be defined in the following. The effective Raman gain is defined as
similar to the Raman gain [
1
G. P. Agrawal, Nonlinear Fiber Optics , Third Edition, (Academic Press, San Diego, 2001).
] except for the forward and backward Raman input powers being replaced with the effective forward and backward Raman input powers
with R̃±
R=R
±
R exp(-α
R
L) the attenuated reflectivies. For long fiber lengths the attenuated reflectivies go to zero and the effective Raman gain equal the Raman gain.
In Eqs. (4) the effective signal and Brillouin input powers are defined as
with the effective reflectivities ρ
±
i
=R
±
i
G
+(L). Furthermore the effective position of the forward and backward Brillouin waves are defined as
where analytical expressions for J
±(z) are given in Appendix B. With this definition, J
+(L) and J
-(L) corresponds to the effective Brillouin fiber lengths of the forward and backward propagating waves respectively. If no Raman gain is introduced these reduce to the conventional effective length
When the Raman gain is introduced J
+(L) differ from J
-(L) since the forward and backward propagating signal in general experience different gain.
With Eqs. (4) the threshold power in the absence of feedback is given, using
Eq. (3), as
where the index R=0 indicates that all the reflectivities are zero.
An implicit equation for the cavity threshold can now be found by taking the ratio between
Eqs. (3) and
(10) giving
where the threshold signal input power is
and the ratio between the Brillouin inputs is
Eq. (11) is the SBS oscillation threshold condition for Raman amplified Brillouin fiber cavities. This implicit equation is the main result of this work. The arbitrary parameter
r has been replaced by the experimentally obtainable
P
thr
R=0, i.e. the threshold in the absence of external feedback. The result,
Eq. (11), reduces to the experimentally verified threshold obtained in Ref. [
9
M. Dämmig, G. Zinner, F. Mitschke, and H. Welling, “Stimulated Brillouin scattering in fibers with and without external feedback,” Phys. Rev. A
48, 3301–3309 (1993).
[CrossRef]
[PubMed]
] if no Raman gain is included and thus also reduces to the resulting equations of Ref. [
8
A. L. Gaeta and R. W. Boyd, “Stimulated Brillouin scattering in the presence of external feedback,” J. Nonlinear Opt. Phys. Mater.
1, 581–594 (1992).
[CrossRef]
] for equal and high reflectivity in both the fiber ends.
Similar results have been achieved for erbium doped amplifiers without taking the feedback into account, but including the temperature dependence on the Brillouin coefficient [
11
V. I. Kovalev and R. G. Harrison, “Continuous wave stimulated Brillouin scattering in optical fibers: New results and applications for high power lasers,” in Topical Problems of Nonlinear Wave Physics
, Proc. SPIE
5975, 59750L-1-13 (2006).
] [
15
A. N. Pilipetskii and V. V. Shkunov, “Calculation of the threshold and of the efficiency of conversion by stimulated scattering in an amplifying medium,” Sov. J. Quantum Elect.
15, 284–286 (1985).
[CrossRef]
].
3. Validation
In the following the signal power threshold is calculated using
Eq. (11) and compared to numerical results. This investigation begin by using realistic parameters for which the oscillation threshold of a discrete and a distributed backward Raman fiber amplification schemes are found. Then the effects of depletion on the threshold is examined and at last the effect of the different parameters of
Eq. (11) are studied. The validation is performed in order to show the applications of the derived threshold condition and point out its limitations. For the investigation the attenuations of all waves are identical, the reflectivities of all waves are identical and the reflections at both fiber ends are equal so that a round-trip reflectivity,
R
+
i
R
-
i
=
R
2, is defined.
The numerical investigation was conducted as described in Ref. [
16
M. E. V. Pedersen, J. R. Ott, and K. Rottwitt, “Self-Pulsation in Raman Fiber Amplifiers,” presented at the Eleventh International Conference on Transparent Optical Networks, Island of São Miguel, Azores, Portugal, 28 June-2 July 2009.
]. In the analytical approach the experimentally obtainable value of
r was chosen so as to match the numerical results. As mentioned in the previous section the value of
r describes the signal threshold power as the reflectivities go to zero.
For the amplification schemes the input power of the Raman amplifiers has been chosen such that they model an optical communication scheme. First a discrete amplifier case with a gain of 20 dB over a fiber of 200 m, e.g to counterbalance a 20 dB loss in a subsequent 100 km fiber. Secondly a distributed amplifier case with a 20 dB gain counterbalancing a simultaneous 20 dB loss in a 100 km fiber. The Raman input powers used to achieve the mentioned gains are calculated in absence of external feedback [
17
K. Rottwitt and A. J. Stentz, “Raman Amplification in Lightwave Communication Systems,” in Optical Fiber Telecommunications
, Chap. 5 , (Academic Press, San Diego, Calif., 2002).
[CrossRef]
]. The coefficients of the two different amplifier schemes are shown in
Table 1. The coefficients correspond to standard silica based optical fibers with different amounts of germanium dopants. If other fibers are used the assumption of instantaneous material response might not be valid.
Table 1. The used values of the amplifier schemes.
|
L
|
g
B
|
g
R
|
α
|
P
in,+
R
|
| Discrete Amp. | 200 m | 0.8 (Wm)-1
| 5 (Wkm)-1
| 0.8 dB=km | 4.7 W |
| Distributed Amp. | 100 km | 0.14 (Wm)-1
| 0.7 (Wkm)-1
| 0.2 dB=km | 0.3 W |
The numerical and analytical results are given in
Figs. 2 and
3 where the oscillation power threshold vs. the round trip reflectivity in dB, defined as 10log
10(
R
2), is shown. The solid line corresponds to the scheme in which Raman amplification is used to counterbalance fiber losses, the dashed line show the threshold without the amplification and the dash-dotted line show the threshold for lossless fibers. The thin lines with errorbars are the numerically obtained results.
When the transmitted signal become unstable, i.e. pulsed, both the Brillouin scattered signal as well as the transmitted input signal exhibit steady state oscillations. If both fiber ends are cleaved at a 90° angle to the direction of the fiber, the reflection at each end equals ≈-15 dB, corresponding to a round trip reflectivity of -30 dB however, typical components and splices have significant lower reflections.
The power threshold for the discrete amplifier scheme is shown in
Fig. 2 versus different round-trip reflectivity. The figure contain three pairs of curves one with a 20 dB gain, one without the gain and one without gain nor loss where a stable output, i.e. not oscillating, occur for input powers to the left of the curves. The thick lines corresponds to the derived threshold while the thin lines with error bars corresponds to numerically found thresholds. The numerical and analytical found thresholds coincide well for high reflectivities, but deviate slightly just before reaching the low reflectivity limit. This deviation is further investigated shortly. The addition of gain is found to reduce the stability region as expected since it increase the signal and Brillouin waves and thus effectively increases the Brillouin gain coefficient. This may also be seen from the derived equation,
Eq. (11), as the addition of Raman amplification increase the effective Brillouin fiber length
J
±(
L).
By comparing the curves for propagation with and without loss it is seen that the stability region increase as attenuation occur. This feature is more pronounced for the distributed amplifier scheme,
Fig. 3, as the fiber is longer and thus more lossy. The increase of the stability region can be understood intuitively as losses will reduce oscillations. This is also seen from
Eq. (11) as increasing the loss decrease the effective Brillouin fiber length which might be most easily seen in the case of no amplification where
J
±(
L)=[1-exp(-
αL)]/
α which goes to zero as
α becomes large.
Another very noticeable feature of
Fig. 3 is that the addition of losses increase the difference between the analytical and numerical curves. The difference for low reflectivities is due to the validity of neglecting depletion of the signal and the Raman waves. The depletion reduces the power in the signal and the Raman waves thus reducing the power transfer to the Brillouin wave thereby increasing the needed signal input power for oscillations. For high round-trip reflectivity the cavity reflectivity is higher than the attenuation keeping the powers at a high level making depletion negligible. Adding the Raman amplification is seen to increase the region in which the approximation of undepletion is valid as the gain counterbalance the losses. In
Fig. 4 the numerical stability calculation of the Raman amplified propagation shown in
Fig. 3 is performed for the three cases of i) depletion included, the thin solid, ii) without depletion on the Raman wave, the thin dashed line, and iii) without depletion on the Raman wave nor the signal wave, thin dash dotted line, all with error bars. In addition the derived analytical solution also shown in
Fig. 4, is shown as the thick solid line. The thin dash dotted line follow the analytical solution within the numerical resolution as it should. Furthermore, it is seen that the dashed line lie between the line with full depletion and the derived analytical solution indicating that the difference is caused by the validity of the assumption of undepleted signal and Raman pump. The assumption of negligible depletion is thus only valid for short fibers, low losses or high reflectivities. The derived
Eq. (11) thus always underestimate the threshold compared to the actual threshold.
Fig. 2. Oscillation power threshold for the discrete Raman amplifier vs. different round trip reflectivity, defined as 10log10(R
2). The blue solid lines corresponds to the scheme with a Raman gain of 20 dB, the green dashed lines show the threshold without the amplification and the red dash-dotted lines show the threshold for lossless fibers. The thick lines are the derived results and the thin lines with error bars are the numerically obtained results. The regions below the curves corresponds to CW output.
Fig. 3. Oscillation power threshold for the distributed Raman amplifier vs. different round trip reflectivity, defined as 10log10(R
2). The blue solid lines corresponds to the scheme in which Raman amplification is used to counterbalance fiber losses as described in the text, the green dashed lines show the threshold without the amplification and the red dash-dotted lines show the threshold for lossless fibers. The thick lines are the derived results and the thin lines with error bars are the numerically obtained results. The regions below the curves corresponds to CW output.
The fact that
Eq. (3) might seem contradictory to the experimental confirmation of the threshold without amplification, derived in Ref. [
9
M. Dämmig, G. Zinner, F. Mitschke, and H. Welling, “Stimulated Brillouin scattering in fibers with and without external feedback,” Phys. Rev. A
48, 3301–3309 (1993).
[CrossRef]
[PubMed]
], which
Eq. (11) reduce to in the limit of no Raman amplification. The reader should though not worry as the experiment of Ref. [
9
M. Dämmig, G. Zinner, F. Mitschke, and H. Welling, “Stimulated Brillouin scattering in fibers with and without external feedback,” Phys. Rev. A
48, 3301–3309 (1993).
[CrossRef]
[PubMed]
] was performed on a short fiber of less than 200 m and did not consider reflections below ~-85 dB.
Fig. 4. The threshold for the distributed Raman amplified scheme, shown in
Fig. 3, which is numerically calculated in the three cases of i) propagation including depletion, the thin blue solid line, ii) propagation without depletion on the Raman wave, the thin green dashed line, and iii) propagation without depletion on the Raman wave nor the signal wave, thin red dash dotted line, all with error bars. The thick black solid line show the derived threshold.
At last the threshold of the discrete Raman amplifier scheme with parameters in
Table 1 is investigated by changing the parameters individually. In
Fig. 5 seven different cases are studied. The thick lines show the derived
Eq. (11) while the thin marked lines corresponds to the numerical calculations. The blue, green and red colored lines show the threshold for backward, forward and bidirectional Raman amplifier schemes respectively, all with a gain of 10 dB in absence of external feedback. The three schemes yield exactly the same threshold in the analytical model. In the numerical calculations the backward amplification scheme is slightly more stable than the forward amplification scheme and the differences become more pronounced at the low reflectivity limit. This is because the signal wave for the backward amplification schemes in case of no external feedback has lower path average power than the forward amplification scheme and thus decrease the energy transfer to the Brillouin wave compared to the forward amplification scheme. Simultaneously the lower path average power increase the effects of depletion thereby making the analytical threshold differ slightly more from the numerical results for the backward amplification scheme compared to the forward and bidirectional schemes. The stability of the bidirectional scheme lie between the two others as it is a combination of these. As the feedback increases the backward and forward schemes approach the bidirectional scheme as the reflections in effect cause bidirectional amplification.
The light-blue line show that the threshold for the backward amplification scheme with only 5 dB gain has a higher threshold compared to the blue line with the same scheme, but 10 dB gain. This is expected as decreasing the gain increase the stability and is also apparent from the effective Brillouin cavity length,
J
±(
L), which decrease as the Raman input decrease. The opposite case of increasing the gain to 20 dB was shown in
Fig. 2 yielding a decrease of the stability region. The yellow and black lines show the threshold for no amplification and including losses and excluding losses respectively. These curves are also shown in
Fig. 2. The pink line show the backward amplification scheme with 10 dB gain, but a 1 km cavity. Increasing the cavity length increase the effective Brillouin cavity length and thus decrease the power threshold.
It is be noted that the analytical threshold yield qualitatively if not quantitatively the same results as the numerical calculations. Calculation times are though in the order of seconds for the implicit analytical equation where the numerical calculations take in the order of hours even days.
Fig. 5. Comparison of forward, backward and bidirectional amplification and an investigation of the effect of propagation length, loss and amount of amplification. The general values of L=200 m, g
B=0.8 (Wm)-1, g
R=5 (Wm)-1 and α=0.8 dB/km when not explicitly stated otherwise. The thick lines indicates the analytically found threshold while the thin lines with error bars show the corresponding numerical calculations. Color codes are as follows: Blue: Backward amplification, 10 dB gain. Green: Forward amplification, 10 dB gain. Red: Bidirectional amplification, 10 dB gain. Light-blue: Backward amplification, 5 dB gain. Pink: Backward amplification, 10 dB gain, 1 km cavity. Yellow: No amplification. Black: No amplification nor loss.