Optomechanical Self-Channeling of Light in a Suspended Planar Dual-Nanoweb Waveguide

  • A. Butsch, C. Conti, F. Biancalana, and P. St. J. Russell, Max Planck Institute for the Science of Light, Guenther-Scharowsky-Straße 1, 91058 Erlangen, Germany and Department of Physics and Institute for Complex Systems-CNR, Sapienza University, Piazzale A. Moro 2, 00185 Rome, Italy
Phys. Rev. Lett. 108, 093903
Sketch of the dual-nanoweb fiber structure. Light propagates in the z direction. The electric field of the TE mode and the magnetic field of the TM mode point parallel to the x axis.

Attractive and repulsive optical forces can appear between two coupled waveguides when light is launched into them [1–4]. Recently, experimental observations of transverse gradient forces and mechanical bending have been reported for a nanostructured waveguide beam [5], stacked ring microcavities [6], a periodically patterned “zipper” cavity [7], and coupled waveguide beams [8]. Optical forces can cause transverse deformations in coupled waveguide systems, which in turn change the effective refractive index, giving rise to the so-called “mechanical Kerr effect” [9,10].

In this Letter, we present a novel optomechanical phenomenon— the self-channeling of light, which we study theoretically in a system of two glass waveguides (“nanowebs”) suspended between the walls of a capillary fiber. When an optical mode is launched into the nanowebs, optical gradient forces cause them to bend inwards or outwards depending on the symmetry of the optical mode profile perpendicular to the nanoweb plane [1]. This mechanical deflection increases the effective refractive index of the mode, causing self-focusing of light within the nanoweb plane and altering the radiation-induced pressure distribution. We show that this can lead to stable self-channeling at power levels as low as a few mW in experimentally realistic structures [11]. The refractive index of the self-channeled mode turns out to be highly sensitive to small changes in the optical power, indeed the resulting differential optomechanical nonlinear coefficient can exceed the γ coefficient in conventional fibers by up to 7 orders of magnitude [12]. The magnitude of the optomechanical nonlinearity also remains high even when the nanowebs have quite dissimilar thicknesses.

Since a point force applied at one position causes a deflection that extends over the whole nanoweb, the optomechanical elastic response is highly nonlocal—a characteristic that is known to favor stable self-trapping [13]. Indeed, self-channeled modes in systems with a nonlocal response are an example of spatial “accessible” solitons [14]. A temporal analogue of such accessible solitons was recently proposed in [15].

Figure 1 shows the generic structure under study: the z axis points along the fiber axis and the y axis is perpendicular to the nanoweb plane. The webs have width L, thickness w and the air gap width is h. The structure supports both TE and TM polarized modes.

For the general case of an asymmetric dual waveguide we define a parameter s, so that the web thicknesses become w+s and w-s. In this case the dispersion relation can be written in the following implicit analytical form:

where p=k0ng2-neff2 and g=k0neff2-1 are the wave vector components in dielectric slabs with refractive index ng and in air, respectively, and A and B are given by A=2gpξ/(g2-p2ξ2) and B=(g2+p2ξ2)/(g2-p2ξ2). The parameter ξ=1 for TE and 1/ng2 for TM polarization, allowing both polarizations to be conveniently treated in a single analysis. For each value of integer m0 two modes exist, one [ + sign in Eq. (1)] corresponding to a mode with 2m field nodes and the other ( - sign) to a mode with 2m+1 nodes. For s=0 Eq. (1) reduces to expressions for the symmetric and antisymmetric modes of a dual-slab waveguide derived in [3].

Normalized Poynting vector distributions of selected TE self-channeled modes at a power of 100 mW. (a) Even m=0 mode and (b) odd m=0 mode for w=200nm, h=300nm, L=70μm, λ=800nm, and P=100mW. (c) Higher-order self-channeled TE mode with two lobes in the x direction (w=200nm). (d) Higher-order self-channeled TE mode with two lobes in each nanoweb in the y direction (w=400nm, h=300nm, λ=600nm). (e) Even and (f) odd m=0 mode for the asymmetric structure with w1=220nm and w2=180nm, the upper web being thicker.Normalized Poynting vector distributions of selected TE self-channeled modes at a power of 100 mW. (a) Even m=0 mode and (b) odd m=0 mode for w=200nm, h=300nm, L=70μm, λ=800nm, and P=100mW. (c) Higher-order self-channeled TE mode with two l... Show more

The electromagnetic force acting on each slab can be calculated by means of the Maxwell stress tensor. Since it is proportional to the intensity of light in the waveguide mode, the field amplitudes must be related to the optical power through the flux of the z component of the time-averaged Poynting vector S¯z. The total optical power P and the power density p per unit length in the x direction are calculated as follows:

