Numerical Simulation of the Photobleaching Process in Laser-Induced Fluorescence Photobleaching Anemometer

At present, a novel flow diagnostic technique for micro/nanofluidics velocity measurement—laser-induced fluorescence photobleaching anemometer (LIFPA)—has been developed and successfully applied in broad areas, e.g., electrokinetic turbulence in micromixers and AC electroosmotic flow. Nevertheless, in previous investigations, to qualitatively reveal the dynamics of the photobleaching process of LIFPA, an approximation of uniform laser distribution was applied. This is different from the actual condition where the laser power density distribution is normally Gaussian. In this investigation, we numerically studied the photobleaching process of fluorescent dye in the laser focus region, according to the convection–diffusion reaction equation. The profiles of effective dye concentration and fluorescence were elucidated. The relationship between the commonly used photobleaching time constant obtained by experiments and the photochemical reaction coefficient is revealed. With the established model, we further discuss the effective spatial resolution of LIFPA and study the influence of the detection region of fluorescence on the performance of the LIFPA system. It is found that at sufficiently high excitation laser power density, LIFPA can even achieve a super-resolution that breaks the limit of optical diffraction. We hope the current investigation can reveal the photobleaching process of fluorescent dye under high laser power density illumination, to enhance our understanding of fluorescent dynamics and photochemistry and develop more powerful photobleaching-related flow diagnostic techniques.


Introduction
Laser-induced fluorescence photobleaching anemometer (LIFPA) is a velocity measurement technique developed by Wang in 2005 [1] and primarily for micro/nanofluidics. It measures flow velocity due to the photobleaching nature of fluorescent dye under a high laser power density. When fluorescent dye molecules are excited by a wavelength-matched laser, either the structures of molecules could be damaged due to the photo instability of the dye molecules [2], or the excited states could be quenched by dye-dye or solvent interactions, when the power density of the excitation beam exceeds a certain threshold. The photobleaching results in a decrease in the effective concentration (C) of fluorescent dye and a weakening fluorescent signal with increasing irradiation time [3]. The faster the fluorescent molecules pass through the excitation region in a flow field with uniform fluorescent molecules, the higher the effective concentration of fluorescent dye, and the higher the fluorescent intensity collected [4]. By establishing the monotonic relationship [5] between fluorescence intensity and flow velocity, the velocity of the flow field can be calculated by detecting the fluorescence signal intensity in the spot area [6]. This technique has been successfully applied to the velocity measurement of complex flow fields such as linearly and nonlinearly oscillating electroosmotic flow [7][8][9][10], and microelectrokinetic turbulence [11][12][13][14][15].
where I f represents the overall fluorescence intensity of the laser focal area, I f 0 is the initial I f at t = 0, t is the bleaching time that the fluorescent molecules pass through the region of laser focus diameter d f , τ is the photobleaching time constant, U represents flow velocity. Then, Zhao et al. [16] studied the photobleaching process and established a theoretical formula to estimate the temporal resolution of LIFPA. Uτ where Ψ s is the weight function of U, and I f ,b is the background fluorescence intensity when the dye is approximately fully photobleached. However, to simplify the analysis, they assumed the local fluorescence intensity (I f ,local ) distribution in the focal area to be uniform, which could lead to deviation from the actual case, where I f ,local normally has a Gaussian distribution. Moreover, the relationships among the complex fluorescent dynamics, the photobleaching reaction process and the superficial and measurable experimental parameters (e.g., τ) have not been established. This could be a major obstacle to the development of the LIFPA technique. In this investigation, the photobleaching process is studied as a photochemical reaction, which can be further numerically analyzed by a convection-diffusion reaction process. We first established a relationship between the photobleaching time constant and photochemical reaction coefficient, by comparing experiments and numerical simulations, to reveal the relationships between the photobleaching reaction process and experimental parameters. Then, the influence of flow velocity and laser intensity on the effective concentration and fluorescent intensity was numerically investigated. At last, the spatial resolution of LIFPA was revisited to better understand the photobleaching phenomenon in flowing liquids and to develop new flow diagnostic techniques.

