Analysis of radiation pressure induced nonlinear optofluidics

We analyze two nonlinear optofluidic processes where nonlinearity is induced by the interplay between optical field and liquid interface. Specifically, guided optical waves generate radiation pressure on the liquid interface, which can in turn distort the liquid interface and modify the properties of the optical field. In the first example, we discuss the feasibility of nonlinear optofluidic solitons, where optical field is governed by the nonlinear Schrödinger equation and nonlinearity is effectively determined by liquid properties. Then, we analyze a nonlinear optofluidic process associated with a high quality (Q) factor whispering gallery mode (WGM) in a liquid droplet. Similar to Kerr effects, the WGM can produce a frequency shift proportional to the WGM power. Using liquid properties that are experimentally attainable, we find that it may only take a few photons to generate measurable WGM resonance shift. Such a possibility may eventually lead to nonlinear optics at single photon energy level. ©2014 Optical Society of America OCIS codes: (190.0190) Nonlinear optics; (190.3270) Kerr effect; (190.6135) Spatial solitons; (230.5750) Resonators. References and links 1. D. Psaltis, S. R. Quake, and C. Yang, “Developing optofluidic technology through the fusion of microfluidics and optics,” Nature 442(7101), 381–386 (2006). 2. U. Levy and R. Shamai, “Tunable optofluidic devices,” Microfluid Nanofluidics 4(1-2), 97–105 (2008). 3. C. Monat, P. Domachuk, and B. J. Eggleton, “Integrated optofluidics: A new river of light,” Nat. Photonics 1(2), 106–114 (2007). 4. S.-H. Kim, J.-H. Choi, S.-K. Lee, S.-H. Kim, S.-M. Yang, Y.-H. Lee, C. Seassal, P. Regrency, and P. Viktorovitch, “Optofluidic integration of a photonic crystal nanolaser,” Opt. Express 16(9), 6515–6527 (2008). 5. Z. Li and D. Psaltis, “Optofluidic dye lasers,” Microfluid Nanofluidics 4(1-2), 145–158 (2007). 6. H. Zhu, I. M. White, J. D. Suter, P. S. Dale, and X. Fan, “Analysis of biomolecule detection with optofluidic ring resonator sensors,” Opt. Express 15(15), 9139–9146 (2008). 7. H. Rokhsari, T. J. Kippenberg, T. Carmon, and K. J. Vahala, “Theoretical and experimental study of radiation pressure-induced mechanical oscillations (parametric instability) in optical microcavities,” IEEE J. Sel. Top. Quantum Electron. 12(1), 96–107 (2006). 8. T. J. Kippenberg and K. J. Vahala, “Cavity opto-mechanics,” Opt. Express 15(25), 17172–17205 (2007). 9. T. J. Kippenberg and K. J. Vahala, “Cavity optomechanics: Back-action at the mesoscale,” Science 321(5893), 1172–1176 (2008). 10. J. Hofer, A. Schliesser, and T. J. Kippenberg, “Cavity optomechanics with ultrahigh-Q crystalline microresonators,” Phys. Rev. A 82(3), 031804 (2010). 11. S. Tallur, S. Sridaran, and S. A. Bhave, “A monolithic radiation-pressure driven, low phase noise silicon nitride opto-mechanical oscillator,” Opt. Express 19(24), 24522–24529 (2011). 12. A. Cho, “Putting light’s light touch to work as optics meets mechanics,” Science 328(5980), 812–813 (2010). 13. A. Ashkin and J. M. Dziedzic, “Radiation pressure on a free liquid surface,” Phys. Rev. Lett. 30(4), 139–142 (1973). 14. I. I. Komissarovak, G. V. Ostrovskaya, and E. N. Shedova, “Light pressure induced deformations of a free liquid surface,” Opt. Commun. 66(1), 15–20 (1988). 15. J.-Z. Zhang and R. K. Chang, “Shape distortion of a single water droplet by laser-induced electrostriction,” Opt. Lett. 13(10), 916–918 (1988). 16. H. M. Lai, P. T. Leung, K. L. Poon, and K. Young, “Electrostrictive distortion of a micrometer-sized droplet by a laser pulse,” J. Opt. Soc. Am. B 6(12), 2430–2437 (1989). #223092 $15.00 USD Received 15 Sep 2014; revised 28 Oct 2014; accepted 3 Nov 2014; published 12 Nov 2014 (C) 2014 OSA 17 November 2014 | Vol. 22, No. 23 | DOI:10.1364/OE.22.028875 | OPTICS EXPRESS 28875 17. A. Casner and J.-P. Delville, “Adaptative lensing driven by the radiation pressure of a continuous-wave laser wave upon a near-critical liquid liquid interface,” Opt. Lett. 26(18), 1418–1420 (2001). 18. A. Casner, J.-P. Delville, and I. Brevik, “Asymmetric optical radiation pressure effects on liquid interfaces under intense illumination,” J. Opt. Soc. Am. B 20(11), 2355–2362 (2003). 19. J.-P. Delville, M. Robert de Saint Vincent, R. D. Schroll, H. Chraïbi, B. Issenmann, R. Wunenburger, D. Lasseux, W. W. Zhang, and E. Brasselet, “Laser microfluidics: fluid actuation by light,” J. Opt. A, Pure Appl. Opt. 11(3), 034015 (2009). 20. G. Bahl, K. H. Kim, W. Lee, J. Liu, X. Fan, and T. Carmon, “Brillouin cavity optomechanics with microfluidic devices,” Nat Commun 4, 1994 (2013). 21. M. Hossein-Zadeh and K. J. Vahala, “Fiber-taper coupling to Whispering-Gallery modes of fluidic resonators embedded in a liquid medium,” Opt. Express 14(22), 10800–10810 (2006). 22. A. Jonáš, Y. Karadag, M. Mestre, and A. Kiraz, “Probing of ultrahigh optical Q-factors of individual liquid microdroplets on superhydrophobic surfaces using tapered optical fiber waveguides,” J. Opt. Soc. Am. 29(12), 3240–3247 (2012). 23. A. Yariv and P. Yeh, Photonics: Optical Electronics in Modern Communications (Oxford University, 2007). 24. R. W. Boyd, Nonlinear Optics, 3rd ed. (Academic, Amsterdam, 2008). 25. J. D. Jackson, Classical Electrodynamics, 3rd ed. (John Wiley & Sons, 1998). 26. H. M. Lai, P. T. Leung, K. Young, P. W. Barber, and S. C. Hill, “Time-independent perturbation for leaking electromagnetic modes in open systems with application to resonances in microdroplets,” Phys. Rev. A 41(9), 5187–5198 (1990). 27. J. R. Buck and H. J. Kimble, “Optimal sizes of dielectric microspheres for cavity QED with strong coupling,” Phys. Rev. A 67(3), 033806 (2003). 28. A. Datta, S. Kundu, M. K. Sanyal, J. Daillant, D. Luzet, C. Blot, and B. Struth, “Dramatic enhancement of capillary wave fluctuations of a decorated water surface,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 71(4), 041604 (2005). 29. H. Leitão, A. M. Somoza, M. M. Telo da Gama, T. Sottmann, and R. Strey, “Scaling of the interfacial tension of microemulsions: A phenomenological description,” J. Chem. Phys. 105(7), 2875 (1996). 30. H. Chraibi, D. Lasseux, R. Wunenburger, E. Arquis, and J.-P. Delville, “Optohydrodynamics of soft fluid interfaces: optical and viscous nonlinear effects,” Eur Phys J E Soft Matter 32(1), 43–52 (2010).


