Optical nanofibers (ONFs) represent versatile nanophotonic platforms for important photonic applications such as optical sensing and quantum and nonlinear optics. The attractiveness of ONFs arises from the tight optical confinement, their wide evanescent field in the subwavelength limit, their surface acoustic properties, and their high tensile strength. Here we investigate Brillouin light scattering in silica-glass ONFs under high tensile strain and show that the fundamental properties of elastic waves dramatically change due to elastic anisotropy and nonlinear elasticity for strain larger than 2%. This yields to unexpected Brillouin strain coefficients for all Brillouin resonances including surface and hybrid waves, followed by a nonlinear evolution at high tensile strength. We further provide a complete theoretical analysis based on third-order nonlinear elasticity of silica that agrees well with our experimental data. These new regimes open the way to the development of compact tensile strain optical sensors based on nanofibers.

Optical micro- and nanofibers (ONFs) are long and uniform ultrathin fibers manufactured by heating and tapering standard optical fibers down to the submicrometer scale.1–5 In addition to providing tight optical confinement, they exhibit a strong evanescent field in the subwavelength limit, which is very attractive for applications such as optical sensing and quantum photonics.3–8 They also possess notable mechanical and elastic properties, with large extensibility and high tensile strength.9 From an acoustic viewpoint, it has been recently shown that ONFs support a new class of acoustic waves compared to standard fibers owing to the strong coupling between shear and longitudinal waves.10,11 These waves include hybrid acoustic waves (HAWs) and surface acoustic waves (SAWs) that move at a lower speed than the longitudinal acoustic velocity. This in turn gives rise to a multipeaked Brillouin spectrum with a Brillouin frequency shift ranging from 5 GHz for SAWs and up to 10 GHz for HAWs.

In this work, we investigate the tensile strain dependence of optical nanofibers using Brillouin spectroscopy and we show that the fundamental elastic properties dramatically change because of strain-induced elastic anisotropy. More specifically, we find that the Brillouin strain coefficients associated with SAWs and HAWs are different from those observed in standard optical fibers ranging from 30% to 110% of a single-mode fiber. We report the observation of Brillouin frequency shifts featuring different slopes and crossings. Unlike in standard fibers, they do not scale with the Brillouin frequency shift and even show a strong nonlinear deviation at high tensile strain. For further understanding, we develop a theoretical model based on third-order elasticity of silica to predict the strain dependence of all acoustic waves in tapered optical fibers including the nanofiber section and the transition regions. We find a very good agreement with experimental results.

This article is divided into four sections. We briefly review Brillouin scattering (BS) in optical fibers and the dependence of the Brillouin frequency shift with respect to strain in Sec. II. Section III describes the experimental measurements of Brillouin strain coefficients in optical nanofibers with different diameters and tensile strain conditions. In Sec. IV, we provide the analytic model based on a nonlinear Hooke’s law to clearly interpret our observations. Finally, our results are further compared with finite element method (FEM)-based numerical simulations to highlight both the nanofiber and the transition region contributions to the Brillouin spectrum.

In 1925, two years after discovering the scattering of light through elastic waves in homogeneous transparent materials,12 Léon Brillouin theoretically predicted that the speed of elastic waves depends on strain because of the nonlinear elasticity and later provided the magnitude and sign of the third-order elastic constants.13,14 These predictions were experimentally verified a few years later with ultrasonic velocity measurements,15 and the first measurement of all third-order elastic constants of fused silica was reported in 1965,16 followed by other measurements still using ultrasonic techniques.17 Brillouin light scattering was later exploited for measuring the longitudinal and transverse velocities in many bulk materials.18–20 In optical fibers, tensile strain measurements lead to the third-order elastic constant C111.21 All these fundamental research works have opened up the development of long range and high-resolution strain and temperature Brillouin distributed optical fiber sensors that are widely and commercially used today for integrity and security in civil engineering and the petroleum industry.22 Brillouin-based fiber optical sensors principally exploit the sensitivity of the Brillouin gain spectrum (BGS) in single-mode fibers to temperature and strain.23–25 However, given the dimension of a standard single-mode optical fiber (SMF) and the acoustic frequency involved in Brillouin scattering, only the bulk longitudinal acoustic wave is probed by Brillouin sensors. The first consequence is that the derivative of the Brillouin frequency shift νB with respect to the tensile strain is given by

(1)

where C11 and C111 are the longitudinal second-order and third-order elastic constants, respectively (see Table I), ϵ is the tensile strain, usually defined in percentage of fiber elongation, and the refractive index dependence to ϵ is neglected. For instance, in standard optical silica fibers, the Brillouin frequency shift is about 11 GHz and thus the strain coefficient is Cϵ ≃ 0.04 MHz/μϵ at a wavelength of 1550 nm, which can also be written as 400 MHz/%.22 Equation (1) generally assumes that only longitudinal mechanical properties are sensed by the apparatus. In standard optical fibers, the maximal tensile strain is limited to 2%, due to the polymer cladding, but naturally the elasticity of silica glass allows going up to 6%.26 