LIFPA Photobleaching Model
The photobleaching process of the fluorescent solution in a steady microchannel flow can be described by a convection-diffusion reaction equation [17][18][19] as ∂C ∂t In Equation (3), C represents the concentration of the fluorescent solution, → u is the flow velocity vector, D is the diffusivity of the dye molecule in the solution, and R c represents the reaction term according to photobleaching, which is determined by the laser intensity distribution, C and a factor k b . By first-order approximation, R c can be expressed as where k b represents the photobleaching rate of dye molecules per unit excitation intensity at a specific excitation wavelength, with a unit of cm 2 /(W·s). It represents how fast the fluorescent dye can be photobleached. The higher the k b , the faster the photobleaching is. I is the intensity of the excitation laser; in general, I has a Gaussian distribution, shown as I max is the peak value of I, and G(x, y, z) is a 3D Gaussian function, shown as where σ x , σ y and σ z represent the standard deviations of laser intensity in x, y and z directions, respectively. Substituting Equations (4)-(6) into Equation (3), with dimensional analysis by where C 0 , L and U 0 are characteristic concentration, length scale and velocity, respectively, we have where P e = U 0 L/D is the Péclet number which represents the ratio of the convection rate to the diffusion rate. In microfluidic applications, normally the disturbance of microflow has the smallest time scale of O(10 −4 s) or higher. In contrast, to obtain good sensitivity and resolution, LIFPA has a fast photobleaching process where τ is of O(10 −6 s), which is approaching the lower limit of photobleaching lifetime, which is normally in the range of nanoseconds to seconds [20][21][22][23][24], or even longer. Meanwhile, the resident time of dye in the laser focus is O(10 −4 s) or lower. Therefore, the time scale of flow is equivalent to or larger than the resident time and τ. A quasi-steady state of the photobleaching process of dye can be approximated in the laser focus region. Thus, ∂C * /∂t * = 0. The LIFPA photobleaching model under a quasi-steady state is obtained as where Z pc = KL/U 0 is a dimensionless number to evaluate the ratio between the photobleaching rate and convective transport rate of dye molecules. The larger the Z pc , the stronger the local photobleaching is. Otherwise, the convective transport of the concentration of the fluorescent solution is more important. K = k b I max is the photochemical reaction coefficient, which represents the photobleaching rate of dye molecules in this investigation. It is a crucial parameter to characterize the photobleaching process through the convection-diffusion reaction equation in COMSOL simulation. Although K has a unit of 1/s, which is the same as 1/τ, K is not 1/τ. One of the focuses of this investigation is to find the relationship between K and τ. The characteristic concentration C 0 was eliminated during the simplification process. This indicates that increasing or decreasing the dye concentration C does not affect the first-order photobleaching process, if in the absence of light absorption. Finally, the local fluorescent intensity I f ,local can be calculated accordingly as [25] with q f being a fluorescence quantum yield factor.