Introduction
Optofluidics [1][2][3][4][5][6] and optomechanics [7][8][9] have recently emerged as two important areas of research.Optofluidics aims to synergize microfluidics and optics to achieve novel functionalities including reconfigurable optical systems [1,2], integrated optics [3,4], lasers [5], and sensing [6].Optomechanics, on the other hand, involves the dynamic interplay between optical field and mechanical motion [7][8][9].To date, most optomechanics-related research utilizes solid resonators [10][11][12].Only a limited number of studies investigated the mechanical interaction between optical field and fluids [13][14][15][16][17][18].The earliest example is perhaps the classic work in [13], where the authors used a focused high power laser beam to create a bulge over a flat air-liquid interface.Similarly, in [15,16], the authors used a high power laser beam focused onto a liquid droplet to distort its interface.More recently, optical radiation pressure was used in [17] to distort a flat liquid interface and form a tunable lens.By using a liquid mixture with extremely small interfacial tension, which can occur near the critical temperature of phase transition, it is possible to significantly reduce the optical power required for large interface distortion [17][18][19].Additionally, stimulated Brillouin scattering in a hollow capillary tube filled with liquid has been reported recently [20].The existence of high-Q whispering gallery mode (WGM) in an all-liquid droplet has also been experimentally confirmed [21,22].The focus of this work, however, is distinct from existing studies.Specifically, our goal is to demonstrate that under appropriate conditions, the dynamic interplay between optical force and liquid interface can lead to processes that are very similar to classical third-order nonlinear processes such as Kerr effects and optical solitons.As will be made clear in this paper, a defining feature of these nonlinear optofluidic processes is that the nonlinearity arises from the distortion of liquid interfaced induced by optical radiation pressure.
To analyze nonlinear optofluidic processes, it is necessary to solve the optical equations (e.g., Maxwell's equations) and the fluidic equations (e.g., Navier-Stokes equations) in a selfconsistent manner, where optical radiation pressure must be balanced by fluidic forces such as surface tension and buoyancy.Consequently, by tuning fluid parameters such as surface tension and density, one should be able to modulate and control the effective strength of optofluidic nonlinearity.This is in sharp contrast with traditional nonlinear optics, where nonlinear susceptibility (e.g., ( )   2 χ or ( ) 3 χ ) is an intrinsic material property and cannot be easily tuned.Furthermore, by using liquids with low surface tension, it should be possible to achieve large nonlinearity with low optical power.In fact, in this paper, we find that by reducing surface tension to a low but experimentally attainable level, the effective strength of nonlinear optofluidics can be several orders of magnitude stronger than the traditional Kerr effects.To the best of our knowledge, this possibility has never been discussed in existing literature.To quantitatively analyze nonlinear optofluidic effects, we consider two examples in this paper.Figure 1(a) shows a liquid structure which contains a high index liquid as the waveguiding core, air as the top cladding, and a lower index liquid as the bottom cladding.Assuming low interfacial tension for the bottom liquid-liquid interface, the presence of optical field should deform the interface and create a bulge.At the same time, this bulge can enhance optical confinement along the transverse direction.Under certain conditions, we find that the propagation of optical field within the liquid bulge can be described by the nonlinear Schrödinger equation, the solutions to which are in the familiar form of optical solitons.This example illustrates a key feature of nonlinear optofluidics, namely that the presence of the optical field distorts the liquid interface; simultaneously, the change in the liquid interface shape also modifies the optical field.As a result, analyses of nonlinear optofluidics often require solving the coupled system of optical waves and fluids.Specifically, the optical field must satisfy the Maxwell's equations plus the boundary conditions imposed by the deformable liquid interface.Meanwhile, the liquid system is governed by the relevant fluid equations, where the impact of optical radiation pressure must also be included.Clearly, the most important feature of nonlinear optofluidics is the coupling between the optical field and the liquids, where the coupling strength is effectively determined by the radiation pressure and the interfacial tension of the liquid interface.
The optofluidic soliton depicted in Fig. 1(a) may not be easy to implement experimentally.A more practical example is the system shown in Fig. 1(b), which is based on a liquid droplet that supports high-Q WGMs.In this example, the optical field of the WGM exerts radiation pressure on the droplet surface.Since the direction of radiation pressure always points from the high index core to the low index cladding [16], a bulge should form along the droplet equator.As a result, given sufficiently high optical power, both the effective optical path of the resonator and the corresponding resonance frequency should change.This phenomenon is very similar to the Kerr effect, where high optical power also shifts the effective cavity length by changing the refractive index of the liquid.In this paper, we derive a closed-form formula that can provide an order of magnitude estimate for this nonlinear optofluidic process.Using common liquid parameters, we find that this nonlinear optofluidic effect can be significantly stronger than the Kerr effect.In fact, for liquids with low but experimentally achievable surface tension, it may even be possible to produce measurable change in WGM resonance frequency at single photon energy level.Such a possibility may ultimately enable us to demonstrate nonlinear interaction between two single photons, which is obviously important for quantum information technology.The optofluidic solitons can exist in the asymmetric liquid waveguide shown in Fig. 2(a).The refractive indices of the three layers are represented by 1 n , 2 n , and 3 n , while their fluid densities are given by 1 ρ , 2 ρ , and 3 ρ , respectively.For simplicity, we assume layer 3 is air and 3 1 n = .In the absence of optical signals, the waveguide core is the planar middle layer with thickness 0 h and has the highest refractive index (i.e., 2

Optofluidic soliton
n n > and 2 3 n n > ).By using the water-in-oil microemulsion system described in [17], it is possible to ensure that the core liquid layer (layer 2) possesses a higher refractive index but a lower mass density than the cladding liquid layer (layer 1), i.e., 2 n n > and 2 1 ρ ρ < .Thus the configuration in Fig. 2 should be hydrodynamically stable for small deformation.
After coupling light into layer 2, the presence of the optical field generates radiation pressure on the dielectric boundaries.We assume that the interfacial tension between layer 1 and 2 is much smaller than that of layer 2 and 3.As a result, we can ignore the deformation of the air-liquid interface (the interface between layer 2 and 3) and focus on the liquid-liquid interface instead.We note that the direction of radiation pressure always points from the high refractive index medium towards the low refractive index medium [16].Consequently, the presence of optical field should generate a liquid bulge as illustrated in Fig. 2(a).In turn, the geometrical deformation can lead to better optical confinement within the bulge.Qualitatively, such a nonlinear optofluidic process is very similar to the self-focusing of spatial solitons.In the following analysis, we establish an analytical framework and confirm that under certain conditions, the aforementioned optofluidic process can indeed be described by the nonlinear Schrödinger equation.
Our first step is to apply the effective index theory [23], a widely used theory that allows us to reduce the original three-dimensional (3D) problem into a two-dimensional (2D) one.In our case, similar to the treatment of structures such as ridge waveguides, we assume that the optical field depends only on x and z and use an effective index ( ) , eff n x z to account for the field's variation along the y direction.For the specific example in Fig. 2(a), we assume that the optical field's distribution along the y axis is simply given by the fundamental y-polarized mode in a planar waveguide.Then, we approximate the x and z dependence of the optical field as: where ( ) Given the effective index in Eq. ( 2), the electric field in Eq. ( 1) should satisfy: ( ) Note that Eq. ( 3) no longer depends on y .Then, assuming a slow varying envelop ( , ) Equation ( 4) shows that the shape of the liquid bulge, i.e., ( ) , can significantly change the properties of the guided optical wave.
In order to determine ( ) , we must first derive its governing equation.For the electric field polarized along the y axis, the optical radiation pressure at the layer 1 and 2 interface is given by . Since we look for a solution that is self-guided and propagates along the x direction, its intensity, i.e., ( ) 2 , A x z , should not depend on x .Therefore, the shape of the bulge ( ) should not depend on x either.
Then, in equilibrium, interfacial tension on a deformed interface is balanced with buoyancy ), this balance leads to: where ρ ρ ρ Δ = − , and σ is the interfacial tension between layer 1 and 2. In Eq. ( 5), we ignore the deformation of the air-liquid interface, since the air-liquid surface tension is assumed to be much higher than the liquid-liquid interfacial tension.To obtain a closed form solution, we further assume that the interfacial tension term in Eq. ( 5) is even smaller than gravity and is therefore neglected.Under this scenario, the height of the liquid bulge is proportional to optical intensity, i.e., ( ) Substituting this result into Eq.( 4), we obtain: where .
Equation (6a) takes the familiar form of the nonlinear Schrödinger equation [24].Its fundamental solution is: Our assumptions in deriving the form of optofluidic solitons are either the standard approach or can be justified using carefully controlled experimental conditions.For example, the slow varying envelope assumption and the effective index theory are commonly used in nonlinear optics and integrated optics.Our assumptions on liquid properties such as density, refractive index, and surface tension, however, warrant additional discussion.For example, two key requirements are: 1) in comparison with the cladding medium (layer 3), the core medium (layer 2) should possess lower density and higher refractive index; 2) the surface tension between layer 1 and 2 should be small.Both requirements can be satisfied by using a water-in-oil microemulsion system that contains a mixture of water, sodium dodecyl sulfate, toluene, and n-butanol-1 [17].At a temperature slightly higher than the critical temperature ( c T = 35 °C) for phase transition, the mixture phase separates into two different micellar phases, with interfacial tension given by ( ) , where  , the Bond number is in the order of O (10).Such a large Bond number indicates that the interfacial tension effect is negligible compared to buoyancy effect.When the Bond number is small, the interfacial tension effect becomes comparable or even stronger than gravity, then we must solve the two coupled equations (Eqs.( 4) and ( 5)) without any approximations.For such cases, analytical solutions would unfortunately be unattainable and numerical methods must be applied.
An attractive feature of the optofluidic soliton is that its nonlinear behaviors can be tuned by varying the parameters of the liquid systems.In particular, an interface with weak interfacial tension should possess high nonlinear efficiency.As an example, consider the waveguide structure analyzed in Fig. 2 10 / V − m [24].Therefore, by using liquid systems with weak interfacial tension, we can effectively increase the third order nonlinearity by several orders of magnitude.
We can use Eq. ( 7) to estimate the optical power required to form the aforementioned optofluidic soliton.For this calculation, we assume 0 10 , and use .Using parameters given above, we find that the optical power of this optofluidic soliton is only 41 µW .We can also use the value for 0 A to estimate the bulge height.Applying the discussion immediately above Eq.( 6), we can relate the magnitude of bulge height to the optical field amplitude as The discussion in this section is mainly theoretical, where our primary aim is to show that radiation pressure induced interface deformation can lead to nonlinear processes very similar to traditional third order nonlinear processes such as optical solitons.In terms of experiment, it is perhaps easier to investigate nonlinear optofluidic processes using a high-Q optical resonator based on a liquid droplet, as will be discussed in the next section.12).The dimensionless constant σ Γ is extracted using Eq. ( 14) and least square fitting (dashed red line).

