Abstract

We investigate the Brillouin dynamic grating generation and detection process in polarization-maintaining fibers for the case of continuous wave operation both theoretically and experimentally. The four interacting light waves couple together through the material density variation due to stimulated Brillouin scattering. The four coupled equations describing this process are derived and solved analytically for two cases: moving fiber Bragg grating approximation and undepleted pump and probe waves approximation. We show that the conventional grating model oversimplifies the Brillouin dynamic grating generation and detection process, since it neglects the coupling between all the four waves, while the four-wave mixing model clearly demonstrates this coupling process and it is verified experimentally by measuring the reflection of the Brillouin dynamic grating. The trends of the theoretical calculation and experimental results agree well with each other confirming that the Brillouin dynamic grating generation and detection process is indeed a four-wave mixing process.

© 2011 OSA

1. Introduction

Since stimulated Brillouin scattering (SBS) was first observed in 1964 [1], it has remained an intensive research topic. Although SBS limits the channel power in optical fiber transmission systems since it could occur at a low input power in optical fibers, many applications could benefit from this nonlinear phenomenon, such as light storage [2], slow and fast light [3], Brillouin lasers and amplifiers [4], as well as distributed fiber sensors [5]. More recently, dynamic gratings generated in polarization-maintaining fibers (PMFs) via SBS, called Brillouin dynamic gratings (BDGs) [6], have attracted much attention. For example, it could be used to realize high-spatial resolution distributed sensing [710], to completely discriminate temperature and strain for Brillouin based sensors [1114], to achieve tunable optical delays [15], and to measure birefringence of PMFs in a distributed manner [16].

SBS occurs when two pump waves, whose frequency difference equals to Brillouin frequency, counter-propagate in an optical fiber. The two pump waves beat together giving rise to density variation associated with an acoustic wave through electrostriction effect [17]. The acoustic wave introduces a moving periodically modulated refractive index acting like a moving fiber Bragg grating (FBG) which diffracts more Stokes wave reinforcing the SBS process. The property of the moving grating could be measured or detected by monitoring the diffracted wave from a third probe wave which is used to illuminate the grating. In order to discriminate the probe and the pump waves, a direct way is to decouple the polarization state of the pump and probe beams, so that most of the experiments are carried out in PMFs. This moving FBG was given the name of BDG in Ref. [6] where separation of the grating generation and diffracting processes was achieved: two pump waves with frequency difference of a Brillouin frequency polarized along the same axis of the PMF generates the BDG, while a probe wave polarized along the orthogonal axis is diffracted strongly when phase matching condition is satisfied. The frequency difference between the diffracted and the probe waves also equals to the Brillouin frequency. So far, it is widely accepted that the fundamental nature of the BDG is similar to that of an FBG so that conventional FBG theory [18] could be directly applied. Recent experiments [19, 20] also tend to confirm this by rigorously measuring the intrinsic spectra of the BDG showing that the linewidth of the BDG appears inversely proportional to the BDG length. Note that BDG spectrum observed at the beginning with a broad bandwidth of ∼ 80 MHz [6] or ∼ 320 MHz [11] were contributed mainly due to the birefringence non-uniformity of PMFs as well as the relative frequency fluctuations between the multiple sources [14].

However, as discussed above, the BDG generation and detection process actually involves four waves interacting with each other through the material density variations; moreover, the diffracted wave has a Brillouin frequency difference with respect to the probe wave. Obviously, this would introduce a second SBS process between diffracted and probe waves, since these two waves have the same polarization state. More specifically, in most of the experiments [6,11,20], the diffracted wave would experience an SBS amplification since its frequency is downshifted by a Brillouin frequency from that of the probe beam whose propagation direction is the same as that of the BDG. Since the power level of the probe wave is relatively high in practice [11,19,20] in order to have strong diffracted wave for detection, this second SBS process is expected to play an important role. Physically speaking, this is a four-wave mixing (FWM) process that the longitudinal acoustic wave couples the four optical waves with different polarization states together. In this work, we will theoretically study this FWM process which is equivalent to the BDG generation and detection process in the case of continuous wave (CW) operation by deriving the four coupled light wave equations. These equations could be solved analytically for two cases: moving FBG model approximation and undepleted pump 1 (see Fig. 1) and probe waves approximation; showing that the FBG model clearly deviates from some of the reported experimental results [6]. In Section 3, experimental results will be provided to compare with our theoretical calculation confirming that the BDG generation and detection process is indeed an FWM process, and the moving FBG model oversimplifies this process. It is worth mentioning that even though this work is intent on studying the CW case, we believe that in the transit case (probe pulse follows one of the pump pulses) [7, 16, 21], when the power of the probe is high, the second SBS process is also expected to have influences on the reflectivity of the BDG characteristics.

 

Fig. 1 Diagram of the theoretical model. BDG generation and detection scheme: an FWM process.

Download Full Size | PPT Slide | PDF

2. Theoretical model and discussion

The diagram of the BDG generation and detection process is shown in Fig. 1. Two pump waves 1 and 2, whose frequencies are locked near the Brillouin frequency, are propagating oppositely inside a PMF with the same polarization state along slow () axis. The frequency of pump 1 wave is assumed to be higher than that of the pump 2 wave, so that the generated BDG through SBS is propagating along +z direction. A probe wave 3 polarized along fast (ŷ) axis has the same propagation direction as pump 1 wave; when phase-matching condition is satisfied, a diffracted wave 4 would be generated whose frequency is downshifted from the frequency of 3 by an amount of about a Brillouin frequency which is determined by the frequency difference between the two pump waves. This diffracted wave would experience a second SBS amplification by 3 along fast axis, since it has the same polarization state with the probe wave. Therefore, inside the PMF, the material density variation (acoustic wave) is contributed by two SBS processes which in turn couple the four optical waves together which could be considered as an FWM process.

In order to obtain a clear view of four waves interacting with each other in a PMF, we need to solve the Maxwell’s equations with the following material density equation [17]:

2ρ˜t2ΓA2ρ˜tvA22ρ˜=12ɛ0γe2E˜tot2,
where ρ̃ = ρρ0 is a change in the local density from its average value ρ0; ΓA is the damping coefficient; γe = ρ0(/)ρ=ρ0 is the electrostrictive constant; vA is the sound velocity inside the fiber; ɛ0 is the vacuum permittivity. The total electric field, tot, inside the PMF could be expressed as tot = x + y = (1 + 2) + ŷ(3 + 4), where
E˜x=12x^[F(x,y)A˜1(z,t)ei(k1zω1t)+F(x,y)A˜2(z,t)ei(k2zω2t)]+c.c.,
E˜y=12y^[F(x,y)A˜3(z,t)ei(k3zω3t)+F(x,y)A˜4(z,t)ei(k4zω4t)]+c.c.,
where Ãj(z,t) (j = 1 – 4) and F(x,y) are the slowly varying fields and dimensionless fundamental mode profile for the four interacting waves, respectively. Note that the mode profiles are assumed to be the same for all the four interacting waves for simplicity. ωj and kj are, respectively, frequencies and propagation constants of optical waves; c.c. represents complex conjugate.

As discussed above, the material density variation is contributed by two SBS processes, which could be expressed as ρ̃ = ρ̃1 + ρ̃2. When the frequency difference of the two pump lasers is locked at the Brillouin frequency, ΩB, under steady-state condition with slowly varying envelope approximation, the coupled equations for all the four optical waves could be expressed as (a more detailed derivation could be found in Appendix A),