The optical gradient pressure is obtained by evaluating the time-averaged y component of the Maxwell stress tensor (i=x,y,z):

The pressure-induced deflection δ(x) of a slab (thickness w), pinned rigidly at each edge ( x=±L/2) to a solid glass wall, is found using the standard plate deflection equation:

where E is Young’s modulus, ν Poisson’s ratio, and σ(x) the pressure distribution. The boundary conditions at x=±L/2 are δ=0 and dδ/dx=0. Since the mechanical deflection δ(x) varies very slowly across the breadth of the structure, the distribution of effective refractive index neff(x) can be accurately calculated using a local mode approach, i.e., solving Eq. (1) for neff(x) using the local value of h(x)=h+δu(x)+δl(x), δu and δl being the deflection profiles of the upper and lower webs ( δ<0 for attractive forces). Once neff(x) is known, the transverse distribution of the field [ e(x) or h(x) for TE or TM] can then be found by numerically solving the Helmholtz equation in the x coordinate:

which at the same time yields the modal index nm of the self-guided mode.

Making use of these connections between radiation-induced pressure, mechanical deformation, and transverse electromagnetic field distribution, an iterative cycle of numerical calculations can be established to seek self-consistent self-channeled solutions. Starting with a trial pressure profile σ0(x), we calculate the deflection δ(x) using (4). By solving (5) we obtain the mode index nm and the transverse field distribution. A new pressure profile σ1(x) is then calculated using (3) and the cycle restarted. Keeping the total optical power P constant, we iterate this procedure until it converges to a self-consistent steady-state solution. For the parameters used in this Letter, convergence is typically reached after less than ten iterations.

Power dependence of the modal index and the differential optomechanical nonlinearity. The modal indices converge to the same value at zero power. (a) Even m=0 TE mode in a structure with w=200nm, h=300nm, λ=800nm, and three different values of width L: 60μm (dashed line), 70μm (solid line), and 80μm (dash-dotted line). (b) is the same as (a) for the odd m=0 TE mode.Power dependence of the modal index and the differential optomechanical nonlinearity. The modal indices converge to the same value at zero power. (a) Even m=0 TE mode in a structure with w=200nm, h=300nm, λ=800nm, and three different values of ... Show more

We now consider a structure with w=200nm, h=300nm, and L=70μm, made from fused silica ( ng=1.45). Figures 2(a) and 2(b) show the self-channeling of the m=0 even and odd TE modes. The wavelength of the light is 800 nm, Young’s modulus 72.5kN/mm2, Poisson’s ratio 0.17, and the optical power was set to 100 mW. Self-channeling causes the nanowebs to be attracted for even modes and repelled for odd modes. The maximum deflection amplitude is 1.4 nm for the even mode (webs bent inward) and 2.5 nm for the odd mode (webs pushed outward). These relatively small deflections are sufficient to create a guiding index profile in the x direction.

In this context, it is necessary to estimate the amplitude of the vibrations induced thermally in a single nanoweb supported inside a fiber of length Lf. The number of vibrational modes of a rectangular plate [16] supported by a mechanical resonance of quality factor Q can be estimated as NLf/(LQ). Since kBThνvib at ambient temperature in these structures [typical resonant frequencies νvib are in the (sub-) MHz range], thermal energy can be expressed as Eth=kBTLf/(2LQ) where kB is Boltzmann’s constant. The deflection δmax at the center of the nanoweb, for simplicity subjected to a force F at its center (evenly distributed over the fiber length), is δmax=FL3(1-ν2)/(16LfEw3). Using this to calculate the stiffness of the nanoweb and its vibrational energy Evib=8Lfw3Eδ2/[L3(1-ν2)] and equating the result to Eth yields a thermal deflection amplitude:

for silica at 293 K. For Q=100 (a value measured in a similar structure [11]), L=70μm, and w=200nm this yields δth15pm, indicating (since we have shown above that deflections of order 1 nm occur at an optical power of 100 mW) that the self-channeling process might be disrupted by thermal vibrations at powers in the range of 1 mW and below.

Although optically induced heating of the material could in principle also cause self-focusing effects, extensive experiments on Kerr-related self-focusing in single nanoweb fibers, at similar drive powers, have shown no evidence that heating plays any role [17]. This is because fused silica has a very low absorption coefficient and a very low value of dn/dT ( <10-5K-1).

Note that when the air gap gets smaller, the index increases for even modes and reduces for odd modes. Thus both modes can be localized (channeled) in the central region of the waveguide. For the parameters chosen, the lateral extent of the m=0 TE mode is 20μm (even) and 16μm (odd). Higher-order self-channeled modes can also exist [Fig. 2(c)]. Changing the parameters to w=400nm, h=300nm, and λ=600nm, a self-channeled m=1 mode appears with two lobes across each nanoweb [Fig. 2(d)].