TABLE I.

Optical and elastic parameters of fused silica used for nonlinear elasticity model and numerical simulations. Pij are electrostrictive tensor components. Cαβ and Cαβγ are second and third order elastic tensor components, respectively. Data are reproduced with permission from E. H. Bogardus, J. Appl. Phys. 36, 2504 (1965). Copyright 1965 AIP Publishing LLC; B. E. Powell and M. J. Skove, J. Appl. Phys. 41, 4913 (1970). Copyright 1970 AIP Publishing LLC; and W. T. Yost and M. A. Breazeale, J. Appl. Phys. 44, 1909 (1973). Copyright 1973 AIP Publishing LLC.

Refractive index n 1.444 C11 78 GPa 
Density ρ 2203 kg/m3 C44 31 GPa 
P11 0.12 C111 578 GPa 
P12 0.27 C112 215 GPa 
P44 −0.073 C123 43 GPa 
Refractive index n 1.444 C11 78 GPa 
Density ρ 2203 kg/m3 C44 31 GPa 
P11 0.12 C111 578 GPa 
P12 0.27 C112 215 GPa 
P44 −0.073 C123 43 GPa 

In this work, we show that this fundamental assumption is no longer valid in optical nanofibers that carry both surface acoustic waves and hybrid acoustic waves.10 More specifically, we show that, when applying a tensile strain on a tapered optical fiber, the nonlinear elasticity of the nanofiber induces a strong elastic anisotropy which in turn affects the acoustic wave propagation and the Brillouin strain coefficients. As we will see, these dramatic changes in the mechanical properties of the nanofibers significantly impact the Brillouin strain coefficients that no longer scale with the Brillouin frequency shift νB.

Figure 1 shows the experimental setup for making, stretching, and characterizing the optical nanofibers using Brillouin spectroscopy. The ONFs were fabricated by heating and tapering uncoated standard single-mode fibers (SMFs) using the heat-brush technique.27 This technique uses a stabilized flame and two motion-controlled translation stages to heat and stretch the fibers. The nanofiber length and taper transition shapes were fully controlled by the computed trajectories of the two translation stages, while keeping the butane flame motionless. Using this technique, we were able to achieve ultrathin fiber waists down to 600 nm over a uniform length up to 80 mm. The transition regions were adiabatically drawn to ensure single-mode conversion from the SMF to the nanofiber with insertion loss down to 1.2 dB at 1550 nm. To probe the Brillouin backscattering signal, we used the heterodyne coherent detection technique shown in Fig. 1. This is an all-fiber setup where the light coming from a narrow linewidth continuous-wave laser at a wavelength of 1550 nm is split into two beams using a 90/10 fiber coupler. One beam is sent to an erbium-doped fiber amplifier (EDFA) and used as the pump wave and the other beam is used as a reference. The pump is injected through an optical circulator into the nanofiber under test, and the backscattered Brillouin signal from the nanofiber is then mixed with the reference using a 50/50 coupler. The two frequency-detuned beams give rise to an optical beat note that is detected in the radiofrequency (RF) domain by a fast photodiode. Note that the polarization state of the reference was adjusted with a polarization controller (not depicted in the experimental scheme) to allow for the best optical beat note visibility. The RF signal is then amplified, and the resulting Brillouin spectrum is recorded with an electrical spectrum analyzer (ESA) once the pulling process is ended and the flame heat has been dissipated. This prevents any pollution, stress, and temperature effects that could shift the Brillouin resonances and alter the measurement results. To further investigate the ONF nonlinear elastic properties, the Brillouin spectrum was measured by applying different tensile strains, still using the motion-controlled translation stages. An axial extensibility up to 6% was achieved before failure, a significant increase on standard SMFs that are usually limited to 2%.22 

FIG. 1.

Scheme of experimental setup for measuring the backward Brillouin spectrum of tapered optical fibers using heterodyne detection. Laser: coherent distributed-feedback (DFB) laser at 1550 nm, amplifier: 33 dBm EDFA, C: optical circulator, ESA: electrical spectrum analyzer. As the translation stages move apart, the tapered optical fiber is stretched by ΔL.

FIG. 1.

Scheme of experimental setup for measuring the backward Brillouin spectrum of tapered optical fibers using heterodyne detection. Laser: coherent distributed-feedback (DFB) laser at 1550 nm, amplifier: 33 dBm EDFA, C: optical circulator, ESA: electrical spectrum analyzer. As the translation stages move apart, the tapered optical fiber is stretched by ΔL.

Close modal