A1z=κ(A1|A2|2+A2A3A4*eiΔkz),
A2z=κ(|A1|2A2+A1A3*A4eiΔkz),
A3z=κ(A3|A4|2+A1A2*A4eiΔkz),
A4z=κ(|A3|2A4+A1*A2A3eiΔkz),
where the coupling constant is
κ=8π3γe2cρ0λ3ΩBΓBAeffao,
with
Aeffao=[F2(x,y)F2(x,y)FA(x,y)]2FA2(x,y),
being the acousto-optic effective area [22, 23]; FA(x, y) is the dimensionless mode profile of the acoustic wave (only the fundamental acoustic mode is considered); angular brackets denote averaging over transverse cross section of the fiber; c is the velocity of light in vacuum; λ is the operation wavelength; ΓB is the Brillouin linewidth. Δk = k3 + k4k1k2 is the phase mismatch which is directly related to the frequency difference, Δν, of pump 1 and probe waves [6]. Note that Eq. (3) have been normalized as that |Aj|2 (j = 1 – 4) represents optical power.

Generally speaking, Eq. (3) needs to be solved numerically. However, with the approximation that pump 1 and probe waves are undepleted, these equations could be solved analytically which would give a clear physical insight for the FWM process. The analysis will be given in Sec. 2.2. Before that, we would like to discuss briefly the moving FBG model which is broadly accepted in all the reports so far, and show that such model oversimplifies the BDG generation and detection process.

2.1. Moving FBG model

So far, most reports assumed the fundamental nature of the BDG to be similar to that of a moving FBG with the experimental confirmation [19,20]. Under this approximation, the frequencies of the two pump waves are locked at Brillouin frequency, the depletion of pump 1 wave and the amplification of pump 2 wave are both ignored. SBS process generates an acoustic wave (density variation) which is considered as a moving FBG. The probe wave whose frequency satisfies the phase-matching condition is diffracted by this grating, and there is no SBS process between probe and diffracted waves. Physically speaking, the diffraction of the probe wave by this moving FBG is actually an acousto-optic process under these assumptions. Mathematically, it could be considered as that Eq. (3) reduce to

A3z=κA1A2*A4eiΔkz,
A4z=κA1*A2A3eiΔkz,
with A1 and A2 being constants. Equations (6a) and (6b) could then be solved analytically for the diffracted wave along with the boundary conditions, A3|z=0 = A3(0) which is a known quantity and A4|z=L = A4(L) = 0, where L is the length of the PMF [18],
A4(z)=2κexp(iΔkz/2)sinh[g(Lz)/2]A1*A2A3(0)gcosh(gL/2)iΔksinh(gL/2),
where g=4κ2|A1|2|A2|2(Δk)2. Under phase-matching condition, Δk = 0, the reflectivity of the BDG could be expressed as
RFBG=|A4(0)A3(0)|2=tanh2(κ|A1||A2|L).
Note that |A1| and |A2| are the square roots of the pump 1 and pump 2 powers, respectively. Figure 2 shows the properties of the functions of tanh2(x) and tanh2(x). It is clearly seen that the trends for tanh2(x) and tanh2(x) are very different. More specifically, when pump powers |A1|2 and |A2|2 are fixed, the reflectivity of the BDG varies with the length of the PMF following the trend showing as the solid curve. This has been confirmed by a recent experiment observing the intrinsic BDG spectrum [20]. On the other hand, when the length L and pump 2 power |A2|2 are kept at constant, the BDG reflectivity versus pump 1 power |A1|2 should follow the trend showing as the dashed curve. However, in the early experiment [6], an exponential growth of the BDG reflectivity with respect to pump 1 power has been observed, which clearly deviates from the result expressed by Eq. (8). Therefore, this approximated model condition oversimplifies the problem and needs to be modified by considering all the four coupling waves.

 

Fig. 2 Properties of the tanh2(x) and tanh2(x) functions.

Download Full Size | PPT Slide | PDF

2.2. FWM process with undepleted pump 1 and probe waves

In most of the BDG experiments, the pump 1 power is much higher than pump 2 power [6, 16, 20]. Moreover, in order to have strong diffracted signal from the BDG, probe power is also chosen to be on the same order of the pump 1 power. Therefore, a more reasonable approximation compared to moving FBG model is that pump 1 and probe power are undepleted, since the length of PMF used in the experiments is short, usually less than tens of meters. In addition, since the probe power is also high in most of experiments, the diffracted signal would also experience an SBS amplification. Under these assumptions, Eq. (3) could be reduced to

A2z=κ(|A1|2A2+A1A3*A4eiΔkz),
A4z=κ(|A3|2A4+A1*A2A3eiΔkz),
with A1 and A3 being constants. Equations (9a) and (9b) could be solved analytically along with the boundary conditions, A2|z=L = A2(L) which is a known quantity and A4|z=L = 0, as
A4(z)=2κA1*A2(L)A3sinh[Φ2+4Ψ(Lz)/2]Φ2+4ΨeΦ(Lz)/2+iΔkL,
where Φ = κ(|A1|2 +|A3|2) – iΔk and Ψ = |A3|2Δk. Under phase-matching condition, Δk = 0, the reflectivity of the BDG could be expressed as (a more detailed derivation could be found in Appendix B),
RFWM=|A4(0)A3|2=|A1|2|A2(L)|2[eκ(|A1|2+|A3|2)L1|A1|2+|A3|2]2.

We plot Eq. (11) with fixed length of the PMF, L, by varying pump or probe powers. The results are shown in Fig. 3. The parameters used to plot these curves are as follows: γe = 0.902, c = 3 × 108 m/s; ρ0 = 2210 kg/m3, Aeffao=80μm2, ΩB/2π = 10875 MHz, ΓB/2π = 30 MHz, λ = 1550 nm, L = 13 m. From Fig. 3(a), we could see that when the probe and pump 2 powers are fixed, when the pump 1 power is increased, the reflectivity increases in a way that is very different from the moving FBG model. The trend could be approximately considered as a cubic curve by expanding the exponential term in Eq. (11) and keeping to the second order within the power range used for calculation, since the gain is not too high. When the pump 1 power is further increased, from Eq. (11), one could expect that the exponential growth would dominate until strong pump depletion occurs where numerical solution is needed for completely solving Eq. (3). However, this similar trend has already been observed in Ref. [6], and we will also show this in our experiment in the next section.

 

Fig. 3 Theoretically calculated BDG reflectivity with respect to (a) pump 1 power |A1|2 for different pump 2 and probe powers (b) pump 2 power |A2|2 for different probe powers and fixed pump 1 power (c) probe power |A3|2 for different pump 1 powers and fixed pump 2 power.

Download Full Size | PPT Slide | PDF

Figure 3(b) shows that when pump 1 and probe powers are fixed, the reflectivity of the BDG should show a linear increasing trend according to Eq. (11). Moveover, when the probe power is higher, the linear slope becomes steeper as well. The BDG reflectivity as a function of probe power is plotted in Fig. 3(c). For higher probe power, the diffracted wave undergoes a stronger SBS amplification which should result in a higher reflectivity. In the next section, we will demonstrate the experimental results (the power levels of the pump and probe waves are chosen to be consistent with the undepleted pump 1 and probe waves approximation), and show that all these trends are observed confirming that the BDG generation and detection process is indeed consistent with an FWM process.

3. Experiment validation