It is interesting to inquire how sensitive the self-channeling is to fabrication imperfections, such as unequal nanoweb thicknesses. In Figs. 2(e) and 2(f) the mode profiles are plotted for a structure with web thicknesses 220 and 180 nm. It can be seen that the self-channeling is robust against quite large degrees of asymmetry.

Contours of constant optomechanical nonlinear coefficient γom as a function of gap width h and wavelength (in units of 106W1km1). The nonlinear coefficient based on the Kerr effect is only 1W1km1. (a) Nonlinearity of the even m=0 TE mode of the dual-nanoweb waveguide with web thickness w=200nm, width L=70μm, and optical bias power P=300mW. (b) Nonlinearity of the odd m=0 TE mode, red curve showing the cutoff condition for the odd mode at given parameters.Contours of constant optomechanical nonlinear coefficient γom as a function of gap width h and wavelength (in units of 106W1km1). The nonlinear coefficient based on the Kerr effect is only 1W1km1. (a) Nonlinearity of the even m=0 TE mode o... Show more

The power dependence of the modal index for the m=0 even and odd self-channeled TE modes ( w=200nm, h=300nm) is presented in Fig. 3 for three different values of L. The lower limit of P in the calculations is determined by the need for the self-channeled mode to be confined tightly enough in the x direction so as not to be affected by the presence of the supporting capillary. By appropriate choice of parameters, powers as small as 10 mW can be sufficient to achieve optomechanical self-channeling in a realistic structure.

For both even and odd modes the modal index increases monotonically with power, the slope being steeper for wider, i.e., mechanically more compliant, structures. The differential change of modal index with variations in power can be regarded as a measure of nonlinearity of the system. We therefore define a differential optomechanical nonlinearity γom:

and plot its value as a function of the optical bias power in Fig. 3(a) for the even and Fig. 3(b) for the odd m=0 TE modes.

The intrinsic (i.e., electronic) γ coefficient can be calculated from the nonlinear refractive index of silica and the effective nonlinear mode area [18]. For the dual-web structure analyzed above it is 1W-1km-1, lying between the value typical of a small-core photonic crystal fiber (PCF) ( 0.2×103W-1km-1) and that of a hollow-core PCF ( 2×10-2W-1km-1) [19]. With peak values of order 108W-1km-1, γom is more than 10×106 times larger than the Kerr nonlinearity.

Although asymmetry causes dephasing and weakens the interaction between the waveguides, it turns out that even when it is strong (e.g., w1=250nm, w2=150nm) the nonlinear coefficient γom remains very high, dropping to roughly one third its value in the symmetric structure for the even and to two thirds for the odd m=0 TE mode, which means that fabrication tolerances for these structures can be quite relaxed.

The even and odd modes exhibit different types of self-channeling behavior. In Fig. 4 γom is plotted as a function of h and wavelength for L=70μm, w=200nm, and P=300mW. The maximum value of the nonlinearity for the given range of parameters reaches 109W-1km-1 for the even m=0 TE mode [Fig. 4(a)] and 5×108W-1km-1 for the odd mode [Fig. 4(b)]. Furthermore, the optical gradient force increases as the web separation decreases [1]. The m=0 odd mode, however, cuts off when h falls below a certain value, for example, at w=200nm and λ=1.3μm the odd TE mode disappears at h250nm. This cutoff spacing decreases with wavelength, as indicated by the red curve in Fig. 4(b). A structure in which only the even mode is self-channeled is potentially very interesting experimentally, since the system will become single mode, making measurements much more straightforward. Because of weaker confinement in the y direction, the odd mode shows maximum nonlinearity at longer wavelength compared to the even mode. Note that at shorter wavelengths the modes are more tightly confined in the y direction, causing the nonlinear coefficient to fall off more rapidly with increasing h.

The TM modes show qualitatively similar behavior, except that at very small spacings ( h<100nm) the m=0 even TM mode experiences higher radiation pressure than the even TE mode [3]. When the distance between the webs is decreased further, or alternatively the optical bias power is increased, the webs will touch (in case of even modes), causing a discontinuity in the optomechanical nonlinearity.

As mentioned above, the stability of the self-channeled modes and the very large differential nonlinear coefficient is a consequence of the extreme nonlocality of the mechanical response. Although the self-channeling nonlinearity is much stronger than the electronic Kerr nonlinearity of the glass, it is, of course, much slower. The time constant of the optomechanical nonlinear response is limited by the mechanical resonant frequencies of the dual-nanoweb [(sub-) MHz range, response time 110μs], which are, however, considerably faster than in other nonlocal media such as liquid crystals ( 1ms, [20]) and photorefractive materials ( 1s, [21]). In addition, the nonlinear response can be further enhanced by modulating the light at the acoustic resonant frequency, as demonstrated in [11].