LIFPA System
The experimental LIFPA system is developed based on a confocal microscope, as shown in Figure 1a. The excitation laser is a 405 nm continuous wave laser (MDL-III-405-500, CNI). The excitation beam is first controlled by an acousto-optic modulator (AOM, 1206C-2-1002, Isomet) for on-off switch. Then, it passes through a spatial pinhole filter (SLF, SFB-16DM, OptoSigma) to improve the beam quality. A diaphragm is used to filter out the zero-order and high-order diffraction spots and retain the first-order spot in the center.
After passing through a collimation lens, a collimated excitation beam with a Gaussian beam profile is obtained.
Subsequently, the excitation beam passes through the dichroic mirror (DM 1 ), which transmits the 405 nm laser and reflects the fluorescence around 480 nm. The excitation beam is then reflected into the objective lens (OL, Olympus PlanApo, 100 × NA 1.4 oil immersion Objective lens) by a mirror (M 1 ). The excitation beam is finally focused in the microchannel flow and illuminates the fluorescent dye. The fluorescent signal passes the optical path along M 1 , DM 1 , and a second mirror (M 2 ), and then passes through a band-pass filter (BP, 470/10 nm, OptoSigma) to eliminate background light noise. The fluorescent signal is collected by a lens (L 2 ) and focused on a pinhole PH 1 for spatial filtering. The filtered fluorescent signal is detected by a photon counter (H7421-40, Hamamatsu). The exposure time of fluorescent measurement, i.e., the sampling duration, is 0.5 ms. A time interval of 0.5 ms is applied between two samples. Therefore, the sampling rate is 1 kHz. Since the photo counter has a saturation count rate of 1.5 × 10 6 photons per second for linear measurement (equivalent to a maximum of 750 photons in the 0.5 ms sampling duration), and the fluorescence of LIFPA is not weak, we used a neutral density filter of OD2 in front of PH 1 , to reduce the fluorescent intensity and guarantee there is no nonlinear saturation of photon counter. In the experiments, the maximum photon number in a 0.5 ms sampling duration is below 533, which is in the linear range of the photon counter. center. After passing through a collimation lens, a collimated excitation beam with a Gaussian beam profile is obtained. Subsequently, the excitation beam passes through the dichroic mirror (DM1), which transmits the 405 nm laser and reflects the fluorescence around 480 nm. The excitation beam is then reflected into the objective lens (OL, Olympus PlanApo, 100 × NA 1.4 oil immersion Objective lens) by a mirror (M1). The excitation beam is finally focused in the microchannel flow and illuminates the fluorescent dye. The fluorescent signal passes the optical path along M1, DM1, and a second mirror (M2), and then passes through a bandpass filter (BP, 470/10 nm, OptoSigma) to eliminate background light noise. The fluorescent signal is collected by a lens (L2) and focused on a pinhole PH1 for spatial filtering. The filtered fluorescent signal is detected by a photon counter (H7421-40, Hamamatsu). The exposure time of fluorescent measurement, i.e., the sampling duration, is 0.5 ms. A time interval of 0.5 ms is applied between two samples. Therefore, the sampling rate is 1 kHz. Since the photo counter has a saturation count rate of 1.5 10 photons per second for linear measurement (equivalent to a maximum of 750 photons in the 0.5 ms sampling duration), and the fluorescence of LIFPA is not weak, we used a neutral density filter of OD2 in front of PH1, to reduce the fluorescent intensity and guarantee there is no nonlinear saturation of photon counter. In the experiments, the maximum photon number in a 0.5 ms sampling duration is below 533, which is in the linear range of the photon counter.

Microchannel and Solution Preparation
The experiment is carried out at the center of the microchannel. For precise positioning of the microchannel with a large travel distance, a 2D translation stage (PI M-521.DG, 1 μm accuracy) and a high-precision 3D piezo nanostage (PI 562.3CD, 1 nm accuracy) are

Microchannel and Solution Preparation
The experiment is carried out at the center of the microchannel. For precise positioning of the microchannel with a large travel distance, a 2D translation stage (PI M-521.DG, 1 µm accuracy) and a high-precision 3D piezo nanostage (PI 562.3CD, 1 nm accuracy) are applied together, as shown in Figure 1a. The structure of the microfluidic chip is shown in Figure 1b. It has three layers. The cover layer is made of 2 mm-thick acrylic with good transparency and mechanical strength. The middle layer, which is also the channel layer, is made of plastic sheets. The bottom plate is a low-fluorescence glass slide with a thickness of 130 µm. The chip is assembled layer-by-layer, to realize a microchannel with a rectangular cross section. The length, width, and height of the microchannel are 5 mm, 360 µm, and 90 µm, respectively. During experiments, the fluorescent solution of Coumarin 102 (C102, Sigma Aldrich) at a concentration of 0.1 mM is injected into the microchannel with a syringe pump to generate a laminar flow. The C102 solution is prepared by dissolving 2.5 mg of C102 powder with 95 mL of deionized water and 5 mL of methanol solution (Analytical Reagent, concentration > 99.5%). The C102 solution has a pH value of 7.82. Finally, after excitation, the photobleached fluorescent solution is drained from the outlet.

Numerical Simulation by COMSOL
In this paper, the photobleaching process of LIFPA is numerically investigated with COMSOL Multiphysics 4.3. The simulation is carried out in a 3D region with rectangular cross sections, as shown in Figure 2a. outlet.

