Separation efficiency maximization in acoustofluidic systems: study of the sample launch-position

The development of lab-on-chip microfluidic systems based on acoustic actuation, and in particular on the acoustophoretic force, has recently attracted significant attention from the scientific community thanks, in part, to the possibility of sample sorting on the basis of both geometrical and mechanical properties. It is commonly recognized that sample prefocusing and launch-position optimization have a substantial effect on the performance of these systems but a clear explanation of how these two parameters influence the system efficiency is still missing. In this manuscript we discuss the impact of both the sample launch position and the sample distribution at the input by the theoretical analysis of a simplified system and by numerical simulations of realistic configurations. The results show that the system performance can be greatly improved by selecting the proper microchannel dimensions and sample-launch position, offering relevant guidelines for the design of micro-acoustofluidic lab-on-chip devices.


Introduction
Over the past twenty years, microuidic and lab-on-chip techniques have generated signicant research interest and successfully permeated into many different elds, with particular attention being paid to lab-on-chip systems for cytology applications. 1 Their inherent micro-size characteristic provides a natural environment to detect, and manipulate cells, even for analysis at a single cell level. 2 Multiple functionalities have already been integrated within a single chip, thanks to the inclusion of different sensors and sample-actuation mechanisms, such as those based on dielectrophoresis, 3,4 optical forces, 5-7 cavitation bubbles 8,9 and magnetic forces. 10,11 An actuation system which is currently attracting considerable attention is that based on the interaction between the sample and acoustic waves. On-chip acoustouidics, which combines the use of ultrasonic acoustic waves with the advantages of microuidic systems, has become an extremely active eld and several review papers, 12,13 tutorials 14 and books 15 have been dedicated to this eld.
Acoustouidics has been successfully applied to many different research studies, ranging from micro-droplet production and manipulation to micro-particle (or cell) sorting, focusing, separation, mixing and arraying. [16][17][18][19][20][21][22][23][24][25][26][27] Among these possible applications, cell separation has attracted a lot of attention including substantial effort devoted towards enabling the isolation of circulating tumor cells from human blood samples. [28][29][30][31][32][33][34][35][36] Several groups have demonstrated the possibility of isolating target cells from a given sample containing a mix of cells with different characteristics using acoustouidics. 37,38 A theoretical analysis of particle separation efficiency in acoustophoretic devices was recently reported and it showed that both intrinsic factors, related to the sample itself, and extrinsic factors, related to the microuidic system, can strongly affect the separation result. Among all the factors, the sample launch position into the active region (i.e. where the acoustophoretic force is present) plays a critical role. 39 The aim of this paper is to investigate how the launch position of the sample inside the channel changes the separation efficiency, and to demonstrate how the best launch position depends on different parameters: the aspect ratio of the channel, the cross-section occupied by the sample distribution at the channel input and the radius of the target bead/cell sample. In the following we focus our attention on investigating the challenging situation where the micro-objects to be selected and separated show a small deviation of their properties from the other beads/cells owing along the channel.

Background knowledge
This section reports some of the basic information about system geometry along with the fundamental equations that determine the performance of the congurations considered in the paper.

