Liquid Gradient Refractive Index Microlens for Dynamically Adjusting the Beam Focusing

An in-plane liquid gradient index (L-GRIN) microlens is designed for dynamically adjusting the beam focusing. The ethylene glycol solution (core liquid) withde-ionized (DI) water (cladding liquid) is co-injected into the lens chamber to form a gradient refractive index profile. The influences of the diffusion coefficient, mass fraction of ethylene glycol and flow rate of liquids on the refractive index profile of L-GRIN microlens are analyzed, and the finite element method and ray tracing method are used to simulate the convection-diffusion process and beam focusing process, which is helpful for the prediction of focusing effects and manipulation of the device. It is found that not only the focal length but the focal spot of the output beam can be adjusted by the diffusion coefficient, mass fraction and flow rate of liquids. The focal length of the microlens varies from 942 to 11 μm when the mass fraction of the ethylene glycol solution varies from 0.05 to 0.4, and the focal length changes from 127.1 to 8 μm by varying the flow rate of the core liquid from 0.5 ˆ 103 to 5 ˆ 103 pL/s when there is no slip between the core and cladding inlet. The multiple adjustable microlens with a simple planar microfluidic structure can be used in integrated optics and lab-on-chip systems.


Introduction
Tunable microlenses are widely used in microfluidic or lab-on-chip systems [1] due to their fine tuning of the light.A tunable microlens adaptively reshapes the input light including adjusting its focal position, intensity, beam profile, and even propagation direction [2].They can also be used for ordinary coupling, collimating and focusing, and in wide range of lab-on-a-chip applications such as flow cytometry [3][4][5], on-chip optical tweezers and single molecule detection [6,7].
Several kinds of tunable microlenses are developed for providing adaptive focusing with adjustable curved refractive surfaces [2], such as a micromeniscus surface actuated by electro-wetting [8], liquid microlenses based on the dielectrophoretic effect [9][10][11], microlenses based on the thermal effect [12], or hydraulically actuated polydimethylsiloxane (PDMS) membranes [13,14].Among them, electro-wetting and dielectrophoresis lenses are electrically driven, while other microlenses whose focal lengths are variable are either based on pressure-induced systems using an elastic membrane, or on thermal effect.A mechanical lens whose focal length is controlled by pumping liquid in and out of the lens chamber caused great interest.The operation mechanism of this lens is simple because it only requires a fluid pumping system.There are two working mechanisms which are often used to implement these tunable liquid microlenses.One is to achieve the refraction of light at curved, optically smooth, liquid-liquid interfaces between two flowing fluids by controlling laminar flows.Recently, the Whitesides group has developed a dynamically reconfigurable liquid-liquid lens [15] with the convex shape of liquid-liquid interfaces in a microfluidic expansion chamber.The curvatures in the liquid-liquid interfaces were controlled by adjusting hydrodynamic flow conditions, and then adjustable focal length was produced.However, in the liquid-liquid lens, a higher speed of the laminar flows is necessary, which means large amounts of fluid have to be supplied discontinuously for keeping the microlens working.Utilization of so-called liquid gradient refractive index (L-GRIN) microlenses [16] is the other approach and also the preferable way to achieve a tunable microlens.Similarly, with the traditional solid self-focusing lens, the refractive index of the L-GRIN microlens increases continuously from the cladding region to the core optical axis.However, a gradient refractive index of the L-GRIN microlens is achieved by dynamically adjusting hydrodynamic flow conditions rather than by changing the geometrical structure of the lens for the traditional one.In contrast to the liquid-liquid interface-based microlens, the L-GRIN lens operates through the diffusion in multiple flows rather than by relying on a clearly defined, curved liquid-liquid interface.Therefore, the flow consumption rate is much lower in the L-GRIN lens [17].These lenses focus a light beam through a liquid medium with a two-dimensional (2D) refractive index gradient.The two-dimensional (2D) refractive index profile and, subsequently, the focal length of the L-GRIN lens can be tuned by changing the ratio of the flow rates of the two inlets [17].However, most L-GRIN microlenses are only capable of out-of-plane focusing.In-plane-focusing tunable microlens are, however, proven to be preferable in integrative devices, which can be fabricated and seamlessly integrated with other on-chip fluidic and optical components such as a lab-on-a-chip lasers [18] and optical waveguides [19].In particular, it has inspired the creation of a variety of innovative devices controlled by flow rates and liquid compositions, and variable light focusing was measured and shown, but the quantitative relation between liquid parameters and focal length has not been demonstrated clearly.It will lead to the unpredictability of the focusing effect in the device.
In this paper, an in-plane tunable L-GRIN microlens was designed for dynamically adjusting the beam focusing, which can be more readily integrated for lab-on-a-chip applications.The convection-diffusion process of liquids in the microfluidic chamber is firstly simulated with the finite element method (FEM) when a high-refractive-index solution is injected side-by-side into a low-refractive-index solution.During the convection and diffusing process, the diffusion profile of the liquids, and hence the refractive index profile within the L-GRIN microlens, varies with the hydrodynamic flow conditions.Thus, the refractive index profile of the L-GRIN microlens is calculated and discussed numerically under different flow conditions, and the beam transmission and focusing process in the L-GRIN microlens is simulated using the ray tracing method.The effects of the diffusion coefficient, mass fraction and flow rate of liquids on the input beam focusing effect, including the focal length and the size of the focal spot, were demonstrated.Therefore, an in-plane tunable beam focusing is achieved using a simple planar microfluidic structure; in addition, with the simulation of the convection and diffusing process, the ray tracing method gives us approaches to predict the focusing effect of an L-GRIN microlens.