Figure 4 shows the experiment setup. Two narrow linewidth (∼ 3 kHz) fiber lasers operating at 1550 nm are used to provide pumps 1 and 2, and their frequency difference is locked by a phase-locking loop in a frequency counter via a photodetector (PD1). The frequency of Laser 1 is higher than Laser 2. A tunable laser (Agilent 81940A) whose linewidth is ∼ 100 kHz with a continuous sweep mode is used as the probe wave. Three polarization controllers (PCs) are used to adjust the polarization states of the three waves. A variable attenuator is used to adjust the power level of the pump 2 wave, while two erbium-doped fiber amplifiers (EDFAs) are employed to adjust power levels of pump 1 and probe waves. Pump 1 and pump 2 waves are launched into the PMF along slow axis in opposite direction to each other through polarization beam splitter/combiner (PBS/PBC), and the probe wave is injected along the orthogonal axis. Power meter 2 (PM2) is used to monitor pump 2 after passing through the PMF in order to determine whether a strong pump depletion happens. In the detection end, a 125 MHz photodetector (PD2) and a high performance oscilloscope (Agilent Infiniium DSO81204B) are used to detect the back diffracted wave via a narrow bandwidth (0.1 nm) tunable filter (TF), whose function is to eliminate the leakage of pump 1 and 2 waves along the slow axis. After the frequency of the two pump waves are locked at the Brillouin frequency of the PMF to generate the BDG, the probe wave is fast swept and the diffracted spectrum of the BDG is recorded by the oscilloscope (A trigger is given to the oscilloscope when the sweep starts for synchronization).

 

Fig. 4 Experiment Setup. PC: polarization controller; PD: photodetector; EDFA: erbium-doped fiber amplifier; ISO: isolator; PBS/PBC: polarization beam splitter/combiner; PM: power meter; TF: tunable filer; PMF: polarization-maintaining fiber.

Download Full Size | PPT Slide | PDF

Two kinds of standard PMFs, 13-m Panda and 10-m Bow-Tie fibers whose Brillouin frequencies are 10875 MHz and 10814 MHz at room temperature respectively, are tested in the experiment. Both fibers are under loose condition. When the frequency difference of the tunable laser and pump 1 wave, Δν, is near 50 GHz for Panda fiber and 43 GHz for Bow-Tie fiber, strong diffracted signal could be observed. Typical spectra of the corresponding BDGs generated in these two kinds of fibers are shown in Figs. 5(a) and 5(b). We intently set the position of highest peak as 0 MHz for comparison. The sweep speed of the tunable laser is 40 nm/s corresponding to ∼ 5 MHz/μs. The noise is reduced in the way that high resolution sampling mode of the oscilloscope is chosen; therefore, each BDG spectrum could be obtained within 200 μs for a 1 GHz sweep range. The frequency resolution of the measured BDG spectrum is ∼ 0.1 MHz (1 GHz/8003 points). Both spectra show multi-peaks due to the non-uniformity of the birefringence inside the PMF. This similar spectrum has been reported recently in a 5-m long PMF (Ref. [14]). In the other previous reports, the BDG spectrum shows only one broad peak, which is attributed to the relative frequency fluctuation between different laser sources. However, in our experiment, the sweep speed is fast enough to avoid such fluctuation. Actually, the sweep speed could be changed from 50 nm/s down to 10 nm/s, the spectrum remains the same for each sweep. However, when the sweep speed is further reduced, the spectrum starts to be different for each sweep, because the frequency fluctuation of the pump lasers cannot be ignored when the sweeping speed of the tunable laser is low. It is worth mentioning that the birefringence non-uniformity of the Panda fiber used in our experiment is much smaller than that of the Bow-Tie fiber, which could be confirmed by Fig. 5, since the overall width of spectrum in Fig. 5(a) is much narrower than that in Fig. 5(b). In the following experiment, we use the Panda fiber to measure the reflectivity of the BDG. As an example, the BDG spectrum changing with respect to the variation of the pump 1 power is shown in Fig. 6, where the pump 2 and probe powers are kept at 10 and 100 mW, respectively. It is clearly seen that the shape of the spectra remain unchanged except that the reflectivity is increasing when the pump 1 power increases. Note that although the birefringence non-uniformity could be suppressed by applying longitudinal strain to the PMF with the coating stripped completely [20], the spectrum obtained in our experiment is good enough to confirm the FWM process analyzed theoretically in last section by finding the maximum reflectivity of the BDG which is obtained by locating the highest peak in the spectrum.

 

Fig. 5 Typical BDG spectra for (a) 13-m Panda and (b) 10-m Bow-Tie fiber.

Download Full Size | PPT Slide | PDF

 

Fig. 6 BDG spectra for different pump 1 power, with the pump 2 power of 10 mW and probe power of 100 mW.

Download Full Size | PPT Slide | PDF

Figure 7(a) shows the BDG reflectivity with respect to pump 1 power for different pump 2 and probe powers. All the trends clearly differ from a linear increasing (the curves are the cubic fits of the data points obtained experimentally) which is consistent with the theoretical calculation. We believe that further increasing the pump 1 power would show a more obvious exponential trend as that in Ref. [6], but the power level is limited by the maximum handling power of our optical components in CW operation. The reflectivity obtained experimentally is much smaller than calculated results because the birefringence of the PMF is nonuniform which results in phase mismatch. In addition, we did not account for other power losses of the signal in the configuration, especially for the TF which causes at least a 3-dB loss. However, this trend already provides a clear evidence that the BDG generation and detection is actually an FWM process, especially for the case that the probe power is comparable to the pump 1 power.

 

Fig. 7 Experimentally obtained BDG reflectivity with respect to (a) pump 1 power |A1|2 for different pump 2 and probe powers (b) pump 2 power |A2|2 for different probe powers and fixed pump 1 power (c) probe power |A3|2 for different pump 1 powers and fixed pump 2 power.

Download Full Size | PPT Slide | PDF

Figure 7(b) shows the BDG reflectivity with respect to pump 2 power for different probe powers with pump 1 power fixed. Both curves which are linear fits of the data points obtained experimentally show a well linear increase of the BDG reflectivity versus pump 2 power, which is also observed in Ref. [6]. For higher probe power, the reflectivity is higher and the slope of the linear increase is steeper as well which is consistent with the theoretical calculation. The reflectivity of the BDG with respect to the probe power is shown in Fig. 7(c). When the probe power is increasing, the BDG reflectivity is higher according to the theoretical analysis which is due to the stronger FWM process. For the moving FBG model, the reflectivity would be constant with respect to the variation of probe power. However, since the BDG generation and detection process in a PMF is actually an FWM process, the reflection of BDG is dependent on the probe power which was ignored in all the previous reports. Even though the present work is focusing on CW case, we believe that for transit case, e.g., time-domain measurement of BDG spectrum (probe pulse follows one of the pump pulses) [7, 16, 21], the reflectivity is also dependent on probe power, since a strong probe pulse would result in an SBS amplification of the diffracted wave.

4. Conclusion

We have theoretically investigated the BDG generation and detection process in a PMF for the case of CW operation. The moving FBG model oversimplifies this process since the results derived from this model deviates from some of the experimental reports. By deriving the four coupled light wave equations, we could obtain analytical solutions under the undepleted pump and probe waves approximation. The trends of both the theoretical calculation and the experiment agree with each other very well showing that the reflectivity of the BDG is dependent on the probe power. This is a clear evidence indicating that the BDG generation and detection process is indeed an FWM process.

Appendix A: derivation of the four coupled wave equations

The theoretical model diagram is shown in Fig. 1. We start the analysis with the material density equation, Eq. (1). The total electric field, tot, on the right-hand side of Eq. (1) is expressed as Eqs. (2) and (3). Since acoustic wave is predominantly longitudinal, we neglect the transversal components of the acoustic wave for SBS process. Then, the right-hand of Eq. (1) takes the form,