Description of analyzed conguration
We take into consideration the simple situation of a rectangular microchannel, with a cross-sectional area of 0.09 mm 2 and with variable aspect ratio (dened as the ratio between width and height, w/h), where an acoustic standing-wave (with a single pressure node and resonating along the microchannel width 40 ) is introduced. We consider a sample-injection area different from the main inlet of the buffer uid, positioned in the "righthalf" of the microchannel, as illustrated graphically in Fig. 1, and we analyze the sample movement across the channel (towards the channel's center) due to the acoustic radiation force caused by the scattering of the acoustic wave by the particles. In real-world situations the sample-injection area yields a certain statistical distribution of the sample starting position (y 0 ,z 0 ); we thus want to analyze the impact of the initial position uncertainty (Dy 0 ,Dz 0 ) on the achievable separation between two sample populations with slightly different properties.
It is worth underlining that in our analysis we neglect the effect of gravity, which can be helpfully used to produce a sample separation in the vertical direction, 41 and the acoustic streaming-induced drag force (generally relevant for particles much smaller than cells).

Fundamental equations and parameter denitions
The reference system used throughout the paper is shown in Fig. 1: the uid ows in the x direction, the acoustic wave resonates along the y-direction (corresponding to the channel width) and the microchannel height is in the z-direction.
One of the main parameters affecting the movement of particles exposed to acoustic waves is the so called acoustic contrast factor (4), which is given by the following equation where r p and r f are the densities of the suspended microparticles and uid, respectively, and b p and b f are the corresponding compressibilities. 40 The acoustic contrast factor is one of the main quantities appearing in the expression of the acoustophoretic force applied to a microsphere in a plane standing wave, see eqn (2), where R is the beads/particles radius, y is its position in the transverse direction, E ac is the acoustic energy density in the microchannel, and k y represents the acoustic wave number.
By imposing the Stokes drag force to be equal to the acoustic force it is possible to calculate the transverse coordinate of the particle as a function of time, as shown in eqn (3). 40 In that equation we identify as y 0 the position, along the y-axis, occupied by the particle when it enters the area of the microchannel where the acoustic wave is present. (3) Regarding the particle movement in the x-direction, we assume that it moves together with the uid, whose velocity v x (y,z) in the microchannel is given by eqn (5). 39 To derive that equation we assume a rectangular microchannel of height h, width w, we identify with VP the pressure gradient and with v c,N a constant factor, dened by eqn (4).
Regarding the z-direction, conversely, we assume that no signicant movement occurs during the time that the microbeads ow along the microchannel. We also initially set Dz 0 ¼ 0, even if we will remove this assumption in Section 4.2.
The above equations allow us to calculate the movement in 3D of a particle owing in the considered microuidic system, and thus they also allow us to dene the bandwidth (BW), displacement (D) and separation-efficiency (SE) parameters, as already reported in ref. 39. To briey recall those denitions, we identify as D i (x) the distance traveled by beads belonging to the i-th population along the y-direction with respect to the starting position, y 0 , and we name BW i (x) the spreading, in the ydirection, of the i-th population.
Starting from these two quantities, the SE parameter is calculated, at each position in the x direction, by the ratio of two distances as shown by eqn (6). As a numerator we use the difference between the distances traveled by two beads, with different characteristics, and injected in the microchannel at Fig. 1 3D scheme of the particle separation mechanism considered in this study. The blue and the green volumes represent the positions occupied by beads with different acoustic contrast factors (80% and 100% of that associated with polystyrene beads in water) while they flow along the x-direction. It should be noted that this figure is only used to introduce the analyzed configuration and the reference system. The z-dependence of the fluid velocity profile has not been taken into account in the calculation of the shown trajectories. the same position y 0 . Conversely at the denominator we sum the two single-side bandwidths with the bead radius. As a result the SE parameter calculated by eqn (6) is greater than 1 when the two populations are completely separated by a distance larger than one bead radius (see Fig. 2).
An analysis of the intrinsic (i.e. sample-related) and extrinsic (system-related) parameters affecting the SE can be found in ref.

39.
3 Theoretical analysis: uniform flowspeed We start our analysis by noticing from eqn (5) that every time the distance from the border is larger than twice the channel's height (i.e. y > 2h) the ratio between the two hyperbolic cosine functions in eqn (5) becomes negligible. As a consequence, in this condition the uid speed through the channel cross-section does not depend on the y-coordinate. When a large aspect-ratio microchannel (i.e. having h ( w) is considered, the area where the ow speed depends on the y-coordinate (y < 2h) represents a negligible portion of the cross-section and it is thus possible to assume that the beads movement along the x-direction does not depend on the position along the y-direction. As a consequence, the sample separation depends on eqn (3) and we can thus consider such a conguration as a 1D-system. To simplify the writing of eqn (3), but without losing important dependences, we can introduce the following denitions: This allows us to rewrite eqn (3) as: We use this compact form to develop an analytical study of the simplied 1D-system which, although an over-simplication, provides useful insights into the impact of different parameters on the separation efficiency (SE), as dened in eqn (6). Using this notation, a change in the particle properties is reected by a change of the A parameter, and it is thus possible to dene an average A and a deviation DA. Similarly, as anticipated in Section 2, Dy 0 represents the uncertainty of the input position, i.e. the maximum variation from the desired value y 0 .

Performance analysis of 1D-systems
In order to derive the SE dependence on the A and y 0 parameters we consider two different populations (a "fast" one and a "slow" one; F/S) and for each of them we consider three different "injection points": the center of the injection area (C), the border of the injection area closer to the microchannel wall (W) and the border closer to the node of the acoustic wave (N). We can thus write 6 equations similar to eqn (8), where we modify A in A AE DA to include the effect of the "fast" and "slow" population, and we use y 0 , y 0 À Dy 0 and y 0 + Dy 0 to take into account the C, W and N starting positions respectively (see inset of Using this denition, it is possible to rewrite eqn (6) as: If we then approximate the differences appearing in the eqn (10) with the corresponding rst-order differential terms we obtain: 100% of that associated with polystyrene beads in water) while they flow along the x-direction. The inset shows a zoomed-in view of the injection region ("C" stands for the center of the injection area, "W" stands for the border of the injection area closer to the microchannel wall and "N" stands for the border closer to the node of the acoustic wave). Right scale: the SE calculated according to eqn (6), and as a function of the x coordinate, shows a peak at about 800 mm where the two populations are well separated.
The above equation can be rewritten in a more useful form by calculating the two derivatives so as to obtain an explicit expression for the SE evolution as a function of time, or of own-distance, as we consider uniform ow speed in the channel. If we then set DA ¼ DA r Â A, where DA r is the relative variation of the A parameter we obtain: In the following we identify with an asterisk (*) the values yielding the best performance, thus y * 0 and t* are the launchposition and time-instant corresponding to the maximum separation efficiency (SE*) achievable between two populations. In particular, thanks to eqn (12), it is possible to derive some non-trivial considerations about y * 0 , t* and SE*: C.(1) The role of DA r : With the given hypothesis the SE(t) (and hence SE*) is directly proportional to DA r , but DA r has no effect on the y 0 value maximizing the SE. It thus means that the optimal launch position of the sample in the channel ðy * 0 Þ doesn't depend on how large the sample property variations are.
C. (2) The maximum-separation instant (t*): The time instant giving the maximum SE cannot be explicitly calculated as the resulting equation is transcendental, but we observe that in eqn (12) the time always appears in the product At. This implies that the maximum-separation instant t* is inversely proportional to A, i.e. the product At* does not depend on A. It is thus useful to rewrite eqn (12) in a more compact way by dening the following terms: The SE achieved at t* can thus be calculated as shown in eqn (14), which has two important characteristics: E is a constant while both F and G only depend on the value of Dy 0 /R.
C. (3) The ratio Dy 0 /R: Following the above considerations we observe that if both Dy 0 and R are multiplied by the same quantity, SE (t*) is le unmodied and thus also the optimal launch position ðy * 0 Þ doesn't change. Nevertheless, it is important to highlight that even if we keep Dy 0 /R constant, a variation of R affects both the value of t* and that of SE(t*) as they are proportional to R À2 and R À1 .
C.(4) The y * 0 value: In this simplied situation y * 0 depends only on the ratio Dy 0 / R. By analyzing Fig.(3) it is interesting to notice that the ideal injection position is close to the border when the ratio Dy 0 /R is almost zero, it rapidly increases when Dy 0 /R grows from 0 to 1 (going from 0% to 12% of w) and then it almost saturates around 20% of w when Dy 0 /R gets much larger than 1.
Aer this study of the 1D system it becomes relevant to analyze, by means of numerical simulations, how the system properties are affected by moving to a 2D and 3D system, so as to verify if the above considerations are still valid.

Numerical simulations: results and discussion
The comparison with a more realistic situation, by means of numerical simulations, is carried out in a two-step process. As a rst step (Section 4.1) we introduce the presence of a nonuniform ow-speed prole across the channel width (y-direction using the reference system shown in Fig. 1) and we compare the results by varying the aspect-ratio (w/h) of the microchannel cross-section. In this case the analysis of the results obtained at large aspect ratio can also be used to assess the validity of the results obtained analytically in the previous section. In this "2D-case" we neglect the dependence of the ow speed (v x ) on the z-direction which is equivalent to assuming that the sample is injected at the half-height of the microchannel and neglecting the vertical sample dispersion. Subsequently, as the second step (Section 4.2), we also include in our model the dependence of v x on the vertical direction v x (y,z) and the vertical-dispersion of the sample at the microchannel input (Dz 0 ), so as to mimic a realistic situation. This allows us to investigate the impact on the achievable SE of the vertical position and spreading of the sample inlet, allowing us to derive some interesting design rules. In both cases we carry out the simulations using typical parameter values corresponding to a water-suspension of polystyrene (PS) microbeads. The parameter values are reported in Table 1. All the simulations were carried out using a custom MATLAB script based on the equations reported in Section 2.2. The beads trajectory at the x-y (2D systems) and x-y-z (3D systems) coordinates as a function of time were calculated combining eqn (3) and (5), thanks to the use of "ODE45" function of MATLAB. The trajectories were calculated over a time vector, composed by 2000 uniformly spaced values, whose time-step depends on microbeads properties and on the channel geometry, as they impact on the time required by beads to reach the microchannel center. As a general indication, the calculation of each point appearing in the gures from 4 to 6 required about half an hour of computation time on a 4-cores processor at 3.50 GHz and with 16 GB of RAM.
It is possible to investigate different microchannel aspectratios in two different ways: keeping one dimension xed and varying the other one (e.g. xed height and variable width), or by simultaneously modifying both dimensions so as to keep the cross-sectional area constant. In the following, unless otherwise specied (in the nal part of Section 4.2), we keep the microchannel area xed at 9 Â 10 À2 mm 2 , corresponding to a square microchannel with a side of 300 mm. The dependencies and trends observed for this specic area are of general validity.

2D-systems approximation
The rst phase of these simulations always involved the iden-tication of the proper time-interval to be considered and of the required spatial-and temporal-resolution. For each parameterset of the microchannel-beads system, we considered different launch positions (y 0 ) and for each of them we derived the corresponding SE, thus allowing us to identify the optimal launch position ðy * 0 Þ. The numerical analysis was initially carried out to verify whether or not the considerations expressed in Section 3 are still valid in the 2D case. In order to verify the previously reported consideration C.(1), we considered in our numerical simulations distinct beads populations having a 4 (see eqn (1)) different from that of standard PS. In Fig. 4 we show the y * 0 values as a function of the microchannel aspect-ratio (w/h) and considering four different bead-population pairs.
The rst pair is composed by beads having an acoustic contrast factor (4) equal to 98% and 102% of the nominal PS value, and it is thus indicated as AE2% in the legend (corresponding to the DA r value in Section 3). In an analogous way, we considered population-pairs with an increased difference of 4, up to the AE20% case, which corresponds to bead-populations having an acoustic contrast factor equal to 80% and 120% respectively of that of PS. As reported in Section 3, it is possible to notice that the y * 0 position doesn't depend on how large the sample properties variations are, provided that the difference is not too big. As it is evident considering the AE20% case, if the 4 variation becomes too large the y * 0 position may start to vary, as the rst-order approximation used to derive eqn (11) is not sufficient anymore.
Nevertheless, as the most critical separation situation is when small differences are present between the sample populations, this limitation is not particularly relevant for our study. It is also interesting to notice that the y * 0 value depends on the aspect ratio of the microchannel cross-section (indicated as w/h in the gures), but it becomes almost constant when w > 10h, as it approaches the 1D-situation theoretically analyzed in Section 3. As a comparison it is interesting to notice that data used to create Fig. 4 yield a ratio Dy 0 /R z 2.7 which corresponds in Fig. 3 to y * 0 z 15%, exactly matching the numerically calculated position for large aspect-ratios.
Subsequently we also veried that the y * 0 value does not depend on the absolute value of 4. To analyze this aspect, we considered three different bead-population pairs, with largely different values of the nominal 4 factor (0.05, 0.5 and 5), while keeping D4 ¼ AE5%. The results show that no change of the y * 0 value is induced by modifying the nominal 4 of the populations, independently of the channel aspect ratio (see Fig. 5).
We then moved to verify the dependence of y * 0 on the Dy 0 and R parameters (i.e. consideration C.(3) of Section 3). According to what previously reported we expect the y * 0 value to depend on  Dy 0 /R, but not on Dy 0 and R separately. To assess this dependence, we compared the ideal y * 0 for ve different congurations of Dy 0 and R, while keeping constant the acoustic contrast factor of the two populations (4 ¼ 0.5; D4 ¼ AE5%). The ve parameterset considered in the numerical simulations are schematically reported in Table 2. We included three different combinations of Dy 0 and R yielding the same Dy 0 /R ratio (8/3), and two different combinations yielding a four-times increase and decrease of the Dy 0 /R value (32/3 and 2/3 respectively).
The results clearly highlight that even in this case, as in the 1D situation previously considered, the optimal launch position depends on Dy 0 /R, and is thus unmodied if both values are multiplied by the same factor (see Fig. 6).
To complete the 2D-approximation analysis of the launch position we created a gure to show the overall dependence of the y * 0 parameter on Dy 0 /R and w/h, which are the only two parameters affecting the y * 0 value. The result of the numerical simulations is reported in Fig. 7 (le panel) as a color-map. Calculations were carried out considering a nominal 4 ¼ 0.5, D4 ¼ AE5% and R ¼ 5 mm, but the reported results have a much more general validity as derived by the above reported analysis.
In the same conditions we also calculated the SE* value, dened according to eqn (6), achievable by proper selection of the injection point, as a function of Dy 0 /R and w/h. The obtained results, reported in right panel of Fig. 7, show the benet of using large aspect ratios and the advantages given by a reduction of Dy 0 , which can be achieved, as an example, by using a prefocusing section. 42 The results of this analysis demonstrate two important aspects: (i) the importance of optimizing the launch position ðy * 0 Þ to improve the SE and (ii) that the y * 0 value, expressed as a percentage of the microchannel width, depends only on microchannel aspect ratio (w/h) and on the Dy 0 /R ratio, but not on other sample factors. Anyway, it is important to highlight that the whole analysis reported up to this section completely neglects the vertical dimension of the microchannel.

Extension to 3D-systems
In this section we analyze a realistic 3D-system, including also the role played by the microchannel vertical dimension, and we thus add in our model two important effects: the dependence of ow speed on the z-coordinate and the fact that the sample distribution has a non-zero dimension also along the vertical direction. As a rst step, we need to properly re-dene the SE parameter: as beads owing at different heights have different velocities, even if we assume to have a rectangular distribution of beads at the input, the positions occupied while owing along the microchannel produce a curved distribution of beads at any other section, as it can be seen in Fig. 8. Since we are interested in having a realistic evaluation of the system performance, we assume the cross-section of the "sampleextraction" port to have a rectangular shape, independently of the curved areas corresponding to the beads positions. We thus keep the denition of the SE unmodied, as in eqn (6), and we include the effect of the curved beads-distribution by re-dening the BW values and |D 1 À D 2 | (see Fig. 8).
In particular we dene the BW parameter of each population as the maximum distance (in the y-direction) between two beads, i.e. considering at each section along the x-axis the bead closer to the microchannel wall (at half of the channel height) and the one closer to the microchannel center (and closer to microchannel oor). The distance |D 1 À D 2 |, which represents    Table 2. Other simulation parameters: 4 ¼ 0.5; D4 ¼ AE5%.
the distance between the "centers" of the two beads distributions, is calculated as the distance between the centers of the two population-bands.
It is important to notice that, as it is evident by Fig. 8, choosing an injection height different from the center of the channel can only worsen the system performance. The owspeed gradient in fact becomes larger as we move away from the middle-height position and thus the beads distribution becomes wider. For this reason, we considered in our analysis the impact of a vertical spreading Dy 0 while keeping the center of the injection channel xed at half-height of the microchannel. Regarding the Dy 0 parameter it should be noticed that it is possible to dene it as a given percentage of the microchannel height, or by its own value (in mm).
To keep consistency with the previous analysis we initially consider the case of Dy 0 dened as a xed percentage of the microchannel height. We decided to start our analysis considering a Dz 0 equal to 5% of the channel height, while keeping all the other parameters set as for the nal 2D simulations: nominal 4 ¼ 0.5, D4 ¼ AE5% and R ¼ 5 mm. As in the previous case we calculated the optimal launch position y * 0 , and the corresponding separation efficiency (SE*) as a function of the microchannel aspect ratio and of the Dy 0 /R value.
The data reported in Fig. 9 (le panel and right panel) show two partially surprising results: the y * 0 obtained in the 3D case exactly matches that obtained in the 2D approximation and also the SE* gure matches that obtained in the 2D case, once rescaled by a constant factor. The reason for these results is that the presence of a vertical spreading implies the presence of beads owing at a different height, where the ow speed is simply scaled (by a factor smaller than 1) with respect to the ow speed at half-height of the channel. As no distortion of the speed prole is introduced, the y * 0 value for beads owing at half height and for those owing at any distance from the channel bottom is the same, provided that beads' interaction with the bottom surface can be neglected. A direct consequence of this is that even analysis carried out using a larger vertical spreading (e.g. Dz 0 equal to 10% or 15%) would yield the same results and thus do not bring any additional information.
The obtained values show that even a minor vertical spreading can have a signicant impact on the achievable SE in case of microchannels with high aspect ratio: as an example, a 5% vertical spreading in the microchannel with aspect ratio 25  corresponds to a Dz 0 as small as AE3 mm, and yields a SE* reduction almost by a factor of 2. On the other side, the use of microchannels with a smaller aspect ratio, although yielding a lower SE* value in the ideal case of Dz 0 ¼ 0, is expected to be signicantly more tolerant to the vertical spreading. We thus investigated the performance of microchannels with different aspect-ratios while xing Dz 0 equal to AE5 mm and AE10 mm.
The y * 0 color-maps obtained in these conditions do not add any relevant information with respect to Fig. 9 and they are thus not reported in the manuscript. Conversely, it is interesting to analyze the data reported in Fig. 10, showing the achievable SE in the above described conditions and with a xed value of Dz 0 . The reported maps highlight that, once Dy 0 /R and Dz 0 are given, it is possible to identify the ideal channel cross-section and then the achievable SE.
It is interesting to notice that while in the 2D case an aspect ratio as large as possible was desirable (see right panel of Fig. 7), in the 3D case the presence of a non-negligible Dz 0 suggests the use of higher channels, so as to mitigate the effect of the vertical spreading of the sample. As a consequence of the necessity to Fig. 9 Left: consider the sample having a certain vertical distribution around the half-height of the channel. The vertical spread is equal to AE5% of the channel height. Optimal injection position y * 0 as a function of the ratio Dy 0 /R and of the channel aspect ratio in the 3D case. Right: corresponding maximum value of SE* as a function of the ratio Dy 0 /R and of the channel aspect ratio in the 3D case. Fig. 10 Left: consider the sample having a certain vertical distribution around the half-height of the channel. The vertical spread is equal to AE5% mm. SE* as a function of the ratio Dy 0 /R and of the channel aspect ratio in the 3D case. Right: consider the sample having a certain vertical distribution around the half-height of the channel. The vertical spread is equal to AE10 mm. SE* as a function of the ratio Dy 0 /R and of the channel aspect ratio in the 3D case.
nd a trade-off between the mitigation of horizontal and vertical spreading, the ideal aspect ratio has a non-obvious dependence on both Dy 0 /R and Dz 0 .
Additionally, it is worth mentioning that in the above reported discussion we exclusively focused our attention on the acoustic radiation force, applied on the owing particles because of the sound-waves scattering, while we neglected the acoustic streaming effect and the related drag-force. Following the analysis reported by Muller et al. 43 it is possible to show that, in order to neglect the acoustic streaming effect, the particle diameter must exceed by a few times the boundary-layer thickness d.
where J is a parameter related to the microchannel geometry and that in case of planar walls is equal to 3/8. With some simple passages, this condition can be written as a limitation on the microchannel width (w) expressed by the below reported equation, where h represents the dynamic viscosity of the uid, r is the uid density and c 0 is the sound-wave speed in the uid: As a reference value, considering cells with 4 ¼ 0.15 and R ¼ 3 mm the maximum value of w allowing to neglect the acoustic streaming effect is z1000 mm, which may thus impose a limitation on the achievable microchannel aspect-ratio (w/h).

Conclusions
In this work we described the results of an analytical and numerical investigation regarding the acoustic separation of microbeads. In particular, we focused on the study and optimization of the sample-launch position in order to maximize the separation efficiency in the challenging situation where the micro-objects to be separated are characterized by a small deviation of their properties from the other micro-particles owing along the channel. We showed that the best sample injection position depends on different factors (the aspect ratio of the channel, the cross section occupied by the beads distribution at the sample inlet and the radius of the bead). The optimization method presented in this study allowed us to derive important design rules which can be applied to free-ow microuidic separation systems independently on the presence or absence of pre-focusing stages and strategies.
It is interesting to notice that thanks to a careful optimization of the injection position of the sample, high SE values can be obtained even in case of no prefocusing techniques. This allows largely simplifying the design and the operation of the microuidic systems.

Conflicts of interest
There are no conicts of interest to declare.