Numerical Simulation by COMSOL
In this paper, the photobleaching process of LIFPA is numerically investigated with COMSOL Multiphysics 4.3. The simulation is carried out in a 3D region with rectangular cross sections, as shown in Figure 2a.
A laser focus region approximately simulating the laser beam is located at the center ( 0) of the computational region and has an axisymmetric structure with 3 parts. At the middle part is a cylindrical region with diameter 203 nm and height ℎ 800 nm, whose diameter is coincident with the full width at half maximum (FWHM) of the laser beam. The top and bottom parts are symmetric cones with diameter 5 μm and height ℎ 1 μm, respectively, to be coincident with the light field according to the high NA of the objective. In the 3D simulation, free tetrahedral elements with variable sizes are applied, as shown in Figure 2b. Since the size of the laser focus region of the excitation beam is only 203 nm, to have sufficient spatial resolution and smooth distribution of concentration field, the minimum size of the elements in the region is 0.5 nm. In the other regions, the minimum and maximum sizes of the elements are 100 nm and 530 nm, respectively. If we simulated the convection-diffusion reaction equation in the entire experimental microchannel, the computation burden is extraordinarily heavy. To save computation resources, in the simulation, we reduce the size of the computation region to 30 μm long, 10 μm wide and 10 μm high. The total number of free tetrahedral elements is 1311608. The photobleaching process is simulated with laminar flow and dilute matter transfer modules in COMSOL. In the laminar flow module, the flow is assumed to be steady, uni-directional and fully developed. Therefore, we have ⃗ for simplification. A noslip boundary condition [26,27] is been applied on the side walls ( ) as Since the sizes of the experimental microchannel are different from those of the computational region, to guarantee that the flow velocities at the laser focus regions of both experiments and simulations are the same, the flow rate in numerical simulation is different from that in experiments. In this investigation, we first calculate the flow velocity at the laser focus region (the third row of Table 1) through experimental flow rate (the first row of Table 1), according to Equation (12)  A laser focus region approximately simulating the laser beam is located at the center (x = y = z = 0) of the computational region and has an axisymmetric structure with 3 parts. At the middle part is a cylindrical region with diameter d f = 203 nm and height h f = 800 nm, whose diameter is coincident with the full width at half maximum (FWHM) of the laser beam. The top and bottom parts are symmetric cones with diameter d c ≈ 5 µm and height h c = 1 µm, respectively, to be coincident with the light field according to the high NA of the objective. In the 3D simulation, free tetrahedral elements with variable sizes are applied, as shown in Figure 2b. Since the size of the laser focus region of the excitation beam is only 203 nm, to have sufficient spatial resolution and smooth distribution of concentration field, the minimum size of the elements in the region is 0.5 nm. In the other regions, the minimum and maximum sizes of the elements are 100 nm and 530 nm, respectively. If we simulated the convection-diffusion reaction equation in the entire experimental microchannel, the computation burden is extraordinarily heavy. To save computation resources, in the simulation, we reduce the size of the computation region to 30 µm long, 10 µm wide and 10 µm high. The total number of free tetrahedral elements is 1,311,608.
The photobleaching process is simulated with laminar flow and dilute matter transfer modules in COMSOL. In the laminar flow module, the flow is assumed to be steady, unidirectional and fully developed. Therefore, we have → u = Ux for simplification. A no-slip boundary condition [26,27] is been applied on the side walls (Ω) as Since the sizes of the experimental microchannel are different from those of the computational region, to guarantee that the flow velocities at the laser focus regions of both experiments and simulations are the same, the flow rate in numerical simulation is different from that in experiments. In this investigation, we first calculate the flow velocity at the laser focus region (the third row of Table 1) through experimental flow rate (the first row of Table 1), according to Equation (12) [8,28,29]: where Q is the flow rate. −w/2 ≤ y ≤ w/2 and −h/2 ≤ z ≤ h/2. Then, the inlet flow rate in the simulation (the second row of Table 1) is calculated based on the flow velocity at the laser focus region (the third row of Table 1) in simulation, which is the same as in experiments. The outlet pressure is the circumstantial pressure. The flow is incompressible with a constant temperature.
In the dilute matter transfer module, the diffusion coefficient of molecules in the solution is D = 1 × 10 −9 m 2 /s. The initial fluorescent solution concentration C = 0.1 mM. At the side walls, no flux is present. The photobleaching of fluorescent molecules occurs only in the laser focus region in Figure 2. The photochemical reaction is stimulated by applying a reaction term based on Equations (4)-(6).