2E˜tot212F2[A˜1A˜2*q12ei(q1zΩt)+A˜3A˜4*q22ei(q2zΩt)]+c.c.,
where Ω = ω1ω2 = ω3ω4, q1 = k1 + k2 and q2 = k3 + k4. Then, we are looking for the solution of the Eq. (1) in the form,
ρ˜=ρ˜1+ρ˜2=12FA(x,y)[Q1(z,t)ei(q1zΩt)+Q2(z,t)ei(q2zΩt)]+c.c.
where Q1(z,t) and Q2(z,t) are the slowly varying amplitudes of the acoustic waves corresponding to the two SBS processes. FA(x, y) is the dimensionless acoustic wave mode profile which is the solution of the unperturbed Eq. (1) and thus satisfies the eigenvalue equation [23],
2FA(x,y)+(ΩBm2vA2qm2)FA(x,y)=0,(m=1,2)
where 2=x2+y2 is the transverse Laplacian operator; ΩBm is the eigenfrequency of the solution of the modal acoustic equation that for a given wave vector qm which satisfies Eq. (A.3). Assuming steady-state condition with slowly varying envelope approximation, substituting Eqs. (A.1) and (A.2) into Eq. (1) and grouping the results by the exponential factors ei(q1z−Ωt) and ei(q2z−Ωt), one could obtain two equations for Q1 and Q2, respectively,
12{iΩΓB1FAQ1+vA2Q1[2FA+(Ω2vA2q12)FA]}=14ɛ0γeF2q12A˜1A˜2*,
12{iΩΓB2FAQ2+vA2Q2[2FA+(Ω2vA2q22)FA]}=14ɛ0γeF2q22A˜3A˜4*,
where ΓBm=qm2ΓA. Two additional approximations have been used in order to obtain Eq. (A.4): the term containing ∂Qm/∂z has been dropped since the hypersonic phonons are strongly damped and thus propagate only over very shot distances before being absorbed [17]; since the ratio ΩΓA/vA2103 is very small, iΩΓA+vA2 could be approximated by its real part [23], vA2 in the term which contains the derivative 2FA. Inserting Eq. (A.3) into Eq. (A.4), multiplying both sides of each resulting equation with FA(x, y) and integrating over the transverse plane, one could obtain that,
Q1(z,t)=ɛ0γeq12A˜1(z,t)A˜2*(z,t)2(ΩB12Ω2iΩΓB1)F2(x,y)FA(x,y)FA2(x,y),
Q2(z,t)=ɛ0γeq22A˜3(z,t)A˜4*(z,t)2(ΩB22Ω2iΩΓB2)F2(x,y)FA(x,y)FA2(x,y),
where angular brackets denote averaging over transverse cross section of the fiber. Note that the eigenfrequency of the acoustic mode ΩBm corresponds to the full resonance occurring both for frequencies and for wave vectors. Strictly speaking, it is the solution of the set of equations
q1=ω1neff,x(ω1)c+(ω1Ω)neff,x(ω1)c,
q2=ω3neff,y(ω3)c+(ω3Ω)neff,y(ω3)c,
and
2FA(x,y)+(Ω2vA2qm2)FA(x,y)=0,(m=1,2)
where neff,x and neff,y are the effective refractive indices of the optical mode along and ŷ axis, respectively. Moreover, approximations have been used as neff,x(ω2) ≈ neff,x(ω1) and neff,y(ω4) ≈ neff,y(ω3), since Ω ≪ ω1, ω3. Thus, for a fixed frequency ω1 or ω3, ΩBm could be calculated as the intersection of the two curves of qm(Ω) obtained from Eqs. (A.6) and (A.7). However, since qm(Ω) is a very weak function because of Ω ≪ ω1, ω3,
q1neff,x(ω1)c(2ω1Ω)2ω1neff,x(ω1)c2ω1nxc,
q2neff,y(ω3)c(2ω3Ω)2ω3neff,y(ω3)c2ω3nyc,
with nx and ny being the refractive index of the PMF along and ŷ axis, respectively. The obtained qm value could be substituted into Eq. (A.7) to find ΩBm directly. Moreover, as discussed before, for BDG generation and detection process, phase matching is required as that Δk = 0 which is equivalent to q1 = q2 = q; therefore, we could take ΩB1 = ΩB2 = ΩB, ΓB1 = ΓB2 = ΓB for further calculations. From Eqs. (A.2) and (A.5), the full material density variation could be expressed as,
ρ˜(x,y,z,t)=14FAF2FAFA2ɛ0γeq2(A˜1A˜2*+A˜3A˜4*)ΩB2Ω2iΩΓBei(qzΩt)+c.c..

Next, we start to derive equations for the evolution of the four optical waves. Maxwell’s equations can be deduced to for the electric field j(j = 1 − 4) inside the PMF as follows:

2E˜jn2c22E˜jt2=1ɛ0c22P˜jt2,
where n = nx for j = 1, 2 and n = ny for j = 3, 4. The total nonlinear polarization induced by the acousto-optic interaction which gives rise to the source term in this equation is given by
P˜tot=ɛ0γeρ0ρ˜E˜tot=ɛ0ɛNLE˜tot,
where
ɛNL=FAU{A˜1A˜2*ei[(k1+k2)zΩt]+A˜3A˜4*ei[(k3+k4)zΩt]}+c.c.
with
U=ɛ0γe2q24ρ0(ΩB2Ω2iΩΓB)F2FAFA2.
By determining different parts of tot that can act as phase-matched source terms for the four waves, one could obtain
P˜1=12ɛ0UFFA(A˜1|A˜2|2+A˜2A˜3A˜4*eiΔkz)ei(k1zω1t)+c.c.,
P˜2=12ɛ0U*FFA(|A˜1|2A˜2+A˜1A˜3*A˜4eiΔkz)ei(k2zω2t)+c.c.,
P˜3=12ɛ0UFFA(A˜3|A˜4|2+A˜1A˜2*A˜4eiΔkz)ei(k3zω3t)+c.c.,
P˜4=12ɛ0U*FFA(|A˜3|2A˜4+A˜1*A˜2A˜3eiΔkz)ei(k4zω4t)+c.c.,
where Δk = k3 + k4k1k2 is the phase mismatch. Neglecting the dependence of the optical modal profile on nonlinear effects, we can substitute the optical modal equation
2F(x,y)+(ωj2n2c2kj2)F(x,y)=0,
along with Eqs. (A.14) into Eq. (A.10), using slowly varying envelope approximation, and multiplying F(x,y) on both sides of each resultant equation, integrating over the transverse cross section, to obtain four equations for four interacting optical waves under steady state,
A1z=iη1(A1|A2|2+A2A3A4*eiΔkz),
A2z=iη1*(|A1|2A2+A1A3*A4eiΔkz),
A3z=iη2(A3|A4|2+A1A2*A4eiΔkz),
A4z=iη2*(|A3|2A4+A1*A2A3eiΔkz),
where
η1=γe2ω13ρ0c4(ΩB2Ω2iΩΓB)Aeffao,
η2=γe2ω33ρ0c4(ΩB2Ω2iΩΓB)Aeffao,
with Aeffao being given by Eq. (5). Note that Eq. (A.16) has been normalized as that |Aj|2 represents optical power by using |Aj|2 = ɛ0cn|Ãj|2F2〉/2. In addition, ω1ω2 and ω3ω4 are used to obtain Eqs. (A.16) and (A.17), and nxny is assumed for normalization. When the frequency difference of the two pump waves has been locked at Brillouin frequency, e.g., Ω = ΩB, Eq. (A.16) reduces to Eq. (3) where a further approximation of λ1λ3 (wavelengths of pump 1 and probe waves) is considered.

Appendix B: analytical solution of the four coupled wave equations under undepleted pump 1 and probe waves approximation

We calculate z derivative of Eq. (9b) with the help of (9a), considering that A1 and A3 are constants, and then the resultant equation for A4 could be expressed as follows:

2A4z2+ΦA4zΨA4=0,
where Φ = κ(|A1|2 + |A3|2) – iΔk and Ψ = |A3|2Δk. The general solution for it is
A4(z)=Cep+z+Depz,
with p±=12Φ±12Φ2+4Ψ and two constants C and D which could be determined by boundary conditions:
A4|z=L=0,
A4z|z=L=κA1*A2(L)A3eiΔkL.
Equation (B.3b) is obtained from Eq. (9b) after insertion of Eq. (B.3a). Substituting Eq. (B.2) into Eq. (B.3), the two constants could be determined as
C=κA1*A2(L)A3eiΔkLep+Lpp+,
D=κA1*A2(L)A3eiΔkLepLpp+.
Therefore, substituting Eq. (B.4) into Eq. (B.2), the diffracted wave A4(z) could be obtained as expressed in Eq. (10). Under phase-matching condition, Δk = 0, so that Φ = κ(|A1|2 + |A3|2) and Ψ = 0; the reflectivity could be expressed in Eq. (11). Note that the expression for pump 2 wave A2 could also be obtained in the similar way.

Acknowledgments

The authors would like to thank Natural Science and Engineering Research Council of Canada (NSERC) Discovery Grants and Canada Research Chair Program for the financial support. D. P. Zhou would like to acknowledge the Province of Ontario Ministry of Research and Innovation and the University of Ottawa for the financial support of the Vision 2020 Postdoctoral Fellowship. The authors would like to thank Dr. Wenhai Li for helpful discussions.

References and links

1. R. Y. Chiao, C. H. Townes, and B. P. Stoicheff, “Stimulated Brillouin scattering and coherent generation of intense hypersonic waves,” Phys. Rev. Lett. 12, 592–595 (1964). [CrossRef]  

2. Z. Zhu, D. J. Gauthier, and R. W. Boyd, “Stored light in an optical fiber via stimulated Brillouin scattering,” Science 318, 1748–1750 (2007). [CrossRef]   [PubMed]  

3. L. Thévenaz, “Slow and fast light in optical fibres,” Nat. Photonics 2, 474–481 (2008). [CrossRef]  

4. G. P. Agrawal, “Nonlinear Fiber Optics,” 4th edition, Academic Press2007.

5. X. Bao and L. Chen, “Recent progress in Brillouin scattering based fiber sensors,” Sensors 11, 4152–4187 (2011). [CrossRef]  

6. K. Y. Song, W. Zou, Z. He, and K. Hotate, “All-optical dynamic grating generation based on Brillouin scattering in polarization-maintaining fiber,” Opt. Lett. 33, 926–928 (2008). [CrossRef]   [PubMed]  

7. Y. Dong, X. Bao, and L. Chen, “Distributed temperature sensing based on birefringence effect on transient Brillouin grating in a polarization-maintaining photonic crystal fiber,” Opt. Lett. 34, 2590–2592 (2009). [CrossRef]   [PubMed]  

8. K. Y. Song and H. J. Yoon, “High-resolution Brillouin optical time domain analysis based on Brillouin dynamic grating,” Opt. Lett. 35, 52–54 (2010). [CrossRef]   [PubMed]  

9. K. Y. Song, S. Chin, N. Primerov, and L. Thévenaz, “Time-domain distributed fiber sensor with 1 cm spatial resolution based on Brillouin dynamic grating,” J. Lightwave Technol. 28, 2062–2067 (2010). [CrossRef]  

10. S. Chin, N. Primerov, and L. Thévenaz, “Sub-centimetre spatial resolution in distributed fibre sensing, based on dynamic Brillouin grating in optical fibers,” to appear in IEEE Sens. J. (2011). [CrossRef]  

11. W. Zou, Z. He, and K. Hotate, “Complete discrimination of strain and temperature using Brillouin frequency shift and birefringence in a polarization-maintaining fiber,” Opt. Express 17, 1248–1255 (2009). [CrossRef]   [PubMed]  

12. W. Zou, Z. He, and K. Hotate, “Demonstration of Brillouin distributed discrimination of strain and temperature using a polarization-maintaining optical fiber,” IEEE Photon. Technol. Lett. 22, 526–528 (2010). [CrossRef]  

13. Y. Dong, L. Chen, and X. Bao, “High-spatial-resolution time-domain simultaneous strain and temperature sensor using Brillouin scattering and birefringence in a polarization-maintaining fiber,” IEEE Photon. Technol. Lett. 22, 1364–1366 (2010). [CrossRef]  

14. W. Zou, Z. He, and K. Hotate, “One-laser-based generation/detection of Brillouin dynamic grating and its application to distributed discrimination of strain and temperature,” Opt. Express 19, 2363–2370 (2011). [CrossRef]   [PubMed]  

15. K. Y. Song, K. Lee, and S. B. Lee, “Tunable optical delays based on Brillouin dynamic grating in optical fibers,” Opt. Express 17, 10344–10349 (2009). [CrossRef]   [PubMed]  

16. Y. Dong, L. Chen, and X. Bao, “Truly distributed birefringence measurement of polarization-maintaining fibers based on transient Brillouin grating,” Opt. Lett. 35, 193–195 (2010). [CrossRef]   [PubMed]  

17. R. W. Boyd, Nonlinear Optics, 3rd ed. (Academic Press, 2008).

18. T. Erdogan, “Fiber grating spectra,” J. Lightwave Technol. 15, 1277–1294 (1997). [CrossRef]  

19. Y. Dong, L. Chen, and X. Bao, “Characterization of the Brillouin grating spectra in a polarization-maintaining fiber,” Opt. Express 18, 18960–18967 (2010). [CrossRef]   [PubMed]  

20. K. Y. Song and H. J. Yoon, “Observation of narrowband intrinsic spectra of Brillouin dynamic gratings,” Opt. Lett. 35, 2958–2960 (2010). [CrossRef]   [PubMed]  

21. K. Y. Song, W. Zou, Z. He, and K. Hotate, “Optical time-domain measurement of Brillouin dynamic grating spectrum in a polarization-maintaining fiber,” Opt. Lett. 34, 1381–1383 (2009). [CrossRef]   [PubMed]  

22. A. Kobyakov, S. Kumar, D. Q. Chowdhury, A. B. Ruffin, M. Sauer, S. R. Bichham, and R. Mishra, “Design concept for optical fibers with enhaced SBS threshold,” Opt. Express 13, 5338–5346 (2005). [CrossRef]   [PubMed]  

23. A. Kobyakov, M. Sauer, and D. Chowdhury, “Stimulated Brillouin scattering in optical fibers,” Adv. Opt. Photon. 2, 1–59 (2010). [CrossRef]  