WGM induced droplet deformation
Figure 3(a) illustrates the second example of radiation pressure induced nonlinear optofluidic processes, where the key element is the high-Q WGM circulating along the equator of a liquid droplet in air.Since the direction of the radiation pressure always points from the high refractive index material to the low index material, the presence of a high-Q WGM would push the spherical droplet surface outwards and consequently enlarge the circumference of the droplet's equator.Since radiation pressure is proportional to the optical power carried by the WGM, under the limit of small droplet deformation, we expect that the equator circumference should increase linearly as a function of WGM power and shift the WGM resonance frequency as a direct result.This phenomenon is very similar to the Kerr effect.In the following analysis, we establish an analytical framework and provide an order of magnitude estimate for the frequency shift associated with this nonlinear optofluidic process.Our analysis suggests that the effective strength of the aforementioned nonlinear optofluidic process can be several orders of magnitude larger than the traditional Kerr effect.The starting point of our analysis is the Young-Laplace equation.We assume a small Bond number, i.e., buoyant force is much smaller than surface-tension force and is ignored.As a result, we find: where p Δ represents the pressure difference inside and outside of the droplet, opt P is the radiation pressure generated by the WGM, σ represents the surface tension of the droplet, and κ is the mean curvature of the droplet.Since our main purpose is to provide an order of magnitude estimate, we take the following steps to simplify our analysis.First, we restrict our considerations to the fundamental transverse electric (TE) mode.Mathematically, this means that the optical field of the WGM can be written as [25]: where the angular mode number l and m are integers and satisfy 0 l > and l m l − ≤ ≤ , ω is angular frequency.The subscript "q" is either "co" if it represents core parameters of the WGM resonator or "cl" if it denotes cladding parameters.For example, co k refers to the optical wave vector within the droplet, and cl k is the cladding wave vector.Similarly, co Z and cl Z represent the impedance of the droplet core and cladding material, respectively.In Eq. ( 9), lm X  is the vector spherical harmonic function [25] and represents the angular variation of the WGM.As to the radial profile of the WGM, the function ( ) l q g k r is either the l-th order spherical Bessel function (for core) or the spherical Hankel function of the first kind (for cladding), i.e.,: where a is the radius of the un-deformed sphere.For the fundamental WGM, the m l = mode and the m l = − mode are identical except for the direction of light circulation.For convenience, we only consider the m l = mode and denote it as | ll .Given Eq. (9a), the WGM radiation pressure opt P on the droplet interface is given by [16]: ( ) where 0 ε is free space permittivity, co n and cl n represent core and cladding refractive indices, and surf E  is the electric field on the droplet surface.Note that the expression for the optical force in Eq. ( 11) is slightly different from the one we used in Eq. ( 5), which is due to the difference in electric field direction [16].From the expression in Eq. (9a), we can verify Numerically solving Eq. ( 9) to ( 11) can give us the form of droplet deformation.However, the mathematical process is complex and can obscure the physics.In this paper, we emphasize the physical reality and utilize several assumptions to obtain an intuitive estimate of droplet deformation.First, we assume that the droplet deformation can be approximated by an oblate spheroid, as shown in Fig. 3(a).This assumption simplifies our analysis considerably.For example, assuming a spheroid shape, the mean curvature of the droplet interface is simply: where e x and p x are defined as ( ) assumption, the WGM resonance shift can be readily estimated using the perturbation theory in [26].The next step is to expand both sides of Eq. ( 8) using spherical harmonic function amount given in Eq. (13a).In other words, the reduction in fluid pressure is a constant term that balances the average magnitude of the radiation force (integrated over the entire droplet surface).After subtracting this constant component (i.e., p Δ ), it is the θ dependence of the radiation pressure that determines droplet deformation.In applying Eq. (13b), we essentially use the first non-zero alternating component to estimate the magnitude of droplet deformation.Obviously, to satisfy Eq. ( 8) exactly, the mean curvature must contain higher "frequency" components.However, using the curvature to calculate droplet shape requires double integration, which can be regarded as a low-pass filter that significantly suppresses these higher "frequency" components.Therefore, it is reasonable to discard the higher order spherical harmonic terms ( 4 L ≥ ), and use Eq. ( 13) to estimate the magnitude of droplet deformation.
To further simplify Eq. (13b), we note that the right hand side of Eq. ( 13b) is a function of .
Through least square fitting, we determine the value of Γ σ to be 1.01.
To obtain a closed form formula, we express 2 surf E  in Eq. ( 11) as: , Equation ( 18) relates droplet deformation to the peak electric field intensity on the droplet interface and is the key result of this paper.In the next section, we consider several examples of liquid droplets and calculate their deformation as a function of optical power associated with the circulating WGM.