Direct Comparison between Experiments and Numerical Simulations
The relationships between U and fluorescence are established by experiments at three different laser powers (P) of 6.9, 11.8 and 18.2 mW. The corresponding I max are 2.96 × 10 7 , 5.05 × 10 7 and 7.80 × 10 7 W/cm 2 , respectively, as shown in Figure 3a. Each data point was calculated by averaging over 2 × 10 4 fluorescent signals, i.e., a 20 s time sequence of fluorescent signals. It is obvious that I f increases with the increasing U. Moreover, when the laser power density is small (e.g., I max = 2.96 × 10 7 W/cm 2 ), I f shows an obvious nonlinear increment with U. When I max is increased to 5.05 × 10 7 and 7.80 × 10 7 W/cm 2 , the nonlinear behavior gradually decreases and the curve becomes more linear, accompanied by a decreasing slope of the U ∼ I f curve. In these cases, velocity fluctuations can be more clearly distinguished by larger fluorescent fluctuations, as shown in Figure 3a. Thus, high laser power density can provide a better sensitivity of LIFPA measurement. three different laser powers ( ) of 6.9, 11.8 and 18.2 mW. The corresponding are 2.96 10 , 5.05 10 and 7.80 10 W/cm 2 , respectively, as shown in Figure 3a. Each data point was calculated by averaging over 2 10 fluorescent signals, i.e., a 20 s time sequence of fluorescent signals. It is obvious that increases with the increasing . Moreover, when the laser power density is small (e.g., 2.96 10 W/cm 2 ), shows an obvious nonlinear increment with . When is increased to 5.05 10 and 7.80 10 W/cm 2 , the nonlinear behavior gradually decreases and the curve becomes more linear, accompanied by a decreasing slope of the ~ curve. In these cases, velocity fluctuations can be more clearly distinguished by larger fluorescent fluctuations, as shown in Figure 3a. Thus, high laser power density can provide a better sensitivity of LIFPA measurement.  In order to compare ~ relationships at different laser power densities, we calculate the normalized fluorescence intensity ( ) with the maximum ( , ) and minimum ( , ) values as [2] , , , In order to compare U ∼ I f relationships at different laser power densities, we calculate the normalized fluorescence intensity ( I f ) with the maximum (I f ,max ) and minimum (I f ,min ) values as [2] The U ∼ I f curves are shown in Figure 3b. It can be clearly seen that the lower the laser power density of the excitation laser, the stronger the curve bending. After normalization, the U ∼ I f curves under each I max only lie in the different bending. In the experiments, we use the photobleaching time constant (τ) to evaluate the strength of the photobleaching. τ can be approximately calculated by nonlinear fitting on U ∼ I f through the equation below [2] where a is an amplification factor which dominates the slope of the U ∼ I f curve and b represents the initial fluorescence intensity of I f at U = 0. τ determines the curvature of the curve. According to Equation (14), a set of U ∼ I f curves are nonlinearly fitted (as shown in Figure 3b) to calculate the corresponding τ. The results of τ calculated from the experiments are listed in Table 2. It can be seen, as P is increased, that τ apparently decreases to as low as 2.43 µs. Subsequently, the numerically calculated U ∼ I f curve is compared with that of the experiments. In the simulation, I f = D I f ,local dV, where D is a cylindrical integration region (i.e., detection region) for fluorescent collection, with a height of h f and diameter of d cl . In this section, D is the cylindrical part of the inset of Figure 2a and d cl = d f , to be coincident with the experimental system we applied. As shown in Figure 3c,d, both the original and normalized numerically calculated fluorescence curves show high consistency with the experimental ones. Since the numerical calculated U ∼ I f curves are merely determined by the photochemical reaction coefficient K, accordingly, a relationship between τ and K can be established, as shown in Figure 4 and Table 2. When I max is increased, we find K increases with I max in an approximately linear manner, while K decreases with τ as K ∼ τ −2.28 , as shown in the inset of Figure 4.
This indicates that the fluorescent photobleaching time shows a nonlinear relationship with K under strong laser illumination. To achieve a higher temporal resolution with smaller τ, K should be significantly increased.  This indicates that the fluorescent photobleaching time shows a nonlinear relationship with under strong laser illumination. To achieve a higher temporal resolution with smaller , should be significantly increased.