Figure 2(a) shows a typical experimental Brillouin spectrum in red, measured for an unstrained ONF with a diameter of 660 nm and a length of 80 mm. Note that these parameters were estimated using the recent technique described by Godet et al.27 The spectrum exhibits many frequency peaks in the frequency range 5–12 GHz, which are signature of all the Brillouin acoustic waves from both the tapered and untapered fiber sections. As shown in Ref. 10, these waves include the surface (SAWs) and hybrid acoustic waves (HAWs) in the uniform waist section, as well as the strong Brillouin resonance around 11 GHz because of the long untapered fiber sections. To get better insight, we have labeled in Fig. 2 the main acoustic modes from the ONF section only. They are denoted by TR(2,1), L(0,1), L(0,2), and L(0,3), respectively, where TR and L stand for torsoradial and longitudinal acoustic modes. As we will see later in the theoretical Sec. IV, the TR(2,1) and L(0,1) are waves with a strong surface character while the L(0,2) and L(0,3) are closer to bulk longitudinal waves. We note that, at this micrometer scale, the acoustic waves propagate as modes, and therefore, they result from the coherent superposition of both the longitudinal and the shear waves. Furthermore, these acoustic modes possess their own dispersion relation giving rise to such multipeaked Brillouin spectrum, compared to standard single-mode fibers that usually exhibit a single peak around 11 GHz, as that observed just above the L(0,3).28,29 It is also important to emphasize that the surface and hybrid acoustic waves travel at acoustic velocities lower than the pure longitudinal acoustic wave. It results in lower Brillouin frequency shifts going down to 5.5 GHz, the minimum frequency being bound by the Rayleigh limit.10 

FIG. 2.

(a) Experimental Brillouin spectrum (red curve) of the optical nanofiber of diameter dMF,1 = 660 nm and a length of 80 mm without axial strain. The main acoustic modes in the nanofiber section are indicated by black arrows. (b) Color map of the experimental Brillouin spectra measured in an optical nanofiber as a function of fiber tensile strain ΔL (vertical axis). (c) Brillouin frequency shifts for each acoustic mode as a function of axial strain in percentage.

FIG. 2.

(a) Experimental Brillouin spectrum (red curve) of the optical nanofiber of diameter dMF,1 = 660 nm and a length of 80 mm without axial strain. The main acoustic modes in the nanofiber section are indicated by black arrows. (b) Color map of the experimental Brillouin spectra measured in an optical nanofiber as a function of fiber tensile strain ΔL (vertical axis). (c) Brillouin frequency shifts for each acoustic mode as a function of axial strain in percentage.

Close modal

Once the ONF has been fabricated, the flame is blown out. They are then characterized at room temperature and stretched by moving apart the translation stages by a distance ΔL. The stretching strongly affects all the Brillouin acoustic resonances, shifting them toward higher frequencies. This is shown by the color map in Fig. 2(b) as a function of ΔL. The Brillouin strain coefficients for all acoustic modes were further extracted from the slopes using the method detailed in supplementary material, Sec. I. The resulting data reported in Fig. 2(c) exhibit several striking features. Perhaps the most striking one is that the acoustic resonances do not scale with the Brillouin frequency νB, as usually expected from Eq. (1). For instance, the acoustic mode L(0,2) has a steeper slope (350 MHz/%, red crosses) than the L(0,3) modes (250 MHz/%, orange crosses). Yet, what is most surprising is that the L(0,1) and TR(2,1) resonances cross each other for a strain larger than 5%, just before failure. We would rather expect increasing slopes of the lines from bottom to top in Fig. 2(c) but this is not the case. This proves that Eq. (1) is no longer valid for tapered fibers. To confirm these observations, we performed further experiments using other ONFs featuring different waists and lengths. Figures 3(a) and 3(b) show experimental results for a 930-nm diameter and 40-mm long nanofiber. Once again, the L(0,3) mode has a greater slope than the L(0,4), and the L(0,2) mode now features a strong nonlinear deviation vs the axial strain as shown in Fig. 3(c).

FIG. 3.

(a) Experimental Brillouin spectrum (red curve) of the optical nanofiber of diameter dMF,2 = 930 nm and a length of 40 mm without axial strain. The main acoustic modes in the nanofiber section are indicated by black arrows. (b) Color map of the experimental Brillouin spectra measured in an optical nanofiber as a function of fiber tensile strain ΔL (vertical axis). (c) Brillouin frequency shifts for each acoustic mode as a function of axial strain in percentage.

FIG. 3.

(a) Experimental Brillouin spectrum (red curve) of the optical nanofiber of diameter dMF,2 = 930 nm and a length of 40 mm without axial strain. The main acoustic modes in the nanofiber section are indicated by black arrows. (b) Color map of the experimental Brillouin spectra measured in an optical nanofiber as a function of fiber tensile strain ΔL (vertical axis). (c) Brillouin frequency shifts for each acoustic mode as a function of axial strain in percentage.

Close modal