References

  • View by:
  • |
  • |
  • |

  1. R. Y. Chiao, C. H. Townes, and B. P. Stoicheff, “Stimulated Brillouin scattering and coherent generation of intense hypersonic waves,” Phys. Rev. Lett. 12, 592–595 (1964).
    [Crossref]
  2. Z. Zhu, D. J. Gauthier, and R. W. Boyd, “Stored light in an optical fiber via stimulated Brillouin scattering,” Science 318, 1748–1750 (2007).
    [Crossref] [PubMed]
  3. L. Thévenaz, “Slow and fast light in optical fibres,” Nat. Photonics 2, 474–481 (2008).
    [Crossref]
  4. G. P. Agrawal, “Nonlinear Fiber Optics,” 4th edition, Academic Press2007.
  5. X. Bao and L. Chen, “Recent progress in Brillouin scattering based fiber sensors,” Sensors 11, 4152–4187 (2011).
    [Crossref]
  6. K. Y. Song, W. Zou, Z. He, and K. Hotate, “All-optical dynamic grating generation based on Brillouin scattering in polarization-maintaining fiber,” Opt. Lett. 33, 926–928 (2008).
    [Crossref] [PubMed]
  7. Y. Dong, X. Bao, and L. Chen, “Distributed temperature sensing based on birefringence effect on transient Brillouin grating in a polarization-maintaining photonic crystal fiber,” Opt. Lett. 34, 2590–2592 (2009).
    [Crossref] [PubMed]
  8. K. Y. Song and H. J. Yoon, “High-resolution Brillouin optical time domain analysis based on Brillouin dynamic grating,” Opt. Lett. 35, 52–54 (2010).
    [Crossref] [PubMed]
  9. K. Y. Song, S. Chin, N. Primerov, and L. Thévenaz, “Time-domain distributed fiber sensor with 1 cm spatial resolution based on Brillouin dynamic grating,” J. Lightwave Technol. 28, 2062–2067 (2010).
    [Crossref]
  10. S. Chin, N. Primerov, and L. Thévenaz, “Sub-centimetre spatial resolution in distributed fibre sensing, based on dynamic Brillouin grating in optical fibers,” to appear in IEEE Sens. J. (2011).
    [Crossref]
  11. W. Zou, Z. He, and K. Hotate, “Complete discrimination of strain and temperature using Brillouin frequency shift and birefringence in a polarization-maintaining fiber,” Opt. Express 17, 1248–1255 (2009).
    [Crossref] [PubMed]
  12. W. Zou, Z. He, and K. Hotate, “Demonstration of Brillouin distributed discrimination of strain and temperature using a polarization-maintaining optical fiber,” IEEE Photon. Technol. Lett. 22, 526–528 (2010).
    [Crossref]
  13. Y. Dong, L. Chen, and X. Bao, “High-spatial-resolution time-domain simultaneous strain and temperature sensor using Brillouin scattering and birefringence in a polarization-maintaining fiber,” IEEE Photon. Technol. Lett. 22, 1364–1366 (2010).
    [Crossref]
  14. W. Zou, Z. He, and K. Hotate, “One-laser-based generation/detection of Brillouin dynamic grating and its application to distributed discrimination of strain and temperature,” Opt. Express 19, 2363–2370 (2011).
    [Crossref] [PubMed]
  15. K. Y. Song, K. Lee, and S. B. Lee, “Tunable optical delays based on Brillouin dynamic grating in optical fibers,” Opt. Express 17, 10344–10349 (2009).
    [Crossref] [PubMed]
  16. Y. Dong, L. Chen, and X. Bao, “Truly distributed birefringence measurement of polarization-maintaining fibers based on transient Brillouin grating,” Opt. Lett. 35, 193–195 (2010).
    [Crossref] [PubMed]
  17. R. W. Boyd, Nonlinear Optics, 3rd ed. (Academic Press, 2008).
  18. T. Erdogan, “Fiber grating spectra,” J. Lightwave Technol. 15, 1277–1294 (1997).
    [Crossref]
  19. Y. Dong, L. Chen, and X. Bao, “Characterization of the Brillouin grating spectra in a polarization-maintaining fiber,” Opt. Express 18, 18960–18967 (2010).
    [Crossref] [PubMed]
  20. K. Y. Song and H. J. Yoon, “Observation of narrowband intrinsic spectra of Brillouin dynamic gratings,” Opt. Lett. 35, 2958–2960 (2010).
    [Crossref] [PubMed]
  21. K. Y. Song, W. Zou, Z. He, and K. Hotate, “Optical time-domain measurement of Brillouin dynamic grating spectrum in a polarization-maintaining fiber,” Opt. Lett. 34, 1381–1383 (2009).
    [Crossref] [PubMed]
  22. A. Kobyakov, S. Kumar, D. Q. Chowdhury, A. B. Ruffin, M. Sauer, S. R. Bichham, and R. Mishra, “Design concept for optical fibers with enhaced SBS threshold,” Opt. Express 13, 5338–5346 (2005).
    [Crossref] [PubMed]
  23. A. Kobyakov, M. Sauer, and D. Chowdhury, “Stimulated Brillouin scattering in optical fibers,” Adv. Opt. Photon. 2, 1–59 (2010).
    [Crossref]

2011 (3)

X. Bao and L. Chen, “Recent progress in Brillouin scattering based fiber sensors,” Sensors 11, 4152–4187 (2011).
[Crossref]

S. Chin, N. Primerov, and L. Thévenaz, “Sub-centimetre spatial resolution in distributed fibre sensing, based on dynamic Brillouin grating in optical fibers,” to appear in IEEE Sens. J. (2011).
[Crossref]

W. Zou, Z. He, and K. Hotate, “One-laser-based generation/detection of Brillouin dynamic grating and its application to distributed discrimination of strain and temperature,” Opt. Express 19, 2363–2370 (2011).
[Crossref] [PubMed]

2010 (8)

W. Zou, Z. He, and K. Hotate, “Demonstration of Brillouin distributed discrimination of strain and temperature using a polarization-maintaining optical fiber,” IEEE Photon. Technol. Lett. 22, 526–528 (2010).
[Crossref]

Y. Dong, L. Chen, and X. Bao, “High-spatial-resolution time-domain simultaneous strain and temperature sensor using Brillouin scattering and birefringence in a polarization-maintaining fiber,” IEEE Photon. Technol. Lett. 22, 1364–1366 (2010).
[Crossref]

A. Kobyakov, M. Sauer, and D. Chowdhury, “Stimulated Brillouin scattering in optical fibers,” Adv. Opt. Photon. 2, 1–59 (2010).
[Crossref]

K. Y. Song and H. J. Yoon, “High-resolution Brillouin optical time domain analysis based on Brillouin dynamic grating,” Opt. Lett. 35, 52–54 (2010).
[Crossref] [PubMed]

Y. Dong, L. Chen, and X. Bao, “Truly distributed birefringence measurement of polarization-maintaining fibers based on transient Brillouin grating,” Opt. Lett. 35, 193–195 (2010).
[Crossref] [PubMed]

K. Y. Song, S. Chin, N. Primerov, and L. Thévenaz, “Time-domain distributed fiber sensor with 1 cm spatial resolution based on Brillouin dynamic grating,” J. Lightwave Technol. 28, 2062–2067 (2010).
[Crossref]

Y. Dong, L. Chen, and X. Bao, “Characterization of the Brillouin grating spectra in a polarization-maintaining fiber,” Opt. Express 18, 18960–18967 (2010).
[Crossref] [PubMed]

K. Y. Song and H. J. Yoon, “Observation of narrowband intrinsic spectra of Brillouin dynamic gratings,” Opt. Lett. 35, 2958–2960 (2010).
[Crossref] [PubMed]

2009 (4)

2008 (2)

2007 (1)

Z. Zhu, D. J. Gauthier, and R. W. Boyd, “Stored light in an optical fiber via stimulated Brillouin scattering,” Science 318, 1748–1750 (2007).
[Crossref] [PubMed]

2005 (1)

1997 (1)

T. Erdogan, “Fiber grating spectra,” J. Lightwave Technol. 15, 1277–1294 (1997).
[Crossref]

1964 (1)

R. Y. Chiao, C. H. Townes, and B. P. Stoicheff, “Stimulated Brillouin scattering and coherent generation of intense hypersonic waves,” Phys. Rev. Lett. 12, 592–595 (1964).
[Crossref]

Agrawal, G. P.

G. P. Agrawal, “Nonlinear Fiber Optics,” 4th edition, Academic Press2007.

Bao, X.

Bichham, S. R.

Boyd, R. W.

Z. Zhu, D. J. Gauthier, and R. W. Boyd, “Stored light in an optical fiber via stimulated Brillouin scattering,” Science 318, 1748–1750 (2007).
[Crossref] [PubMed]

R. W. Boyd, Nonlinear Optics, 3rd ed. (Academic Press, 2008).

Chen, L.

Chiao, R. Y.

R. Y. Chiao, C. H. Townes, and B. P. Stoicheff, “Stimulated Brillouin scattering and coherent generation of intense hypersonic waves,” Phys. Rev. Lett. 12, 592–595 (1964).
[Crossref]

Chin, S.

S. Chin, N. Primerov, and L. Thévenaz, “Sub-centimetre spatial resolution in distributed fibre sensing, based on dynamic Brillouin grating in optical fibers,” to appear in IEEE Sens. J. (2011).
[Crossref]