Estimate of nonlinear optofluidic effects based on high-Q WGMs
Using Eq. ( 18), we can calculate the radiation pressure induced WGM frequency shift and quantitatively compare nonlinear optofluidics with the traditional Kerr effect.We choose the following parameters in our calculations.The core index is 1.
The droplet radius varies from 10 m μ to 400 m μ .All WGMs are the fundamental TE mode | ll , and their angular momenta l are chosen to ensure that the WGM wavelengths λ are in the vicinity of 1.56 m μ .Table 1 lists the l and λ of the TE | ll mode in such droplets with different radii a .The resonance frequencies are obtained by matching transverse field components across the droplet interface [27], which is assumed to be a perfect sphere.Equation ( 18) links droplet deformation with peak surface field intensity, which is unfortunately not directly measurable.Therefore, we need to further relate the peak field intensity in Eq. ( 18) with the optical power or energy of the circulating WGM.To accomplish this, we take advantage of the fact that the deformation is linearly proportional to .By integrating the Poynting vector across the 0 φ = plane, we can calculate the total optical power carried by the WGM.Similarly, through volume integration of energy density, we can obtain the total energy stored within the WGM.Then, we multiply the electric field intensity with an appropriate normalization factor such that the total power or the energy of the WGM becomes the desired value.This normalization factor is then used to calculate 2 peak surf E associated with the desired WGM power or energy.Once 2 peak surf E is obtained, the only unknown parameter in Eq. ( 18) is the dimensionless factor lm Γ θ , which can be obtained through numerical integration of Eq. ( 17) and is shown in Fig. 4(d).Using the procedure described above, we set the power associated with the circulating WGM to be 1 W and use Eq. ( 18) to calculate the radiation induced interface deformation for different droplet radius.The results are shown in Fig. 5(a).The red circles give the normalized droplet deformation ( ΔR / a ) induced by the WGM radiation pressure, where we assume the surface tension is that of water ( 72 σ = mN / m).The blue and black diamonds represent the estimated Kerr effect in water (blue) and CS 2 (black), produced by the same peak surface field . These estimates are simply obtained by using ( ) where ( ) χ is the third order nonlinear susceptibility of water ( ).Since both ΔR / a and Δn represent the same physical effects, i.e., the relative change in the optical path of the high-Q resonator, they can be shown in the same figure for direct comparison.Clearly, the nonlinear optofluidic effect is three to five orders of magnitude stronger than the Kerr effect.Since the nonlinear optofluidic process can be significantly stronger than the Kerr effect, it is worth asking the question: Is it possible to use only a few photons to generate measurable WGM frequency shift?To answer this question, we first calculate the value of that the total WGM energy equals to single photon energy ω  .Then, we substitute this value into Eq.( 18) and calculate the normalized radius change ΔR / a induced by the presence of a single photon.Finally, using the perturbation theory results in [26], the WGM frequency shift associated with the deformed spheroid is: Combining Eqs. ( 18) and ( 19), we can now calculate the radiation pressure induced frequency shift.The results are shown in Fig. 5 Given the fact that WGMs with Q factors as high as 2.3 × 10 6 have been observed using a liquid resonator [22], Fig. 5(b) suggests that one should be able to use only a few photons to produce experimentally observable resonance shift.If one can reduce the surface tension further down to the level of 0.1 / mN m σ = , then the presence of even a single photon can significantly change the characteristics of the liquid droplet resonator.
It is worth remarking that liquid systems with ultralow surface tensions have been experimentally attained using surfactants [28][29][30].For example, by introducing a bimolecular layer of preformed ferric stearate, one can reduce the surface tension of water down to 1 / mN m [28].In emulsion or microemulsion systems containing water (a polar fluid) and oil (a nonpolar fluid), surfactants can lower surface tension as small as 1 / N m μ at an optimal concentration [29].Therefore, for these liquid systems with ultralow surface tension, nonlinear optics at single photon energy level should be experimentally feasible.
In our WGM analyses, the effect of gravity is ignored.To justify this choice, we can estimate the Bond number of the liquid droplet.For a sphere with radius a , the gravitational and the surface tension effects can be estimated as ).Such a resonator is similar to the one investigated in [21].As a specific example, we estimate the deformation induced by 1W WGM power circulating within an oil-in-water droplet with radius .This result suggests that the conclusions we obtained using water-in-air droplets are broadly applicable to other types of liquid resonators.
Equation (18) suggests that droplet deformation depends on interfacial tension but not on liquid viscosity.This is to be expected, since our discussion is based on static analysis, i.e., the Young-Laplace equation in Eq. ( 8), which does not involve viscosity.We do, however, expect that viscosity should play a role in the dynamics of droplet deformation.
Finally, we point out that both the spherical WGM system analyzed here and the planar soliton case discussed in section 2 can be connected to liquid-based spheroid resonators.On the one hand, the spherical droplet analyzed in section 3 is a special case of the liquid spheroid resonators.On the other hand, a prolate spheroid resonator with an extremely large major axis can essentially be regarded as an infinitely long dielectric cylinder, which can also support WGM that circulates along the perimeter of the cylinder.And if we "unwrap" the infinitely long cylinder, the WGM traveling along the cylinder perimeter becomes similar to the self-guided soliton waves analyzed in section 2.