Structure and Principle of L-GRIN Micrlens
The schematic of the L-GRIN microlens designed is shown in Figure 1.The structure of the microlens as shown in Figure 1a consists of the microfluid chamber, core inlet, cladding inlet and outlet.The core liquid and the cladding liquid are injected through the core inlet and cladding inlet, respectively, which flow out through the outlet.H in_core , W in_core , H in_clad , W in_clad , H out , W out represent the height and width of the core inlet, cladding inlet and outlet, respectively.The main part of the L-GRIN microlens is a micro-cylindrical chamber, where the diffusion and convection process of liquids occurs and the gradient refractive index profile appears.The xoy cross-sectional view of the L-GRIN microlens together with the substrate is also illustrated in Figure 1a, where d core which is designed to be 50 µm stands for the diameter of the core inlet.Similarly, d clad is the diameter of the cladding inlet and is designed to be 150 µm.In the chamber, ethylene glycolsolution (core liquid) is injected side-by-side into de-ionized (DI) water (cladding liquid) from the same direction as shown in Figure 1b.On account of the well-known manner in which the refractive index of liquids varies with the wavelength of incident light, in our simulation, the wavelength of incident light is set to be 500nm without considering the dispersion of light.The diffusion and convection process between ethylene glycol solution and DI water results in a 2D refractive index profile in the lens chamber, as shown in the left graph of Figure 1c.An axis-symmetric refractive index gradient is accordingly achieved in the xoy plane and the simulation result is given in the right graph of Figure 1c, which shows that the refractive index of the L-GRIN microlens increases continuously from the cladding region to the core optical axis.The maximum refractive index comes near the inlet, and under the influence of diffusion, the refractive index decreases along the flowing direction of the liquids.
2015, 6ines2015, 6, page-page 3 with the wavelength of incident light, in our simulation, the wavelength of incident light is set to be 500nm without considering the dispersion of light.The diffusion and convection process between ethylene glycol solution and DI water results in a 2D refractive index profile in the lens chamber, as shown in the left graph of Figure 1c.An axis-symmetric refractive index gradient is accordingly achieved in the xoy plane and the simulation result is given in the right graph of Figure 1c, which shows that the refractive index of the L-GRIN microlens increases continuously from the cladding region to the core optical axis.The maximum refractive index comes near the inlet, and under the influence of diffusion, the refractive index decreases along the flowing direction of the liquids.To simulate and optimize the refractive index profile and light propagation in the L-GRIN microlens under different conditions, the FEM and optical ray tracing method are adopted.The refractive index profile can be calculated by simulating the diffusion and convection process of liquids in the microfluidic chamber.When the refractive index of liquids is mainly affected by diffusion, the process can be described by Fick's second law: In which C is the solution concentration and t is the diffusion time.D represents the diffusion coefficient that may determine the diffusion speed between the core liquid and the cladding liquid.As shown in Equation (1), the solution concentration changes with time and space.However, in the To simulate and optimize the refractive index profile and light propagation in the L-GRIN microlens under different conditions, the FEM and optical ray tracing method are adopted.The refractive index profile can be calculated by simulating the diffusion and convection process of liquids in the microfluidic chamber.When the refractive index of liquids is mainly affected by diffusion, the process can be described by Fick's second law: In which C is the solution concentration and t is the diffusion time.D represents the diffusion coefficient that may determine the diffusion speed between the core liquid and the cladding liquid.As shown in Equation (1), the solution concentration changes with time and space.However, in the microfluidic chamber, the diffusion and the convection have a combined effect on the refractive index profile of the L-GRIN microlens.The convective-diffusive process in the chamber is therefore rewritten as [20]: where U is defined to be the average flow velocity in the chamber.For a steady state flow, the distribution of concentration in the device does not vary with time, which means that BC{Bt " 0.
Considering the symmetry of the lens structure, the convective-diffusive process can be simplified as: The concentration distribution for the full, developed, steady-state flow can be expressed as Equation ( 4) [21].
In the normalized coordinate system, x 1 " x{R, y 1 " y{R, where R is the radius of the microfluidic chamber.Then c 1 = c/C 0 is the normalized concentration, and C 0 stands for the initial concentration.Pe is called the Peclet number, which is a dimensionless number and represents the relative proportion of convection to diffusion.Pe is proportional to the laminar velocity across the transversal section of the microfluidic chamber and the length of the microfluidic chamber, and Pe is also inversely proportional to the diffusion coefficient.Pe is usually used to describe the coefficient involved in the convection and diffusion process for a common analysis.However, the effect caused by convection and diffusion is analyzed separately in the present paper in order to dynamically reshape the light beam of the L-GRIN microlens.
For the structure we designed, U " pQ core `Qclad q{R 2 π is defined as the average velocity, where Q core is the flow rate of the core liquid and Q core is the flow rate of the cladding liquid; r is the normalized initial boundary ratio between the core and cladding liquids and is expressed as [22]: αβ 2 pr ´1q 4 `r pr ´1q rr 2 `3r ´αpr ´1qpr ´4qsβ ´r4 " 0 where α " Q core {Q clad , β " µ core {µ clad , and µ core is the liquid viscosity of the core liquid and µ clad is the liquid viscosity of the cladding liquid.The initial parameters are set as c 1 | x 1 "0,0ďy 1 ăr"1 " 1, c 1 ˇˇx1 "0,y 1 "r " r and c 1 ˇˇx1 "0,rďy 1 ď1 " 0 , in which Bc 1 By 1 ˇˇˇx "8,0ďy 1 ď1 " 0 stands for the full mixing state at the outlet and Bc 1 By 1 ˇˇˇy "1 " 0 represents the full mixing state at the chamber walls.
Once the liquids for the core and the cladding liquids are determined, the liquid viscosity and diffusion coefficient are considered to be constant.The position of the focal point can be tuned along the x-direction by changing the refractive index profile, which can be achieved by adjusting the flow rate of the core or cladding liquid.According to Equation (2), if the average velocity U is large, the diffusion in the xoy plane can be neglected and only the refractive index profile in the yoz plane needs to be calculated, where the plane is perpendicular to the liquid flow direction.On the contrary, when the average velocity U is small, the diffusion in both x and y directions has to be discussed.The refractive index profile can be obtained according to the concentration profile [23].
where n(x', y') is the refractive index profile of the L-GRIN microlens in the x'oy' plane, and c 1 core and c 1 clad are the normalized concentration profile of the core and cladding liquids in the x'oy' plane, respectively.The focusing of the input light beam through the microlens is determined by the refractive index profile.
Because of the decisive effect of the convection-diffusion process on the refractive index profile of the L-GRIN microlens, the average velocity U and the diffusion coefficient D, which are also found to vary with the concentration C and temperature, are important parameters for tunable focusing.For example, the diffusion coefficients between DI water and ethylene glycol are 3.75 ˆ10 ´10 and 1.17 ˆ10 ´9 m 2 /s with different mass fractions of 0.025 and 0.95 for ethylene glycol [24].In addition, when the temperature varies from 30 to 50 ˝C, the diffusion coefficient changes from 3.15 ˆ10 ´10 to 6.45 ˆ10 ´10 m 2 /s with a fixed mass fraction of 0.8 for the ethylene glycol [25].Due to the temperature of liquids usually being stationary during the convection-diffusion process, thediffusion coefficient D, concentration C and average velocity U of liquids are the main influence factors for adjusting the focusing performance of the L-GRIN microlens.It was found by simulations as follows that an adjustable focal length could be achieved under a relatively lower average velocity while the size of the focal spot had to be tuned at a high enough average velocity.Therefore, both the focal length and the focal spot of the output beam can be adjusted by carefully controlling the average velocity of the liquids.Suppose that the DI water and ethylene glycol solution are chosen as the cladding and the core liquids, respectively, and two cases of the average velocity with relative slip and without relative slip between the core and cladding inlet are disscused.The subsequent discussion will be focused on how the adjustment of the output beam will be implemented by adjusting these influence factors.
In order to validate our numerical simulation, the comparison between the observed light propagation in Reference [23] and our simulation results was illustrated in Figure 2. Figure 2a shows light propagation in the optofluidic waveguides under different fluidic conditions, and Figure 2b shows our simulation results at the same fluidic conditions.When driven at a highflow rate of Q clad = 50 µL¨min ´1 (Pe = 0.07), the light propagates in a straight line and shows no focusing effect.In contrast, at a low flow rate of Q clad = 0.5 µL¨min ´1, the light exhibits curved-ray trajectories and converges repeatedly for different core widths, r = 1/6, 1/3 and 1/2 (Pe = 0.0005, 0.0007 and 0.001, respectively) [23].
2015, 6ines2015, 6, page-page 5 where n(x', y') is the refractive index profile of the L-GRIN microlens in the x'oy' plane, and core c′ and clad c′ are the normalized concentration profile of the core and cladding liquids in the x'oy' plane, respectively.The focusing of the input light beam through the microlens is determined by the refractive index profile.
Because of the decisive effect of the convection-diffusion process on the refractive index profile of the L-GRIN microlens, the average velocity U and the diffusion coefficient D, which are also found to vary with the concentration C and temperature, are important parameters for tunable focusing.For example, the diffusion coefficients between DI water and ethylene glycol are 3.75 × 10 −10 and 1.17 × 10 −9 m 2 /s with different mass fractions of 0.025 and 0.95 for ethylene glycol [24].In addition, when the temperature varies from 30 to 50 °C, the diffusion coefficient changes from 3.15 × 10 −10 to 6.45 × 10 −10 m 2 /s with a fixed mass fraction of 0.8 for the ethylene glycol [25].Due to the temperature of liquids usually being stationary during the convection-diffusion process, thediffusion coefficient D, concentration C and average velocity U of liquids are the main influence factors for adjusting the focusing performance of the L-GRIN microlens.It was found by simulations as follows that an adjustable focal length could be achieved under a relatively lower average velocity while the size of the focal spot had to be tuned at a high enough average velocity.Therefore, both the focal length and the focal spot of the output beam can be adjusted by carefully controlling the average velocity of the liquids.Suppose that the DI water and ethylene glycol solution are chosen as the cladding and the core liquids, respectively, and two cases of the average velocity with relative slip and without relative slip between the core and cladding inlet are disscused.The subsequent discussion will be focused on how the adjustment of the output beam will be implemented by adjusting these influence factors.
In order to validate our numerical simulation, the comparison between the observed light propagation in Reference [23] and our simulation results was illustrated in Figure 2. Figure 2a shows light propagation in the optofluidic waveguides under different fluidic conditions, and Figure 2b shows our simulation results at the same fluidic conditions.When driven at a highflow rate of Qclad = 50 μL•min −1 (Pe = 0.07), the light propagates in a straight line and shows no focusing effect.In contrast, at a low flow rate of Qclad = 0.5 μL•min −1 , the light exhibits curved-ray trajectories and converges repeatedly for different core widths, r = 1/6, 1/3 and 1/2 (Pe = 0.0005, 0.0007 and 0.001, respectively) [23].[23].At a high flow rate of Q clad = 50 µL¨min ´1, the light is confined in the core due to the step-index distribution.At Q clad = 0.5 µL¨min ´1, the diffusion-induced gradient of the refractive index causes the light to repeatedly merge.With increasing r, the length of the first section becomes larger, for example P 1/6 = 270 µm, P 1/3 = 300 µm and P 1/2 = 340 µm.The focusing period also increases; for example, the first period for P 1/3 is 300 µm, the second is extended to 490 µm and the third is 590 µm.(Scale bar equals 300 µm.)(b)Simulated light propagation under the same flow conditions.With increasing r, the length of the first section becomes larger, for example, P 1/6 = 300 µm, P 1/3 = 360 µm and P 1/2 = 410 µm.The focusing period also increases; for example, the first period for P 1/3 is 360 µm, the second is extended to 540 µm and the third is 680 µm.
Similarly, the length of the first section increases with the core width and decreases almost linearly with the flow rate of the core fluid in both figures.However, when examining the light propagation pattern in detail, it is found that the self-focusing period is chirped such that the focusing period increases.For instance, as shown in Figure 2a, the first period for P 1/3 is 300 µm, the second is extended to 490 µm, and the third is 590 µm, owing to the diminishing bidirectional gradient contrast downstream as a result of diffusion [23].Furthermore, the first period for P 1/3 is 360 µm, the second is extended to 540 µm, and the third is 680 µm in Figure 2b.The small difference in focus length maybe due to the difference between two structures, and the change trends are the same.The match between the experimental results in Reference [23] and our simulation results illustrated that the follow-up simulations could be trusted.

TheRefractive Index Profile of the L-GRIN Microlens
In order to form the gradient refractive index profile in the lens chamber, the ethylene glycolsolution (n core = 1.432) and DI water (n clad = 1.332) are co-injected into the lens chamber from the same direction.Once contacting DI water, the ethylene glycol starts to diffuse from the ethylene glycol solution into the DI water.The average velocity U is calculated according to U " pQ core `Qclad q{R 2 π.For the purpose of maintaining the same flow rate of the ethylene glycol solution and the DI water (supposingthere is norelative slip between the core and cladding liquids), Q clad {Q core " 8 should be satisfied because the cross-sectional area of the cladding inlet is eight times as large as that of the core inlet.In order to illustrate the influence of the diffusion coefficient D on the refractive index profile of the L-GRIN microlens, the refractive index profiles in the xoy plane were simulated at two different diffusion coefficients (D = 1 ˆ10 ´9 m 2 /s and D = 4 ˆ10 ´10 m 2 /s) as shown in Figure 3. Figure 3a,b indicate that the variation of the diffusion coefficient resulting from the variation of the environmental conditions such as temperature has a strong influence on the refractive index profile, and the diffusion phenomenon becomes noticeable with the increase of the diffusion coefficient.This illustrates the fact that the diffusion is a highly activated process for the refractive index profile, provided the flow rates of the liquids are the same.
2015, 6ines2015, 6, page-page 6 Similarly, the length of the first section increases with the core width and decreases almost linearly with the flow rate of the core fluid in both figures.However, when examining the light propagation pattern in detail, it is found that the self-focusing period is chirped such that the focusing period increases.For instance, as shown in Figure 2a, the first period for P1/3 is 300 μm, the second is extended to 490 μm, and the third is 590 μm, owing to the diminishing bidirectional gradient contrast downstream as a result of diffusion [23].Furthermore, the first period for P1/3 is 360 μm, the second is extended to 540 μm, and the third is 680 μm in Figure 2b.The small difference in focus length maybe due to the difference between two structures, and the change trends are the same.The match between the experimental results in Reference [23] and our simulation results illustrated that the follow-up simulations could be trusted.

TheRefractive Index Profile of the L-GRIN Microlens
In order to form the gradient refractive index profile in the lens chamber, the ethylene glycolsolution (ncore = 1.432) and DI water (nclad = 1.332) are co-injected into the lens chamber from the same direction.Once contacting DI water, the ethylene glycol starts to diffuse from the ethylene glycol solution into the DI water.The average velocity U is calculated according to . For the purpose of maintaining the same flow rate of the ethylene glycol solution and the DI water (supposingthere is norelative slip between the core and cladding liquids), clad core 8 Q Q = should be satisfied because the cross-sectional area of the cladding inlet is eight times as large as that of the core inlet.In order to illustrate the influence of the diffusion coefficient D on the refractive index profile of the L-GRIN microlens, the refractive index profiles in the xoy plane were simulated at two different diffusion coefficients (D = 1×10 −9 m 2 /s and D = 4×10 −10 m 2 /s) as shown in Figure 3. Figure 3a,b indicate that the variation of the diffusion coefficient resulting from the variation of the environmental conditions such as temperature has a strong influence on the refractive index profile, and the diffusion phenomenon becomes noticeable with the increase of the diffusion coefficient.This illustrates the fact that the diffusion is a highly activated process for the refractive index profile, provided the flow rates of the liquids are the same.

The Adjustment of the Focal Lengthof the L-GRIN Microlens
It is known that a relatively lower average velocity of liquids offers remarkable diffusion effect, and the focus with an adjustable focal length of the L-GRIN microlens forms by changing the diffusion process.The diffusion process of the liquid is influenced by the average velocity U of liquids and the diffusion coefficient D which is affected simultaneously by the concentration C of the solution.Therefore, the focal length of the L-GRIN microlens can be tuned by adjusting the concentration and the flow rateof the liquids.Hereinafter, the concentration of the ethylene glycol solution is altered through changing the mass fraction of the ethylene glycol solution.

The Adjustment of the Focal Lengthof the L-GRIN Microlens
It is known that a relatively lower average velocity of liquids offers remarkable diffusion effect, and the focus with an adjustable focal length of the L-GRIN microlens forms by changing the diffusion process.The diffusion process of the liquid is influenced by the average velocity U of liquids and the diffusion coefficient D which is affected simultaneously by the concentration C of the solution.Therefore, the focal length of the L-GRIN microlens can be tuned by adjusting the concentration and the flow rateof the liquids.Hereinafter, the concentration of the ethylene glycol solution is altered through changing the mass fraction of the ethylene glycol solution.

Tuning Mass Fraction of Ethylene Glycol Solution
Because the diffusion coefficient changes with the concentration of the core liquid, the refractive index profile of the L-GRIN microlens varies with the mass fraction of the ethylene glycol solution.As the diffusion progresses, the solution flowing into different regions has a different concentration, and the diffusion coefficient is, therefore, also different.However, according to Equation ( 4), compared with the concentrated ethylene glycol solution, the diluted ethylene glycol solution will have an obvious gradient refractive index profile whose diffusion coefficient is larger and can be considered approximately constant during the diffusion process.To facilitate the analysis, smaller mass fractions of ethylene glycol ranging from 0.05 to 0.4 were selected and the focal length of the L-GRIN microlens with different mass fractions was calculated while the flow rate of the liquids was kept the same (Q core = 1 ˆ10 3 pL/s, Q clad = 8 ˆ10 3 pL/s).Firstly, we simulated the refractive index profiles of the L-GRIN microlens with the mass fraction increasing from 0.05 to 0.4 by 0.05, and simulation results showed that the sharper peak profile of the refractive index could be achieved with the increase of the mass fraction.Figure 4a shows the refractive index profile with the mass fraction of 0.2, and the cross-sectional refractive index profiles at five different locations (x = 50, 100, 150, 200 and 250µm) are shown in Figure 4b.The beam transmission and focusing process can be simulated by using the ray tracing method based on the refractive index profile, which is shown in Figure 4c.The focal length is the distance from the outlet to the focal point, and Figure 4d shows that the focal length of the L-GRIN microlens can be tuned with different mass fractions of the ethylene glycol solution.When the flow rates of the core and cladding liquid are fixed, the focal length of the L-GRIN microlens varies from 942 to 11 µm with a mass fraction of ethylene glycol ranging from 0.05 to 0.4.Therefore, controlling the mass fraction of the core liquid is a useful way to adjust the focal length of the device.

Tuning Mass Fraction of Ethylene Glycol Solution
Because the diffusion coefficient changes with the concentration of the core liquid, the refractive index profile of the L-GRIN microlens varies with the mass fraction of the ethylene glycol solution.As the diffusion progresses, the solution flowing into different regions has a different concentration, and the diffusion coefficient is, therefore, also different.However, according to Equation ( 4), compared with the concentrated ethylene glycol solution, the diluted ethylene glycol solution will have an obvious gradient refractive index profile whose diffusion coefficient is larger and can be considered approximately constant during the diffusion process.To facilitate the analysis, smaller mass fractions of ethylene glycol ranging from 0.05 to 0.4 were selected and the focal length of the L-GRIN microlens with different mass fractions was calculated while the flow rate of the liquids was kept the same (Qcore = 1 × 10 3 pL/s, Qclad = 8 × 10 3 pL/s).Firstly, we simulated the refractive index profiles of the L-GRIN microlens with the mass fraction increasing from 0.05 to 0.4 by 0.05, and simulation results showed that the sharper peak profile of the refractive index could be achieved with the increase of the mass fraction.Figure 4a shows the refractive index profile with the mass fraction of 0.2, and the cross-sectional refractive index profiles at five different locations (x = 50, 100, 150, 200 and 250μm) are shown in Figure 4b.The beam transmission and focusing process can be simulated by using the ray tracing method based on the refractive index profile, which is shown in Figure 4c.The focal length is the distance from the outlet to the focal point, and Figure 4d shows that the focal length of the L-GRIN microlens can be tuned with different mass fractions of the ethylene glycol solution.When the flow rates of the core and cladding liquid are fixed, the focal length of the L-GRIN microlens varies from 942 to 11 μm with a mass fraction of ethylene glycol ranging from 0.05 to 0.4.Therefore, controlling the mass fraction of the core liquid is a useful way to adjust the focal length of the device.

Tuning the Flow Rate of Liquids
The refractive index profile of the L-GRIN microlens can also be tuned by changing the flow rate of the liquids.In order to analyze the influence of the flow rate of the liquids on the focal length of the L-GRIN microlens, the ethylene glycol solution is co-injected with DI water into the lens chamber at the same flow rate ranging from 0.5 × 10 3 to 5 × 10 3 pL/s.In the analysis, we suppose there is no relative slip between the ethylene glycol solution and the DI water.Additionally, for

Tuning the Flow Rate of Liquids
The refractive index profile of the L-GRIN microlens can also be tuned by changing the flow rate of the liquids.In order to analyze the influence of the flow rate of the liquids on the focal length of the L-GRIN microlens, the ethylene glycol solution is co-injected with DI water into the lens chamber at the same flow rate ranging from 0.5 ˆ10 3 to 5 ˆ10 3 pL/s.In the analysis, we suppose there is no relative slip between the ethylene glycol solution and the DI water.Additionally, for simplifying the calculation, simulations were started by setting D " 8 ˆ10 ´10 m 2 {s and µ = 1 ˆ10 ´3 Pa¨s with a fixed mass fraction of 0.3.
We simulated the refractive index profiles with the flow rate increasing from 0.5 ˆ10 3 to 5 ˆ10 3 pL/s by 0.5 ˆ10 3 pL/s, and the results demonstrated that the diffusion was no longer dominant with the increase of the flow rate.Figure 5a shows the refractive index profile with the flow rate of 4 ˆ10 3 pL/s, and the cross-sectional refractive index profiles at five different locations (x = 50, 100, 150, 200 and 250 µm) are shown in Figure 5b. Figure 5c shows the focal length of the L-GRIN microlens is adjusted with the different flow rate.According to our calculations, the focal length decreases observably from 128 to 8 µm when the flow rate varies from 0.5 ˆ10 3 to 5 ˆ10 3 pL/s.However, when the flow rate is larger than 5 ˆ10 3 pL/s, the converging of the light beam caused by focusing becomes more and more unobvious.Therefore, continuously tuning the flow rate of the liquids in a certain range provides a tunable focal length when the mass fraction of the ethylene glycol is kept constant.In general, the focal length of the lens is related to the lens radius, the refractive index difference between the lens center and border, and the lens thickness; in addition, it is dependent on the wavelength of the incident light which will change the refractive index profile of the lens.Thus, a shorter wavelength of incident light (when other parameters are constant) results in larger refractive index contrast, which causes light to bend toward the lens axis more significantly and leads to the decreased focal length.

The Adjustment of the Focal Spot of the L-GRIN Microlens
The focus becomes a tunable focal spot at a relatively higher average velocity.The adjustment of the focal spot of the L-GRIN microlens can also be achieved by changing the flow rate of the liquids when the average velocity is higher.The case of no relative slipbetween the core and cladding liquids with a lower average velocity has been analyzed above.Analyzing a more sophisticated case will be helpful to better understand the influence of flow rate on the refractive index profile.For instance, the ethylene glycol solution is injected in different flow rates with DI water, andthere is a relative slip between the core and cladding liquids.For analyzing the influence of the relative slip between the core and cladding liquids on the refractive index profile, the cladding flow rate was kept constant (40 ˆ10 3 pL/s) and the flow rate of the core inlet varied from 2 ˆ10 3 pL/s to 50 ˆ10 3 pL/s by 5 ˆ10 3 pL/s.Figure 6a shows the refractive index profile with the core inlet flow rate of 25 ˆ10 3 pL/s, and the cross-sectional refractive index profiles at five different locations (x = 50, 100, 150, 200 and 250 µm) are shown in Figure 6b.It is found that there is little variation of the refractive index in central regions with different cross-sectional lengths along the x axis.The width in the central region along the y axis, in which the refractive index keeps almost constant, is defined as the core width.When the core inlet flow rate is Q core = 25 ˆ10 3 pL/s, the core width of the refractive index is 24 µm, as shown in Figure 6b.The constant refractive index in the central region hasno obvious effect on the focusing process;thus, the size of the output beam spot can be determined mainly by the core width.The curve with circle in Figure 6c shows the relationship between the core inlet flow rate and the core width.Additionally, the curve width dot shows the changing trend of half the width of the refrative index profile with different Q core .When the core inlet flow rate remainsconstant, with a relatively low core inlet flow rate (lower than 10 ˆ10 3 pL/s),the core width approachs zero.In that case, the L-GRIN microlens can theoretically focus the beam onto one point when the ratio of the core inlet flow rate to the cladding inlet flow rate is lower than 0.25.However, the core width increases along with the increaseof the core inlet flow rate (higher than 10 ˆ10 3 pL/s); that is, when increasing the ratio of the flow rates of the ethylene glycol solution to the DI water, the size of the output spot may continue to grow.Thus, controlling the ratio of the core inlet flow rate to the cladding inlet flow rate is a useful way to dynamically adjust the size of the light beam in microscale.
2015, 6ines2015, 6, page-page of the relative slip between the core and cladding liquids on the refractive index profile, the cladding flow rate was kept constant (40 × 10 3 pL/s) and the flow rate of the core inlet varied from 2 × 10 3 pL/s to 50 × 10 3 pL/s by 5 × 10 3 pL/s.Figure 6a shows the refractive index profile with the core inlet flow rate of 25 × 10 3 pL/s, and the cross-sectional refractive index profiles at five different locations (x = 50, 100, 150, 200 and 250 μm) are shown in Figure 6b.It is found that there is little variation of the refractive index in central regions with different cross-sectional lengths along the x axis.The width in the central region along they axis, in which the refractive index keeps almost constant, is defined as the core width.When the core inlet flow rate is Qcore = 25×10 3 pL/s, the core width of the refractive index is 24 μm, as shown in Figure 6b.The constant refractive index in the central region hasno obvious effect on the focusing process;thus, the size of the output beam spot can be determined mainly by the core width.The curve with circle in Figure 6c shows the relationship between the core inlet flow rate and the core width.Additionally,the curve width dot shows the changing trend of half the width of the refrative index profile with different Qcore.When the core inlet flow rate remainsconstant, with a relatively low core inlet flow rate (lower than 10 × 10 3 pL/s),the core width approachs zero.In that case, the L-GRIN microlens can theoretically focus the beam onto one point when the ratio of the core inlet flow rate to the cladding inlet flow rate is lower than 0.25.However, the core width increases along with the increaseof the core inlet flow rate (higher than 10 × 10 3 pL/s); that is, when increasing the ratio of the flow rates of the ethylene glycol solution to the DI water,the size of the output spot may continue to grow.Thus, controlling the ratio of the core inlet flow rate to the cladding inlet flow rate is a useful way to dynamically adjust the size of the light beam in microscale.Figure 6b reveals that the cross-sectional refractive index profile in the yoz plane changes slightly even though x is different, which is fitted using six-degree polynomial functions.The optimal parameters of the functions areselected through the test of the partial fitting square sum on the basis of the variance analysis of the fitting equation.For example, the fitted refractive index profile at x = 125 μm with Qcore = 2.5 × 10 4 pL/s and Qcore = 4 × 10 4 pL/s can be expressed as: Figure 6b reveals that the cross-sectional refractive index profile in the yoz plane changes slightly even though x is different, which is fitted using six-degree polynomial functions.The optimal parameters of the functions areselected through the test of the partial fitting square sum on the basis of the variance analysis of the fitting equation.For example, the fitted refractive index profile at x = 125 µm with Q core = 2.5 ˆ10 4 pL/s and Q core = 4 ˆ10 4 pL/s can be expressed as: n " 1.1 ˆ10 ´12 s 6 ´2.2 ˆ10 ´13 s 5 ´6.45 ˆ10 ´9s 4 ´2.73 ˆ10 ´10 s 3 `3.62 ˆ10 ´6s 2 `1.99 ˆ10 ´6s `1.3533 (7) where s " a y 2 `z2 This is a six-order polynomial of s, which represents the non-linear relation between the refractive index n of the liquids and the chamber position in the x direction.Figure 7a shows that the fitting refractive index profile matches well with the simulated data.For analyzing the beam focusing effect of the L-GRIN microlens, the beam transmission and focusing process was simulated using the ray tracing method according to the refractive index profile, which is shown in Figure 7b.A beam spot with the diameter of 23.5 µm and the focal length of 235.3 µm was achieved.The ability of the tunable focusing of the beam is significant for a wide range of applications in integrated optics and lab-on-a-chip systems.Figure 7a shows that the fitting refractive index profile matches well with the simulated data.For analyzing the beam focusing effect of the L-GRIN microlens, the beam transmission and focusing process was simulated using the ray tracing method according to the refractive index profile, which is shown in Figure 7b.A beam spot with the diameter of 23.5 μm and the focal length of 235.3 μm was achieved.The ability of the tunable focusing of the beam is significant for a wide range of applications in integrated optics and lab-on-a-chip systems.In the L-GRIN microlens, the diameter of the focal spot is affected by the focal length and the wavelength of incident light.For a given incident light, one can conclude that the focal spot is the smallest when the focal length is the shortest.As the wavelength increases, the focal length is gradually increased, and the focal spot is increased.
These discussions above are all for a steady-state flow.However, the response time of the focal length and the focal spot change are important parameters for the stability of an adaptive lens operation, which is mainly dependent on the diffusion speed of two inlet flows.According to Equation (1), the response time varies with the concentration C and the diffusion coefficient D of the liquids, which are also found to be affected by the temperature, the average velocity U and the viscosity of the liquids, the structure and materials of the chamber, etc.For the fixed inlet liquids and the chamber, as the temperature increases, the response time decreases.However, increasing the average velocity of the liquids would lend to an unfavorable slow response time because the contacting time of the two liquids is shortened.Therefore, for practical applications we need to choose the suitable liquids and optimize the microlens parameters; meanwhile, the microlens material should be carefully chosen.

Conclusions
We designed an in-plane L-GRIN microlens that can realize two-dimensional light beam In the L-GRIN microlens, the diameter of the focal spot is affected by the focal length and the wavelength of incident light.For a given incident light, one can conclude that the focal spot is the smallest when the focal length is the shortest.As the wavelength increases, the focal length is gradually increased, and the focal spot is increased.
These discussions above are all for a steady-state flow.However, the response time of the focal length and the focal spot change are important parameters for the stability of an adaptive lens operation, which is mainly dependent on the diffusion speed of two inlet flows.According to Equation (1), the response time varies with the concentration C and the diffusion coefficient D of the liquids, which are also found to be affected by the temperature, the average velocity U and the viscosity of the liquids, the structure and materials of the chamber, etc.For the fixed inlet liquids and the chamber, as the temperature increases, the response time decreases.However, increasing the average velocity of the liquids would lend to an unfavorable slow response time because the

Figure 2 .
Figure 2. (a) Observed light propagation under different flow conditions[23].At a high flow rate of Qclad = 50 μL•min −1 , the light is confined in the core due to the step-index distribution.At Qclad = 0.5 μL•min −1 , the diffusion-induced gradient of the refractive index causes the light to repeatedly merge.With increasing r, the length of the first section becomes larger, for example P1/6 = 270 μm, P1/3 = 300 μm and P1/2 = 340 μm.The focusing period also increases; for example, the first period for P1/3 is 300 μm, the second is extended to 490 μm and the third is 590 μm.(Scale bar equals 300 μm.)(b) Simulated light propagation under the same flow conditions.With increasing r, the length of the first section becomes larger, for example, P1/6 = 300 μm, P1/3 = 360 μm and P1/2 = 410 μm.The focusing period also increases; for example, the first period for P1/3 is 360 μm, the second is extended to 540 μm and the third is 680 μm.

Figure 2 .
Figure 2. (a) Observed light propagation under different flow conditions[23].At a high flow rate of Q clad = 50 µL¨min ´1, the light is confined in the core due to the step-index distribution.At Q clad = 0.5 µL¨min ´1, the diffusion-induced gradient of the refractive index causes the light to repeatedly merge.With increasing r, the length of the first section becomes larger, for example P 1/6 = 270 µm, P 1/3 = 300 µm and P 1/2 = 340 µm.The focusing period also increases; for example, the first period for P 1/3 is 300 µm, the second is extended to 490 µm and the third is 590 µm.(Scale bar equals 300 µm.)(b)Simulated light propagation under the same flow conditions.With increasing r, the length of the first section becomes larger, for example, P 1/6 = 300 µm, P 1/3 = 360 µm and P 1/2 = 410 µm.The focusing period also increases; for example, the first period for P 1/3 is 360 µm, the second is extended to 540 µm and the third is 680 µm.

Figure 4 .
Figure 4. (a) The refractive index profile in the xoy plane with the mass fraction of 0.2; (b) The cross-sectional refractive index profiles at five different locations (x = 50, 100, 150, 200 and 250 μm) with Qcore = 1 × 10 3 pL/s, Q clad =8 × 10 3 pL/s; (c) Beam transmission and focusing process; (d) The focal length of the L-GRIN microlens with different mass fractions of the ethylene glycol solution.

Figure 4 .
Figure 4. (a) The refractive index profile in the xoy plane with the mass fraction of 0.2; (b) The cross-sectional refractive index profiles at five different locations (x = 50, 100, 150, 200 and 250 µm) with Q core = 1 ˆ10 3 pL/s, Q clad = 8 ˆ10 3 pL/s; (c) Beam transmission and focusing process; (d) The focal length of the L-GRIN microlens with different mass fractions of the ethylene glycol solution.

Figure 5 .
Figure 5. (a) The refractive index profile in the xoy plane with the flow rate of 4 ˆ10 3 pL/s; (b) The cross-sectional refractive index profiles at five different locations (x = 50, 100, 150, 200 and 250 µm) with Q clad = 8 ˆ10 3 pL/s, Q core = 1 ˆ10 3 pL/s; (c) The focal length of the L-GRIN microlens with different flow rates of solutions.

Figure 6 .
Figure 6.(a) The refractive index profile in the xoy plane with Q core = 25 ˆ10 3 pL/s; (b) The cross-sectional refractive index profiles at five different locations (x = 50, 100, 150, 200 and 250 µm) at Q core = 25 ˆ10 3 pL/s; (c) The variation of core width and half width with different Q core .
This is a six-order polynomial of s, which represents the non-linear relation between refractive index n of the liquids and the chamber position in the x direction.

Figure 7 .
Figure 7. (a) The fitting refractive index profile of the simulated data and six-degree polynomial at x = 125 µm with Q core = 2.5 ˆ10 4 and Q core = 4 ˆ10 4 pL/s; (b) Beam transmission and focusing process.