K. Y. Song, S. Chin, N. Primerov, and L. Thévenaz, “Time-domain distributed fiber sensor with 1 cm spatial resolution based on Brillouin dynamic grating,” J. Lightwave Technol. 28, 2062–2067 (2010).
[Crossref]

Chowdhury, D.

Chowdhury, D. Q.

Dong, Y.

Erdogan, T.

T. Erdogan, “Fiber grating spectra,” J. Lightwave Technol. 15, 1277–1294 (1997).
[Crossref]

Gauthier, D. J.

Z. Zhu, D. J. Gauthier, and R. W. Boyd, “Stored light in an optical fiber via stimulated Brillouin scattering,” Science 318, 1748–1750 (2007).
[Crossref] [PubMed]

He, Z.

Hotate, K.

Kobyakov, A.

Kumar, S.

Lee, K.

Lee, S. B.

Mishra, R.

Primerov, N.

S. Chin, N. Primerov, and L. Thévenaz, “Sub-centimetre spatial resolution in distributed fibre sensing, based on dynamic Brillouin grating in optical fibers,” to appear in IEEE Sens. J. (2011).
[Crossref]

K. Y. Song, S. Chin, N. Primerov, and L. Thévenaz, “Time-domain distributed fiber sensor with 1 cm spatial resolution based on Brillouin dynamic grating,” J. Lightwave Technol. 28, 2062–2067 (2010).
[Crossref]

Ruffin, A. B.

Sauer, M.

Song, K. Y.

Stoicheff, B. P.

R. Y. Chiao, C. H. Townes, and B. P. Stoicheff, “Stimulated Brillouin scattering and coherent generation of intense hypersonic waves,” Phys. Rev. Lett. 12, 592–595 (1964).
[Crossref]

Thévenaz, L.

S. Chin, N. Primerov, and L. Thévenaz, “Sub-centimetre spatial resolution in distributed fibre sensing, based on dynamic Brillouin grating in optical fibers,” to appear in IEEE Sens. J. (2011).
[Crossref]

K. Y. Song, S. Chin, N. Primerov, and L. Thévenaz, “Time-domain distributed fiber sensor with 1 cm spatial resolution based on Brillouin dynamic grating,” J. Lightwave Technol. 28, 2062–2067 (2010).
[Crossref]

L. Thévenaz, “Slow and fast light in optical fibres,” Nat. Photonics 2, 474–481 (2008).
[Crossref]

Townes, C. H.

R. Y. Chiao, C. H. Townes, and B. P. Stoicheff, “Stimulated Brillouin scattering and coherent generation of intense hypersonic waves,” Phys. Rev. Lett. 12, 592–595 (1964).
[Crossref]

Yoon, H. J.

Zhu, Z.

Z. Zhu, D. J. Gauthier, and R. W. Boyd, “Stored light in an optical fiber via stimulated Brillouin scattering,” Science 318, 1748–1750 (2007).
[Crossref] [PubMed]

Zou, W.

Adv. Opt. Photon. (1)

IEEE Photon. Technol. Lett. (2)

W. Zou, Z. He, and K. Hotate, “Demonstration of Brillouin distributed discrimination of strain and temperature using a polarization-maintaining optical fiber,” IEEE Photon. Technol. Lett. 22, 526–528 (2010).
[Crossref]

Y. Dong, L. Chen, and X. Bao, “High-spatial-resolution time-domain simultaneous strain and temperature sensor using Brillouin scattering and birefringence in a polarization-maintaining fiber,” IEEE Photon. Technol. Lett. 22, 1364–1366 (2010).
[Crossref]

IEEE Sens. J. (1)

S. Chin, N. Primerov, and L. Thévenaz, “Sub-centimetre spatial resolution in distributed fibre sensing, based on dynamic Brillouin grating in optical fibers,” to appear in IEEE Sens. J. (2011).
[Crossref]

J. Lightwave Technol. (2)

Nat. Photonics (1)

L. Thévenaz, “Slow and fast light in optical fibres,” Nat. Photonics 2, 474–481 (2008).
[Crossref]

Opt. Express (5)

Opt. Lett. (6)

Phys. Rev. Lett. (1)

R. Y. Chiao, C. H. Townes, and B. P. Stoicheff, “Stimulated Brillouin scattering and coherent generation of intense hypersonic waves,” Phys. Rev. Lett. 12, 592–595 (1964).
[Crossref]

Science (1)

Z. Zhu, D. J. Gauthier, and R. W. Boyd, “Stored light in an optical fiber via stimulated Brillouin scattering,” Science 318, 1748–1750 (2007).
[Crossref] [PubMed]

Sensors (1)

X. Bao and L. Chen, “Recent progress in Brillouin scattering based fiber sensors,” Sensors 11, 4152–4187 (2011).
[Crossref]

Other (2)

G. P. Agrawal, “Nonlinear Fiber Optics,” 4th edition, Academic Press2007.

R. W. Boyd, Nonlinear Optics, 3rd ed. (Academic Press, 2008).

Cited By

OSA participates in Crossref's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1
Fig. 1

Diagram of the theoretical model. BDG generation and detection scheme: an FWM process.

Fig. 2
Fig. 2

Properties of the tanh2(x) and tanh 2 ( x ) functions.

Fig. 3
Fig. 3

Theoretically calculated BDG reflectivity with respect to (a) pump 1 power |A1|2 for different pump 2 and probe powers (b) pump 2 power |A2|2 for different probe powers and fixed pump 1 power (c) probe power |A3|2 for different pump 1 powers and fixed pump 2 power.

Fig. 4
Fig. 4

Experiment Setup. PC: polarization controller; PD: photodetector; EDFA: erbium-doped fiber amplifier; ISO: isolator; PBS/PBC: polarization beam splitter/combiner; PM: power meter; TF: tunable filer; PMF: polarization-maintaining fiber.

Fig. 5
Fig. 5

Typical BDG spectra for (a) 13-m Panda and (b) 10-m Bow-Tie fiber.

Fig. 6
Fig. 6

BDG spectra for different pump 1 power, with the pump 2 power of 10 mW and probe power of 100 mW.

Fig. 7
Fig. 7

Experimentally obtained BDG reflectivity with respect to (a) pump 1 power |A1|2 for different pump 2 and probe powers (b) pump 2 power |A2|2 for different probe powers and fixed pump 1 power (c) probe power |A3|2 for different pump 1 powers and fixed pump 2 power.

Equations (51)

Equations on this page are rendered with MathJax. Learn more.