The main purpose of this Letter has been to show for the first time that steady-state self-channeled modes can exist in this dual-nanoweb system. We now briefly indicate how the analysis can be extended to describe the dynamics of the self-channeling process. A paraxial model can be constructed for the propagation of an optical beam with intensity |A(x,z)|2:

where k=ωne/c, ne being the effective modal index of the undistorted nanoweb pair and Δn=neff(x)-ne the perturbation due to optomechanical forces. At equilibrium, for a given power P, the pressure-induced deformation in the vicinity of x=0 can be approximated by a parabolic function, leading to Δn(x)n0(P)-n2(P)x2. Explicit expressions for n0(P) and n2(P) can then be obtained from a second order expansion in x of the Green’s function of Eq. (4): n0(P)=(L/2)3 P/24, and n2(P)=LP/16 (further details will be reported elsewhere). The fact that Δn(x) is a function of power, and not of local field intensity, allows us to consider the optomechanical solitons as a unique realization of the so-called accessible solitons [14], which are self-trapped beams in the presence of a highly nonlocal response. A whole family of Hermite-Gaussian solutions, including fundamental and higher-order spatial solitons can be retrieved. A detailed study of the dynamics of these nonlocal self-trapped beams, and the investigation of other features, like collisions and modulational instability, will be reported in a future publication.

In conclusion, dual-nanoweb fibers are a promising vehicle for attaining ultrahigh nonlinear coefficients in optical fibers. In contrast to previously reported nonlocal nonlinear systems [20,21], they allow the formation of stable self-trapped modes that can propagate over many meters along a flexible pathway. This may permit the observation of spatiotemporal solitons (light bullets), and related phenomena such as supercontinuum generation inside spatially self-trapped beams.


  1. M. L. Povinelli et al., Opt. Lett. 30, 3042 (2005).
  2. A. Mizrahi and L. Schächter, Opt. Lett. 32, 692 (2007).
  3. F. Riboli et al., Eur. Phys. J. D 46, 157 (2007).
  4. W. H. P. Pernice et al., Opt. Express 17, 16 032 (2009).
  5. M. Li et al., Nature (London) 456, 480 (2008).
  6. G. S. Wiederhecker et al., Nature (London) 462, 633 (2009).
  7. M. Eichenfield et al., Nature (London) 459, 550 (2009).
  8. M. Li, W. H. Pernice and H. X. Tang, Nature Photon. 3, 464 (2009).
  9. W. H. P. Pernice, M. Li, and H. X. Tang, Appl. Phys. Lett. 95, 123507 (2009).
  10. J. Ma and M. L. Povinelli, Opt. Express 19, 10 102 (2011).
  11. A. Butsch et al., in Proceedings of the International FiO Conference San Jose, October 2011 (Optical Society of America, 2011), paper PDPC3.
  12. The γ coefficient is usually defined as the nonlinear phase change of a guided mode per W of optical power per km of propagation.

  13. O. Bang et al., Phys. Rev. E 66, 046619 (2002).
  14. A. W. Snyder and D. J. Mitchell, Science 276, 1538 (1997).
  15. C. Conti et al., Phys. Rev. Lett. 105, 263902 (2010).
  16. A. W. Leissa, J. Sound Vib. 31, 257 (1973).
  17. C. Kreuzer, A. Podlipensky, and P. St. J. Russell, Opt. Lett. 35, 2816 (2010).
  18. J. Laegsgaard, N. A. Mortensen, and A. Bjarklev, J. Opt. Soc. Am. B 20, 2046 (2003).
  19. P. St. J. Russell, J. Lightwave Technol. 24, 4729 (2006).
  20. C. Conti, M. Peccianti, and G. Assanto, Phys. Rev. Lett. 92, 113902 (2004).
  21. C. Rotschild et al., Phys. Rev. Lett. 95, 213904 (2005).

About the Authors

Image of A. Butsch
Image of C. Conti
Image of F. Biancalana
Image of P. St. J. Russell

Related Articles

Focus: Surface Waves Store Bouncing Droplet’s History
Nonlinear Dynamics

Focus: Surface Waves Store Bouncing Droplet’s History

A droplet bouncing and wandering across a liquid surface can produce waves that store the history of its chaotic motion. Read More »

Synopsis: Get the Hydrogen Out
Industrial Physics

Synopsis: Get the Hydrogen Out

Using clean helium made with a two-step purification process could prevent damaging blockages in cryostats. Read More »

Synopsis: Emptiness Constrains the Universe

Synopsis: Emptiness Constrains the Universe

The distribution of galaxies around regions of relatively empty space can be used to constrain cosmological parameters. Read More »

More Articles