We then repeated the same experiment using many taper waists ranging from 0.6 μm up to 1.9 μm. Figure 4 shows the normalized Brillouin strain coefficients measured for various acoustic resonances, defined as CN=1νBdνBdϵ. For comparison, the normalized strain coefficient of an SMF-28 is plotted as a horizontal black line (normalized coefficient, CSMF = 4.25). According to Eq. (1), it corresponds to a pure longitudinal acoustic wave while the other bottom black line corresponds to the strain coefficient for a pure shear wave.30 As can be seen, the coefficients of most nanofibers (colored crosses) are highly dispersed between these two limits and even beyond for taper diameters smaller than 750 nm. These experimental data obviously show the limits of the standard model. In Sec. IV, we will present a modified theory of tensile strain based on third-order nonlinear elastic constants that allows for predicting the Brillouin strain coefficients observed in nanofibers.

FIG. 4.

Normalized Brillouin strain coefficient CN of acoustic resonances for many fiber taper diameters ranging from 0.6 μm up to 1.9 μm (tensile strain <2%). CSMF is the coefficient for SMF-28 related to longitudinal waves. CS refers to the computed coefficient using the acoustic shear velocity.

FIG. 4.

Normalized Brillouin strain coefficient CN of acoustic resonances for many fiber taper diameters ranging from 0.6 μm up to 1.9 μm (tensile strain <2%). CSMF is the coefficient for SMF-28 related to longitudinal waves. CS refers to the computed coefficient using the acoustic shear velocity.

Close modal

There are four main reasons that can explain the unexpected Brillouin strain coefficients observed in Figs. 2–4. The first reason comes from the large extensibility of nanofibers. Tensile strain going up to 6% leads to dramatic changes of their dimensions that can significantly affects the Brillouin resonances. Next, the hybrid nature of the acoustic modes generated in ONFs gives rise to a strong nonlinear dispersion of the Brillouin resonances. The third reason relies on the nonlinear elasticity of silica which requires considering all the third-order elastic tensor coefficients (TOEC). This is fundamentally different from Eq. (1) in which only the longitudinal third-order elastic coefficient C111 was taken in to account. At last, when tensile strain is applied onto a material, it acquires an elastic anisotropy that strongly affects the acoustic wave propagation. To elucidate this question, let us first recall the basics of backward Brillouin scattering (BS) in optical fibers, which arises from the interaction between two counterpropagating and frequency-detuned optical waves, and a longitudinal elastic wave, as shown schematically in Fig. 5.31 This inelastic scattering is usually modeled using three coupled-amplitude equations for the pump wave (blue), the backward Stokes wave (red), and the elastic wave (purple).32 Assuming the latter as a plane wave, this 3-wave interaction satisfies a phase-matching condition that reads as βa = βpβS. Neglecting the very small dispersion between the pump and Stokes waves, as |βp| ≃ |βS|, the phase-matching relation can be approximated by |βa| ≃ |2βp|, leading to the Brillouin frequency shift νBvβp/π, where v is the velocity of the acoustic wave. It is worth emphasizing here that the acoustic velocity v represents the main contribution to the strain effect.21 

FIG. 5.

Scheme of backward Brillouin scattering in an optical nanofiber. Blue arrow: pump beam, red arrow: Stokes beam, and purple: acoustic wave. When applying tensile strain, the acoustic wave behavior dramatically changes because of the elastic anisotropy.

FIG. 5.

Scheme of backward Brillouin scattering in an optical nanofiber. Blue arrow: pump beam, red arrow: Stokes beam, and purple: acoustic wave. When applying tensile strain, the acoustic wave behavior dramatically changes because of the elastic anisotropy.

Close modal

Unlike single-mode optical fibers, optical nanofibers have an air clad and a very small core at the scale of both the optical and acoustical wavelengths. Hence, the elastic waves involved in Brillouin scattering can no longer be considered as plane waves, but as acoustic modes that include both longitudinal and shear elastic components. A direct consequence is that the acoustic velocity v must be computed for each acoustic mode. To do so, we assume the displacement field as the superposition of a longitudinal wave of velocity vL and a shear wave of velocity vs. Then, we calculate the stress field using the Hooke’s law. In air, all stress components applied to the boundary of the ONF must be zero. Those boundary conditions eventually lead to a dispersion equation whose solutions are the acoustic velocities. However, the computation of the acoustic velocities is here quite complex, since the tensile strain can reach up to a few percent, which is far beyond the linear elasticity regime.

If the medium is strained so that infinitesimal deformations are no longer considered, nonlinear acoustic effects occur and the third-order elastic constants must be taken into account. The nonlinear Hooke’s law can be written as29 

(2)

where σα is the α component of the stress tensor, ϵβ is the β component of the strain tensor, Cαβ is the (α, β) component second-order elastic constant, and Cαβγ the (α, β, γ) third-order elastic constant, with α, β, γ ∈ (rr, θθ, zz, rz, θz, ) ↔ (1, 2, 3, 4, 5, 6) in cylindrical coordinates. The nonlinear regime has to be modeled with much attention, as we now have to solve jointly the static and dynamic solutions. The static solution is the amount of strain ϵ¯zz induced by the tensile strain, while the dynamic solutions are the acoustic waves ϵ̃β that propagate along the strained optical microfibers. However, because of acoustic nonlinearity, the dynamic solutions are affected by the static solution at high tensile strain. In the following we will first solve the static problem by computing ϵ¯zz and we will further address the dynamic problem by computing the acoustic velocities for a given solution ϵ¯zz.

