We introduce tunable optofluidic microlasers based on active optical resonant cavities formed by optically stretched, dye-doped emulsion droplets confined in a dual-beam optical trap. To achieve tunable dye lasing, optically pumped droplets of oil dispersed in water are stretched by light in the dual-beam trap. Subsequently, resonant path lengths of whispering gallery modes (WGMs) propagating in the droplet are modified, leading to shifts in the microlaser emission wavelengths. Using this technique, we present all-optical, almost reversible spectral tuning of the lasing WGMs and show that the direction of tuning depends on the position of the pump beam focus on the droplet. In addition, we study the effects of temperature changes on the spectral position of lasing WGMs and demonstrate that droplet heating leads to red-tuning of the droplet lasing wavelength.
© 2013 Optical Society of America
Liquid droplets have been long recognized and studied for their peculiar optical properties . At the microscale, shape of liquids is dictated by the liquid surface tension that tends to minimize liquid surface area, thus transforming a small liquid parcel into a perfect sphere. Such spherical droplets with exceptionally smooth surfaces can act as optical cavities hosting morphology-dependent resonances called whispering gallery modes (WGMs) that are confined to small volumes near the droplet rim . Consequently, very high electric field intensities can be achieved in microdroplets that can mediate strong coupling between light and quantum systems such as atoms, molecules, or quantum dots .
Droplet-based optical resonators have found a number of applications in optofluidics - an exciting new research field that exploits unique properties of fluids for creating optical components and systems . Dye lasing , stimulated Raman scattering , and Raman lasing  in microdroplets have been demonstrated at low threshold pump powers. Apart from their signifi-cance for the fundamental physics research, microdroplets are well suited to another important line of applications that benefit from the sensitivity of the WGMs to the microdroplet size as well as the refractive indices of the droplet liquid and the surrounding medium. These properties naturally make microdroplets very attractive in biological and chemical sensing, and in developing optical devices that necessitate a tunable optical microcavity, e.g. tunable optical switches or light sources .
Tunability of shape, size, and composition is an important advantage of reconfigurable optofluidic components based on microdroplets. The most straightforward way of tuning spectral position of WGMs of airborne droplets relies on controlled changes of the liquid cavity size using evaporation/condensation , electrowetting , or photothermal effect . Adjustment of resonant frequency of microdroplet optical cavity can also be achieved by controlled deformation of the microdroplet which lifts the degeneracy of WGMs with different azimuthal mode numbers. Frequency splitting of the azimuthal WGMs then depends on the extent of the microdroplet deformation [11, 12]. For microdroplets surrounded by a solid elastomer, tuning has been achieved by mechanical shape deformation  or by changing refractive index of liquid crystal droplets in external electric field . However, for the case of microdroplets surrounded by another liquid in microfluidic channels, reversible spectral tuning has not been shown yet, despite different demonstrations of optical functionality of emulsion droplets including dye lasing [15, 16]. In order to address the issue of spectral tuning, researchers have demonstrated multicolor dye lasing in microfluidic channels by fast switching between consecutive droplets doped with different dye molecules . Recently, emission wavelength tuning by continuous decrease of the droplet size due to dissolving of benzyl alcohol microdroplets in water has also been demonstrated . Nevertheless, such dissolution-based tuning of the lasing emission wavelength is irreversible.
As shown by Guck et al., transparent dielectric microobjects suspended in a liquid can be deformed by the forces of light . At the object surface where the material properties (in particular, the dielectric constant) change abruptly, transfer of momentum from light to the object occurs. Such surface momentum transfer then leads to surface stresses that can cause deformation of the illuminated object, with the amount of deformation directly proportional to the incident light power. Optically induced deformation of soft microobjects (lipid vesicles, living cells) has been successfully demonstrated using a pair of counter-propagating laser beams in the so-called optical stretcher geometry . Optical stretcher uses the dual-beam optical trap configuration [20–22] which consists of two precisely aligned, slightly divergent counter-propagating laser beams emitted, for example, from optical fibers. Because of the Gaussian transversal intensity profile of both beams, high refractive index objects (e.g. oil droplets suspended in water) are pulled toward the beams axis and, subsequently, move along this axis to the location where the radiation pressures of both beams just balance each other. At this position, the objects are stably confined and their stretching in the optical field of the two beams can be studied.
Optical stretcher represents a convenient, all-optical method for tuning the emission spectra of lasing microdroplets via controlled droplet deformation that leads to frequency-splitting of degenerate WGMs of originally spherical droplets. However, the magnitude of optical stretching forces is rather small in comparison with the typical interfacial tension forces acting between two immiscible liquids in emulsions. In order to achieve sufficient droplet deformations without using prohibitively high stretching laser powers, interfacial tension has to be lowered by several orders of magnitude. This can be accomplished by using various mixtures of surfactant molecules and inorganic salts for the emulsion preparation. Emulsions with ultra-low interfacial tension have been used by Møller et al.  and Ward et al.  in the single beam and multiple beam optical tweezers, respectively, to observe large deformations of the emulsion droplets with moderate trapping powers. Despite the above demonstrations, optical stretcher geometry has not been employed yet for light-induced tuning of the emission spectra of lasing microdroplets.
In this paper, we introduce novel tunable optofluidic dye microlasers that use optical stretching forces for the adjustment of the laser emission wavelength. Our microlasers are based on droplets of immersion oil emulsified in water and confined in a dual-beam optical trap. Since the droplets are stained with a fluorescent dye, they function as active ultrahigh-Q optical resonant cavities that are optically pumped with an independent pulsed green laser. In order to reduce the stretching power required for microlaser tuning, we lower the interfacial tension between the droplets and aqueous host medium to ∼1 mN/m using a surfactant and sodium chloride. We demonstrate almost reversible spectral tuning of the lasing WGMs with tuning slopes as large as ∼0.5 nm/W (nm per Watt of the total stretching power). We also show that the direction of wavelength tuning depends on the position of the pump beam focus on the droplet. In addition, we characterize the effects of changes of the droplet and immersion medium temperature on the spectral position of lasing WGMs and show that these thermally induced spectral shifts are caused by simultaneous changes of the droplet size and refractive index of both the droplet and the host medium resulting from the droplet heating by the stretching laser beam.
2. Experimental setup for microdroplet lasing and stretching
Experimental setup for the observation of tunable lasing from optically trapped and stretched emulsion microdroplets is shown in Fig. 1. Studied oil-in-water emulsions were contained within home-made PDMS microfluidic chips featuring a single channel with square cross-section of 160 × 160 μm2. These chips supported two multi-mode optical fibers with 50 μm core and 125 μm cladding diameters. The faces of the fibers were parallel to the opposite side walls of the channel and their optical axes were precisely aligned to form the dual-beam trap (see inset in Fig. 1(a)). The microfluidic channel had one inlet and one outlet with a syringe connected to the inlet via tubing to inject previously prepared microdroplet emulsions and washing solutions. Detailed manufacturing process for the chips is given in the next subsection. In order to observe the droplets, excite and collect the fluorescence emission from the droplets, and implement single beam optical trap, PDMS chips were mounted on an xyz stage of an inverted optical microscope equipped with a water immersion objective (NA=1.2, 60×; Nikon).
Laser light source for the dual-beam trapping and optical stretching was a 1070 nm Ytterbium fiber laser (YLM-10-LP-SC; IPG Photonics) with 10 W maximum output power. Its output beam was directed to a λ/2 waveplate (WP1) mounted on a stepper-motor-driven rotary stage and a polarizing beam splitter (PBS1) which together served to adjust the total stretching power coupled into the fibers. Following PBS1, another polarizing beam splitter PBS2 was used to divide the beam into two parts, with λ/2 waveplate WP2 mounted on a manual rotary stage adjusting the division ratio of the total beam power between the fibers. Subsequently, both beams were coupled to respective single-mode patch cables (P5-980A-PCAPC-1; Thorlabs) using compact fiber couplers (PAF-X-11-C; Thorlabs). The output ports of the patch cables were then connected to the free ends of two multi-mode fibers embedded in the microfluidic chip. In order to measure the total actual power coupled to the fibers, a glass plate with an angle of 45 degrees relative to the laser beam was inserted in the beam path after PBS1; this plate reflected a small fraction of the total power (∼5 %) towards a power meter. Thus measured power fraction was calibrated by measuring the total output power from the bare fiber ends before inserting them into the chip.
Even though the power balance between the two stretching beams could be adjusted by WP2 for a constant total injected laser power, the power ratio displayed small fluctuations when the total injected power was modified by WP1. Consequently, a droplet confined in the dual-beam optical trap moved in all three directions with respect to the observation focal point when the total stretching power was altered. This phenomenon was further augmented by slight imperfections in the coaxial alignment of the two fibers. In addition, immersion oil used in the experiments was slightly denser than water and, thus, the oil droplets tended to sediment on the bottom of the microchannel. Since the transversal gradient optical force of the dual beam trap that counteracted gravity in our experimental arrangement (gravity pointing in negative z-direction) depends linearly on the total beam power, changes of the dual beam trapping power led to additional vertical displacement of the trapped droplet from the observation focal plane. Such unwanted droplet displacements changed the efficiency of the droplet pumping by the pulsed laser beam (see below) and, subsequently, the intensity of lasing emission from the trapped droplet. To overcome these issues and fix the position of the droplet in the axial and transversal directions during the droplet stretching, an additional infrared laser beam from a cw solid state laser (wavelength 1064 nm, maximum power 300 mW; CrystaLaser) was introduced to create an independent single beam trap (optical tweezers) perpendicular to the dual beam trap (see inset in Fig. 1(a)). This laser beam was sent through a 3× beam expander and reflected from two dichroic mirrors D1, D2 before being focused into the sample chamber by the water immersion microscope objective. The overall transmission of all optical elements in the single beam trapping path was measured to be about 20%, resulting in the average trapping power of ∼ 60 mW at the specimen. Using the single beam trap, the studied droplet could be brought exactly to the position of the dual beam trap by adjusting the xyz stage of the inverted microscope. Furthermore, in order to change the position of the single beam trap in the objective focal plane xy, the first lens of the beam expander was mounted on an xy stage that could be steered in x and y directions (perpendicular to the beam axis).
A 532 nm pulsed laser beam generated by frequency-doubling the output of a home-built passively Q-switched Nd – YVO4 laser (20 ns pulse width, 33 kHz repetition rate, maximal average output power 20 mW) was used as the pump light source to excite fluorescently-stained trapped droplets. The pump laser beam was focused by the same microscope objective as the optical tweezers laser beam, after being transmitted through the dichroic mirror D1 and reflected from the dichroic mirror D2. The diameter of the pump beam focal spot was estimated from the beam spot images acquired after the beam reflection from the cover glass to be 1.24 μm. A 2× beam expander with the first lens mounted on an xy stage was added to the optical path of the pump beam to steer the beam focus position in 3-D (xy stage was used to steer the focus of the beam in xy direction and the distance between the two lenses of the beam expander was used to adjust its focal point position relative to the objective lens in the axial direction). This way, focusing points of the optical tweezers beam and the pump beam could be independently adjusted relative to each other and also relative to the imaging focal plane of the objective. Before the experiments, the droplet position was adjusted by moving the optical tweezers in the focal plane so that the focal spot of the pump laser lied near the droplet’s rim for the maximal efficiency of the WGM excitation.
Fluorescence from a trapped microdroplet was collected using the same microscope objective as used for focusing the optical tweezers and pump beams, and dispersed by a monochromator (focal length 500 mm; Acton Research) on the chip of a cooled CCD detector (Pixis 100; Princeton Instruments) controlled by WinSpec application software (Princeton Instruments). In Fig. 1(c), a typical fluorescence spectrum of a trapped microdroplet is shown. The high-intensity WGMs located in the spectral interval 620 – 630 nm indicate droplet lasing. In order to prevent fast photobleaching of dye molecules inside the droplets, a shutter was added to the pump laser beam path; this shutter was only opened during the spectrum acquisition time and provided trigger signal for the spectroscopic detector to synchronize the data recording. To observe and take an image of the sample, the light in the imaging path of the microscope could be diverted to an ordinary CCD camera via a flippable mirror (FM).
In order to automate the change of the stretching power and synchronize it with the collection of the droplet fluorescence spectra at each level of the stretching power, a LabVIEW code was developed to rotate the stepper motor, open the shutter for 8 ms, and record the stretching laser power while, at the same time, externally triggered WinSpec software was collecting the droplet emission spectra with an exposure time of 8 ms. This way, the droplet deformation cycle could be carried out as fast as possible to minimize the effects of the droplet heating during the spectrum recording interval.
3. Microfluidic chip fabrication
Microfluidic chips for fiber-based dual-beam optical trapping and stretching were manufactured from PDMS elastomer by casting liquid PDMS precursor into a home-made reusable mold (master). The chip master was prepared in the following manner: first, a piece of optical fiber stripped of the plastic jacket with a diameter of 125 μm was glued onto a cleaned glass substrate with a UV activated adhesive (Loctite 3922 or 3494). This fiber fragment served to form insertion ports for the dual-beam trap fibers in the chip. In the middle of the fiber length, two square capillary tubes (external width of 160 μm) were aligned coaxially with each other and perpendicularly to the fiber, one from each side of the fiber, and glued to the substrate. These capillary tubes defined the liquid channel in the chip. Central openings of the capillary tubes adjacent to the fiber were connected by UV activated adhesive. Resultant cross-like structure with a stalk formed by the capillary tubes and arms in the middle formed by the stripped optical fiber was placed into a plexiglass casting frame that was filled with liquid PDMS precursor and cured at the temperature of ∼ 85°C for 2 hours. After solidification, the resulting PDMS block was peeled off the mold and the PDMS surface was activated together with a 150 μm thick microscope cover slip in an oxygen-argon microwave-excited plasma (oxygen-argon 1:1, total pressure 500 mTorr, microwave power 50 W). Subsequently, the activated PDMS chip was pressed against the activated cover slip, thus sealing the chip permanently. The resulting chip featured one straight channel with square cross-section of 160 × 160 μm2 and two perfectly aligned cylindrical slots for optical fibers in the middle. The liquid inlet and outlet of the channel were realized during PDMS casting with a wire inserted into the mold or after PDMS solidification by a hole-punch.
In order to achieve optimal stretching forces for droplets with diameters around 50 μm, beam diameters at the droplet stable trapping position need to be comparable with the droplet size. This requirement was addressed by using multimode optical fibers (AFS50/125Y, NA=0.22; Thorlabs) with core diameters of 50 μm. These fibers were striped of their plastic jacket, cleaved, and inserted into the slots in the chip with the help of two small 3-D translational stages, one for each fiber. When the optimal positions of the fibers were reached, connections between the chip and the fibers were sealed using epoxy glue applied to the fiber slots from outside. In order to enable cleaning of the fibers, separation distance between the front faces of the fibers was chosen a little smaller than the width of microfluidic channel, equal to 134 μm.
4. Immersion oil/water emulsion system
The emulsion system selected for our study consisted of microdroplets of dye-doped microscope immersion oil (n = 1.515, ρ = 1.02 g/cm3, Merck) dispersed in deionized water (n = 1.334, ρ = 1.00 g/cm3). This emulsion system of virtually immiscible liquids provides a high refractive index contrast together with a low density contrast between the droplet and the surrounding medium. This implies a good confinement of WGMs inside the trapped micro-droplet, strong optical trapping force, and a small buoyancy that further facilitates the optical trapping of the droplets. The dye DiI(3) (1,1’-Dioctadecyl-3,3,3’,3’-tetramethylindo carbocyanine perchlorate, Sigma Aldrich) that has a hydrophilic chromophore and hydrophobic side chains was selected as the laser gain medium. This dye has its excitation peak at ∼ 560 nm and emission peak at ∼ 575 nm, and is well suited for the efficient excitation of the WGMs due to its surfactant nature that ensures its localization at the oil-water interface where WGMs also reside [16,25]. In our experiments DiI(3) was dissolved in chloroform at 10 mM concentration and the dye solution was subsequently added to the immersion oil in 2 to 100 ratio. This solution was put in the oven at ∼ 80°C and evacuated for 10 min to remove the chloroform content. Hence, the final solution contained only immersion oil with the dye concentration of 200 μM.
5. Characterization and adjustment of interfacial tension of emulsion droplets
Interfacial tension (γ) between two immiscible liquids can be determined using Wilhelmy plate method . In this method, contacting planar layers of studied liquids are formed in a container and a solid plate from a material that is wetted by the denser liquid is brought into contact with the liquid-liquid interface, perpendicularly to the interface. Subsequently, vertical force acting on the plate is measured that is directly proportional to γ. Since the denser phase in our experimental system was immersion oil, we used Wilhelmy plate made of Teflon that is wetted well by oil. Using Wilhelmy plate method, we measured the value of γ ≈ 13.1mN/m for the dye-doped immersion oil/water interface.
Tuning of the lasing emission spectra of dye-doped microdroplets by optical stretching requires lowering of surface forces to the level comparable with optical forces. In order to achieve sufficiently low interfacial tensions, we added AOT surfactant (Docusate sodium salt; Sigma-Aldrich) and sodium chloride to the aqueous phase. This surfactant/salt system has been shown to decrease γ to extremely low values (down to 0.001 mN/m) depending on the concentration of both the surfactant and NaCl [24, 27]. However, such low values of γ are problematic in experiments where the droplets serve as optical resonant cavities since their shape can be easily disturbed by thermal fluctuations which cause degradation of the droplet cavity Q-factor. Here, we used 5 mM AOT surfactant and 40 mM NaCl to decrease the interfacial tension to γ = (1.5 ± 0.1)mN/m. We measured the interfacial tension in the presence of surfactant and salt in a range of temperatures between 25°C and 50°C and found a negligible variation of γ with temperature. Thus, heating of droplets due to residual absorption of the stretching laser light does not cause any appreciable changes of surface forces.
Emulsion microdroplets were obtained by mixing the dye-doped immersion oil with the surfactant-salt-water solution in 1 to 50 ratio, and shaking the obtained mixture by hand. Droplets produced in this way had a large size distribution with diameters ranging from < 1 μm up to ∼ 100 μm and the emulsions were stable over ∼ 10 hours. Prepared emulsions were subsequently loaded into the PDMS chip using liquid flow controlled by gravitational hydraulic pressure built by the height of water column in a 20 cm long piece of tubing connected to the inlet of the microchip. After loading, microdroplets were observed to float freely within the channel, thanks to their small buoyancy, without sticking to the walls of the sample chamber.
6. Simulations of WGM tuning by droplet stretching
In order to calculate the WGM tuning as a function of the stretching power and find out the effects of interfacial tension, droplet size, beam profiles, fiber distance, and refractive index of the used liquids on the tuning range, a simulation of the droplet deformation by optical stretcher was carried out based on a 3-D ray-tracing method developed in . In this method, an incident beam is divided into geometrical rays of variable intensities perpendicular to the wavefront. Each of these rays is then traced through the droplet interface using Snell’s law and for each refraction or reflection, the transmitted and reflected powers are calculated for TE and TM polarization by related Fresnel coefficients. Subsequently, the forces acting on the interface as a result of ray reflection or refraction can be calculated from the momentum conservation principle. To calculate the total stretching force (Fstretch) acting on the droplet, optical forces of all individual rays forming the incident beam are summed. This ray tracing is continued until the remaining power in the reflected/refracted rays becomes sufficiently small.
If the fibers are assumed to be coaxial and their output beam profiles and powers to be identical, droplet deformation will be symmetric relative to its equatorial plane xz shown in Fig. 2(a). Since the optical stretching forces in moderately focused beams point predominantly along the beam propagation axis y, the droplet will elongate in the y-direction and assume the shape that can be approximated by a prolate spheroid for small deformations. This shape change increases the droplet surface area and, thus, its surface energy. Consequently, a restoring surface force Fsurf arises that grows with increasing droplet deformation until it exactly balances the optical stretching force Fstretch. It should be born in mind that Fstretch depends on the droplet shape and, thus, it changes with the droplet deformation.
Let us describe the deformed droplet shape by a prolate spheroid , where a is the polar semi-major axis along the dual beam trap axis y, and b is the equatorial semi-minor axis. The droplet volume is conserved after stretching so and, thus, , where r0 is the radius of the unperturbed drop and a = r0 + Δa. In equilibrium, Fstretch equals to Fsurf that depends on the droplet elongation Δa along the y-axis as :29]: Figure 2(b) shows the results of a simulation for a 50 μm diameter oil droplet placed in aqueous surfactant solution with γ = 1.5 mN/m and located at the center of the dual beam trap. The trap is generated between two multimode fibers with cores of diameter 50 μm and numerical aperture NA = 0.22, and fiber facets are separated by a distance 134 μm. In this figure we plot the wavelength tuning of WGMs with mean wavelength λ = 625 nm circulating in the equatorial plane (off-axis excitation) and the polar plane (on-axis excitation) as a function of the stretching laser power (Pstretch) which corresponds to the total power delivered from the two multimode fibers. For both fibers we assume beam waist diameter equals 50 μm and the beam waists are located at the fiber end facets. At the trap center located 67 μm away from the fiber end facets, beam diameter of 55 μm was found assuming an embedded Gaussian beam with 50 μm beam waist diameter propagating with a correction M-factor of 4.0 given by the asymptotic divergence angle of 9.4° that corresponds to the numerical aperture of the fiber output. As we can see, for WGMs propagating in the equatorial plane, tuning slope is negative and equal to −0.96 nm/W, while WGMs propagating in the polar plane are tuned with a positive slope of 0.48 nm/W. Thus, the direction of the lasing wavelength tuning depends on the orientation of the modal path of particular WGM with respect to the droplet planes of symmetry. Due to approximate description of the stretching laser beam profiles and assumption of droplet deforming precisely towards a prolate spheroid, calculated values of the WGM tuning slopes and the tuning slope ratio may differ from the true ones. However, the direction of WGM tuning by optical stretching is predicted correctly.
7. Observation of tunable dye lasing in optically stretched microdroplets
Experiments with tunable dye lasing in optically stretched microdroplets were carried out using multimode fibers with parameters identical to those used in the simulations. In order to excite different WGMs in the stretched droplets, the focus of the pump beam was moved to different locations on the droplet.
First, the pump beam was focused at the rim of the droplet in the equatorial plane perpendicular to the dual beam trap axis (off-axis excitation in Fig. 1(b)). Figure 3(a) illustrates spectral tuning of three consecutive lasing WGMs with off-axis excitation for a droplet with the diameter of 47 μm and γ ≈ 1.5mN/m. Pstretch was changed from 0 W to 0.88 W in two cycles and the total of 38 spectra were recorded with the minimal values of Pstretch corresponding to frame numbers 1, 19, and 38 and the maximal values of Pstretch corresponding to frame numbers 10 and 29. Time interval between acquiring consecutive spectra was 50 ms and the droplet was exposed with the pump beam fluence of ∼ 490mJ/cm2 at the specimen for 8 ms. From Fig. 3(a), it is clear that the equatorial plane WGMs experience power-dependent blue tuning consistent with the simulation results presented in Fig. 2. As indicated by the white vertical line located near the middle WGM in Fig. 3(a), the WGMs do not return back to their initial spectral positions after completing a single tuning cycle and a small blue drift of ∼ 0.04nm is observed in all lasing WGMs. We interpreted this drift as a result of gradual dissolution of the droplet medium in water triggered by heating of the droplet by the stretching laser beams. This explanation was supported by the fact that the drift amount was directly proportional to the total duration of the tuning cycle and the maximal Pstretch (data not shown). During data recording, the effect of the droplet dissolution was minimized by performing the experiments as fast as possible. In order to take the effect of residual droplet dissolution into consideration, we assumed a constant dissolution rate during each tuning cycle resulting in a linear drift of the WGM positions with time. Hence, WGM positions were corrected by subtracting a line connecting the first and last raw WGM position in the tuning cycle corresponding to Pstretch = 0W. This correction was applied to all of the tuning slope calculations presented in this paper. Figure 3(b) shows the tuning of the lasing WGM wavelength for the central peak in Fig. 3(a) corrected for droplet dissolution as a function of Pstretch for two subsequent tuning cycles. The slope of a linear fit to the experimental data is −0.52 nm/W which is ∼ 54% of the predicted tuning slope of −0.97 nm/W. We address this discrepancy between the experiment and the simulation in the next section.
In the second set of experiments, pump beam was focused at the droplet rim in the polar plane containing the dual beam trap axis (on-axis excitation in Fig. 1(b)) and all other experimental parameters were kept fixed. The experiments in both configurations of the pump beam focus were carried out with different droplets since it was not possible to hold the droplet in the trap during the change of the pump beam position. Figure 3(c) shows spectral tuning of three consecutive lasing WGMs with on-axis excitation for a droplet with the diameter of 49 μm. Pstretch was changed from 0 W to 0.75 W in two cycles and the total of 38 spectra were recorded with the minimal Pstretch corresponding to frame numbers 1, 19, and 38 and the maximal Pstretch corresponding to frame numbers 10 and 29. As in the previous experiment, time interval between acquiring consecutive spectra was 50 ms and the droplet was exposed with the pump beam fluence of ∼ 490 mJ/cm2 at the specimen for 8 ms. As illustrated in Fig. 3(d), for on-axis excitation, lasing WGMs experience red tuning that increases linearly with Pstretch. Here, a linear fit gives positive tuning slope of 0.49 nm/W that is ∼ 2% larger than the predicted value of 0.48 nm/W.
Tunable lasing experiments with different configurations of the pump beam reveal a systematic shift of experimental tuning slopes with respect to the predictions towards the red end of the spectrum. As argued in the next section, this shift can be explained by simultaneous changes of the droplet size and refractive index of both the droplet and the host medium resulting from heating by the stretching laser beams.
8. Influence of thermal effects on WGM tuning
Both the simulations shown in Fig. 2 and the asymptotic formula giving the spectral positions of WGMs with different azimuthal mode numbers as a function of the droplet ellipticity [11,12] predict the ratio of the spectral tuning slopes of the off-axis WGMs (azimuthal mode number = angular momentum index) and the on-axis WGMs (azimuthal mode number = 0) to be −2.0 for droplet sizes studied in this paper. However, in our experiments this ratio was observed to be −1.1. In order to elucidate this difference, we carried out control experiments with dye-doped immersion oil droplets suspended in pure water without surfactant. As mentioned earlier, interfacial tension for oil/water system γ ≈ 13.1mN/m is almost nine times larger than interfacial tension in the presence of surfactant γ ≈ 1.5mN/m. Consequently, one expects nine times smaller tuning slope for the surfactant-free emulsion system in comparison to the data presented in Figs. 3(b) and 3(d).
In Fig. 4(a), spectral tuning of lasing WGMs with off-axis excitation is shown for a 44 μm diameter immersion oil droplet suspended in pure water. The excitation power, exposure time, and time between subsequent frames were identical to the experiments performed with droplets containing surfactant (Fig. 3). In the surfactant-free experiments, similarly to the previous cases with surfactant, a monotonous blue drift was observed resulting from the droplet dissolution. In addition, a small cyclic tuning of the WGM positions was also detected with changing Pstretch. However, unlike in the previous case with surfactant, spectral tuning for off-axis excitation was towards the red end of the spectrum. A line fit to the dissolution-corrected tuning data shown in Fig. 4(b) gives a positive tuning slope of 0.085nm/W, in contrast to our prediction of a negative slope of −0.11nm/W from the simulations. In Fig. 4(c), spectral tuning of lasing WGMs with on-axis excitation is shown for a 55 μm diameter immersion oil droplet suspended in pure water. Again, a line fit to the dissolution-corrected tuning values in Fig. 4(d) shows a positive tuning slope of 0.138nm/W that has a proper sign but is 2.6 times larger than the prediction of 0.054nm/W.
Observation of red-tuning independent of the pump beam position suggests that in the surfactant-free case, WGM tuning is not dominated by droplet stretching which is associated with opposite tuning directions for off-axis and on-axis WGMs. Instead, as we discuss below, heating of the emulsion system by the optical stretching laser beams becomes significant, leading to temperature-induced red shifts observed in WGM spectral positions. Relatively large water absorption at the stretching laser wavelength of 1070 nm leads to an increase in the temperature of water and the droplet by more than 10°C per 1 W of Pstretch . For viscous liquids such as immersion oil and water, the dominant mechanism of heat transfer on the scale of hundreds of microns is conduction. For the heat conductivity of water and immersion oil (∼ 0.5Wm−1K−1), steady state is reached within 10 – 20 ms after temperature jumps of 10°C. Since the duration of individual steps of the droplet tuning cycle is > 50 ms, we can assume that the temperature distribution within the droplet reaches equilibrium before each spectral acquisition. Immediate results of such temperature changes are the droplet expansion and the change of the refractive index of the droplet and water. As mentioned in Section 5, interfacial tension does not change appreciably within the temperature range of ∼ 25° C above the room temperature. Therefore, temperature-induced changes of surface forces that could cause additional spectral tuning can be ruled out in the first approximation.
In order to predict the contributions of thermal expansion and refractive index changes to the WGM tuning slopes, we adopted an asymptotic formula which gives the spectral positions of WGMs in terms of oil/water relative refractive index n1/n2 and droplet radius r0. This asymptotic formula holds for droplets with large size parameters x = 2πr0n2/λ and its use is justified in our case where typical r0 ≈ 25 μm and λ ≈ 625nm give x ≈ 335. First derivative of this asymptotic formula with respect to temperature gives the following expression for thermal tuning of WGM positions of spherical droplets:
The above simple model of thermally induced tuning shows that in the off-axis excitation geometry, there are two competing mechanisms that tune the WGM spectral positions in opposite directions. While droplet stretching leads to a blue-tuning of WGMs, heating effects tune WGMs to the red end of the spectrum. The overall tuning slopes we observe in the experiments then result from the superposition of these two effects. For emulsions with surfactant, optical stretching dominates and the total wavelength tuning slope is negative (see Figs. 3(a) and 3(b)) whereas in the surfactant-free emulsions, the total wavelength tuning slope is positive due to the dominance of heating effects (see Figs. 4(a) and 4(b)). In the on-axis excitation geometry, both thermal and stretching tuning mechanisms lead to positive tuning slopes which explains why the measured tuning slope in Figs. 4(c) and 4(d) is more than the amount expected solely from the droplet deformation due to optical stretching.
In all experiments in which the excitation beam position was either on-axis or off-axis, we studied only WGMs propagating in a single plane of symmetry of an optically stretched droplet (polar or equatorial, respectively). Thus, it was not possible to monitor simultaneously the spectral positions of both blue-shifting and red-shifting WGMs for a single stretched droplet. However, such information is crucial for determining quantitatively the relative contribution of droplet stretching and thermal effects to the observed WGM tuning. In order to excite simultaneously WGMs propagating along different paths in a single droplet, the excitation beam position was moved to the droplet center (see Fig. 1(b)). Figure 5 summarizes experimental results obtained with on-center excitation of a 42 μm diameter droplet for γ = 1.5mN/m. With this excitation geometry, WGMs in the droplet fluorescence spectra display a spectral split δλ that grows with increasing Pstretch (see Fig. 5(a)). Within the single frequency-split WGM group, the shortest-wavelength WGM with λB propagates in the equatorial plane whereas the longest-wavelength WGM with λR propagates in the polar plane. As derived by Chen et al. , wavelength λ0 corresponding to the degenerate WGMs propagating in an undeformed spherical droplet is given by λ0 = λB + 2δλ/3. This formula allows calculating the value of λ0 for each stretching power from the experimental data. In case of pure stretching deformation of the droplet, λ0 must be independent of Pstretch. Non-negligible thermal effects then induce changes of λ0 with changing Pstretch. In Fig. 5(b), values of λB, λR, and λ0 are shown for WGM1, WGM2 labeled in Fig. 5(a). Details of changes of λ0 of the two WGMs with changing Pstretch presented in Fig. 5(c) then clearly illustrate that the thermal tuning effects cannot be neglected. A linear fit of λ0 as a function of Pstretch given in Fig. 5(d) yields the average thermal tuning slope of 0.085 nm/W. This value agrees well with our estimation of 0.082 nm/W obtained from the above presented model of thermal tuning effects for the droplet studied in Fig. 5. Measured thermal tuning slope can be used to estimate pure stretching-induced tuning observed in the case of emulsions with surfactants (Fig. 3). For off-axis and on-axis excitations, the pure stretching tuning slopes given by the difference between total measured tuning slopes and the average thermal tuning slope become −0.61 nm/W and 0.41 nm/W, respectively. The experimental ratio of the stretching tuning slopes is then −1.5 which is close to the value of −2.0 predicted by the simulations and the analytical formula for deformed droplet WGMs. The remaining difference between the experiments and the prediction may be attributed to several factors. One of them is the fact that individual tuning slopes measured in various excitation geometries were obtained with droplets of different diameters. As the values of both the stretching and thermal tuning slopes depend on the droplet size, experimental ratio of the corrected tuning slopes calculated from measurements carried out with different droplets is generally biased. Also, surfactant concentration in individual droplets is likely not identical which results in variations of the droplet surface tension and, thus, different droplet deformations under otherwise identical experimental conditions. In addition, non-uniform temperature distribution within the droplets heated during stretching  induces local variations of refractive index around the droplet circumference. These variations differ for the polar and equatorial WGMs and cause deviations of thermally-induced WGM spectral tuning with respect to the case of a uniformly heated droplet. Finally, despite small magnitude of optical stretching deformations, the exact shape of stretched droplets may differ from a prolate spheroid. For such general droplet shapes, ratio of tuning slopes different from −2.0 is expected.
In this article, we have presented tunable optofluidic dye microlasers based on active optical resonant cavities formed by optically trapped microdroplets of fluorescently stained immersion oil emulsified in water. All-optical tuning of the laser emission wavelength was achieved by a controlled deformation of the droplet shape using light-induced forces generated by dual-beam optical trap. By lowering interfacial tension between the droplets and the aqueous host medium with a surfactant and sodium chloride, we have demonstrated almost reversible spectral tuning of the lasing droplet WGMs with tuning slopes up to ∼0.5 nm per Watt of the total stretching laser power. These spectral tuning slopes can be increased by further reduction of the droplet interfacial tension. Droplet liquid was observed to slightly dissolve in the host liquid during each tuning cycle, resulting in an irreversible blue-shift in the WGM spectral positions. This can be prevented by further optimization of the emulsion system in order to obtain fully reversible spectral tuning of the WGMs. Tuning direction of our liquid-cavity microlasers could be easily reversed by changing the position of the optical pumping beam on the droplet. Comparison of experimental data with simulations of the WGM tuning revealed an additional mechanism of reversible spectral tuning triggered by thermally induced changes of the droplet size and refractive index of both the droplet and the host medium. Such thermally induced tuning is a consequence of the absorption of the near-infrared optical stretcher light in water; it can be avoided by changing the wavelength of the used laser beam or replacing aqueous host medium with deuterium oxide that has minimal absorption in the near-infrared spectral region. Our droplet-based tunable and disposable optofluidic lasers can be used with a wide choice of fluorescent dyes and they can be easily integrated into microfluidic lab-on-a-chip platforms. Thus, they represent an attractive alternative for tunable sources of light for on-chip applications in chemical and biological analysis and sensing.
The authors acknowledge support from TÜBİTAK (Grant No. 111T059), Czech Science Foundation ( P205/11/P294), Technology Agency of the Czech Republic ( TA03010642), Ministry of Education, Youth and Sports of the Czech Republic together with the European Commission (ALISI CZ.1.05/2.1.00/01.0017). The authors thank Martin Šiler for help with thermal analysis and A. L. Demirel for help with the surface tension measurements.
References and links
1. V. V. Datsyuk, “Optics of microdroplets,” J. Mol. Liq. 93, 159–175 (2001). [CrossRef]
2. G. C. Righini, Y. Dumeige, P. Feron, M. Ferrari, G. N. Conti, D. Ristic, and S. Soria, “Whispering gallery mode microresonators: Fundamentals and applications,” Riv. Nuovo Cimento 34, 435–488 (2011).
6. J. B. Snow, S.-X. Qian, and R. K. Chang, “Stimulated Raman scattering from individual water and ethanol droplets at morphology-dependent resonances,” Opt. Lett. 10, 37–39 (1985). [CrossRef] [PubMed]
7. A. Sennaroglu, A. Kiraz, M. A. Dündar, A. Kurt, and A. L. Demirel, “Raman lasing near 630 nm from stationary glycerol-water microdroplets on a superhydrophobic surface,” Opt. Lett. 32, 2197–2199 (2007). [CrossRef] [PubMed]
8. A. Kiraz, A. Kurt, M. A. Dündar, and A. L. Demirel, “Simple largely tunable optical microcavity,” Appl. Phys. Lett. 89, 081118 (2006). [CrossRef]
9. A. Kiraz, Y. Karadağ, and A. F. Coskun, “Spectral tuning of liquid microdroplets standing on a superhydrophobic surface using electrowetting,” Appl. Phys. Lett. 92, 191104 (2008). [CrossRef]
11. G. Chen, M. M. Mazumder, Y. R. Chemla, A. Serpengüzel, R. K. Chang, and S. C. Hill, “Wavelength variation of laser emission along the entire rim of slightly deformed microdroplets,” Opt. Lett. 18, 1993–1995 (1993). [CrossRef] [PubMed]
12. S. C. Yorulmaz, M. Mestre, M. Muradoglu, B. E. Alaca, and A. Kiraz, “Controlled observation of nondegenerate cavity modes in a microdroplet on a superhydrophobic surface,” Opt. Commun. 282, 3024–3027 (2009). [CrossRef]
14. M. Humar, M. Ravnik, S. Pajk, and I. Muševič, “Electrically tunable liquid crystal optical microresonators,” Nat. Photonics 3, 595–600 (2009). [CrossRef]
16. M. Aas, A. Jonáš, and A. Kiraz, “Lasing in optically manipulated, dye-doped emulsion microdroplets,” Opt. Commun. 290, 183–187 (2013). [CrossRef]
17. S. K. Y. Tang, Z. Li, A. R. Abate, J. J. Agresti, D. A. Weitz, D. Psaltis, and G. M. Whitesides, “A multi-color fast-switching microfluidic droplet dye laser,” Lab Chip 9, 2767–2771 (2009). [CrossRef] [PubMed]
19. J. Guck, R. Ananthakrishnan, H. Mahmood, T. J. Moon, C. C. Cunningham, and J. Käs, “The optical stretcher: A novel laser tool to micromanipulate cells,” Biophys. J. 81, 767–784 (2001). [CrossRef] [PubMed]
20. A. Ashkin, “Acceleration and trapping of particles by radiation pressure,” Phys. Rev. Lett. 24, 156–159 (1970). [CrossRef]
22. T. Čižmár, O. Brzobohatý, K. Dholakia, and P. Zemánek, “The holographic optical micro-manipulation system based on counter-propagating beams,” Laser Phys. Lett. 8, 50–56 (2011). [CrossRef]
23. P. C. F. Møller and L. B. Oddershede, “Quantification of droplet deformation by electromagnetic trapping,” EPL 88, 48005 (2009). [CrossRef]
24. A. D. Ward, M. G. Berry, C. D. Mellor, and C. D. Bain, “Optical sculpture: controlled deformation of emulsion droplets with ultralow interfacial tensions using optical tweezers,” Chem. Commun. 2006, 4515–4517 (2006). [CrossRef]
25. S. Holler, N. L. Goddard, and S. Arnold, “Spontaneous emission spectra from microdroplets,” J. Chem. Phys. 108, 6545–6547 (1998). [CrossRef]
26. J. Drelich, C. Fang, and C. L. White, The Encyclopedia of Surface and Colloid Science: Measurement of Interfacial Tension in Fluid/Fluid Systems(Marcel- Dekker, 2002).
27. R. Aveyard, B. P. Binks, S. Clark, and J. Mead, “Interfacial tension minima in oil-water-surfactant systems. Behaviour of alkaneaqueous NaCl systems containing aerosol OT,” J. Chem. Soc., Faraday Trans. 1 , 82, 125–142 (1986). [CrossRef]
28. H. Sosa-Martínez and J. C. Gutiérrez-Vega, “Optical forces on a Mie spheroidal particle arbitrarily oriented in a counterpropagating trap,” J. Opt. Soc. Am. B 26, 2109–2116 (2009). [CrossRef]
29. S. Ramanujan, Collected papers of Srinivasa Ramanujan (Cambridge University Press, 1927).
31. C. C. Lam, P. T. Leung, and K. Young, “Explicit asymptotic formulas for the positions, widths, and strengths of resonances in Mie scattering,” J. Opt. Soc. Am. B 9, 1585–1592 (1992). [CrossRef]