2 ρ ˜ t 2 Γ A 2 ρ ˜ t v A 2 2 ρ ˜ = 1 2 ɛ 0 γ e 2 E ˜ tot 2 ,
E ˜ x = 1 2 x ^ [ F ( x , y ) A ˜ 1 ( z , t ) e i ( k 1 z ω 1 t ) + F ( x , y ) A ˜ 2 ( z , t ) e i ( k 2 z ω 2 t ) ] + c . c . ,
E ˜ y = 1 2 y ^ [ F ( x , y ) A ˜ 3 ( z , t ) e i ( k 3 z ω 3 t ) + F ( x , y ) A ˜ 4 ( z , t ) e i ( k 4 z ω 4 t ) ] + c . c . ,
A 1 z = κ ( A 1 | A 2 | 2 + A 2 A 3 A 4 * e i Δ k z ) ,
A 2 z = κ ( | A 1 | 2 A 2 + A 1 A 3 * A 4 e i Δ kz ) ,
A 3 z = κ ( A 3 | A 4 | 2 + A 1 A 2 * A 4 e i Δ k z ) ,
A 4 z = κ ( | A 3 | 2 A 4 + A 1 * A 2 A 3 e i Δ k z ) ,
κ = 8 π 3 γ e 2 c ρ 0 λ 3 Ω B Γ B A eff ao ,
A eff ao = [ F 2 ( x , y ) F 2 ( x , y ) F A ( x , y ) ] 2 F A 2 ( x , y ) ,
A 3 z = κ A 1 A 2 * A 4 e i Δ kz ,
A 4 z = κ A 1 * A 2 A 3 e i Δ kz ,
A 4 ( z ) = 2 κ exp ( i Δ kz / 2 ) sinh [ g ( L z ) / 2 ] A 1 * A 2 A 3 ( 0 ) g cosh ( gL / 2 ) i Δ k sinh ( gL / 2 ) ,
R FBG = | A 4 ( 0 ) A 3 ( 0 ) | 2 = tanh 2 ( κ | A 1 | | A 2 | L ) .
A 2 z = κ ( | A 1 | 2 A 2 + A 1 A 3 * A 4 e i Δ k z ) ,
A 4 z = κ ( | A 3 | 2 A 4 + A 1 * A 2 A 3 e i Δ kz ) ,
A 4 ( z ) = 2 κ A 1 * A 2 ( L ) A 3 sinh [ Φ 2 + 4 Ψ ( L z ) / 2 ] Φ 2 + 4 Ψ e Φ ( L z ) / 2 + i Δ kL ,
R FWM = | A 4 ( 0 ) A 3 | 2 = | A 1 | 2 | A 2 ( L ) | 2 [ e κ ( | A 1 | 2 + | A 3 | 2 ) L 1 | A 1 | 2 + | A 3 | 2 ] 2 .
2 E ˜ tot 2 1 2 F 2 [ A ˜ 1 A ˜ 2 * q 1 2 e i ( q 1 z Ω t ) + A ˜ 3 A ˜ 4 * q 2 2 e i ( q 2 z Ω t ) ] + c . c . ,
ρ ˜ = ρ ˜ 1 + ρ ˜ 2 = 1 2 F A ( x , y ) [ Q 1 ( z , t ) e i ( q 1 z Ω t ) + Q 2 ( z , t ) e i ( q 2 z Ω t ) ] + c . c .
2 F A ( x , y ) + ( Ω Bm 2 v A 2 q m 2 ) F A ( x , y ) = 0 , ( m = 1 , 2 )
1 2 { i Ω Γ B 1 F A Q 1 + v A 2 Q 1 [ 2 F A + ( Ω 2 v A 2 q 1 2 ) F A ] } = 1 4 ɛ 0 γ e F 2 q 1 2 A ˜ 1 A ˜ 2 * ,
1 2 { i Ω Γ B 2 F A Q 2 + v A 2 Q 2 [ 2 F A + ( Ω 2 v A 2 q 2 2 ) F A ] } = 1 4 ɛ 0 γ e F 2 q 2 2 A ˜ 3 A ˜ 4 * ,
Q 1 ( z , t ) = ɛ 0 γ e q 1 2 A ˜ 1 ( z , t ) A ˜ 2 * ( z , t ) 2 ( Ω B 1 2 Ω 2 i Ω Γ B 1 ) F 2 ( x , y ) F A ( x , y ) F A 2 ( x , y ) ,
Q 2 ( z , t ) = ɛ 0 γ e q 2 2 A ˜ 3 ( z , t ) A ˜ 4 * ( z , t ) 2 ( Ω B 2 2 Ω 2 i Ω Γ B 2 ) F 2 ( x , y ) F A ( x , y ) F A 2 ( x , y ) ,
q 1 = ω 1 n eff , x ( ω 1 ) c + ( ω 1 Ω ) n eff , x ( ω 1 ) c ,
q 2 = ω 3 n eff , y ( ω 3 ) c + ( ω 3 Ω ) n eff , y ( ω 3 ) c ,
2 F A ( x , y ) + ( Ω 2 v A 2 q m 2 ) F A ( x , y ) = 0 , ( m = 1 , 2 )
q 1 n eff , x ( ω 1 ) c ( 2 ω 1 Ω ) 2 ω 1 n eff , x ( ω 1 ) c 2 ω 1 n x c ,
q 2 n eff , y ( ω 3 ) c ( 2 ω 3 Ω ) 2 ω 3 n eff , y ( ω 3 ) c 2 ω 3 n y c ,
ρ ˜ ( x , y , z , t ) = 1 4 F A F 2 F A F A 2 ɛ 0 γ e q 2 ( A ˜ 1 A ˜ 2 * + A ˜ 3 A ˜ 4 * ) Ω B 2 Ω 2 i Ω Γ B e i ( qz Ω t ) + c . c ..
2 E ˜ j n 2 c 2 2 E ˜ j t 2 = 1 ɛ 0 c 2 2 P ˜ j t 2 ,
P ˜ tot = ɛ 0 γ e ρ 0 ρ ˜ E ˜ tot = ɛ 0 ɛ NL E ˜ tot ,
ɛ NL = F A U { A ˜ 1 A ˜ 2 * e i [ ( k 1 + k 2 ) z Ω t ] + A ˜ 3 A ˜ 4 * e i [ ( k 3 + k 4 ) z Ω t ] } + c . c .
U = ɛ 0 γ e 2 q 2 4 ρ 0 ( Ω B 2 Ω 2 i Ω Γ B ) F 2 F A F A 2 .
P ˜ 1 = 1 2 ɛ 0 U F F A ( A ˜ 1 | A ˜ 2 | 2 + A ˜ 2 A ˜ 3 A ˜ 4 * e i Δ k z ) e i ( k 1 z ω 1 t ) + c . c . ,
P ˜ 2 = 1 2 ɛ 0 U * F F A ( | A ˜ 1 | 2 A ˜ 2 + A ˜ 1 A ˜ 3 * A ˜ 4 e i Δ k z ) e i ( k 2 z ω 2 t ) + c . c . ,
P ˜ 3 = 1 2 ɛ 0 U F F A ( A ˜ 3 | A ˜ 4 | 2 + A ˜ 1 A ˜ 2 * A ˜ 4 e i Δ kz ) e i ( k 3 z ω 3 t ) + c . c . ,
P ˜ 4 = 1 2 ɛ 0 U * F F A ( | A ˜ 3 | 2 A ˜ 4 + A ˜ 1 * A ˜ 2 A ˜ 3 e i Δ kz ) e i ( k 4 z ω 4 t ) + c . c . ,
2 F ( x , y ) + ( ω j 2 n 2 c 2 k j 2 ) F ( x , y ) = 0 ,
A 1 z = i η 1 ( A 1 | A 2 | 2 + A 2 A 3 A 4 * e i Δ kz ) ,
A 2 z = i η 1 * ( | A 1 | 2 A 2 + A 1 A 3 * A 4 e i Δ kz ) ,
A 3 z = i η 2 ( A 3 | A 4 | 2 + A 1 A 2 * A 4 e i Δ kz ) ,
A 4 z = i η 2 * ( | A 3 | 2 A 4 + A 1 * A 2 A 3 e i Δ kz ) ,
η 1 = γ e 2 ω 1 3 ρ 0 c 4 ( Ω B 2 Ω 2 i Ω Γ B ) A eff ao ,
η 2 = γ e 2 ω 3 3 ρ 0 c 4 ( Ω B 2 Ω 2 i Ω Γ B ) A eff ao ,
2 A 4 z 2 + Φ A 4 z Ψ A 4 = 0 ,
A 4 ( z ) = C e p + z + D e p z ,
A 4 | z = L = 0 ,
A 4 z | z = L = κ A 1 * A 2 ( L ) A 3 e i Δ kL .
C = κ A 1 * A 2 ( L ) A 3 e i Δ kL e p + L p p + ,
D = κ A 1 * A 2 ( L ) A 3 e i Δ kL e p L p p + .

Metrics