For the static problem, we can assume that only the longitudinal stress is significant, and therefore, we can write the nonlinear Hooke’s law as σ¯zz=c11ϵ¯zz+12c111ϵ¯zz2. Integrating the equation along the tapered fiber yields the expression of strain ϵ¯zz as a function of the total elongation9,33 (see supplementary material, Sec. I, for more details).

Before addressing the nonlinear dynamic problem, let us recall that ordinary linear elasticity and Hooke’s law describe the physics of springs. It states that stress is proportional to deformation via the second-order elastic constants. One of the direct consequences of elasticity is the ability of a vibration to propagate as a wave. The wave velocity is then related to the second-order elastic constants. In an isotropic medium, the twelve components of the elastic tensor result from a linear combination of two linearly independent coefficients, namely, the Young and shear moduli, or the two Lamé constants (λ, μ). As a result, in the linear regime, the second term of the right-hand side of Eq. (2) is negligible and the only term left is related to the stiffness matrix Cαβ. It is a symmetric tensor for an isotropic material, which can be written in Voigt notation as

(3)

where λ = 16 GPa and μ = 31 GPa are the Lamé constants of fused silica. The colors in the matrix indicate the symmetry. Solving the linear dynamic problem leads to two types of waves: the shear or longitudinal elastic waves.

In the nonlinear case, the tensile strain ϵ¯zz is strong enough to induce an acoustic nonlinearity. In the second term of the right-hand side of Eq. (2), all terms involving ϵ¯zz are no longer negligible. The Hooke’s law for the dynamic stress then reduces to

(4)

where

(5)

Next, to solve the dynamic problem, we need to consider the strain-induced elastic anisotropy. In an isotropic medium, the matrix coefficients come from a linear combination of three independent coefficients C111, C112, and C123, respectively. As previously discussed, only the tensile strain ϵ¯zz is strong enough to provide an effective acoustic nonlinearity that breaks the isotropic symmetry of the tensor. Therefore, Eq. (5) possesses only Cαβ3 third-order coefficients, and the effective elastic tensor can be written as

(6)

where the data have been calculated from Table I and ○ denotes the Hadamard product.16,34,35 The tensor matrix is now typical of a medium with transverse isotropic symmetry with respect to the z axis. The coefficients highlighted in blue and green give rise to a strong asymmetry in the other directions. For example, when the ONF is stretched at the maximum strain ϵ¯zz6%, the elastic coefficient C12 significantly increases by more than 90%.

Computing the acoustic mode velocity in a silica rod with transverse isotropic symmetry is somewhat complicated.36 In our case, however, the strong symmetry of our problem and the limited solutions we are interested in allow the problem to be simplified. For example, the longitudinal and shear waves approximation is still valid. The elastic coefficients related to shear waves are similarly affected by the acoustic nonlinearity, as shown by the lower-right submatrix coefficients, in which the pink and purple coefficients are quite close. Furthermore, although the upper-left submatrix diagonal coefficients are no longer the same, the phase-matching condition requires that the acoustic wave propagation occurs along the z direction. Therefore, the effect of the acoustic nonlinearity on C33=vz2ρ=(λ+2μ)1+C111λ+2μϵ¯zz is one of major sources of perturbation, where ρ is the material density. In Eq. (6), the coefficient affecting the longitudinal velocity along z is C111λ+2μ=6.71 (highlighted in blue). It increases the longitudinal velocity vz of a plane wave propagating along z up to 20% under high tensile strain. As a comparison, the longitudinal velocity for the other directions (vx and vy) increases by 9% only. The strong asymmetry of the elastic properties does not prevent to perform the standard computation of acoustic modes since the phase-matching condition for Brillouin backscattering concerns waves propagating essentially along the z direction. Given these observations, we assume that the displacement field is still the superposition of a pure longitudinal wave ϕ propagating along z with the velocity vz=C33ρ and of a shear wave Ψ of velocity vs=C44ρ, and therefore, it can read as

(7)

where in cylinder coordinates

with p=Ω2vz2k2 and q=Ω2vs2k2, n is the azimuthal order, Ω the acoustic wave pulsation, k is the acoustic longitudinal wavevector, A, B, and C are amplitudes constants.37 From Eq. (4), we then compute the stress component involved in the boundary conditions, and this yields for the Hooke’s law

(8)
(9)
(10)

The boundary conditions impose that these three components are zero at the surface of the nanofiber. This leads to a system of 3 equations with nonzero solutions if the determinant is null. Solving the resulting dispersion equation we find the acoustic modes wavevector βa from which we can get the Brillouin frequency shift. Comparing the coefficients highlighted in colors with their detailed expression in Eq. (6), we find that the acoustic nonlinearity has a contribution of about 3 times ϵ¯zz for most tensor coefficients, except one coefficient in the off diagonal green component C12=λ1+C112λϵ¯zz with C112λ=15.05. This coefficient alone has a big impact on the dispersion equation. As the boundary condition requires Eq. (8) to be null, it strongly affects how the z-longitudinal strain component ϵzz=uzz is coupled to the transverse strains ϵrr=urr and ϵθθ=1ruθθ+urr.