Effective Concentration Distribution
The influence of photobleaching on the effective concentration C of fluorescent dye can be directly observed in Figure 5. The photobleaching generates a cometlike region of C. When the flow velocity U at the focus region is small (e.g., U = 4.6 mm/s), as shown in Figure 5a, a strongly photobleached region can be found in the focus region of the laser. The cometlike region of C has a large width attributed to the lateral diffusion. When U is increased, the fluorescent dye is less photobleached and C in the focus region is significantly higher than the counterparts of smaller U, as can be seen from Figure 5b,c. The tails of the cometlike regions exhibit larger length and smaller width. Downstream of the focus area, C gradually recovers due to the molecular diffusion. Similar results can also be found in the second and third rows of Figure 5.

Fluorescence Intensity Distribution
The distributions of fluorescence under different flow velocities and excited b ferent laser power densities are calculated according to Equation (10). The resu shown in Figure 6, where the fluorescent intensity is normalized by the maximum It is obvious that the fluorescent intensity in the focus area is symmetric on the x ax decreases along the flow direction as expected.
As shown in Figure 6a-c, when is increased with fixed 2.96 10 W is increased, while is decreased. The streamwise positions of the peak fluor intensities move downstream and toward the center of the focus area. Similar resu also be found on the second and third rows. When is increased with fixed  The distribution of C is dominated by three dimensionless numbers, i.e., P e and Z pc , in Equation (9), and Z pd = Z pc P e , which is the ratio between photobleaching rate and diffusion rate. When P e is increased with constant Z pc , the influence of diffusivity becomes smaller and the width of the cometlike region is decreased. When Z pc is increased with constant P e , the influence of photobleaching is enhanced and the value of C becomes smaller in the focus area. When Z pd is fixed, the tail length of the cometlike region of C increases with P e .

Fluorescence Intensity Distribution
The distributions of fluorescence under different flow velocities and excited by different laser power densities are calculated according to Equation (10). The results are shown in Figure 6, where the fluorescent intensity is normalized by the maximum value. It is obvious that the fluorescent intensity in the focus area is symmetric on the x axis, and decreases along the flow direction as expected.
As shown in Figure 6a-c, when U is increased with fixed I max = 2.96 × 10 7 W/cm 2 , P e is increased, while Z pc is decreased. The streamwise positions of the peak fluorescent intensities move downstream and toward the center of the focus area. Similar results can also be found on the second and third rows. When I max is increased with fixed U = 41.1 mm/s, P e is constant, while Z pc is increased. From Figure 6c,f,i as examples, it can be seen that the streamwise positions of the peak fluorescent intensities move upstream and even leave the integration region. The fluorescent intensity at the center of focus area continuously decreases.  Figure 6d,g indicate, as is increased, the major body of fluorescence mov stream and even out of the integration region. This could lead to a "waste" of fluore for calculating flow velocity. This reminds us to rethink the influence of the integ region on the performance, especially the spatial resolution of the LIFPA system. fore, we revisited the spatial resolution of LIFPA first in this section. Then, we di possible method to break the optical diffraction limit and realize a super-resolution ity measurement.

Spatial Resolution of Effective Velocity Measurement with LIFPA
To investigate the details of the photobleaching process of LIFPA, we calcula fluorescent intensity profile , ( , , / max , ) through the ce  Figure 6d,g indicate, as Z pc is increased, the major body of fluorescence moves upstream and even out of the integration region. This could lead to a "waste" of fluorescence for calculating flow velocity. This reminds us to rethink the influence of the integration region on the performance, especially the spatial resolution of the LIFPA system. Therefore, we revisited the spatial resolution of LIFPA first in this section. Then, we discuss a possible method to break the optical diffraction limit and realize a super-resolution velocity measurement.