Summary
In this paper we analyze the possibility of using all liquid systems to achieve nonlinear optical effects with extremely low power threshold.The defining feature of the proposed nonlinear optofluidic processes is the interface deformation due to optical radiation pressure.In particular, we find that through the formation of a bulged interface, optical waves could selffocus within a liquid slab waveguide.Under appropriate conditions, these self-guided optical waves are governed by the nonlinear Schrödinger equation and can be described by the familiar soliton solutions.We then consider a spherical liquid droplet that supports high-Q WGMs.With sufficiently high power, the optical force associated with the circulating WGM could deform the liquid droplet and induce a frequency shift of the WGM resonance.By applying spherical harmonic expansion, we estimate that the radiation pressure induced nonlinearity is several orders of magnitude stronger than the traditional Kerr effect.In fact, we find that by using a liquid system with low but experimentally achievable surface tension, it is possible to produce measurable frequency shift at the energy level of a few photons.Such effect may ultimately lead to nonlinear interactions between two single photons.

Fig. 1 .
Fig. 1.(a) An optofluidic soliton formed from a self-guided optical wave that is confined within a liquid bulge formed through radiation pressure.(b) A liquid droplet that contains a high-Q WGM circulating along the equator.The radiation pressure of the WGM forms the bulge, which in turn shifts the WGM resonance frequency.