Once the dispersion equation has been solved, we further need to combine the static and dynamic solutions to get the acoustic wavevector and to compute the Brillouin frequency shifts as a function of tensile strain applied on the nanofiber (see supplementary material, Sec. II for detailed calculations). Figures 6(a) and 6(b) shows the results of these calculations as white dots and white diamonds superimposed on the experimental Brillouin spectra. The white dots describe the TR(2,1) Brillouin acoustic mode as a function of tensile strain while the white diamonds correspond to the L(0,x) longitudinal modes, with x = 1–3. As can be seen, the agreement between theory and experiment is excellent. In particular, all the Brillouin frequency shifts and the crossing of TR(2,1) and L(0,1) modes are very well fitted by our model. To further illustrate the influence of the induced anisotropy, the same computation was made without taking into account the asymmetry induced by C12 green component. The results are reported by the blue diamonds in Fig. 6(b). Comparing the white and blue diamonds shows that the anisotropy coming from C12 is the main reason for the mode crossing between L(0,1) and TR(2,1). This mode crossing can be easily explained by looking at Fig. 6(b) that shows the four main acoustic modes propagating in the nanofiber. Note that the deformation has been magnified for better visibility and that the color stands for the longitudinal deformation ϵzz. When comparing the four modes in Fig. 6(b), L(0,1) definitely shows the most hybrid nature. It has the strongest combination of both transverse deformation and longitudinal strain. In contrast, the TR(2,1) mode, highlighted in yellow, is almost a pure surface mode. L(0,2) and L(0,3), although being hybrids, are rather surface-like or longitudinal-like modes. The fact that C12 affects the coupling between the transverse deformations and the longitudinal strain explains why the L(0,1) velocity is so strongly altered.

FIG. 6.

(a) Comparison between our analytic model (white dots and diamonds) and the experimental Brillouin spectra in false color for a nanofiber of diameter dMF,1 = 660 nm vs tensile strain from 0% to 6%. The blue diamonds show the analytic model results with C112 = 0. When C112 is zero, the anisotropy vanishes and the theory does not predict TR(2,1), L(0,1) Brillouin resonances crossing anymore. (b) Experimental Brillouin spectrum of the fiber without strain. The lower insets show a simulation of both the transverse deformation and the longitudinal strain of the nanofiber for four acoustic modes including the TR(2,1), L(0,1), L(0,2), and L(0,3), respectively. The vertical scale is normalized to roughly one acoustic wavelength.

FIG. 6.

(a) Comparison between our analytic model (white dots and diamonds) and the experimental Brillouin spectra in false color for a nanofiber of diameter dMF,1 = 660 nm vs tensile strain from 0% to 6%. The blue diamonds show the analytic model results with C112 = 0. When C112 is zero, the anisotropy vanishes and the theory does not predict TR(2,1), L(0,1) Brillouin resonances crossing anymore. (b) Experimental Brillouin spectrum of the fiber without strain. The lower insets show a simulation of both the transverse deformation and the longitudinal strain of the nanofiber for four acoustic modes including the TR(2,1), L(0,1), L(0,2), and L(0,3), respectively. The vertical scale is normalized to roughly one acoustic wavelength.

Close modal

Although our analytic model provides a good agreement with experimental data for these four modes, it does not find the small frequency peak around 9 GHz in Fig. 6(b) and its weak slope under tensile strain. Once again, this behavior is in complete contradiction with the standard model described by Eq. (1) that scales with the Brillouin frequency. This is actually because these small peaks come from the two taper transitions and not from the nanofiber itself (see Fig. 5). However, the two adiabatic taper transitions experience less tensile strain than the nanofiber and thus their contribution must be separately computed. To that end, we used a finite-element method (FEM), which allows for including the electrostrictive force and for modeling the exact solutions.38 Attention must be paid, however, because the transition regions have wider and varying diameters. They support a great number of acoustic modes, most of which are not excited because of a weak transverse overlap with the optical modes. By including the electrostrictive force in the simulations, we rule out most of the nonexcited acoustic modes.

The numerical results reported in Fig. 7(a) show the Brillouin spectrum as a function of tensile strain for a nanofiber of 660 nm, including both the tapered and untapered fiber sections. The comparison with the experimental results shown in Fig. 7(b) is very good and even better than with the analytical model. The computed Brillouin spectrum now clearly shows some peaks in the 9–10 GHz range at zero strain because of the tapered fiber transitions, as those experimentally observed. Unlike the analytic model, the FEM does not make any assumption on the solutions and allows for predicting the Brillouin spectrum of nanofibers for every diameter. To get further insight, Fig. 7(c) shows the simulated Brillouin spectrum of the optical taper transition regions only. One can clearly see the lines starting from 9 to 10 GHz without tensile strain, as experimentally shown in Fig. 7(b).