Spatial Resolution of Effective Velocity Measurement with LIFPA
To investigate the details of the photobleaching process of LIFPA, we calculated the fluorescent intensity profile I f ,local ( I f ,local = I f ,local / max y=z=0 I f ,local ) through the center of the laser focus spot (−406 ≤ x ≤ 406 nm), as shown in Figure 7. It can be seen, when K = 1.85 × 10 5 1/s (Figure 7a), the I f ,local curves according to different U show a small difference. As K is further increased, as shown in Figure 7b,c, the fluorescent intensity at the center part becomes smaller and the streamwise positions of the peak I f ,local move downstream at larger U. One noteworthy phenomenon is that when K is sufficiently large (see Figure 7c,d as an example), there is a valley with flat I f ,local , where the fluorescent dye molecules have been strongly photobleached. When measuring velocity fluctuations, only the part marked with a width of d eff (i.e., actual resolution, which is defined as the distance from the left edge of the integration region to the position where I f ,local = 0.05 I f ,local,1 , I f ,local,1 is marked in Figure 7c and represents the I f ,local at the left edge of the integration region) in Figure 7d has an apparent contribution to the fluorescent variation. The other part in the integration region has a negligible contribution to the fluorescent variation, and accordingly, the velocity variation, especially at higher velocity. In other words, the velocity measurement by LIFPA is only sensitive in the marked part, which can be considered as the point spreading function (PSF) of the fluorescent intensity. Therefore, the width of PSF, i.e., d eff , actually determines the spatial resolution of LIFPA. other words, the velocity measurement by LIFPA is only sensitive in the marked part, which can be considered as the point spreading function (PSF) of the fluorescent intensity. Therefore, the width of PSF, i.e., , actually determines the spatial resolution of LIFPA. This is intriguing, since in conventional LIFPA techniques, the spatial resolution of LIFPA is believed to be equivalent to the spatial resolution of the optical system (i.e., the FWHM of the point spread function). From the current research, it can be seen when is sufficiently large, e.g., = 8 × 10 1/s, the averaged over the four curves shown in Figure 7d is around 58 nm. The actual spatial resolution of LIFPA can be apparently smaller than that of the optical system (203 nm in this investigation) and the diffraction limit, i.e., the super-resolution of velocity measurement can be realized. Specifically, this can be achieved by simply increasing in a confocal microscope, or properly selecting a fluorescent dye with large . For instance, if = is increased to 2 × 10 1/s, the averaged can theoretically be as low as 30 nm. Relative to the LIFPA system, developed on the basis of stimulated emission depletion (STED) super-resolution techniques [6] to realize super-resolution velocity measurements, where two beams (one excitation beam and one depletion beam) must be collimated and well aligned, the current investigation shows that LIFPA can be intrinsically super-resolution even with a simple confocal microscope, if is sufficiently large. Additionally, in the STED LIFPA system, the power density of the excitation beam cannot be very high, otherwise the spatial resolution of STED is reduced. Thus, the temporal resolution of the STED LIFPA system is normally low. However, for a sufficiently large , both ultrahigh spatial and temporal resolutions can be achieved simultaneously.

Influence of Integration Region on Velocity Measurement
In this section, we further study the influence of of integration region on the ve- This is intriguing, since in conventional LIFPA techniques, the spatial resolution of LIFPA is believed to be equivalent to the spatial resolution of the optical system (i.e., the FWHM of the point spread function). From the current research, it can be seen when K is sufficiently large, e.g., K = 8 × 10 6 1/s, the averaged d eff over the four curves shown in Figure 7d is around 58 nm. The actual spatial resolution of LIFPA can be apparently smaller than that of the optical system (203 nm in this investigation) and the diffraction limit, i.e., the super-resolution of velocity measurement can be realized. Specifically, this can be achieved by simply increasing I max in a confocal microscope, or properly selecting a fluorescent dye with large k b . For instance, if K = I max k b is increased to 2 × 10 7 1/s, the averaged d eff can theoretically be as low as 30 nm.
Relative to the LIFPA system, developed on the basis of stimulated emission depletion (STED) super-resolution techniques [6] to realize super-resolution velocity measurements, where two beams (one excitation beam and one depletion beam) must be collimated and well aligned, the current investigation shows that LIFPA can be intrinsically superresolution even with a simple confocal microscope, if K is sufficiently large. Additionally, in the STED LIFPA system, the power density of the excitation beam cannot be very high, otherwise the spatial resolution of STED is reduced. Thus, the temporal resolution of the STED LIFPA system is normally low. However, for a sufficiently large K, both ultrahigh spatial and temporal resolutions can be achieved simultaneously.

