Model for a Torsional-Mode Ultrasonic Transducer for an Acousto-Optic In-Fiber Isolator

Torsional ultrasound modes can couple the optical polarization states in a birefringent fiber. Polarization coupling produced by interaction with a higher-order torsional mode slightly above its cutoff may provide a route to producing an in-fiber isolator suitable for use in narrow-band high-power fiber amplifiers. This paper describes a model of a transducer generating torsional modes in a cylindrical fiber. This model predicts that almost all of the power applied to the transducer is radiated into the desired mode. The paper also discusses effects produced by acoustic absorption and the dependence of the acoustic velocity on temperature.


Introduction
By matching the optical polarization beat length L B in a birefringent silica optical fiber to the axial wavelength of a traveling torsional wave of angular frequency Ω propagating in the fiber, it is possible to couple light from one polarization to the other [1,2].When light is coupled from the slow polarization to the fast polarization, it is down-shifted in frequency by Ω if the light and sound are copropagating, but up-shifted if the light and sound are counter-propagating.This nonreciprocal frequency shift can be used to make an infiber acousto-optic (AO) isolator.This requires placing two torsional transducers on the fiber separated by D = πv g /2Ω, where v g is the group velocity of light in the fast polarization.These transducers propagate sound in the same direction.One watt of acoustic power from each transducer is sufficient to couple 50% of the light from one polarization to the other over a few centimeters [2].The polymer jacket is removed from the fiber in the interaction regions to prevent acoustic damping.The relative acoustic phase of the two transducers is adjusted such that in the forward direction all of the light at the wavelength of interest returns to the slow polarization, whereas in the backward direction all of the light at this wavelength goes into the fast polarization.Polarizers at both ends, such as 45-degree tilted fiber Bragg gratings [3], remove light in the fast polarization to make an isolator.To date an isolator based on torsional acoustic waves has not been made, but 100% polarization switching has been demonstrated using the lowest-order (s = 0) torsional mode [1].However, an analogous isolator was demonstrated using the lowestorder flexural mode at Ω/2π = 4.31 MHz [4].In this case coupling was between the two transverse modes of a twomode fiber.Although a 22 dB extinction ratio was reported, the large D = 11 m spacing between the transducers makes this approach unattractive for high-power operation.The key to decreasing this spacing is to increase Ω, while using an acoustic mode for which phase matching is possible at high frequency.Raising the frequency to about 50 MHz reduces the spacing D between transducers to about 1 m.Except for the short regions following the transducers where acoustic coupling occurs, the fiber between the transducers can be coiled and jacketed with polymer.Figure 1 shows the layout of an AO isolator at four different length scales.
The phase-matching acceptance for polarization switching, which occurs over a few centimeters, is much wider than the transmission features produced by the transducer spacing D and the linear polarizers.The transmission in the forward direction as a function of wavelength λ is a comb with periodicity λL B /D, while the backward transmission is an out-of-phase comb.This is illustrated in Figure 2.    The position of these combs depends on the relative acoustic phase of the two transducers.Hwang et al. [4] found that it was necessary to use a feedback circuit to adjust this phase to compensate for shifting of the combs produced by vibrations or temperature changes.The width of the comb maxima and minima is sufficient to prevent backward propagation of light produced by stimulated Brillouin scattering in narrowband Yb-doped fiber amplifiers.The isolator can be used for multiple wavelengths commensurate with the comb period, such asthose that might be produced by spectral beam combination.However, it cannot protect against inadvertent lasing at other wavelengths within the gain bandwidth of Yb.This type of isolator can be tuned for use with other types of fiber amplifier besides Yb, such as Er or Tm.
A transducer spacing of 1 m imposes a power limitation on operation of an AO isolator, since undesirable nonlinear effects can build up over this distance.For example, we estimate that the nonlinear refractive index of fused silica will limit isolator operation to about 200 W in a 10 μm core fiber and about 1 kW in a 25 μm core fiber.By contrast, present fiber-coupled isolators are limited to about 50 W. Bulk isolators can withstand higher power but distort the beam quality at powers above about 100 W.
The frequency for phase-matching of the lowest-order torsional mode is too low for use as a high-power isolator, since D becomes very large and the comb spacing very narrow.However, a transducer producing 100% polarization coupling with this mode can be used as a tunable notch filter [1], when placed between crossed polarizers.This can filter out undesired wavelengths, including wavelengths within the gain bandwidth of the amplifying medium.