FIG. 7.

(a) Numerical simulations of the Brillouin spectrum including the uniform nanofiber, the taper transitions (up to 3 μm diameter), and the single-mode fiber. (b) Experimental Brillouin spectra for a 660 nm-diameter tapered optical fiber vs strain. (c) Numerical simulations of the tapered transitions only for diameter up to 3 μm. The tensile strain for each part of the transition regions was computed using Ref. 9.

FIG. 7.

(a) Numerical simulations of the Brillouin spectrum including the uniform nanofiber, the taper transitions (up to 3 μm diameter), and the single-mode fiber. (b) Experimental Brillouin spectra for a 660 nm-diameter tapered optical fiber vs strain. (c) Numerical simulations of the tapered transitions only for diameter up to 3 μm. The tensile strain for each part of the transition regions was computed using Ref. 9.

Close modal

To confirm the model, comparison for different diameters was performed. In Figs. 8(a)–8(c), we report the measurement and simulation data for 660 nm and 930 nm diameters. The Brillouin resonances were extracted from the simulation data by finding the local maxima of the simulated Brillouin spectrum. The superposition of experimental data and the simulated Brillouin resonances (dots) shows a very good agreement for all diameters, as shown in Figs. 8(a) and 8(b). Furthermore, Fig. 8(b) shows that the Brillouin resonances for a 930 nm-diameter nanofiber nonlinearly shift with the tensile strain. This behavior is not due to the fourth-order elastic constant but to the complex acoustic dispersion. This nonlinear evolution is well predicted by both analytic and FEM models when TOECs are included. The nonlinear shift with respect to tensile strain is here due to the fact that the two branches around 8 and 9 GHz correspond to elastic waves strongly coupled at the nanofiber surface.27 This strong coupling leads to an avoided crossing that dramatically affects the Brillouin resonances.

FIG. 8.

Comparison between the experimental Brillouin spectra in color map vs finite-element method-based numerical simulations for two different diameter optical nanofibers: (a) 660 nm (b) 930 nm.

FIG. 8.

Comparison between the experimental Brillouin spectra in color map vs finite-element method-based numerical simulations for two different diameter optical nanofibers: (a) 660 nm (b) 930 nm.

Close modal

To conclude, we have reported a detailed investigation of backward Brillouin scattering in tapered optical nanofibers under tensile strain. Our results revealed that the fundamental properties of acoustic waves dramatically change because of elastic anisotropy. We observed, in particular, unexpected Brillouin strain coefficients ranging from 30% to 110% of the single-mode fiber’s value, followed by a nonlinear evolution at high tensile strength. We further provided a complete theoretical model based on third-order nonlinear elasticity of silica that agrees well with our experimental data. Our model is valid for a wide range of fiber taper diameters and it was further checked using FEM-based numerical simulations. The complex dispersion relation of acoustic modes is further responsible for the nonlinear behavior with respect to tensile strain. Finally, this work contributes to a further understanding of the complex light-sound interaction in optical waveguides and these new behaviors could open the way to the development of compact tensile strength optical sensors based on optical nanofibers.

The supplementary material includes two sections. The first section details the solving of the static problem and the computation of the inhomogeneous distribution of tensile strain along a tapered fiber. It provides the nanofiber static strain ϵ¯ for a tapered fiber elongation ΔL.

The second section details the combination of both static and dynamic strains which is not straightforward in a nonlinear acoustic problem. It shows how an effective acoustic wavelength is computed to find the acoustic wavevector.

This work was funded by the French National Research Agency under Grant Agreement Nos. ANR-16-CE24-0010-03, ANR-17-EURE-0002, and ANR-15-IDEX-0003. Adrien Godet is thankful to the Conseil Régional de Bourgogne Franche-Comté for student scholarship.