Influence of Integration Region on Velocity Measurement
In this section, we further study the influence of d cl of integration region on the velocity calibration curve, and attempt to show another way to improve LIFPA's spatial resolution. The results are plotted in Figure 8. Relative to Figure 3a where d cl = d f , the velocity calibration curves calculated at d cl = 0.5d f show clearly smaller values and the value ranges of I f , as plotted in Figure 8a, in line with our expectations. Although it could lead to a negative influence on velocity measurement due to worse signal-to-noise ratio (SNR), the spatial resolution could be clearly improved to super-resolution level. From Figures 3b and 8b,c, it can be seen that changing d cl will not affect the result that the curves show more bending at lower I max . As d cl is increased, the influence of I max on the bending of velocity calibration curves becomes smaller, as can be seen from Figure 8b,c. The U ∼ I f plots become indistinguishable between I max = 5.05 × 10 7 and 7.80 × 10 7 W/cm 2 in Figure 8c. We plot four U ∼ I f curves of different d cl in Figure 8d at I max = 7.80 × 10 7 W/cm 2 . Despite the smaller fluorescent signal at smaller d cl , the velocity calibration curve becomes more linear at smaller d cl . It again indicates the temporal resolution is also improved, accompanied by the spatial resolution. Figures 3b and 8b,c, it can be seen that changing will not affect the result that the curves show more bending at lower . As is increased, the influence of on the bending of velocity calibration curves becomes smaller, as can be seen from Figure  8b,c. The ~ plots become indistinguishable between 5.05 10 and 7.80 10 W/cm 2 in Figure 8c. We plot four ~ curves of different in Figure 8d at 7.80 10 W/cm 2 . Despite the smaller fluorescent signal at smaller , the velocity calibration curve becomes more linear at smaller . It again indicates the temporal resolution is also improved, accompanied by the spatial resolution. In addition, we plotted the , ~ under 0.5 in Figure 9. It can be seen that as is reduced, the corresponding is also reduced. At 8.00 10 1/s, the averaged over the four curves is around 55 nm, which is slightly smaller than that at shown in Figure 7d. Thus, by decreasing the size of the integration region, the actual spatial resolution of LIFPA can be further enhanced. In addition, we plotted the I f ,local ∼ x under d cl = 0.5d f in Figure 9. It can be seen that as d cl is reduced, the corresponding d eff is also reduced. At K = 8.00 × 10 6 1/s, the averaged d eff over the four curves is around 55 nm, which is slightly smaller than that at d cl = d f shown in Figure 7d. Thus, by decreasing the size of the integration region, the actual spatial resolution of LIFPA can be further enhanced. In addition, we plotted the , ~ under 0.5 in Figure 9. It can be seen that as is reduced, the corresponding is also reduced. At 8.00 10 1/s, the averaged over the four curves is around 55 nm, which is slightly smaller than that at shown in Figure 7d. Thus, by decreasing the size of the integration region, the actual spatial resolution of LIFPA can be further enhanced.

Conclusions
In this investigation, we studied the photobleaching process of fluorescent dye in the laser focus region, according to the model based on the convection-diffusion reaction equation. The profiles of effective dye concentration and fluorescence were elucidated. The relationship between the commonly used photobleaching time constant obtained by experiments and the photochemical reaction coefficient was revealed. We further studied the influence of the detection region of fluorescence on the performance of the LIFPA system, and found that at sufficiently high excitation laser power density, LIFPA can achieve super-resolution and break the limit of optical diffraction, even in a confocal microscope. We hope the current investigation can promote the development of LIFPA, and reveal the photobleaching process of fluorescent dye under high laser power density illumination, to enhance our understanding of fluorescent dynamics and photochemistry and develop more powerful photobleaching-related flow diagnostic techniques.