Matching of Transducer Modes and Free-Oscillation Modes
In the analysis below we consider the fiber to be a homogeneous elastic cylindrical glass rod, neglecting perturbations which might arise because of different elastic properties of the core or the stress rods.Torsional waves in a cylindrical rod are characterized by a displacement u(r, z) which is entirely in the azimuthal direction, but independent of the azimuthal coordinate θ.If the cylinder radius is a, we can introduce a dimensionless frequency γ = Ωa/V , where V is the transverse sound velocity.This velocity is variously reported by Engan to be 3750 m/s [2] or 3764 m/s [5] in fused silica.Vacher et al. [6] found that V at 5 MHz and 16 GHz are the same to within 0.1% at room temperature, suggesting that dispersion is small.However, V was found to increase with temperature at 0.3 m/s kelvin.It is convenient to express r and z in units of a. Then the torsional waves are of the form u = (1/2)r exp(±iγz) or u = J 1 (αr) exp(±iβz), where The torsional modes of a free elastic cylinder are a subset of all the possible vibrational modes [7].In general the displacement u has an azimuthal dependence cos(nθ) or sin(nθ).The modes with n = 0 can be further subdivided into torsional and axial-radial.The modes with n = 1 or higher are called flexural.The dispersion relations for the modes, except for the torsional modes, depend on the Poisson ratio for the material.Figure 3 shows the dispersion relations for n = 0 (red) and n = 1 (green) for the Poisson ratio 0.17 of fused silica.The lowest and highest red curves shown are the two lowest torsional modes.The other red curves are axial-radial.The lowest torsional, radial-axial, and flexural modes are the only ones with no cutoff.For torsional oscillations the two nonzero components of strain (with proportional components of stress) are 2ε rθ = r∂(u/r)∂r and 2ε zθ = ∂u/∂z [8].Our model for the transducer is a section of fiber between z = −L and z = 0 where u at r = 1 is a singlefrequency oscillation independent of z.The model does not specify how such oscillation is produced.However, at the end of the paper we cite some methods that have been used to excite torsional vibrations in fibers and rods.ε zθ = 0 at the surface of the fiber in the transducer, whereas ε rθ = 0 at the surface of the free fiber.The displacement in the transducer is of the form where , and κ n is positive real for propagating modes or positive imaginary for cutoff modes.D n and H n are the modal amplitudes.We could multiply the first term above by a constant, but we choose to normalize in this way for now.Outside the transducer the fiber is uncontrained at the surface; so ε rθ should vanish at r = 1.Since ((rd/dr))[(1/r)J 1 (αr)] = −αJ 2 (αr), we choose α = Z n , where Z n is the n'th zero of J 2 , and write the displacement for z > 0 in the form where n and β n is positive real for propagating modes or positive imaginary for cutoff modes.For positive z the propagating modes with modal amplitudes C n carry power only in the positive z direction.The mode with a cutoff at γ = Z n we call the s = n mode.Since Z 1 = 5.136, the Figure 3: Dispersion relations for torsional, radial-axial, and flexural modes for an elastic cylinder.Lowest and highest red curves are torsional; other red curves are axial-radial; green curves are flexural.β is longitudinal wave vector in units of 1/a.γ is angular frequency in units of V/a.s = 1 cutoff frequency for a fiber of diameter 2a = 125 μm is Ω/2π = VZ 1 /2πa = 49 MHz.By operating just above the cutoff for the s = n mode, we can match the axial wavelength 2πa/β n to the optical polarization beat length L B (typically about 3 mm).The phase velocity Vγ/β n approaches infinity at cutoff, while the sound group velocity Vβ n /γ approaches zero.The s = 0 mode, corresponding to the C 0 term, has no cutoff, propagates at speed V without dispersion, and can only be phase matched at much lower frequency.Since Z 1 is greater than Y 1 = 3.832, a propagating mode inside the transducer already exists when γ passes through the s = 1 cutoff.
At z < −L we write the displacement as ( Here F n are the modal amplitudes for sound propagating in the −z direction.The last two terms with amplitudes G 0 or G p are inserted to describe reflections of one of the modes above cutoff back into the transducer.For example, if the s = 1 mode is above cutoff, we take G 1 to be nonzero and G 0 = 0.The case of nonzero G 0 is not of much practical interest but is included for completeness.When the G coefficients are zero, equal acoustic power is radiated from the transducer in both directions.The backward power is undesirable.It can be damped by leaving the jacket on the fiber at z < −L.However, more efficient operation might be possible by reflecting the power.For a mode near cutoff, reflection can be produced by raising the fiber temperature behind the transducer.This increases V , which puts the mode below cutoff in this region.
If the reflection is total, then |G p | = |F p |, but the effects of the reflection depend on the phase Θ of G p .This phase depends on the location of the thermal gradient producing the reflection and cycles through one period when this location shifts by half the longitudinal sound wavelength in the unheated fiber.Under phase-matched conditions this is also half the optical polarization beat length L B .
Reflection should be produced close to the transducer in order to minimize AO coupling by light going in the −z direction.At z = −L and z = 0 both u and ∂u/∂z should be continuous.This leads to the following set of four equations: Here we have introduced Note that in the absence of reflection (G 0 = 0, G p = 0) we can take B n = 0, F 0 = C 0 , and F n = C n , so that emission from the transducer is symmetrical.
The torsional modes outside the transducer are orthogonal, as can be seen from the following identities: Multiplying ( 4)-( 7) by either r or J 1 (Z m r), integrating, and using the above identities, we obtain We can eliminate F 0 , F m , C 0 , and C m from these equations to obtain a set of equations to be obeyed by A n and another set of equations to be obeyed by B n .In fact, if we define Z 0 = 0 and β 0 = γ, the equations derived from ( 12)-( 15) can be treated on the same footing as those derived from ( 16)- (19).
We obtain If one first specifies a value of G p (or G 0 ), these equations constitute a set of linear inhomogeneous equations for the amplitudes A n and B n .In order to solve these equations we choose an integer N (typically 200) and truncate by letting n run from 1 to N and m run from 0 to N − 1.Then we obtain A n and B n by numerical matrix inversion.The coefficients F 0 , F m , C 0 , and C m needed for the modal expansions of u outside the transducer are then obtained from ( 12), ( 14), (16), and To check the convergence of our numerical procedure, we can plot the real and imaginary parts of the left and right sides of ( 4)- (7) as a function of r to see if they match.Equations ( 4) and ( 6) represent the displacements u(−L, r) and u(0, r), while (5) and ( 7) are proportional to the strain ε zθ .We can also calculate −r∂(u/r)∂r from ( 4) and ( 6) to obtain expressions proportional to the strain ε rθ .The real parts correspond to times when the transducer amplitude (the displacement at r = 1) is maximal, while the imaginary parts correspond to times when the transducer amplitude is zero.In making these plots the number of terms we retain in the sums (typically 140) is less than N, since the final terms can be corrupted by the matrix truncation.
We find that the matching of the external and internal displacements and strains is good.The differences appear as low-amplitude, high-frequency oscillations.The differences are most apparent in ε zθ , which is zero at r = 1 in the transducer but appears to be discontinuous as a function of z at z = −L and z = 0 when r = 1.The displacements and strains are insensitive to the choice of N if N is not too small.

Energy Considerations
Since our model neglects sound absorption in the fiber, the power radiated in the modes which are above cutoff should equal the power P T applied to the transducer.Modes below cutoff are exponentially attenuated away from the transducer and do not propagate power.Let U be the amplitude of the transducer motion at r = 1.Then the modal coefficients calculated in the previous section should be multiplied by U/|J 1 (γ)| to determine the physical displacement amplitudes.By comparison with the power expressions given in [2], we see that the power propagating in the s = 0 mode at z > 0 is P s=0,+ = [γ|C 0 |/2J 1 (γ)] 2 P , where P = (1/4)πρV 3 U 2 and ρ = 2200 kg/m 2 is the density of fused silica.If, for example, U = 90 nm, then P = 0.738 W. The power in higher-order modes which are above cutoff is P s=n,+ = 2[|C n |J 1 (Z n )/J 1 (γ)] 2 β n γP .Analogous expressions apply to the power propagating at z < −L.
Division by |J 1 (γ)| is not possible when γ = Y n for some n.In this case there is no steady state in the absence of damping.When Y 1 < γ < Z 1 , there is one propagating mode inside the transducer.We call this the t = 1 mode.This mode can radiate only to the s = 0 mode outside the transducer, and resonances occur as a function of L. Internal values of u are typically large compared to u at the surface, and the t = 1 power within the transducer can be much larger than the power driving the transducer.When γ exceeds Z 1 , the s = 1 channel opens, sound vibrations are no longer trapped in the transducer, the t = 1 power is small, and the generated s = 1 power is almost independent of L. In general the t = n power propagating in the forward and backward direction is given for modes above cutoff by To calculate the time-averaged power driving the transducer we use the stress evaluated at the surface of the fiber.The time-averaged power dP T driving the transducer over the surface area a 2 dθdz is obtained by multiplying this stress by the surface velocity −UΩ sin(Ωt) and averaging over time, +H n exp(−iκ n z) . (24) Integrating over the surface of the transducer and using ( 8) and ( 10) yields It is often more meaningful to express powers relative to P T than in terms of P .

Numerical Results
Let us begin by considering a transducer with L = 10 and γ = 5.13767, slightly above the s = 1 cutoff.For a fiber with a diameter of 2a = 125 μm this corresponds to a frequency Ω/2π = 49.06MHz and an axial s = 1 wavelength 2πa/β 1 = 2.7 mm, on the order of the polarization beat length L B in polarization-maintaining optical fibers.One watt of power in this mode corresponds to an amplitude at r = 1 of 85.7 nm.At this frequency the s = 0 wavelength 2πa/γ = 76.4 μm is much too small for phase matching.We wish to compare the situation where the sound is emitted equally in both directions to the case where the s = 1 mode is perfectly reflected at z = −L with phase Θ = 0 or 180 • .U is held constant in Figures 4 and 5, but P T depends on whether there is reflection.The red curves are at z = 0 with reflection, the green curves are at z = −L with reflection, and the blue curves are at either end with no reflection.
In the case Θ = 0, P T = 1.47P with reflection, but P T = 2.96P with no reflection.Thus, reflection acts to reduce the power needed to drive the transducer to a given amplitude U but does not greatly affect the way the transducer vibrates.The modal efficiencies with reflection are P s=0,+ /P T = 0.56%, P s=1,+ /P T = 98.64%, and P s=0,− /P T = 0.07%.P s=1,− /P T = 25.76%propagates in each direction between the point of reflection and z = −L.Within the transducer P t=1,+ /P T = 0.028% and P t=1,− /P T = 0.776%.The amplitude at r = 1 of the s = 1 mode is 0.985U.
In the case Θ = 180 • , P T = 120.60P; so it becomes much harder to drive the transducer to a given amplitude U.The modal efficiencies are P s=0,+ /P T = 58.76%,P s=1,+ /P T = 4.94%, and P s=0,− /P T = 36.18%.Thus, the transducer becomes a poor emitter of s = 1 ultrasound, although a large power P s=1,− /P T = 1943.6%propagates in each direction in the region z < −L between the point of reflection and the left end of the transducer.Within the transducer P t=1,+ /P T = 63.71% and P t=1,− /P T = 2.02%.The transducer motion is qualitatively quite different because of the reflection.
In Figure 6 we plot the efficiency P s=1,+ /P T as a function of Θ.We see that the efficiency is high over a broad range.On the other hand, without reflection 49.73% of the applied power P T is emitted into the s = 1 mode in each direction, and the undesirable backward sound must be damped to prevent interaction with light in the fiber core.
The transducer simulations described below are for no reflections at z < −L.Above the s = 1 cutoff the applied power P T and the radiated power are only weakly dependent on L. Figure 7 shows the s = 1 efficiency for each direction as a function of L at γ = 5.13767.The s = 0 and internal t = 1 powers are much smaller than the s = 1 power.The situation is different below the s = 1 cutoff, where acoustic energy can leave the transducer only in the s = 0 channel.This is illustrated in Figure 8 for γ = 5.12, such that 2π/κ 1 = 1.85.While the efficiency for generating the s = 0 mode in each direction remains at 1/2 (red), P T goes to zero at resonances when L is a multiple of 2π/κ 1 (blue).Consequently the ratio P t=1 /P T becomes large near these resonances (green).We expect that damping effects not included in our model will act to suppress these resonances, but it may be best to avoid them by a suitable choice of L.
In order to increase the frequency Ω, so as to shorten the distance between the two transducers in an isolator, or to use a larger-diameter fiber, it is of interest to consider increasing γ above the cutoff for higher-order modes.Figure 9 shows the s = 2 efficiency as a function of L at γ = 8.41849, slightly above the Z 2 = 8.41724 cutoff (red).This figure also shows the s = 3 efficiency at γ = 11.62075,slightly above the Z 3 = 11.61984cutoff (green).We see that, as desired, most of the acoustic power goes into the mode that is just above cutoff.Scaling to large Ω will eventually be limited by the increased attenuation of ultrasound at very high frequencies.
At γ = 8.416, just below the s = 2 cutoff, there is no dominance of s = 1 over s = 0. Figure 10 shows P s=0 /P T (red) and P s=1 /P T (green), as well as the internal resonances P t=1 /P T (blue) and P t=2 /P T (violet).The s = 0 and s = 1 efficiencies in each direction sum to 1/2.

Acousto-Optical Coupling, Acoustic Attenuation, and Thermal Effects
According to Engan [2] the coupling of the optical polarization states by phase-matched torsional acoustic waves is sinusoidal.That is, one polarization amplitude increases as sin(K n z), while the other decreases as cos(K n z).Here z is position (not in units of a) along the fiber.We call K n z the coupling angle.The coupling constant K n for the s = n mode, provided that mode is above cutoff, is where P is the acoustic power injected at z = 0 into the phase-matched mode.Coupling of 50% of the light from one polarization to the other occurs over a distance L c = π/(4K n ).For 1 W of power in a fiber with a 2.7 mm polarization beat length and a 125 μm diameter, the coupling length for the s = 1 mode is L c = 3.54 cm and the phase-matching frequency is 49.06 MHz.For the same power and beat length, but a 250 μm diameter, the coupling length for the s = 4 mode is L c = 3.58 cm and the phase-matching frequency is 70.66 MHz.
When one takes damping of the sound waves into account, the power decreases exponentially after the sound leaves the transducer at z = 0.This means that β n has an imaginary part even above cutoff.The coupling angle at position This means that a coupling angle of π/4 is only possible if L c < L a , where the absorption length is defined as L a = a/ Im(β n ).The acoustic power decays as exp(−2z/L a ).When there is damping the distance over which sufficient coupling occurs is greater than L c .
We assume that the damping of the s = 0 mode is the same as that found for transverse sound waves in bulk fused Advances in OptoElectronics silica.Room-temperature measurements by Krischer [9] in the 200 to 1000 MHz range found the power attenuation to be 5.7 × 10 −18 dB/μsec( f /Hz) 2 , where f = Ω/2π.This agrees well with the value given by Kul'bitskaya and Shutilov [10] but is smaller than the value 1.4 × 10 −17 dB/μsec( f /Hz) 2 reported by Fraser et al. [11].Measurements by Krause [12] at 590 MHz showed a decrease in attenuation with increasing temperature with a minimum at about 375 C. For the present discussion we assume that Krischer's result can be extrapolated to lower frequencies.Dividing by V and neglecting dispersion in the speed of sound, we find that Im(β 0 )/a = 4.43 × 10 −18 Ω 2 sec 2 /m or β 0 = γ + i(a 0 /a)γ 2 , where a 0 = 6.63 × 10 −11 m.
For the s = 0 mode β 0 /a is the total wave vector, and the transverse wave vector is zero.For higher-order modes we take β n = (β 2 0 − Z 2 n ) 1/2 , since the transverse wave vector is Z n /a.Note that damping smears out the cutoff, and that β n has both real and imaginary parts on either side of the cutoff.Suppose that the beat length is L B = 2.7 mm.For the example above of 2a = 125 μm and f = 49.06MHz, the absorption length for the s = 1 mode is L a = 6.32 cm.The minimum power required at z = 0 to achieve 50% optical coupling is P = 0.313 W. The s = 2 mode in this fiber at 80.39 MHz has L a = 1.44 cm and requires at least P = 2.36 W. For the example 2a = 250 μm and f = 70.66MHz, the absorption length for the s = 4 mode is L a = 2.12 cm and at least 2.87 W in this mode is required.Absorbed acoustic power can produce a temperature gradient along the fiber.This can cause a chirp in β n , adversely affecting phase matching with the optical field.On the other hand, the phase-matching acceptance increases when the AO coupling length decreases.Heating can also be caused by optical power in the cladding.Cladding light can be removed before the isolator by embedding a section of fiber in high-index material.Air flow past the fiber can also reduce temperature gradients.
Even in the absence of temperature gradients and acoustic damping, the temperature dependence of V [6] affects   the cutoff frequencies and the phase-matching conditions.Figure 11  the distance over which the power decreases by 1/e 2 .If the fiber diameter were to be slightly different at the locations of the two transducers in an isolator, phase matching for both transducers at the same frequency could be maintained by having the transducers at different temperatures.Backward acoustic power produced by a transducer can be reflected by heating the fiber behind the transducer over a sufficient length to prevent tunneling through the region where the reflected mode is below cutoff.The phase Θ of the reflected light depends on how and where this heating is done.

Discussion
Using higher-order torsional acoustic modes near cutoff to couple optical polarization states in a birefringent fiber allows phase matching to the polarization beat length at high frequencies in the 50 MHz range.Two torsionalmode transducers producing 50% polarization coupling and separated on the fiber by about 1 m can be used as an in-fiber optical isolator suitable for use with narrow-band fiber amplifiers of arbitrary wavelength.We have presented an analytical and numerical model of a torsional-mode transducer.This model indicates that the mode just above cutoff radiates most of the acoustic power applied to the  transducer, irrespective of the transducer length.In the absence of acoustic reflection the power is radiated equally in both directions, and undesirable backward radiation must be absorbed.Reflecting the backward radiation towards the transducer with a suitable phase can enhance the efficiency in the forward direction to almost 100%.Various types of transducers have been developed for exciting torsional vibrations in fibers or rods, but not at the high frequencies needed for an AO isolator.Transducers for exciting the s = 0 mode in optical fibers have used a conical glass horn with the fiber mounted along its axis and shear-mode piezoelectric transducers attached to its base [1,[13][14][15].Other torsional transducers have been developed which use paired [16] or segmented [17,18] piezoelectric elements.Torsional oscillations in metallic rods have been produced by electromagnetic-acoustic transduction [19].Silicon-based micro-electromechanical systems (MEMSs) technology may provide a route to making higher-frequency Advances in OptoElectronics torsional transducers.We are working with Sandia National Laboratories to explore this possibility.

Figure 1 :
Figure 1: Components of an AO isolator are shown schematically at four different length scales.Sizes are approximate.Forward propagating light enters from the left with slow polarization.Polarizers (possibly 45 • tilted fiber Bragg gratings) are oriented to pass the slow polarization.Copropagating torsional acoustic waves generated by transducers act as 50% polarization couplers.Transducer 2 returns forward light to the slow polarization, while transducer 1 shifts backward light entirely to the fast polarization.

Figure 2 :
Figure 2: Calculated forward (red) and backward (green) transmission through an isolator using fiber with beat length L B = 3 mm and transducer separation D = 105 cm.
. If, for example, one is interested in the situation where the p mode is totally reflected at z < −L, we repeat the calculation for different values of |G p | in order to home in on the value of |G p | for which |G p | = |F p |.The result depends on the phase Θ of G p , which in turn depends on how and where the reflection occurs.

7 −Figure 4 :
Figure 4: Displacement u(r) (a) and (d), strain ε rθ (r) (b) and (e), and strain ε zθ (r) (c) and (f) are shown at the two ends of the transducer.Transverse coordinate r is relative to the fiber radius a. Scaling is arbitrary.The normalized frequency is γ = 5.13767.Red (z = 0) and green (z = −L) curves are for unit reflectivity of the s = 1 mode back towards z = −L with phase Θ = 0. Blue curve is for zero reflectivity and the same transducer amplitude U. (a), (b), (c), and (d) show the real parts, corresponding to times when u(1) = U. (d), (e), and (f) show the imaginary parts, corresponding to times when u(1) = 0. High-frequency oscillations in some of the curves are a numerical artifact.Since ε zθ (r) in the absence of reflection has opposite signs at z = 0 and z = −L, the right-hand green plots have been multiplied by minus 1 for better comparison with the blue plots.

Figure 5 :
Figure 5: Displacement u(r) (a) and (d), strain ε rθ (r) (b) and (e), and strain ε zθ (r) (c) and (f) are shown at the two ends of the transducer.This is for the same conditions as in Figure4, except that the phase of the reflected s = 1 mode is 180 • .

Figure 6 :
Figure 6: Efficiency P s=1,+ /P T as a function of the phase Θ of the reflected s = 1 mode at z = −L.

Figure 7 :
Figure 7: Efficiency for s = 1 generation at γ = 5.13767.Transducer length L is in units of a.

Figure 8 :Figure 9 :
Figure 8: Efficiency for s = 0 generation (red) at γ = 5.12 remains at 1/2, but P T (blue) becomes small near resonances.As a result the internal t = 1 efficiency P t=1 /P T (green) is large at the resonances.

Figure 11 :
Figure 11: Axial wavelength (left of cutoff resonances) and 1/e 2 attenuation length (right of resonances) is shown as a function of temperature for several fiber diameters.The frequency is chosen to produce an axial wavelength of 3 mm at 20 C for a 125.0 μmdiameter fiber.