1.
T. A.
Birks
and
Y. W.
Li
,
J. Lightwave Technol.
10
,
432
(
1992
).
2.
L.
Tong
,
R. R.
Gattass
,
J. B.
Ashcom
,
S.
He
,
J.
Lou
,
M.
Shen
,
I.
Maxwell
, and
E.
Mazur
,
Nature
426
,
816
(
2003
).
3.
L.
Tong
and
M.
Sumetsky
, “
Subwavelength and nanometer diameter optical fibers
,” in
Advanced Topics in Science and Technology in China
(
Springer-Verlag
,
Berlin, Heidelberg
,
2010
).
5.
M. A.
Foster
,
A. C.
Turner
,
M.
Lipson
, and
A. L.
Gaeta
,
Opt. Express
16
,
1300
(
2008
).
6.
C.
Sayrin
,
C.
Clausen
,
B.
Albrecht
,
P.
Schneeweiss
, and
A.
Rauschenbeutel
,
Optica
2
,
353
(
2015
).
7.
B.
Gouraud
,
D.
Maxein
,
A.
Nicolas
,
O.
Morin
, and
J.
Laurat
,
Phys. Rev. Lett.
114
,
180503
(
2015
).
8.
E.
Meyer-Scott
,
A.
Dot
,
R.
Ahmad
,
L.
Li
,
M.
Rochette
, and
T.
Jennewein
,
Appl. Phys. Lett.
106
,
081111
(
2015
).
9.
S.
Holleis
,
T.
Hoinkes
,
C.
Wuttke
,
P.
Schneeweiss
, and
A.
Rauschenbeutel
,
Appl. Phys. Lett.
104
,
163109
(
2014
).
10.
J.-C.
Beugnot
,
S.
Lebrun
,
G.
Pauliat
,
H.
Maillotte
,
V.
Laude
, and
T.
Sylvestre
,
Nat. Commun.
5
,
5242
(
2014
).
11.
O.
Florez
,
P. F.
Jarschel
,
Y. A. V.
Espinel
,
C. M. B.
Cordeiro
,
T. P. M.
Alegre
,
G. S.
Wiederhecker
, and
P.
Dainese
,
Nat. Commun.
7
,
11759
(
2016
).
14.
L.
Brillouin
,
Les Tenseurs En Mécanique Et En Élasticité
(
Dover Publications
,
1946
).
15.
D. S.
Hughes
and
J. L.
Kelly
,
Phys. Rev.
92
,
1145
(
1953
).
16.
E. H.
Bogardus
,
J. Appl. Phys.
36
,
2504
(
1965
).
17.
Q.
Wang
,
G. A.
Saunders
,
H. B.
Senin
, and
E. F.
Lambson
,
J. Non-Cryst. Solids
143
,
65
(
1992
).
18.
19.
R.
Vacher
and
J.
Pelous
,
Phys. Lett. A
53
,
233
(
1975
).
20.
J.
Pelous
and
R.
Vacher
,
Solid State Commun.
18
,
657
(
1976
).
21.
T.
Horiguchi
,
T.
Kurashima
, and
M.
Tateda
,
IEEE Photonics Technol. Lett.
1
,
107
(
1989
).
22.
C. A.
Galindez-Jamioy
and
J. M.
López-Higuera
, “
Brillouin distributed fiber sensors: An overview and applications
,” https://www.hindawi.com/journals/js/2012/204121/abs/ (
2012
).
23.
D.
Culverhouse
,
F.
Farahi
,
C. N.
Pannell
, and
D. A.
Jackson
,
Electron. Lett.
25
,
913
(
1989
).
24.
T.
Horiguchi
,
T.
Kurashima
, and
M.
Tateda
,
IEEE Photonics Technol. Lett.
2
,
352
(
1990
).
25.
M.
Niklès
,
L.
Thévenaz
, and
P. A.
Robert
,
Opt. Lett.
21
,
758
(
1996
).
26.
M.
Guerette
,
C. R.
Kurkjian
,
S.
Semjonov
, and
L.
Huang
,
J. Am. Ceram. Soc.
99
,
841
(
2016
).
27.
A.
Godet
,
A.
Ndao
,
T.
Sylvestre
,
V.
Pecheur
,
S.
Lebrun
,
G.
Pauliat
,
J.-C.
Beugnot
, and
K.
Phan Huy
,
Optica
4
,
1232
(
2017
).
28.
N. L.
Rowell
and
G. I.
Stegeman
,
Phys. Rev. Lett.
41
,
970
(
1978
).
29.
A. D.
Volkov
,
A. I.
Kokshaiskii
,
A. I.
Korobov
, and
V. M.
Prokhorov
,
Acoust. Phys.
61
,
651
(
2015
).
30.
Y.
Tanaka
and
K.
Ogusu
,
IEEE Photonics Technol. Lett.
11
,
865
(
1999
).
31.
A.
Kobyakov
,
M.
Sauer
, and
D.
Chowdhury
,
Adv. Opt. Photonics
2
,
1
(
2010
).
32.
R. W.
Boyd
,
Nonlinear Optics
, 3rd ed. (
Academic Press
,
Amsterdam, Boston
,
2008
).
34.
B. E.
Powell
and
M. J.
Skove
,
J. Appl. Phys.
41
,
4913
(
1970
).
35.
W. T.
Yost
and
M. A.
Breazeale
,
J. Appl. Phys.
44
,
1909
(
1973
).
36.
F.
Honarvar
,
E.
Enjilela
,
A. N.
Sinclair
, and
S. A.
Mirnezami
,
Int. J. Solids Struct.
44
,
5236
(
2007
).
37.
D.
Royer
and
E.
Dieulesaint
,
Elastic Waves in Solids I: Free and Guided Propagation
(
Springer
,
2000
).
38.
J.-C.
Beugnot
and
V.
Laude
,
Phys. Rev. B
86
,
224304
(
2012
).

Supplementary Material