Fig. 2 .
Fig. 2. Illustration of an optofluidic soliton.The structure contains two liquids (refractive indices 1 n and 2 n ) and air ( 3 1 n = ).The thickness of the waveguide is 0 h in the absence of optical field.The radiation pressure of the guided optical signal produces the bulge shown in the figure.The thickness of the bulge that serves as the waveguide core is denoted as ( ) ( ) 0 h x h h x = +Δ .(b) The effective index of the asymmetric dielectric waveguide defined in (a) as a function of the core layer thickness ( ) h x .The waveguide parameters are 1 1.5 n = , v = .The two micellar phases possess different densities and refractive indices, with the density difference given by .The relative magnitude of buoyancy to interfacial paragraph above.Then, based on Eq. (7b)assume that the optical beam is confined in a cross-section of width 10 W µm = (i.e., soliton width) and height 3 H µm = (i.e., waveguide thickness), the optical power of such a beam can be simply estimated as

Fig. 3 .
Fig. 3. (a) A liquid droplet with a high-Q WGM circulating near its equator.The radiation pressure of the WGM deforms the original spherical droplet (the blue circle) and generates the bulge (the solid black line), which is approximated as an oblate spheroid (the dashed black line).The normalized equator radius e x is defined as the ratio of the spheroid radius at the equator ( ) a R + Δ and the radius of the original sphere a .(b) The integral ( ) ( ) 20 0 ( ) , sin e e F x a x Y d π κ θ θ θ θ = −  opt P , depends on θ but not on φ .Additionally, for the | ll mode,2 surf E is symmetric with respect to the equator plane (i.e., .) Furthermore, under the spheroid #223092 -$15.00USD Received 15 Sep 2014; revised 28 Oct 2014; accepted 3 Nov 2014; published 12 Nov 2014 (C) 2014 OSA 17 November 2014 | Vol.22, No. 23 | DOI:10.1364/OE.22.028875| OPTICS EXPRESS 28883

ex
only.Multiplying the integral by droplet radius a , we can define a dimensionless parameter and numerically evaluate it as a function of e x .The result, which is shown in Fig. 3(b), suggests that this integral is almost a perfect linear function of ex .Thus we introduce a dimensionless parameter Γ σ as:

Fig. 4 . 2 E 2 E 2 E
Fig. 4. (a) The radial dependence of 2 E  of a fundamental TE mode ( l = 257) in a spherical

Fig. 5 .
Fig. 5. (a) Radiation pressure induced droplet deformation ( / R a Δ ) in droplets with different radii (red circles).The total power of the WGM that circulates along the droplet equator is fixed at 1 W.For comparison, the changes in refractive index ( n Δ ) due to the Kerr effect are shown in the same figure, which are estimated using

1 W WGM power should induce a relative radius change of 4 / 4 /
procedure described above, we find that with the results shown in Fig.5(a), for an water-in-air droplet with identical radius, interfacial tension and WGM power, the deformation is

Table 1 .
The lm Γ θ values (blue circles) are obtained numerically using Eq.(17) and simply connected together using the dashed line.